1
|
Rivoire O. A role for conformational changes in enzyme catalysis. Biophys J 2024; 123:1563-1578. [PMID: 38704639 PMCID: PMC11213973 DOI: 10.1016/j.bpj.2024.04.030] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/24/2024] [Revised: 04/15/2024] [Accepted: 04/29/2024] [Indexed: 05/06/2024] Open
Abstract
The role played by conformational changes in enzyme catalysis is controversial. In addition to examining specific enzymes, studying formal models can help identify the conditions under which conformational changes promote catalysis. Here, we present a model demonstrating how conformational changes can break a generic trade-off due to the conflicting requirements of successive steps in catalytic cycles, namely high specificity for the transition state to accelerate the chemical transformation and low affinity for the products to favor their release. The mechanism by which the trade-off is broken is a transition between conformations with different affinities for the substrate. The role of the effector that induces the transition is played by a substrate "handle," a part of the substrate that is not chemically transformed but whose interaction with the enzyme is nevertheless essential to rapidly complete the catalytic cycle. A key element of the model is the formalization of the constraints causing the trade-off that the presence of multiple states breaks, which we attribute to the strong chemical similarity between successive reaction states-substrates, transition states, and products. For the sake of clarity, we present our model for irreversible one-step unimolecular reactions. In this context, we demonstrate how the different forms that chemical similarities between reaction states can take impose limits on the overall catalytic turnover. We first analyze catalysts without internal degrees of freedom and then show how two-state catalysts can overcome their limitations. Our results recapitulate previous proposals concerning the role of conformational changes and substrate handles in a formalism that makes explicit the constraints that elicit these features. In addition, our approach establishes links with studies in the field of heterogeneous catalysis, where the same trade-offs are observed and where overcoming them is a well-recognized challenge.
Collapse
|
2
|
Lei YK, Yagi K, Sugita Y. Learning QM/MM potential using equivariant multiscale model. J Chem Phys 2024; 160:214109. [PMID: 38828815 DOI: 10.1063/5.0205123] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/24/2024] [Accepted: 05/09/2024] [Indexed: 06/05/2024] Open
Abstract
The machine learning (ML) method emerges as an efficient and precise surrogate model for high-level electronic structure theory. Its application has been limited to closed chemical systems without considering external potentials from the surrounding environment. To address this limitation and incorporate the influence of external potentials, polarization effects, and long-range interactions between a chemical system and its environment, the first two terms of the Taylor expansion of an electrostatic operator have been used as extra input to the existing ML model to represent the electrostatic environments. However, high-order electrostatic interaction is often essential to account for external potentials from the environment. The existing models based only on invariant features cannot capture significant distribution patterns of the external potentials. Here, we propose a novel ML model that includes high-order terms of the Taylor expansion of an electrostatic operator and uses an equivariant model, which can generate a high-order tensor covariant with rotations as a base model. Therefore, we can use the multipole-expansion equation to derive a useful representation by accounting for polarization and intermolecular interaction. Moreover, to deal with long-range interactions, we follow the same strategy adopted to derive long-range interactions between a target system and its environment media. Our model achieves higher prediction accuracy and transferability among various environment media with these modifications.
Collapse
Affiliation(s)
- Yao-Kun Lei
- Theoretical Molecular Science Laboratory, RIKEN Cluster for Pioneering Research, Wako, Saitama 351-0198, Japan
- Computational Biophysics Research Team, RIKEN Center for Computational Science, Kobe, Hyogo 650-0047, Japan
- RIKEN Interdisciplinary Theoretical and Mathematical Sciences Program (iTHEMS), Wako, Saitama 351-0198, Japan
| | - Kiyoshi Yagi
- Theoretical Molecular Science Laboratory, RIKEN Cluster for Pioneering Research, Wako, Saitama 351-0198, Japan
- Computational Biophysics Research Team, RIKEN Center for Computational Science, Kobe, Hyogo 650-0047, Japan
| | - Yuji Sugita
- Theoretical Molecular Science Laboratory, RIKEN Cluster for Pioneering Research, Wako, Saitama 351-0198, Japan
- Computational Biophysics Research Team, RIKEN Center for Computational Science, Kobe, Hyogo 650-0047, Japan
- RIKEN Interdisciplinary Theoretical and Mathematical Sciences Program (iTHEMS), Wako, Saitama 351-0198, Japan
- Laboratory for Biomolecular Function Simulation, RIKEN Center for Biosystems Dynamics Research, Kobe, Hyogo 650-0047, Japan
| |
Collapse
|
3
|
Nam K, Shao Y, Major DT, Wolf-Watz M. Perspectives on Computational Enzyme Modeling: From Mechanisms to Design and Drug Development. ACS OMEGA 2024; 9:7393-7412. [PMID: 38405524 PMCID: PMC10883025 DOI: 10.1021/acsomega.3c09084] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 11/14/2023] [Revised: 01/15/2024] [Accepted: 01/19/2024] [Indexed: 02/27/2024]
Abstract
Understanding enzyme mechanisms is essential for unraveling the complex molecular machinery of life. In this review, we survey the field of computational enzymology, highlighting key principles governing enzyme mechanisms and discussing ongoing challenges and promising advances. Over the years, computer simulations have become indispensable in the study of enzyme mechanisms, with the integration of experimental and computational exploration now established as a holistic approach to gain deep insights into enzymatic catalysis. Numerous studies have demonstrated the power of computer simulations in characterizing reaction pathways, transition states, substrate selectivity, product distribution, and dynamic conformational changes for various enzymes. Nevertheless, significant challenges remain in investigating the mechanisms of complex multistep reactions, large-scale conformational changes, and allosteric regulation. Beyond mechanistic studies, computational enzyme modeling has emerged as an essential tool for computer-aided enzyme design and the rational discovery of covalent drugs for targeted therapies. Overall, enzyme design/engineering and covalent drug development can greatly benefit from our understanding of the detailed mechanisms of enzymes, such as protein dynamics, entropy contributions, and allostery, as revealed by computational studies. Such a convergence of different research approaches is expected to continue, creating synergies in enzyme research. This review, by outlining the ever-expanding field of enzyme research, aims to provide guidance for future research directions and facilitate new developments in this important and evolving field.
Collapse
Affiliation(s)
- Kwangho Nam
- Department
of Chemistry and Biochemistry, University
of Texas at Arlington, Arlington, Texas 76019, United States
| | - Yihan Shao
- Department
of Chemistry and Biochemistry, University
of Oklahoma, Norman, Oklahoma 73019-5251, United States
| | - Dan T. Major
- Department
of Chemistry and Institute for Nanotechnology & Advanced Materials, Bar-Ilan University, Ramat-Gan 52900, Israel
| | | |
Collapse
|
4
|
Hill TD, Basnet S, Lepird HH, Rightnowar BW, Moran SD. Anisotropic dynamics of an interfacial enzyme active site observed using tethered substrate analogs and ultrafast 2D IR spectroscopy. J Chem Phys 2023; 159:165101. [PMID: 37870142 PMCID: PMC10597647 DOI: 10.1063/5.0167991] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/16/2023] [Accepted: 09/29/2023] [Indexed: 10/24/2023] Open
Abstract
Enzymes accelerate the rates of biomolecular reactions by many orders of magnitude compared to bulk solution, and it is widely understood that this catalytic effect arises from a combination of polar pre-organization and electrostatic transition state stabilization. A number of recent reports have also implicated ultrafast (femtosecond-picosecond) timescale motions in enzymatic activity. However, complications arising from spatially-distributed disorder, the occurrence of multiple substrate binding modes, and the influence of hydration dynamics on solvent-exposed active sites still confound many experimental studies. Here we use ultrafast two-dimensional infrared (2D IR) spectroscopy and covalently-tethered substrate analogs to examine dynamical properties of the promiscuous Pyrococcus horikoshii ene-reductase (PhENR) active site in two binding configurations mimicking proposed "inactive" and "reactive" Michaelis complexes. Spectral diffusion measurements of aryl-nitrile substrate analogs reveal an end-to-end tradeoff between fast (sub-ps) and slow (>5 ps) motions. Fermi resonant aryl-azide analogs that sense interactions of coupled oscillators are described. Lineshape and quantum beat analyses of these probes reveal characteristics that correlate with aryl-nitrile frequency fluctuation correlation functions parameters, demonstrating that this anisotropy is an intrinsic property of the water-exposed active site, where countervailing gradients of fast dynamics and disorder in the reactant ground state are maintained near the hydration interface. Our results suggest several plausible factors leading to state-selective rate enhancement and promiscuity in PhENR. This study also highlights a strategy to detect perturbations to vibrational modes outside the transparent window of the mid-IR spectrum, which may be extended to other macromolecular systems.
Collapse
Affiliation(s)
| | - Sunil Basnet
- School of Chemical and Biomolecular Sciences, Southern Illinois University Carbondale, 1245 Lincoln Drive MC 4409, Carbondale, Illinois 62901, USA
| | - Hannah H. Lepird
- School of Chemical and Biomolecular Sciences, Southern Illinois University Carbondale, 1245 Lincoln Drive MC 4409, Carbondale, Illinois 62901, USA
| | - Blaze W. Rightnowar
- School of Chemical and Biomolecular Sciences, Southern Illinois University Carbondale, 1245 Lincoln Drive MC 4409, Carbondale, Illinois 62901, USA
| | - Sean D. Moran
- School of Chemical and Biomolecular Sciences, Southern Illinois University Carbondale, 1245 Lincoln Drive MC 4409, Carbondale, Illinois 62901, USA
| |
Collapse
|
5
|
Kulkarni M, Söderhjelm P. Free-Energy Landscape and Rate Estimation of the Aromatic Ring Flips in Basic Pancreatic Trypsin Inhibitors Using Metadynamics. J Chem Theory Comput 2023; 19:6605-6618. [PMID: 37698852 PMCID: PMC10569046 DOI: 10.1021/acs.jctc.3c00460] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/01/2023] [Indexed: 09/13/2023]
Abstract
Aromatic side chains (phenylalanine and tyrosine) of a protein flip by 180° around the Cβ-Cγ axis (χ2 dihedral of the side chain), producing two symmetry-equivalent states. The study of ring flip dynamics with nuclear magnetic resonance (NMR) experiments helps to understand local conformational fluctuations. Ring flips are categorized as slow (milliseconds and onward) or fast (nanoseconds to near milliseconds) based on timescales accessible to NMR experiments. In this study, we investigated the ability of the infrequent metadynamics approach to estimate the flip rate and discriminate between slow and fast ring flips for eight individual aromatic side chains (F4, Y10, Y21, F22, Y23, F33, Y35, and F45) of the basic pancreatic trypsin inhibitor. Well-tempered metadynamics simulations were performed to estimate the ring-flipping free-energy surfaces for all eight aromatic residues. The results indicate that χ2 as a standalone collective variable (CV) is not sufficient to obtain computationally consistent results. Inclusion of a complementary CV, such as χ1(Cα-Cβ), solved the problem for most residues and enabled us to classify fast and slow ring flips. This indicates the importance of librational motions in ring flips. Multiple pathways and mechanisms were observed for residues F4, Y10, and F22. Recrossing events were observed for residues F22 and F33, indicating a possible role of friction effects in ring flipping. The results demonstrate the successful application of infrequent metadynamics to estimate ring flip rates and identify certain limitations of the approach.
Collapse
Affiliation(s)
- Mandar Kulkarni
- Division of Biophysical Chemistry, Lund University, Chemical Center, 22100 Lund, Sweden
| | - Pär Söderhjelm
- Division of Biophysical Chemistry, Lund University, Chemical Center, 22100 Lund, Sweden
| |
Collapse
|
6
|
Martins-Costa MTC, Ruiz-López MF. Electrostatics and Chemical Reactivity at the Air-Water Interface. J Am Chem Soc 2023; 145:1400-1406. [PMID: 36622259 DOI: 10.1021/jacs.2c12089] [Citation(s) in RCA: 20] [Impact Index Per Article: 20.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/10/2023]
Abstract
It has been recently discovered that chemical reactions at aqueous interfaces can be orders of magnitude faster compared to conventional bulk phase reactions, but despite its wide-ranging implications, which extend from atmospheric to synthetic chemistry or technological applications, the phenomenon is still incompletely understood. The role of strong electric fields due to space asymmetry and the accumulation of ions at the interface has been claimed as a possible cause from some experiments, but the reorganization of the solvent around the reactive system should provide even greater additional electrostatic contributions that have not yet been analyzed. In this study, with the help of first-principles molecular dynamics simulations, we go deeper into this issue by a careful assessment of solvation electrostatics at the air-water interface. Our simulations confirm that electrostatic forces can indeed be a key factor in rate acceleration compared to bulk solution. Remarkably, the study reveals that the effect cannot simply be attributed to the magnitude of the local electric field and that the fluctuations of the full electrostatic potential resulting from unique dynamical behavior of the solvation shells at the interface must be accounted for. This finding paves the way for future applications of the phenomenon in organic synthesis, especially for charge transfer or redox reactions in thin films and microdroplets.
Collapse
Affiliation(s)
- Marilia T C Martins-Costa
- Laboratoire de Physique et Chimie Théoriques, UMR CNRS 7019, University of Lorraine, CNRS, BP 70239, 54506 Vandoeuvre-lès-Nancy, France
| | - Manuel F Ruiz-López
- Laboratoire de Physique et Chimie Théoriques, UMR CNRS 7019, University of Lorraine, CNRS, BP 70239, 54506 Vandoeuvre-lès-Nancy, France
| |
Collapse
|
7
|
Rapp C, Nidetzky B. Hydride Transfer Mechanism of Enzymatic Sugar Nucleotide C2 Epimerization Probed with a Loose-Fit CDP-Glucose Substrate. ACS Catal 2022; 12:6816-6830. [PMID: 35747200 PMCID: PMC9207888 DOI: 10.1021/acscatal.2c00257] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/14/2022] [Revised: 05/12/2022] [Indexed: 11/29/2022]
Abstract
![]()
Transient oxidation–reduction
through hydride transfer with
tightly bound NAD coenzyme is used by a large class of sugar nucleotide
epimerases to promote configurational inversion of carbon stereocenters
in carbohydrate substrates. A requirement for the epimerases to coordinate
hydride abstraction and re-addition with substrate rotation in the
binding pocket poses a challenge for dynamical protein conformational
selection linked to enzyme catalysis. Here, we studied the thermophilic
C2 epimerase from Thermodesulfatator atlanticus (TaCPa2E) in combination with a slow CDP-glucose
substrate (kcat ≈ 1.0 min–1; 60 °C) to explore the sensitivity of the enzymatic hydride
transfer toward environmental fluctuations affected by temperature
(20–80 °C). We determined noncompetitive primary kinetic
isotope effects (KIE) due to 2H at the glucose C2 and showed
that a normal KIE on the kcat (Dkcat) reflects isotope sensitivity of
the hydrogen abstraction to enzyme-NAD+ in a rate-limiting
transient oxidation. The Dkcat peaked at 40 °C was 6.1 and decreased to 2.1 at low (20 °C)
and 3.3 at high temperature (80 °C). The temperature profiles
for kcat with the 1H and 2H substrate showed a decrease in the rate below a dynamically
important breakpoint (∼40 °C), suggesting an equilibrium
shift to an impaired conformational landscape relevant for catalysis
in the low-temperature region. Full Marcus-like model fits of the
rate and KIE profiles provided evidence for a high-temperature reaction
via low-frequency conformational sampling associated with a broad
distribution of hydride donor–acceptor distances (long-distance
population centered at 3.31 ± 0.02 Å), only poorly suitable
for quantum mechanical tunneling. Collectively, dynamical characteristics
of TaCPa2E-catalyzed hydride transfer during transient
oxidation of CDP-glucose reveal important analogies to mechanistically
simpler enzymes such as alcohol dehydrogenase and dihydrofolate reductase.
A loose-fit substrate (in TaCPa2E) resembles structural
variants of these enzymes by extensive dynamical sampling to balance
conformational flexibility and catalytic efficiency.
Collapse
Affiliation(s)
- Christian Rapp
- Institute of Biotechnology and Biochemical Engineering, Graz University of Technology, Petersgasse 10-12/1, 8010 Graz, Austria
| | - Bernd Nidetzky
- Institute of Biotechnology and Biochemical Engineering, Graz University of Technology, Petersgasse 10-12/1, 8010 Graz, Austria
- Austrian Centre of Industrial Biotechnology (ACIB), Petersgasse 14, 8010 Graz, Austria
| |
Collapse
|
8
|
Prah A, Mavri J, Stare J. An electrostatic duel: subtle differences in the catalytic performance of monoamine oxidase A and B isoenzymes elucidated at the residue level using quantum computations. Phys Chem Chem Phys 2021; 23:26459-26467. [PMID: 34806105 DOI: 10.1039/d1cp03993h] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Abstract
The origin of the immense catalytic power of enzymes remains one of the biggest unresolved questions in biochemistry, with electrostatics being one of the main contenders. Herein, we report results that not only confirm that electrostatics is the driving force behind enzyme catalysis, but also that it is capable of tuning subtle differences in the catalytic performance between structurally similar enzymes, as demonstrated using the example of isoenzymes, monoamine oxidases A and B. Using our own computationally efficient multiscale model [A. Prah, et al., ACS Catal., 2019, 9, 1231] we analyzed the rate-limiting step of the reaction between phenylethylamine and both isoenzymes and deduced that the electrostatic environment provided by isoenzyme B has a perceivably higher catalytic influence on all the considered parameters of the reaction (energy barrier, charge transfer, dipole moment, and HOMO-LUMO gap). This is in full agreement with the available experimental kinetic data and with our own simulations of the reaction in question. In-depth analysis of individual amino acid contributions of both isoenzymes to the barrier (based on the interaction between the electric field provided by the enzyme and the dipole moment of the reacting moiety) shows that the majority of the difference between the isoenzymes can be attributed to a small number of sizable differences between the aligned amino acid pairs, whereas in most of the pairs the difference in contribution to the barrier is vanishingly small. These results suggest that electrostatics largely controls the substrate selectivity of enzymes and validates our approach as being capable of discerning fine nuances in the selectivity of structurally related isoenzymes.
Collapse
Affiliation(s)
- Alja Prah
- Theory Department, National Institute of Chemistry, Slovenia. .,University of Ljubljana, Faculty of Pharmacy, Slovenia
| | - Janez Mavri
- Theory Department, National Institute of Chemistry, Slovenia.
| | - Jernej Stare
- Theory Department, National Institute of Chemistry, Slovenia.
| |
Collapse
|
9
|
Maffucci I, Laage D, Stirnemann G, Sterpone F. Differences in thermal structural changes and melting between mesophilic and thermophilic dihydrofolate reductase enzymes. Phys Chem Chem Phys 2020; 22:18361-18373. [DOI: 10.1039/d0cp02738c] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/25/2023]
Abstract
The thermal resistance of two homolog enzymes is investigated, with an emphasis on their local stability and flexibility, and on the possible implications regarding their reactivity.
Collapse
Affiliation(s)
- Irene Maffucci
- CNRS Laboratoire de Biochimie Théorique
- Institut de Biologie Physico-Chimique
- PSL University
- Paris
- France
| | - Damien Laage
- PASTEUR
- Département de chimie
- École Normale Supérieure
- PSL University
- Sorbonne Université
| | - Guillaume Stirnemann
- CNRS Laboratoire de Biochimie Théorique
- Institut de Biologie Physico-Chimique
- PSL University
- Paris
- France
| | - Fabio Sterpone
- CNRS Laboratoire de Biochimie Théorique
- Institut de Biologie Physico-Chimique
- PSL University
- Paris
- France
| |
Collapse
|
10
|
Lei YK, Zhang J, Zhang Z, Gao YQ. Dynamic Electric Field Complicates Chemical Reactions in Solutions. J Phys Chem Lett 2019; 10:2991-2997. [PMID: 31094529 DOI: 10.1021/acs.jpclett.9b01038] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/09/2023]
Abstract
Chemical reactions can be strongly influenced by an external electric field (EEF), but because the EEF is often time-dependent and in case it does not adapt quickly enough to the reaction progress, especially during fast barrier crossing processes, dynamic effects could be important. Here we find that electrostatic interactions can reduce the height of the reaction barrier for a Claissen rearrangement reaction and accelerate the key motions for bonding. Meanwhile, strong electrostatic interactions can modify the barrier into an effective potential well, confining the system into the barrier until solvents adjust themselves to provide an appropriate EEF for charge redistribution. In this case, the otherwise concerted mechanism becomes a stepwise one. Consequently, the motion of solvents modulates the reaction dynamics and leads to heterogeneous reaction paths, even in a seemingly homogeneous aqueous solution. In addition, an excessive stabilization of transition state retards the barrier crossing process, making the thermodynamically favorable pathway less favored dynamically.
Collapse
Affiliation(s)
- Yao Kun Lei
- Institute of Theoretical and Computational Chemistry, College of Chemistry and Molecular Engineering , Peking University , Beijing 100871 , China
| | - Jun Zhang
- Institute of Theoretical and Computational Chemistry, College of Chemistry and Molecular Engineering , Peking University , Beijing 100871 , China
| | - Zhen Zhang
- Department of Physics , Tangshan Normal University , Tangshan 063000 , China
| | - Yi Qin Gao
- Institute of Theoretical and Computational Chemistry, College of Chemistry and Molecular Engineering , Peking University , Beijing 100871 , China
| |
Collapse
|
11
|
Delgado M, Görlich S, Longbotham JE, Scrutton NS, Hay S, Moliner V, Tuñón I. Convergence of theory and experiment on the role of preorganization, quantum tunneling and enzyme motions into flavoenzyme-catalyzed hydride transfer. ACS Catal 2019; 7:3190-3198. [PMID: 31157122 DOI: 10.1021/acscatal.7b00201] [Citation(s) in RCA: 28] [Impact Index Per Article: 5.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Hydride transfer is one of the most common reactions catalyzed by enzymatic systems and it has become an object of study due to possible significant quantum tunneling effects. In the present work, we provide a combination of theoretical QM/MM simulations and experimental measurements of the rate constants and kinetic isotopic effects (KIEs) for the hydride transfer reaction catalyzed by morphinone reductase, MR. Quantum mechanical tunneling coefficients, computed in the framework of variational transition-state theory, play a significant role in this reaction, reaching values of 23.8 ± 5.5 for the lightest isotopologue; one of the largest values reported for enzymatic systems. This prediction is supported by the agreement between the theoretically predicted rate constants and the corresponding experimental values. Simulations indicate that the role of protein motions can be satisfactorily described as equilibrium fluctuations along the reaction coordinate, in line with a high degree of preorganization displayed by this enzyme.
Collapse
Affiliation(s)
- Manuel Delgado
- Department
of Physical and Analytical Chemistry, University Jaume I, 12071 Castelló de la Plana, Spain
| | - Stefan Görlich
- Manchester
Institute of Biotechnology and School of Chemistry, The University of Manchester, 131 Princess Street, Manchester M1 7DN, United Kingdom
| | - James E. Longbotham
- Manchester
Institute of Biotechnology and School of Chemistry, The University of Manchester, 131 Princess Street, Manchester M1 7DN, United Kingdom
| | - Nigel S. Scrutton
- Manchester
Institute of Biotechnology and School of Chemistry, The University of Manchester, 131 Princess Street, Manchester M1 7DN, United Kingdom
| | - Sam Hay
- Manchester
Institute of Biotechnology and School of Chemistry, The University of Manchester, 131 Princess Street, Manchester M1 7DN, United Kingdom
| | - Vicent Moliner
- Department
of Physical and Analytical Chemistry, University Jaume I, 12071 Castelló de la Plana, Spain
- Department
of Chemistry, University of Bath, Bath BA2 7AY, United Kingdom
| | - Iñaki Tuñón
- Department
of Physical Chemistry, University of València, 46100 Burjassot, Spain
| |
Collapse
|
12
|
Prah A, Frančišković E, Mavri J, Stare J. Electrostatics as the Driving Force Behind the Catalytic Function of the Monoamine Oxidase A Enzyme Confirmed by Quantum Computations. ACS Catal 2019. [DOI: 10.1021/acscatal.8b04045] [Citation(s) in RCA: 25] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/23/2023]
Affiliation(s)
- Alja Prah
- Theory Department, National Institute of Chemistry, 1001 Ljubljana, Slovenia
- Faculty of Pharmacy, University of Ljubljana, 1000 Ljubljana, Slovenia
| | - Eric Frančišković
- Theory Department, National Institute of Chemistry, 1001 Ljubljana, Slovenia
| | - Janez Mavri
- Theory Department, National Institute of Chemistry, 1001 Ljubljana, Slovenia
| | - Jernej Stare
- Theory Department, National Institute of Chemistry, 1001 Ljubljana, Slovenia
| |
Collapse
|
13
|
Tuñón I, Williams IH. The transition state and cognate concepts. ADVANCES IN PHYSICAL ORGANIC CHEMISTRY 2019. [DOI: 10.1016/bs.apoc.2019.09.001] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/01/2022]
|
14
|
Garcia-Meseguer R, Carpenter BK. Re-Evaluating the Transition State for Reactions in Solution. European J Org Chem 2018. [DOI: 10.1002/ejoc.201800841] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022]
Affiliation(s)
| | - Barry K. Carpenter
- School of Chemistry; Cardiff University; CF10 3AT Cardiff United Kingdom
| |
Collapse
|
15
|
Roca M, Ruiz-Pernía JJ, Castillo R, Oliva M, Moliner V. Temperature dependence of dynamic, tunnelling and kinetic isotope effects in formate dehydrogenase. Phys Chem Chem Phys 2018; 20:25722-25737. [PMID: 30280169 DOI: 10.1039/c8cp04244f] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Abstract
The origin of the catalytic power of enzymes has been a question of debate for a long time. In this regard, the possible contribution of protein dynamics in enzymatic catalysis has become one of the most controversial topics. In the present work, the hydride transfer step in the formate dehydrogenase (FDH EC 1.2.1.2) enzyme is studied by means of molecular dynamic (MD) simulations with quantum mechanics/molecular mechanics (QM/MM) potentials in order to explore any correlation between dynamics, tunnelling effects and the rate constant. The temperature dependence of the kinetic isotope effects (KIEs), which is one of the few tests that can be studied by experiments and simulations to shed light on this debate, has been computed and the results have been compared with previous experimental data. The classical mechanical free energy barrier and the number of recrossing trajectories seem to be temperature-independent while the quantum vibrational corrections and the tunnelling effects are slightly temperature-dependent over the interval of 5-45 °C. The computed primary KIEs are in very good agreement with previous experimental data, being almost temperature-independent within the standard deviations. The modest dependence on the temperature is due to just the quantum vibrational correction contribution. These results, together with the analysis of the evolution of the collective variables such as the electrostatic potential or the electric field created by the protein on the key atoms involved in the reaction, confirm that while the protein is well preorganised, some changes take place along the reaction that favour the hydride transfer and the product release. Coordinates defining these movements are, in fact, part of the real reaction coordinate.
Collapse
Affiliation(s)
- Maite Roca
- Departament de Química Física i Analítica, Universitat Jaume I, 12071 Castellón, Spain.
| | | | - Raquel Castillo
- Departament de Química Física i Analítica, Universitat Jaume I, 12071 Castellón, Spain.
| | - Mónica Oliva
- Departament de Química Física i Analítica, Universitat Jaume I, 12071 Castellón, Spain.
| | - Vicent Moliner
- Departament de Química Física i Analítica, Universitat Jaume I, 12071 Castellón, Spain.
| |
Collapse
|
16
|
Affiliation(s)
- Valerie Vaissier Welborn
- Kenneth S. Pitzer Center for Theoretical Chemistry and Department of Chemistry, University of California, Berkeley, California 94720, United States
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| | - Teresa Head-Gordon
- Kenneth S. Pitzer Center for Theoretical Chemistry and Department of Chemistry, University of California, Berkeley, California 94720, United States
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
- Department of Chemical and Biomolecular Engineering and Department of Bioengineering, University of California, Berkeley, California 94720, United States
| |
Collapse
|
17
|
Molecular modeling of conformational dynamics and its role in enzyme evolution. Curr Opin Struct Biol 2018; 52:50-57. [PMID: 30205262 DOI: 10.1016/j.sbi.2018.08.004] [Citation(s) in RCA: 26] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/16/2018] [Accepted: 08/20/2018] [Indexed: 12/19/2022]
Abstract
With increasing computational power, biomolecular simulations have become an invaluable tool for understanding enzyme mechanisms and the origins of enzyme catalysis. More recently, computational studies have started to focus on understanding how enzyme activity itself evolves, both in terms of enhancing the native or new activities on existing enzyme scaffolds, or completely de novo on previously non-catalytic scaffolds. In this context, both experiment and molecular modeling provided strong evidence for an important role of conformational dynamics in the evolution of enzyme functions. This contribution will present a brief overview of the current state of the art for computationally exploring enzyme conformational dynamics in enzyme evolution, and, using several showcase studies, illustrate the ways molecular modeling can be used to shed light on how enzyme function evolves, at the most fundamental molecular level.
Collapse
|
18
|
Petrović D, Risso VA, Kamerlin SCL, Sanchez-Ruiz JM. Conformational dynamics and enzyme evolution. J R Soc Interface 2018; 15:20180330. [PMID: 30021929 PMCID: PMC6073641 DOI: 10.1098/rsif.2018.0330] [Citation(s) in RCA: 103] [Impact Index Per Article: 17.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/10/2018] [Accepted: 06/27/2018] [Indexed: 12/21/2022] Open
Abstract
Enzymes are dynamic entities, and their dynamic properties are clearly linked to their biological function. It follows that dynamics ought to play an essential role in enzyme evolution. Indeed, a link between conformational diversity and the emergence of new enzyme functionalities has been recognized for many years. However, it is only recently that state-of-the-art computational and experimental approaches are revealing the crucial molecular details of this link. Specifically, evolutionary trajectories leading to functional optimization for a given host environment or to the emergence of a new function typically involve enriching catalytically competent conformations and/or the freezing out of non-competent conformations of an enzyme. In some cases, these evolutionary changes are achieved through distant mutations that shift the protein ensemble towards productive conformations. Multifunctional intermediates in evolutionary trajectories are probably multi-conformational, i.e. able to switch between different overall conformations, each competent for a given function. Conformational diversity can assist the emergence of a completely new active site through a single mutation by facilitating transition-state binding. We propose that this mechanism may have played a role in the emergence of enzymes at the primordial, progenote stage, where it was plausibly promoted by high environmental temperatures and the possibility of additional phenotypic mutations.
Collapse
Affiliation(s)
- Dušan Petrović
- Department of Chemistry, BMC, Uppsala University, Box 576, 751 23 Uppsala, Sweden
| | - Valeria A Risso
- Departamento de Quimica Fisica, Facultad de Ciencias, University of Granada, 18071 Granada, Spain
| | | | - Jose M Sanchez-Ruiz
- Departamento de Quimica Fisica, Facultad de Ciencias, University of Granada, 18071 Granada, Spain
| |
Collapse
|
19
|
Zhang Q, Pan Z, Zhang L, Zhang R, Chen Z, Jin T, Wu T, Chen X, Zhuang W. Ion effect on the dynamics of water hydrogen bonding network: A theoretical and computational spectroscopy point of view. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2018. [DOI: 10.1002/wcms.1373] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/08/2022]
Affiliation(s)
- Qiang Zhang
- State Key Laboratory of Structural Chemistry, Fujian Institute of Research on the Structure of MatterChinese Academy of SciencesFuzhouChina
- Department of ChemistryBohai UniversityJinzhouChina
| | - Zhijun Pan
- State Key Laboratory of Structural Chemistry, Fujian Institute of Research on the Structure of MatterChinese Academy of SciencesFuzhouChina
| | - Lu Zhang
- State Key Laboratory of Structural Chemistry, Fujian Institute of Research on the Structure of MatterChinese Academy of SciencesFuzhouChina
| | - Ruiting Zhang
- School of Physics and Optoelectronic EngineeringXidian UniversityXi'anChina
| | - Zhening Chen
- State Key Laboratory of Structural Chemistry, Fujian Institute of Research on the Structure of MatterChinese Academy of SciencesFuzhouChina
| | - Tan Jin
- State Key Laboratory of Structural Chemistry, Fujian Institute of Research on the Structure of MatterChinese Academy of SciencesFuzhouChina
| | - Tianmin Wu
- State Key Laboratory of Structural Chemistry, Fujian Institute of Research on the Structure of MatterChinese Academy of SciencesFuzhouChina
| | - Xian Chen
- State Key Laboratory of Structural Chemistry, Fujian Institute of Research on the Structure of MatterChinese Academy of SciencesFuzhouChina
| | - Wei Zhuang
- State Key Laboratory of Structural Chemistry, Fujian Institute of Research on the Structure of MatterChinese Academy of SciencesFuzhouChina
| |
Collapse
|
20
|
Schramm VL, Schwartz SD. Promoting Vibrations and the Function of Enzymes. Emerging Theoretical and Experimental Convergence. Biochemistry 2018; 57:3299-3308. [PMID: 29608286 DOI: 10.1021/acs.biochem.8b00201] [Citation(s) in RCA: 46] [Impact Index Per Article: 7.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
A complete understanding of enzyme catalysis requires knowledge of both transition state features and the detailed motions of atoms that cause reactant molecules to form and traverse the transition state. The seeming intractability of the problem arises from the femtosecond lifetime of chemical transition states, preventing most experimental access. Computational chemistry is admirably suited to short time scale analysis but can be misled by inappropriate starting points or by biased assumptions. Kinetic isotope effects provide an experimental approach to transition state structure and a method for obtaining transition state analogues but, alone, do not inform how that transition state is reached. Enzyme structures with transition state analogues provide computational starting points near the transition state geometry. These well-conditioned starting points, combined with the unbiased computational method of transition path sampling, provide realistic atomistic motions involved in transition state formation and passage. In many, but not all, enzymatic systems, femtosecond local protein motions near the catalytic site are linked to transition state formation. These motions are not inherently revealed by most approaches of transition state theory, because transition state theory replaces dynamics with the statistics of the transition state. Experimental and theoretical convergence of the link between local catalytic site vibrational modes and catalysis comes from heavy atom ("Born-Oppenheimer") enzymes. Fully labeled and catalytic site local heavy atom labels perturb the probability of finding enzymatic transition states in ways that can be analyzed and predicted by transition path sampling. Recent applications of these experimental and computational approaches reveal how subpicosecond local catalytic site protein modes play important roles in creating the transition state.
Collapse
Affiliation(s)
- Vern L Schramm
- Department of Biochemistry , Albert Einstein College of Medicine , Bronx , New York 10461 , United States
| | - Steven D Schwartz
- Department of Chemistry and Biochemistry , University of Arizona , Tucson , Arizona 85721 , United States
| |
Collapse
|
21
|
Santos-Martins D, Calixto AR, Fernandes PA, Ramos MJ. A Buried Water Molecule Influences Reactivity in α-Amylase on a Subnanosecond Time Scale. ACS Catal 2018. [DOI: 10.1021/acscatal.7b04400] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Affiliation(s)
- Diogo Santos-Martins
- UCIBIO@REQUIMTE, Departamento de Química e Bioquímica, Faculdade de Ciências, Universidade do Porto, Rua do Campo Alegre s/n, 4169-007 Porto, Portugal
| | - Ana R. Calixto
- UCIBIO@REQUIMTE, Departamento de Química e Bioquímica, Faculdade de Ciências, Universidade do Porto, Rua do Campo Alegre s/n, 4169-007 Porto, Portugal
| | - Pedro A. Fernandes
- UCIBIO@REQUIMTE, Departamento de Química e Bioquímica, Faculdade de Ciências, Universidade do Porto, Rua do Campo Alegre s/n, 4169-007 Porto, Portugal
| | - Maria J. Ramos
- UCIBIO@REQUIMTE, Departamento de Química e Bioquímica, Faculdade de Ciências, Universidade do Porto, Rua do Campo Alegre s/n, 4169-007 Porto, Portugal
| |
Collapse
|
22
|
Roca M, Moliner V, Tuñón I. Origin of Enzymatic Kinetic Isotope Effects in Human Purine Nucleoside Phosphorylase. ACS Catal 2017. [DOI: 10.1021/acscatal.7b04199] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- M. Roca
- Departament
de Química Física i Analítica, Universitat Jaume I, 12071 Castelló, Spain
- Departament
de Química Física, Universitat de València, 46100 Burjassot, Spain
| | - V. Moliner
- Departament
de Química Física i Analítica, Universitat Jaume I, 12071 Castelló, Spain
| | - I. Tuñón
- Departament
de Química Física, Universitat de València, 46100 Burjassot, Spain
| |
Collapse
|
23
|
Zinovjev K, Tuñón I. Quantifying the limits of transition state theory in enzymatic catalysis. Proc Natl Acad Sci U S A 2017; 114:12390-12395. [PMID: 29101125 PMCID: PMC5703300 DOI: 10.1073/pnas.1710820114] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
While being one of the most popular reaction rate theories, the applicability of transition state theory to the study of enzymatic reactions has been often challenged. The complex dynamic nature of the protein environment raised the question about the validity of the nonrecrossing hypothesis, a cornerstone in this theory. We present a computational strategy to quantify the error associated to transition state theory from the number of recrossings observed at the equicommittor, which is the best possible dividing surface. Application of a direct multidimensional transition state optimization to the hydride transfer step in human dihydrofolate reductase shows that both the participation of the protein degrees of freedom in the reaction coordinate and the error associated to the nonrecrossing hypothesis are small. Thus, the use of transition state theory, even with simplified reaction coordinates, provides a good theoretical framework for the study of enzymatic catalysis.
Collapse
Affiliation(s)
- Kirill Zinovjev
- Departament de Química Física, Universitat de València, 46100 Burjassot, Spain
| | - Iñaki Tuñón
- Departament de Química Física, Universitat de València, 46100 Burjassot, Spain
| |
Collapse
|
24
|
Lee YV, Choi SB, Wahab HA, Choong YS. Active Site Flexibility of Mycobacterium tuberculosis Isocitrate Lyase in Dimer Form. J Chem Inf Model 2017; 57:2351-2357. [PMID: 28820943 DOI: 10.1021/acs.jcim.7b00265] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Tuberculosis (TB) still remains a global threat due to the emergence of a drug-resistant strain. Instead of focusing on the drug target of active stage TB, we are highlighting the isocitrate lyase (ICL) at the dormant stage TB. ICL is one of the persistent factors for Mycobacterium tuberculosis (MTB) to survive during the dormant phase. In addition, the absence of ICL in human has made ICL a potential drug target for TB therapy. However, the dynamic details of ICL which could give insights to the ICL-ligand interaction have yet to be solved. Therefore, a series of ICL dimer dynamics studies through molecular dynamics simulation were performed in this work. The ICL active site entrance gate closure is contributed to by hydrogen bonding and electrostatic interactions with the C-terminal. Analysis suggested that the open-closed behavior of the ICL active site entrance depends on the type of ligand present in the active site. We also observed four residues (Ser91, Asp108, Asp153, and Cys191) which could possibly be the nucleophiles for nucleophilic attack on the cleavage of isocitrate at the C2-C3 bond. We hope that the elucidation of ICL dynamics can benefit future works such as lead identification or antibody design against ICL for TB therapeutics.
Collapse
Affiliation(s)
- Yie-Vern Lee
- Institute for Research in Molecular Medicine (INFORMM), Universiti Sains Malaysia , 11800 Minden, Penang, Malaysia
| | - Sy Bing Choi
- Pharmaceutical Design and Simulation Laboratory, School of Pharmaceutical Sciences, Universiti Sains Malaysia , 11800 Minden, Penang, Malaysia
| | - Habibah A Wahab
- Pharmaceutical Design and Simulation Laboratory, School of Pharmaceutical Sciences, Universiti Sains Malaysia , 11800 Minden, Penang, Malaysia
| | - Yee Siew Choong
- Institute for Research in Molecular Medicine (INFORMM), Universiti Sains Malaysia , 11800 Minden, Penang, Malaysia
| |
Collapse
|
25
|
Duboué-Dijon E, Pluhařová E, Domin D, Sen K, Fogarty AC, Chéron N, Laage D. Coupled Valence-Bond State Molecular Dynamics Description of an Enzyme-Catalyzed Reaction in a Non-Aqueous Organic Solvent. J Phys Chem B 2017; 121:7027-7041. [PMID: 28675789 DOI: 10.1021/acs.jpcb.7b03102] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Enzymes are widely used in nonaqueous solvents to catalyze non-natural reactions. While experimental measurements showed that the solvent nature has a strong effect on the reaction kinetics, the molecular details of the catalytic mechanism in nonaqueous solvents have remained largely elusive. Here we study the transesterification reaction catalyzed by the paradigm subtilisin Carlsberg serine protease in an organic apolar solvent. The rate-limiting acylation step involves a proton transfer between active-site residues and the nucleophilic attack of the substrate to form a tetrahedral intermediate. We design the first coupled valence-bond state model that simultaneously describes both reactions in the enzymatic active site. We develop a new systematic procedure to parametrize this model on high-level ab initio QM/MM free energy calculations that account for the molecular details of the active site and for both substrate and protein conformational fluctuations. Our calculations show that the reaction energy barrier changes dramatically with the solvent and protein conformational fluctuations. We find that the mechanism of the tetrahedral intermediate formation during the acylation step is similar to that determined under aqueous conditions, and that the proton transfer and nucleophilic attack reactions occur concertedly. We identify the reaction coordinate to be mostly due to the rearrangement of some residual water molecules close to the active site.
Collapse
Affiliation(s)
- Elise Duboué-Dijon
- École Normale Supérieure - PSL Research University , Département de Chimie, Sorbonne Universités - UPMC Univ Paris 06, CNRS UMR 8640 PASTEUR, 24 rue Lhomond, 75005 Paris, France
| | - Eva Pluhařová
- École Normale Supérieure - PSL Research University , Département de Chimie, Sorbonne Universités - UPMC Univ Paris 06, CNRS UMR 8640 PASTEUR, 24 rue Lhomond, 75005 Paris, France
| | - Dominik Domin
- École Normale Supérieure - PSL Research University , Département de Chimie, Sorbonne Universités - UPMC Univ Paris 06, CNRS UMR 8640 PASTEUR, 24 rue Lhomond, 75005 Paris, France
| | - Kakali Sen
- École Normale Supérieure - PSL Research University , Département de Chimie, Sorbonne Universités - UPMC Univ Paris 06, CNRS UMR 8640 PASTEUR, 24 rue Lhomond, 75005 Paris, France
| | - Aoife C Fogarty
- École Normale Supérieure - PSL Research University , Département de Chimie, Sorbonne Universités - UPMC Univ Paris 06, CNRS UMR 8640 PASTEUR, 24 rue Lhomond, 75005 Paris, France
| | - Nicolas Chéron
- École Normale Supérieure - PSL Research University , Département de Chimie, Sorbonne Universités - UPMC Univ Paris 06, CNRS UMR 8640 PASTEUR, 24 rue Lhomond, 75005 Paris, France
| | - Damien Laage
- École Normale Supérieure - PSL Research University , Département de Chimie, Sorbonne Universités - UPMC Univ Paris 06, CNRS UMR 8640 PASTEUR, 24 rue Lhomond, 75005 Paris, France
| |
Collapse
|
26
|
Ojeda-May P, Nam K. Acceleration of Semiempirical QM/MM Methods through Message Passage Interface (MPI), Hybrid MPI/Open Multiprocessing, and Self-Consistent Field Accelerator Implementations. J Chem Theory Comput 2017. [PMID: 28628742 DOI: 10.1021/acs.jctc.7b00322] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
The strategy and implementation of scalable and efficient semiempirical (SE) QM/MM methods in CHARMM are described. The serial version of the code was first profiled to identify routines that required parallelization. Afterward, the code was parallelized and accelerated with three approaches. The first approach was the parallelization of the entire QM/MM routines, including the Fock matrix diagonalization routines, using the CHARMM message passage interface (MPI) machinery. In the second approach, two different self-consistent field (SCF) energy convergence accelerators were implemented using density and Fock matrices as targets for their extrapolations in the SCF procedure. In the third approach, the entire QM/MM and MM energy routines were accelerated by implementing the hybrid MPI/open multiprocessing (OpenMP) model in which both the task- and loop-level parallelization strategies were adopted to balance loads between different OpenMP threads. The present implementation was tested on two solvated enzyme systems (including <100 QM atoms) and an SN2 symmetric reaction in water. The MPI version exceeded existing SE QM methods in CHARMM, which include the SCC-DFTB and SQUANTUM methods, by at least 4-fold. The use of SCF convergence accelerators further accelerated the code by ∼12-35% depending on the size of the QM region and the number of CPU cores used. Although the MPI version displayed good scalability, the performance was diminished for large numbers of MPI processes due to the overhead associated with MPI communications between nodes. This issue was partially overcome by the hybrid MPI/OpenMP approach which displayed a better scalability for a larger number of CPU cores (up to 64 CPUs in the tested systems).
Collapse
Affiliation(s)
| | - Kwangho Nam
- Department of Chemistry and Biochemistry, University of Texas at Arlington , Arlington, Texas 76019-0065, United States
| |
Collapse
|
27
|
Warshel A, Bora RP. Perspective: Defining and quantifying the role of dynamics in enzyme catalysis. J Chem Phys 2017; 144:180901. [PMID: 27179464 DOI: 10.1063/1.4947037] [Citation(s) in RCA: 134] [Impact Index Per Article: 19.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022] Open
Abstract
Enzymes control chemical reactions that are key to life processes, and allow them to take place on the time scale needed for synchronization between the relevant reaction cycles. In addition to general interest in their biological roles, these proteins present a fundamental scientific puzzle, since the origin of their tremendous catalytic power is still unclear. While many different hypotheses have been put forward to rationalize this, one of the proposals that has become particularly popular in recent years is the idea that dynamical effects contribute to catalysis. Here, we present a critical review of the dynamical idea, considering all reasonable definitions of what does and does not qualify as a dynamical effect. We demonstrate that no dynamical effect (according to these definitions) has ever been experimentally shown to contribute to catalysis. Furthermore, the existence of non-negligible dynamical contributions to catalysis is not supported by consistent theoretical studies. Our review is aimed, in part, at readers with a background in chemical physics and biophysics, and illustrates that despite a substantial body of experimental effort, there has not yet been any study that consistently established a connection between an enzyme's conformational dynamics and a significant increase in the catalytic contribution of the chemical step. We also make the point that the dynamical proposal is not a semantic issue but a well-defined scientific hypothesis with well-defined conclusions.
Collapse
Affiliation(s)
- Arieh Warshel
- Department of Chemistry, University of Southern California, SGM 418, 3620 McClintock Avenue, Los Angeles, California 90089, USA
| | - Ram Prasad Bora
- Department of Chemistry, University of Southern California, SGM 418, 3620 McClintock Avenue, Los Angeles, California 90089, USA
| |
Collapse
|
28
|
Krzemińska A, Moliner V, Świderek K. Dynamic and Electrostatic Effects on the Reaction Catalyzed by HIV-1 Protease. J Am Chem Soc 2016; 138:16283-16298. [PMID: 27935692 DOI: 10.1021/jacs.6b06856] [Citation(s) in RCA: 58] [Impact Index Per Article: 7.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/11/2023]
Abstract
HIV-1 Protease (HIV-1 PR) is one of the three enzymes essential for the replication process of HIV-1 virus, which explains why it has been the main target for design of drugs against acquired immunodeficiency syndrome (AIDS). This work is focused on exploring the proteolysis reaction catalyzed by HIV-1 PR, with special attention to the dynamic and electrostatic effects governing its catalytic power. Free energy surfaces for all possible mechanisms have been computed in terms of potentials of mean force (PMFs) within hybrid QM/MM potentials, with the QM subset of atoms described at semiempirical (AM1) and DFT (M06-2X) level. The results suggest that the most favorable reaction mechanism involves formation of a gem-diol intermediate, whose decomposition into the product complex would correspond to the rate-limiting step. The agreement between the activation free energy of this step with experimental data, as well as kinetic isotope effects (KIEs), supports this prediction. The role of the protein dynamic was studied by protein isotope labeling in the framework of the Variational Transition State Theory. The predicted enzyme KIEs, also very close to the values measured experimentally, reveal a measurable but small dynamic effect. Our calculations show how the contribution of dynamic effects to the effective activation free energy appears to be below 1 kcal·mol-1. On the contrary, the electric field created by the protein in the active site of the enzyme emerges as being critical for the electronic reorganization required during the reaction. These electrostatic properties of the active site could be used as a mold for future drug design.
Collapse
Affiliation(s)
- Agnieszka Krzemińska
- Institute of Applied Radiation Chemistry, Faculty of Chemistry, Lodz University of Technology , Zeromskiego 116, 90-924 Lodz, Poland
| | - Vicent Moliner
- Departament de Química Física i Analítica, Universitat Jaume I , 12071 Castelló, Spain
| | - Katarzyna Świderek
- Institute of Applied Radiation Chemistry, Faculty of Chemistry, Lodz University of Technology , Zeromskiego 116, 90-924 Lodz, Poland.,Departament de Química Física i Analítica, Universitat Jaume I , 12071 Castelló, Spain
| |
Collapse
|
29
|
Zoi I, Suarez J, Antoniou D, Cameron SA, Schramm VL, Schwartz SD. Modulating Enzyme Catalysis through Mutations Designed to Alter Rapid Protein Dynamics. J Am Chem Soc 2016; 138:3403-9. [PMID: 26927977 PMCID: PMC4794390 DOI: 10.1021/jacs.5b12551] [Citation(s) in RCA: 41] [Impact Index Per Article: 5.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
The relevance of sub-picosecond protein motions to the catalytic event remains a topic of debate. Heavy enzymes (isotopically substituted) provide an experimental tool for bond-vibrational links to enzyme catalysis. A recent transition path sampling study with heavy purine nucleoside phosphorylase (PNP) characterized the experimentally observed mass-dependent slowing of barrier crossing (Antoniou, D.; Ge, X.; Schramm, V. L.; Schwartz, S. D. J. Phys. Chem. Lett. 2012, 3, 3538). Here we computationally identify second-sphere amino acid residues predicted to influence the freedom of the catalytic site vibrational modes linked to heavy enzyme effects in PNP. We mutated heavy and light PNPs to increase the catalytic site vibrational freedom. Enzymatic barrier-crossing rates were converted from mass-dependent to mass-independent as a result of the mutations. The mutagenic uncoupling of femtosecond motions between catalytic site groups and reactants decreased transition state barrier crossing by 2 orders of magnitude, an indication of the femtosecond dynamic contributions to catalysis.
Collapse
Affiliation(s)
- Ioanna Zoi
- Department of Chemistry and Biochemistry, University of Arizona, 1306 East University Blvd, Tucson, Arizona 85721, United States
| | - Javier Suarez
- Department of Biochemistry, Albert Einstein College of Medicine, 1300 Morris Park Ave, Bronx, New York 10461, United States
| | - Dimitri Antoniou
- Department of Chemistry and Biochemistry, University of Arizona, 1306 East University Blvd, Tucson, Arizona 85721, United States
| | - Scott A. Cameron
- Department of Biochemistry, Albert Einstein College of Medicine, 1300 Morris Park Ave, Bronx, New York 10461, United States
| | - Vern L. Schramm
- Department of Biochemistry, Albert Einstein College of Medicine, 1300 Morris Park Ave, Bronx, New York 10461, United States
| | - Steven D. Schwartz
- Department of Chemistry and Biochemistry, University of Arizona, 1306 East University Blvd, Tucson, Arizona 85721, United States
| |
Collapse
|
30
|
Antoniou D, Schwartz SD. Phase Space Bottlenecks in Enzymatic Reactions. J Phys Chem B 2016; 120:433-9. [PMID: 26756622 PMCID: PMC4734068 DOI: 10.1021/acs.jpcb.5b11157] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/13/2015] [Revised: 01/11/2016] [Indexed: 01/09/2023]
Abstract
The definition of a transition state on an individual reactive trajectory is made via a committor analysis. In the past, the bottleneck definition has often been applied in configuration space. This is an approximation, and in order to expand this definition, we are revisiting an enzyme in which we had identified a fast subpicosecond motion that makes the reaction possible. First we used a time-series analysis method to identify the exact time when this motion initiates donor-acceptor compression. Then we modified the standard committor analysis of transition path sampling to identify events in phase space and found that there is a dividing surface in phase space significantly earlier than the configurationally defined transition-state crossing.
Collapse
Affiliation(s)
- Dimitri Antoniou
- Department of Chemistry and Biochemistry, University of Arizona, 1306 East University Boulevard, Tucson, Arizona 85721, United States
| | - Steven D. Schwartz
- Department of Chemistry and Biochemistry, University of Arizona, 1306 East University Boulevard, Tucson, Arizona 85721, United States
| |
Collapse
|
31
|
Buß O, Jager S, Dold SM, Zimmermann S, Hamacher K, Schmitz K, Rudat J. Statistical Evaluation of HTS Assays for Enzymatic Hydrolysis of β-Keto Esters. PLoS One 2016; 11:e0146104. [PMID: 26730596 PMCID: PMC4711668 DOI: 10.1371/journal.pone.0146104] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/28/2015] [Accepted: 12/14/2015] [Indexed: 01/15/2023] Open
Abstract
β-keto esters are used as precursors for the synthesis of β-amino acids, which are building blocks for some classes of pharmaceuticals. Here we describe the comparison of screening procedures for hydrolases to be used for the hydrolysis of β-keto esters, the first step in the preparation of β-amino acids. Two of the tested high throughput screening (HTS) assays depend on coupled enzymatic reactions which detect the alcohol released during ester hydrolysis by luminescence or absorption. The third assay detects the pH shift due to acid formation using an indicator dye. To choose the most efficient approach for screening, we assessed these assays with different statistical methods—namely, the classical Z’-factor, standardized mean difference (SSMD), the Kolmogorov-Smirnov-test, and t-statistics. This revealed that all three assays are suitable for HTS, the pH assay performing best. Based on our data we discuss the explanatory power of different statistical measures. Finally, we successfully employed the pH assay to identify a very fast hydrolase in an enzyme-substrate screening.
Collapse
Affiliation(s)
- O. Buß
- Karlsruhe Institute of Technology, Technical Biology, Karlsruhe, Germany
- * E-mail: (OB); (SJ)
| | - S. Jager
- Technische Universität Darmstadt, Computational Biology and Simulation, Darmstadt, Germany
- * E-mail: (OB); (SJ)
| | - S. -M. Dold
- Karlsruhe Institute of Technology, Technical Biology, Karlsruhe, Germany
| | - S. Zimmermann
- Karlsruhe Institute of Technology, Biomolecular Separation Engineering, Karlsruhe, Germany
| | - K. Hamacher
- Technische Universität Darmstadt, Computational Biology and Simulation, Darmstadt, Germany
| | - K. Schmitz
- Technische Universität Darmstadt, Biological Chemistry, Darmstadt, Germany
| | - J. Rudat
- Karlsruhe Institute of Technology, Technical Biology, Karlsruhe, Germany
| |
Collapse
|
32
|
Matyushov DV. Protein electron transfer: is biology (thermo)dynamic? JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2015; 27:473001. [PMID: 26558324 DOI: 10.1088/0953-8984/27/47/473001] [Citation(s) in RCA: 43] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/05/2023]
Abstract
Simple physical mechanisms are behind the flow of energy in all forms of life. Energy comes to living systems through electrons occupying high-energy states, either from food (respiratory chains) or from light (photosynthesis). This energy is transformed into the cross-membrane proton-motive force that eventually drives all biochemistry of the cell. Life's ability to transfer electrons over large distances with nearly zero loss of free energy is puzzling and has not been accomplished in synthetic systems. The focus of this review is on how this energetic efficiency is realized. General physical mechanisms and interactions that allow proteins to fold into compact water-soluble structures are also responsible for a rugged landscape of energy states and a broad distribution of relaxation times. Specific to a protein as a fluctuating thermal bath is the protein-water interface, which is heterogeneous both dynamically and structurally. The spectrum of interfacial fluctuations is a consequence of protein's elastic flexibility combined with a high density of surface charges polarizing water dipoles into surface nanodomains. Electrostatics is critical to the protein function and the relevant questions are: (i) What is the spectrum of interfacial electrostatic fluctuations? (ii) Does the interfacial biological water produce electrostatic signatures specific to proteins? (iii) How is protein-mediated chemistry affected by electrostatics? These questions connect the fluctuation spectrum to the dynamical control of chemical reactivity, i.e. the dependence of the activation free energy of the reaction on the dynamics of the bath. Ergodicity is often broken in protein-driven reactions and thermodynamic free energies become irrelevant. Continuous ergodicity breaking in a dense spectrum of relaxation times requires using dynamically restricted ensembles to calculate statistical averages. When applied to the calculation of the rates, this formalism leads to the nonergodic activated kinetics, which extends the transition-state theory to dynamically dispersive media. Releasing the grip of thermodynamics in kinetic calculations through nonergodicity provides the mechanism for an efficient optimization between reaction rates and the spectrum of relaxation times of the protein-water thermal bath. Bath dynamics, it appears, play as important role as the free energy in optimizing biology's performance.
Collapse
Affiliation(s)
- Dmitry V Matyushov
- Department of Physics and School of Molecular Sciences, Arizona State University, PO Box 871504, Tempe, AZ 85287-1504, USA
| |
Collapse
|