An official website of the United States government
Genetic information is encoded in the DNA double helix, which, in its physiological milieu, is characterized by the iconical Watson-Crick nucleo-base pairing. Recent NMR relaxation experiments revealed the transient presence of an alternative, Hoogsteen (HG) base pairing pattern in naked DNA duplexes, and estimated its relative stability and lifetime. In contrast with DNA, such structures were not observed in RNA duplexes. Understanding HG base pairing is important because the underlying "breathing" motion between the two conformations can significantly modulate protein binding. However, a detailed mechanistic insight into the transition pathways and kinetics is still missing. We performed enhanced sampling simulation (with combined metadynamics and adaptive force-bias method) and Markov state modeling to obtain accurate free energy, kinetics, and the intermediates in the transition pathway between Watson-Crick and HG base pairs for both naked B-DNA and A-RNA duplexes. The Markov state model constructed from our unbiased MD simulation data revealed previously unknown complex extrahelical intermediates in the seemingly simple process of base flipping in B-DNA. Extending our calculation to A-RNA, for which HG base pairing is not observed experimentally, resulted in relatively unstable, single-hydrogen-bonded, distorted Hoogsteen-like bases. Unlike B-DNA, the transition pathway primarily involved base paired and intrahelical intermediates with transition timescales much longer than that of B-DNA. The seemingly obvious flip-over reaction coordinate (i.e., the glycosidic torsion angle) is unable to resolve the intermediates. Instead, a multidimensional picture involving backbone dihedral angles and distance between hydrogen bond donor and acceptor atoms is required to gain insight into the molecular mechanism.
more »
« less
Reaction Coordinate and Thermodynamics of Base Flipping in RNA
Base flipping is a key biophysical event involved in recognition of various ligands by ribonucleic acid (RNA) molecules. However, the mechanism of base flipping in RNA remains poorly understood, in part due to the lack of atomistic details on complex rearrangements in neighboring bases. In this work, we applied transition path sampling (TPS) methods to study base flipping in a double-stranded RNA (dsRNA) molecule that is known to interact with RNA-editing enzymes through this mechanism. We obtained an ensemble of 1000 transition trajectories to describe the base-flipping process. We used the likelihood maximization method to determine the refined reaction coordinate (RC) consisting of two collective variables (CVs), a distance and a dihedral angle between nucleotides that form stacking interactions with the flipping base. The free energy profile projected along the refined RC revealed three minima, two corresponding to the initial and final states and one for a metastable state. We suggest that the metastable state likely represents a wobbled conformation of nucleobases observed in NMR studies that is often characterized as the flipped state. The analyses of reactive trajectories further revealed that the base flipping is coupled to a global conformational change in a stem-loop of dsRNA.
more »
« less
- PAR ID:
- 10215014
- Date Published:
- Journal Name:
- Journal of Chemical Theory and Computation
- ISSN:
- 1549-9618
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
RNA oligonucleotides are crucial for a range of biological functions and in many biotechnological applications. Herein, we measured, for the first time, the conductance of individual double-stranded (ds)RNA molecules and compared it with the conductance of single DNA : RNA hybrids. The average conductance values are similar for both biomolecules, but the distribution of conductance values shows an order of magnitude higher variability for dsRNA, indicating higher molecular flexibility of dsRNA. Microsecond Molecular Dynamics simulations explain this difference and provide structural insights into the higher stability of DNA : RNA duplex with atomic level of detail. The rotations of 2′-OH groups of the ribose rings and the bases in RNA strands destabilize the duplex structure by weakening base stacking interactions, affecting charge transport, and making single-molecule conductance of dsRNA more variable (dynamic disorder). The results demonstrate that a powerful combination of state-of-the-art biomolecular electronics techniques and computational approaches can provide valuable insights into biomolecules’ biophysics with unprecedented spatial resolution.more » « less
-
Abstract RNA turnover pathways ensure appropriate gene expression levels by eliminating unwanted transcripts. Dis3-like 2 (Dis3L2) is a 3′–5′ exoribonuclease that plays a critical role in human development. Dis3L2 independently degrades structured substrates, including coding and noncoding 3′ uridylated RNAs. While the basis for Dis3L2’s substrate recognition has been well characterized, the mechanism of structured RNA degradation by this family of enzymes is unknown. We characterized the discrete steps of the degradation cycle by determining cryogenic electron microscopy structures representing snapshots along the RNA turnover pathway and measuring kinetic parameters for RNA processing. We discovered a dramatic conformational change that is triggered by double-stranded RNA (dsRNA), repositioning two cold shock domains by 70 Å. This movement exposes a trihelix linker region, which acts as a wedge to separate the two RNA strands. Furthermore, we show that the trihelix linker is critical for dsRNA, but not single-stranded RNA, degradation. These findings reveal the conformational plasticity of Dis3L2 and detail a mechanism of structured RNA degradation.more » « less
-
null (Ed.)The transfer of a β-hydrogen from a metal-alkyl group to ethylene is a fundamental organometallic transformation. Previously proposed mechanisms for this transformation involve either a two-step β-hydrogen elimination and migratory insertion sequence with a metal hydride intermediate or a one-step concerted pathway. Here, we report density functional theory (DFT) quasiclassical direct dynamics trajectories that reveal new dynamical mechanisms for the β-hydrogen transfer of [Cp*Rh III (Et)(ethylene)] + . Despite the DFT energy landscape showing a two-step mechanism with a Rh–H intermediate, quasiclassical trajectories commencing from the β-hydrogen elimination transition state revealed complete dynamical skipping of this intermediate. The skipping occurred either extremely fast (typically <100 femtoseconds (fs)) through a dynamically ballistic mechanism or slower through a dynamically unrelaxed mechanism. Consistent with trajectories begun at the transition state, all trajectories initiated at the Rh–H intermediate show continuation along the reaction coordinate. All of these trajectory outcomes are consistent with the Rh–H intermediate <1 kcal mol −1 stabilized relative to the β-hydrogen elimination and migratory insertion transition states. For Co, which on the energy landscape is a one-step concerted mechanism, trajectories showed extremely fast traversing of the transition-state zone (<50 fs), and this concerted mechanism is dynamically different than the Rh ballistic mechanism. In contrast to Rh, for Ir, in addition to dynamically ballistic and unrelaxed mechanisms, trajectories also stopped at the Ir–H intermediate. This is consistent with an Ir–H intermediate that is stabilized by ∼3 kcal mol −1 relative to the β-hydrogen elimination and migratory insertion transition states. Overall, comparison of Rh to Co and Ir provides understanding of the relationship between the energy surface shape and resulting dynamical mechanisms of an organometallic transformation.more » « less
-
Mycoviruses are widely distributed across fungi, including the yeasts of the Saccharomycotina subphylum. This manuscript reports the first double-stranded RNA (dsRNA) virus isolated from Pichia membranifaciens. This novel virus has been named Pichia membranifaciens virus L-A (PmV-L-A) and is a member of the Totiviridae. PmV-L-A is 4579 bp in length, with RNA secondary structures similar to the packaging, replication, and frameshift signals of totiviruses that infect Saccharomycotina yeasts. PmV-L-A was found to be part of a monophyletic group within the I-A totiviruses, implying a shared ancestry between mycoviruses isolated from the Pichiaceae and Saccharomycetaceae yeasts. Energy-minimized AlphaFold2 molecular models of the PmV-L-A Gag protein revealed structural conservation with the Gag protein of Saccharomyces cerevisiae virus L-A (ScV-L-A). The predicted tertiary structure of the PmV-L-A Pol and other homologs provided a possible mechanism for totivirus RNA replication due to structural similarities with the RNA-dependent RNA polymerases of mammalian dsRNA viruses. Insights into the structure, function, and evolution of totiviruses gained from yeasts are essential because of their emerging role in animal disease and their parallels with mammalian viruses.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1021/acs.jctc.0c01199
