On the dynamics of some small structural motifs in rRNA upon ligand binding

Abstract

The present study characterizes using molecular dynamics simulations the behavior of the GAA (1186-1188) hairpin triloops with their closing c-g base pairs in large ribonucleoligand complexes (PDB IDs: 1njn, 1nwy, 1jzx). The relative energies of the motifs in the complexes with respect to that in the reference structure (unbound form of rRNA; PDB ID: 1njp) display the trends that agree with those of the conformational parameters reported in a previous study1 utilizing the de novo pseudotorsional (ηθ) approach. The RNA regions around the actual RNA-ligand contacts, which experience the most substantial conformational changes upon formation of the complexes were identified. The thermodynamic parameters, based on a two-state conformational model of RNA sequences containing 15, 21 and 27 nucleotides in the immediate vicinity of the particular binding sites, were evaluated. From a more structural standpoint, the strain of a triloop, being far from the specific contacts and interacting primarily with other parts of the ribosome, was established as a structural feature which conforms to the trend of the average values of the thermodynamic variables corresponding to the three motifs defined by the 15-, 21-and 27-nucleotide sequences. From a more functional standpoint, RNA-ligand recognition is suggested to be presumably dictated by the types of ligands in the complexes.