The structure and deformability of double-stranded DNA and RNA depend on the sequence of bases, affecting biological processes and nanostructure design, but this dependence is incompletely understood. Here we present mechanical properties of DNA and RNA duplexes inferred from atomic-resolution, explicit-solvent molecular dynamics (MD) simulations of 107 DNA and 107 RNA oligomers containing all hexanucleotide sequences. In addition to the level of rigid bases, minor and major grooves, we probe the length and sequence dependence of global material constants such as persistence lengths, stretching and twisting rigidities. We propose a simple model to predict sequence-dependent shape and nonlocal, harmonic stiffness for an arbitrary sequence, validate it on an independent set of MD simulations for DNA and RNA duplexes containing all pentamers, and demonstrate its utility in various applications. The large amount of the simulated data enabled us to study rare events, such as base-pair opening, or flips of the A-RNA sugar pucker into the B domain and the related dynamics of the 2'-OH group. Together, this work provides a comprehensive sequence-specific description of DNA and RNA duplex mechanics, forming a baseline for further research and allowing for a broad range of applications.

• MeSH
• DNA * chemistry   MeSH
• Nucleic Acid Conformation MeSH
• RNA * chemistry   MeSH
• Base Sequence MeSH
• Molecular Dynamics Simulation MeSH
• Publication type
• Journal Article MeSH
• Names of Substances
• DNA * MeSH
• RNA * MeSH

Molecular dynamics (MD) simulations are an important and well-established tool for investigating RNA structural dynamics, but their accuracy relies heavily on the quality of the employed force field (ff). In this work, we present a comprehensive evaluation of widely used pair-additive and polarizable RNA ffs using the challenging UUCG tetraloop (TL) benchmark system. Extensive standard MD simulations, initiated from the NMR structure of the 14-mer UUCG TL, revealed that most ffs did not maintain the native state, instead favoring alternative loop conformations. Notably, three very recent variants of pair-additive ffs, OL3CP-gHBfix21, DES-Amber, and OL3R2.7, successfully preserved the native structure over a 10 × 20 μs time scale. To further assess these ffs, we performed enhanced sampling folding simulations of the shorter 8-mer UUCG TL, starting from the single-stranded conformation. Estimated folding free energies (ΔG°fold) varied significantly among these three ffs, with values of 0.0 ± 0.6, 2.4 ± 0.8, and 7.4 ± 0.2 kcal/mol for OL3CP-gHBfix21, DES-Amber, and OL3R2.7, respectively. The ΔG°fold value predicted by the OL3CP-gHBfix21 ff was closest to experimental estimates, ranging from -1.6 to -0.7 kcal/mol. In contrast, the higher ΔG°fold values obtained using DES-Amber and OL3R2.7 were unexpected, suggesting that key interactions are inaccurately described in the folded, unfolded, or misfolded ensembles. These discrepancies led us to further test DES-Amber and OL3R2.7 ffs on additional RNA and DNA systems, where further performance issues were observed. Our results emphasize the complexity of accurately modeling RNA dynamics and suggest that creating an RNA ff capable of reliably performing across a wide range of RNA systems remains extremely challenging. In conclusion, our study provides valuable insights into the capabilities of current RNA ffs and highlights key areas for future ff development.

• MeSH
• Nucleic Acid Conformation MeSH
• RNA * chemistry   MeSH
• Molecular Dynamics Simulation * MeSH
• Thermodynamics MeSH
• Publication type
• Journal Article MeSH
• Names of Substances
• RNA * MeSH

Lipid-mediated delivery of active pharmaceutical ingredients (API) opened new possibilities in advanced therapies. By encapsulating an API into a lipid nanocarrier (LNC), one can safely deliver APIs not soluble in water, those with otherwise strong adverse effects, or very fragile ones such as nucleic acids. However, for the rational design of LNCs, a detailed understanding of the composition-structure-function relationships is missing. This review presents currently available computational methods for LNC investigation, screening, and design. The state-of-the-art physics-based approaches are described, with the focus on molecular dynamics simulations in all-atom and coarse-grained resolution. Their strengths and weaknesses are discussed, highlighting the aspects necessary for obtaining reliable results in the simulations. Furthermore, a machine learning, i.e., data-based learning, approach to the design of lipid-mediated API delivery is introduced. The data produced by the experimental and theoretical approaches provide valuable insights. Processing these data can help optimize the design of LNCs for better performance. In the final section of this Review, state-of-the-art of computer simulations of LNCs are reviewed, specifically addressing the compatibility of experimental and computational insights.

• Keywords
• ionizable lipid, lipid nanocarrier, lipid nanoparticle, liposome, molecular simulation, vesicle,
• MeSH
• Pharmaceutical Preparations chemistry   administration & dosage   MeSH
• Drug Delivery Systems * methods   MeSH
• Humans MeSH
• Lipids * chemistry   MeSH
• Nanoparticles chemistry   MeSH
• Drug Carriers * chemistry   MeSH
• Computer Simulation MeSH
• Molecular Dynamics Simulation MeSH
• Machine Learning MeSH
• Check Tag
• Humans MeSH
• Publication type
• Journal Article MeSH
• Review MeSH
• Names of Substances
• Pharmaceutical Preparations MeSH
• Lipids * MeSH
• Drug Carriers * MeSH

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 * chemistry   MeSH
• Nucleic Acid Conformation * MeSH
• Base Pairing MeSH
• RNA * chemistry   MeSH
• Molecular Dynamics Simulation * MeSH
• Publication type
• Journal Article MeSH
• Names of Substances
• 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.

