1
|
Prasad VK, Otero-de-la-Roza A, DiLabio GA. Small-Basis Set Density-Functional Theory Methods Corrected with Atom-Centered Potentials. J Chem Theory Comput 2022; 18:2913-2930. [PMID: 35412817 DOI: 10.1021/acs.jctc.2c00036] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Abstract
Density functional theory (DFT) is currently the most popular method for modeling noncovalent interactions and thermochemistry. The accurate calculation of noncovalent interaction energies, reaction energies, and barrier heights requires choosing an appropriate functional and, typically, a relatively large basis set. Deficiencies of the density-functional approximation and the use of a limited basis set are the leading sources of error in the calculation of noncovalent and thermochemical properties in molecular systems. In this article, we present three new DFT methods based on the BLYP, M06-2X, and CAM-B3LYP functionals in combination with the 6-31G* basis set and corrected with atom-centered potentials (ACPs). ACPs are one-electron potentials that have the same form as effective-core potentials, except they do not replace any electrons. The ACPs developed in this work are used to generate energy corrections to the underlying DFT/basis-set method such that the errors in predicted chemical properties are minimized while maintaining the low computational cost of the parent methods. ACPs were developed for the elements H, B, C, N, O, F, Si, P, S, and Cl. The ACP parameters were determined using an extensive training set of 118655 data points, mostly of complete basis set coupled-cluster level quality. The target molecular properties for the ACP-corrected methods include noncovalent interaction energies, molecular conformational energies, reaction energies, barrier heights, and bond separation energies. The ACPs were tested first on the training set and then on a validation set of 42567 additional data points. We show that the ACP-corrected methods can predict the target molecular properties with accuracy close to complete basis set wavefunction theory methods, but at a computational cost of double-ζ DFT methods. This makes the new BLYP/6-31G*-ACP, M06-2X/6-31G*-ACP, and CAM-B3LYP/6-31G*-ACP methods uniquely suited to the calculation of noncovalent, thermochemical, and kinetic properties in large molecular systems.
Collapse
Affiliation(s)
- Viki Kumar Prasad
- Department of Chemistry, University of British Columbia, Okanagan, 3247 University Way, Kelowna, British Columbia V1V 1V7, Canada
| | - Alberto Otero-de-la-Roza
- Departamento de Química Física y Analítica, Facultad de Química, Universidad de Oviedo, MALTA Consolider Team, Oviedo E-33006, Spain
| | - Gino A DiLabio
- Department of Chemistry, University of British Columbia, Okanagan, 3247 University Way, Kelowna, British Columbia V1V 1V7, Canada
| |
Collapse
|
2
|
Prasad VK, Otero-de-la-Roza A, DiLabio GA. Fast and Accurate Quantum Mechanical Modeling of Large Molecular Systems Using Small Basis Set Hartree-Fock Methods Corrected with Atom-Centered Potentials. J Chem Theory Comput 2022; 18:2208-2232. [PMID: 35313106 DOI: 10.1021/acs.jctc.1c01128] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Abstract
There has been significant interest in developing fast and accurate quantum mechanical methods for modeling large molecular systems. In this work, by utilizing a machine learning regression technique, we have developed new low-cost quantum mechanical approaches to model large molecular systems. The developed approaches rely on using one-electron Gaussian-type functions called atom-centered potentials (ACPs) to correct for the basis set incompleteness and the lack of correlation effects in the underlying minimal or small basis set Hartree-Fock (HF) methods. In particular, ACPs are proposed for ten elements common in organic and bioorganic chemistry (H, B, C, N, O, F, Si, P, S, and Cl) and four different base methods: two minimal basis sets (MINIs and MINIX) plus a double-ζ basis set (6-31G*) in combination with dispersion-corrected HF (HF-D3/MINIs, HF-D3/MINIX, HF-D3/6-31G*) and the HF-3c method. The new ACPs are trained on a very large set (73 832 data points) of noncovalent properties (interaction and conformational energies) and validated additionally on a set of 32 048 data points. All reference data are of complete basis set coupled-cluster quality, mostly CCSD(T)/CBS. The proposed ACP-corrected methods are shown to give errors in the tenths of a kcal/mol range for noncovalent interaction energies and up to 2 kcal/mol for molecular conformational energies. More importantly, the average errors are similar in the training and validation sets, confirming the robustness and applicability of these methods outside the boundaries of the training set. In addition, the performance of the new ACP-corrected methods is similar to complete basis set density functional theory (DFT) but at a cost that is orders of magnitude lower, and the proposed ACPs can be used in any computational chemistry program that supports effective-core potentials without modification. It is also shown that ACPs improve the description of covalent and noncovalent bond geometries of the underlying methods and that the improvement brought about by the application of the ACPs is directly related to the number of atoms to which they are applied, allowing the treatment of systems containing some atoms for which ACPs are not available. Overall, the ACP-corrected methods proposed in this work constitute an alternative accurate, economical, and reliable quantum mechanical approach to describe the geometries, interaction energies, and conformational energies of systems with hundreds to thousands of atoms.
Collapse
Affiliation(s)
- Viki Kumar Prasad
- Department of Chemistry, University of British Columbia, Okanagan, 3247 University Way, Kelowna, British Columbia, Canada V1V 1V7
| | - Alberto Otero-de-la-Roza
- MALTA Consolider Team, Departamento de Química Física y Analítica, Facultad de Química, Universidad de Oviedo, E-33006 Oviedo, Spain
| | - Gino A DiLabio
- Department of Chemistry, University of British Columbia, Okanagan, 3247 University Way, Kelowna, British Columbia, Canada V1V 1V7
| |
Collapse
|
3
|
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
|
4
|
Gutten O, Jurečka P, Aliakbar Tehrani Z, Buděšínský M, Řezáč J, Rulíšek L. Conformational energies and equilibria of cyclic dinucleotides in vacuo and in solution: computational chemistry vs. NMR experiments. Phys Chem Chem Phys 2021; 23:7280-7294. [PMID: 33876088 DOI: 10.1039/d0cp05993e] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
Performance of computational methods in modelling cyclic dinucleotides - an important and challenging class of compounds - has been evaluated by two different benchmarks: (1) gas-phase conformational energies and (2) qualitative agreement with NMR observations of the orientation of the χ-dihedral angle in solvent. In gas-phase benchmarks, where CCSD(T) and DLPNO-CCSD(T) methods have been used as the reference, most of the (dispersion corrected) density functional approximations are accurate enough to justify prioritizing computational cost and compatibility with other modelling options as the criterion of choice. NMR experiments of 3'3'-c-di-AMP, 3'3'-c-GAMP, and 3'3'-c-di-GMP show the overall prevalence of the anti-conformation of purine bases, but some population of syn-conformations is observed for guanines. Implicit solvation models combined with quantum-chemical methods struggle to reproduce this behaviour, probably due to a lack of dynamics and explicitly modelled solvent, leading to structures that are too compact. Molecular dynamics simulations overrepresent the syn-conformation of guanine due to the overestimation of an intramolecular hydrogen bond. Our combination of experimental and computational benchmarks provides "error bars" for modelling cyclic dinucleotides in solvent, where such information is generally difficult to obtain, and should help gauge the interpretability of studies dealing with binding of cyclic dinucleotides to important pharmaceutical targets. At the same time, the presented analysis calls for improvement in both implicit solvation models and force-field parameters.
Collapse
Affiliation(s)
- Ondrej Gutten
- Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences, Flemingovo náměstí 2, 166 10, Praha 6, Czech Republic.
| | | | | | | | | | | |
Collapse
|
5
|
Mráziková K, Mlýnský V, Kührová P, Pokorná P, Kruse H, Krepl M, Otyepka M, Banáš P, Šponer J. UUCG RNA Tetraloop as a Formidable Force-Field Challenge for MD Simulations. J Chem Theory Comput 2020; 16:7601-7617. [PMID: 33215915 DOI: 10.1021/acs.jctc.0c00801] [Citation(s) in RCA: 16] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/23/2023]
Abstract
Explicit solvent atomistic molecular dynamics (MD) simulations represent an established technique to study structural dynamics of RNA molecules and an important complement for diverse experimental methods. However, performance of molecular mechanical (MM) force fields (ff's) remains far from satisfactory even after decades of development, as apparent from a problematic structural description of some important RNA motifs. Actually, some of the smallest RNA molecules belong to the most challenging systems for MD simulations and, among them, the UUCG tetraloop is saliently difficult. We report a detailed analysis of UUCG MD simulations, depicting the sequence of events leading to the loss of the UUCG native state during MD simulations. The total amount of MD simulation data analyzed in this work is close to 1.3 ms. We identify molecular interactions, backbone conformations, and substates that are involved in the process. Then, we unravel specific ff deficiencies using diverse quantum mechanical/molecular mechanical (QM/MM) and QM calculations. Comparison between the MM and QM methods shows discrepancies in the description of the 5'-flanking phosphate moiety and both signature sugar-base interactions. Our work indicates that poor behavior of the UUCG tetraloop in simulations is a complex issue that cannot be attributed to one dominant and straightforwardly correctable factor. Instead, there is a concerted effect of multiple ff inaccuracies that are coupled and amplifying each other. We attempted to improve the simulation behavior by some carefully tailored interventions, but the results were still far from satisfactory, underlying the difficulties in development of accurate nucleic acid ff's.
Collapse
Affiliation(s)
- Klaudia Mráziková
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 65 Brno, Czech Republic.,National Centre for Biomolecular Research, Faculty of Science, Masaryk University, Kamenice 5, 625 00 Brno, Czech Republic
| | - Vojtěch Mlýnský
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 65 Brno, Czech Republic
| | - Petra Kührová
- Regional Centre of Advanced Technologies and Materials, Faculty of Science, Palacký University, Šlechtitelů 27, 783 71 Olomouc, Czech Republic
| | - Pavlína Pokorná
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 65 Brno, Czech Republic.,National Centre for Biomolecular Research, Faculty of Science, Masaryk University, Kamenice 5, 625 00 Brno, Czech Republic
| | - Holger Kruse
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 65 Brno, Czech Republic
| | - Miroslav Krepl
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 65 Brno, Czech Republic.,Regional Centre of Advanced Technologies and Materials, Faculty of Science, Palacký University, Šlechtitelů 27, 783 71 Olomouc, Czech Republic
| | - Michal Otyepka
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 65 Brno, Czech Republic.,Regional Centre of Advanced Technologies and Materials, Faculty of Science, Palacký University, Šlechtitelů 27, 783 71 Olomouc, Czech Republic
| | - Pavel Banáš
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 65 Brno, Czech Republic.,Regional Centre of Advanced Technologies and Materials, Faculty of Science, Palacký University, Šlechtitelů 27, 783 71 Olomouc, Czech Republic
| | - Jiří Šponer
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 65 Brno, Czech Republic
| |
Collapse
|
6
|
Kumar N, Marx D. Deciphering the Self-Cleavage Reaction Mechanism of Hairpin Ribozyme. J Phys Chem B 2020; 124:4906-4918. [PMID: 32453954 DOI: 10.1021/acs.jpcb.0c03768] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Hairpin ribozyme catalyzes the reversible self-cleavage of phosphodiester bonds which plays prominent roles in key biological processes involving RNAs. Despite impressive advances on ribozymatic self-cleavage, critical aspects of its molecular reaction mechanism remain controversially debated. Here, we generate and analyze the multidimensional free energy landscape that underlies the reaction using extensive QM/MM metadynamics simulations to investigate in detail the full self-cleavage mechanism. This allows us to answer several pertinent yet controversial questions concerning activation of the 2'-OH group, the mechanistic role of water molecules present in the active site, and the full reaction pathway including the structures of transition states and intermediates. Importantly, we find that a sufficiently unrestricted reaction subspace must be mapped using accelerated sampling methods in order to compute the underlying free energy landscape. It is shown that lower-dimensional sampling where the bond formation and cleavage steps are coupled does not allow the system to sufficiently explore the landscape. On the basis of a three-dimensional free energy surface spanned by flexible generalized coordinates, we find that 2'-OH is indirectly activated by adjacent G8 nucleobase in conjunction with stabilizing H-bonding involving water. This allows the proton of the 2'-OH group to directly migrate toward the 5'-leaving group via a nonbridging oxygen of the phosphodiester link. At variance with similar enzymatic processes where water wires connected to protonable side chains of the protein matrix act as transient proton shuttles, no such de/reprotonation events of water molecules are found to be involved in this ribozymatic transesterification. Overall, our results support an acid-catalyzed reaction mechanism where A38 nucleobase directly acts as an acid whereas G8, in stark contrast, participates only indirectly via stabilizing the nascent nucleophile for subsequent attack. Moreover, we conclude that self-cleavage of hairpin ribozyme follows an AN + DN two-step associative pathway where the rate-determining step is the cleavage of the phosphodiester bond. These results provide a major advancement in our understanding of the unique catalytic mechanism of hairpin ribozyme which will fruitfully impact on the design of synthetic ribozymes.
Collapse
Affiliation(s)
- Narendra Kumar
- Lehrstuhl für Theoretische Chemie, Ruhr-Universität Bochum, 44780 Bochum, Germany
| | - Dominik Marx
- Lehrstuhl für Theoretische Chemie, Ruhr-Universität Bochum, 44780 Bochum, Germany
| |
Collapse
|
7
|
Řezáč J, Bím D, Gutten O, Rulíšek L. Toward Accurate Conformational Energies of Smaller Peptides and Medium-Sized Macrocycles: MPCONF196 Benchmark Energy Data Set. J Chem Theory Comput 2018; 14:1254-1266. [DOI: 10.1021/acs.jctc.7b01074] [Citation(s) in RCA: 55] [Impact Index Per Article: 9.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/23/2023]
Affiliation(s)
- Jan Řezáč
- Gilead Sciences Research Center and The Institute of Organic Chemistry and Biochemistry, Czech Academy of Sciences, Flemingovo náměstí 2, 166 10 Praha 6, Czech Republic
| | - Daniel Bím
- Gilead Sciences Research Center and The Institute of Organic Chemistry and Biochemistry, Czech Academy of Sciences, Flemingovo náměstí 2, 166 10 Praha 6, Czech Republic
| | - Ondrej Gutten
- Gilead Sciences Research Center and The Institute of Organic Chemistry and Biochemistry, Czech Academy of Sciences, Flemingovo náměstí 2, 166 10 Praha 6, Czech Republic
| | - Lubomír Rulíšek
- Gilead Sciences Research Center and The Institute of Organic Chemistry and Biochemistry, Czech Academy of Sciences, Flemingovo náměstí 2, 166 10 Praha 6, Czech Republic
| |
Collapse
|
8
|
Dzowo YK, Wolfbrandt C, Resendiz MJE, Wang H. Modeling of canonical and C2′- O-thiophenylmethyl modified hexamers of RNA. Insights into the nature of structural changes and thermal stability. NEW J CHEM 2018. [DOI: 10.1039/c8nj01739e] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Abstract
Modification of the C2′-O-position with thiophenylmethyl groups on both strands leads to thermal stabilization of the duplex. Predicting the effects that modifications will have on structure of RNA is of importance in the development of new RNA technologies.
Collapse
Affiliation(s)
| | | | | | - Haobin Wang
- Department of Chemistry
- University of Colorado Denver
- Denver
- USA
| |
Collapse
|
9
|
Harikrishna S, Pradeepkumar PI. Probing the Binding Interactions between Chemically Modified siRNAs and Human Argonaute 2 Using Microsecond Molecular Dynamics Simulations. J Chem Inf Model 2017; 57:883-896. [DOI: 10.1021/acs.jcim.6b00773] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022]
Affiliation(s)
- S. Harikrishna
- Department of Chemistry, Indian Institute of Technology Bombay, Mumbai−400076, India
| | - P. I. Pradeepkumar
- Department of Chemistry, Indian Institute of Technology Bombay, Mumbai−400076, India
| |
Collapse
|
10
|
Darré L, Ivani I, Dans PD, Gómez H, Hospital A, Orozco M. Small Details Matter: The 2′-Hydroxyl as a Conformational Switch in RNA. J Am Chem Soc 2016; 138:16355-16363. [DOI: 10.1021/jacs.6b09471] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Affiliation(s)
- Leonardo Darré
- Institute
for Research in Biomedicine (IRB Barcelona), The Barcelona Institute of Science and Technology, 08028 Barcelona, Spain
- Joint
BSC-IRB Program in Computational Biology, Institute for Research in Biomedicine, 08028 Barcelona, Spain
| | - Ivan Ivani
- Institute
for Research in Biomedicine (IRB Barcelona), The Barcelona Institute of Science and Technology, 08028 Barcelona, Spain
- Joint
BSC-IRB Program in Computational Biology, Institute for Research in Biomedicine, 08028 Barcelona, Spain
| | - Pablo D. Dans
- Institute
for Research in Biomedicine (IRB Barcelona), The Barcelona Institute of Science and Technology, 08028 Barcelona, Spain
- Joint
BSC-IRB Program in Computational Biology, Institute for Research in Biomedicine, 08028 Barcelona, Spain
| | - Hansel Gómez
- Institute
for Research in Biomedicine (IRB Barcelona), The Barcelona Institute of Science and Technology, 08028 Barcelona, Spain
- Joint
BSC-IRB Program in Computational Biology, Institute for Research in Biomedicine, 08028 Barcelona, Spain
| | - Adam Hospital
- Institute
for Research in Biomedicine (IRB Barcelona), The Barcelona Institute of Science and Technology, 08028 Barcelona, Spain
- Joint
BSC-IRB Program in Computational Biology, Institute for Research in Biomedicine, 08028 Barcelona, Spain
| | - Modesto Orozco
- Institute
for Research in Biomedicine (IRB Barcelona), The Barcelona Institute of Science and Technology, 08028 Barcelona, Spain
- Joint
BSC-IRB Program in Computational Biology, Institute for Research in Biomedicine, 08028 Barcelona, Spain
- Department
of Biochemistry and Biomedicine, Faculty of Biology, University of Barcelona, 08028 Barcelona, Spain
| |
Collapse
|
11
|
Dans PD, Walther J, Gómez H, Orozco M. Multiscale simulation of DNA. Curr Opin Struct Biol 2016; 37:29-45. [DOI: 10.1016/j.sbi.2015.11.011] [Citation(s) in RCA: 99] [Impact Index Per Article: 12.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/01/2015] [Revised: 11/23/2015] [Accepted: 11/25/2015] [Indexed: 01/05/2023]
|
12
|
Gkionis K, Kruse H, Šponer J. Derivation of Reliable Geometries in QM Calculations of DNA Structures: Explicit Solvent QM/MM and Restrained Implicit Solvent QM Optimizations of G-Quadruplexes. J Chem Theory Comput 2016; 12:2000-16. [DOI: 10.1021/acs.jctc.5b01025] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Affiliation(s)
- Konstantinos Gkionis
- Institute of Biophysics,
Academy of Sciences of the Czech Republic, Královopolská 135, 612
65 Brno, Czech Republic
| | - Holger Kruse
- Institute of Biophysics,
Academy of Sciences of the Czech Republic, Královopolská 135, 612
65 Brno, 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, Masaryk University, Campus Bohunice, Kamenice 5, 625 00 Brno, Czech Republic
| |
Collapse
|
13
|
Havrila M, Zgarbová M, Jurečka P, Banáš P, Krepl M, Otyepka M, Šponer J. Microsecond-Scale MD Simulations of HIV-1 DIS Kissing-Loop Complexes Predict Bulged-In Conformation of the Bulged Bases and Reveal Interesting Differences between Available Variants of the AMBER RNA Force Fields. J Phys Chem B 2015; 119:15176-90. [DOI: 10.1021/acs.jpcb.5b08876] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Marek Havrila
- Institute
of Biophysics, Academy of Sciences of the Czech Republic, Královopolská
135, 612 65 Brno, Czech Republic
| | - Marie Zgarbová
- Regional
Centre of Advanced Technologies and Materials, Department of Physical
Chemistry, Faculty of Science, Palacký University, tř.
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, Palacký University, tř.
17 listopadu 12, 771 46 Olomouc, Czech Republic
| | - Pavel Banáš
- Regional
Centre of Advanced Technologies and Materials, Department of Physical
Chemistry, Faculty of Science, Palacký University, tř.
17 listopadu 12, 771 46 Olomouc, Czech Republic
| | - Miroslav Krepl
- Institute
of Biophysics, Academy of Sciences of the Czech Republic, Královopolská
135, 612 65 Brno, Czech Republic
| | - Michal Otyepka
- Regional
Centre of Advanced Technologies and Materials, Department of Physical
Chemistry, Faculty of Science, Palacký University, tř.
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, Masaryk University, Campus Bohunice, Kamenice 5, 625 00 Brno, Czech Republic
| |
Collapse
|
14
|
Kruse H, Mladek A, Gkionis K, Hansen A, Grimme S, Sponer J. Quantum chemical benchmark study on 46 RNA backbone families using a dinucleotide unit. J Chem Theory Comput 2015; 11:4972-91. [PMID: 26574283 DOI: 10.1021/acs.jctc.5b00515] [Citation(s) in RCA: 80] [Impact Index Per Article: 8.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/21/2023]
Abstract
We have created a benchmark set of quantum chemical structure-energy data denoted as UpU46, which consists of 46 uracil dinucleotides (UpU), representing all known 46 RNA backbone conformational families. Penalty-function-based restrained optimizations with COSMO TPSS-D3/def2-TZVP ensure a balance between keeping the target conformation and geometry relaxation. The backbone geometries are close to the clustering-means of their respective RNA bioinformatics family classification. High-level wave function methods (DLPNO-CCSD(T) as reference) and a wide-range of dispersion-corrected or inclusive DFT methods (DFT-D3, VV10, LC-BOP-LRD, M06-2X, M11, and more) are used to evaluate the conformational energies. The results are compared to the Amber RNA bsc0χOL3 force field. Most dispersion-corrected DFT methods surpass the Amber force field significantly in accuracy and yield mean absolute deviations (MADs) for relative conformational energies of ∼0.4-0.6 kcal/mol. Double-hybrid density functionals represent the most accurate class of density functionals. Low-cost quantum chemical methods such as PM6-D3H+, HF-3c, DFTB3-D3, as well as small basis set calculations corrected for basis set superposition errors (BSSEs) by the gCP procedure are also tested. Unfortunately, the presently available low-cost methods are struggling to describe the UpU conformational energies with satisfactory accuracy. The UpU46 benchmark is an ideal test for benchmarking and development of fast methods to describe nucleic acids, including force fields.
Collapse
Affiliation(s)
- Holger Kruse
- 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
| | - Arnost Mladek
- Institute of Biophysics, Academy of Sciences of the Czech Republic , Královopolská 135, 612 65 Brno, Czech Republic
| | - Konstantinos Gkionis
- 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
| | - Andreas Hansen
- Mulliken Center for Theoretical Chemistry, Institut für Physikalische und Theoretische Chemie der Universität Bonn , Beringstr. 4, D-53115 Bonn, Germany
| | - Stefan Grimme
- Mulliken Center for Theoretical Chemistry, Institut für Physikalische und Theoretische Chemie der Universität Bonn , Beringstr. 4, D-53115 Bonn, Germany
| | - Jiri Sponer
- 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
|
Gresh N, Sponer JE, Devereux M, Gkionis K, de Courcy B, Piquemal JP, Sponer J. Stacked and H-Bonded Cytosine Dimers. Analysis of the Intermolecular Interaction Energies by Parallel Quantum Chemistry and Polarizable Molecular Mechanics. J Phys Chem B 2015; 119:9477-95. [DOI: 10.1021/acs.jpcb.5b01695] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/10/2023]
Affiliation(s)
- Nohad Gresh
- Chemistry & Biology, Nucleo(s)tides & Immunology for Therapy (CBNIT), CNRS UMR8601, Université Paris Descartes, PRES Sorbonne Paris Cité, UFR Biomédicale, 45 rue des Saints-Pères, 75270 Paris Cedex 06, France
- Laboratoire
de Chimie Théorique, Sorbonne Universités, UPMC, Paris 6, case courrier
137, 4, place Jussieu, Paris, F75252, France
- Laboratoire
de Chimie Théorique, UMR 7616 CNRS, case courrier 137, 4, place Jussieu, Paris, F75252, France
| | - Judit E. Sponer
- 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
| | - Mike Devereux
- Department
of Chemistry, University of Basel, Klingelbergstrasse 80, Basel CH 4056, Switzerland
| | - Konstantinos Gkionis
- Institute
of Biophysics, Academy of Sciences of the Czech Republic, Kralovopolska,
135, 612 65 Brno, Czech Republic
| | - Benoit de Courcy
- Chemistry & Biology, Nucleo(s)tides & Immunology for Therapy (CBNIT), CNRS UMR8601, Université Paris Descartes, PRES Sorbonne Paris Cité, UFR Biomédicale, 45 rue des Saints-Pères, 75270 Paris Cedex 06, France
- Laboratoire
de Chimie Théorique, Sorbonne Universités, UPMC, Paris 6, case courrier
137, 4, place Jussieu, Paris, F75252, France
- Laboratoire
de Chimie Théorique, UMR 7616 CNRS, case courrier 137, 4, place Jussieu, Paris, F75252, France
| | - Jean-Philip Piquemal
- Laboratoire
de Chimie Théorique, Sorbonne Universités, UPMC, Paris 6, case courrier
137, 4, place Jussieu, Paris, F75252, France
- Laboratoire
de Chimie Théorique, UMR 7616 CNRS, case courrier 137, 4, place Jussieu, Paris, F75252, France
| | - Jiri Sponer
- 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
|
16
|
Šponer JE, Šponer J, Giorgi A, Di Mauro E, Pino S, Costanzo G. Untemplated nonenzymatic polymerization of 3',5'cGMP: a plausible route to 3',5'-linked oligonucleotides in primordia. J Phys Chem B 2015; 119:2979-89. [PMID: 25625780 DOI: 10.1021/acs.jpcb.5b00601] [Citation(s) in RCA: 24] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Abstract
The high-energy 3',5' phosphodiester linkages conserved in 3',5' cyclic GMPs offer a genuine solution for monomer activation required by the transphosphorylation reactions that could lead to the emergence of the first simple oligonucleotide sequences on the early Earth. In this work we provide an in-depth characterization of the effect of the reaction conditions on the yield of the polymerization reaction of 3',5' cyclic GMPs both in aqueous environment as well as under dehydrating conditions. We show that the threshold temperature of the polymerization is about 30 °C lower under dehydrating conditions than in solution. In addition, we present a plausible exergonic reaction pathway for the polymerization reaction, which involves transient formation of anionic centers at the O3' positions of the participating riboses. We suggest that excess Na(+) cations inhibit the polymerization reaction because they block the anionic mechanism via neutralizing the negatively charged O3'. Our experimental findings are compatible with a prebiotic scenario, where gradual desiccation of the environment could induce polymerization of 3',5' cyclic GMPs synthesized in liquid.
Collapse
Affiliation(s)
- Judit E Šponer
- Institute of Biophysics, Academy of Sciences of the Czech Republic , Královopolská 135, CZ-61265 Brno, Czech Republic
| | | | | | | | | | | |
Collapse
|
17
|
Grimme S. A General Quantum Mechanically Derived Force Field (QMDFF) for Molecules and Condensed Phase Simulations. J Chem Theory Comput 2014; 10:4497-514. [DOI: 10.1021/ct500573f] [Citation(s) in RCA: 134] [Impact Index Per Article: 13.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Affiliation(s)
- Stefan Grimme
- Mulliken Center for Theoretical
Chemistry, Institut für Physikalische und Theoretische Chemie der Universität Bonn, Beringstr. 4, D-53115 Bonn, Germany
| |
Collapse
|
18
|
Wilson KA, Kellie JL, Wetmore SD. DNA-protein π-interactions in nature: abundance, structure, composition and strength of contacts between aromatic amino acids and DNA nucleobases or deoxyribose sugar. Nucleic Acids Res 2014; 42:6726-41. [PMID: 24744240 PMCID: PMC4041443 DOI: 10.1093/nar/gku269] [Citation(s) in RCA: 108] [Impact Index Per Article: 10.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022] Open
Abstract
Four hundred twenty-eight high-resolution DNA-protein complexes were chosen for a bioinformatics study. Although 164 crystal structures (38% of those searched) contained no interactions, 574 discrete π-contacts between the aromatic amino acids and the DNA nucleobases or deoxyribose were identified using strict criteria, including visual inspection. The abundance and structure of the interactions were determined by unequivocally classifying the contacts as either π-π stacking, π-π T-shaped or sugar-π contacts. Three hundred forty-four nucleobase-amino acid π-π contacts (60% of all interactions identified) were identified in 175 of the crystal structures searched. Unprecedented in the literature, 230 DNA-protein sugar-π contacts (40% of all interactions identified) were identified in 137 crystal structures, which involve C-H···π and/or lone-pair···π interactions, contain any amino acid and can be classified according to sugar atoms involved. Both π-π and sugar-π interactions display a range of relative monomer orientations and therefore interaction energies (up to -50 (-70) kJ mol(-1) for neutral (charged) interactions as determined using quantum chemical calculations). In general, DNA-protein π-interactions are more prevalent than perhaps currently accepted and the role of such interactions in many biological processes may yet to be uncovered.
Collapse
Affiliation(s)
- Katie A Wilson
- Department of Chemistry and Biochemistry, University of Lethbridge, 4401 University Drive West, Lethbridge, AB, T1K 3M4, Canada
| | - Jennifer L Kellie
- Department of Chemistry and Biochemistry, University of Lethbridge, 4401 University Drive West, Lethbridge, AB, T1K 3M4, Canada
| | - Stacey D Wetmore
- Department of Chemistry and Biochemistry, University of Lethbridge, 4401 University Drive West, Lethbridge, AB, T1K 3M4, Canada
| |
Collapse
|
19
|
Shea JE. Simulations Pave the Way for Exploring New Frontiers in the Biological Sciences. J Phys Chem Lett 2014; 5:1783-1784. [PMID: 26270383 DOI: 10.1021/jz5007934] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/04/2023]
Affiliation(s)
- Joan-Emma Shea
- Department of Chemistry and Biochemistry and Department of Physics, University of California, Santa Barbara, California 93106-9510, United States
| |
Collapse
|
20
|
Šponer J, Banáš P, Jurečka P, Zgarbová M, Kührová P, Havrila M, Krepl M, Stadlbauer P, Otyepka M. Molecular Dynamics Simulations of Nucleic Acids. From Tetranucleotides to the Ribosome. J Phys Chem Lett 2014; 5:1771-82. [PMID: 26270382 DOI: 10.1021/jz500557y] [Citation(s) in RCA: 113] [Impact Index Per Article: 11.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/16/2023]
Abstract
We present a brief overview of explicit solvent molecular dynamics (MD) simulations of nucleic acids. We explain physical chemistry limitations of the simulations, namely, the molecular mechanics (MM) force field (FF) approximation and limited time scale. Further, we discuss relations and differences between simulations and experiments, compare standard and enhanced sampling simulations, discuss the role of starting structures, comment on different versions of nucleic acid FFs, and relate MM computations with contemporary quantum chemistry. Despite its limitations, we show that MD is a powerful technique for studying the structural dynamics of nucleic acids with a fast growing potential that substantially complements experimental results and aids their interpretation.
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, Campus Bohunice, Kamenice 5, 625 00 Brno, Czech Republic
| | - Pavel Banáš
- §Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacký University, tř. 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, Palacký University, tř. 17 listopadu 12, 771 46 Olomouc, Czech Republic
| | - Marie Zgarbová
- §Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacký University, tř. 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, Palacký University, tř. 17 listopadu 12, 771 46 Olomouc, Czech Republic
| | - Marek Havrila
- †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
| | - Miroslav Krepl
- †Institute of Biophysics, Academy of Sciences of the Czech Republic, Královopolská 135, 612 65 Brno, Czech Republic
| | - Petr Stadlbauer
- †Institute of Biophysics, Academy of Sciences of the Czech Republic, Královopolská 135, 612 65 Brno, Czech Republic
| | - Michal Otyepka
- §Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacký University, tř. 17 listopadu 12, 771 46 Olomouc, Czech Republic
| |
Collapse
|
21
|
Kruse H, Havrila M, Šponer J. 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. J Chem Theory Comput 2014; 10:2615-29. [DOI: 10.1021/ct500183w] [Citation(s) in RCA: 30] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
Affiliation(s)
- Holger Kruse
- CEITEC
− Central European Institute of Technology, Campus Bohunice, Kamenice
5, 625 00 Brno, Czech Republic
| | - Marek Havrila
- CEITEC
− Central European Institute of Technology, Campus Bohunice, Kamenice
5, 625 00 Brno, Czech Republic
- Institute
of Biophysics, Academy of Sciences of the Czech Republic, Královopolská
135, 612 65 Brno, Czech Republic
| | - Jiřı́ Šponer
- CEITEC
− Central European Institute of Technology, Campus Bohunice, Kamenice
5, 625 00 Brno, Czech Republic
- Institute
of Biophysics, Academy of Sciences of the Czech Republic, Královopolská
135, 612 65 Brno, Czech Republic
| |
Collapse
|