QM Computations on Complete Nucleic Acids Building Blocks: Analysis of the Sarcin–Ricin RNA Motif Using DFT-D3, HF-3c, PM6-D3H, and MM Approaches
- 6 May 2014
- journal article
- research article
- Published by American Chemical Society (ACS) in Journal of Chemical Theory and Computation
- Vol. 10 (6) , 2615-2629
Abstract
A set of conformations obtained from explicit solvent molecular dynamics (MD) simulations of the Sarcin–Ricin internal loop (SRL) RNA motif is investigated using quantum mechanical (QM, TPSS-D3/def2-TZVP DFT-D3) and molecular mechanics (MM, AMBER parm99bsc0+χol3 force field) methods. Solvent effects are approximated using implicit solvent methods (COSMO for DFT-D3; GB and PB for MM). Large-scale DFT-D3 optimizations of the full 11-nucleotide motif are compared to MM results and reveal a higher flexibility of DFT-D3 over the MM in the optimization procedure. Conformational energies of the SRL motif expose significant differences in the DFT-D3 and MM energy descriptions that explain difficulties in MD simulations of the SRL motif. The TPSS-D3 data are in excellent agreement with results obtained by the hybrid functionals PW6B95-D3 and M06-2X. Computationally more efficient methods such as PM6-D3H and HF-3c show promising but partly inconsistent results. It is demonstrated that large-scale DFT-D3 computations on complete nucleic acids building blocks are a viable tool to complement the picture obtained from MD simulations and can be used as benchmarks for faster computational methods. Methodological challenges of large-scale QM computations on nucleic acids such as missing solvent–solute interactions and the truncation of the studied systems are discussed.Keywords
This publication has 97 references indexed in Scilit:
- Understanding the Sequence Preference of Recurrent RNA Building Blocks Using Quantum Chemistry: The Intrastrand RNA Dinucleotide PlatformJournal of Chemical Theory and Computation, 2011
- Refinement of the Cornell et al. Nucleic Acids Force Field Based on Reference Quantum Chemical Calculations of Glycosidic Torsion ProfilesJournal of Chemical Theory and Computation, 2011
- Assessing the Performance of the MM/PBSA and MM/GBSA Methods. 1. The Accuracy of Binding Free Energy Calculations Based on Molecular Dynamics SimulationsJournal of Chemical Information and Modeling, 2010
- Molecular Dynamics and Quantum Mechanics of RNA: Conformational and Chemical Change We Can Believe InAccounts of Chemical Research, 2009
- Dissection of the high rate constant for the binding of a ribotoxin to the ribosomeProceedings of the National Academy of Sciences, 2009
- Optimization of parameters for semiempirical methods V: Modification of NDDO approximations and application to 70 elementsJournal of Molecular Modeling, 2007
- UCSF Chimera—A visualization system for exploratory research and analysisJournal of Computational Chemistry, 2004
- VMD: Visual molecular dynamicsJournal of Molecular Graphics, 1996
- Molecular dynamics with coupling to an external bathThe Journal of Chemical Physics, 1984
- Comparison of simple potential functions for simulating liquid waterThe Journal of Chemical Physics, 1983