• MeSH
• DNA * chemistry   MeSH
• Nucleic Acid Conformation MeSH
• Elasticity MeSH
• RNA * chemistry   MeSH
• Temperature MeSH
• Publication type
• Journal Article MeSH
• Research Support, Non-U.S. Gov't MeSH
• Names of Substances
• DNA * MeSH
• RNA * MeSH

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.

• MeSH
• DNA chemistry   MeSH
• RNA, Double-Stranded * MeSH
• Nucleic Acid Conformation MeSH
• Magnetic Phenomena MeSH
• RNA chemistry   MeSH
• Molecular Dynamics Simulation * MeSH
• Temperature MeSH
• Publication type
• Journal Article MeSH
• Research Support, Non-U.S. Gov't MeSH
• Names of Substances
• DNA MeSH
• RNA, Double-Stranded * MeSH
• RNA MeSH

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
• Phosphates * MeSH
• Molecular Conformation MeSH
• Nucleotide Motifs MeSH
• RNA * chemistry   MeSH
• Molecular Dynamics Simulation MeSH
• Publication type
• Journal Article MeSH
• Names of Substances
• Phosphates * MeSH
• RNA * MeSH

Recognition of single-stranded RNA (ssRNA) by RNA recognition motif (RRM) domains is an important class of protein-RNA interactions. Many such complexes were characterized using nuclear magnetic resonance (NMR) and/or X-ray crystallography techniques, revealing ensemble-averaged pictures of the bound states. However, it is becoming widely accepted that better understanding of protein-RNA interactions would be obtained from ensemble descriptions. Indeed, earlier molecular dynamics simulations of bound states indicated visible dynamics at the RNA-RRM interfaces. Here, we report the first atomistic simulation study of spontaneous binding of short RNA sequences to RRM domains of HuR and SRSF1 proteins. Using a millisecond-scale aggregate ensemble of unbiased simulations, we were able to observe a few dozen binding events. HuR RRM3 utilizes a pre-binding state to navigate the RNA sequence to its partially disordered bound state and then to dynamically scan its different binding registers. SRSF1 RRM2 binding is more straightforward but still multiple-pathway. The present study necessitated development of a goal-specific force field modification, scaling down the intramolecular van der Waals interactions of the RNA which also improves description of the RNA-RRM bound state. Our study opens up a new avenue for large-scale atomistic investigations of binding landscapes of protein-RNA complexes, and future perspectives of such research are discussed.

• MeSH
• ELAV-Like Protein 1 metabolism   MeSH
• RNA Recognition Motif genetics   MeSH
• RNA-Binding Proteins * metabolism   MeSH
• RNA * chemistry   MeSH
• RNA Recognition Motif Proteins metabolism   MeSH
• Protein Binding MeSH
• Binding Sites MeSH
• Publication type
• Journal Article MeSH
• Research Support, Non-U.S. Gov't MeSH
• Names of Substances
• ELAV-Like Protein 1 MeSH
• RNA-Binding Proteins * MeSH
• RNA * MeSH
• RNA Recognition Motif Proteins MeSH

Kink-turns are highly bent internal loop motifs commonly found in the ribosome and other RNA complexes. They frequently act as binding sites for proteins and mediate tertiary interactions in larger RNA structures. Kink-turns have been a topic of intense research, but their elastic properties in the folded state are still poorly understood. Here we use extensive all-atom molecular dynamics simulations to parameterize a model of kink-turn in which the two flanking helical stems are represented by effective rigid bodies. Time series of the full set of six interhelical coordinates enable us to extract minimum energy shapes and harmonic stiffness constants for kink-turns from different RNA functional classes. The analysis suggests that kink-turns exhibit isotropic bending stiffness but are highly anisotropic with respect to lateral displacement of the stems. The most flexible lateral displacement mode is perpendicular to the plane of the static bend. These results may help understand the structural adaptation and mechanical signal transmission by kink-turns in complex natural and artificial RNA structures.

• MeSH
• Nucleic Acid Conformation MeSH
• Ribosomes metabolism   MeSH
• RNA * chemistry   MeSH
• Molecular Dynamics Simulation * MeSH
• Binding Sites MeSH
• Publication type
• Journal Article MeSH
• Research Support, Non-U.S. Gov't MeSH
• Names of Substances
• RNA * MeSH

The widespread Mn2+-sensing yybP-ykoY riboswitch controls the expression of bacterial Mn2+ homeostasis genes. Here, we first determine the crystal structure of the ligand-bound yybP-ykoY riboswitch aptamer from Xanthomonas oryzae at 2.96 Å resolution, revealing two conformations with docked four-way junction (4WJ) and incompletely coordinated metal ions. In >100 µs of MD simulations, we observe that loss of divalents from the core triggers local structural perturbations in the adjacent docking interface, laying the foundation for signal transduction to the regulatory switch helix. Using single-molecule FRET, we unveil a previously unobserved extended 4WJ conformation that samples transient docked states in the presence of Mg2+. Only upon adding sub-millimolar Mn2+, however, can the 4WJ dock stably, a feature lost upon mutation of an adenosine contacting Mn2+ in the core. These observations illuminate how subtly differing ligand preferences of competing metal ions become amplified by the coupling of local with global RNA dynamics.

• MeSH
• RNA, Bacterial chemistry   genetics   metabolism   MeSH
• Escherichia coli genetics   MeSH
• Magnesium metabolism   MeSH
• Nucleic Acid Conformation MeSH
• Crystallography, X-Ray MeSH
• Lactococcus lactis genetics   metabolism   MeSH
• Ligands MeSH
• Manganese metabolism   MeSH
• Molecular Conformation MeSH
• Models, Molecular MeSH
• Mutation MeSH
• Gene Expression Regulation, Bacterial MeSH
• Riboswitch physiology   MeSH
• Signal Transduction * MeSH
• Molecular Dynamics Simulation MeSH
• Binding Sites MeSH
• Xanthomonas metabolism   MeSH
• Publication type
• Journal Article MeSH
• Research Support, Non-U.S. Gov't MeSH
• Research Support, N.I.H., Extramural MeSH
• Research Support, U.S. Gov't, Non-P.H.S. MeSH
• Names of Substances
• RNA, Bacterial MeSH
• Magnesium MeSH
• Ligands MeSH
• Manganese MeSH
• Riboswitch MeSH

