Distribution of graph-distances in Boltzmann ensembles of RNA secondary structures

Large RNA molecules often carry multiple functional domains whose spatial arrangement is an important determinant of their function. Pre-mRNA splicing, furthermore, relies on the spatial proximity of the splice junctions that can be separated by very long introns. Similar effects appear in the processing of RNA virus genomes. Albeit a crude measure, the distribution of spatial distances in thermodynamic equilibrium therefore provides useful information on the overall shape of the molecule can provide insights into the interplay of its functional domains. Spatial distance can be approximated by the graph-distance in RNA secondary structure. We show here that the equilibrium distribution of graph-distances between arbitrary nucleotides can be computed in polynomial time by means of dynamic programming. A naive implementation would yield recursions with a very high time complexity of O(n^11). Although we were able to reduce this to O(n^6) for many practical applications a further reduction seems difficult. We conclude, therefore, that sampling approaches, which are much easier to implement, are also theoretically favorable for most real-life applications, in particular since these primarily concern long-range interactions in very large RNA molecules.Comment: Peer-reviewed and presented as part of the 13th Workshop on Algorithms in Bioinformatics (WABI2013

Isothermal folding of G-quadruplexes

Thermodynamic studies of G-quadruplex stability are an essential complement to structures obtained by NMR or x-ray crystallography. An understanding of the energetics of quadruplex folding provides a necessary foundation for the physical interpretation of quadruplex formation and reactivity. While thermal denaturation methods are most commonly used to evaluate quadruplex stability, it is also possible to study folding using isothermal titration methods. G-quadruplex folding is tightly coupled to specific cation binding. We describe here protocols for monitoring the cation-driven quadruplex folding transition using circular dichroism or absorbance, and for determination of the distribution of free and bound cation using a fluorescence indicator. Together these approaches provide insight into quadruplex folding at constant temperature, and characterize the linkage between cation binding and folding

Magnesium Fluctuations Modulate RNA Dynamics in the SAM-I Riboswitch

Experiments demonstrate that Mg(2+) is crucial for structure and function of RNA systems, yet the detailed molecular mechanism of Mg(2+) action on RNA is not well understood. We investigate the interplay between RNA and Mg(2+) at atomic resolution through ten 2 microsecond explicit solvent molecular dynamics simulations of the SAM-I riboswitch with varying ion concentrations. The structure, including three stemloops, is very stable on this timescale. Simulations reveal that outer sphere coordinated Mg(2+) ions fluctuate on the same time scale as the RNA, and that their dynamics couple. Locally, Mg(2+) association affects RNA conformation through tertiary bridging interactions; globally, increasing Mg(2+) concentration slows RNA fluctuations. Outer sphere Mg(2+) ions responsible for these effects account for 80% of Mg(2+) in our simulations. These ions are transiently bound to the RNA, maintaining interactions, but shuttled from site to site. Outer sphere Mg(2+) are separated from the RNA by a single hydration shell, occupying a thin layer 3-5Å from the RNA. Distribution functions reveal outer sphere Mg(2+) are positioned by electronegative atoms, hydration layers, and have a preference for the major groove. Diffusion analysis suggests transient outer sphere Mg(2+) dynamics are glassy. Since outer sphere Mg(2+) ions account for most of the Mg(2+) in our simulations, these ions may change the paradigm of Mg(2+)-RNA interactions. Rather than a few inner sphere ions anchoring the RNA structure surrounded by a continuum of diffuse ions, we observe a layer of outer sphere coordinated Mg(2+) that is transiently bound but strongly coupled to the RNA