The transition from B-DNA to A-DNA occurs in many protein-DNA interactions or in DNA/RNA hybrid duplexes, and thus plays a role in many important biomolecular processes that convey the biological function of DNA. However, the stability of A-DNA is severely underestimated in current AMBER force fields such as OL15, OL21 or bsc1, potentially leading to unstable or deformed protein-DNA complexes. In this study, we refine the deoxyribose dihedral potential to increase the stability of the north (N) puckering present in A-DNA. The new parameters, termed OL24, model A/B equilibrium in B-DNA duplexes in water in good agreement with nuclear magnetic resonance (NMR) experiment. They also improve the description of DNA/RNA hybrids and the transition of the DNA duplex to the A-form in concentrated ethanol solutions. These refinements significantly improve the modeling of protein-DNA complexes, increasing their structural stability and A-form population, while maintaining accurate representation of canonical B-DNA duplexes. Overall, the new parameters should allow more reliable modeling of the thermodynamic equilibrium between A- and B-DNA forms and the interactions of DNA with proteins.

• MeSH
• A-DNA * chemie   MeSH
• B-DNA * chemie   MeSH
• deoxyribosa * chemie   MeSH
• DNA * chemie   MeSH
• konformace nukleové kyseliny MeSH
• simulace molekulární dynamiky MeSH
• termodynamika MeSH
• Publikační typ
• časopisecké články MeSH
• Názvy látek
• A-DNA * MeSH
• B-DNA * MeSH
• deoxyribosa * MeSH
• DNA * 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
• katalýza MeSH
• kvantová teorie * MeSH
• kyslík chemie   MeSH
• protony MeSH
• RNA katalytická chemie   metabolismus   MeSH
• simulace molekulární dynamiky * MeSH
• termodynamika 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
• hairpin ribozyme MeSH Prohlížeč
• kyslík MeSH
• protony MeSH
• RNA katalytická 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
• konformace nukleové kyseliny MeSH
• kvantová teorie * MeSH
• nukleové kyseliny chemie   MeSH
• RNA chemie   MeSH
• výpočetní biologie * 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, Non-P.H.S. MeSH
• Názvy látek
• nukleové kyseliny 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
• chemické modely * MeSH
• DNA chemie   ultrastruktura   MeSH
• kinetika MeSH
• konformace nukleové kyseliny MeSH
• kvantová teorie MeSH
• mechanický stres MeSH
• modul pružnosti MeSH
• molekulární modely * MeSH
• počítačová simulace MeSH
• Publikační typ
• časopisecké články MeSH
• práce podpořená grantem MeSH
• Názvy látek
• 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
• biofyzika metody   MeSH
• fosfáty chemie   MeSH
• fosforylace MeSH
• hořčík chemie   MeSH
• katalýza MeSH
• konformace nukleové kyseliny MeSH
• kvantová teorie MeSH
• lidé MeSH
• molekulární modely MeSH
• počítačová simulace MeSH
• ribozomy chemie   MeSH
• RNA katalytická chemie   MeSH
• RNA virová chemie   MeSH
• RNA chemie   MeSH
• software MeSH
• Check Tag
• lidé MeSH
• Publikační typ
• časopisecké články MeSH
• práce podpořená grantem MeSH
• přehledy MeSH
• Research Support, N.I.H., Extramural MeSH
• Názvy látek
• fosfáty MeSH
• hammerhead ribozyme MeSH Prohlížeč
• hořčík MeSH
• RNA katalytická MeSH
• RNA virová MeSH
• RNA MeSH

The hepatitis delta virus (HDV) ribozyme is an RNA motif embedded in human pathogenic HDV RNA. Previous experimental studies have established that the active-site nucleotide C75 is essential for self-cleavage of the ribozyme, although its exact catalytic role in the process remains debated. Structural data from X-ray crystallography generally indicate that C75 acts as the general base that initiates catalysis by deprotonating the 2'-OH nucleophile at the cleavage site, while a hydrated magnesium ion likely protonates the 5'-oxygen leaving group. In contrast, some mechanistic studies support the role of C75 acting as general acid and thus being protonated before the reaction. We report combined quantum chemical/molecular mechanical calculations for the C75 general base pathway, utilizing the available structural data for the wild type HDV genomic ribozyme as a starting point. Several starting configurations differing in magnesium ion placement were considered and both one-dimensional and two-dimensional potential energy surface scans were used to explore plausible reaction paths. Our calculations show that C75 is readily capable of acting as the general base, in concert with the hydrated magnesium ion as the general acid. We identify a most likely position for the magnesium ion, which also suggests it acts as a Lewis acid. The calculated energy barrier of the proposed mechanism, approximately 20 kcal/mol, would lower the reaction barrier by approximately 15 kcal/mol compared with the uncatalyzed reaction and is in good agreement with experimental data.

• MeSH
• cytosin metabolismus   MeSH
• genomika MeSH
• katalýza MeSH
• kvantová teorie * MeSH
• molekulární modely * MeSH
• RNA katalytická chemie   metabolismus   MeSH
• sekundární struktura proteinů MeSH
• studie proveditelnosti MeSH
• vazebná místa MeSH
• virus hepatitidy delta enzymologie   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
• cytosin MeSH
• RNA katalytická MeSH