Nejvíce citovaný článek - PubMed ID 36989143
Sensitivity of the RNA Structure to Ion Conditions as Probed by Molecular Dynamics Simulations of Common Canonical RNA Duplexes
Mixed double helices formed by RNA and DNA strands, commonly referred to as hybrid duplexes or hybrids, are essential in biological processes like transcription and reverse transcription. They are also important for their applications in CRISPR gene editing and nanotechnology. Yet, despite their significance, the hybrid duplexes have been seldom modeled by atomistic molecular dynamics methodology, and there is no benchmark study systematically assessing the force-field performance. Here, we present an extensive benchmark study of polypurine tract (PPT) and Dickerson-Drew dodecamer hybrid duplexes using contemporary and commonly utilized pairwise additive and polarizable nucleic acid force fields. Our findings indicate that none of the available force-field choices accurately reproduces all the characteristic structural details of the hybrid duplexes. The AMBER force fields are unable to populate the C3'-endo (north) pucker of the DNA strand and underestimate inclination. The CHARMM force field accurately describes the C3'-endo pucker and inclination but shows base pair instability. The polarizable force fields struggle with accurately reproducing the helical parameters. Some force-field combinations even demonstrate a discernible conflict between the RNA and DNA parameters. In this work, we offer a candid assessment of the force-field performance for mixed DNA/RNA duplexes. We provide guidance on selecting utilizable force-field combinations and also highlight potential pitfalls and best practices for obtaining optimal performance.
- MeSH
- DNA * chemie MeSH
- konformace nukleové kyseliny * MeSH
- párování bází MeSH
- RNA * chemie MeSH
- simulace molekulární dynamiky * MeSH
- Publikační typ
- časopisecké články MeSH
- Názvy látek
- DNA * MeSH
- RNA * MeSH
Nucleic acid double helices in their DNA, RNA, and DNA-RNA hybrid form play a fundamental role in biology and are main building blocks of artificial nanostructures, but how their properties depend on temperature remains poorly understood. Here, we report thermal dependence of dynamic bending persistence length, twist rigidity, stretch modulus, and twist-stretch coupling for DNA, RNA, and hybrid duplexes between 7°C and 47°C. The results are based on all-atom molecular dynamics simulations using different force field parameterizations. We first demonstrate that unrestrained molecular dynamics can reproduce experimentally known mechanical properties of the duplexes at room temperature. Beyond experimentally known features, we also infer the twist rigidity and twist-stretch coupling of the hybrid duplex. As for the temperature dependence, we found that increasing temperature softens all the duplexes with respect to bending, twisting, and stretching. The relative decrease of the stretch moduli is 0.003-0.004/°C, similar for all the duplex variants despite their very different stretching stiffness, whereas RNA twist stiffness decreases by 0.003/°C, and smaller values are found for the other elastic moduli. The twist-stretch couplings are nearly unaffected by temperature. The stretching, bending, and twisting stiffness all include an important entropic component. Relation of our results to the two-state model of DNA flexibility is discussed. Our work provides temperature-dependent elasticity of nucleic acid duplexes at the microsecond scale relevant for initial stages of protein binding.
RNA plays critical roles in the transmission and regulation of genetic information and is increasingly used in biomedical and biotechnological applications. Functional RNAs contain extended double-stranded regions, and the structure of double-stranded RNA (dsRNA) has been revealed at high resolution. However, the dependence of the properties of the RNA double helix on environmental effects, notably temperature, is still poorly understood. Here, we use single-molecule magnetic tweezer measurements to determine the dependence of the dsRNA twist on temperature. We find that dsRNA unwinds with increasing temperature, even more than DNA, with ΔTwRNA = -14.4 ± 0.7°/(°C·kbp), compared to ΔTwDNA = -11.0 ± 1.2°/(°C·kbp). All-atom molecular dynamics (MD) simulations using a range of nucleic acid force fields, ion parameters, and water models correctly predict that dsRNA unwinds with rising temperature but significantly underestimate the magnitude of the effect. These MD data, together with additional MD simulations involving DNA and DNA-RNA hybrid duplexes, reveal a linear correlation between the twist temperature decrease and the helical rise, in line with DNA but at variance with RNA experimental data. We speculate that this discrepancy might be caused by some unknown bias in the RNA force fields tested or by as yet undiscovered transient alternative structures in the RNA duplex. Our results provide a baseline to model more complex RNA assemblies and to test and develop new parametrizations for RNA simulations. They may also inspire physical models of the temperature-dependent dsRNA structure.
Molecular dynamics (MD) simulations represent an established tool to study RNA molecules. The outcome of MD studies depends, however, on the quality of the force field (ff). Here we suggest a correction for the widely used AMBER OL3 ff by adding a simple adjustment of the nonbonded parameters. The reparameterization of the Lennard-Jones potential for the -H8···O5'- and -H6···O5'- atom pairs addresses an intranucleotide steric clash occurring in the type 0 base-phosphate interaction (0BPh). The nonbonded fix (NBfix) modification of 0BPh interactions (NBfix0BPh modification) was tuned via a reweighting approach and subsequently tested using an extensive set of standard and enhanced sampling simulations of both unstructured and folded RNA motifs. The modification corrects minor but visible intranucleotide clash for the anti nucleobase conformation. We observed that structural ensembles of small RNA benchmark motifs simulated with the NBfix0BPh modification provide better agreement with experiments. No side effects of the modification were observed in standard simulations of larger structured RNA motifs. We suggest that the combination of OL3 RNA ff and NBfix0BPh modification is a viable option to improve RNA MD simulations.
- MeSH
- fosfáty * MeSH
- molekulární konformace MeSH
- nukleotidové motivy MeSH
- RNA * chemie MeSH
- simulace molekulární dynamiky MeSH
- Publikační typ
- časopisecké články MeSH
- Názvy látek
- fosfáty * MeSH
- RNA * MeSH