Molecular dynamics (MD) simulations became a leading tool for investigation of structural dynamics of nucleic acids. Despite recent efforts to improve the empirical potentials (force fields, ffs), RNA ffs have persisting deficiencies, which hamper their utilization in quantitatively accurate simulations. Previous studies have shown that at least two salient problems contribute to difficulties in the description of free-energy landscapes of small RNA motifs: (i) excessive stabilization of the unfolded single-stranded RNA ensemble by intramolecular base-phosphate and sugar-phosphate interactions and (ii) destabilization of the native folded state by underestimation of stability of base pairing. Here, we introduce a general ff term (gHBfix) that can selectively fine-tune nonbonding interaction terms in RNA ffs, in particular, the H bonds. The gHBfix potential affects the pairwise interactions between all possible pairs of the specific atom types, while all other interactions remain intact; i.e., it is not a structure-based model. In order to probe the ability of the gHBfix potential to refine the ff nonbonded terms, we performed an extensive set of folding simulations of RNA tetranucleotides and tetraloops. On the basis of these data, we propose particular gHBfix parameters to modify the AMBER RNA ff. The suggested parametrization significantly improves the agreement between experimental data and the simulation conformational ensembles, although our current ff version still remains far from being flawless. While attempts to tune the RNA ffs by conventional reparametrizations of dihedral potentials or nonbonded terms can lead to major undesired side effects, as we demonstrate for some recently published ffs, gHBfix has a clear promising potential to improve the ff performance while avoiding introduction of major new imbalances.

• MeSH
• RNA chemistry   MeSH
• Molecular Dynamics Simulation * MeSH
• Hydrogen Bonding MeSH
• Publication type
• Journal Article MeSH
• Names of Substances
• RNA MeSH

We have carried out an extended set of standard and enhanced-sampling MD simulations (for a cumulative simulation time of 620 μs) with the aim to study folding landscapes of the rGGGUUAGGG and rGGGAGGG parallel G-hairpins (PH) with propeller loop. We identify folding and unfolding pathways of the PH, which is bridged with the unfolded state via an ensemble of cross-like structures (CS) possessing mutually tilted or perpendicular G-strands interacting via guanine-guanine H-bonding. The oligonucleotides reach the PH conformation from the unfolded state via a conformational diffusion through the folding landscape, i.e. as a series of rearrangements of the H-bond interactions starting from compacted anti-parallel hairpin-like structures. Although isolated PHs do not appear to be thermodynamically stable we suggest that CS and PH-types of structures are sufficiently populated during RNA guanine quadruplex (GQ) folding within the context of complete GQ-forming sequences. These structures may participate in compact coil-like ensembles that involve all four G-strands and already some bound ions. Such ensembles can then rearrange into the fully folded parallel GQs via conformational diffusion. We propose that the basic atomistic folding mechanism of propeller loops suggested in this work may be common for their formation in RNA and DNA GQs.

• MeSH
• G-Quadruplexes * MeSH
• Guanine chemistry   metabolism   MeSH
• Kinetics MeSH
• RNA chemistry   metabolism   MeSH
• RNA Folding * MeSH
• Base Sequence MeSH
• Molecular Dynamics Simulation MeSH
• Thermodynamics MeSH
• Hydrogen Bonding MeSH
• Publication type
• Journal Article MeSH
• Research Support, Non-U.S. Gov't MeSH
• Names of Substances
• Guanine MeSH
• RNA MeSH

With both catalytic and genetic functions, ribonucleic acid (RNA) is perhaps the most pluripotent chemical species in molecular biology, and its functions are intimately linked to its structure and dynamics. Computer simulations, and in particular atomistic molecular dynamics (MD), allow structural dynamics of biomolecular systems to be investigated with unprecedented temporal and spatial resolution. We here provide a comprehensive overview of the fast-developing field of MD simulations of RNA molecules. We begin with an in-depth, evaluatory coverage of the most fundamental methodological challenges that set the basis for the future development of the field, in particular, the current developments and inherent physical limitations of the atomistic force fields and the recent advances in a broad spectrum of enhanced sampling methods. We also survey the closely related field of coarse-grained modeling of RNA systems. After dealing with the methodological aspects, we provide an exhaustive overview of the available RNA simulation literature, ranging from studies of the smallest RNA oligonucleotides to investigations of the entire ribosome. Our review encompasses tetranucleotides, tetraloops, a number of small RNA motifs, A-helix RNA, kissing-loop complexes, the TAR RNA element, the decoding center and other important regions of the ribosome, as well as assorted others systems. Extended sections are devoted to RNA-ion interactions, ribozymes, riboswitches, and protein/RNA complexes. Our overview is written for as broad of an audience as possible, aiming to provide a much-needed interdisciplinary bridge between computation and experiment, together with a perspective on the future of the field.

• MeSH
• DNA chemistry   MeSH
• Catalysis MeSH
• Nucleic Acid Conformation * MeSH
• Computer Simulation MeSH
• RNA chemistry   MeSH
• Molecular Dynamics Simulation * MeSH
• Publication type
• Journal Article MeSH
• Research Support, Non-U.S. Gov't MeSH
• Review MeSH
• Research Support, N.I.H., Extramural MeSH
• Names of Substances
• DNA MeSH
• RNA MeSH

