1
|
Brezina K, Beck H, Marsalek O. Reducing the Cost of Neural Network Potential Generation for Reactive Molecular Systems. J Chem Theory Comput 2023; 19:6589-6604. [PMID: 37747971 PMCID: PMC10569056 DOI: 10.1021/acs.jctc.3c00391] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/06/2023] [Indexed: 09/27/2023]
Abstract
Although machine learning potentials have recently had a substantial impact on molecular simulations, the construction of a robust training set can still become a limiting factor, especially due to the requirement of a reference ab initio simulation that covers all the relevant geometries of the system. Recognizing that this can be prohibitive for certain systems, we develop the method of transition tube sampling that mitigates the computational cost of training set and model generation. In this approach, we generate classical or quantum thermal geometries around a transition path describing a conformational change or a chemical reaction using only a sparse set of local normal mode expansions along this path and select from these geometries by an active learning protocol. This yields a training set with geometries that characterize the whole transition without the need for a costly reference trajectory. The performance of the method is evaluated on different molecular systems with the complexity of the potential energy landscape increasing from a single minimum to a double proton-transfer reaction with high barriers. Our results show that the method leads to training sets that give rise to models applicable in classical and path integral simulations alike that are on par with those based directly on ab initio calculations while providing the computational speedup we have come to expect from machine learning potentials.
Collapse
Affiliation(s)
- Krystof Brezina
- Charles University, Faculty of Mathematics
and Physics, Ke Karlovu
3, 121 16, Prague
2, Czech Republic
| | - Hubert Beck
- Charles University, Faculty of Mathematics
and Physics, Ke Karlovu
3, 121 16, Prague
2, Czech Republic
| | - Ondrej Marsalek
- Charles University, Faculty of Mathematics
and Physics, Ke Karlovu
3, 121 16, Prague
2, Czech Republic
| |
Collapse
|
2
|
Shukla S, Jakowski J, Kadian S, Narayan RJ. Computational approaches to delivery of anticancer drugs with multidimensional nanomaterials. Comput Struct Biotechnol J 2023; 21:4149-4158. [PMID: 37675288 PMCID: PMC10477808 DOI: 10.1016/j.csbj.2023.08.010] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/09/2023] [Revised: 08/14/2023] [Accepted: 08/14/2023] [Indexed: 09/08/2023] Open
Abstract
Functionalized nanotubes (NTs), nanosheets, nanorods, and porous organometallic scaffolds are potential in vivo carriers for cancer therapeutics. Precise delivery through these agents depends on factors like hydrophobicity, payload capacity, bulk/surface adsorption, orientation of molecules inside the host matrix, bonding, and nonbonding interactions. Herein, we summarize advances in simulation techniques, which are extremely valuable in initial geometry optimization and evaluation of the loading and unloading behavior of encapsulated drug molecules. Computational methods broadly involve the use of quantum and classical mechanics for studying the behavior of molecular properties. Combining theoretical processes with experimental techniques, such as X-ray crystallography, NMR spectroscopy, and bioassays, can provide a more comprehensive understanding of the structure and function of biological molecules. This integrated approach has led to numerous breakthroughs in drug discovery, enzyme design, and the study of complex biological processes. This short review provides an overview of results and challenges described from erstwhile investigations on the molecular interaction of anticancer drugs with nanocarriers of different aspect ratios.
Collapse
Affiliation(s)
- Shubhangi Shukla
- Joint Department of Biomedical Engineering, North Carolina State University, Raleigh, NC 27695-7907, United States
| | - Jacek Jakowski
- Center for Nanophase Materials Sciences, Oak Ridge National Laboratory, Oak Ridge, TN, United States
- Computational Sciences and Engineering Division, Oak Ridge National Laboratory, Oak Ridge, TN, United States
| | - Sachin Kadian
- Joint Department of Biomedical Engineering, North Carolina State University, Raleigh, NC 27695-7907, United States
| | - Roger J. Narayan
- Joint Department of Biomedical Engineering, North Carolina State University, Raleigh, NC 27695-7907, United States
| |
Collapse
|
3
|
Alibakhshi A, Steffen J, Pinilla C, Hartke B. Comparison of Implicit and Explicit Solvent Approaches in Ab Initio Evaluation of Thermochemistry in Solution: Application in Studying Boron Isotope Fractionation in Water. J Phys Chem A 2023; 127:2503-2510. [PMID: 36917555 DOI: 10.1021/acs.jpca.3c00154] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/16/2023]
Abstract
Evaluation of thermochemistry in solution plays a key role in numerous fields. For this task, the solvent effects are commonly included in theoretical computations based on either implicit or explicit solvent approaches. In the present study, we evaluate and compare the performance of some of the most widely applied methods based on these two approaches. For studying the solvent effect on thermochemistry with an explicit solvent, we demonstrate that partial normal mode analysis with frozen geometry of solvent molecules for multiple solute-solvent configurations can yield quite accurate and reliable results for a drastically reduced computational cost. As a case study, we consider the evaluation of the equilibrium constant for the boron isotope exchange between boric acid and borate (k3-4) in pure and saline water which is of high geochemical importance. Employing three different rigorous and high-precision theoretical approaches, we provide a reliable estimation of k3-4 which is a value within 1.028 to 1.030 for both pure and saline water which is in excellent agreement with experimental data.
Collapse
Affiliation(s)
- Amin Alibakhshi
- Theoretical Chemistry, Institute for Physical Chemistry, Christian-Albrechts-University, Olshausenstr. 40, 24118 Kiel, Germany
- Theoretical Chemistry, Ruhr-Universitaet Bochum, D-44780 Bochum, Germany
| | - Julien Steffen
- Theoretical Chemistry, Institute for Physical Chemistry, Christian-Albrechts-University, Olshausenstr. 40, 24118 Kiel, Germany
| | - Carlos Pinilla
- Departamento de Física y Geociencias, Universidad del Norte, Km 5 via Puerto Colombia, Barranquilla 080020, Colombia
- School of Chemistry, University of Bristol, Cantock's Close Road, BS8 1TS Bristol, U.K
| | - Bernd Hartke
- Theoretical Chemistry, Institute for Physical Chemistry, Christian-Albrechts-University, Olshausenstr. 40, 24118 Kiel, Germany
| |
Collapse
|
4
|
Korol R, Rosa-Raíces JL, Bou-Rabee N, Miller TF. Dimension-free path-integral molecular dynamics without preconditioning. J Chem Phys 2020; 152:104102. [PMID: 32171215 DOI: 10.1063/1.5134810] [Citation(s) in RCA: 16] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Convergence with respect to imaginary-time discretization (i.e., the number of ring-polymer beads) is an essential part of any path-integral-based molecular dynamics (MD) calculation. However, an unfortunate property of existing non-preconditioned numerical integration schemes for path-integral molecular dynamics-including essentially all existing ring-polymer molecular dynamics (RPMD) and thermostatted RPMD (T-RPMD) methods-is that for a given MD time step, the overlap between the exact ring-polymer Boltzmann-Gibbs distribution and that sampled using MD becomes zero in the infinite-bead limit. This has clear implications for hybrid Metropolis Monte Carlo/MD sampling schemes, and it also causes the divergence with bead number of the primitive path-integral kinetic-energy expectation value when using standard RPMD or T-RPMD. We show that these and other problems can be avoided through the introduction of "dimension-free" numerical integration schemes for which the sampled ring-polymer position distribution has non-zero overlap with the exact distribution in the infinite-bead limit for the case of a harmonic potential. Most notably, we introduce the BCOCB integration scheme, which achieves dimension freedom via a particular symmetric splitting of the integration time step and a novel implementation of the Cayley modification [R. Korol et al., J. Chem. Phys. 151, 124103 (2019)] for the free ring-polymer half-steps. More generally, we show that dimension freedom can be achieved via mollification of the forces from the external physical potential. The dimension-free path-integral numerical integration schemes introduced here yield finite error bounds for a given MD time step, even as the number of beads is taken to infinity; these conclusions are proven for the case of a harmonic potential and borne out numerically for anharmonic systems that include liquid water. The numerical results for BCOCB are particularly striking, allowing for nearly three-fold increases in the stable time step for liquid water with respect to the Bussi-Parrinello (OBABO) and Leimkuhler (BAOAB) integrators, while introducing negligible errors in the calculated statistical properties and absorption spectrum. Importantly, the dimension-free, non-preconditioned integration schemes introduced here preserve ergodicity and global second-order accuracy, and they remain simple, black-box methods that avoid additional computational costs, tunable parameters, or system-specific implementations.
Collapse
Affiliation(s)
- Roman Korol
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, USA
| | - Jorge L Rosa-Raíces
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, USA
| | - Nawaf Bou-Rabee
- Department of Mathematical Sciences, Rutgers University Camden, Camden, New Jersey 08102, USA
| | - Thomas F Miller
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, USA
| |
Collapse
|
5
|
Poltavsky I, Kapil V, Ceriotti M, Kim KS, Tkatchenko A. Accurate Description of Nuclear Quantum Effects with High-Order Perturbed Path Integrals (HOPPI). J Chem Theory Comput 2020; 16:1128-1135. [PMID: 31913625 DOI: 10.1021/acs.jctc.9b00881] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Imaginary time path-integral (PI) simulations that account for nuclear quantum effects (NQE) beyond the harmonic approximation are increasingly employed together with modern electronic-structure calculations. Existing PI methods are applicable to molecules, liquids, and solids; however, the computational cost of such simulations increases dramatically with decreasing temperature. To address this challenge, here, we propose to combine high-order PI factorization with perturbation theory (PT). Already for conventional second-order PI simulations, the PT ansatz increases the accuracy 2-fold compared to fourth-order schemes with the same settings. In turn, applying PT to high-order path integrals (HOPI) further improves the efficiency of simulations for molecular and condensed matter systems especially at low temperatures. We present results for bulk liquid water, the aspirin molecule, and the CH5+ molecule. Perturbed HOPI simulations remain both efficient and accurate down to 20 K and provide a convenient method to estimate the convergence of quantum-mechanical observables.
Collapse
Affiliation(s)
- Igor Poltavsky
- Physics and Materials Science Research Unit , University of Luxembourg , L-1511 Luxembourg City , Luxembourg
| | - Venkat Kapil
- Laboratory of Computational Science and Modelling, Institute of Materials , Ecole Polytechnique Fédérale de Lausanne , Lausanne , Switzerland
| | - Michele Ceriotti
- Laboratory of Computational Science and Modelling, Institute of Materials , Ecole Polytechnique Fédérale de Lausanne , Lausanne , Switzerland
| | - Kwang S Kim
- Department of Chemistry, School of Natural Science , Ulsan National Institute of Science and Technology , Ulsan 44919 , Korea
| | - Alexandre Tkatchenko
- Physics and Materials Science Research Unit , University of Luxembourg , L-1511 Luxembourg City , Luxembourg
| |
Collapse
|
6
|
Kannath S, Adamczyk P, Ferro-Costas D, Fernández-Ramos A, Major DT, Dybala-Defratyka A. Role of Microsolvation and Quantum Effects in the Accurate Prediction of Kinetic Isotope Effects: The Case of Hydrogen Atom Abstraction in Ethanol by Atomic Hydrogen in Aqueous Solution. J Chem Theory Comput 2020; 16:847-859. [PMID: 31904954 PMCID: PMC7588029 DOI: 10.1021/acs.jctc.9b00774] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
![]()
Hydrogen abstraction from ethanol
by atomic hydrogen in aqueous
solution is studied using two theoretical approaches: the multipath
variational transition state theory (MP-VTST) and a path-integral
formalism in combination with free-energy perturbation and umbrella
sampling (PI-FEP/UM). The performance of the models is compared to
experimental values of H kinetic isotope effects (KIE). Solvation
models used in this study ranged from purely implicit, via mixed–microsolvation
treated quantum mechanically via the density functional theory (DFT)
to fully explicit representation of the solvent, which was incorporated
using a combined quantum mechanical-molecular mechanical (QM/MM) potential.
The effects of the transition state conformation and the position
of microsolvating water molecules interacting with the solute on the
KIE are discussed. The KIEs are in good agreement with experiment
when MP-VTST is used together with a model that includes microsolvation
of the polar part of ethanol by five or six water molecules, emphasizing
the importance of explicit solvation in KIE calculations. Both, MP-VTST
and PI-FEP/UM enable detailed characterization of nuclear quantum
effects accompanying the hydrogen atom transfer reaction in aqueous
solution.
Collapse
Affiliation(s)
- Suraj Kannath
- Institute of Applied Radiation Chemistry, Faculty of Chemistry , Lodz University of Technology , Zeromskiego 116 , 90-924 Lodz , Poland
| | - Paweł Adamczyk
- Institute of Applied Radiation Chemistry, Faculty of Chemistry , Lodz University of Technology , Zeromskiego 116 , 90-924 Lodz , Poland
| | - David Ferro-Costas
- LAQV@REQUIMTE, Department of Chemistry & Biochemistry, Faculty of Sciences , University of Porto , Rua do Campo Alegre , 4169-007 Porto , Portugal.,Center for Research in Biological Chemistry and Molecular Materials (CIQUS) , Universidade de Santiago de Compostela , 15782 Santiago de Compostela , Spain
| | - Antonio Fernández-Ramos
- Center for Research in Biological Chemistry and Molecular Materials (CIQUS) , Universidade de Santiago de Compostela , 15782 Santiago de Compostela , Spain
| | - Dan Thomas Major
- Department of Chemistry and Institute for Nanotechnology & Advanced Materials , Bar-Ilan University , Ramat-Gan 52900 , Israel
| | - Agnieszka Dybala-Defratyka
- Institute of Applied Radiation Chemistry, Faculty of Chemistry , Lodz University of Technology , Zeromskiego 116 , 90-924 Lodz , Poland
| |
Collapse
|
7
|
Karandashev K, Vaníček J. Accelerating equilibrium isotope effect calculations. II. Stochastic implementation of direct estimators. J Chem Phys 2019; 151:134116. [PMID: 31594323 DOI: 10.1063/1.5124995] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Path integral calculations of equilibrium isotope effects and isotopic fractionation are expensive due to the presence of path integral discretization errors, statistical errors, and thermodynamic integration errors. Whereas the discretization errors can be reduced by high-order factorization of the path integral and statistical errors by using centroid virial estimators, two recent papers proposed alternative ways to completely remove the thermodynamic integration errors: Cheng and Ceriotti [J. Chem. Phys. 141, 244112 (2015)] employed a variant of free-energy perturbation called "direct estimators," while Karandashev and Vaníček [J. Chem. Phys. 143, 194104 (2017)] combined the thermodynamic integration with a stochastic change of mass and piecewise-linear umbrella biasing potential. Here, we combine the former approach with the stochastic change in mass in order to decrease its statistical errors when applied to larger isotope effects and perform a thorough comparison of different methods by computing isotope effects first on a harmonic model and then on methane and methanium, where we evaluate all isotope effects of the form CH4-xDx/CH4 and CH5-xDx +/CH5 +, respectively. We discuss the reasons for a surprising behavior of the original method of direct estimators, which performed well for a much larger range of isotope effects than what had been expected previously, as well as some implications of our work for the more general problem of free energy difference calculations.
Collapse
Affiliation(s)
- Konstantin Karandashev
- Laboratory of Theoretical Physical Chemistry, Institut des Sciences et Ingénierie Chimiques, Ecole Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| | - Jiří Vaníček
- Laboratory of Theoretical Physical Chemistry, Institut des Sciences et Ingénierie Chimiques, Ecole Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| |
Collapse
|
8
|
Liu M, Youmans KN, Gao J. Dual QM and MM Approach for Computing Equilibrium Isotope Fractionation Factor of Organic Species in Solution. Molecules 2018; 23:E2644. [PMID: 30326599 PMCID: PMC6222756 DOI: 10.3390/molecules23102644] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/17/2018] [Revised: 10/10/2018] [Accepted: 10/11/2018] [Indexed: 11/17/2022] Open
Abstract
A dual QM and MM approach for computing equilibrium isotope effects has been described. In the first partition, the potential energy surface is represented by a combined quantum mechanical and molecular mechanical (QM/MM) method, in which a solute molecule is treated quantum mechanically, and the remaining solvent molecules are approximated classically by molecular mechanics. In the second QM/MM partition, differential nuclear quantum effects responsible for the isotope effect are determined by a statistical mechanical double-averaging formalism, in which the nuclear centroid distribution is sampled classically by Newtonian molecular dynamics and the quantum mechanical spread of quantized particles about the centroid positions is treated using the path integral (PI) method. These partitions allow the potential energy surface to be properly represented such that the solute part is free of nuclear quantum effects for nuclear quantum mechanical simulations, and the double-averaging approach has the advantage of sampling efficiency for solvent configuration and for path integral convergence. Importantly, computational precision is achieved through free energy perturbation (FEP) theory to alchemically mutate one isotope into another. The PI-FEP approach is applied to model systems for the 18O enrichment found in cellulose of trees to determine the isotope enrichment factor of carbonyl compounds in water. The present method may be useful as a general tool for studying isotope fractionation in biological and geochemical systems.
Collapse
Affiliation(s)
- Meiyi Liu
- Laboratory of Theoretical and Computational Chemistry, Theoretical Chemistry Institute, Jilin University, Changchun 130023, China.
| | - Katelyn N Youmans
- Department of Chemistry and Supercomputing Institute, University of Minnesota, Minneapolis, MN 55455, USA.
| | - Jiali Gao
- Laboratory of Theoretical and Computational Chemistry, Theoretical Chemistry Institute, Jilin University, Changchun 130023, China.
- Department of Chemistry and Supercomputing Institute, University of Minnesota, Minneapolis, MN 55455, USA.
| |
Collapse
|
9
|
Vasquez L, Rostkowski M, Gelman F, Dybala-Defratyka A. Can Path Integral Molecular Dynamics Make a Good Approximation for Vapor Pressure Isotope Effects Prediction for Organic Solvents? A Comparison to ONIOM QM/MM and QM Cluster Calculation. J Phys Chem B 2018; 122:7353-7364. [PMID: 29961315 DOI: 10.1021/acs.jpcb.8b03444] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/17/2023]
Abstract
Isotopic fractionation of volatile organic compounds (VOCs), which are under strict measures of control because of their potential harm to the environment and humans, has an important ecological aspect, as the isotopic composition of compounds may depend on the conditions in which such compounds are distributed in Nature. Therefore, detailed knowledge on isotopic fractionation, not only experimental but also based on theoretical models, is crucial to follow conditions and pathways within which these contaminants are spread throughout the ecosystems. In this work, we present carbon and, for the first time, bromine vapor pressure isotope effect (VPIE) on the evaporation process from pure-phase systems-dibromomethane and bromobenzene, the representatives of aliphatic and aromatic brominated VOCs. We combine isotope effects measurements with their theoretical prediction using three computational techniques, namely path integral molecular dynamics, QM cluster, and hybrid ONIOM models. While evaporation of both compounds resulted in normal bromine VPIEs, the difference in the direction of carbon isotopic fractionation is observed for the aliphatic and aromatic compounds, where VPIEs are inverse and normal, respectively. Even though theoretical models tested here turned out to be insufficient for quantitative agreement with the experimental values, cluster electronic structure calculations, as well as two-layer ONIOM computations, provided better reproduction of experimental trends.
Collapse
Affiliation(s)
- Luis Vasquez
- Institute of Applied Radiation Chemistry, Faculty of Chemistry , Lodz University of Technology , Żeromskiego 116 , 90-924 Łódź , Poland
| | - Michal Rostkowski
- Institute of Applied Radiation Chemistry, Faculty of Chemistry , Lodz University of Technology , Żeromskiego 116 , 90-924 Łódź , Poland
| | - Faina Gelman
- Geological Survey of Israel , Malkhei Israel Street 30 , 95501 Jerusalem , Israel
| | - Agnieszka Dybala-Defratyka
- Institute of Applied Radiation Chemistry, Faculty of Chemistry , Lodz University of Technology , Żeromskiego 116 , 90-924 Łódź , Poland
| |
Collapse
|
10
|
Abstract
Ab initio molecular dynamics is an irreplaceable technique for the realistic simulation of complex molecular systems and processes from first principles. This paper proposes a comprehensive and self-contained review of ab initio molecular dynamics from a computational perspective and from first principles. Quantum mechanics is presented from a molecular dynamics perspective. Various approximations and formulations are proposed, including the Ehrenfest, Born–Oppenheimer, and Hartree–Fock molecular dynamics. Subsequently, the Kohn–Sham formulation of molecular dynamics is introduced as well as the afferent concept of density functional. As a result, Car–Parrinello molecular dynamics is discussed, together with its extension to isothermal and isobaric processes. Car–Parrinello molecular dynamics is then reformulated in terms of path integrals. Finally, some implementation issues are analysed, namely, the pseudopotential, the orbital functional basis, and hybrid molecular dynamics.
Collapse
|
11
|
Poltavsky I, DiStasio RA, Tkatchenko A. Perturbed path integrals in imaginary time: Efficiently modeling nuclear quantum effects in molecules and materials. J Chem Phys 2018; 148:102325. [DOI: 10.1063/1.5006596] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/01/2023] Open
Affiliation(s)
- Igor Poltavsky
- Physics and Materials Science Research Unit, University of Luxembourg, Luxembourg L-1511, Luxembourg
| | - Robert A. DiStasio
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853, USA
| | - Alexandre Tkatchenko
- Physics and Materials Science Research Unit, University of Luxembourg, Luxembourg L-1511, Luxembourg
| |
Collapse
|
12
|
|
13
|
Karandashev K, Xu ZH, Meuwly M, Vaníček J, Richardson JO. Kinetic isotope effects and how to describe them. STRUCTURAL DYNAMICS (MELVILLE, N.Y.) 2017; 4:061501. [PMID: 29282447 PMCID: PMC5729036 DOI: 10.1063/1.4996339] [Citation(s) in RCA: 27] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/16/2017] [Accepted: 08/23/2017] [Indexed: 06/01/2023]
Abstract
We review several methods for computing kinetic isotope effects in chemical reactions including semiclassical and quantum instanton theory. These methods describe both the quantization of vibrational modes as well as tunneling and are applied to the ⋅H + H2 and ⋅H + CH4 reactions. The absolute rate constants computed with the semiclassical instanton method both using on-the-fly electronic structure calculations and fitted potential-energy surfaces are also compared directly with exact quantum dynamics results. The error inherent in the instanton approximation is found to be relatively small and similar in magnitude to that introduced by using fitted surfaces. The kinetic isotope effect computed by the quantum instanton is even more accurate, and although it is computationally more expensive, the efficiency can be improved by path-integral acceleration techniques. We also test a simple approach for designing potential-energy surfaces for the example of proton transfer in malonaldehyde. The tunneling splittings are computed, and although they are found to deviate from experimental results, the ratio of the splitting to that of an isotopically substituted form is in much better agreement. We discuss the strengths and limitations of the potential-energy surface and based on our findings suggest ways in which it can be improved.
Collapse
Affiliation(s)
- Konstantin Karandashev
- Laboratory of Theoretical Physical Chemistry, Institut des Sciences et Ingénierie Chimiques, Ecole Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| | - Zhen-Hao Xu
- Department of Chemistry, University of Basel, Klingelbergstrasse 80, CH-4056 Basel, Switzerland
| | - Markus Meuwly
- Department of Chemistry, University of Basel, Klingelbergstrasse 80, CH-4056 Basel, Switzerland
| | - Jiří Vaníček
- Laboratory of Theoretical Physical Chemistry, Institut des Sciences et Ingénierie Chimiques, Ecole Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| | - Jeremy O Richardson
- Laboratory of Physical Chemistry, Department of Chemistry and Applied Biosciences, Eidgenössische Technische Hochschule Zürich (ETHZ), CH-8093 Zürich, Switzerland
| |
Collapse
|
14
|
Buchowiecki M, Prątnicki F. High temperature partition function – a key role of ro-vibrational coupling and inflection points. Mol Phys 2017. [DOI: 10.1080/00268976.2017.1334132] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/19/2022]
Affiliation(s)
- Marcin Buchowiecki
- Institute of Physics, Department of Mathematics and Physics, University of Szczecin, Szczecin, Poland
| | - Filip Prątnicki
- Institute of Physics, Department of Mathematics and Physics, University of Szczecin, Szczecin, Poland
| |
Collapse
|
15
|
Dupuis R, Benoit M, Tuckerman ME, Méheut M. Importance of a Fully Anharmonic Treatment of Equilibrium Isotope Fractionation Properties of Dissolved Ionic Species As Evidenced by Li +(aq). Acc Chem Res 2017. [PMID: 28644616 DOI: 10.1021/acs.accounts.6b00607] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Equilibrium fractionation of stable isotopes is critically important in fields ranging from chemistry, including medicinal chemistry, electrochemistry, geochemistry, and nuclear chemistry, to environmental science. The dearth of reliable estimates of equilibrium fractionation factors, from experiment or from natural observations, has created a need for accurate computational approaches. Because isotope fractionation is a purely quantum mechanical phenomenon, exact calculation of fractionation factors is nontrivial. Consequently, a severe approximation is often made, in which it is assumed that the system can be decomposed into a set of independent harmonic oscillators. Reliance on this often crude approximation is one of the primary reasons that theoretical prediction of isotope fractionation has lagged behind experiment. A class of problems for which one might expect the harmonic approximation to perform most poorly is the isotopic fractionation between solid and solution phases. In order to illustrate the errors associated with the harmonic approximation, we have considered the fractionation of Li isotopes between aqueous solution and phyllosilicate minerals, where we find that the harmonic approximation overestimates isotope fractionation factors by as much as 30% at 25 °C. Lithium is a particularly interesting species to examine, as natural lithium isotope signatures provide information about hydrothermal processes, carbon cycle, and regulation of the Earth's climate by continental alteration. Further, separation of lithium isotopes is of growing interest in the nuclear industry due to a need for pure 6Li and 7Li isotopes. Moving beyond the harmonic approximation entails performing exact quantum calculations, which can be achieved using the Feynman path integral formulation of quantum statistical mechanics. In the path integral approach, a system of quantum particles is represented as a set of classical-like ring-polymer chains, whose interparticle interactions are determined by the rules of quantum mechanics. Because a classical isomorphism exists between the true quantum system and the system of ring-polymers, classical-like methods can be applied. Recent developments of efficient path integral approaches for the exact calculation of isotope fractionation now allow the case of the aforementioned dissolved Li fractionation properties to be studied in detail. Applying this technique, we find that the calculations yield results that are in good agreement with both experimental data and natural observations. Importantly, path integral methods, being fully atomistic, allow us to identify the origins of anharmonic effects and to make reliable predictions at temperatures that are experimentally inaccessible yet are, nevertheless, relevant for natural phenomena.
Collapse
Affiliation(s)
- Romain Dupuis
- DIPC, Paseo Manuel de Lardizabal, 4, 20018 San Sebastiàn, Spain
| | - Magali Benoit
- CEMES CNRS UPR8011, 29 rue Jeanne Marvig, 31055 Cedex Toulouse, France
| | - Mark E. Tuckerman
- Department of Chemistry, New York University, New York, New York 10003, United States
- Courant Institute of Mathematical Sciences, New York University, New York, New York 10012, United States
- NYU-ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai 200062, China
| | - Merlin Méheut
- GET, CNRS UMR 5563, IRD UR 154, Université Toulouse III, OMP, 14 avenue Edouard Belin, 31400 Toulouse, France
| |
Collapse
|
16
|
Dupuis R, Dolado JS, Benoit M, Surga J, Ayuela A. Quantum Nuclear Dynamics of Protons within Layered Hydroxides at High Pressure. Sci Rep 2017; 7:4842. [PMID: 28687787 PMCID: PMC5501843 DOI: 10.1038/s41598-017-04080-2] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/03/2017] [Accepted: 05/08/2017] [Indexed: 11/09/2022] Open
Abstract
Studies of the structure of hydroxides under pressure using neutron diffraction reveal that the high concentration of hydrogen is distributed in a disordered network. The disorder in the hydrogen-bond network and possible phase transitions are reported to occur at pressures within the range accessible to experiments for layered calcium hydroxides, which are considered to be exemplary prototype materials. In this study, the static and dynamical properties of these layered hydroxides are investigated using a quantum approach describing nuclear motion, shown herein to be required particularly when studying diffusion processes involving light hydrogen atoms. The effect of high-pressure on the disordered hydrogen-bond network shows that the protons tunnel back and forth across the barriers between three potential minima around the oxygen atoms. At higher pressures the structure has quasi two-dimensional layers of hydrogen atoms, such that at low temperatures this causes the barrier crossing of the hydrogen to be significantly rarefied. Furthermore, for moderate values of both temperature and pressure this process occurs less often than the usual mechanism of proton transport via vacancies, limiting global proton diffusion within layers at high pressure.
Collapse
Affiliation(s)
- Romain Dupuis
- DIPC, Paseo Manuel de Lardizabal, 4, 20018, San Sebastian, Spain.
| | - Jorge S Dolado
- TECNALIA, Parque Científico y Tecnológico de Bizkaia 700 Edificio, 48160, Elexalde Derio, Spain
| | - Magali Benoit
- CEMES, 29 Rue Jeanne Marvig, 31055, Toulouse Cedex 4, France
| | - Jose Surga
- INTEVEP, Urb. Santa Rosa, Los Teques, 1201, Venezuela
| | - Andrés Ayuela
- DIPC, Paseo Manuel de Lardizabal, 4, 20018, San Sebastian, Spain.,CFM-MPC CSIC-UPV/EHU, Paseo Manuel de Lardizabal, 5, 20018, San Sebastian, Spain
| |
Collapse
|
17
|
Unravelling the influence of quantum proton delocalization on electronic charge transfer through the hydrogen bond. Chem Phys Lett 2017. [DOI: 10.1016/j.cplett.2017.04.034] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/24/2022]
|
18
|
Karandashev K, Vaníček J. Accelerating equilibrium isotope effect calculations. I. Stochastic thermodynamic integration with respect to mass. J Chem Phys 2017. [DOI: 10.1063/1.4981260] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023] Open
Affiliation(s)
- Konstantin Karandashev
- Laboratory of Theoretical Physical Chemistry, Institut des Sciences et Ingénierie Chimiques, École Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| | - Jiří Vaníček
- Laboratory of Theoretical Physical Chemistry, Institut des Sciences et Ingénierie Chimiques, École Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| |
Collapse
|
19
|
Moustafa SG, Schultz AJ, Kofke DA. Harmonically Assisted Methods for Computing the Free Energy of Classical Crystals by Molecular Simulation: A Comparative Study. J Chem Theory Comput 2017; 13:825-834. [DOI: 10.1021/acs.jctc.6b01082] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Sabry G. Moustafa
- Department of Chemical and
Biological Engineering, University at Buffalo, The State University of New York, Buffalo, New York 14260-4200, United States
| | - Andrew J. Schultz
- Department of Chemical and
Biological Engineering, University at Buffalo, The State University of New York, Buffalo, New York 14260-4200, United States
| | - David A. Kofke
- Department of Chemical and
Biological Engineering, University at Buffalo, The State University of New York, Buffalo, New York 14260-4200, United States
| |
Collapse
|
20
|
Affiliation(s)
- Venkat Kapil
- Laboratory of Computational Science and Modelling, Institute of Materials, Ecole Polytechnique Fédérale de Lausanne, Lausanne, Switzerland
| | - Jörg Behler
- Lehrstuhl für Theoretische Chemie, Ruhr-Universität Bochum, Bochum, Germany
| | - Michele Ceriotti
- Laboratory of Computational Science and Modelling, Institute of Materials, Ecole Polytechnique Fédérale de Lausanne, Lausanne, Switzerland
| |
Collapse
|
21
|
Litman YE, Videla PE, Rodriguez J, Laria D. Positional Isotope Exchange in HX·(H2O)n (X = F, I) Clusters at Low Temperatures. J Phys Chem A 2016; 120:7213-24. [DOI: 10.1021/acs.jpca.6b06681] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Yair E. Litman
- Departamento
de Química Inorgánica, Analítica y Química-Física
and INQUIMAE-CONICET, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires Ciudad Universitaria, Pabellón II, 1428 Buenos Aires, Argentina
| | - Pablo E. Videla
- Departamento
de Química Inorgánica, Analítica y Química-Física
and INQUIMAE-CONICET, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires Ciudad Universitaria, Pabellón II, 1428 Buenos Aires, Argentina
| | - Javier Rodriguez
- Departamento
de Física de la Materia Condensada, Comisión Nacional de Energía Atómica, Avenida Libertador 8250, 1429 Buenos Aires, Argentina
- ECyT, UNSAM, Martín
de Irigoyen 3100, 1650 San Martín, Pcia. de Buenos Aires, Argentina
| | - Daniel Laria
- Departamento
de Química Inorgánica, Analítica y Química-Física
and INQUIMAE-CONICET, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires Ciudad Universitaria, Pabellón II, 1428 Buenos Aires, Argentina
- Departamento
de Física de la Materia Condensada, Comisión Nacional de Energía Atómica, Avenida Libertador 8250, 1429 Buenos Aires, Argentina
| |
Collapse
|
22
|
Uhl F, Marx D, Ceriotti M. Accelerated path integral methods for atomistic simulations at ultra-low temperatures. J Chem Phys 2016; 145:054101. [DOI: 10.1063/1.4959602] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Felix Uhl
- 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
| | - Michele Ceriotti
- Laboratory of Computational Science and Modelling, Institute of Materials, École Polytechnique Fédérale de Lausanne, 1015 Lausanne, Switzerland
| |
Collapse
|
23
|
Abstract
Path-integral free energy perturbation (PI-FEP) theory is presented to directly determine the ratio of quantum mechanical partition functions of different isotopologs in a single simulation. Furthermore, a double averaging strategy is used to carry out the practical simulation, separating the quantum mechanical path integral exactly into two separate calculations, one corresponding to a classical molecular dynamics simulation of the centroid coordinates, and another involving free-particle path-integral sampling over the classical, centroid positions. An integrated centroid path-integral free energy perturbation and umbrella sampling (PI-FEP/UM, or simply, PI-FEP) method along with bisection sampling was summarized, which provides an accurate and fast convergent method for computing kinetic isotope effects for chemical reactions in solution and in enzymes. The PI-FEP method is illustrated by a number of applications, to highlight the computational precision and accuracy, the rule of geometrical mean in kinetic isotope effects, enhanced nuclear quantum effects in enzyme catalysis, and protein dynamics on temperature dependence of kinetic isotope effects.
Collapse
Affiliation(s)
- J Gao
- Theoretical Chemistry Institute, Jilin University, Changchun, Jilin Province, PR China; University of Minnesota, Minneapolis, MN, United States.
| |
Collapse
|
24
|
Fang W, Chen J, Rossi M, Feng Y, Li XZ, Michaelides A. Inverse Temperature Dependence of Nuclear Quantum Effects in DNA Base Pairs. J Phys Chem Lett 2016; 7:2125-31. [PMID: 27195654 PMCID: PMC4933496 DOI: 10.1021/acs.jpclett.6b00777] [Citation(s) in RCA: 27] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/08/2023]
Abstract
Despite the inherently quantum mechanical nature of hydrogen bonding, it is unclear how nuclear quantum effects (NQEs) alter the strengths of hydrogen bonds. With this in mind, we use ab initio path integral molecular dynamics to determine the absolute contribution of NQEs to the binding in DNA base pair complexes, arguably the most important hydrogen-bonded systems of all. We find that depending on the temperature, NQEs can either strengthen or weaken the binding within the hydrogen-bonded complexes. As a somewhat counterintuitive consequence, NQEs can have a smaller impact on hydrogen bond strengths at cryogenic temperatures than at room temperature. We rationalize this in terms of a competition of NQEs between low-frequency and high-frequency vibrational modes. Extending this idea, we also propose a simple model to predict the temperature dependence of NQEs on hydrogen bond strengths in general.
Collapse
Affiliation(s)
- Wei Fang
- Thomas Young Centre, London Centre for Nanotechnology, Department of Chemistry, and Department of
Physics and Astronomy, University College
London, London WC1E 6BT, United Kingdom
| | - Ji Chen
- Thomas Young Centre, London Centre for Nanotechnology, Department of Chemistry, and Department of
Physics and Astronomy, University College
London, London WC1E 6BT, United Kingdom
| | - Mariana Rossi
- Physical
and Theoretical Chemistry Lab, University
of Oxford, South Parks
Road, OX1 3QZ Oxford, United Kingdom
| | - Yexin Feng
- School
of Physics and Electronics, Hunan University, Changsha 410082, People’s Republic of China
| | - Xin-Zheng Li
- International
Center for Quantum Materials, School of Physics and Collaborative
Innovation Center of Quantum Matter, Peking
University, Beijing 100871, People’s Republic of China
- E-mail: (X.-Z.L.)
| | - Angelos Michaelides
- Thomas Young Centre, London Centre for Nanotechnology, Department of Chemistry, and Department of
Physics and Astronomy, University College
London, London WC1E 6BT, United Kingdom
- E-mail: (A.M.)
| |
Collapse
|
25
|
Ceriotti M, Fang W, Kusalik PG, McKenzie RH, Michaelides A, Morales MA, Markland TE. Nuclear Quantum Effects in Water and Aqueous Systems: Experiment, Theory, and Current Challenges. Chem Rev 2016; 116:7529-50. [DOI: 10.1021/acs.chemrev.5b00674] [Citation(s) in RCA: 339] [Impact Index Per Article: 42.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Affiliation(s)
- Michele Ceriotti
- Laboratory
of Computational Science and Modeling, Institute of Materials, École Polytechnique Fédérale de Lausanne, 1015 Lausanne, Switzerland
| | - Wei Fang
- Thomas
Young Centre, London Centre for Nanotechnology and Department of Physics
and Astronomy, University College London, London WC1E 6BT, United Kingdom
| | - Peter G. Kusalik
- Department
of Chemistry, University of Calgary, 2500 University Drive NW, Calgary, Alberta T2N 1N4, Canada
| | - Ross H. McKenzie
- School
of Mathematics and Physics, University of Queensland, Brisbane, 4072 Queensland Australia
| | - Angelos Michaelides
- Thomas
Young Centre, London Centre for Nanotechnology and Department of Physics
and Astronomy, University College London, London WC1E 6BT, United Kingdom
| | - Miguel A. Morales
- Lawrence Livermore National Laboratory, Livermore, California 94550, United States
| | - Thomas E. Markland
- Department
of Chemistry, Stanford University, 333 Campus Drive, Stanford, California 94305, United States
| |
Collapse
|
26
|
|
27
|
Karandashev K, Vaníček J. Accelerating quantum instanton calculations of the kinetic isotope effects. J Chem Phys 2015; 143:194104. [DOI: 10.1063/1.4935701] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Affiliation(s)
- Konstantin Karandashev
- Laboratory of Theoretical Physical Chemistry, Institut des Sciences et Ingénierie Chimiques, Ecole Polytechnique Fédérale de Lausanne (EPFL), Lausanne CH-1015, Switzerland
| | - Jiří Vaníček
- Laboratory of Theoretical Physical Chemistry, Institut des Sciences et Ingénierie Chimiques, Ecole Polytechnique Fédérale de Lausanne (EPFL), Lausanne CH-1015, Switzerland
| |
Collapse
|
28
|
Rossi M, Fang W, Michaelides A. Stability of Complex Biomolecular Structures: van der Waals, Hydrogen Bond Cooperativity, and Nuclear Quantum Effects. J Phys Chem Lett 2015; 6:4233-8. [PMID: 26722963 DOI: 10.1021/acs.jpclett.5b01899] [Citation(s) in RCA: 29] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/16/2023]
Abstract
Biomolecules are complex systems stabilized by a delicate balance of weak interactions, making it important to assess all energetic contributions in an accurate manner. However, it is a priori unclear which contributions make more of an impact. Here, we examine stacked polyglutamine (polyQ) strands, a peptide repeat often found in amyloid aggregates. We investigate the role of hydrogen bond (HB) cooperativity, van der Waals (vdW) dispersion interactions, and quantum contributions to free energies, including anharmonicities through density functional theory and ab initio path integral simulations. Of these various factors, we find that the largest impact on structural stabilization comes from vdW interactions. HB cooperativity is the second largest contribution as the size of the stacked chain grows. Competing nuclear quantum effects make the net quantum contribution small but very sensitive to anharmonicities, vdW, and the number of HBs. Our results suggest that a reliable treatment of these systems can only be attained by considering all of these components.
Collapse
Affiliation(s)
- Mariana Rossi
- Physical and Theoretical Chemistry Lab, University of Oxford , South Parks Road, OX1 3QZ Oxford, United Kingdom
- St. Edmund Hall , Queen's Lane, OX1 4AR Oxford, United Kingdom
| | - Wei Fang
- Thomas Young Centre, London Centre for Nanotechnology, and Department of Chemistry, University College London , 17-19 Gordon Street, WC1H 0AH London, United Kingdom
| | - Angelos Michaelides
- Thomas Young Centre, London Centre for Nanotechnology, and Department of Chemistry, University College London , 17-19 Gordon Street, WC1H 0AH London, United Kingdom
| |
Collapse
|
29
|
Poltavsky I, Tkatchenko A. Modeling quantum nuclei with perturbed path integral molecular dynamics. Chem Sci 2015; 7:1368-1372. [PMID: 29910893 PMCID: PMC5975916 DOI: 10.1039/c5sc03443d] [Citation(s) in RCA: 29] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/12/2015] [Accepted: 10/28/2015] [Indexed: 11/25/2022] Open
Abstract
Here we combine perturbation theory with the Feynman–Kac imaginary-time path integral approach to quantum mechanics for modeling quantum nuclear effects.
The quantum nature of nuclear motions plays a vital role in the structure, stability, and thermodynamics of molecules and materials. The standard approach to model nuclear quantum fluctuations in chemical and biological systems is to use path-integral molecular dynamics. Unfortunately, conventional path-integral simulations can have an exceedingly large computational cost due to the need to employ an excessive number of coupled classical subsystems (beads) for quantitative accuracy. Here, we combine perturbation theory with the Feynman–Kac imaginary-time path integral approach to quantum mechanics and derive an improved non-empirical partition function and estimators to calculate converged quantum observables. Our perturbed path-integral (PPI) method requires the same ingredients as the conventional approach, but increases the accuracy and efficiency of path integral simulations by an order of magnitude. Results are presented for the thermodynamics of fundamental model systems, an empirical water model containing 256 water molecules within periodic boundary conditions, and ab initio simulations of nitrogen and benzene molecules. For all of these examples, PPI simulations with 4 to 8 classical beads recover the nuclear quantum contribution to the total energy and heat capacity at room temperature within a 3% accuracy, paving the way toward seamless modeling of nuclear quantum effects in realistic molecules and materials.
Collapse
Affiliation(s)
- Igor Poltavsky
- Fritz-Haber-Institut der Max-Planck-Gesellschaft , Faradayweg 4-6 , 14195 Berlin , Germany .
| | - Alexandre Tkatchenko
- Fritz-Haber-Institut der Max-Planck-Gesellschaft , Faradayweg 4-6 , 14195 Berlin , Germany .
| |
Collapse
|
30
|
Sinitskiy AV, Voth GA. A reductionist perspective on quantum statistical mechanics: Coarse-graining of path integrals. J Chem Phys 2015; 143:094104. [DOI: 10.1063/1.4929790] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Anton V. Sinitskiy
- Department of Chemistry, James Franck Institute, Institute for Biophysical Dynamics, and Computation Institute, The University of Chicago, 5735 S. Ellis Ave., Chicago, Illinois 60637, USA
| | - Gregory A. Voth
- Department of Chemistry, James Franck Institute, Institute for Biophysical Dynamics, and Computation Institute, The University of Chicago, 5735 S. Ellis Ave., Chicago, Illinois 60637, USA
| |
Collapse
|
31
|
Cheng B, Ceriotti M. Direct path integral estimators for isotope fractionation ratios. J Chem Phys 2014; 141:244112. [DOI: 10.1063/1.4904293] [Citation(s) in RCA: 25] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
32
|
Quantum delocalization of protons in the hydrogen-bond network of an enzyme active site. Proc Natl Acad Sci U S A 2014; 111:18454-9. [PMID: 25503367 DOI: 10.1073/pnas.1417923111] [Citation(s) in RCA: 89] [Impact Index Per Article: 8.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
Enzymes use protein architectures to create highly specialized structural motifs that can greatly enhance the rates of complex chemical transformations. Here, we use experiments, combined with ab initio simulations that exactly include nuclear quantum effects, to show that a triad of strongly hydrogen-bonded tyrosine residues within the active site of the enzyme ketosteroid isomerase (KSI) facilitates quantum proton delocalization. This delocalization dramatically stabilizes the deprotonation of an active-site tyrosine residue, resulting in a very large isotope effect on its acidity. When an intermediate analog is docked, it is incorporated into the hydrogen-bond network, giving rise to extended quantum proton delocalization in the active site. These results shed light on the role of nuclear quantum effects in the hydrogen-bond network that stabilizes the reactive intermediate of KSI, and the behavior of protons in biological systems containing strong hydrogen bonds.
Collapse
|