1
|
Wang P, Zhang J, Zhang S, Lu D, Zhu Y. Using High-Throughput Molecular Dynamics Simulation to Enhance the Computational Design of Kemp Elimination Enzymes. J Chem Inf Model 2023; 63:1323-1337. [PMID: 36782360 DOI: 10.1021/acs.jcim.3c00002] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/15/2023]
Abstract
Computational enzyme design has been successfully applied to identify new alternatives to natural enzymes for the biosynthesis of important compounds. However, the moderate catalytic activities of de novo designed enzymes indicate that the modeling accuracy of current computational enzyme design methods should be improved. Here, high-throughput molecular dynamics simulations were used to enhance computational enzyme design, thus allowing the identification of variants with higher activities in silico. Different time schemes of high-throughput molecular dynamics simulations were tested to identify the catalytic features of evolved Kemp eliminases. The 20 × 1 ns molecular dynamics simulation scheme was sufficiently accurate and computationally viable to screen the computationally designed massive variants of Kemp elimination enzymes. The developed hybrid computational strategy was used to redesign the most active Kemp eliminase, HG3.17, and five variants were generated and experimentally confirmed to afford higher catalytic efficiencies than that of HG3.17, with one double variant (D52Q/A53S) exhibiting a 55% increase. The hybrid computational enzyme design strategy is general and computationally economical, with which we anticipate the efficient creation of practical enzymes for industrial biocatalysis.
Collapse
Affiliation(s)
- Pengyu Wang
- College of Life Science and Technology, Beijing University of Chemical Technology, Beijing 100029, China.,Department of Chemical Engineering, Tsinghua University, Beijing 100084, China
| | - Jun Zhang
- Department of Chemical Engineering, Tsinghua University, Beijing 100084, China
| | - Shengyu Zhang
- College of Life Science and Technology, Beijing University of Chemical Technology, Beijing 100029, China
| | - Diannan Lu
- Department of Chemical Engineering, Tsinghua University, Beijing 100084, China
| | - Yushan Zhu
- College of Life Science and Technology, Beijing University of Chemical Technology, Beijing 100029, China
| |
Collapse
|
2
|
Tomasini M, Caporaso L, Trouvé J, Poater J, Gramage‐Doria R, Poater A. Unravelling Enzymatic Features in a Supramolecular Iridium Catalyst by Computational Calculations. Chemistry 2022; 28:e202201970. [PMID: 35788999 PMCID: PMC9804516 DOI: 10.1002/chem.202201970] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/26/2022] [Indexed: 01/05/2023]
Abstract
Non-biological catalysts following the governing principles of enzymes are attractive systems to disclose unprecedented reactivities. Most of those existing catalysts feature an adaptable molecular recognition site for substrate binding that are prone to undergo conformational selection pathways. Herein, we present a non-biological catalyst that is able to bind substrates via the induced fit model according to in-depth computational calculations. The system, which is constituted by an inflexible substrate-recognition site derived from a zinc-porphyrin in the second coordination sphere, features destabilization of ground states as well as stabilization of transition states for the relevant iridium-catalyzed C-H bond borylation of pyridine. In addition, this catalyst appears to be most suited to tightly bind the transition state rather than the substrate. Besides these features, which are reminiscent of the action modes of enzymes, new elementary catalytic steps (i. e. C-B bond formation and catalyst regeneration) have been disclosed owing to the unique distortions encountered in the different intermediates and transition states.
Collapse
Affiliation(s)
- Michele Tomasini
- Institut de Química Computacional i CatàlisiDepartament de QuímicaUniversitat de Gironac/Mª Aurèlia Capmany 6917003GironaCataloniaSpain,Department of ChemistryUniversity of SalernoVia Ponte Don Melillo84084FiscianoItaly
| | - Lucia Caporaso
- Department of ChemistryUniversity of SalernoVia Ponte Don Melillo84084FiscianoItaly
| | | | - Jordi Poater
- Departament de Química Inorgànica i Orgànica & IQTCUBUniversitat de Barcelona08028BarcelonaSpain,ICREA08010BarcelonaSpain
| | | | - Albert Poater
- Institut de Química Computacional i CatàlisiDepartament de QuímicaUniversitat de Gironac/Mª Aurèlia Capmany 6917003GironaCataloniaSpain
| |
Collapse
|
3
|
Natural Evolution Provides Strong Hints about Laboratory Evolution of Designer Enzymes. Proc Natl Acad Sci U S A 2022; 119:e2207904119. [PMID: 35901204 PMCID: PMC9351539 DOI: 10.1073/pnas.2207904119] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/30/2023] Open
Abstract
Laboratory evolution combined with computational enzyme design provides the opportunity to generate novel biocatalysts. Nevertheless, it has been challenging to understand how laboratory evolution optimizes designer enzymes by introducing seemingly random mutations. A typical enzyme optimized with laboratory evolution is the abiological Kemp eliminase, initially designed by grafting active site residues into a natural protein scaffold. Here, we relate the catalytic power of laboratory-evolved Kemp eliminases to the statistical energy ([Formula: see text]) inferred from their natural homologous sequences using the maximum entropy model. The [Formula: see text] of designs generated by directed evolution is correlated with enhanced activity and reduced stability, thus displaying a stability-activity trade-off. In contrast, the [Formula: see text] for mutants in catalytic-active remote regions (in which remote residues are important for catalysis) is strongly anticorrelated with the activity. These findings provide an insight into the role of protein scaffolds in the adaption to new enzymatic functions. It also indicates that the valley in the [Formula: see text] landscape can guide enzyme design for abiological catalysis. Overall, the connection between laboratory and natural evolution contributes to understanding what is optimized in the laboratory and how new enzymatic function emerges in nature, and provides guidance for computational enzyme design.
Collapse
|
4
|
Evolution of dynamical networks enhances catalysis in a designer enzyme. Nat Chem 2021; 13:1017-1022. [PMID: 34413499 DOI: 10.1038/s41557-021-00763-6] [Citation(s) in RCA: 34] [Impact Index Per Article: 11.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/17/2020] [Accepted: 06/30/2021] [Indexed: 02/08/2023]
Abstract
Activation heat capacity is emerging as a crucial factor in enzyme thermoadaptation, as shown by the non-Arrhenius behaviour of many natural enzymes. However, its physical origin and relationship to the evolution of catalytic activity remain uncertain. Here we show that directed evolution of a computationally designed Kemp eliminase reshapes protein dynamics, which gives rise to an activation heat capacity absent in the original design. These changes buttress transition-state stabilization. Extensive molecular dynamics simulations show that evolution results in the closure of solvent-exposed loops and a better packing of the active site. Remarkably, this gives rise to a correlated dynamical network that involves the transition state and large parts of the protein. This network tightens the transition-state ensemble, which induces a negative activation heat capacity and non-linearity in the activity-temperature dependence. Our results have implications for understanding enzyme evolution and suggest that selectively targeting the conformational dynamics of the transition-state ensemble by design and evolution will expedite the creation of novel enzymes.
Collapse
|
5
|
Zhang J, Kong K, Li X, Zhang Q. Kemp-type elimination of 1-arylsulfonyl-3-iodo-1 H-indazoles. SYNTHETIC COMMUN 2021. [DOI: 10.1080/00397911.2021.1939053] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/21/2022]
Affiliation(s)
- Jiawei Zhang
- Department of Medicinal Chemistry, School of Pharmacy, Fudan University, Shanghai, China
| | - Kaimin Kong
- Department of Medicinal Chemistry, School of Pharmacy, Fudan University, Shanghai, China
| | - Xiangjun Li
- Department of Medicinal Chemistry, School of Pharmacy, Fudan University, Shanghai, China
| | - Qian Zhang
- Department of Medicinal Chemistry, School of Pharmacy, Fudan University, Shanghai, China
| |
Collapse
|
6
|
Gerlt JA. Evolution of Enzyme Function and the Development of Catalytic Efficiency: Triosephosphate Isomerase, Jeremy R. Knowles, and W. John Albery. Biochemistry 2021; 60:3529-3538. [PMID: 34015914 DOI: 10.1021/acs.biochem.1c00211] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Every reader knows that an enzyme accelerates a reaction by reducing the activation-energy barrier. However, understanding how this is achieved by the structure of the enzyme and its interactions with stable complexes and transition states and, then, using this to (re)design enzymes to catalyze novel reactions remain the "holy grail" of mechanistic enzymology. The necessary foundation is the free-energy profile that specifies the energies of the bound substate, product, and intervening intermediates as well as the transition states by which they are interconverted. When this free-energy profile is compared to that for the uncatalyzed reaction, strategies for establishing and enhancing catalysis can be identified. This Perspective reminds readers that the first free-energy profile determined for an enzyme-catalyzed reaction, that for triosephosphate isomerase, was published in Biochemistry in 1976 by Jeremy R. Knowles, W. John Albery, and co-workers. They used the profile to propose three steps of increasing "subtlety" that can be influenced by evolutionary pressure to increase the flux through the reaction coordinate: (1) "uniform binding" of the substrate, product, and intermediates; (2) "differential binding" of complexes so that these are isoenergetic (to minimize the energy of the intervening transition states); and (3) "catalysis of an elementary step" in which the transition state for the kinetically significant chemical step is stabilized so that flux can be determined by the rate of substrate binding or product dissociation. These papers continue to guide mechanistic studies of enzyme-catalyzed reactions and provide principles for the (re)design of novel enzymes.
Collapse
Affiliation(s)
- John A Gerlt
- Departments of Biochemistry and Chemistry and Institute for Genomic Biology, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, United States
| |
Collapse
|
7
|
Bunzel HA, Anderson JLR, Mulholland AJ. Designing better enzymes: Insights from directed evolution. Curr Opin Struct Biol 2021; 67:212-218. [PMID: 33517098 DOI: 10.1016/j.sbi.2020.12.015] [Citation(s) in RCA: 33] [Impact Index Per Article: 11.0] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/12/2020] [Revised: 12/03/2020] [Accepted: 12/28/2020] [Indexed: 12/18/2022]
Abstract
De novo enzymes can be created by computational design and directed evolution. Here, we review recent insights into the origins of catalytic power in evolved designer enzymes to pinpoint opportunities for next-generation designs: Evolution precisely organizes active sites, introduces catalytic H-bonding networks, invokes electrostatic catalysis, and creates dynamical networks embedding the active site in a reactive protein scaffold. Such insights foster our fundamental knowledge of enzyme catalysis and fuel the future design of tailor-made enzymes.
Collapse
Affiliation(s)
- H Adrian Bunzel
- School of Biochemistry, University of Bristol, Bristol BS8 1TD, UK; Centre for Computational Chemistry, School of Chemistry, University of Bristol, Bristol BS8 1TS, UK
| | | | - Adrian J Mulholland
- Centre for Computational Chemistry, School of Chemistry, University of Bristol, Bristol BS8 1TS, UK.
| |
Collapse
|
8
|
Li WL, Head-Gordon T. Catalytic Principles from Natural Enzymes and Translational Design Strategies for Synthetic Catalysts. ACS CENTRAL SCIENCE 2021; 7:72-80. [PMID: 33532570 PMCID: PMC7844850 DOI: 10.1021/acscentsci.0c01556] [Citation(s) in RCA: 30] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/18/2020] [Indexed: 05/19/2023]
Abstract
As biocatalysts, enzymes are characterized by their high catalytic efficiency and strong specificity but are relatively fragile by requiring narrow and specific reactive conditions for activity. Synthetic catalysts offer an opportunity for more chemical versatility operating over a wider range of conditions but currently do not reach the remarkable performance of natural enzymes. Here we consider some new design strategies based on the contributions of nonlocal electric fields and thermodynamic fluctuations to both improve the catalytic step and turnover for rate acceleration in arbitrary synthetic catalysts through bioinspired studies of natural enzymes. With a focus on the enzyme as a whole catalytic construct, we illustrate the translational impact of natural enzyme principles to synthetic enzymes, supramolecular capsules, and electrocatalytic surfaces.
Collapse
Affiliation(s)
- Wan-Lu Li
- Chemical
Sciences Division, Lawrence Berkeley National
Laboratory, Berkeley, California 94720, United States
- Kenneth
S. Pitzer Center for Theoretical Chemistry, University of California Berkeley, Berkeley, California 94720, United States
| | - Teresa Head-Gordon
- Chemical
Sciences Division, Lawrence Berkeley National
Laboratory, Berkeley, California 94720, United States
- Kenneth
S. Pitzer Center for Theoretical Chemistry, University of California Berkeley, Berkeley, California 94720, United States
- Department of Chemistry, Department of Chemical and Biomolecular Engineering, and Department of
Bioengineering, University of California
Berkeley, Berkeley, California 94720, United States
| |
Collapse
|
9
|
Biler M, Crean RM, Schweiger AK, Kourist R, Kamerlin SCL. Ground-State Destabilization by Active-Site Hydrophobicity Controls the Selectivity of a Cofactor-Free Decarboxylase. J Am Chem Soc 2020; 142:20216-20231. [PMID: 33180505 PMCID: PMC7735706 DOI: 10.1021/jacs.0c10701] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/08/2020] [Indexed: 01/11/2023]
Abstract
Bacterial arylmalonate decarboxylase (AMDase) and evolved variants have become a valuable tool with which to access both enantiomers of a broad range of chiral arylaliphatic acids with high optical purity. Yet, the molecular principles responsible for the substrate scope, activity, and selectivity of this enzyme are only poorly understood to date, greatly hampering the predictability and design of improved enzyme variants for specific applications. In this work, empirical valence bond and metadynamics simulations were performed on wild-type AMDase and variants thereof to obtain a better understanding of the underlying molecular processes determining reaction outcome. Our results clearly reproduce the experimentally observed substrate scope and support a mechanism driven by ground-state destabilization of the carboxylate group being cleaved by the enzyme. In addition, our results indicate that, in the case of the nonconverted or poorly converted substrates studied in this work, increased solvent exposure of the active site upon binding of these substrates can disturb the vulnerable network of interactions responsible for facilitating the AMDase-catalyzed cleavage of CO2. Finally, our results indicate a switch from preferential cleavage of the pro-(R) to the pro-(S) carboxylate group in the CLG-IPL variant of AMDase for all substrates studied. This appears to be due to the emergence of a new hydrophobic pocket generated by the insertion of the six amino acid substitutions, into which the pro-(S) carboxylate binds. Our results allow insight into the tight interaction network determining AMDase selectivity, which in turn provides guidance for the identification of target residues for future enzyme engineering.
Collapse
Affiliation(s)
- Michal Biler
- Department
of Chemistry−BMC, Uppsala University, BMC Box 576, S-751 23 Uppsala, Sweden
| | - Rory M. Crean
- Department
of Chemistry−BMC, Uppsala University, BMC Box 576, S-751 23 Uppsala, Sweden
| | - Anna K. Schweiger
- Institute
of Molecular Biotechnology, Graz University
of Technology, NAWI Graz,
Petersgasse 14, 8010 Graz, Austria
| | - Robert Kourist
- Institute
of Molecular Biotechnology, Graz University
of Technology, NAWI Graz,
Petersgasse 14, 8010 Graz, Austria
| | | |
Collapse
|
10
|
Xia YL, Li YP, Fu YX, Liu SQ. The Energetic Origin of Different Catalytic Activities in Temperature-Adapted Trypsins. ACS OMEGA 2020; 5:25077-25086. [PMID: 33043186 PMCID: PMC7542600 DOI: 10.1021/acsomega.0c02401] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/22/2020] [Accepted: 09/09/2020] [Indexed: 05/08/2023]
Abstract
Psychrophilic enzymes were always observed to have higher catalytic activity (k cat) than their mesophilic homologs at room temperature, while the origin of this phenomenon remains obscure. Here, we used two different temperature-adapted trypsins, the psychrophilic Atlantic cod trypsin (ACT) and the mesophilic bovine trypsin (BT), as a model system to explore the energetic origin of their different catalytic activities using computational methods. The results reproduce the characteristic changing trends in the activation free energy, activation enthalpy, and activation entropy between the psychrophilic and mesophilic enzymes, where, in particular, the slightly decreased activation free energy of ACT is determined by its considerably reduced activation enthalpy rather than by its more negative activation entropy compared to BT. The calculated electrostatic contributions to the solvation free energies in the reactant state/ground sate (RS/GS) and transition state (TS) show that, going from BT to ACT, the TS stabilization has a predominant effect over the RS stabilization on lowering the activation enthalpy of ACT. Comparison between the solvation energy components reveals a more optimized electrostatic preorganization to the TS in ACT, which provides a larger stabilization to the TS through reducing the reorganization energy, thus resulting in the lower activation enthalpy and hence lower activation free energy of ACT. Thus, it can be concluded that it is the difference in the protein electrostatic environment, and hence its different stabilizing effects on the TS, that brings about the different catalytic activities of different temperature-adapted trypsins.
Collapse
Affiliation(s)
- Yuan-Ling Xia
- State
Key Laboratory for Conservation and Utilization of Bio-Resources in
Yunnan, Yunnan University, Kunming 650091, Yunnan, China
- Editorial
Office of Journal of Yunnan University (Natural Sciences Edition), Yunnan University, Kunming 650091, Yunnan, China
| | - Yong-Ping Li
- School
of Agriculture, Yunnan University, Kunming 650091, Yunnan, China
| | - Yun-Xin Fu
- State
Key Laboratory for Conservation and Utilization of Bio-Resources in
Yunnan, Yunnan University, Kunming 650091, Yunnan, China
- Human
Genetics Center and Division of Biostatistics, School of Public Health, The University of Texas Health Science Center, Houston, Texas 77030, United States
| | - Shu-Qun Liu
- State
Key Laboratory for Conservation and Utilization of Bio-Resources in
Yunnan, Yunnan University, Kunming 650091, Yunnan, China
| |
Collapse
|
11
|
Crean RM, Gardner JM, Kamerlin SCL. Harnessing Conformational Plasticity to Generate Designer Enzymes. J Am Chem Soc 2020; 142:11324-11342. [PMID: 32496764 PMCID: PMC7467679 DOI: 10.1021/jacs.0c04924] [Citation(s) in RCA: 49] [Impact Index Per Article: 12.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/05/2020] [Indexed: 02/08/2023]
Abstract
Recent years have witnessed an explosion of interest in understanding the role of conformational dynamics both in the evolution of new enzymatic activities from existing enzymes and in facilitating the emergence of enzymatic activity de novo on scaffolds that were previously non-catalytic. There are also an increasing number of examples in the literature of targeted engineering of conformational dynamics being successfully used to alter enzyme selectivity and activity. Despite the obvious importance of conformational dynamics to both enzyme function and evolvability, many (although not all) computational design approaches still focus either on pure sequence-based approaches or on using structures with limited flexibility to guide the design. However, there exist a wide variety of computational approaches that can be (re)purposed to introduce conformational dynamics as a key consideration in the design process. Coupled with laboratory evolution and more conventional existing sequence- and structure-based approaches, these techniques provide powerful tools for greatly expanding the protein engineering toolkit. This Perspective provides an overview of evolutionary studies that have dissected the role of conformational dynamics in facilitating the emergence of novel enzymes, as well as advances in computational approaches that allow one to target conformational dynamics as part of enzyme design. Harnessing conformational dynamics in engineering studies is a powerful paradigm with which to engineer the next generation of designer biocatalysts.
Collapse
Affiliation(s)
- Rory M. Crean
- Department of Chemistry -
BMC, Uppsala University, Box 576, 751 23 Uppsala, Sweden
| | - Jasmine M. Gardner
- Department of Chemistry -
BMC, Uppsala University, Box 576, 751 23 Uppsala, Sweden
| | - Shina C. L. Kamerlin
- Department of Chemistry -
BMC, Uppsala University, Box 576, 751 23 Uppsala, Sweden
| |
Collapse
|
12
|
Risso VA, Romero-Rivera A, Gutierrez-Rus LI, Ortega-Muñoz M, Santoyo-Gonzalez F, Gavira JA, Sanchez-Ruiz JM, Kamerlin SCL. Enhancing a de novo enzyme activity by computationally-focused ultra-low-throughput screening. Chem Sci 2020; 11:6134-6148. [PMID: 32832059 PMCID: PMC7407621 DOI: 10.1039/d0sc01935f] [Citation(s) in RCA: 17] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/05/2020] [Accepted: 05/18/2020] [Indexed: 01/02/2023] Open
Abstract
Directed evolution has revolutionized protein engineering. Still, enzyme optimization by random library screening remains sluggish, in large part due to futile probing of mutations that are catalytically neutral and/or impair stability and folding. FuncLib is a novel approach which uses phylogenetic analysis and Rosetta design to rank enzyme variants with multiple mutations, on the basis of predicted stability. Here, we use it to target the active site region of a minimalist-designed, de novo Kemp eliminase. The similarity between the Michaelis complex and transition state for the enzymatic reaction makes this system particularly challenging to optimize. Yet, experimental screening of a small number of active-site variants at the top of the predicted stability ranking leads to catalytic efficiencies and turnover numbers (∼2 × 104 M-1 s-1 and ∼102 s-1) for this anthropogenic reaction that compare favorably to those of modern natural enzymes. This result illustrates the promise of FuncLib as a powerful tool with which to speed up directed evolution, even on scaffolds that were not originally evolved for those functions, by guiding screening to regions of the sequence space that encode stable and catalytically diverse enzymes. Empirical valence bond calculations reproduce the experimental activation energies for the optimized eliminases to within ∼2 kcal mol-1 and indicate that the enhanced activity is linked to better geometric preorganization of the active site. This raises the possibility of further enhancing the stability-guidance of FuncLib by computational predictions of catalytic activity, as a generalized approach for computational enzyme design.
Collapse
Affiliation(s)
- Valeria A Risso
- Departamento de Química Física, Facultad de Ciencias , Unidad de Excelencia de Química aplicada a Biomedicina y Medioambiente (UEQ) , Universidad de Granada , 18071 Granada , Spain .
| | - Adrian Romero-Rivera
- Science for Life Laboratory , Department of Chemistry-BMC , Uppsala University , BMC Box 576 , S-751 23 Uppsala , Sweden .
| | - Luis I Gutierrez-Rus
- Departamento de Química Física, Facultad de Ciencias , Unidad de Excelencia de Química aplicada a Biomedicina y Medioambiente (UEQ) , Universidad de Granada , 18071 Granada , Spain .
| | - Mariano Ortega-Muñoz
- Departamento de Química Orgánica , Facultad de Ciencias , Unidad de Excelencia de Química aplicada a Biomedicina y Medioambiente (UEQ) , Universidad de Granada , 18071 Granada , Spain
| | - Francisco Santoyo-Gonzalez
- Departamento de Química Orgánica , Facultad de Ciencias , Unidad de Excelencia de Química aplicada a Biomedicina y Medioambiente (UEQ) , Universidad de Granada , 18071 Granada , Spain
| | - Jose A Gavira
- Laboratorio de Estudios Cristalográficos , Instituto Andaluz de Ciencias de la Tierra , CSIC, Unidad de Excelencia de Química aplicada a Biomedicina y Medioambiente (UEQ) , University of Granada , Avenida de las Palmeras 4 , 18100 Armilla , Granada , Spain
| | - Jose M Sanchez-Ruiz
- Departamento de Química Física, Facultad de Ciencias , Unidad de Excelencia de Química aplicada a Biomedicina y Medioambiente (UEQ) , Universidad de Granada , 18071 Granada , Spain .
| | - Shina C L Kamerlin
- Science for Life Laboratory , Department of Chemistry-BMC , Uppsala University , BMC Box 576 , S-751 23 Uppsala , Sweden .
| |
Collapse
|
13
|
Mondal D, Kolev V, Warshel A. Combinatorial Approach for Exploring Conformational Space and Activation Barriers in Computer-Aided Enzyme Design. ACS Catal 2020; 10:6002-6012. [PMID: 34178420 PMCID: PMC8225234 DOI: 10.1021/acscatal.0c01206] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Computer-aided enzyme design is a field of great potential importance for biotechnological applications, medical advances, and a fundamental understanding of enzyme action. However, reaching a predictive ability in this direction is extremely challenging. It requires both the ability to predict quantitatively the activation barriers in cases where the structure and sequence are known and the ability to predict the effect of different mutations. In this work, we propose a protocol for predicting reasonable starting structures of mutants of proteins with known structures and for calculating the activation barriers of the generated mutants. Our approach also allows us to use the predicted structures of the generated mutant to predict structures and activation barriers for subsequent set of mutations. This protocol is used to examine the reliability of the in silico directed evolution of Kemp eliminase and haloalkane dehalogenase. We also used the results of single and double mutations as a base for predicting the effect of transition-state stabilization by multiple concurrent mutations. This strategy seems to be useful in creating an activity funnel that provides a qualitative ranking of the catalytic power of different mutants.
Collapse
Affiliation(s)
- Dibyendu Mondal
- Department of Chemistry, University of Southern California, Los Angeles, California 90089, United States
| | - Vesselin Kolev
- Department of Chemistry, University of Southern California, Los Angeles, California 90089, United States
| | - Arieh Warshel
- Department of Chemistry, University of Southern California, Los Angeles, California 90089, United States
| |
Collapse
|
14
|
Beker W, Sokalski WA. Bottom-Up Nonempirical Approach To Reducing Search Space in Enzyme Design Guided by Catalytic Fields. J Chem Theory Comput 2020; 16:3420-3429. [PMID: 32282205 PMCID: PMC7467639 DOI: 10.1021/acs.jctc.0c00139] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Abstract
Currently developed protocols of theozyme design still lead to biocatalysts with much lower catalytic activity than enzymes existing in nature, and, so far, the only avenue of improvement was the in vitro laboratory-directed evolution (LDE) experiments. In this paper, we propose a different strategy based on "reversed" methodology of mutation prediction. Instead of common "top-down" approach, requiring numerous assumptions and vast computational effort, we argue for a "bottom-up" approach that is based on the catalytic fields derived directly from transition state and reactant complex wave functions. This enables direct one-step determination of the general quantitative angular characteristics of optimal catalytic site and simultaneously encompasses both the transition-state stabilization (TSS) and ground-state destabilization (GSD) effects. We further extend the static catalytic field approach by introducing a library of atomic multipoles for amino acid side-chain rotamers, which, together with the catalytic field, allow one to determine the optimal side-chain orientations of charged amino acids constituting the elusive structure of a preorganized catalytic environment. Obtained qualitative agreement with experimental LDE data for Kemp eliminase KE07 mutants validates the proposed procedure, yielding, in addition, a detailed insight into possible dynamic and epistatic effects.
Collapse
Affiliation(s)
- Wiktor Beker
- Wroclaw University of Science and Technology, Wroclaw, Poland
| | | |
Collapse
|
15
|
Kries H, Bloch JS, Bunzel HA, Pinkas DM, Hilvert D. Contribution of Oxyanion Stabilization to Kemp Eliminase Efficiency. ACS Catal 2020. [DOI: 10.1021/acscatal.0c00575] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/02/2023]
Affiliation(s)
- Hajo Kries
- Laboratory of Organic Chemistry, ETH Zurich, 8093 Zurich, Switzerland
| | - Joël S. Bloch
- Laboratory of Organic Chemistry, ETH Zurich, 8093 Zurich, Switzerland
| | - H. Adrian Bunzel
- Laboratory of Organic Chemistry, ETH Zurich, 8093 Zurich, Switzerland
| | - Daniel M. Pinkas
- Laboratory of Organic Chemistry, ETH Zurich, 8093 Zurich, Switzerland
| | - Donald Hilvert
- Laboratory of Organic Chemistry, ETH Zurich, 8093 Zurich, Switzerland
| |
Collapse
|
16
|
Prejanò M, Medina FE, Ramos MJ, Russo N, Fernandes PA, Marino T. How the Destabilization of a Reaction Intermediate Affects Enzymatic Efficiency: The Case of Human Transketolase. ACS Catal 2020; 10:2872-2881. [PMID: 33828899 PMCID: PMC8016368 DOI: 10.1021/acscatal.9b04690] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/30/2019] [Revised: 02/04/2020] [Indexed: 12/16/2022]
Abstract
![]()
Atomic
resolution X-ray crystallography has shown that an intermediate
(the X5P-ThDP adduct) of the catalytic cycle of transketolase (TK)
displays a significant, putatively highly energetic, out-of-plane
distortion in a sp2 carbon
adjacent to a lytic bond, suggested to lower the barrier of the subsequent
step, and thus was postulated to embody a clear-cut demonstration
of the intermediate destabilization effect. The lytic
bond of the subsequent rate-limiting step was very elongated in the
X-ray structure (1.61 Å), which was proposed to be a consequence
of the out-of-plane distortion. Here we use high-level QM and QM/MM
calculations to study the intermediate destabilization effect. We show that the intrinsic energy penalty for the observed
distortion is small (0.2 kcal·mol–1) and that
the establishment of a favorable hydrogen bond within X5P-ThDP, instead
of enzyme steric strain, was found to be the main cause for the distortion.
As the net energetic effect of the distortion is small, the establishment
of the internal hydrogen bond (−0.6 kcal·mol–1) offsets the associated penalty. This makes the distorted structure
more stable than the nondistorted one. Even though the energy contributions
determined here are close to the accuracy of the computational methods
in estimating penalties for geometric distortions, our data show that
the intermediate destabilization effect provides
a small contribution to the observed reaction rate and does not represent
a catalytic effect that justifies the many orders of magnitude which
enzymes accelerate reaction rates. The results help to understand
the intrinsic enzymatic machinery behind enzyme’s amazing proficiency.
Collapse
Affiliation(s)
- Mario Prejanò
- Dipartimento di Chimica e Tecnologie Chimiche, Università della Calabria, 87036 Arcavacata di Rende (CS), Italy
| | - Fabiola E. Medina
- 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
| | - Nino Russo
- Dipartimento di Chimica e Tecnologie Chimiche, Università della Calabria, 87036 Arcavacata di Rende (CS), Italy
| | - 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
| | - Tiziana Marino
- Dipartimento di Chimica e Tecnologie Chimiche, Università della Calabria, 87036 Arcavacata di Rende (CS), Italy
| |
Collapse
|
17
|
Berta D, Buigues PJ, Badaoui M, Rosta E. Cations in motion: QM/MM studies of the dynamic and electrostatic roles of H + and Mg 2+ ions in enzyme reactions. Curr Opin Struct Biol 2020; 61:198-206. [PMID: 32065923 DOI: 10.1016/j.sbi.2020.01.002] [Citation(s) in RCA: 19] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/07/2019] [Revised: 01/04/2020] [Accepted: 01/06/2020] [Indexed: 12/17/2022]
Abstract
Here we discuss current trends in the simulations of enzymatic reactions focusing on phosphate catalysis. The mechanistic details of the proton transfers coupled to the phosphate cleavage is one of the key challenges in QM/MM calculations of these and other enzyme catalyzed reactions. The lack of experimental information offers both an opportunity for computations as well as often unresolved controversies. We discuss the example of small GTPases including the important human Ras protein. The high dimensionality and chemical complexity of these reactions demand carefully chosen computational techniques both in terms of the underlying quantum chemical theory and the sampling of the conformational ensemble. We also point out the important role of Mg2+ ions, and recent advances in their transient involvement in the catalytic mechanisms.
Collapse
Affiliation(s)
- Dénes Berta
- Department of Chemistry, King's College London, London, SE1 1DB, United Kingdom
| | - Pedro J Buigues
- Department of Chemistry, King's College London, London, SE1 1DB, United Kingdom
| | - Magd Badaoui
- Department of Chemistry, King's College London, London, SE1 1DB, United Kingdom
| | - Edina Rosta
- Department of Chemistry, King's College London, London, SE1 1DB, United Kingdom.
| |
Collapse
|
18
|
Acosta-Silva C, Bertran J, Branchadell V, Oliva A. Kemp Elimination Reaction Catalyzed by Electric Fields. Chemphyschem 2020; 21:295-306. [PMID: 31840917 DOI: 10.1002/cphc.201901155] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/05/2019] [Revised: 12/12/2019] [Indexed: 12/20/2022]
Abstract
The Kemp elimination reaction is the most widely used in the de novo design of new enzymes. The effect of two different kinds of electric fields in the reactions of acetate as a base with benzisoxazole and 5-nitrobenzisoxazole as substrates have been theoretically studied. The effect of the solvent reaction field has been calculated using the SMD continuum model for several solvents; we have shown that solvents inhibit both reactions, the decrease of the reaction rate being larger as far as the dielectric constant is increased. The diminution of the reaction rate is especially remarkable between aprotic organic solvents and protic solvents as water, the electrostatic term of the hydrogen bonds being the main factor for the large inhibitory effect of water. The presence of an external electric field oriented in the direction of the charge transfer (z axis) increases it and, so, the reaction rate. In the reaction of the nitro compound, if the electric field is oriented in an orthogonal direction (x axis) the charge transfer to the NO2 group is favored and there is a subsequent increase of the reaction rate. However, this increase is smaller than the one produced by the field in the z axis. It is worthwhile mentioning that one of the main effects of external electric fields of intermediate intensity is the reorientation of the reactants. Finally, the implications of our results in the de novo design of enzymes are discussed.
Collapse
Affiliation(s)
- Carles Acosta-Silva
- Departament de Química, Universitat Autònoma de Barcelona, 08193, Bellaterra, Spain
| | - Joan Bertran
- Departament de Química, Universitat Autònoma de Barcelona, 08193, Bellaterra, Spain
| | - Vicenç Branchadell
- Departament de Química, Universitat Autònoma de Barcelona, 08193, Bellaterra, Spain
| | - Antoni Oliva
- Departament de Química, Universitat Autònoma de Barcelona, 08193, Bellaterra, Spain
| |
Collapse
|
19
|
Kulkarni Y, Kamerlin SCL. Computational physical organic chemistry using the empirical valence bond approach. ADVANCES IN PHYSICAL ORGANIC CHEMISTRY 2019. [DOI: 10.1016/bs.apoc.2019.07.001] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/01/2022]
|
20
|
Exploring the challenges of computational enzyme design by rebuilding the active site of a dehalogenase. Proc Natl Acad Sci U S A 2018; 116:389-394. [PMID: 30587585 DOI: 10.1073/pnas.1804979115] [Citation(s) in RCA: 22] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
Rational enzyme design presents a major challenge that has not been overcome by computational approaches. One of the key challenges is the difficulty in assessing the magnitude of the maximum possible catalytic activity. In an attempt to overcome this challenge, we introduce a strategy that takes an active enzyme (assuming that its activity is close to the maximum possible activity), design mutations that reduce the catalytic activity, and then try to restore that catalysis by mutating other residues. Here we take as a test case the enzyme haloalkane dehalogenase (DhlA), with a 1,2-dichloroethane substrate. We start by demonstrating our ability to reproduce the results of single mutations. Next, we design mutations that reduce the enzyme activity and finally design double mutations that are aimed at restoring the activity. Using the computational predictions as a guide, we conduct an experimental study that confirms our prediction in one case and leads to inconclusive results in another case with 1,2-dichloroethane as substrate. Interestingly, one of our predicted double mutants catalyzes dehalogenation of 1,2-dibromoethane more efficiently than the wild-type enzyme.
Collapse
|
21
|
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
|
22
|
Barbosa ACC, Neves RPP, Sousa SF, Ramos MJ, Fernandes PA. Mechanistic Studies of a Flavin Monooxygenase: Sulfur Oxidation of Dibenzothiophenes by DszC. ACS Catal 2018. [DOI: 10.1021/acscatal.8b01877] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Ana C. C. Barbosa
- 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
| | - Rui P. P. Neves
- 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
| | - Sérgio F. Sousa
- 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
| | - 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
| |
Collapse
|
23
|
|
24
|
Vaissier V, Sharma SC, Schaettle K, Zhang T, Head-Gordon T. Computational Optimization of Electric Fields for Improving Catalysis of a Designed Kemp Eliminase. ACS Catal 2017. [DOI: 10.1021/acscatal.7b03151] [Citation(s) in RCA: 47] [Impact Index Per Article: 6.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Valerie Vaissier
- Chemical
Sciences Division, Lawrence Berkeley National Laboratories, Berkeley, California 94720, United States
| | | | | | | | - Teresa Head-Gordon
- Chemical
Sciences Division, Lawrence Berkeley National Laboratories, Berkeley, California 94720, United States
| |
Collapse
|
25
|
Abstract
Hepatitis C virus (HCV) currently affects several million people across the globe. One of the major classes of drugs against HCV inhibits the NS3/4A protease of the polyprotein chain. Efficacy of these drugs is severely limited due to the high mutation rate that results in several genetically related quasispecies. The molecular mechanism of drug resistance is frequently deduced from structural studies and binding free energies. However, prediction of new mutations requires the evaluation of both binding free energy of the drug as well as the parameters (kcat and KM) for the natural substrate. The vitality values offer a good approach to investigate and predict mutations that render resistance to the inhibitor. A successful mutation should only affect the binding of the drug and not the catalytic activity and binding of the natural substrate. In this article, we have calculated the vitality values for four known drug inhibitors that are either currently in use or in clinical trials, evaluating binding free energies by the relevant PDLD/S-LRA method and activation barriers by the EVB method. The molecular details pertaining to resistance are also discussed. We show that our calculations are able to reproduce the catalytic effects and binding free energies in a good agreement with the corresponding observed values. Importantly, previous computational approaches have not been able to achieve this task. The trend for the vitality values is in accordance with experimental findings. Finally, we calculate the vitality values for mutations that have either not been studied experimentally or reported for some inhibitors.
Collapse
Affiliation(s)
- Garima Jindal
- Department of Chemistry, University of Southern California , 3620 McClintock Avenue, Los Angeles, California 90089, United States
| | - Dibyendu Mondal
- Department of Chemistry, University of Southern California , 3620 McClintock Avenue, Los Angeles, California 90089, United States
| | - Arieh Warshel
- Department of Chemistry, University of Southern California , 3620 McClintock Avenue, Los Angeles, California 90089, United States
| |
Collapse
|