Abstract:Tetraloops are very abundant structural elements of RNA that are formed by four nucleotides in a hairpin loop which is closed by a double stranded helical stem with some Watson-Crick base pairs. A tetraloop r(GCGAAGGC) was identified from the crystal structure of the central domain of 16S rRNA (727-730) in the Thermus thermophilus 30S ribosomal complex. The crystal structure of the 30S complex includes a total of 104 nucleotides from the central domain of the 16S rRNA and three ribosomal proteins S6, S15 and S18. Independent biochemical experiments have demonstrated that protein S15 plays the role in initiating the formation of the central domain of this complex. In the crystal, the tetraloop interacts with the protein S15 at two sites: one of them is associated with hydrogen bond interactions between residue His50 and nucleotide G730, and the other is associated with the occurrence of residue Arg53 beside A728. This paper uses molecular dynamics (MD) simulations to investigate the protein-dependent structural stability of the tetraloop and demonstrates the folding pathway of this tetraloop via melting denaturation and its subsequent refolding. Three important results are derived from these simulations: (i) The stability of nucleotide A728 appears to be protein dependent. Without the interaction with S15, A728 flips away from stacking with A729. (ii) The melting temperature demonstrated by the simulations is analogous to the results of thermodynamic experiments. In addition, the simulated folding of the tetraloop is stepwise: the native shape of the backbone is formed first; this is then followed by the formation of the Watson- Crick base pairs in the stem; and finally the hydrogen bonds and base stacking in the loop are formed. (iii) The tetraloop structure is similar to the crystal structure at salt concentrations of 0.1 M and 1.0 M used for the simulations, but the refolded structure at 0.1 M salt is more comparable to the crystal structure than at 1.0 M. The results from the simulations using both the Generalized Born continuum model and explicit solvent model (Particle Mesh Ewald) generate a similar pathway for unfolding/refolding of the tetraloop.

RNA secondary structure formation: a solvable model of heteropolymer folding

Statistical mechanics of secondary structures formed by random RNA sequences

Modeling Loop Composition and Ion Concentration Effects in RNA Hairpin Folding Stability

Statistical Physics of RNA-folding

Numerical Performance of Abs Codes for Nonlinear Systems of Equations

RNA Secondary Structure Thermodynamics

Theory for RNA folding, stretching, and melting including loops and salt

Statistical mechanics of RNA folding: a lattice approach

Bistable secondary structures of small RNAs and their structural probing by comparative imino proton NMR spectroscopy

A universality class for RNA-like polymers and double polymers

On the behavior of random RNA secondary structures near the glass transition

RNAlike polymer model: exact calculation on the Bethe lattice

Helix-Based RNA Landscape Partition and Alternative Secondary Structure Determination

A detailed field theory for RNA-like molecules with periodic base sequence

Free Energy Landscape and Multiple Folding Pathways of an H-Type RNA Pseudoknot

Cotranscriptional kinetic folding of RNA secondary structures including pseudoknots

ANNaMo: Coarse-grained modelling for folding and assembly of RNA and DNA systems

Accurately Modeling RNA Stem-Loops in an Implicit Solvent Environment

Critical phenomena of RNA-like polymers -- a synopsis

Molecular Dynamics Simulations of the Denaturation and Refolding of an RNA Tetraloop.

Dynamic Extended Folding: Modeling the RNA Secondary Structures During Co-Transcriptional Folding.