The cyclooxygenase-2 is a pro-inflammatory and cancer marker, whose mRNA stability and translation is regulated by the CUG-binding protein 2 interacting with AU-rich sequences in the 3' untranslated region. Here, we present the solution NMR structure of CUG-binding protein 2 RRM3 in complex with 5'-UUUAA-3' originating from the COX-2 3'-UTR. We show that RRM3 uses the same binding surface and protein moieties to interact with AU- and UG-rich RNA motifs, binding with low and high affinity, respectively. Using NMR spectroscopy, isothermal titration calorimetry and molecular dynamics simulations, we demonstrate that distinct sub-states characterized by different aromatic side-chain conformations at the RNA-binding surface allow for high- or low-affinity binding with functional implications. This study highlights a mechanism for RNA discrimination possibly common to multiple RRMs as several prominent members display a similar rearrangement of aromatic residues upon binding their targets.The RNA Recognition Motif (RRM) is the most ubiquitous RNA binding domain. Here the authors combined NMR and molecular dynamics simulations and show that the RRM RNA binding surface exists in different states and that a conformational switch of aromatic side-chains fine-tunes sequence specific binding affinities.

• MeSH
• 3' Untranslated Regions MeSH
• Amino Acid Motifs MeSH
• CELF Proteins chemistry   genetics   metabolism   MeSH
• Cyclooxygenase 2 genetics   MeSH
• Phenylalanine chemistry   metabolism   MeSH
• Protein Conformation MeSH
• Magnetic Resonance Spectroscopy MeSH
• RNA, Messenger chemistry   metabolism   MeSH
• Nerve Tissue Proteins chemistry   genetics   metabolism   MeSH
• Molecular Dynamics Simulation MeSH
• Amino Acid Substitution MeSH
• AU Rich Elements MeSH
• Binding Sites MeSH
• Publication type
• Journal Article MeSH
• Research Support, Non-U.S. Gov't MeSH
• Names of Substances
• 3' Untranslated Regions MeSH
• CELF Proteins MeSH
• CELF2 protein, human MeSH Browser
• Cyclooxygenase 2 MeSH
• Phenylalanine MeSH
• RNA, Messenger MeSH
• Nerve Tissue Proteins MeSH
• PTGS2 protein, human MeSH Browser

The computer-aided folding of biomolecules, particularly RNAs, is one of the most difficult challenges in computational structural biology. RNA tetraloops are fundamental RNA motifs playing key roles in RNA folding and RNA-RNA and RNA-protein interactions. Although state-of-the-art Molecular Dynamics (MD) force fields correctly describe the native state of these tetraloops as a stable free-energy basin on the microsecond time scale, enhanced sampling techniques reveal that the native state is not the global free energy minimum, suggesting yet unidentified significant imbalances in the force fields. Here, we tested our ability to fold the RNA tetraloops in various force fields and simulation settings. We employed three different enhanced sampling techniques, namely, temperature replica exchange MD (T-REMD), replica exchange with solute tempering (REST2), and well-tempered metadynamics (WT-MetaD). We aimed to separate problems caused by limited sampling from those due to force-field inaccuracies. We found that none of the contemporary force fields is able to correctly describe folding of the 5'-GAGA-3' tetraloop over a range of simulation conditions. We thus aimed to identify which terms of the force field are responsible for this poor description of TL folding. We showed that at least two different imbalances contribute to this behavior, namely, overstabilization of base-phosphate and/or sugar-phosphate interactions and underestimated stability of the hydrogen bonding interaction in base pairing. The first artifact stabilizes the unfolded ensemble, while the second one destabilizes the folded state. The former problem might be partially alleviated by reparametrization of the van der Waals parameters of the phosphate oxygens suggested by Case et al., while in order to overcome the latter effect we suggest local potentials to better capture hydrogen bonding interactions.

• MeSH
• Nucleic Acid Conformation MeSH
• RNA chemistry   metabolism   MeSH
• RNA Folding MeSH
• Molecular Dynamics Simulation * MeSH
• RNA Stability MeSH
• Static Electricity MeSH
• Temperature MeSH
• Hydrogen Bonding MeSH
• Publication type
• Journal Article MeSH
• Names of Substances
• RNA MeSH

RNA recognition motif (RRM) proteins represent an abundant class of proteins playing key roles in RNA biology. We present a joint atomistic molecular dynamics (MD) and experimental study of two RRM-containing proteins bound with their single-stranded target RNAs, namely the Fox-1 and SRSF1 complexes. The simulations are used in conjunction with NMR spectroscopy to interpret and expand the available structural data. We accumulate more than 50 μs of simulations and show that the MD method is robust enough to reliably describe the structural dynamics of the RRM-RNA complexes. The simulations predict unanticipated specific participation of Arg142 at the protein-RNA interface of the SRFS1 complex, which is subsequently confirmed by NMR and ITC measurements. Several segments of the protein-RNA interface may involve competition between dynamical local substates rather than firmly formed interactions, which is indirectly consistent with the primary NMR data. We demonstrate that the simulations can be used to interpret the NMR atomistic models and can provide qualified predictions. Finally, we propose a protocol for 'MD-adapted structure ensemble' as a way to integrate the simulation predictions and expand upon the deposited NMR structures. Unbiased μs-scale atomistic MD could become a technique routinely complementing the NMR measurements of protein-RNA complexes.

