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

Investor logo
Investor logo
Investor logo

Warning

This publication doesn't include Faculty of Arts. It includes Central European Institute of Technology. Official publication website can be found on muni.cz.
Authors

KRUSE Holger HAVRILA Marek ŠPONER Jiří

Year of publication 2014
Type Article in Periodical
Magazine / Source Journal of Chemical Theory and Computation
MU Faculty or unit

Central European Institute of Technology

Citation
Web http://pubs.acs.org/doi/abs/10.1021/ct500183w
Doi http://dx.doi.org/10.1021/ct500183w
Field Physical chemistry and theoretical chemistry
Keywords DENSITY-FUNCTIONAL THEORY; MOLECULAR-DYNAMICS SIMULATIONS; QUANTUM-CHEMICAL COMPUTATIONS; BASE-PHOSPHATE INTERACTIONS; AUXILIARY BASIS-SETS; PARTICLE MESH EWALD; AMBER FORCE-FIELD; NONCOVALENT INTERACTIONS; INTERACTION ENERGIES; STRUCTURAL BIOLOGY
Description 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+chi(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.
Related projects:

You are running an old browser version. We recommend updating your browser to its latest version.