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

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

The hairpin ribozyme is a prominent member of small ribozymes since it does not require metal ions to achieve catalysis. Guanine 8 (G8) and adenine 38 (A38) have been identified as key participants in self-cleavage and -ligation. We have carried out hybrid quantum-mechanical/molecular mechanical (QM/MM) calculations to evaluate the energy along several putative reaction pathways. The error of our DFT description of the QM region was tested and shown to be ~1 kcal/mol. We find that self-cleavage of the hairpin ribozyme may follow several competing microscopic reaction mechanisms, all with calculated activation barriers in good agreement with those from experiment (20-21 kcal/mol). The initial nucleophilic attack of the A-1(2'-OH) group on the scissile phosphate is predicted to be rate-limiting in all these mechanisms. An unprotonated G8(-) (together with A38H(+)) yields a feasible activation barrier (20.4 kcal/mol). Proton transfer to a nonbridging phosphate oxygen also leads to feasible reaction pathways. Finally, our calculations consider thio-substitutions of one or both nonbridging oxygens of the scissile phosphate and predict that they have only a negligible effect on the reaction barrier, as observed experimentally.

• MeSH
• Catalysis MeSH
• Quantum Theory * MeSH
• Oxygen chemistry   MeSH
• Protons MeSH
• RNA, Catalytic chemistry   metabolism   MeSH
• Molecular Dynamics Simulation * 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
• hairpin ribozyme MeSH Browser
• Oxygen MeSH
• Protons MeSH
• RNA, Catalytic MeSH

In this feature article, we provide a side-by-side introduction for two research fields: quantum chemical calculations of molecular interaction in nucleic acids and RNA structural bioinformatics. Our main aim is to demonstrate that these research areas, while largely separated in contemporary literature, have substantial potential to complement each other that could significantly contribute to our understanding of the exciting world of nucleic acids. We identify research questions amenable to the combined application of modern ab initio methods and bioinformatics analysis of experimental structures while also assessing the limitations of these approaches. The ultimate aim is to attain valuable physicochemical insights regarding the nature of the fundamental molecular interactions and how they shape RNA structures, dynamics, function, and evolution.

• MeSH
• Nucleic Acid Conformation MeSH
• Quantum Theory * MeSH
• Nucleic Acids chemistry   MeSH
• RNA chemistry   MeSH
• Computational Biology * 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
• Nucleic Acids MeSH
• RNA MeSH

We combined atomistic molecular-dynamics simulations with quantum-mechanical calculations to investigate the sequence dependence of the stretching behavior of duplex DNA. Our combined quantum-mechanical/molecular-mechanical approach demonstrates that molecular-mechanical force fields are able to describe both the backbone and base-base interactions within the highly distorted nucleic acid structures produced by stretching the DNA from the 5' ends, which include conformations containing disassociated basepairs, just as well as these force fields describe relaxed DNA conformations. The molecular-dynamics simulations indicate that the force-induced melting pathway is sequence-dependent and is influenced by the availability of noncanonical hydrogen-bond interactions that can assist the disassociation of the DNA basepairs. The biological implications of these results are discussed.

• MeSH
• Models, Chemical * MeSH
• DNA chemistry   ultrastructure   MeSH
• Kinetics MeSH
• Nucleic Acid Conformation MeSH
• Quantum Theory MeSH
• Stress, Mechanical MeSH
• Elastic Modulus MeSH
• Models, Molecular * MeSH
• Computer Simulation MeSH
• Publication type
• Journal Article MeSH
• Research Support, Non-U.S. Gov't MeSH
• Names of Substances
• DNA MeSH

Hybrid QM/MM methods combine the rigor of quantum mechanical (QM) calculations with the low computational cost of empirical molecular mechanical (MM) treatment allowing to capture dynamic properties to probe critical atomistic details of enzyme reactions. Catalysis by RNA enzymes (ribozymes) has only recently begun to be addressed with QM/MM approaches and is thus still a field under development. This review surveys methodology as well as recent advances in QM/MM applications to RNA mechanisms, including those of the HDV, hairpin, and hammerhead ribozymes, as well as the ribosome. We compare and correlate QM/MM results with those from QM and/or molecular dynamics (MD) simulations, and discuss scope and limitations with a critical eye on current shortcomings in available methodologies and computer resources. We thus hope to foster mutual appreciation and facilitate collaboration between experimentalists and theorists to jointly advance our understanding of RNA catalysis at an atomistic level.

• MeSH
• Biophysics methods   MeSH
• Phosphates chemistry   MeSH
• Phosphorylation MeSH
• Magnesium chemistry   MeSH
• Catalysis MeSH
• Nucleic Acid Conformation MeSH
• Quantum Theory MeSH
• Humans MeSH
• Models, Molecular MeSH
• Computer Simulation MeSH
• Ribosomes chemistry   MeSH
• RNA, Catalytic chemistry   MeSH
• RNA, Viral chemistry   MeSH
• RNA chemistry   MeSH
• Software MeSH
• Check Tag
• Humans 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
• Phosphates MeSH
• hammerhead ribozyme MeSH Browser
• Magnesium MeSH
• RNA, Catalytic MeSH
• RNA, Viral MeSH
• RNA MeSH

This review provides a critical assessment of the advantages and limitations of modeling methods available for guanine quadruplex (G-DNA) molecules. We characterize the relations of simulations to the experimental techniques and explain the actual meaning and significance of the results. The following aspects are discussed: pair-additive approximation of the empirical force fields, sampling limitations stemming from the simulation time and accuracy of description of base stacking, H-bonding, sugar-phosphate backbone and ions by force fields. Several methodological approaches complementing the classical explicit solvent molecular dynamics simulations are commented on, including enhanced sampling methods, continuum solvent methods, free energy calculations and gas phase simulations. The successes and pitfalls of recent simulation studies of G-DNA are demonstrated on selected results, including studies of cation interactions and dynamics of G-DNA stems, studies of base substitutions (inosine, thioguanine and mixed tetrads), analysis of possible kinetic intermediates in folding pathway of a G-DNA stem and analysis of loop regions of G-DNA molecules.

• MeSH
• DNA chemistry   MeSH
• G-Quadruplexes * MeSH
• Guanine chemistry   MeSH
• Ligands MeSH
• Models, Molecular MeSH
• Computer Simulation MeSH
• Thermodynamics MeSH
• Publication type
• Journal Article MeSH
• Research Support, Non-U.S. Gov't MeSH
• Names of Substances
• DNA MeSH
• Guanine MeSH
• Ligands 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
• Magnesium chemistry   MeSH
• Cations, Divalent chemistry   MeSH
• Cations, Monovalent chemistry   MeSH
• Nucleic Acid Conformation MeSH
• Models, Molecular MeSH
• Molecular Sequence Data MeSH
• RNA, Catalytic chemistry   MeSH
• Base Sequence MeSH
• Sodium chemistry   MeSH
• Binding Sites MeSH
• Hepatitis Delta Virus enzymology   MeSH
• Water chemistry   MeSH
• Hydrogen Bonding MeSH
• Publication type
• Journal Article MeSH
• Research Support, Non-U.S. Gov't MeSH
• Research Support, N.I.H., Extramural MeSH
• Names of Substances
• Magnesium MeSH
• Cations, Divalent MeSH
• Cations, Monovalent MeSH
• RNA, Catalytic MeSH
• Sodium MeSH
• Water MeSH

Explicit solvent molecular dynamics (MD) simulations were carried out for sarcin-ricin domain (SRD) motifs from 23S (Escherichia coli) and 28S (rat) rRNAs. The SRD motif consists of GAGA tetraloop, G-bulged cross-strand A-stack, flexible region and duplex part. Detailed analysis of the overall dynamics, base pairing, hydration, cation binding and other SRD features is presented. The SRD is surprisingly static in multiple 25 ns long simulations and lacks any non-local motions, with root mean square deviation (r.m.s.d.) values between averaged MD and high-resolution X-ray structures of 1-1.4 A. Modest dynamics is observed in the tetraloop, namely, rotation of adenine in its apex and subtle reversible shift of the tetraloop with respect to the adjacent base pair. The deformed flexible region in low-resolution rat X-ray structure is repaired by simulations. The simulations reveal few backbone flips, which do not affect positions of bases and do not indicate a force field imbalance. Non-Watson-Crick base pairs are rigid and mediated by long-residency water molecules while there are several modest cation-binding sites around SRD. In summary, SRD is an unusually stiff rRNA building block. Its intrinsic structural and dynamical signatures seen in simulations are strikingly distinct from other rRNA motifs such as Loop E and Kink-turns.

• MeSH
• Endoribonucleases metabolism   MeSH
• Escherichia coli genetics   MeSH
• Fungal Proteins metabolism   MeSH
• Cations chemistry   MeSH
• Nucleic Acid Conformation MeSH
• Rats MeSH
• Crystallography, X-Ray MeSH
• Models, Molecular * MeSH
• Base Pairing MeSH
• Computer Simulation MeSH
• Ricin metabolism   MeSH
• RNA, Ribosomal, 23S chemistry   metabolism   MeSH
• RNA, Ribosomal, 28S chemistry   metabolism   MeSH
• Carbohydrates chemistry   MeSH
• Binding Sites MeSH
• Water chemistry   MeSH
• Hydrogen Bonding MeSH
• Animals MeSH
• Check Tag
• Rats MeSH
• Animals MeSH
• Publication type
• Journal Article MeSH
• Research Support, Non-U.S. Gov't MeSH
• Research Support, U.S. Gov't, Non-P.H.S. MeSH
• Names of Substances
• alpha-sarcin MeSH Browser
• Endoribonucleases MeSH
• Fungal Proteins MeSH
• Cations MeSH
• Ricin MeSH
• RNA, Ribosomal, 23S MeSH
• RNA, Ribosomal, 28S MeSH
• Carbohydrates MeSH
• Water MeSH