• MeSH
• Protein Conformation MeSH
• Humans MeSH
• Magnetic Resonance Spectroscopy MeSH
• Models, Molecular MeSH
• RNA Recognition Motif genetics   MeSH
• Multiprotein Complexes chemistry   genetics   MeSH
• RNA chemistry   genetics   MeSH
• Amino Acid Sequence genetics   MeSH
• Serine-Arginine Splicing Factors chemistry   genetics   MeSH
• RNA Splicing Factors chemistry   genetics   MeSH
• Molecular Dynamics Simulation MeSH
• Binding Sites MeSH
• Check Tag
• Humans MeSH
• Publication type
• Journal Article MeSH
• Research Support, Non-U.S. Gov't MeSH
• Names of Substances
• Multiprotein Complexes MeSH
• RBFOX1 protein, human MeSH Browser
• RNA MeSH
• Serine-Arginine Splicing Factors MeSH
• RNA Splicing Factors MeSH
• SRSF1 protein, human MeSH Browser

In numerous Gram-positive bacteria, the glmS ribozyme or catalytic riboswitch regulates the expression of glucosamine-6-phosphate (GlcN6P) synthase via site-specific cleavage of its sugar-phosphate backbone in response to GlcN6P ligand binding. Biochemical data have suggested a crucial catalytic role for an active site guanine (G40 in Thermoanaerobacter tengcongensis, G33 in Bacillus anthracis). We used hybrid quantum chemical/molecular mechanical (QM/MM) calculations to probe the mechanism where G40 is deprotonated and acts as a general base. The calculations suggest that the deprotonated guanine G40(-) is sufficiently reactive to overcome the thermodynamic penalty arising from its rare protonation state, and thus is able to activate the A-1(2'-OH) group toward nucleophilic attack on the adjacent backbone. Furthermore, deprotonation of A-1(2'-OH) and nucleophilic attack are predicted to occur as separate steps, where activation of A-1(2'-OH) precedes nucleophilic attack. Conversely, the transition state associated with the rate-determining step corresponds to concurrent nucleophilic attack and protonation of the G1(O5') leaving group by the ammonium moiety of the GlcN6P cofactor. Overall, our calculations help to explain the crucial roles of G40 (as a general base) and GlcN6P (as a general acid) during glmS ribozyme self-cleavage. In addition, we show that the QM/MM description of the glmS ribozyme self-cleavage reaction is significantly more sensitive to the size of the QM region and the quality of the QM-MM coupling than that of other small ribozymes.

• Keywords
• QM/MM, RNA catalysis, glmS, riboswitch, ribozyme,
• MeSH
• Catalysis MeSH
• Riboswitch MeSH
• RNA, Catalytic chemistry   MeSH
• Models, Theoretical MeSH
• Thermodynamics MeSH
• Publication type
• Journal Article MeSH
• Research Support, Non-U.S. Gov't MeSH
• Research Support, N.I.H., Extramural MeSH
• Names of Substances
• Riboswitch MeSH
• RNA, Catalytic MeSH

The hepatitis delta virus (HDV) ribozyme is a catalytic RNA motif embedded in the human pathogenic HDV RNA. It catalyzes self-cleavage of its sugar-phosphate backbone with direct participation of the active site cytosine C75. Biochemical and structural data support a general acid role of C75. Here, we used hybrid quantum mechanical/molecular mechanical (QM/MM) calculations to probe the reaction mechanism and changes in Gibbs energy along the ribozyme's reaction pathway with an N3-protonated C75H(+) in the active site, which acts as the general acid, and a partially hydrated Mg(2+) ion with one deprotonated, inner-shell coordinated water molecule that acts as the general base. We followed eight reaction paths with a distinct position and coordination of the catalytically important active site Mg(2+) ion. For six of them, we observed feasible activation barriers ranging from 14.2 to 21.9 kcal mol(-1), indicating that the specific position of the Mg(2+) ion in the active site is predicted to strongly affect the kinetics of self-cleavage. The deprotonation of the U-1(2'-OH) nucleophile and the nucleophilic attack of the resulting U-1(2'-O(-)) on the scissile phosphodiester are found to be separate steps, as deprotonation precedes the nucleophilic attack. This sequential mechanism of the HDV ribozyme differs from the concerted nucleophilic activation and attack suggested for the hairpin ribozyme. We estimate the pKa of the U-1(2'-OH) group to range from 8.8 to 11.2, suggesting that it is lowered by several units from that of a free ribose, comparable to and most likely smaller than the pKa of the solvated active site Mg(2+) ion. Our results thus support the notion that the structure of the HDV ribozyme, and particularly the positioning of the active site Mg(2+) ion, facilitate deprotonation and activation of the 2'-OH nucleophile.

• MeSH
• Hepatitis D virology   MeSH
• Magnesium chemistry   MeSH
• Catalytic Domain MeSH
• Nucleic Acid Conformation MeSH
• Crystallography, X-Ray MeSH
• Quantum Theory MeSH
• Humans MeSH
• Models, Molecular MeSH
• RNA, Catalytic chemistry   MeSH
• RNA, Viral chemistry   MeSH
• Thermodynamics MeSH
• Hepatitis Delta Virus chemistry   enzymology   MeSH
• Check Tag
• Humans MeSH
• Publication type
• Journal Article MeSH
• Research Support, Non-U.S. Gov't MeSH
• Research Support, N.I.H., Extramural MeSH
• Names of Substances
• hairpin ribozyme MeSH Browser
• Magnesium MeSH
• RNA, Catalytic MeSH
• RNA, Viral MeSH

