Molecular dynamics (MD) simulations were employed to investigate the structure, dynamics, and local base-pair step deformability of the free 16S ribosomal helix 44 from Thermus thermophilus and of a canonical A-RNA double helix. While helix 44 is bent in the crystal structure of the small ribosomal subunit, the simulated helix 44 is intrinsically straight. It shows, however, substantial instantaneous bends that are isotropic. The spontaneous motions seen in simulations achieve large degrees of bending seen in the X-ray structure and would be entirely sufficient to allow the dynamics of the upper part of helix 44 evidenced by cryo-electron microscopic studies. Analysis of local base-pair step deformability reveals a patch of flexible steps in the upper part of helix 44 and in the area proximal to the bulge bases, suggesting that the upper part of helix 44 has enhanced flexibility. The simulations identify two conformational substates of the second bulge area (bottom part of the helix) with distinct base pairing. In agreement with nuclear magnetic resonance (NMR) and X-ray studies, a flipped out conformational substate of conserved 1492A is seen in the first bulge area. Molecular dynamics (MD) simulations reveal a number of reversible alpha-gamma backbone flips that correspond to transitions between two known A-RNA backbone families. The flipped substates do not cumulate along the trajectory and lead to a modest transient reduction of helical twist with no significant influence on the overall geometry of the duplexes. Despite their considerable flexibility, the simulated structures are very stable with no indication of substantial force field inaccuracies.
- MeSH
- bakteriální RNA chemie genetika MeSH
- hořčík chemie MeSH
- konformace nukleové kyseliny MeSH
- krystalografie rentgenová MeSH
- magnetická rezonanční spektroskopie MeSH
- molekulární modely MeSH
- molekulární struktura MeSH
- párování bází genetika MeSH
- počítačová simulace * MeSH
- RNA ribozomální 16S chemie genetika MeSH
- sekvence nukleotidů MeSH
- sodík chemie MeSH
- termodynamika MeSH
- Thermus thermophilus chemie genetika MeSH
- vazebná místa genetika MeSH
- vodíková vazba MeSH
- Publikační typ
- časopisecké články MeSH
- práce podpořená grantem MeSH
- Názvy látek
- bakteriální RNA MeSH
- hořčík MeSH
- RNA ribozomální 16S MeSH
- sodík MeSH
The hepatitis delta virus (HDV) ribozyme is an RNA enzyme from the human pathogenic HDV. Cations play a crucial role in self-cleavage of the HDV ribozyme, by promoting both folding and chemistry. Experimental studies have revealed limited but intriguing details on the location and structural and catalytic functions of metal ions. Here, we analyze a total of approximately 200 ns of explicit-solvent molecular dynamics simulations to provide a complementary atomistic view of the binding of monovalent and divalent cations as well as water molecules to reaction precursor and product forms of the HDV ribozyme. Our simulations find that an Mg2+ cation binds stably, by both inner- and outer-sphere contacts, to the electronegative catalytic pocket of the reaction precursor, in a position to potentially support chemistry. In contrast, protonation of the catalytically involved C75 in the precursor or artificial placement of this Mg2+ into the product structure result in its swift expulsion from the active site. These findings are consistent with a concerted reaction mechanism in which C75 and hydrated Mg2+ act as general base and acid, respectively. Monovalent cations bind to the active site and elsewhere assisted by structurally bridging long-residency water molecules, but are generally delocalized.
- MeSH
- hořčík chemie MeSH
- kationty dvojmocné chemie MeSH
- kationty jednomocné chemie MeSH
- konformace nukleové kyseliny MeSH
- molekulární modely MeSH
- molekulární sekvence - údaje MeSH
- RNA katalytická chemie MeSH
- sekvence nukleotidů MeSH
- sodík chemie MeSH
- vazebná místa MeSH
- virus hepatitidy delta enzymologie MeSH
- voda chemie MeSH
- vodíková vazba MeSH
- Publikační typ
- časopisecké články MeSH
- práce podpořená grantem MeSH
- Research Support, N.I.H., Extramural MeSH
- Názvy látek
- hořčík MeSH
- kationty dvojmocné MeSH
- kationty jednomocné MeSH
- RNA katalytická MeSH
- sodík MeSH
- voda MeSH
The hepatitis delta virus (HDV) ribozyme is a self-cleaving RNA enzyme involved in the replication of a human pathogen, the hepatitis delta virus. Recent crystal structures of the precursor and product of self-cleavage, together with detailed kinetic analyses, have led to hypotheses on the catalytic strategies employed by the HDV ribozyme. We report molecular dynamics (MD) simulations (approximately 120 ns total simulation time) to test the plausibility that specific conformational rearrangements are involved in catalysis. Site-specific self-cleavage requires cytidine in position 75 (C75). A precursor simulation with unprotonated C75 reveals a rather weak dynamic binding of C75 in the catalytic pocket with spontaneous, transient formation of a H-bond between U-1(O2') and C75(N3). This H-bond would be required for C75 to act as the general base. Upon protonation in the precursor, C75H+ has a tendency to move towards its product location and establish a firm H-bonding network within the catalytic pocket. However, a C75H+(N3)-G1(O5') H-bond, which would be expected if C75 acted as a general acid catalyst, is not observed on the present simulation timescale. The adjacent loop L3 is relatively dynamic and may serve as a flexible structural element, possibly gated by the closing U20.G25 base-pair, to facilitate a conformational switch induced by a protonated C75H+. L3 also controls the electrostatic environment of the catalytic core, which in turn may modulate C75 base strength and metal ion binding. We find that a distant RNA tertiary interaction involving a protonated cytidine (C41) becomes unstable when left unprotonated, leading to disruptive conformational rearrangements adjacent to the catalytic core. A Na ion temporarily compensates for the loss of the protonated hydrogen bond, which is strikingly consistent with the experimentally observed synergy between low pH and high Na+ concentrations in mediating residual self-cleavage of the HDV ribozyme in the absence of divalents.
- MeSH
- katalytická doména MeSH
- konformace nukleové kyseliny MeSH
- krystalografie rentgenová MeSH
- lidé MeSH
- molekulární modely MeSH
- prekurzory RNA chemie genetika metabolismus MeSH
- protony MeSH
- RNA katalytická chemie genetika metabolismus MeSH
- RNA virová chemie genetika metabolismus MeSH
- sekvence nukleotidů MeSH
- statická elektřina MeSH
- termodynamika MeSH
- virus hepatitidy delta enzymologie genetika MeSH
- vodíková vazba MeSH
- Check Tag
- lidé MeSH
- Publikační typ
- časopisecké články MeSH
- práce podpořená grantem MeSH
- Research Support, N.I.H., Extramural MeSH
- Research Support, U.S. Gov't, P.H.S. MeSH
- Názvy látek
- prekurzory RNA MeSH
- protony MeSH
- RNA katalytická MeSH
- RNA virová MeSH