Reparameterization of RNA χ Torsion Parameters for the AMBER Force Field and Comparison to NMR Spectra for Cytidine and Uridine
Open Access
- 16 April 2010
- journal article
- research article
- Published by American Chemical Society (ACS) in Journal of Chemical Theory and Computation
- Vol. 6 (5) , 1520-1531
- https://doi.org/10.1021/ct900604a
Abstract
A reparameterization of the torsional parameters for the glycosidic dihedral angle, χ, for the AMBER99 force field in RNA nucleosides is used to provide a modified force field, AMBER99χ. Molecular dynamics simulations of cytidine, uridine, adenosine, and guanosine in aqueous solution using the AMBER99 and AMBER99χ 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 AMBER99χ force field, each molecular dynamics simulation time was extended to 120 ns for convergence purposes. Nuclear magnetic resonance (NMR) spectroscopy, including one-dimensional (1D) 1H, steady-state 1D 1H nuclear Overhauser effect (NOE), and transient 1D 1H 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 AMBER99χ 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° for the base orientation of purines. The AMBER99χ force field prefers anti conformations around 185°, which is more consistent with the quantum mechanical calculations and known 3D structures of folded ribonucleic acids (RNAs). Evidently, the AMBER99χ force field predicts the structural characteristics of ribonucleosides better than the AMBER99 force field and should improve structural and thermodynamic predictions of RNA structures.Keywords
This publication has 62 references indexed in Scilit:
- Electrostatic interactions in the hairpin ribozyme account for the majority of the rate acceleration without chemical participation by nucleobasesRNA, 2008
- Quantum Mechanical/Molecular Mechanical Simulation Study of the Mechanism of Hairpin Ribozyme CatalysisJournal of the American Chemical Society, 2008
- Toward True DNA Base-Stacking Energies: MP2, CCSD(T), and Complete Basis Set CalculationsJournal of the American Chemical Society, 2002
- Global Minimum of the Adenine···Thymine Base Pair Corresponds Neither to Watson−Crick Nor to Hoogsteen Structures. Molecular Dynamic/Quenching/AMBER and ab Initio beyond Hartree−Fock StudiesJournal of the American Chemical Society, 2000
- A hybrid method for solutes in complex solvents: Density functional theory combined with empirical force fieldsThe Journal of Chemical Physics, 1999
- The GROMOS Biomolecular Simulation Program PackageThe Journal of Physical Chemistry A, 1999
- A Second Generation Force Field for the Simulation of Proteins, Nucleic Acids, and Organic MoleculesJournal of the American Chemical Society, 1995
- A Quantitative Model of the Escherichia coli 16 S RNA in the 30 S Ribosomal SubunitJournal of Molecular Biology, 1994
- Free energy perturbation method for chemical reactions in the condensed phase: a dynamic approach based on a combined quantum and molecular mechanics potentialJournal of the American Chemical Society, 1987
- CHARMM: A program for macromolecular energy, minimization, and dynamics calculationsJournal of Computational Chemistry, 1983