The hepatitis delta virus (HDV) ribozyme is a member of the class of small, self-cleaving catalytic RNAs found in a wide range of genomes from HDV to human. Both pre- and post-catalysis (precursor and product) crystal structures of the cis-acting genomic HDV ribozyme have been determined. These structures, together with extensive solution probing, have suggested that a significant conformational change accompanies catalysis. A recent crystal structure of a trans-acting precursor, obtained at low pH and by molecular replacement from the previous product conformation, conforms to the product, raising the possibility that it represents an activated conformer past the conformational change. Here, using fluorescence resonance energy transfer (FRET), we discovered that cleavage of this ribozyme at physiological pH is accompanied by a structural lengthening in magnitude comparable to previous trans-acting HDV ribozymes. Conformational heterogeneity observed by FRET in solution appears to have been removed upon crystallization. Analysis of a total of 1.8 µsec of molecular dynamics (MD) simulations showed that the crystallographically unresolved cleavage site conformation is likely correctly modeled after the hammerhead ribozyme, but that crystal contacts and the removal of several 2'-oxygens near the scissile phosphate compromise catalytic in-line fitness. A cis-acting version of the ribozyme exhibits a more dynamic active site, while a G-1 residue upstream of the scissile phosphate favors poor fitness, allowing us to rationalize corresponding changes in catalytic activity. Based on these data, we propose that the available crystal structures of the HDV ribozyme represent intermediates on an overall rugged RNA folding free-energy landscape.

• Keywords
• conformational change, molecular dynamics simulation, small ribozyme, steady-state FRET, time-resolved FRET,
• MeSH
• Catalytic Domain MeSH
• Catalysis MeSH
• Kinetics MeSH
• Nucleic Acid Conformation MeSH
• Crystallography, X-Ray MeSH
• Models, Molecular MeSH
• Fluorescence Resonance Energy Transfer methods   MeSH
• RNA, Catalytic chemistry   MeSH
• RNA, Small Nuclear chemistry   metabolism   MeSH
• RNA, Viral chemistry   MeSH
• Molecular Dynamics Simulation MeSH
• RNA Cleavage MeSH
• Hepatitis Delta Virus enzymology   genetics   MeSH
• Publication type
• Journal Article MeSH
• Research Support, Non-U.S. Gov't MeSH
• Research Support, N.I.H., Extramural MeSH
• Research Support, U.S. Gov't, Non-P.H.S. MeSH
• Names of Substances
• hammerhead ribozyme MeSH Browser
• RNA, Catalytic MeSH
• RNA, Small Nuclear MeSH
• RNA, Viral MeSH
• U1 small nuclear RNA MeSH Browser

The sarcin-ricin RNA motif (SR motif) is one of the most prominent recurrent RNA building blocks that occurs in many different RNA contexts and folds autonomously, that is, in a context-independent manner. In this study, we combined bioinformatics analysis with explicit-solvent molecular dynamics (MD) simulations to better understand the relation between the RNA sequence and the evolutionary patterns of the SR motif. A SHAPE probing experiment was also performed to confirm the fidelity of the MD simulations. We identified 57 instances of the SR motif in a nonredundant subset of the RNA X-ray structure database and analyzed their base pairing, base-phosphate, and backbone-backbone interactions. We extracted sequences aligned to these instances from large rRNA alignments to determine the frequency of occurrence for different sequence variants. We then used a simple scoring scheme based on isostericity to suggest 10 sequence variants with a highly variable expected degree of compatibility with the SR motif 3D structure. We carried out MD simulations of SR motifs with these base substitutions. Nonisosteric base substitutions led to unstable structures, but so did isosteric substitutions which were unable to make key base-phosphate interactions. The MD technique explains why some potentially isosteric SR motifs are not realized during evolution. We also found that the inability to form stable cWW geometry is an important factor in the case of the first base pair of the flexible region of the SR motif. A comparison of structural, bioinformatics, SHAPE probing, and MD simulation data reveals that explicit solvent MD simulations neatly reflect the viability of different sequence variants of the SR motif. Thus, MD simulations can efficiently complement bioinformatics tools in studies of conservation patterns of RNA motifs and provide atomistic insight into the role of their different signature interactions.

• MeSH
• Nucleic Acid Conformation MeSH
• Nucleotide Motifs MeSH
• Base Pairing MeSH
• RNA, Ribosomal chemistry   metabolism   MeSH
• RNA chemistry   metabolism   MeSH
• Solvents chemistry   MeSH
• Molecular Dynamics Simulation MeSH
• Hydrogen Bonding MeSH
• Computational Biology MeSH
• Publication type
• Journal Article MeSH
• Research Support, Non-U.S. Gov't MeSH
• Research Support, N.I.H., Extramural MeSH
• Names of Substances
• RNA, Ribosomal MeSH
• RNA MeSH
• Solvents MeSH

Explicit solvent molecular dynamics simulations have been used to complement preceding experimental and computational studies of folding of guanine quadruplexes (G-DNA). We initiate early stages of unfolding of several G-DNAs by simulating them under no-salt conditions and then try to fold them back using standard excess salt simulations. There is a significant difference between G-DNAs with all-anti parallel stranded stems and those with stems containing mixtures of syn and anti guanosines. The most natural rearrangement for all-anti stems is a vertical mutual slippage of the strands. This leads to stems with reduced numbers of tetrads during unfolding and a reduction of strand slippage during refolding. The presence of syn nucleotides prevents mutual strand slippage; therefore, the antiparallel and hybrid quadruplexes initiate unfolding via separation of the individual strands. The simulations confirm the capability of G-DNA molecules to adopt numerous stable locally and globally misfolded structures. The key point for a proper individual folding attempt appears to be correct prior distribution of syn and anti nucleotides in all four G-strands. The results suggest that at the level of individual molecules, G-DNA folding is an extremely multi-pathway process that is slowed by numerous misfolding arrangements stabilized on highly variable timescales.

• MeSH
• DNA chemistry   MeSH
• G-Quadruplexes * MeSH
• DNA, Single-Stranded chemistry   MeSH
• Humans MeSH
• Molecular Dynamics Simulation * MeSH
• Telomere chemistry   MeSH
• Check Tag
• Humans MeSH
• Publication type
• Journal Article MeSH
• Research Support, Non-U.S. Gov't MeSH
• Names of Substances
• DNA MeSH
• DNA, Single-Stranded MeSH

