1
|
Yuan Y, Mills MJL, Zhang Z, Ma Y, Zhao C, Su W. A general RNA force field: comprehensive analysis of energy minima of molecular fragments of RNA. J Mol Model 2021; 27:137. [PMID: 33903935 DOI: 10.1007/s00894-021-04746-9] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/20/2020] [Accepted: 04/14/2021] [Indexed: 11/29/2022]
Abstract
Force fields are actively used to study RNA. Development of accurate force fields relies on a knowledge of how the variation of properties of molecules depends on their structure. Detailed scrutiny of RNA's conformational preferences is needed to guide such development. Towards this end, minimum energy structures for each of a set of 16 small RNA-derived molecules were obtained by geometry optimization at the HF/6-31G(d,p), B3LYP/apc-1, and MP2/cc-pVDZ levels of theory. The number of minima computed for a given fragment was found to be related to both its size and flexibility. Atomic electrostatic multipole moments of atoms occurring in the [HO-P(O3)-CH2-] fragment of 30 sugar-phosphate-sugar geometries were calculated at the HF/6-31G(d,p) and B3LYP/apc-1 levels of theory, and the transferability of these properties between different conformations was investigated. The atomic multipole moments were found to be highly transferable between different conformations with small standard deviations. These results indicate necessary elements of the development of accurate RNA force fields.
Collapse
Affiliation(s)
- Yongna Yuan
- School of Information Science & Engineering, Lanzhou University, No. 222 South Tianshui Road, Lanzhou, 730000, China.
| | - Matthew J L Mills
- 3M Corporate Research Analytical Laboratory, Saint Paul, MN, 55114, USA
| | - Zhuangzhuang Zhang
- School of Information Science & Engineering, Lanzhou University, No. 222 South Tianshui Road, Lanzhou, 730000, China.,Xi'an Microelectronic Technology Institute, No.198 Taibai South Road, Xi'an, 710000, China
| | - Yan Ma
- School of Information Science & Engineering, Lanzhou University, No. 222 South Tianshui Road, Lanzhou, 730000, China
| | - Chunyan Zhao
- School of Pharmacy, Lanzhou University, No. 222 South Tianshui Road, Lanzhou, 730000, China.
| | - Wei Su
- School of Information Science & Engineering, Lanzhou University, No. 222 South Tianshui Road, Lanzhou, 730000, China
| |
Collapse
|
2
|
Ortiz de Luzuriaga I, Lopez X, Gil A. Learning to Model G-Quadruplexes: Current Methods and Perspectives. Annu Rev Biophys 2021; 50:209-243. [PMID: 33561349 DOI: 10.1146/annurev-biophys-060320-091827] [Citation(s) in RCA: 14] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Abstract
G-quadruplexes have raised considerable interest during the past years for the development of therapies against cancer. These noncanonical structures of DNA may be found in telomeres and/or oncogene promoters, and it has been observed that the stabilization of such G-quadruplexes may disturb tumor cell growth. Nevertheless, the mechanisms leading to folding and stabilization of these G-quadruplexes are still not well established, and they are the focus of much current work in this field. In seminal works, stabilization was observed to be produced by cations. However, subsequent studies showed that different kinds of small molecules, from planar and nonplanar organic molecules to square-planar and octahedral metal complexes, may also lead to the stabilization of G-quadruplexes. Thus, the comprehension and rationalization of the interaction of these small molecules with G-quadruplexes are also important topics of current interest in medical applications. To shed light on the questions arising from the literature on the formation of G-quadruplexes, their stabilization, and their interaction with small molecules, synergies between experimental studies and computational works are needed. In this review, we mainly focus on in silico approaches and provide a broad compilation of different leading studies carried out to date by different computational methods. We divide these methods into twomain categories: (a) classical methods, which allow for long-timescale molecular dynamics simulations and the corresponding analysis of dynamical information, and (b) quantum methods (semiempirical, quantum mechanics/molecular mechanics, and density functional theory methods), which allow for the explicit simulation of the electronic structure of the system but, in general, are not capable of being used in long-timescale molecular dynamics simulations and, therefore, give a more static picture of the relevant processes.
Collapse
Affiliation(s)
- Iker Ortiz de Luzuriaga
- CIC nanoGUNE BRTA, 20018 Donostia, Euskadi, Spain; .,Polimero eta Material Aurreratuak: Fisika, Kimika eta Teknologia, Kimika Fakultatea, Euskal Herriko Uniberstitatea, UPV/EHU, 20080 Donostia, Euskadi, Spain
| | - Xabier Lopez
- Polimero eta Material Aurreratuak: Fisika, Kimika eta Teknologia, Kimika Fakultatea, Euskal Herriko Uniberstitatea, UPV/EHU, 20080 Donostia, Euskadi, Spain.,Donostia International Physics Center, 20018 Donostia, Spain
| | - Adrià Gil
- CIC nanoGUNE BRTA, 20018 Donostia, Euskadi, Spain; .,BioISI-Biosystems and Integrative Sciences Institute, Faculdade de Ciências, Universidade de Lisboa, 1749-016 Lisboa, Portugal;
| |
Collapse
|
3
|
Abstract
Our current knowledge on the unique roles of RNA in cells makes it vital to investigate the properties of RNA systems using computational methods because of the potential pharmaceutical applications. With the continuous advancement of computer technology, it is now possible to study RNA folding. Molecular mechanics calculations are useful in discovering the structural and thermodynamic properties of RNA systems. Yet, the predictions depend on the quality of the RNA force field, which is a set of parameters describing the potential energy of the system. Torsional parameters are one of the terms in a force field that can be revised using physics-based approaches. This chapter focuses on improvements provided by revisions of torsional parameters of the AMBER (Assisted Model Building with Energy Refinement) RNA force field. The theory behind torsional revisions and re-parameterization of several RNA torsions is briefly described. Applications of the revised torsional parameters to study RNA nucleosides, single-stranded RNA tetramers, and RNA repeat expansions are described in detail. It is concluded that RNA force fields require constant revisions and should be benchmarked against diverse RNA systems such as single strands and internal loops in order to test their qualities.
Collapse
|
4
|
Šponer J, Bussi G, Krepl M, Banáš P, Bottaro S, Cunha RA, Gil-Ley A, Pinamonti G, Poblete S, Jurečka P, Walter NG, Otyepka M. RNA Structural Dynamics As Captured by Molecular Simulations: A Comprehensive Overview. Chem Rev 2018; 118:4177-4338. [PMID: 29297679 PMCID: PMC5920944 DOI: 10.1021/acs.chemrev.7b00427] [Citation(s) in RCA: 336] [Impact Index Per Article: 56.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/14/2017] [Indexed: 12/14/2022]
Abstract
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.
Collapse
Affiliation(s)
- Jiří Šponer
- Institute of Biophysics of the Czech Academy of Sciences , Kralovopolska 135 , Brno 612 65 , Czech Republic
| | - Giovanni Bussi
- Scuola Internazionale Superiore di Studi Avanzati , Via Bonomea 265 , Trieste 34136 , Italy
| | - Miroslav Krepl
- Institute of Biophysics of the Czech Academy of Sciences , Kralovopolska 135 , Brno 612 65 , Czech Republic
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science , Palacky University Olomouc , 17. listopadu 12 , Olomouc 771 46 , Czech Republic
| | - Pavel Banáš
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science , Palacky University Olomouc , 17. listopadu 12 , Olomouc 771 46 , Czech Republic
| | - Sandro Bottaro
- Structural Biology and NMR Laboratory, Department of Biology , University of Copenhagen , Copenhagen 2200 , Denmark
| | - Richard A Cunha
- Scuola Internazionale Superiore di Studi Avanzati , Via Bonomea 265 , Trieste 34136 , Italy
| | - Alejandro Gil-Ley
- Scuola Internazionale Superiore di Studi Avanzati , Via Bonomea 265 , Trieste 34136 , Italy
| | - Giovanni Pinamonti
- Scuola Internazionale Superiore di Studi Avanzati , Via Bonomea 265 , Trieste 34136 , Italy
| | - Simón Poblete
- Scuola Internazionale Superiore di Studi Avanzati , Via Bonomea 265 , Trieste 34136 , Italy
| | - Petr Jurečka
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science , Palacky University Olomouc , 17. listopadu 12 , Olomouc 771 46 , Czech Republic
| | - Nils G Walter
- Single Molecule Analysis Group and Center for RNA Biomedicine, Department of Chemistry , University of Michigan , Ann Arbor , Michigan 48109 , United States
| | - Michal Otyepka
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science , Palacky University Olomouc , 17. listopadu 12 , Olomouc 771 46 , Czech Republic
| |
Collapse
|
5
|
Negi I, Kathuria P, Sharma P, Wetmore SD. How do hydrophobic nucleobases differ from natural DNA nucleobases? Comparison of structural features and duplex properties from QM calculations and MD simulations. Phys Chem Chem Phys 2018; 19:16365-16374. [PMID: 28657627 DOI: 10.1039/c7cp02576a] [Citation(s) in RCA: 23] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
Computational (DFT and MD simulation) methods are employed to systematically characterize the structural and energetic properties of five hydrophobic nucleobases (FEMO, MMO2, NaM, 5SICS and TPT3) that constitute four unnatural base pairs (FEMO:5SICS, MMO2:5SICS, NaM:5SICS and TPT3:NaM). These hydrophobic bases have been recently shown to be replicated when present between natural bases in DNA duplexes, with the highest replication fidelity and efficiency occuring for the TPT3:NaM pair. Our QM calculations suggest that the preferred (anti) glycosidic orientations of nucleosides containing hydrophobic bases are similar to the natural DNA nucleosides despite differences in their chemical structures. However, due to the inability to form interbase hydrogen bonds, hydrophobic base pairs intrinsically prefer nonplanar, distorted geometries, many of which are stabilized through π-π stacking interactions. Furthermore, the intrinsic stacking potential between a hydrophobic and a natural base is similar to that between two natural bases, indicating that the strength of stacking interactions in DNA duplexes containing hydrophobic bases is likely comparable to natural DNA. However, in contrast to the isolated base-pair geometries, our MD simulations suggest that the hydrophobic base pairs adopt variable geometries within DNA, which range from stacked (5SICS:FEMO) to nearly planar (5SICS:NaM and SICS:MMO2) to planar (TPT3:NaM). As a result, the duplex structural features at the site of modification depend on the identity of the hydrophobic base pair, where the TPT3:NaM pair causes the least structural changes compared to natural DNA. Overall, the structural insight obtained from our calculations on DNA containing hydrophobic base pairs explains the experimentally-observed higher fidelity and efficiency during replication of TPT3:NaM compared to other hydrophobic nucleobase pairs. By providing valuable structural information that explains the intrinsic and duplex properties of this class of unnatural nucleobases, the present work may aid the future design of improved hydrophobic analogues.
Collapse
Affiliation(s)
- Indu Negi
- Computational Biochemistry Laboratory, Department of Chemistry and Centre for Advanced Studies in Chemistry, Panjab University, Chandigarh, 160014, India.
| | | | | | | |
Collapse
|
6
|
Deb I, Sarzynska J, Nilsson L, Lahiri A. Rapid communication capturing the destabilizing effect of dihydrouridine through molecular simulations. Biopolymers 2016; 101:985-91. [PMID: 24729441 DOI: 10.1002/bip.22495] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/16/2014] [Revised: 04/04/2014] [Accepted: 04/07/2014] [Indexed: 01/08/2023]
Abstract
The structural effects of the commonly occurring modified nucleoside dihydrouridine (D) observed experimentally in model oligonucleotides include a strong destabilization of the C3'-endo sugar conformation of D, the disruption of stacking interactions of neighboring residues with D and a possible destabilization of the C3'-endo sugar pucker of the 5'-neighboring nucleoside. Our simulations with a combination of a set of parameters for modified RNA residues with the recently developed AMBER FF99χ force field having reoptimized glycosidic torsion angle parameters for standard nucleosides was found to reproduce the destabilizing effect of dihydrouridine better than with the AMBER FF99 force field for nucleic acids for which the parameters for the modified residues were originally developed.
Collapse
Affiliation(s)
- Indrajit Deb
- Department of Biophysics, Molecular Biology and Bioinformatics, University of Calcutta, Kolkata, 700009, West Bengal, India
| | | | | | | |
Collapse
|
7
|
Deb I, Pal R, Sarzynska J, Lahiri A. Reparameterizations of theχTorsion and Lennard-JonesσParameters Improve the Conformational Characteristics of Modified Uridines. J Comput Chem 2016; 37:1576-88. [DOI: 10.1002/jcc.24374] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/15/2016] [Accepted: 03/05/2016] [Indexed: 01/12/2023]
Affiliation(s)
- Indrajit Deb
- Department of Biophysics, Molecular Biology and Bioinformatics; University of Calcutta; 92 APC Road Kolkata West Bengal 700009 India
- Institute of Bioorganic Chemistry, Polish Academy of Sciences; Noskowskiego 12/14 Poznan 61-704 Poland
| | - Rupak Pal
- Department of Biophysics, Molecular Biology and Bioinformatics; University of Calcutta; 92 APC Road Kolkata West Bengal 700009 India
| | - Joanna Sarzynska
- Institute of Bioorganic Chemistry, Polish Academy of Sciences; Noskowskiego 12/14 Poznan 61-704 Poland
| | - Ansuman Lahiri
- Department of Biophysics, Molecular Biology and Bioinformatics; University of Calcutta; 92 APC Road Kolkata West Bengal 700009 India
| |
Collapse
|
8
|
Nganou C, Kennedy SD, McCamant DW. Disagreement Between the Structure of the dTpT Thymine Pair Determined by NMR and Molecular Dynamics Simulations Using Amber 14 Force Fields. J Phys Chem B 2016; 120:1250-8. [PMID: 26836489 DOI: 10.1021/acs.jpcb.6b00191] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
We report a disagreement between the predicted structures of the dTpT thymine pair (thymidylyl(3' → 5')thymidine) using nuclear magnetic resonance (NMR) spectroscopy and molecular dynamics (MD) simulations using the AMBER ff14SB and ff14 + ε/ζOL1 + χOL4 force fields for DNA. The NMR structure was determined using NOE couplings to thymine's H6 and J(HH) couplings between sugar protons. The MD simulation used replica exchange methods to produce converged statistics in a 500 ns trajectory. NMR data indicate that both thymine nucleotides in the pair display an anti conformation of B-DNA, while the MD simulations predict a structure in which the 5'-thymine is flipped into a syn conformation and the 3'-thymine is in an anti conformation. The syn conformation of the 5'-thymine predicted by MD appears by a ∼ 180-deg flip of the glycosidic angle in comparison to the B-form anti structure. Differences in the distortion of the sugar pucker between 5'-thymine and 3'-thymine further highlighted the surprisingly different conformation of the 5'- and 3'-ends. While both MD and NMR indicate the deoxyribose sugars to be primarily in the 2'-endo conformation typical of B-form DNA, the MD simulations predict a more twisted conformation (2'-endo/1'-exo) for the 5'-sugar and significant flexibility of C3' of the 3'-sugar. We conclude that the current AMBER force field does not accurately predict the conformation of single-stranded thymine, in agreement with previous work investigating single-stranded DNA.
Collapse
Affiliation(s)
- Collins Nganou
- Department of Chemistry, University of Rochester , Rochester, New York 14627, United States
| | - Scott D Kennedy
- Department of Biochemistry and Biophysics, University of Rochester , Rochester, New York 14642, United States
| | - David W McCamant
- Department of Chemistry, University of Rochester , Rochester, New York 14627, United States
| |
Collapse
|
9
|
Schrodt MV, Andrews CT, Elcock AH. Large-Scale Analysis of 48 DNA and 48 RNA Tetranucleotides Studied by 1 μs Explicit-Solvent Molecular Dynamics Simulations. J Chem Theory Comput 2015; 11:5906-17. [PMID: 26580891 PMCID: PMC4806854 DOI: 10.1021/acs.jctc.5b00899] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
An understanding of how the conformational behavior of single-stranded DNAs and RNAs depend on sequence is likely to be important for attempts to rationalize the thermodynamics of nucleic acid folding. In an attempt to further our understanding of such sequence dependences, we report here the results of 192 (1 μs) explicit-solvent molecular dynamics (MD) simulations of 48 DNA and 48 RNA tetranucleotide sequences performed using recently reported modifications to the AMBER force field. Each tetranucleotide was simulated starting from two different conformations, a fully natively stacked and a completely unstacked conformation, and populations of the various possible base stacking arrangements were analyzed. The simulations indicate that, for both DNA and RNA, the populations of fully natively stacked conformations increase linearly with increasing numbers of purines in the sequence, while the conformational entropies, computed by two complementary methods, decrease. Despite the comparatively short simulation times, the computed free energies of stacking of the 16 possible combinations of bases in the middle of the sequences are found to be in good correspondence with values reported recently from simulations of dinucleoside monophosphates using the same force field. Finally, consistent with recent reports from other groups, non-native stacking interactions, i.e., between bases that are not adjacent in sequence, are shown to be a recurring feature of the simulations; in particular, stacking interactions of bases in a i:i+2 relationship are shown to occur significantly more frequently when the intervening base is a pyrimidine. Given that the high prevalence of non-native stacking interactions is thought to be unrealistic, it appears that further parametrization work will be required before accurate conformational descriptions of single-stranded nucleic acids can be obtained with current force fields.
Collapse
Affiliation(s)
| | - Casey T. Andrews
- Department of Biochemistry, University of Iowa, Iowa City, IA 52242
| | - Adrian H. Elcock
- Department of Biochemistry, University of Iowa, Iowa City, IA 52242
| |
Collapse
|
10
|
Bochicchio A, Rossetti G, Tabarrini O, Krauβ S, Carloni P. Molecular view of ligands specificity for CAG repeats in anti-Huntington therapy. J Chem Theory Comput 2015; 11:4911-22. [PMID: 26574279 DOI: 10.1021/acs.jctc.5b00208] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Abstract
Huntington's disease is a fatal and devastating neurodegenerative genetic disorder for which there is currently no cure. It is characterized by Huntingtin protein's mRNA transcripts with 36 or more CAG repeats. Inhibiting the formation of pathological complexes between these expanded transcripts and target proteins may be a valuable strategy against the disease. Yet, the rational design of molecules specifically targeting the expanded CAG repeats is limited by the lack of structural information. Here, we use well-tempered metadynamics-based free energy calculations to investigate pose and affinity of two ligands targeting CAG repeats for which affinities have been previously measured. The first consists of two 4-guanidinophenyl rings linked by an ester group. It is the most potent ligand identified so far, with Kd = 60(30) nM. The second consists of a 4-phenyl dihydroimidazole and 4-1H-indole dihydroimidazole connected by a C-C bond (Kd = 700(80) nM). Our calculations reproduce the experimental affinities and uncover the recognition pattern between ligands' and their RNA target. They also provide a molecular basis for the markedly different affinity of the two ligands for CAG repeats as observed experimentally. These findings may pave the way for a structure-based hit-to-lead optimization to further improve ligand selectivity toward CAG repeat-containing mRNAs.
Collapse
Affiliation(s)
| | - Giulia Rossetti
- Department of Oncology, Hematology and Stem Cell Transplantation, RWTH Aachen University , D-52074 Aachen, North Rhine-Westphalia, Germany
| | - Oriana Tabarrini
- Department of Pharmaceutical Sciences, Università di Perugia , Via del Liceo 1, I-06123 Perugia, Perugia, Italy
| | - Sybille Krauβ
- German Center for Neurodegenerative Diseases (DZNE), Sigmund-Freud-Strasse 25, D-53127 Bonn, North Rhine-Westphalia, Germany
| | | |
Collapse
|
11
|
Brown RF, Andrews CT, Elcock AH. Stacking Free Energies of All DNA and RNA Nucleoside Pairs and Dinucleoside-Monophosphates Computed Using Recently Revised AMBER Parameters and Compared with Experiment. J Chem Theory Comput 2015; 11:2315-28. [PMID: 26574427 PMCID: PMC4651843 DOI: 10.1021/ct501170h] [Citation(s) in RCA: 74] [Impact Index Per Article: 8.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Abstract
We report the results of a series of 1-μs-long explicit-solvent molecular dynamics (MD) simulations performed to compare the free energies of stacking (ΔGstack) of all possible combinations of DNA and RNA nucleoside (NS) pairs and dinucleoside-monophosphates (DNMPs). For both NS pairs and DNMPs, we show that the computed stacking free energies are in reasonable qualitative agreement with experimental measurements and appear to provide the closest correspondence with experimental data yet found among computational studies; in all cases, however, the computed stacking free energies are too favorable relative to experimental data. Comparisons of NS-pair systems indicate that stacking interactions are very similar in RNA and DNA systems except when a thymine or uracil base is involved: the presence of a thymine base favors stacking by ∼0.3 kcal/mol relative to a uracil base. One exception is found in the self-stacking of cytidines, which are found to be significantly more favorable for the DNA form; an analysis of the rotational orientations sampled during stacking events suggests that this is likely to be due to more favorable sugar-sugar interactions in stacked complexes of deoxycytidines. Comparisons of the DNMP systems indicate that stacking interactions are more favorable in RNA than in DNA except, again, when thymine or uracil bases are involved. Finally, additional simulations performed using a previous generation of the AMBER force field-in which the description of glycosidic bond rotations was less than optimal-produce computed stacking free energies that are in poorer agreement with experimental data. Overall, the simulations provide a comprehensive view of stacking thermodynamics in NS pairs and in DNMPs as predicted by a state-of-the-art MD force field.
Collapse
Affiliation(s)
- Reid F. Brown
- Department of Biochemistry, University of Iowa, Iowa City, IA 52242
| | - Casey T. Andrews
- Department of Biochemistry, University of Iowa, Iowa City, IA 52242
| | - Adrian H. Elcock
- Department of Biochemistry, University of Iowa, Iowa City, IA 52242
| |
Collapse
|
12
|
Fujiwara SI, Sawada K, Amisaki T. Molecular dynamics study on conformational differences between dGMP and 8-oxo-dGMP: Effects of metal ions. J Mol Graph Model 2014; 51:158-67. [PMID: 24929814 DOI: 10.1016/j.jmgm.2014.05.007] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/11/2013] [Revised: 05/09/2014] [Accepted: 05/22/2014] [Indexed: 01/05/2023]
Abstract
The modified nucleotide base 7,8-dihydro-8-oxo-guanine (8-oxo-G) is one of the major sources of spontaneous mutagenesis. Nucleotide-sanitizing enzymes, such as the MutT homolog-1 (MTH1) and nudix-type motif 5 (NUDT5), selectively remove 8-oxo-G from the cellular pool of nucleotides. Previous studies showed that, although the syn conformation generally predominates in purine nucleotides with a bulky substituent at the 8-position, 8-oxo-dGMP binds to both MTH1 and NUDT5 in the anti conformation. This study was initiated to investigate the possibility that 8-oxo-dGMP itself may adopt the anti conformation. Molecular dynamics simulations of mononucleotides (dGMP, 8-oxo-dGMP) in aqueous solution were performed. 8-oxo-dGMP adopted the anti conformation as well as the syn conformation, and the proportion of adopting the anti conformation increased in the presence of metal ions. When 8-oxo-dGMP was in the anti conformation, a metal ion was located between the oxygen atom of phosphate and the oxygen atom at the 8-position of 8-oxo-G. The types of stable anti conformations of 8-oxo-dGMP differed, depending on the ionic radii and charges of coexisting ions. These data suggested a role for metal ions, other than as cofactors for the hydrolysis of the di- and tri-phosphate forms of mononucleotides; that the metal ions help retain the anti conformation of the N-glycosidic torsion angle of 8-oxo-dGMP to promote the binding between the 8-oxo-G deoxynucleotide and the nucleotide-sanitizing enzymes.
Collapse
Affiliation(s)
- Shin-Ichi Fujiwara
- Department of Biological Regulation, Faculty of Medicine, Tottori University, 86 Nishi-cho, Yonago 683-8503, Japan.
| | - Kenichiro Sawada
- Department of Biological Regulation, Faculty of Medicine, Tottori University, 86 Nishi-cho, Yonago 683-8503, Japan
| | - Takashi Amisaki
- Department of Biological Regulation, Faculty of Medicine, Tottori University, 86 Nishi-cho, Yonago 683-8503, Japan
| |
Collapse
|
13
|
Deb I, Sarzynska J, Nilsson L, Lahiri A. Conformational preferences of modified uridines: comparison of AMBER derived force fields. J Chem Inf Model 2014; 54:1129-42. [PMID: 24697757 DOI: 10.1021/ci400582a] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The widespread occurrence of modified residues in RNA sequences necessitates development of accurate parameters for these modifications for reliable modeling of RNA structure and dynamics. A comprehensive set of parameters for the 107 naturally occurring RNA modifications was proposed by Aduri et al. (J. Chem. Theory Comput. 2007, 3, 1464-1475) for the AMBER FF99 force field. In this work, we tested these parameters on a set of modified uridine residues, namely, dihydrouridine, 2-thiouridine, 4-thiouridine, pseudouridine, and uridine-5-oxyacetic acid, by performing molecular dynamics and replica exchange molecular dynamics simulations of these nucleosides. Although our simulations using the FF99 force field did not, in general, reproduce the experimentally observed conformational characteristics well, combination of the parameter set with recent revisions of the FF99 force field for RNA showed noticeable improvement for some of the nucleosides.
Collapse
Affiliation(s)
- Indrajit Deb
- Department of Biophysics, Molecular Biology & Bioinformatics, University of Calcutta , Kolkata 700009, West Bengal, India
| | | | | | | |
Collapse
|
14
|
Mládek A, Banáš P, Jurečka P, Otyepka M, Zgarbová M, Šponer J. Energies and 2'-Hydroxyl Group Orientations of RNA Backbone Conformations. Benchmark CCSD(T)/CBS Database, Electronic Analysis, and Assessment of DFT Methods and MD Simulations. J Chem Theory Comput 2013; 10:463-80. [PMID: 26579924 DOI: 10.1021/ct400837p] [Citation(s) in RCA: 23] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Abstract
Sugar-phosphate backbone is an electronically complex molecular segment imparting RNA molecules high flexibility and architectonic heterogeneity necessary for their biological functions. The structural variability of RNA molecules is amplified by the presence of the 2'-hydroxyl group, capable of forming multitude of intra- and intermolecular interactions. Bioinformatics studies based on X-ray structure database revealed that RNA backbone samples at least 46 substates known as rotameric families. The present study provides a comprehensive analysis of RNA backbone conformational preferences and 2'-hydroxyl group orientations. First, we create a benchmark database of estimated CCSD(T)/CBS relative energies of all rotameric families and test performance of dispersion-corrected DFT-D3 methods and molecular mechanics in vacuum and in continuum solvent. The performance of the DFT-D3 methods is in general quite satisfactory. The B-LYP-D3 method provides the best trade-off between accuracy and computational demands. B3-LYP-D3 slightly outperforms the new PW6B95-D3 and MPW1B95-D3 and is the second most accurate density functional of the study. The best agreement with CCSD(T)/CBS is provided by DSD-B-LYP-D3 double-hybrid functional, although its large-scale applications may be limited by high computational costs. Molecular mechanics does not reproduce the fine energy differences between the RNA backbone substates. We also demonstrate that the differences in the magnitude of the hyperconjugation effect do not correlate with the energy ranking of the backbone conformations. Further, we investigated the 2'-hydroxyl group orientation preferences. For all families, we conducted a QM and MM hydroxyl group rigid scan in gas phase and solvent. We then carried out set of explicit solvent MD simulations of folded RNAs and analyze 2'-hydroxyl group orientations of different backbone families in MD. The solvent energy profiles determined primarily by the sugar pucker match well with the distribution data derived from the simulations. The QM and MM energy profiles predict the same 2'-hydroxyl group orientation preferences. Finally, we demonstrate that the high energy of unfavorable and rarely sampled 2'-hydroxyl group orientations can be attributed to clashes between occupied orbitals.
Collapse
Affiliation(s)
- Arnošt Mládek
- Institute of Biophysics, Academy of Sciences of the Czech Republic , Královopolská 135, 612 65 Brno, Czech Republic
| | - Pavel Banáš
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacky University , tr. 17. listopadu 12, 771 46, Olomouc, Czech Republic
| | - Petr Jurečka
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacky University , tr. 17. listopadu 12, 771 46, Olomouc, Czech Republic
| | - Michal Otyepka
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacky University , tr. 17. listopadu 12, 771 46, Olomouc, Czech Republic
| | - Marie Zgarbová
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacky University , tr. 17. listopadu 12, 771 46, Olomouc, Czech Republic
| | - Jiří Šponer
- Institute of Biophysics, Academy of Sciences of the Czech Republic , Královopolská 135, 612 65 Brno, Czech Republic.,CEITEC, Central European Institute of Technology , Campus Bohunice, Kamenice 5, 625 00 Brno, Czech Republic
| |
Collapse
|
15
|
Song J, Ji C, Zhang JZH. The critical effect of polarization on the dynamical structure of guanine quadruplex DNA. Phys Chem Chem Phys 2013; 15:3846-54. [PMID: 23399949 DOI: 10.1039/c2cp44100d] [Citation(s) in RCA: 30] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/20/2023]
Abstract
Guanine quadruplex DNA (G-DNA), found in eukaryotic telomeres and recently in non-telomeric genomic DNA, plays important biological roles and their structures are being explored as potential targets for therapeutic intervention. Since the quadruplex structure of G-DNA is stabilized by cations, electrostatic interaction is expected to play important roles in the dynamical structure of G-DNA. In current work, MD simulation was carried out to study the dynamical structure of a special G-DNA (with sequence d(G(4)T(4)G(4))) complexed with five K(+) ions. In order to properly include polarization in MD simulation, a new set of polarized nucleic acid specific charge based on fragment quantum chemistry calculation was developed for G-DNA. Our study showed that polarization of the nucleobases by K(+) enhanced electrostatic attraction between the base and ions. This increased attractive interaction is critical to stabilizing the stem-loop junction ions in G-DNA. Without this polarization effect, as in MD simulation using a standard (nonpolarizable) force field, the top and bottom cations escaped into the solvent within just a few nanoseconds. Furthermore, an incorrect bifurcated bonding geometry of G-DNA, found in previous MD simulation study under a standard force field but not observed in experiments, disappeared in the present stimulation using the new polarized force field. The current study bridged an important gap between the simulation study and experimental observation on the dynamical structure of G-DNA.
Collapse
Affiliation(s)
- Jianing Song
- State Key Laboratory of Precision Spectroscopy, East China Normal University, Shanghai 200062, China
| | | | | |
Collapse
|
16
|
Kara M, Zacharias M. Theoretical studies of nucleic acids folding. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2013. [DOI: 10.1002/wcms.1146] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/14/2023]
Affiliation(s)
- Mahmut Kara
- Physics Department T38, Technical University Munich, Garching, Germany
| | - Martin Zacharias
- Martin Zacharias, Physics Department T38, Technical University Munich, Garching, Germany
| |
Collapse
|
17
|
Zgarbová M, Luque FJ, Šponer J, Cheatham TE, Otyepka M, Jurečka P. Toward Improved Description of DNA Backbone: Revisiting Epsilon and Zeta Torsion Force Field Parameters. J Chem Theory Comput 2013; 9:2339-2354. [PMID: 24058302 PMCID: PMC3775469 DOI: 10.1021/ct400154j] [Citation(s) in RCA: 221] [Impact Index Per Article: 20.1] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Abstract
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.
Collapse
Affiliation(s)
- Marie Zgarbová
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacky University, 17 listopadu 12, 771 46, Olomouc, Czech Republic
| | - F. Javier Luque
- Departament de Fisicoquímica and Institut de Biomedicina (IBUB), Facultad de Farmacia, Universidad de Barcelona, Campus Torribera, Prat de la Riba 171, Edifici Verdaguer, 080921 Santa Coloma de Gramanet, Spain
| | - Jiří Šponer
- Institute of Biophysics, Academy of Sciences of the Czech Republic, Královopolská 135, 612 65 Brno, Czech Republic
- CEITEC – Central European Institute of Technology, Campus Bohunice, Kamenice 5, 625 00 Brno, Czech Republic
| | - Thomas E. Cheatham
- Departments of Medicinal Chemistry, University of Utah, 2000 East 30 South Skaggs 105, Salt Lake City, UT 84112, USA
| | - Michal Otyepka
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacky University, 17 listopadu 12, 771 46, Olomouc, Czech Republic
| | - Petr Jurečka
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacky University, 17 listopadu 12, 771 46, Olomouc, Czech Republic
| |
Collapse
|
18
|
Kührová P, Banáš P, Best RB, Šponer J, Otyepka M. Computer Folding of RNA Tetraloops? Are We There Yet? J Chem Theory Comput 2013; 9:2115-25. [DOI: 10.1021/ct301086z] [Citation(s) in RCA: 76] [Impact Index Per Article: 6.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Affiliation(s)
- Petra Kührová
- Regional Centre of Advanced Technologies
and Materials, Department of Physical Chemistry, Faculty of Science, Palacky University Olomouc, 17. listopadu 12, 771 46 Olomouc, Czech Republic
| | - Pavel Banáš
- Regional Centre of Advanced Technologies
and Materials, Department of Physical Chemistry, Faculty of Science, Palacky University Olomouc, 17. listopadu 12, 771 46 Olomouc, Czech Republic
- Institute
of Biophysics, Academy of Sciences of the Czech Republic, Kralovopolska
135, 612 65 Brno, Czech Republic
| | - Robert B. Best
- Laboratory of Chemical Physics,
National Institute of Diabetes and Digestive and Kidney Diseases, National Institutes of Health, Bethesda, Maryland 20892-0520,
United States
| | - Jiří Šponer
- Institute
of Biophysics, Academy of Sciences of the Czech Republic, Kralovopolska
135, 612 65 Brno, Czech Republic
- CEITEC − Central European
Institute of Technology, Masaryk University, Campus Bohunice, Kamenice 5, 625 00 Brno, Czech Republic
| | - Michal Otyepka
- Regional Centre of Advanced Technologies
and Materials, Department of Physical Chemistry, Faculty of Science, Palacky University Olomouc, 17. listopadu 12, 771 46 Olomouc, Czech Republic
- Institute
of Biophysics, Academy of Sciences of the Czech Republic, Kralovopolska
135, 612 65 Brno, Czech Republic
| |
Collapse
|
19
|
Mládek A, Krepl M, Svozil D, Čech P, Otyepka M, Banáš P, Zgarbová M, Jurečka P, Šponer J. Benchmark quantum-chemical calculations on a complete set of rotameric families of the DNA sugar–phosphate backbone and their comparison with modern density functional theory. Phys Chem Chem Phys 2013; 15:7295-310. [DOI: 10.1039/c3cp44383c] [Citation(s) in RCA: 28] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
20
|
Sharma P, Lait LA, Wetmore SD. Exploring the limits of nucleobase expansion: computational design of naphthohomologated (xx-) purines and comparison to the natural and xDNA purines. Phys Chem Chem Phys 2013; 15:15538-49. [DOI: 10.1039/c3cp52656a] [Citation(s) in RCA: 21] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/05/2023]
|
21
|
Šponer J, Mládek A, Šponer JE, Svozil D, Zgarbová M, Banáš P, Jurečka P, Otyepka M. The DNA and RNA sugar-phosphate backbone emerges as the key player. An overview of quantum-chemical, structural biology and simulation studies. Phys Chem Chem Phys 2012; 14:15257-77. [PMID: 23072945 DOI: 10.1039/c2cp41987d] [Citation(s) in RCA: 68] [Impact Index Per Article: 5.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Abstract
Knowledge of geometrical and physico-chemical properties of the sugar-phosphate backbone substantially contributes to the comprehension of the structural dynamics, function and evolution of nucleic acids. We provide a side by side overview of structural biology/bioinformatics, quantum chemical and molecular mechanical/simulation studies of the nucleic acids backbone. We highlight main features, advantages and limitations of these techniques, with a special emphasis given to their synergy. The present status of the research is then illustrated by selected examples which include classification of DNA and RNA backbone families, benchmark structure-energy quantum chemical calculations, parameterization of the dihedral space of simulation force fields, incorporation of arsenate into DNA, sugar-phosphate backbone self-cleavage in small RNA enzymes, and intricate geometries of the backbone in recurrent RNA building blocks. Although not apparent from the current literature showing limited overlaps between the QM, simulation and bioinformatics studies of the nucleic acids backbone, there in fact should be a major cooperative interaction between these three approaches in studies of the sugar-phosphate backbone.
Collapse
Affiliation(s)
- Jiří Šponer
- Institute of Biophysics, Academy Sciences of the Czech Republic, Královopolská 135, 612 65 Brno, Czech Republic.
| | | | | | | | | | | | | | | |
Collapse
|
22
|
Zgarbová M, Luque FJ, Šponer J, Otyepka M, Jurečka P. A Novel Approach for Deriving Force Field Torsion Angle Parameters Accounting for Conformation-Dependent Solvation Effects. J Chem Theory Comput 2012; 8:3232-42. [DOI: 10.1021/ct3001987] [Citation(s) in RCA: 26] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/04/2023]
Affiliation(s)
- Marie Zgarbová
- Regional Centre of Advanced
Technologies and Materials, Department of Physical Chemistry, Faculty
of Science, Palacky University, 17. listopadu 12, 77146 Olomouc, Czech
Republic
| | - F. Javier Luque
- Department de
Fisicoquímica
and Institut de Biomedicina (IBUB), Facultat de Farmàcia, Universitat
de Barcelona, Avgda Diagonal 643, Barcelona 08028, Spain
| | - Jiří Šponer
- Institute of Biophysics, Academy
of Sciences of the Czech Republic, Královopolská 135,
612 65 Brno, Czech Republic
- CEITEC - Central European Institute
of Technology, Masaryk University, Campus Bohunice, Kamenice 5, 625
00 Brno, Czech Republic
| | - Michal Otyepka
- Regional Centre of Advanced
Technologies and Materials, Department of Physical Chemistry, Faculty
of Science, Palacky University, 17. listopadu 12, 77146 Olomouc, Czech
Republic
| | - Petr Jurečka
- Regional Centre of Advanced
Technologies and Materials, Department of Physical Chemistry, Faculty
of Science, Palacky University, 17. listopadu 12, 77146 Olomouc, Czech
Republic
| |
Collapse
|
23
|
Beššeová I, Banáš P, Kührová P, Košinová P, Otyepka M, Šponer J. Simulations of A-RNA Duplexes. The Effect of Sequence, Solute Force Field, Water Model, and Salt Concentration. J Phys Chem B 2012; 116:9899-916. [DOI: 10.1021/jp3014817] [Citation(s) in RCA: 59] [Impact Index Per Article: 4.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Affiliation(s)
- Ivana Beššeová
- Institute
of Biophysics, Academy of Sciences of the Czech Republic, Kralovopolska
135, 612 65 Brno, Czech Republic
| | - Pavel Banáš
- Institute
of Biophysics, Academy of Sciences of the Czech Republic, Kralovopolska
135, 612 65 Brno, Czech Republic
- Regional Centre of Advanced
Technologies and Materials, Department of Physical Chemistry, Faculty
of Science, Palacky University, tr. 17
listopadu 12, 771 46, Olomouc, Czech Republic
| | - Petra Kührová
- Regional Centre of Advanced
Technologies and Materials, Department of Physical Chemistry, Faculty
of Science, Palacky University, tr. 17
listopadu 12, 771 46, Olomouc, Czech Republic
| | - Pavlína Košinová
- Regional Centre of Advanced
Technologies and Materials, Department of Physical Chemistry, Faculty
of Science, Palacky University, tr. 17
listopadu 12, 771 46, Olomouc, Czech Republic
| | - Michal Otyepka
- Institute
of Biophysics, Academy of Sciences of the Czech Republic, Kralovopolska
135, 612 65 Brno, Czech Republic
- Regional Centre of Advanced
Technologies and Materials, Department of Physical Chemistry, Faculty
of Science, Palacky University, tr. 17
listopadu 12, 771 46, Olomouc, Czech Republic
| | - Jiří Šponer
- Institute
of Biophysics, Academy of Sciences of the Czech Republic, Kralovopolska
135, 612 65 Brno, Czech Republic
- CEITEC - Central European Institute of Technology, Campus Bohunice, Kamenice
5, 625 00 Brno, Czech Republic
| |
Collapse
|
24
|
Krepl M, Zgarbová M, Stadlbauer P, Otyepka M, Banáš P, Koča J, Cheatham TE, Jurečka P, Šponer J. Reference simulations of noncanonical nucleic acids with different χ variants of the AMBER force field: quadruplex DNA, quadruplex RNA and Z-DNA. J Chem Theory Comput 2012; 8:2506-2520. [PMID: 23197943 PMCID: PMC3506181 DOI: 10.1021/ct300275s] [Citation(s) in RCA: 200] [Impact Index Per Article: 16.7] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Abstract
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.
Collapse
Affiliation(s)
- Miroslav Krepl
- Institute of Biophysics, Academy of Sciences of the Czech Republic, Kralovopolska 135, 612 65 Brno, Czech Republic
| | - Marie Zgarbová
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacky University, tr. 17 listopadu 12, 771 46, Olomouc, Czech Republic
| | - Petr Stadlbauer
- Institute of Biophysics, Academy of Sciences of the Czech Republic, Kralovopolska 135, 612 65 Brno, Czech Republic
| | - Michal Otyepka
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacky University, tr. 17 listopadu 12, 771 46, Olomouc, Czech Republic
| | - Pavel Banáš
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacky University, tr. 17 listopadu 12, 771 46, Olomouc, Czech Republic
| | - Jaroslav Koča
- CEITEC – Central European Institute of Technology, Campus Bohunice, Kamenice 5, 625 00 Brno, Czech Republic
- National Center for Biomolecular Research, Masaryk University, Campus Bohunice, Kamenice 5, 625 00 Brno, Czech Republic
| | - Thomas E. Cheatham
- Department of Medicinal Chemistry, College of Pharmacy, University of Utah, Salt Lake City, UT
| | - Petr Jurečka
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacky University, tr. 17 listopadu 12, 771 46, Olomouc, Czech Republic
| | - Jiří Šponer
- Institute of Biophysics, Academy of Sciences of the Czech Republic, Kralovopolska 135, 612 65 Brno, Czech Republic
- CEITEC – Central European Institute of Technology, Campus Bohunice, Kamenice 5, 625 00 Brno, Czech Republic
| |
Collapse
|
25
|
šponer J, Cang X, Cheatham TE. Molecular dynamics simulations of G-DNA and perspectives on the simulation of nucleic acid structures. Methods 2012; 57:25-39. [PMID: 22525788 PMCID: PMC3775459 DOI: 10.1016/j.ymeth.2012.04.005] [Citation(s) in RCA: 105] [Impact Index Per Article: 8.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/15/2012] [Revised: 04/04/2012] [Accepted: 04/06/2012] [Indexed: 11/29/2022] Open
Abstract
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.
Collapse
Affiliation(s)
- Jiří šponer
- Institute of Biophysics, Academy of Sciences of the Czech Republic, Královopolská 135, 612 65 Brno, Czech Republic
- CEITEC – Central European Institute of Technology, Masaryk University, Campus Bohunice, Kamenice 5, 625 00 Brno, Czech Republic
| | - Xiaohui Cang
- Drug Discovery and Design Center, State Key Laboratory of Drug Research, Shanghai Institute of Materia Medica, Chinese Academy of Sciences, 555 Zuchongzhi Road, Shanghai 201203, China
| | - Thomas E. Cheatham
- Department of Medicinal Chemistry, College of Pharmacy, Skaggs Hall 201, 2000 East 30 South, University of Utah, Salt Lake City, UT 84112, United States
| |
Collapse
|
26
|
Zgarbová M, Otyepka M, Šponer J, Mládek A, Banáš P, Cheatham TE, Jurečka P. Refinement of the Cornell et al. Nucleic Acids Force Field Based on Reference Quantum Chemical Calculations of Glycosidic Torsion Profiles. J Chem Theory Comput 2011; 7:2886-2902. [PMID: 21921995 PMCID: PMC3171997 DOI: 10.1021/ct200162x] [Citation(s) in RCA: 769] [Impact Index Per Article: 59.2] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/08/2011] [Indexed: 01/19/2023]
Abstract
We report a reparameterization of the glycosidic torsion χ of the Cornell et al. AMBER force field for RNA, χ(OL). The parameters remove destabilization of the anti region found in the ff99 force field and thus prevent formation of spurious ladder-like structural distortions in RNA simulations. They also improve the description of the syn region and the syn-anti balance as well as enhance MD simulations of various RNA structures. Although χ(OL) can be combined with both ff99 and ff99bsc0, we recommend the latter. We do not recommend using χ(OL) for B-DNA because it does not improve upon ff99bsc0 for canonical structures. However, it might be useful in simulations of DNA molecules containing syn nucleotides. Our parametrization is based on high-level QM calculations and differs from conventional parametrization approaches in that it incorporates some previously neglected solvation-related effects (which appear to be essential for obtaining correct anti/high-anti balance). Our χ(OL) force field is compared with several previous glycosidic torsion parametrizations.
Collapse
Affiliation(s)
- Marie Zgarbová
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacky University, 17 listopadu 12, 77146 Olomouc, Czech Republic
| | - Michal Otyepka
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacky University, 17 listopadu 12, 77146 Olomouc, Czech Republic
- Institute of Biophysics, Academy of Sciences of the Czech Republic, Královopolská 135, 612 65 Brno, Czech Republic
| | - Jiří Šponer
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacky University, 17 listopadu 12, 77146 Olomouc, Czech Republic
- Institute of Biophysics, Academy of Sciences of the Czech Republic, Královopolská 135, 612 65 Brno, Czech Republic
| | - Arnošt Mládek
- Institute of Biophysics, Academy of Sciences of the Czech Republic, Královopolská 135, 612 65 Brno, Czech Republic
| | - Pavel Banáš
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacky University, 17 listopadu 12, 77146 Olomouc, Czech Republic
- Institute of Biophysics, Academy of Sciences of the Czech Republic, Královopolská 135, 612 65 Brno, Czech Republic
| | - Thomas E. Cheatham
- Departments of Medicinal Chemistry and Pharmaceutics and Pharmaceutical Chemistry, University of Utah, 2000 East 30 South Skaggs 201, Salt Lake City, Utah 84112, United States
| | - Petr Jurečka
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacky University, 17 listopadu 12, 77146 Olomouc, Czech Republic
| |
Collapse
|
27
|
Yildirim I, Stern HA, Tubbs JD, Kennedy SD, Turner DH. Benchmarking AMBER force fields for RNA: comparisons to NMR spectra for single-stranded r(GACC) are improved by revised χ torsions. J Phys Chem B 2011; 115:9261-70. [PMID: 21721539 PMCID: PMC3140773 DOI: 10.1021/jp2016006] [Citation(s) in RCA: 82] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/17/2011] [Revised: 05/06/2011] [Indexed: 12/27/2022]
Abstract
Accurately modeling unpaired regions of RNA is important for predicting structure, dynamics, and thermodynamics of folded RNA. Comparisons between NMR data and molecular dynamics simulations provide a test of force fields used for modeling. Here, NMR spectroscopy, including NOESY, (1)H-(31)P HETCOR, DQF-COSY, and TOCSY, was used to determine conformational preferences for single-stranded GACC RNA. The spectra are consistent with a conformational ensemble containing major and minor A-form-like structures. In a series of 50 ns molecular dynamics (MD) simulations with the AMBER99 force field in explicit solvent, initial A-form-like structures rapidly evolve to disordered conformations. A set of 50 ns simulations with revised χ torsions (AMBER99χ force field) gives two primary conformations, consistent with the NMR spectra. A single 1.9 μs MD simulation with the AMBER99χ force field showed that the major and minor conformations are retained for almost 68% of the time in the first 700 ns, with multiple transformations from A-form to non-A-form conformations. For the rest of the simulation, random-coil structures and a stable non-A-form conformation inconsistent with NMR spectra were seen. Evidently, the AMBER99χ force field improves structural predictions for single-stranded GACC RNA compared to the AMBER99 force field, but further force field improvements are needed.
Collapse
Affiliation(s)
- Ilyas Yildirim
- Department of Chemistry and Center for RNA Biology, University of Rochester, Rochester, New York 14627, United States
| | - Harry A. Stern
- Department of Chemistry and Center for RNA Biology, University of Rochester, Rochester, New York 14627, United States
| | - Jason D. Tubbs
- Department of Chemistry and Center for RNA Biology, University of Rochester, Rochester, New York 14627, United States
| | - Scott D. Kennedy
- Department of Chemistry and Center for RNA Biology, University of Rochester, Rochester, New York 14627, United States
- Department of Biochemistry and Biophysics, School of Medicine and Dentistry, University of Rochester, Rochester, New York 14642, United States
| | - Douglas H. Turner
- Department of Chemistry and Center for RNA Biology, University of Rochester, Rochester, New York 14627, United States
| |
Collapse
|
28
|
Banáš P, Hollas D, Zgarbová M, Jurečka P, Orozco M, Cheatham TE, Šponer J, Otyepka M. Performance of Molecular Mechanics Force Fields for RNA Simulations: Stability of UUCG and GNRA Hairpins. J Chem Theory Comput 2010; 6:3836-3849. [DOI: 10.1021/ct100481h] [Citation(s) in RCA: 293] [Impact Index Per Article: 20.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
Affiliation(s)
- Pavel Banáš
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacky University Olomouc, tr. 17. listopadu 12, 771 46 Olomouc, Czech Republic, Institute of Biophysics, Academy of Sciences of the Czech Republic, Kralovopolska 135, 612 65 Brno, Czech Republic, Joint Research Program in Computational Biology, Institut de Recerca Biomédica and Barcelona Superocomputing Center, Baldiri i Reixac 10, Barcelona 08028, Spain, Jordi Girona 31, Barcelona 08028, Spain
| | - Daniel Hollas
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacky University Olomouc, tr. 17. listopadu 12, 771 46 Olomouc, Czech Republic, Institute of Biophysics, Academy of Sciences of the Czech Republic, Kralovopolska 135, 612 65 Brno, Czech Republic, Joint Research Program in Computational Biology, Institut de Recerca Biomédica and Barcelona Superocomputing Center, Baldiri i Reixac 10, Barcelona 08028, Spain, Jordi Girona 31, Barcelona 08028, Spain
| | - Marie Zgarbová
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacky University Olomouc, tr. 17. listopadu 12, 771 46 Olomouc, Czech Republic, Institute of Biophysics, Academy of Sciences of the Czech Republic, Kralovopolska 135, 612 65 Brno, Czech Republic, Joint Research Program in Computational Biology, Institut de Recerca Biomédica and Barcelona Superocomputing Center, Baldiri i Reixac 10, Barcelona 08028, Spain, Jordi Girona 31, Barcelona 08028, Spain
| | - Petr Jurečka
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacky University Olomouc, tr. 17. listopadu 12, 771 46 Olomouc, Czech Republic, Institute of Biophysics, Academy of Sciences of the Czech Republic, Kralovopolska 135, 612 65 Brno, Czech Republic, Joint Research Program in Computational Biology, Institut de Recerca Biomédica and Barcelona Superocomputing Center, Baldiri i Reixac 10, Barcelona 08028, Spain, Jordi Girona 31, Barcelona 08028, Spain
| | - Modesto Orozco
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacky University Olomouc, tr. 17. listopadu 12, 771 46 Olomouc, Czech Republic, Institute of Biophysics, Academy of Sciences of the Czech Republic, Kralovopolska 135, 612 65 Brno, Czech Republic, Joint Research Program in Computational Biology, Institut de Recerca Biomédica and Barcelona Superocomputing Center, Baldiri i Reixac 10, Barcelona 08028, Spain, Jordi Girona 31, Barcelona 08028, Spain
| | - Thomas E. Cheatham
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacky University Olomouc, tr. 17. listopadu 12, 771 46 Olomouc, Czech Republic, Institute of Biophysics, Academy of Sciences of the Czech Republic, Kralovopolska 135, 612 65 Brno, Czech Republic, Joint Research Program in Computational Biology, Institut de Recerca Biomédica and Barcelona Superocomputing Center, Baldiri i Reixac 10, Barcelona 08028, Spain, Jordi Girona 31, Barcelona 08028, Spain
| | - Jiří Šponer
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacky University Olomouc, tr. 17. listopadu 12, 771 46 Olomouc, Czech Republic, Institute of Biophysics, Academy of Sciences of the Czech Republic, Kralovopolska 135, 612 65 Brno, Czech Republic, Joint Research Program in Computational Biology, Institut de Recerca Biomédica and Barcelona Superocomputing Center, Baldiri i Reixac 10, Barcelona 08028, Spain, Jordi Girona 31, Barcelona 08028, Spain
| | - Michal Otyepka
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacky University Olomouc, tr. 17. listopadu 12, 771 46 Olomouc, Czech Republic, Institute of Biophysics, Academy of Sciences of the Czech Republic, Kralovopolska 135, 612 65 Brno, Czech Republic, Joint Research Program in Computational Biology, Institut de Recerca Biomédica and Barcelona Superocomputing Center, Baldiri i Reixac 10, Barcelona 08028, Spain, Jordi Girona 31, Barcelona 08028, Spain
| |
Collapse
|
29
|
Sano E, Li W, Yuki H, Liu X, Furihata T, Kobayashi K, Chiba K, Neya S, Hoshino T. Mechanism of the decrease in catalytic activity of human cytochrome P450 2C9 polymorphic variants investigated by computational analysis. J Comput Chem 2010; 31:2746-58. [DOI: 10.1002/jcc.21568] [Citation(s) in RCA: 31] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/08/2022]
|
30
|
Coletta A, Morozzo della Rocca B, Jaisankar P, Majumder HK, Chillemi G, Sanna N, Desideri A. Assignment of UV−vis Spectrum of (3,3′)-Diindolylmethane, a Leishmania donovani Topoisomerase IB Inhibitor and a Candidate DNA Minor Groove Binder. J Phys Chem A 2010; 114:7121-6. [PMID: 20550156 DOI: 10.1021/jp101494d] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Andrea Coletta
- Dipartimento di Biologia, Universitá degli Studi di Roma “Tor Vergata”, Via della Ricerca Scientifica, 00133 Roma, Italy, Departments of Medicinal Chemistry and Parasitology, Indian Institut of Chemical Biology, 4 Raja S.C. Mullick Road, Kolkata-700032, India, and CASPUR Interuniversities Consortium for Supercomputing Applications, Via dei Tizii 6b, 00185 Roma, Italy
| | - Blasco Morozzo della Rocca
- Dipartimento di Biologia, Universitá degli Studi di Roma “Tor Vergata”, Via della Ricerca Scientifica, 00133 Roma, Italy, Departments of Medicinal Chemistry and Parasitology, Indian Institut of Chemical Biology, 4 Raja S.C. Mullick Road, Kolkata-700032, India, and CASPUR Interuniversities Consortium for Supercomputing Applications, Via dei Tizii 6b, 00185 Roma, Italy
| | - Parasuraman Jaisankar
- Dipartimento di Biologia, Universitá degli Studi di Roma “Tor Vergata”, Via della Ricerca Scientifica, 00133 Roma, Italy, Departments of Medicinal Chemistry and Parasitology, Indian Institut of Chemical Biology, 4 Raja S.C. Mullick Road, Kolkata-700032, India, and CASPUR Interuniversities Consortium for Supercomputing Applications, Via dei Tizii 6b, 00185 Roma, Italy
| | - Hemanta K. Majumder
- Dipartimento di Biologia, Universitá degli Studi di Roma “Tor Vergata”, Via della Ricerca Scientifica, 00133 Roma, Italy, Departments of Medicinal Chemistry and Parasitology, Indian Institut of Chemical Biology, 4 Raja S.C. Mullick Road, Kolkata-700032, India, and CASPUR Interuniversities Consortium for Supercomputing Applications, Via dei Tizii 6b, 00185 Roma, Italy
| | - Giovanni Chillemi
- Dipartimento di Biologia, Universitá degli Studi di Roma “Tor Vergata”, Via della Ricerca Scientifica, 00133 Roma, Italy, Departments of Medicinal Chemistry and Parasitology, Indian Institut of Chemical Biology, 4 Raja S.C. Mullick Road, Kolkata-700032, India, and CASPUR Interuniversities Consortium for Supercomputing Applications, Via dei Tizii 6b, 00185 Roma, Italy
| | - Nico Sanna
- Dipartimento di Biologia, Universitá degli Studi di Roma “Tor Vergata”, Via della Ricerca Scientifica, 00133 Roma, Italy, Departments of Medicinal Chemistry and Parasitology, Indian Institut of Chemical Biology, 4 Raja S.C. Mullick Road, Kolkata-700032, India, and CASPUR Interuniversities Consortium for Supercomputing Applications, Via dei Tizii 6b, 00185 Roma, Italy
| | - Alessandro Desideri
- Dipartimento di Biologia, Universitá degli Studi di Roma “Tor Vergata”, Via della Ricerca Scientifica, 00133 Roma, Italy, Departments of Medicinal Chemistry and Parasitology, Indian Institut of Chemical Biology, 4 Raja S.C. Mullick Road, Kolkata-700032, India, and CASPUR Interuniversities Consortium for Supercomputing Applications, Via dei Tizii 6b, 00185 Roma, Italy
| |
Collapse
|
31
|
Mlýnský V, Banáš P, Hollas D, Réblová K, Walter NG, Šponer J, Otyepka M. Extensive molecular dynamics simulations showing that canonical G8 and protonated A38H+ forms are most consistent with crystal structures of hairpin ribozyme. J Phys Chem B 2010; 114:6642-52. [PMID: 20420375 PMCID: PMC2872159 DOI: 10.1021/jp1001258] [Citation(s) in RCA: 74] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/02/2023]
Abstract
The hairpin ribozyme is a prominent member of the group of small catalytic RNAs (RNA enzymes or ribozymes) because it does not require metal ions to achieve catalysis. Biochemical and structural data have implicated guanine 8 (G8) and adenine 38 (A38) as catalytic participants in cleavage and ligation catalyzed by the hairpin ribozyme, yet their exact role in catalysis remains disputed. To gain insight into dynamics in the active site of a minimal self-cleaving hairpin ribozyme, we have performed extensive classical, explicit-solvent molecular dynamics (MD) simulations on time scales of 50-150 ns. Starting from the available X-ray crystal structures, we investigated the structural impact of the protonation states of G8 and A38, and the inactivating A-1(2'-methoxy) substitution employed in crystallography. Our simulations reveal that a canonical G8 agrees well with the crystal structures while a deprotonated G8 profoundly distorts the active site. Thus MD simulations do not support a straightforward participation of the deprotonated G8 in catalysis. By comparison, the G8 enol tautomer is structurally well tolerated, causing only local rearrangements in the active site. Furthermore, a protonated A38H(+) is more consistent with the crystallography data than a canonical A38. The simulations thus support the notion that A38H(+) is the dominant form in the crystals, grown at pH 6. In most simulations, the canonical A38 departs from the scissile phosphate and substantially perturbs the structures of the active site and S-turn. Yet, we occasionally also observe formation of a stable A-1(2'-OH)...A38(N1) hydrogen bond, which documents the ability of the ribozyme to form this hydrogen bond, consistent with a potential role of A38 as general base catalyst. The presence of this hydrogen bond is, however, incompatible with the expected in-line attack angle necessary for self-cleavage, requiring a rapid transition of the deprotonated 2'-oxyanion to a position more favorable for in-line attack after proton transfer from A-1(2'-OH) to A38(N1). The simulations revealed a potential force field artifact, occasional but irreversible formation of "ladder-like", underwound A-RNA structure in one of the external helices. Although it does not affect the catalytic center of the hairpin ribozyme, further studies are under way to better assess possible influence of such force field behavior on long RNA simulations.
Collapse
Affiliation(s)
- Vojtěch Mlýnský
- Department of Physical Chemistry, Faculty of Science, Palacky University Olomouc, tr. 17. listopadu 12, 771 46 Olomouc, Czech Republic
| | - Pavel Banáš
- Department of Physical Chemistry, Faculty of Science, Palacky University Olomouc, tr. 17. listopadu 12, 771 46 Olomouc, Czech Republic
- Institute of Biophysics, Academy of Sciences of the Czech Republic, Kralovopolska 135, 612 65 Brno, Czech Republic
| | - Daniel Hollas
- Department of Physical Chemistry, Faculty of Science, Palacky University Olomouc, tr. 17. listopadu 12, 771 46 Olomouc, Czech Republic
| | - Kamila Réblová
- Institute of Biophysics, Academy of Sciences of the Czech Republic, Kralovopolska 135, 612 65 Brno, Czech Republic
| | - Nils G. Walter
- Department of Chemistry, Single Molecule Analysis Group, University of Michigan, 930 N. University Ave., Ann Arbor, MI 48109-1055, USA
| | - Jiří Šponer
- Department of Physical Chemistry, Faculty of Science, Palacky University Olomouc, tr. 17. listopadu 12, 771 46 Olomouc, Czech Republic
- Institute of Biophysics, Academy of Sciences of the Czech Republic, Kralovopolska 135, 612 65 Brno, Czech Republic
| | - Michal Otyepka
- Department of Physical Chemistry, Faculty of Science, Palacky University Olomouc, tr. 17. listopadu 12, 771 46 Olomouc, Czech Republic
- Institute of Biophysics, Academy of Sciences of the Czech Republic, Kralovopolska 135, 612 65 Brno, Czech Republic
| |
Collapse
|
32
|
Yildirim I, Stern HA, Kennedy SD, Tubbs JD, Turner DH. Reparameterization of RNA chi Torsion Parameters for the AMBER Force Field and Comparison to NMR Spectra for Cytidine and Uridine. J Chem Theory Comput 2010; 6:1520-1531. [PMID: 20463845 PMCID: PMC2867398 DOI: 10.1021/ct900604a] [Citation(s) in RCA: 136] [Impact Index Per Article: 9.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/17/2009] [Indexed: 11/29/2022]
Abstract
A reparameterization of the torsional parameters for the glycosidic dihedral angle, chi, for the AMBER99 force field in RNA nucleosides is used to provide a modified force field, AMBER99chi. Molecular dynamics simulations of cytidine, uridine, adenosine, and guanosine in aqueous solution using the AMBER99 and AMBER99chi force fields are compared with NMR results. For each nucleoside and force field, 10 individual molecular dynamics simulations of 30 ns each were run. For cytidine with AMBER99chi force field, each molecular dynamics simulation time was extended to 120 ns for convergence purposes. Nuclear magnetic resonance (NMR) spectroscopy, including one-dimensional (1D) (1)H, steady-state 1D (1)H nuclear Overhauser effect (NOE), and transient 1D (1)H NOE, was used to determine the sugar puckering and preferred base orientation with respect to the ribose of cytidine and uridine. The AMBER99 force field overestimates the population of syn conformations of the base orientation and of C2'-endo sugar puckering of the pyrimidines, while the AMBER99chi force field's predictions are more consistent with NMR results. Moreover, the AMBER99 force field prefers high anti conformations with glycosidic dihedral angles around 310 degrees for the base orientation of purines. The AMBER99chi force field prefers anti conformations around 185 degrees , which is more consistent with the quantum mechanical calculations and known 3D structures of folded ribonucleic acids (RNAs). Evidently, the AMBER99chi force field predicts the structural characteristics of ribonucleosides better than the AMBER99 force field and should improve structural and thermodynamic predictions of RNA structures.
Collapse
Affiliation(s)
- Ilyas Yildirim
- Department of Chemistry and Center for RNA Biology, University of Rochester, Rochester, New York 14627, and Department of Biochemistry and Biophysics, School of Medicine and Dentistry, University of Rochester, Rochester, New York 14642
| | | | | | | | | |
Collapse
|
33
|
Zgarbová M, Otyepka M, Šponer J, Hobza P, Jurečka P. Large-scale compensation of errors in pairwise-additive empirical force fields: comparison of AMBER intermolecular terms with rigorous DFT-SAPT calculations. Phys Chem Chem Phys 2010; 12:10476-93. [DOI: 10.1039/c002656e] [Citation(s) in RCA: 67] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
34
|
Fadrná E, Špačková N, Sarzyñska J, Koča J, Orozco M, Cheatham TE, Kulinski T, Šponer J. Single Stranded Loops of Quadruplex DNA As Key Benchmark for Testing Nucleic Acids Force Fields. J Chem Theory Comput 2009; 5:2514-30. [DOI: 10.1021/ct900200k] [Citation(s) in RCA: 111] [Impact Index Per Article: 7.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/30/2022]
Affiliation(s)
- Eva Fadrná
- National Centre for Biomolecular Research, Faculty of Science, Masaryk University, Kotlářská 2, 611 37 Brno, Czech Republic, Institute of Biophysics, Academy of Sciences of the Czech Republic, Královopolská 135, 612 65 Brno, Czech Republic, Institute of Bioorganic Chemistry, Polish Academy of Sciences, Noskowskiego 12/14, 61 704 Poznań, Poland, Joint IRB-BSC program on Computational Biology, Institute for Research in Biomedicine, Baldiri Reixac 10-12, 08028 Barcelona, Spain, Barcelona Supercomputing
| | - Nad’a Špačková
- National Centre for Biomolecular Research, Faculty of Science, Masaryk University, Kotlářská 2, 611 37 Brno, Czech Republic, Institute of Biophysics, Academy of Sciences of the Czech Republic, Královopolská 135, 612 65 Brno, Czech Republic, Institute of Bioorganic Chemistry, Polish Academy of Sciences, Noskowskiego 12/14, 61 704 Poznań, Poland, Joint IRB-BSC program on Computational Biology, Institute for Research in Biomedicine, Baldiri Reixac 10-12, 08028 Barcelona, Spain, Barcelona Supercomputing
| | - Joanna Sarzyñska
- National Centre for Biomolecular Research, Faculty of Science, Masaryk University, Kotlářská 2, 611 37 Brno, Czech Republic, Institute of Biophysics, Academy of Sciences of the Czech Republic, Královopolská 135, 612 65 Brno, Czech Republic, Institute of Bioorganic Chemistry, Polish Academy of Sciences, Noskowskiego 12/14, 61 704 Poznań, Poland, Joint IRB-BSC program on Computational Biology, Institute for Research in Biomedicine, Baldiri Reixac 10-12, 08028 Barcelona, Spain, Barcelona Supercomputing
| | - Jaroslav Koča
- National Centre for Biomolecular Research, Faculty of Science, Masaryk University, Kotlářská 2, 611 37 Brno, Czech Republic, Institute of Biophysics, Academy of Sciences of the Czech Republic, Královopolská 135, 612 65 Brno, Czech Republic, Institute of Bioorganic Chemistry, Polish Academy of Sciences, Noskowskiego 12/14, 61 704 Poznań, Poland, Joint IRB-BSC program on Computational Biology, Institute for Research in Biomedicine, Baldiri Reixac 10-12, 08028 Barcelona, Spain, Barcelona Supercomputing
| | - Modesto Orozco
- National Centre for Biomolecular Research, Faculty of Science, Masaryk University, Kotlářská 2, 611 37 Brno, Czech Republic, Institute of Biophysics, Academy of Sciences of the Czech Republic, Královopolská 135, 612 65 Brno, Czech Republic, Institute of Bioorganic Chemistry, Polish Academy of Sciences, Noskowskiego 12/14, 61 704 Poznań, Poland, Joint IRB-BSC program on Computational Biology, Institute for Research in Biomedicine, Baldiri Reixac 10-12, 08028 Barcelona, Spain, Barcelona Supercomputing
| | - Thomas E. Cheatham
- National Centre for Biomolecular Research, Faculty of Science, Masaryk University, Kotlářská 2, 611 37 Brno, Czech Republic, Institute of Biophysics, Academy of Sciences of the Czech Republic, Královopolská 135, 612 65 Brno, Czech Republic, Institute of Bioorganic Chemistry, Polish Academy of Sciences, Noskowskiego 12/14, 61 704 Poznań, Poland, Joint IRB-BSC program on Computational Biology, Institute for Research in Biomedicine, Baldiri Reixac 10-12, 08028 Barcelona, Spain, Barcelona Supercomputing
| | - Tadeusz Kulinski
- National Centre for Biomolecular Research, Faculty of Science, Masaryk University, Kotlářská 2, 611 37 Brno, Czech Republic, Institute of Biophysics, Academy of Sciences of the Czech Republic, Královopolská 135, 612 65 Brno, Czech Republic, Institute of Bioorganic Chemistry, Polish Academy of Sciences, Noskowskiego 12/14, 61 704 Poznań, Poland, Joint IRB-BSC program on Computational Biology, Institute for Research in Biomedicine, Baldiri Reixac 10-12, 08028 Barcelona, Spain, Barcelona Supercomputing
| | - Jiří Šponer
- National Centre for Biomolecular Research, Faculty of Science, Masaryk University, Kotlářská 2, 611 37 Brno, Czech Republic, Institute of Biophysics, Academy of Sciences of the Czech Republic, Královopolská 135, 612 65 Brno, Czech Republic, Institute of Bioorganic Chemistry, Polish Academy of Sciences, Noskowskiego 12/14, 61 704 Poznań, Poland, Joint IRB-BSC program on Computational Biology, Institute for Research in Biomedicine, Baldiri Reixac 10-12, 08028 Barcelona, Spain, Barcelona Supercomputing
| |
Collapse
|