1
|
Iyengar SS, Schlegel HB, Sumner I, Li J. Rare Events Sampling Methods for Quantum and Classical Ab Initio Molecular Dynamics. J Phys Chem A 2024; 128:5386-5397. [PMID: 38951489 DOI: 10.1021/acs.jpca.3c07385] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/03/2024]
Abstract
We provide an approach to sample rare events during classical ab initio molecular dynamics and quantum wavepacket dynamics. For classical AIMD, a set of fictitious degrees of freedom are introduced that may harmonically interact with the electronic and nuclear degrees of freedom to steer the dynamics in a conservative fashion toward energetically forbidden regions. A similar approach when introduced for quantum wavepacket dynamics has the effect of biasing the trajectory of the wavepacket centroid toward the regions of the potential surface that are difficult to sample. The approach is demonstrated for a phenol-amine system, which is a prototypical problem for condensed phase-proton transfer, and for model potentials undergoing wavepacket dynamics. In all cases, the approach yields trajectories that conserve energy while sampling rare events.
Collapse
Affiliation(s)
- Srinivasan S Iyengar
- Department of Chemistry, Department of Physics, and the Indiana University Quantum Science and Engineering Center (IU-QSEC), Indiana University, 800 E. Kirkwood Avenue, Bloomington 47405, Indiana, United States
| | - H Bernhard Schlegel
- Department of Chemistry, Wayne State University, Detroit 48202, Michigan, United States
| | - Isaiah Sumner
- Department of Chemistry and Biochemistry, James Madison University, Harrisonburg 22807, Virginia, United States
| | - Junjie Li
- Texas Advanced Computing Center, The University of Texas at Austin, Austin 78758, Texas, United States
| |
Collapse
|
2
|
Iyengar SS, Zhang JH, Saha D, Ricard TC. Graph-| Q⟩⟨ C|: A Quantum Algorithm with Reduced Quantum Circuit Depth for Electronic Structure. J Phys Chem A 2023; 127:9334-9345. [PMID: 37906738 DOI: 10.1021/acs.jpca.3c04261] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/02/2023]
Abstract
The accurate determination of chemical properties is known to have a critical impact on multiple fundamental chemical problems but is deeply hindered by the steep algebraic scaling of electron correlation calculations and the exponential scaling of quantum nuclear dynamics. With the advent of new quantum computing hardware and associated developments in creating new paradigms for quantum software, this avenue has been recognized as perhaps one way to address exponentially complex challenges in quantum chemistry and molecular dynamics. In this paper, we discuss a new approach to drastically reduce the quantum circuit depth (by several orders of magnitude) and help improve the accuracy in the quantum computation of electron correlation energies for large molecular systems. The method is derived from a graph-theoretic approach to molecular fragmentation and enables us to create a family of projection operators that decompose quantum circuits into separate unitary processes. Some of these processes can be treated on quantum hardware and others on classical hardware in a completely asynchronous and parallel fashion. Numerical benchmarks are provided through the computation of unitary coupled-cluster singles and doubles (UCCSD) energies for medium-sized protonated and neutral water clusters using the new quantum algorithms presented here.
Collapse
Affiliation(s)
- Srinivasan S Iyengar
- Department of Chemistry, Department of Physics, and the Indiana University Quantum Science and Engineering Center (IU-QSEC), Indiana University, 800 E. Kirkwood Avenue, Bloomington, Indiana 47405, United States
| | - Juncheng Harry Zhang
- Department of Chemistry, Department of Physics, and the Indiana University Quantum Science and Engineering Center (IU-QSEC), Indiana University, 800 E. Kirkwood Avenue, Bloomington, Indiana 47405, United States
| | - Debadrita Saha
- Department of Chemistry, Department of Physics, and the Indiana University Quantum Science and Engineering Center (IU-QSEC), Indiana University, 800 E. Kirkwood Avenue, Bloomington, Indiana 47405, United States
| | - Timothy C Ricard
- Department of Chemistry, Department of Physics, and the Indiana University Quantum Science and Engineering Center (IU-QSEC), Indiana University, 800 E. Kirkwood Avenue, Bloomington, Indiana 47405, United States
| |
Collapse
|
3
|
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
|
4
|
Kumar A, DeGregorio N, Ricard T, Iyengar SS. Graph-Theoretic Molecular Fragmentation for Potential Surfaces Leads Naturally to a Tensor Network Form and Allows Accurate and Efficient Quantum Nuclear Dynamics. J Chem Theory Comput 2022; 18:7243-7259. [PMID: 36332133 DOI: 10.1021/acs.jctc.2c00484] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022]
Abstract
Molecular fragmentation methods have revolutionized quantum chemistry. Here, we use a graph-theoretically generated molecular fragmentation method, to obtain accurate and efficient representations for multidimensional potential energy surfaces and the quantum time-evolution operator, which plays a critical role in quantum chemical dynamics. In doing so, we find that the graph-theoretic fragmentation approach naturally reduces the potential portion of the time-evolution operator into a tensor network that contains a stream of coupled lower-dimensional propagation steps to potentially achieve quantum dynamics with reduced complexity. Furthermore, the fragmentation approach used here has previously been shown to allow accurate and efficient computation of post-Hartree-Fock electronic potential energy surfaces, which in many cases has been shown to be at density functional theory cost. Thus, by combining the advantages of molecular fragmentation with the tensor network formalism, the approach yields an on-the-fly quantum dynamics scheme where both the electronic potential calculation and nuclear propagation portion are enormously simplified through a single stroke. The method is demonstrated by computing approximations to the propagator and to potential surfaces for a set of coupled nuclear dimensions within a protonated water wire problem exhibiting the Grotthuss mechanism of proton transport. In all cases, our approach has been shown to reduce the complexity of representing the quantum propagator, and by extension action of the propagator on an initial wavepacket, by several orders, with minimal loss in accuracy.
Collapse
Affiliation(s)
- Anup Kumar
- Department of Chemistry, and the Indiana University Quantum Science and Engineering Center (IU-QSEC), Indiana University, Bloomington, Indiana 47405, United States
| | - Nicole DeGregorio
- Department of Chemistry, and the Indiana University Quantum Science and Engineering Center (IU-QSEC), Indiana University, Bloomington, Indiana 47405, United States
| | - Timothy Ricard
- Department of Chemistry, and the Indiana University Quantum Science and Engineering Center (IU-QSEC), Indiana University, Bloomington, Indiana 47405, United States
| | - Srinivasan S Iyengar
- Department of Chemistry, and the Indiana University Quantum Science and Engineering Center (IU-QSEC), Indiana University, Bloomington, Indiana 47405, United States
| |
Collapse
|
5
|
Zhang JH, Iyengar SS. Graph-| Q⟩⟨ C|, a Graph-Based Quantum/Classical Algorithm for Efficient Electronic Structure on Hybrid Quantum/Classical Hardware Systems: Improved Quantum Circuit Depth Performance. J Chem Theory Comput 2022; 18:2885-2899. [PMID: 35412836 DOI: 10.1021/acs.jctc.1c01303] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
We present a procedure to reduce the depth of quantum circuits and improve the accuracy of results in computing post-Hartree-Fock electronic structure energies in large molecular systems. The method is based on molecular fragmentation where a molecular system is divided into overlapping fragments through a graph-theoretic procedure. This allows us to create a set of projection operators that decompose the unitary evolution of the full system into separate sets of processes, some of which can be treated on quantum hardware and others on classical hardware. Thus, we develop a procedure for an electronic structure that can be asynchronously spawned onto a potentially large ensemble of classical and quantum hardware systems. We demonstrate this method by computing Unitary Coupled Cluster Singles and Doubles (UCCSD) energies for a set of [H2]n clusters, with n ranging from 4 to 128. We implement our methodology using quantum circuits, and when these quantum circuits are processed on a quantum simulator, we obtain energies in agreement with the UCCSD energies in the milli-hartree energy range. We also show that our circuit decomposition approach yields up to 9 orders of magnitude reduction in the number of CNOT gates and quantum circuit depth for the large-sized clusters when compared to a standard quantum circuit implementation available on IBM's Quantum Information Science kit, known as Qiskit.
Collapse
Affiliation(s)
- Juncheng Harry Zhang
- Department of Chemistry and Department of Physics, Indiana University, 800 E. Kirkwood Avenue, Bloomington, Indiana 47405, United States
| | - Srinivasan S Iyengar
- Department of Chemistry and Department of Physics, Indiana University, 800 E. Kirkwood Avenue, Bloomington, Indiana 47405, United States
| |
Collapse
|
6
|
Brémond É, Savarese M, Rega N, Ciofini I, Adamo C. Free Energy Profiles of Proton Transfer Reactions: Density Functional Benchmark from Biased Ab Initio Dynamics. J Chem Theory Comput 2022; 18:1501-1511. [PMID: 35129987 DOI: 10.1021/acs.jctc.1c01002] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
By coupling an enhanced sampling algorithm with an orbital-localized variant of Car-Parrinello molecular dynamics, the so-called atomic centered density matrix propagation model, we reconstruct the free energy profiles along reaction pathways using different density functional approximations (DFAs) ranging from locals to hybrids. In particular, we compare the computed free energy barrier height of proton transfer (PT) reactions to those obtained by a more traditional static approach, based on the intrinsic reaction coordinate (IRC), for two case systems, namely malonaldehyde and formic acid dimer. The obtained results show that both the IRC profiles and the potentials of mean force, derived from biased dynamic trajectories, are very sensitive to the density functional approximation applied. More precisely, we observe that, with the notable exception of M06-L, local density functionals always strongly underestimate the reaction barrier heights. More generally, we find that also the shape of the free energy profile is very sensitive to the density functional choice, thus highlighting the effect, often neglected, that the choice of DFA has also in the case of dynamics simulations.
Collapse
Affiliation(s)
- Éric Brémond
- Université de Paris, ITODYS, CNRS, F-75006 Paris, France
| | - Marika Savarese
- Chimie ParisTech-PSL, CNRS, Institute of Chemistry for Health and Life Sciences, F-75005 Paris, France
| | - Nadia Rega
- Dipartimento di Scienze Chimiche, Università di Napoli Federico II, Complesso Universitario di M.S. Angelo, via Cintia, I-80126 Napoli, Italy.,Scuola Superiore Meridionale, Largo S. Marcellino 10, I-80138 Napoli, Italy.,Centro Interdipartimentale di Ricerca sui Biomateriali (CRIB), Piazzale Tecchio 80, I-80125, Napoli, Italy
| | - Ilaria Ciofini
- Chimie ParisTech-PSL, CNRS, Institute of Chemistry for Health and Life Sciences, F-75005 Paris, France
| | - Carlo Adamo
- Chimie ParisTech-PSL, CNRS, Institute of Chemistry for Health and Life Sciences, F-75005 Paris, France.,Institut Universitaire de France, 103 Boulevard Saint Michel, F-75005 Paris, France
| |
Collapse
|
7
|
Saha D, Iyengar SS, Richerme P, Smith JM, Sabry A. Mapping Quantum Chemical Dynamics Problems to Spin-Lattice Simulators. J Chem Theory Comput 2021; 17:6713-6732. [PMID: 34694820 DOI: 10.1021/acs.jctc.1c00688] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The accurate computational determination of chemical, materials, biological, and atmospheric properties has a critical impact on a wide range of health and environmental problems, but is deeply limited by the computational scaling of quantum mechanical methods. The complexity of quantum chemical studies arises from the steep algebraic scaling of electron correlation methods and the exponential scaling in studying nuclear dynamics and molecular flexibility. To date, efforts to apply quantum hardware to such quantum chemistry problems have focused primarily on electron correlation. Here, we provide a framework that allows for the solution of quantum chemical nuclear dynamics by mapping these to quantum spin-lattice simulators. Using the example case of a short-strong hydrogen-bonded system, we construct the Hamiltonian for the nuclear degrees of freedom on a single Born-Oppenheimer surface and show how it can be transformed to a generalized Ising model Hamiltonian. We then demonstrate a method to determine the local fields and spin-spin couplings needed to identically match the molecular and spin-lattice Hamiltonians. We describe a protocol to determine the on-site and intersite coupling parameters of this Ising Hamiltonian from the Born-Oppenheimer potential and nuclear kinetic energy operator. Our approach represents a paradigm shift in the methods used to study quantum nuclear dynamics, opening the possibility to solve both electronic structure and nuclear dynamics problems using quantum computing systems.
Collapse
Affiliation(s)
- Debadrita Saha
- Department of Chemistry, and the Indiana University Quantum Science and Engineering Center (IU-QSEC), Indiana University, Bloomington, Indiana 47405, United States
| | - Srinivasan S Iyengar
- Department of Chemistry, and the Indiana University Quantum Science and Engineering Center (IU-QSEC), Indiana University, Bloomington, Indiana 47405, United States
| | - Philip Richerme
- Department of Physics and the Indiana University Quantum Science and Engineering Center (IU-QSEC), Indiana University, Bloomington, Indiana 47405, United States
| | - Jeremy M Smith
- Department of Chemistry, Indiana University, 800 E. Kirkwood Avenue, Bloomington, Indiana 47405, United States
| | - Amr Sabry
- Department of Computer Science, School of Informatics, Computing, and Engineering, and the Indiana University Quantum Science and Engineering Center (IU-QSEC), Indiana University, Bloomington, Indiana 47405, United States
| |
Collapse
|
8
|
Trumbore CN. Shear-Induced Amyloid Aggregation in the Brain: V. Are Alzheimer's and Other Amyloid Diseases Initiated in the Lower Brain and Brainstem by Cerebrospinal Fluid Flow Stresses? J Alzheimers Dis 2021; 79:979-1002. [PMID: 33386802 PMCID: PMC7990457 DOI: 10.3233/jad-201025] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 11/16/2020] [Indexed: 12/13/2022]
Abstract
Amyloid-β (Aβ) and tau oligomers have been identified as neurotoxic agents responsible for causing Alzheimer's disease (AD). Clinical trials using Aβ and tau as targets have failed, giving rise to calls for new research approaches to combat AD. This paper provides such an approach. Most basic AD research has involved quiescent Aβ and tau solutions. However, studies involving laminar and extensional flow of proteins have demonstrated that mechanical agitation of proteins induces or accelerates protein aggregation. Recent MRI brain studies have revealed high energy, chaotic motion of cerebrospinal fluid (CSF) in lower brain and brainstem regions. These and studies showing CSF flow within the brain have shown that there are two energetic hot spots. These are within the third and fourth brain ventricles and in the neighborhood of the circle of Willis blood vessel region. These two regions are also the same locations as those of the earliest Aβ and tau AD pathology. In this paper, it is proposed that cardiac systolic pulse waves that emanate from the major brain arteries in the lower brain and brainstem regions and whose pulse waves drive CSF flows within the brain are responsible for initiating AD and possibly other amyloid diseases. It is further proposed that the triggering of these diseases comes about because of the strengthening of systolic pulses due to major artery hardening that generates intense CSF extensional flow stress. Such stress provides the activation energy needed to induce conformational changes of both Aβ and tau within the lower brain and brainstem region, producing unique neurotoxic oligomer molecule conformations that induce AD.
Collapse
Affiliation(s)
- Conrad N. Trumbore
- Department of Chemistry and Biochemistry, University of Delaware, Newark, DE, USA
| |
Collapse
|
9
|
Jaglan R, Mandal D. The role of potential energy surface in quantum mechanical tunneling: A computational perspective. COMPUT THEOR CHEM 2020. [DOI: 10.1016/j.comptc.2020.112920] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/16/2022]
|
10
|
Banik S, Sansi AK, Nandan S, Roy TK. On the Proton Shuttle Motion in Protonated Acetylene: An Electronic Structure Perspective. ChemistrySelect 2020. [DOI: 10.1002/slct.202002524] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022]
Affiliation(s)
- Subrata Banik
- Department of ChemistrySchool of Chemical and BiotechnologySASTRA Deemed University Thanjavur 613401 Tamilnadu India
| | - Ankit Kumar Sansi
- Department of Chemistry and Chemical SciencesCentral University of Jammu, Rahya-Suchani (Bagla) Jammu 180011 India
| | - Shiv Nandan
- Department of Chemistry and Chemical SciencesCentral University of Jammu, Rahya-Suchani (Bagla) Jammu 180011 India
| | - Tapta Kanchan Roy
- Department of Chemistry and Chemical SciencesCentral University of Jammu, Rahya-Suchani (Bagla) Jammu 180011 India
| |
Collapse
|
11
|
Prah A, Ogrin P, Mavri J, Stare J. Nuclear quantum effects in enzymatic reactions: simulation of the kinetic isotope effect of phenylethylamine oxidation catalyzed by monoamine oxidase A. Phys Chem Chem Phys 2020; 22:6838-6847. [DOI: 10.1039/d0cp00131g] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
By using computational techniques for quantizing nuclear motion one can accurately reproduce kinetic isotope effect of enzymatic reactions, as demonstrated for phenylethylamine oxidation catalyzed by the monoamine oxidase A enzyme.
Collapse
Affiliation(s)
- Alja Prah
- Theory Department
- National Institute of Chemistry
- Ljubljana
- Slovenia
- University of Ljubljana
| | - Peter Ogrin
- Theory Department
- National Institute of Chemistry
- Ljubljana
- Slovenia
- University of Ljubljana
| | - Janez Mavri
- Theory Department
- National Institute of Chemistry
- Ljubljana
- Slovenia
| | - Jernej Stare
- Theory Department
- National Institute of Chemistry
- Ljubljana
- Slovenia
| |
Collapse
|
12
|
DeGregorio N, Iyengar SS. Challenges in constructing accurate methods for hydrogen transfer reactions in large biological assemblies: rare events sampling for mechanistic discovery and tensor networks for quantum nuclear effects. Faraday Discuss 2020; 221:379-405. [DOI: 10.1039/c9fd00071b] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
We present two methods that address the computational complexities arising in hydrogen transfer reactions in enzyme active sites.
Collapse
Affiliation(s)
- Nicole DeGregorio
- Department of Chemistry
- Department of Physics
- Indiana University
- Bloomington
- USA
| | | |
Collapse
|
13
|
Trumbore CN. Shear-Induced Amyloid Formation in the Brain: III. The Roles of Shear Energy and Seeding in a Proposed Shear Model. J Alzheimers Dis 2019; 65:47-70. [PMID: 30040710 PMCID: PMC6087447 DOI: 10.3233/jad-171003] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Abstract
If cerebrospinal and interstitial fluids move through very narrow brain flow channels, these restrictive surroundings generate varying levels of fluid shear and different shear rates, and dissolved amyloid monomers absorb different shear energies. It is proposed that dissolved amyloid-β protein (Aβ) and other amyloid monomers undergo shear-induced conformational changes that ultimately lead to amyloid monomer aggregation even at very low brain flow and shear rates. Soluble Aβ oligomers taken from diseased brains initiate in vivo amyloid formation in non-diseased brains. The brain environment is apparently responsible for this result. A mechanism involving extensional shear is proposed for the formation of a seed Aβ monomer molecule that ultimately promotes templated conformational change of other Aβ molecules. Under non-quiescent, non-equilibrium conditions, gentle extensional shear within the brain parenchyma, and perhaps even during laboratory preparation of Aβ samples, may be sufficient to cause subtle conformational changes in these monomers. These result from brain processes that significantly lower the high activation energy predicted for the quiescent Aβ dimerization process. It is further suggested that changes in brain location and changes brought about by aging expose Aβ molecules to different shear rates, total shear, and types of shear, resulting in different conformational changes in these molecules. The consequences of such changes caused by variable shear energy are proposed to underlie formation of amyloid strains causing different amyloid diseases. Amyloid researchers are urged to undertake studies with amyloids, anti-amyloid drugs, and antibodies while all of these are under shear conditions similar to those in the brain.
Collapse
Affiliation(s)
- Conrad N Trumbore
- Department of Chemistry and Biochemistry, University of Delaware, Newark, DE, USA
| |
Collapse
|
14
|
DeGregorio N, Iyengar SS. Adaptive Dimensional Decoupling for Compression of Quantum Nuclear Wave Functions and Efficient Potential Energy Surface Representations through Tensor Network Decomposition. J Chem Theory Comput 2019; 15:2780-2796. [DOI: 10.1021/acs.jctc.8b01113] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Nicole DeGregorio
- Department of Chemistry and Department of Physics, Indiana University, 800 E. Kirkwood Ave, Bloomington, Indiana 47405, United States
| | - Srinivasan S. Iyengar
- Department of Chemistry and Department of Physics, Indiana University, 800 E. Kirkwood Ave, Bloomington, Indiana 47405, United States
| |
Collapse
|
15
|
Ricard TC, Haycraft C, Iyengar SS. Adaptive, Geometric Networks for Efficient Coarse-Grained Ab Initio Molecular Dynamics with Post-Hartree–Fock Accuracy. J Chem Theory Comput 2018; 14:2852-2866. [DOI: 10.1021/acs.jctc.8b00186] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Affiliation(s)
- Timothy C. Ricard
- Department of Chemistry and Department of Physics, Indiana University, 800 East Kirkwood Avenue, Bloomington, Indiana 47405, United States
| | - Cody Haycraft
- Department of Chemistry and Department of Physics, Indiana University, 800 East Kirkwood Avenue, Bloomington, Indiana 47405, United States
| | - Srinivasan S. Iyengar
- Department of Chemistry and Department of Physics, Indiana University, 800 East Kirkwood Avenue, Bloomington, Indiana 47405, United States
| |
Collapse
|
16
|
DeGregorio N, Iyengar SS. Efficient and Adaptive Methods for Computing Accurate Potential Surfaces for Quantum Nuclear Effects: Applications to Hydrogen-Transfer Reactions. J Chem Theory Comput 2017; 14:30-47. [DOI: 10.1021/acs.jctc.7b00927] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022]
Affiliation(s)
- Nicole DeGregorio
- Department of Chemistry and
Department of Physics, Indiana University, 800 E. Kirkwood Avenue, Bloomington, Indiana 47405, United States
| | - Srinivasan S. Iyengar
- Department of Chemistry and
Department of Physics, Indiana University, 800 E. Kirkwood Avenue, Bloomington, Indiana 47405, United States
| |
Collapse
|
17
|
Jevtic S, Anders J. A qualitative quantum rate model for hydrogen transfer in soybean lipoxygenase. J Chem Phys 2017; 147:114108. [PMID: 28938801 DOI: 10.1063/1.4998941] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022] Open
Abstract
The hydrogen transfer reaction catalysed by soybean lipoxygenase (SLO) has been the focus of intense study following observations of a high kinetic isotope effect (KIE). Today high KIEs are generally thought to indicate departure from classical rate theory and are seen as a strong signature of tunnelling of the transferring particle, hydrogen or one of its isotopes, through the reaction energy barrier. In this paper, we build a qualitative quantum rate model with few free parameters that describes the dynamics of the transferring particle when it is exposed to energetic potentials exerted by the donor and the acceptor. The enzyme's impact on the dynamics is modelled by an additional energetic term, an oscillatory contribution known as "gating." By varying two key parameters, the gating frequency and the mean donor-acceptor separation, the model is able to reproduce well the KIE data for SLO wild-type and a variety of SLO mutants over the experimentally accessible temperature range. While SLO-specific constants have been considered here, it is possible to adapt these for other enzymes.
Collapse
Affiliation(s)
- S Jevtic
- Department of Mathematics, Huxley Building, Imperial College, London SW7 2AZ, United Kingdom
| | - J Anders
- CEMPS, Physics and Astronomy, University of Exeter, Exeter EX4 4QL, United Kingdom
| |
Collapse
|
18
|
Hirshberg B, Gerber RB. Mean-Field Methods for Time-Dependent Quantum Dynamics of Many-Atom Systems. ADVANCES IN QUANTUM CHEMISTRY 2017. [DOI: 10.1016/bs.aiq.2017.01.002] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/20/2022]
|
19
|
Alborzpour JP, Tew DP, Habershon S. Efficient and accurate evaluation of potential energy matrix elements for quantum dynamics using Gaussian process regression. J Chem Phys 2016; 145:174112. [DOI: 10.1063/1.4964902] [Citation(s) in RCA: 44] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Jonathan P. Alborzpour
- Department of Chemistry and Centre for Scientific Computing, University of Warwick, Coventry CV4 7AL, United Kingdom
| | - David P. Tew
- School of Chemistry, University of Bristol, Bristol BS8 1TS, United Kingdom
| | - Scott Habershon
- Department of Chemistry and Centre for Scientific Computing, University of Warwick, Coventry CV4 7AL, United Kingdom
| |
Collapse
|
20
|
Li J, Iyengar SS. Ab Initio Molecular Dynamics Using Recursive, Spatially Separated, Overlapping Model Subsystems Mixed within an ONIOM-Based Fragmentation Energy Extrapolation Technique. J Chem Theory Comput 2015; 11:3978-91. [DOI: 10.1021/acs.jctc.5b00433] [Citation(s) in RCA: 21] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Junjie Li
- Department of Chemistry and
Department of Physics, Indiana University, 800 E. Kirkwood Ave, Bloomington, Indiana 47405, United States
| | - Srinivasan S. Iyengar
- Department of Chemistry and
Department of Physics, Indiana University, 800 E. Kirkwood Ave, Bloomington, Indiana 47405, United States
| |
Collapse
|
21
|
Phatak P, Venderley J, Debrota J, Li J, Iyengar SS. Active Site Dynamical Effects in the Hydrogen Transfer Rate-limiting Step in the Catalysis of Linoleic Acid by Soybean Lipoxygenase-1 (SLO-1): Primary and Secondary Isotope Contributions. J Phys Chem B 2015; 119:9532-46. [PMID: 26079999 DOI: 10.1021/acs.jpcb.5b02385] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Using ab initio molecular dynamics (AIMD) simulations that facilitate the treatment of rare events, we probe the active site participation in the rate-determining hydrogen transfer step in the catalytic oxidation of linoleic acid by soybean lipoxygenase-1 (SLO-1). The role of two different active site components is probed. (a) On the hydrogen atom acceptor side of the active site, the hydrogen bonding propensity between the acceptor side hydroxyl group, which is bound to the iron cofactor, and the backbone carboxyl group of isoleucine (residue number 839) is studied toward its role in promoting the hydrogen transfer event. Primary and secondary (H/D) isotope effects are also probed and a definite correlation with subtle secondary H/D isotope effects is found. With increasing average nuclear kinetic energy, the increase in transfer probability is enhanced due to the presence of the hydrogen bond between the backbone carbonyl of I839 and the acceptor oxygen. Further increase in average nuclear kinetic energy reduces the strength of this secondary hydrogen bond which leads to a deterioration in hydrogen transfer rates and finally embrances an Arrhenius-like behavior. (b) On the hydrogen atom donor side, the coupling between vibrational modes predominantly localized on the donor-side linoleic acid group and the reactive mode is probed. There appears to be a qualitative difference in the coupling between modes that belong to linoleic acid and the hydrogen transfer mode, for hydrogen and deuterium transfer. For example, the donor side secondary hydrogen atom is much more labile (by nearly a factor of 5) during deuterium transfer as compared to the case for hydrogen transfer. This appears to indicate a greater coupling between the modes belonging to the linoleic acid scaffold and the deuterium transfer mode and also provides a new rationalization for the abnormal (nonclassical) secondary isotope effect results obtained by Knapp, Rickert, and Klinman in J. Am. Chem. Soc. , 2002 , 124 , 3865 . To substantiate our findings noted in point a above, we have suggested an I839 → A839 or I839 → V839 mutation. This will modify the bulkiness of hydrogen the bonding residue, allowing greater flexibility in the secondary hydrogen bond formation highlighted above and adversely affecting the reaction rate.
Collapse
Affiliation(s)
- Prasad Phatak
- Department of Chemistry and Department of Physics, Indiana University, 800 East Kirkwood Avenue, Bloomington, Indiana 47405, United States
| | - Jordan Venderley
- Department of Chemistry and Department of Physics, Indiana University, 800 East Kirkwood Avenue, Bloomington, Indiana 47405, United States
| | - John Debrota
- Department of Chemistry and Department of Physics, Indiana University, 800 East Kirkwood Avenue, Bloomington, Indiana 47405, United States
| | - Junjie Li
- Department of Chemistry and Department of Physics, Indiana University, 800 East Kirkwood Avenue, Bloomington, Indiana 47405, United States
| | - Srinivasan S Iyengar
- Department of Chemistry and Department of Physics, Indiana University, 800 East Kirkwood Avenue, Bloomington, Indiana 47405, United States
| |
Collapse
|
22
|
Brela MZ, Wójcik MJ, Boczar M, Witek Ł, Yasuda M, Ozaki Y. Car–Parrinello Molecular Dynamics Simulations of Infrared Spectra of Crystalline Vitamin C with Analysis of Double Minimum Proton Potentials for Medium-Strong Hydrogen Bonds. J Phys Chem B 2015; 119:7922-30. [DOI: 10.1021/acs.jpcb.5b02777] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Affiliation(s)
- Mateusz Z. Brela
- Research
Group of Molecular Modelling of Catalytic Processes, Faculty of Chemistry, Jagiellonian University, Ingardena 3, 30-060 Krakow, Poland
| | - Marek J. Wójcik
- Laboratory
of Molecular Spectroscopy, Faculty of Chemistry, Jagiellonian University, Ingardena 3, 30-060 Krakow, Poland
| | - Marek Boczar
- Laboratory
of Molecular Spectroscopy, Faculty of Chemistry, Jagiellonian University, Ingardena 3, 30-060 Krakow, Poland
| | - Łukasz Witek
- Photochemistry
and Luminescence Research Group, Faculty of Chemistry, Jagiellonian University, Ingardena 3, 30-060 Krakow, Poland
| | - Mitsuru Yasuda
- Department
of Chemistry, School of Science and Technology, Kwansei Gakuin University, Sanda, Hyogo 669-1337, Japan
| | - Yukihiro Ozaki
- Department
of Chemistry, School of Science and Technology, Kwansei Gakuin University, Sanda, Hyogo 669-1337, Japan
| |
Collapse
|
23
|
Engel H, Eitan R, Azuri A, Major DT. Nuclear quantum effects in chemical reactions via higher-order path-integral calculations. Chem Phys 2015. [DOI: 10.1016/j.chemphys.2015.01.001] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/25/2022]
|
24
|
Saller MAC, Habershon S. Basis Set Generation for Quantum Dynamics Simulations Using Simple Trajectory-Based Methods. J Chem Theory Comput 2014; 11:8-16. [DOI: 10.1021/ct500657f] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Maximilian A. C. Saller
- Department
of Chemistry and
Centre for Scientific Computing, University of Warwick, Coventry, CV4 7AL, United Kingdom
| | - Scott Habershon
- Department
of Chemistry and
Centre for Scientific Computing, University of Warwick, Coventry, CV4 7AL, United Kingdom
| |
Collapse
|
25
|
Audran G, Brémond P, Marque SR, Siri D, Santelli M. Energetics of the biosynthesis of cyclopentenones from unsaturated fatty acids. Tetrahedron 2014. [DOI: 10.1016/j.tet.2014.09.056] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/24/2022]
|
26
|
Trakhtenberg LI. Tunneling transfer of atomic particles in chemical and biological reactions: The role of intermolecular vibrations and media reorganization. RUSSIAN JOURNAL OF PHYSICAL CHEMISTRY A 2014. [DOI: 10.1134/s003602441411020x] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/23/2022]
|
27
|
Mazzuca JW, Garashchuk S, Jakowski J. The effect of local substrate motion on quantum hydrogen transfer in soybean lipoxygenase-1 modeled with QTES-DFTB dynamics. Chem Phys Lett 2014. [DOI: 10.1016/j.cplett.2014.08.006] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
|
28
|
Prociuk AH, Iyengar SS. A Multiwavelet Treatment of the Quantum Subsystem in Quantum Wavepacket Ab Initio Molecular Dynamics through an Hierarchical Partitioning of Momentum Space. J Chem Theory Comput 2014; 10:2950-63. [DOI: 10.1021/ct5003016] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Alexander H. Prociuk
- Department of Chemistry and
Department of Physics, Indiana University, Bloomington, Indiana 47405, United States
| | - Srinivasan S. Iyengar
- Department of Chemistry and
Department of Physics, Indiana University, Bloomington, Indiana 47405, United States
| |
Collapse
|
29
|
Li J, Li X, Iyengar SS. Vibrational Properties of Hydrogen-Bonded Systems Using the Multireference Generalization to the “On-the-Fly” Electronic Structure within Quantum Wavepacket ab Initio Molecular Dynamics (QWAIMD). J Chem Theory Comput 2014; 10:2265-80. [DOI: 10.1021/ct5002347] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Junjie Li
- Department of Chemistry and
Department of Physics, Indiana University, 800 East Kirkwood Avenue, Bloomington, Indiana 47405, United States
| | - Xiaohu Li
- Department of Chemistry and
Department of Physics, Indiana University, 800 East Kirkwood Avenue, Bloomington, Indiana 47405, United States
| | - Srinivasan S. Iyengar
- Department of Chemistry and
Department of Physics, Indiana University, 800 East Kirkwood Avenue, Bloomington, Indiana 47405, United States
| |
Collapse
|
30
|
Hirao H, Thellamurege N, Zhang X. Applications of density functional theory to iron-containing molecules of bioinorganic interest. Front Chem 2014; 2:14. [PMID: 24809043 PMCID: PMC4010748 DOI: 10.3389/fchem.2014.00014] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/11/2013] [Accepted: 03/10/2014] [Indexed: 12/29/2022] Open
Abstract
The past decades have seen an explosive growth in the application of density functional theory (DFT) methods to molecular systems that are of interest in a variety of scientific fields. Owing to its balanced accuracy and efficiency, DFT plays particularly useful roles in the theoretical investigation of large molecules. Even for biological molecules such as proteins, DFT finds application in the form of, e.g., hybrid quantum mechanics and molecular mechanics (QM/MM), in which DFT may be used as a QM method to describe a higher prioritized region in the system, while a MM force field may be used to describe remaining atoms. Iron-containing molecules are particularly important targets of DFT calculations. From the viewpoint of chemistry, this is mainly because iron is abundant on earth, iron plays powerful (and often enigmatic) roles in enzyme catalysis, and iron thus has the great potential for biomimetic catalysis of chemically difficult transformations. In this paper, we present a brief overview of several recent applications of DFT to iron-containing non-heme synthetic complexes, heme-type cytochrome P450 enzymes, and non-heme iron enzymes, all of which are of particular interest in the field of bioinorganic chemistry. Emphasis will be placed on our own work.
Collapse
Affiliation(s)
- Hajime Hirao
- Division of Chemistry and Biological Chemistry, School of Physical and Mathematical Sciences, Nanyang Technological UniversitySingapore, Singapore
| | | | | |
Collapse
|
31
|
Garashchuk S, Jakowski J, Rassolov VA. Approximate quantum trajectory dynamics for reactive processes in condensed phase. MOLECULAR SIMULATION 2014. [DOI: 10.1080/08927022.2014.907493] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/25/2022]
|
32
|
Migliore A, Polizzi NF, Therien M, Beratan DN. Biochemistry and theory of proton-coupled electron transfer. Chem Rev 2014; 114:3381-465. [PMID: 24684625 PMCID: PMC4317057 DOI: 10.1021/cr4006654] [Citation(s) in RCA: 340] [Impact Index Per Article: 34.0] [Reference Citation Analysis] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/19/2013] [Indexed: 02/01/2023]
Affiliation(s)
- Agostino Migliore
- Department
of Chemistry, Department of Biochemistry, and Department of Physics, Duke University, Durham, North Carolina 27708, United States
| | - Nicholas F. Polizzi
- Department
of Chemistry, Department of Biochemistry, and Department of Physics, Duke University, Durham, North Carolina 27708, United States
| | - Michael
J. Therien
- Department
of Chemistry, Department of Biochemistry, and Department of Physics, Duke University, Durham, North Carolina 27708, United States
| | - David N. Beratan
- Department
of Chemistry, Department of Biochemistry, and Department of Physics, Duke University, Durham, North Carolina 27708, United States
| |
Collapse
|
33
|
Quantum mechanical modeling: a tool for the understanding of enzyme reactions. Biomolecules 2013; 3:662-702. [PMID: 24970187 PMCID: PMC4030948 DOI: 10.3390/biom3030662] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/06/2013] [Revised: 09/17/2013] [Accepted: 09/19/2013] [Indexed: 01/16/2023] Open
Abstract
Most enzyme reactions involve formation and cleavage of covalent bonds, while electrostatic effects, as well as dynamics of the active site and surrounding protein regions, may also be crucial. Accordingly, special computational methods are needed to provide an adequate description, which combine quantum mechanics for the reactive region with molecular mechanics and molecular dynamics describing the environment and dynamic effects, respectively. In this review we intend to give an overview to non-specialists on various enzyme models as well as established computational methods and describe applications to some specific cases. For the treatment of various enzyme mechanisms, special approaches are often needed to obtain results, which adequately refer to experimental data. As a result of the spectacular progress in the last two decades, most enzyme reactions can be quite precisely treated by various computational methods.
Collapse
|
34
|
Abstract
Quantum tunneling and protein dynamics have emerged as important components of enzyme function. This review focuses on soybean lipoxygenase-1, to illustrate how the properties of enzymatic C-H bond activation link protein motions to the fundamental bond making-breaking processes.
Collapse
Affiliation(s)
- Judith P. Klinman
- Department of Chemistry, Department of Molecular and Cell Biology, and California Institute for Quantitative Biosciences (QB3), University of California, Berkeley, California 94720-3220, United States
| |
Collapse
|
35
|
Bushnell EAC, Jamil R, Gauld JW. Gaining insight into the chemistry of lipoxygenases: a computational investigation into the catalytic mechanism of (8R)-lipoxygenase. J Biol Inorg Chem 2013; 18:343-55. [PMID: 23361122 DOI: 10.1007/s00775-013-0978-4] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/01/2012] [Accepted: 01/14/2013] [Indexed: 10/27/2022]
Abstract
Lipoxygenases (LOXs) are ubiquitous in nature and catalyze a range of life-essential reactions within organisms. In particular they are critical to the formation of eicosanoids, which are critical for normal cell function. However, a number of important questions about the reactivity and mechanism of these enzymes still remain. Specifically, although the initial step in the mechanism of LOXs has been well studied, little is known of subsequent steps. Thus, with use of a quantum mechanical/molecular mechanical approach, the complete catalytic mechanism of (8R)-LOX was investigated. The results have provided a better understanding of the general chemistry of LOXs as a whole. In particular, from comparisons with soybean LOX-1, it appears that the initial proton-coupled electron transfer may be very similar among all LOXs. Furthermore, LOXs appear to undergo multistate reactivity where potential spin inversion of an electron may occur either in the attack of O(2) or in the regeneration of the active site. Lastly, it is shown that with the explicit modeling of the environment, the regeneration of the active center likely occurs via the rotation of the intermediate followed by an outer-sphere [Formula: see text] transfer as opposed to the formation of a "purple intermediate" complex.
Collapse
Affiliation(s)
- Eric A C Bushnell
- Department of Chemistry and Biochemistry, University of Windsor, Windsor, ON, N9B 3P4, Canada
| | | | | |
Collapse
|
36
|
Alonso JL, Castro A, Clemente-Gallardo J, Echenique P, Mazo JJ, Polo V, Rubio A, Zueco D. Non-adiabatic effects within a single thermally averaged potential energy surface: Thermal expansion and reaction rates of small molecules. J Chem Phys 2012; 137:22A533. [DOI: 10.1063/1.4747699] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
37
|
Schrauben JN, Cattaneo M, Day TC, Tenderholt AL, Mayer JM. Multiple-site concerted proton-electron transfer reactions of hydrogen-bonded phenols are nonadiabatic and well described by semiclassical Marcus theory. J Am Chem Soc 2012; 134:16635-45. [PMID: 22974135 PMCID: PMC3476473 DOI: 10.1021/ja305668h] [Citation(s) in RCA: 71] [Impact Index Per Article: 5.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/27/2022]
Abstract
Photo-oxidations of hydrogen-bonded phenols using excited-state polyarenes are described to derive fundamental understanding of multiple-site concerted proton-electron transfer reactions (MS-CPET). Experiments have examined phenol bases having -CPh(2)NH(2), -Py, and -CH(2)Py groups ortho to the phenol hydroxyl group and tert-butyl groups in the 4,6-positions for stability (HOAr-NH(2), HOAr-Py, and HOAr-CH(2)Py, respectively; Py = pyridyl; Ph = phenyl). The photo-oxidations proceed by intramolecular proton transfer from the phenol to the pendent base concerted with electron transfer to the excited polyarene. For comparison, 2,4,6-(t)Bu(3)C(6)H(2)OH, a phenol without a pendent base and tert-butyl groups in the 2,4,6-positions, has also been examined. Many of these bimolecular reactions are fast, with rate constants near the diffusion limit. Combining the photochemical k(CPET) values with those from prior thermal stopped-flow kinetic studies gives data sets for the oxidations of HOAr-NH(2) and HOAr-CH(2)Py that span over 10(7) in k(CPET) and nearly 0.9 eV in driving force (ΔG(o)'). Plots of log(k(CPET)) vs ΔG(o)', including both excited-state anthracenes and ground state aminium radical cations, define a single Marcus parabola in each case. These two data sets are thus well described by semiclassical Marcus theory, providing a strong validation of the use of this theory for MS-CPET. The parabolas give λ(CPET) ≅ 1.15-1.2 eV and H(ab) ≅ 20-30 cm(-1). These experiments represent the most direct measurements of H(ab) for MS-CPET reactions to date. Although rate constants are available only up to the diffusion limit, the parabolas clearly peak well below the adiabatic limit of ca. 6 × 10(12) s(-1). Thus, this is a very clear demonstration that the reactions are nonadiabatic. The nonadiabatic character slows the reactions by a factor of ~45. Results for the oxidation of HOAr-Py, in which the phenol and base are conjugated, and for oxidation of 2,4,6-(t)Bu(3)C(6)H(2)OH, which lacks a base, show that both have substantially lower λ and larger pre-exponential terms. The implications of these results for MS-CPET reactions are discussed.
Collapse
Affiliation(s)
| | | | - Thomas C. Day
- Department of Chemistry, University of Washington, Seattle WA 98195
| | | | - James M. Mayer
- Department of Chemistry, University of Washington, Seattle WA 98195
| |
Collapse
|
38
|
Trakhtenberg LI, Nadtochenko VA. Tunneling proton transfer in biological systems. Role of temperature and pressure. RUSSIAN JOURNAL OF PHYSICAL CHEMISTRY A 2012. [DOI: 10.1134/s0036024412090142] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/22/2022]
|
39
|
Phatak P, Sumner I, Iyengar SS. Gauging the flexibility of the active site in soybean lipoxygenase-1 (SLO-1) through an atom-centered density matrix propagation (ADMP) treatment that facilitates the sampling of rare events. J Phys Chem B 2012; 116:10145-64. [PMID: 22838384 PMCID: PMC3558621 DOI: 10.1021/jp3015047] [Citation(s) in RCA: 25] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
We present a computational methodology to sample rare events in large biological enzymes that may involve electronically polarizing, reactive processes. The approach includes simultaneous dynamical treatment of electronic and nuclear degrees of freedom, where contributions from the electronic portion are computed using hybrid density functional theory and the computational costs are reduced through a hybrid quantum mechanics/molecular mechanics (QM/MM) treatment. Thus, the paper involves a QM/MM dynamical treatment of rare events. The method is applied to probe the effect of the active site elements on the critical hydrogen transfer step in the soybean lipoxygenase-1 (SLO-1) catalyzed oxidation of linoleic acid. It is found that the dynamical fluctuations and associated flexibility of the active site are critical toward maintaining the electrostatics in the regime where the reactive process can occur smoothly. Physical constraints enforced to limit the active site flexibility are akin to mutations and, in the cases studied, have a detrimental effect on the electrostatic fluctuations, thus adversely affecting the hydrogen transfer process.
Collapse
Affiliation(s)
- Prasad Phatak
- Department of Chemistry and Department of Physics, Indiana University, 800 E. Kirkwood Ave, Bloomington, IN-47405
| | - Isaiah Sumner
- Department of Chemistry and Department of Physics, Indiana University, 800 E. Kirkwood Ave, Bloomington, IN-47405
| | - Srinivasan S. Iyengar
- Department of Chemistry and Department of Physics, Indiana University, 800 E. Kirkwood Ave, Bloomington, IN-47405
| |
Collapse
|
40
|
Mazzuca J, Garashchuk S, Jakowski J. Description of proton transfer in soybean lipoxygenase-1 employing approximate quantum trajectory dynamics. Chem Phys Lett 2012. [DOI: 10.1016/j.cplett.2012.06.019] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/01/2022]
|
41
|
Weinberg DR, Gagliardi CJ, Hull JF, Murphy CF, Kent CA, Westlake BC, Paul A, Ess DH, McCafferty DG, Meyer TJ. Proton-Coupled Electron Transfer. Chem Rev 2012; 112:4016-93. [DOI: 10.1021/cr200177j] [Citation(s) in RCA: 1125] [Impact Index Per Article: 93.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Affiliation(s)
- David R. Weinberg
- Department
of Chemistry, University
of North Carolina at Chapel Hill, Chapel Hill, North Carolina 27599-3290,
United States
- Department of Physical and Environmental
Sciences, Colorado Mesa University, 1100 North Avenue, Grand Junction,
Colorado 81501-3122, United States
| | - Christopher J. Gagliardi
- Department
of Chemistry, University
of North Carolina at Chapel Hill, Chapel Hill, North Carolina 27599-3290,
United States
| | - Jonathan F. Hull
- Department
of Chemistry, University
of North Carolina at Chapel Hill, Chapel Hill, North Carolina 27599-3290,
United States
| | - Christine Fecenko Murphy
- Department
of Chemistry, B219
Levine Science Research Center, Box 90354, Duke University, Durham,
North Carolina 27708-0354, United States
| | - Caleb A. Kent
- Department
of Chemistry, University
of North Carolina at Chapel Hill, Chapel Hill, North Carolina 27599-3290,
United States
| | - Brittany C. Westlake
- The American Chemical Society,
1155 Sixteenth Street NW, Washington, District of Columbia 20036,
United States
| | - Amit Paul
- Department
of Chemistry, University
of North Carolina at Chapel Hill, Chapel Hill, North Carolina 27599-3290,
United States
| | - Daniel H. Ess
- Department
of Chemistry, University
of North Carolina at Chapel Hill, Chapel Hill, North Carolina 27599-3290,
United States
| | - Dewey Granville McCafferty
- Department
of Chemistry, B219
Levine Science Research Center, Box 90354, Duke University, Durham,
North Carolina 27708-0354, United States
| | - Thomas J. Meyer
- Department
of Chemistry, University
of North Carolina at Chapel Hill, Chapel Hill, North Carolina 27599-3290,
United States
| |
Collapse
|
42
|
Pacheco AB, Dietrick SM, Stevens PS, Iyengar SS. "Pump-probe" atom-centered density matrix propagation studies to gauge anharmonicity and energy repartitioning in atmospheric reactive adducts: case study of the OH + isoprene and OH + butadiene reaction intermediates. J Phys Chem A 2012; 116:4108-28. [PMID: 22401490 DOI: 10.1021/jp212330e] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Time-resolved "pump-probe" ab initio molecular dynamics studies are constructed to probe the stability of reaction intermediates, the mechanism of energy transfer, and energy repartitioning, for moieties involved during the interaction of volatile organic compunds with hydroxyl radical. These systems are of prime importance in the atmosphere. Specifically, the stability of reaction intermediates of hydroxyl radical adducts to isoprene and butadiene molecules is used as a case study to develop novel computational techniques involving "pump-probe" ab initio molecular dynamics. Starting with the various possible hydroxyl radical adducts to isoprene and butadiene, select vibrational modes of each of the adducts are populated with excess energy to mimic the initial conditions of an experiment. The flow of energy into the remaining modes is then probed by subjecting the excited adducts to ab initio molecular dynamics simulations. It is found that the stability of the adducts arises directly due to the anhormonically driven coupling of the modes to facilitate repartitioning of the excess vibrational energy. This kind of vibrational repartitioning has a critical influence on the energy density.
Collapse
Affiliation(s)
- Alexander B Pacheco
- Department of Chemistry, Indiana University, Bloomington, Indiana 47405, USA
| | | | | | | |
Collapse
|
43
|
Brela M, Stare J, Pirc G, Sollner-Dolenc M, Boczar M, Wójcik MJ, Mavri J. Car-Parrinello simulation of the vibrational spectrum of a medium strong hydrogen bond by two-dimensional quantization of the nuclear motion: application to 2-hydroxy-5-nitrobenzamide. J Phys Chem B 2012; 116:4510-8. [PMID: 22429110 DOI: 10.1021/jp2094559] [Citation(s) in RCA: 52] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023]
Abstract
The nature of medium strong intra- and intermolecular hydrogen bonding in 2-hydroxy-5-nitrobenzamide in the crystal phase was examined by infrared spectroscopy and Car-Parrinello molecular dynamics simulation. The focal point of our study was the part of the infrared spectra associated with the O-H and N-H stretching modes that are very sensitive to the strength of hydrogen bonding. For spectra calculations we used an isolated dimer and the fully periodic crystal environment. We calculated the spectra by using harmonic approximation, the time course of the dipole moment function as obtained from the Car-Parrinello simulation, and the quantization of the nuclear motion of the proton for an instantaneous snapshot of the structures in one and two dimensions. Although quantitative assessment of the agreement between the computed and experimental band contour is difficult due to the fact that the experimental band is very broad, we feel that the most reasonable qualitative agreement with the experiment is obtained from snapshot structures and two-dimensional quantization of the proton motion. We have also critically examined the methods of constructing the one-dimensional proton potential. Perspectives are given for the treatment of nuclear quantum effects in biocatalysis.
Collapse
Affiliation(s)
- Mateusz Brela
- Laboratory for Biocomputing and Bioinformatics, National Institute of Chemistry, Ljubljana, Slovenia
| | | | | | | | | | | | | |
Collapse
|
44
|
Markle TF, Tenderholt AL, Mayer JM. Probing quantum and dynamic effects in concerted proton-electron transfer reactions of phenol-base compounds. J Phys Chem B 2012; 116:571-84. [PMID: 22148459 PMCID: PMC3974916 DOI: 10.1021/jp2091736] [Citation(s) in RCA: 34] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
The oxidation of three phenols, which contain an intramolecular hydrogen bond to a pendent pyridine or amine group, has been shown, in a previous experimental study, to undergo concerted proton-electron transfer (CPET). In this reaction, the electron is transferred to an outer-sphere oxidant, and the proton is transferred from the oxygen to nitrogen atom. In the present study, this reaction is studied computationally using a version of Hammes-Schiffer's multistate continuum theory where CPET is formulated as a transmission frequency between neutral and cation vibrational-electronic states. The neutral and cation proton vibrational wave functions are computed from one-dimensional potential energy surfaces (PESs) for the transferring proton in a fixed heavy atom framework. The overlap integrals for these neutral/cation wave functions, considering several initial (i.e., neutral) and final (i.e., cation) vibrational states, are used to evaluate the relative rates of oxidation. The analysis is extended to heavy atom configurations with various proton donor-acceptor (i.e., O-N) distances to assess the importance of heavy atom "gating". Such changes in d(ON) dramatically affect the nature of the proton PESs and wave functions. Surprisingly, the most reactive configurations have similar donor-acceptor distances despite the large (~0.2 Å) differences in the optimized structures. These theoretical results qualitatively reproduce the experimental faster reactivity of the reaction of the pyridyl derivative 1 versus the CH(2)-pyridyl 2, but the computed factor of 5 is smaller than the experimental 10(2). The amine derivative is calculated to react similarly to 1, which does not agree with the experiments, likely due to some of the simplifying assumptions made in applying the theory. The computed kinetic isotope effects (KIEs) and their temperature dependence are in agreement with experimental results.
Collapse
Affiliation(s)
| | - Adam L. Tenderholt
- Department of Chemistry, University of Washington, Box 351700 Seattle, WA 98195-1700
| | - James M. Mayer
- Department of Chemistry, University of Washington, Box 351700 Seattle, WA 98195-1700
| |
Collapse
|
45
|
Dietrick SM, Pacheco AB, Phatak P, Stevens PS, Iyengar SS. Influence of Water on Anharmonicity, Stability, and Vibrational Energy Distribution of Hydrogen-Bonded Adducts in Atmospheric Reactions: Case Study of the OH + Isoprene Reaction Intermediate Using Ab Initio Molecular Dynamics. J Phys Chem A 2011; 116:399-414. [DOI: 10.1021/jp204511v] [Citation(s) in RCA: 25] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Scott M. Dietrick
- Department of Chemistry, Indiana University, Bloomington, Indiana 47405, United States
| | - Alexander B. Pacheco
- Department of Chemistry, Indiana University, Bloomington, Indiana 47405, United States
| | - Prasad Phatak
- Department of Chemistry, Indiana University, Bloomington, Indiana 47405, United States
| | - Philip S. Stevens
- Department of Chemistry and Center for Research in Environmental Science, School of Public and Environmental Affairs, Indiana University, Bloomington, Indiana 47405, United States
| | - Srinivasan S. Iyengar
- Departments of Chemistry and Physics, Indiana University, Bloomington, Indiana 47405, United States
| |
Collapse
|
46
|
Li X, Iyengar SS. Quantum Wavepacket Ab Initio Molecular Dynamics for Extended Systems. J Phys Chem A 2011; 115:6269-84. [DOI: 10.1021/jp112389m] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Xiaohu Li
- Department of Chemistry and Department of Physics, Indiana University, 800 East Kirkwood Avenue, Bloomington, IN 47405
| | - Srinivasan S. Iyengar
- Department of Chemistry and Department of Physics, Indiana University, 800 East Kirkwood Avenue, Bloomington, IN 47405
| |
Collapse
|
47
|
Azuri A, Engel H, Doron D, Major DT. Path-Integral Calculations of Nuclear Quantum Effects in Model Systems, Small Molecules, and Enzymes via Gradient-Based Forward Corrector Algorithms. J Chem Theory Comput 2011; 7:1273-86. [DOI: 10.1021/ct100716c] [Citation(s) in RCA: 29] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Asaf Azuri
- Department of Chemistry and the Lise Meitner-Minerva Center of Computational Quantum Chemistry, Bar-Ilan University, Ramat-Gan 52900, Israel
| | - Hamutal Engel
- Department of Chemistry and the Lise Meitner-Minerva Center of Computational Quantum Chemistry, Bar-Ilan University, Ramat-Gan 52900, Israel
| | - Dvir Doron
- Department of Chemistry and the Lise Meitner-Minerva Center of Computational Quantum Chemistry, Bar-Ilan University, Ramat-Gan 52900, Israel
| | - Dan Thomas Major
- Department of Chemistry and the Lise Meitner-Minerva Center of Computational Quantum Chemistry, Bar-Ilan University, Ramat-Gan 52900, Israel
| |
Collapse
|
48
|
Li X, Iyengar SS. Quantum wavepacket ab initio molecular dynamics: generalizations using an extended Lagrangian treatment of diabatic states coupled through multireference electronic structure. J Chem Phys 2011; 133:184105. [PMID: 21073211 DOI: 10.1063/1.3504167] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We present a generalization to our previously developed quantum wavepacket ab initio molecular dynamics (QWAIMD) method by using multiple diabatic electronic reduced single particle density matrices, propagated within an extended Lagrangian paradigm. The Slater determinantal wavefunctions associated with the density matrices utilized may be orthogonal or nonorthogonal with respect to each other. This generalization directly results from an analysis of the variance in electronic structure with quantum nuclear degrees of freedom. The diabatic electronic states are treated here as classical parametric variables and propagated simultaneously along with the quantum wavepacket and classical nuclei. Each electronic density matrix is constrained to be N-representable. Consequently two sets of new methods are derived: extended Lagrangian-QWAIMD (xLag-QWAIMD) and diabatic extended Lagrangian-QWAIMD (DxLag-QWAIMD). In both cases, the instantaneous potential energy surface for the quantum nuclear degrees of freedom is constructed from the diabatic states using an on-the-fly nonorthogonal multireference formalism. By introducing generalized grid-based electronic basis functions, we eliminate the basis set dependence on the quantum nucleus. Subsequent reuse of the two-electron integrals during the on-the-fly potential energy surface computation stage yields a substantial reduction in computational costs. Specifically, both xLag-QWAIMD and DxLag-QWAIMD turn out to be about two orders of magnitude faster than our previously developed time-dependent deterministic sampling implementation of QWAIMD. Energy conservation properties, accuracy of the associated potential surfaces, and vibrational properties are analyzed for a family of hydrogen bonded systems.
Collapse
Affiliation(s)
- Xiaohu Li
- Department of Chemistry, Indiana University, 800 E. Kirkwood Ave, Bloomington, Indiana 47405, USA
| | | |
Collapse
|
49
|
Hocker D, Li X, Iyengar SS. Shannon Entropy Based Time-Dependent Deterministic Sampling for Efficient “On-the-Fly” Quantum Dynamics and Electronic Structure. J Chem Theory Comput 2011; 7:256-68. [DOI: 10.1021/ct1005856] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Affiliation(s)
- David Hocker
- Department of Chemistry and Department of Physics, Indiana University, 800 E. Kirkwood Ave, Bloomington, Indiana 47405, United States
| | - Xiaohu Li
- Department of Chemistry and Department of Physics, Indiana University, 800 E. Kirkwood Ave, Bloomington, Indiana 47405, United States
| | - Srinivasan S. Iyengar
- Department of Chemistry and Department of Physics, Indiana University, 800 E. Kirkwood Ave, Bloomington, Indiana 47405, United States
| |
Collapse
|
50
|
McGeagh JD, Ranaghan KE, Mulholland AJ. Protein dynamics and enzyme catalysis: insights from simulations. BIOCHIMICA ET BIOPHYSICA ACTA-PROTEINS AND PROTEOMICS 2010; 1814:1077-92. [PMID: 21167324 DOI: 10.1016/j.bbapap.2010.12.002] [Citation(s) in RCA: 52] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/20/2010] [Revised: 11/25/2010] [Accepted: 12/03/2010] [Indexed: 10/18/2022]
Abstract
The role of protein dynamics in enzyme catalysis is one of the most active and controversial areas in enzymology today. Some researchers claim that protein dynamics are at the heart of enzyme catalytic efficiency, while others state that dynamics make no significant contribution to catalysis. What is the biochemist - or student - to make of the ferocious arguments in this area? Protein dynamics are complex and fascinating, as molecular dynamics simulations and experiments have shown. The essential question is: do these complex motions have functional significance? In particular, how do they affect or relate to chemical reactions within enzymes, and how are chemical and conformational changes coupled together? Biomolecular simulations can analyse enzyme reactions and dynamics in atomic detail, beyond that achievable in experiments: accurate atomistic modelling has an essential part to play in clarifying these issues. This article is part of a Special Issue entitled: Protein Dynamics: Experimental and Computational Approaches.
Collapse
Affiliation(s)
- John D McGeagh
- Centre for Computational Chemistry, School of Chemistry, University of Bristol, Cantock's Close, BS8 1TS, United Kingdom
| | | | | |
Collapse
|