We provide theoretical predictions of the intrinsic stability of different arrangements of guanine quadruplex (G-DNA) stems. Most computational studies of nucleic acids have applied Molecular Mechanics (MM) approaches using simple pairwise-additive force fields. The principle limitation of such calculations is the highly approximate nature of the force fields. In this study, we for the first time apply accurate QM computations (DFT-D3 with large atomic orbital basis sets) to essentially complete DNA building blocks, seven different folds of the cation-stabilized two-quartet G-DNA stem, each having more than 250 atoms. The solvent effects are approximated by COSMO continuum solvent. We reveal sizable differences between MM and QM descriptions of relative energies of different G-DNA stems, which apparently reflect approximations of the DNA force field. Using the QM energy data, we propose correction to earlier free energy estimates of relative stabilities of different parallel, hybrid, and antiparallel G-stem folds based on classical simulations. The new energy ranking visibly improves the agreement between theory and experiment. We predict the 5'-anti-anti-3' GpG dinucleotide step to be the most stable one, closely followed by the 5'-syn-anti-3' step. The results are in good agreement with known experimental structures of 2-, 3-, and 4-quartet G-DNA stems. Besides providing specific results for G-DNA, our study highlights basic limitations of force field modeling of nucleic acids. Although QM computations have their own limitations, mainly the lack of conformational sampling and the approximate description of the solvent, they can substantially improve the quality of calculations currently relying exclusively on force fields.

• MeSH
• DNA chemistry   MeSH
• G-Quadruplexes * MeSH
• Guanine chemistry   MeSH
• Quantum Theory * MeSH
• Models, Molecular MeSH
• Publication type
• Journal Article MeSH
• Research Support, Non-U.S. Gov't MeSH
• Research Support, N.I.H., Extramural MeSH
• Research Support, U.S. Gov't, Non-P.H.S. MeSH
• Names of Substances
• DNA MeSH
• Guanine MeSH

We present a refinement of the backbone torsion parameters ε and ζ of the Cornell et al. AMBER force field for DNA simulations. The new parameters, denoted as εζOL1, were derived from quantum-mechanical calculations with inclusion of conformation-dependent solvation effects according to the recently reported methodology (J. Chem. Theory Comput. 2012, 7(9), 2886-2902). The performance of the refined parameters was analyzed by means of extended molecular dynamics (MD) simulations for several representative systems. The results showed that the εζOL1 refinement improves the backbone description of B-DNA double helices and G-DNA stem. In B-DNA simulations, we observed an average increase of the helical twist and narrowing of the major groove, thus achieving better agreement with X-ray and solution NMR data. The balance between populations of BI and BII backbone substates was shifted towards the BII state, in better agreement with ensemble-refined solution experimental results. Furthermore, the refined parameters decreased the backbone RMS deviations in B-DNA MD simulations. In the antiparallel guanine quadruplex (G-DNA) the εζOL1 modification improved the description of non-canonical α/γ backbone substates, which were shown to be coupled to the ε/ζ torsion potential. Thus, the refinement is suggested as a possible alternative to the current ε/ζ torsion potential, which may enable more accurate modeling of nucleic acids. However, long-term testing is recommended before its routine application in DNA simulations.

• Publication type
• Journal Article MeSH

Riboswitches often occur in the 5'-untranslated regions of bacterial mRNA where they regulate gene expression. The preQ(1) riboswitch controls the biosynthesis of a hypermodified nucleoside queuosine in response to binding the queuosine metabolic intermediate. Structures of the ligand-bound and ligand-free states of the preQ(1) riboswitch from Thermoanaerobacter tengcongensis were determined recently by X-ray crystallography. We used multiple, microsecond-long molecular dynamics simulations (29 μs in total) to characterize the structural dynamics of preQ(1) riboswitches in both states. We observed different stabilities of the stem in the bound and free states, resulting in different accessibilities of the ribosome-binding site. These differences are related to different stacking interactions between nucleotides of the stem and the associated loop, which itself adopts different conformations in the bound and free states. We suggest that the loop not only serves to bind preQ(1) but also transmits information about ligand binding from the ligand-binding pocket to the stem, which has implications for mRNA accessibility to the ribosome. We explain functional results obscured by a high salt crystallization medium and help to refine regions of disordered electron density, which demonstrates the predictive power of our approach. Besides investigating the functional dynamics of the riboswitch, we have also utilized this unique small folded RNA system for analysis of performance of the RNA force field on the μs time scale. The latest AMBER parmbsc0χ(OL3) RNA force field is capable of providing stable trajectories of the folded molecule on the μs time scale. On the other hand, force fields that are not properly balanced lead to significant structural perturbations on the sub-μs time scale, which could easily lead to inappropriate interpretation of the simulation data.

• MeSH
• RNA, Bacterial chemistry   MeSH
• Crystallography, X-Ray MeSH
• Models, Molecular MeSH
• Riboswitch * MeSH
• Molecular Dynamics Simulation * MeSH
• Thermoanaerobacter chemistry   MeSH
• Publication type
• Journal Article MeSH
• Research Support, Non-U.S. Gov't MeSH
• Research Support, N.I.H., Extramural MeSH
• Names of Substances
• RNA, Bacterial MeSH
• Riboswitch * MeSH

