1
|
Brieuc F, Schran C, Uhl F, Forbert H, Marx D. Converged quantum simulations of reactive solutes in superfluid helium: The Bochum perspective. J Chem Phys 2020; 152:210901. [DOI: 10.1063/5.0008309] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/18/2023] Open
Affiliation(s)
- Fabien Brieuc
- Lehrstuhl für Theoretische Chemie, Ruhr-Universität Bochum, 44780 Bochum, Germany
| | - Christoph Schran
- Lehrstuhl für Theoretische Chemie, Ruhr-Universität Bochum, 44780 Bochum, Germany
| | - Felix Uhl
- Lehrstuhl für Theoretische Chemie, Ruhr-Universität Bochum, 44780 Bochum, Germany
| | - Harald Forbert
- Center for Solvation Science ZEMOS, Ruhr-Universität Bochum, 44780 Bochum, Germany
| | - Dominik Marx
- Lehrstuhl für Theoretische Chemie, Ruhr-Universität Bochum, 44780 Bochum, Germany
| |
Collapse
|
2
|
Rosales AR, Quinn TR, Wahlers J, Tomberg A, Zhang X, Helquist P, Wiest O, Norrby PO. Application of Q2MM to predictions in stereoselective synthesis. Chem Commun (Camb) 2018; 54:8294-8311. [PMID: 29971313 DOI: 10.1039/c8cc03695k] [Citation(s) in RCA: 31] [Impact Index Per Article: 5.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Abstract
Quantum-Guided Molecular Mechanics (Q2MM) can be used to derive transition state force fields (TSFFs) that allow the fast and accurate predictions of stereoselectivity for a wide range of catalytic enantioselective reactions. The basic ideas behind the derivation of TSFFs using Q2MM are discussed and the steps involved in obtaining a TSFF using the Q2MM code, publically available at github.com/q2mm, are shown. The applicability for a range of reactions, including several non-standard applications of Q2MM, is demonstrated. Future developments of the method are also discussed.
Collapse
Affiliation(s)
- Anthony R Rosales
- Department of Chemistry and Biochemistry, University of Notre Dame, Notre Dame, IN 46556, USA.
| | | | | | | | | | | | | | | |
Collapse
|
3
|
Okumura H, Higashi M, Yoshida Y, Sato H, Akiyama R. Theoretical approaches for dynamical ordering of biomolecular systems. Biochim Biophys Acta Gen Subj 2017; 1862:212-228. [PMID: 28988931 DOI: 10.1016/j.bbagen.2017.10.001] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/06/2017] [Revised: 09/30/2017] [Accepted: 10/04/2017] [Indexed: 01/21/2023]
Abstract
BACKGROUND Living systems are characterized by the dynamic assembly and disassembly of biomolecules. The dynamical ordering mechanism of these biomolecules has been investigated both experimentally and theoretically. The main theoretical approaches include quantum mechanical (QM) calculation, all-atom (AA) modeling, and coarse-grained (CG) modeling. The selected approach depends on the size of the target system (which differs among electrons, atoms, molecules, and molecular assemblies). These hierarchal approaches can be combined with molecular dynamics (MD) simulation and/or integral equation theories for liquids, which cover all size hierarchies. SCOPE OF REVIEW We review the framework of quantum mechanical/molecular mechanical (QM/MM) calculations, AA MD simulations, CG modeling, and integral equation theories. Applications of these methods to the dynamical ordering of biomolecular systems are also exemplified. MAJOR CONCLUSIONS The QM/MM calculation enables the study of chemical reactions. The AA MD simulation, which omits the QM calculation, can follow longer time-scale phenomena. By reducing the number of degrees of freedom and the computational cost, CG modeling can follow much longer time-scale phenomena than AA modeling. Integral equation theories for liquids elucidate the liquid structure, for example, whether the liquid follows a radial distribution function. GENERAL SIGNIFICANCE These theoretical approaches can analyze the dynamic behaviors of biomolecular systems. They also provide useful tools for exploring the dynamic ordering systems of biomolecules, such as self-assembly. This article is part of a Special Issue entitled "Biophysical Exploration of Dynamical Ordering of Biomolecular Systems" edited by Dr. Koichi Kato.
Collapse
Affiliation(s)
- Hisashi Okumura
- Research Center for Computational Science, Institute for Molecular Science, Okazaki, Aichi 444-8585, Japan; Department of Structural Molecular Science, The Graduate University for Advanced Studies, Okazaki, Aichi 444-8585, Japan.
| | - Masahiro Higashi
- Department of Chemistry, Biology and Marine Science, University of the Ryukyus, Nishihara, Okinawa 903-0213, Japan
| | - Yuichiro Yoshida
- Department of Molecular Engineering, Graduate School of Engineering, Kyoto University, Kyoto 615-8510, Japan
| | - Hirofumi Sato
- Department of Molecular Engineering, Graduate School of Engineering, Kyoto University, Kyoto 615-8510, Japan; Elements Strategy Initiative for Catalysts and Batteries, Kyoto University, Japan
| | - Ryo Akiyama
- Department of Chemistry, Kyushu University, Fukuoka 819-0395, Japan
| |
Collapse
|
4
|
Higashi M, Saito S. Quantitative Evaluation of Site Energies and Their Fluctuations of Pigments in the Fenna–Matthews–Olson Complex with an Efficient Method for Generating a Potential Energy Surface. J Chem Theory Comput 2016; 12:4128-37. [DOI: 10.1021/acs.jctc.6b00516] [Citation(s) in RCA: 29] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2023]
Affiliation(s)
- Masahiro Higashi
- Department
of Chemistry, Biology and Marine Science, University of the Ryukyus, 1 Senbaru, Nishihara, Okinawa 903-0213, Japan
| | - Shinji Saito
- Department
of Theoretical and Computational Molecular Science, Institute for Molecular Science, 38 Nishigo-Naka, Myodaiji, Okazaki, Aichi 444-8585, Japan
- The Graduate University for Advanced Studies (SOKENDAI), 38 Nishigo-Naka, Myodaiji, Okazaki, Aichi 444-8585, Japan
| |
Collapse
|
5
|
Higashi M, Truhlar DG. Combined Electrostatically Embedded Multiconfiguration Molecular Mechanics and Molecular Mechanical Method: Application to Molecular Dynamics Simulation of a Chemical Reaction in Aqueous Solution with Hybrid Density Functional Theory. J Chem Theory Comput 2015; 4:1032-9. [PMID: 26636356 DOI: 10.1021/ct8000816] [Citation(s) in RCA: 30] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Abstract
We here combine the electrostatically embedded multiconfiguration molecular mechanics (EE-MCMM) method for generating global potential energy surfaces in the presence of an electrostatic potential with molecular mechanics (MM). The resulting EE-MCMM/MM method is illustrated by applying it to carry out a molecular dynamics simulation for the symmetric bimolecular reaction Cl(-) + CH3Cl' → ClCH3 + Cl'(-) in aqueous solution with hybrid density functional theory as the quantum mechanical level. The potential of mean force is calculated, and the free energy barrier is found to be 25.3 kcal/mol, which is in good agreement with previous work. The advantage of the combined EE-MCMM and MM method is that the number of quantum mechanical calculations required for the active subsystem is very small compared to straight direct dynamics.
Collapse
Affiliation(s)
- Masahiro Higashi
- Department of Chemistry and Supercomputing Institute, 207 Pleasant Street SE, University of Minnesota, Minneapolis, Minnesota 55455-0431
| | - Donald G Truhlar
- Department of Chemistry and Supercomputing Institute, 207 Pleasant Street SE, University of Minnesota, Minneapolis, Minnesota 55455-0431
| |
Collapse
|
6
|
Xu X, Yu T, Papajak E, Truhlar DG. Multistructural variational transition state theory: kinetics of the hydrogen abstraction from carbon-2 of 2-methyl-1-propanol by hydroperoxyl radical including all structures and torsional anharmonicity. J Phys Chem A 2012; 116:10480-7. [PMID: 23020791 DOI: 10.1021/jp307504p] [Citation(s) in RCA: 21] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
We calculated the forward and reverse rate constants of the hydrogen abstraction reaction from carbon-2 of 2-methyl-1-propanol by hydroperoxyl radical over the temperature range 250-2400 K by using multistructural canonical variational transition state theory (MS-CVT) including both multiple-structure and torsional potential anharmonicity effects by the multistructural torsional anharmonicity (MS-T) method. In these calculations, multidimensional tunneling (MT) probabilities used to compute the tunneling transmission coefficients were evaluated by the small-curvature tunneling (SCT) approximation. Comparison with the rate constants obtained by the single-structural harmonic oscillator (SS-HO) approximation shows that multistructural anharmonicity increases the forward rate constants for all temperatures, but the reverse rate constants are reduced for temperatures lower than 430 K and increased for higher temperatures. The neglect of multistructural torsional anharmonicity would lead to errors of factors of 1.5, 8.8, and 13 at 300, 1000, and 2400 K, respectively, for the forward reaction, and would lead to errors of factors of 0.76, 3.0, and 6.0, respectively, at these temperatures for the reverse reaction.
Collapse
Affiliation(s)
- Xuefei Xu
- Department of Chemistry and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| | | | | | | |
Collapse
|
7
|
Xu X, Papajak E, Zheng J, Truhlar DG. Multi-structural variational transition state theory: kinetics of the 1,5-hydrogen shift isomerization of the 1-butoxyl radical including all structures and torsional anharmonicity. Phys Chem Chem Phys 2012; 14:4204-16. [DOI: 10.1039/c2cp23692c] [Citation(s) in RCA: 29] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
8
|
Tishchenko O, Truhlar DG. Gradient-based multiconfiguration Shepard interpolation for generating potential energy surfaces for polyatomic reactions. J Chem Phys 2010; 132:084109. [PMID: 20192292 DOI: 10.1063/1.3310296] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
This paper describes and illustrates a way to construct multidimensional representations of reactive potential energy surfaces (PESs) by a multiconfiguration Shepard interpolation (MCSI) method based only on gradient information, that is, without using any Hessian information from electronic structure calculations. MCSI, which is called multiconfiguration molecular mechanics (MCMM) in previous articles, is a semiautomated method designed for constructing full-dimensional PESs for subsequent dynamics calculations (classical trajectories, full quantum dynamics, or variational transition state theory with multidimensional tunneling). The MCSI method is based on Shepard interpolation of Taylor series expansions of the coupling term of a 2 x 2 electronically diabatic Hamiltonian matrix with the diagonal elements representing nonreactive analytical PESs for reactants and products. In contrast to the previously developed method, these expansions are truncated in the present version at the first order, and, therefore, no input of electronic structure Hessians is required. The accuracy of the interpolated energies is evaluated for two test reactions, namely, the reaction OH+H(2)-->H(2)O+H and the hydrogen atom abstraction from a model of alpha-tocopherol by methyl radical. The latter reaction involves 38 atoms and a 108-dimensional PES. The mean unsigned errors averaged over a wide range of representative nuclear configurations (corresponding to an energy range of 19.5 kcal/mol in the former case and 32 kcal/mol in the latter) are found to be within 1 kcal/mol for both reactions, based on 13 gradients in one case and 11 in the other. The gradient-based MCMM method can be applied for efficient representations of multidimensional PESs in cases where analytical electronic structure Hessians are too expensive or unavailable, and it provides new opportunities to employ high-level electronic structure calculations for dynamics at an affordable cost.
Collapse
Affiliation(s)
- Oksana Tishchenko
- Department of Chemistry and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, USA.
| | | |
Collapse
|
9
|
Zhang Y, Lin H. Quantum tunneling in testosterone 6beta-hydroxylation by cytochrome P450: reaction dynamics calculations employing multiconfiguration molecular-mechanical potential energy surfaces. J Phys Chem A 2010; 113:11501-8. [PMID: 19480428 DOI: 10.1021/jp901850c] [Citation(s) in RCA: 12] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022]
Abstract
Testosterone hydroxylation is a prototypical reaction of human cytochrome P450 3A4, which metabolizes about 50% of oral drugs on the market. Reaction dynamics calculations were carried out for the testosterone 6beta-hydrogen abstraction and the 6beta-d(1)-testosterone 6beta-duterium abstraction employing a model that consists of the substrate and the active oxidant compound I. The calculations were performed at the level of canonical variational transition state theory with multidimensional tunneling and were based on a semiglobal full-dimensional potential energy surface generated by the multiconfiguration molecular mechanics technique. The tunneling coefficients were found to be around 3, indicating substantial contributions by quantum tunneling. However, the tunneling made only modest contributions to the kinetic isotope effects. The kinetic isotope effects were computed to be about 2 in the doublet spin state and about 5 in the quartet spin state.
Collapse
Affiliation(s)
- Yan Zhang
- Chemistry Department, University of Colorado Denver, Denver, Colorado 80217-3364, USA
| | | |
Collapse
|
10
|
Tishchenko O, Truhlar DG. Non-Hermitian Multiconfiguration Molecular Mechanics. J Chem Theory Comput 2009; 5:1454-61. [DOI: 10.1021/ct900077g] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Affiliation(s)
- Oksana Tishchenko
- Department of Chemistry and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431
| | - Donald G. Truhlar
- Department of Chemistry and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431
| |
Collapse
|
11
|
Tishchenko O, Truhlar DG. Efficient global representations of potential energy functions: Trajectory calculations of bimolecular gas-phase reactions by multiconfiguration molecular mechanics. J Chem Phys 2009; 130:024105. [DOI: 10.1063/1.3042145] [Citation(s) in RCA: 13] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
12
|
Higashi M, Truhlar DG. Electrostatically Embedded Multiconfiguration Molecular Mechanics Based on the Combined Density Functional and Molecular Mechanical Method. J Chem Theory Comput 2008; 4:790-803. [DOI: 10.1021/ct800004y] [Citation(s) in RCA: 49] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/24/2023]
Affiliation(s)
- Masahiro Higashi
- Department of Chemistry and Supercomputing Institute, 207 Pleasant Street SE, University of Minnesota, Minneapolis, Minnesota 55455-0431
| | - Donald G. Truhlar
- Department of Chemistry and Supercomputing Institute, 207 Pleasant Street SE, University of Minnesota, Minneapolis, Minnesota 55455-0431
| |
Collapse
|
13
|
Albu TV, Espinosa-García J, Truhlar DG. Computational Chemistry of Polyatomic Reaction Kinetics and Dynamics: The Quest for an Accurate CH5 Potential Energy Surface. Chem Rev 2007; 107:5101-32. [DOI: 10.1021/cr078026x] [Citation(s) in RCA: 58] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Titus V. Albu
- Department of Chemistry, Box 5055, Tennessee Technological University, Cookeville, Tennessee 38505
| | | | - Donald G. Truhlar
- Department of Chemistry and Supercomputer Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431
| |
Collapse
|