Refinement of empirical force fields for nucleic acids requires their extensive testing using as wide range of systems as possible. However, finding unambiguous reference data is not easy. In this paper, we analyze four systems which we suggest should be included in standard portfolio of molecules to test nucleic acids force fields, namely, parallel and antiparallel stranded DNA guanine quadruplex stems, RNA quadruplex stem, and Z-DNA. We highlight parameters that should be monitored to assess the force field performance. The work is primarily based on 8.4 μs of 100-250 ns trajectories analyzed in detail followed by 9.6 μs of additional selected back up trajectories that were monitored to verify that the results of the initial analyses are correct. Four versions of the Cornell et al. AMBER force field are tested, including an entirely new parmχ(OL4) variant with χ dihedral specifically reparametrized for DNA molecules containing syn nucleotides. We test also different water models and ion conditions. While improvement for DNA quadruplexes is visible, the force fields still do not fully represent the intricate Z-DNA backbone conformation.

• Publication type
• Journal Article MeSH

The L1 stalk is a key mobile element of the large ribosomal subunit which interacts with tRNA during translocation. Here, we investigate the structure and mechanical properties of the rRNA H76/H75/H79 three-way junction at the base of the L1 stalk from four different prokaryotic organisms. We propose a coarse-grained elastic model and parameterize it using large-scale atomistic molecular dynamics simulations. Global properties of the junction are well described by a model in which the H76 helix is represented by a straight, isotropically flexible elastic rod, while the junction core is represented by an isotropically flexible spherical hinge. Both the core and the helix contribute substantially to the overall H76 bending fluctuations. The presence of wobble pairs in H76 does not induce any increased flexibility or anisotropy to the helix. The half-closed conformation of the L1 stalk seems to be accessible by thermal fluctuations of the junction itself, without any long-range allosteric effects. Bending fluctuations of H76 with a bulge introduced in it suggest a rationale for the precise position of the bulge in eukaryotes. Our elastic model can be generalized to other RNA junctions found in biological systems or in nanotechnology.

• MeSH
• Biomechanical Phenomena MeSH
• Nucleic Acid Conformation MeSH
• Ribosomal Proteins chemistry   MeSH
• RNA, Ribosomal, 23S chemistry   MeSH
• Molecular Dynamics Simulation MeSH
• Ribosome Subunits, Large, Archaeal chemistry   MeSH
• Ribosome Subunits, Large, Bacterial chemistry   MeSH
• Publication type
• Journal Article MeSH
• Research Support, Non-U.S. Gov't MeSH
• Names of Substances
• ribosomal protein L1 MeSH Browser
• Ribosomal Proteins MeSH
• RNA, Ribosomal, 23S MeSH

The article reviews the application of biomolecular simulation methods to understand the structure, dynamics and interactions of nucleic acids with a focus on explicit solvent molecular dynamics simulations of guanine quadruplex (G-DNA and G-RNA) molecules. While primarily dealing with these exciting and highly relevant four-stranded systems, where recent and past simulations have provided several interesting results and novel insight into G-DNA structure, the review provides some general perspectives on the applicability of the simulation techniques to nucleic acids.

• MeSH
• DNA chemistry   MeSH
• G-Quadruplexes * MeSH
• Guanine chemistry   MeSH
• Nucleic Acid Conformation MeSH
• Ligands MeSH
• RNA chemistry   MeSH
• Solvents chemistry   MeSH
• Molecular Dynamics Simulation * MeSH
• Telomere chemistry   MeSH
• Hydrogen Bonding MeSH
• Publication type
• Journal Article MeSH
• Research Support, Non-U.S. Gov't MeSH
• Review MeSH
• Research Support, N.I.H., Extramural MeSH
• Research Support, U.S. Gov't, Non-P.H.S. MeSH
• Names of Substances
• DNA MeSH
• Guanine MeSH
• Ligands MeSH
• RNA MeSH
• Solvents MeSH

Folded RNA molecules are shaped by an astonishing variety of highly conserved noncanonical molecular interactions and backbone topologies. The dinucleotide platform is a widespread recurrent RNA modular building submotif formed by the side-by-side pairing of bases from two consecutive nucleotides within a single strand, with highly specific sequence preferences. This unique arrangement of bases is cemented by an intricate network of noncanonical hydrogen bonds and facilitated by a distinctive backbone topology. The present study investigates the gas-phase intrinsic stabilities of the three most common RNA dinucleotide platforms - 5'-GpU-3', ApA, and UpC - via state-of-the-art quantum-chemical (QM) techniques. The mean stability of base-base interactions decreases with sequence in the order GpU > ApA > UpC. Bader's atoms-in-molecules analysis reveals that the N2(G)…O4(U) hydrogen bond of the GpU platform is stronger than the corresponding hydrogen bonds in the other two platforms. The mixed-pucker sugar-phosphate backbone conformation found in most GpU platforms, in which the 5'-ribose sugar (G) is in the C2'-endo form and the 3'-sugar (U) in the C3'-endo form, is intrinsically more stable than the standard A-RNA backbone arrangement, partially as a result of a favorable O2'…O2P intra-platform interaction. Our results thus validate the hypothesis of Lu et al. (Lu Xiang-Jun, et al. Nucleic Acids Res. 2010, 38, 4868-4876), that the superior stability of GpU platforms is partially mediated by the strong O2'…O2P hydrogen bond. In contrast, ApA and especially UpC platform-compatible backbone conformations are rather diverse and do not display any characteristic structural features. The average stabilities of ApA and UpC derived backbone conformers are also lower than those of GpU platforms. Thus, the observed structural and evolutionary patterns of the dinucleotide platforms can be accounted for, to a large extent, by their intrinsic properties as described by modern QM calculations. In contrast, we show that the dinucleotide platform is not properly described in the course of atomistic explicit-solvent simulations. Our work also gives methodological insights into QM calculations of experimental RNA backbone geometries. Such calculations are inherently complicated by rather large data and refinement uncertainties in the available RNA experimental structures, which often preclude reliable energy computations.

• Publication type
• Journal Article MeSH