1
|
Shandilya E, Rallabandi B, Maiti S. In situ enzymatic control of colloidal phoresis and catalysis through hydrolysis of ATP. Nat Commun 2024; 15:3603. [PMID: 38684662 PMCID: PMC11059368 DOI: 10.1038/s41467-024-47912-2] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/21/2023] [Accepted: 04/16/2024] [Indexed: 05/02/2024] Open
Abstract
The ability to sense chemical gradients and respond with directional motility and chemical activity is a defining feature of complex living systems. There is a strong interest among scientists to design synthetic systems that emulate these properties. Here, we realize and control such behaviors in a synthetic system by tailoring multivalent interactions of adenosine nucleotides with catalytic microbeads. We first show that multivalent interactions of the bead with gradients of adenosine mono-, di- and trinucleotides (AM/D/TP) control both the phoretic motion and a proton-transfer catalytic reaction, and find that both effects are diminished greatly with increasing valence of phosphates. We exploit this behavior by using enzymatic hydrolysis of ATP to AMP, which downregulates multivalent interactivity in situ. This produces a sudden increase in transport of the catalytic microbeads (a phoretic jump), which is accompanied by increased catalytic activity. Finally, we show how this enzymatic activity can be systematically tuned, leading to simultaneous in situ spatial and temporal control of the location of the microbeads, as well as the products of the reaction that they catalyze. These findings open up new avenues for utilizing multivalent interaction-mediated programming of complex chemo-mechanical behaviors into active systems.
Collapse
Affiliation(s)
- Ekta Shandilya
- Department of Chemical Sciences, Indian Institute of Science Education and Research (IISER), Mohali, Knowledge City, Manauli, 140306, India
| | - Bhargav Rallabandi
- Department of Mechanical Engineering, University of California, Riverside, CA, 92521, USA.
| | - Subhabrata Maiti
- Department of Chemical Sciences, Indian Institute of Science Education and Research (IISER), Mohali, Knowledge City, Manauli, 140306, India.
| |
Collapse
|
2
|
Jiang Y, Ding N, Shao Q, Stull SL, Cheng Z, Yang ZJ. Substrate Positioning Dynamics Involves a Non-Electrostatic Component to Mediate Catalysis. J Phys Chem Lett 2023; 14:11480-11489. [PMID: 38085952 DOI: 10.1021/acs.jpclett.3c02444] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2023]
Abstract
Substrate positioning dynamics (SPD) orients the substrate in the active site, thereby influencing catalytic efficiency. However, it remains unknown whether SPD effects originate primarily from electrostatic perturbation inside the enzyme or can independently mediate catalysis with a significant non-electrostatic component. In this work, we investigated how the non-electrostatic component of SPD affects transition state (TS) stabilization. Using high-throughput enzyme modeling, we selected Kemp eliminase variants with similar electrostatics inside the enzyme but significantly different SPD. The kinetic parameters of these mutants were experimentally characterized. We observed a valley-shaped, two-segment linear correlation between the TS stabilization free energy (converted from kinetic parameters) and substrate positioning index (a metric to quantify SPD). The energy varies by approximately 2 kcal/mol. Favorable SPD was observed for the distal mutant R154W, increasing the proportion of reactive conformations and leading to the lowest activation free energy. These results indicate the substantial contribution of the non-electrostatic component of SPD to enzyme catalytic efficiency.
Collapse
Affiliation(s)
- Yaoyukun Jiang
- Department of Chemistry, Vanderbilt University, Nashville, Tennessee 37235, United States
| | - Ning Ding
- Department of Chemistry, Vanderbilt University, Nashville, Tennessee 37235, United States
| | - Qianzhen Shao
- Department of Chemistry, Vanderbilt University, Nashville, Tennessee 37235, United States
| | - Sebastian L Stull
- Department of Chemistry, Vanderbilt University, Nashville, Tennessee 37235, United States
| | - Zihao Cheng
- Department of Chemistry, Vanderbilt University, Nashville, Tennessee 37235, United States
| | - Zhongyue J Yang
- Department of Chemistry, Vanderbilt University, Nashville, Tennessee 37235, United States
- Center for Structural Biology, Vanderbilt University, Nashville, Tennessee 37235, United States
- Vanderbilt Institute of Chemical Biology, Vanderbilt University, Nashville, Tennessee 37235, United States
- Data Science Institute, Vanderbilt University, Nashville, Tennessee 37235, United States
- Department of Chemical and Biomolecular Engineering, Vanderbilt University, Nashville, Tennessee 37235, United States
| |
Collapse
|
3
|
Jiang Y, Stull SL, Shao Q, Yang ZJ. Convergence in determining enzyme functional descriptors across Kemp eliminase variants. ELECTRONIC STRUCTURE (BRISTOL, ENGLAND) 2022; 4:044007. [PMID: 37425623 PMCID: PMC10327861 DOI: 10.1088/2516-1075/acad51] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 07/11/2023]
Abstract
Molecular simulations have been extensively employed to accelerate biocatalytic discoveries. Enzyme functional descriptors derived from molecular simulations have been leveraged to guide the search for beneficial enzyme mutants. However, the ideal active-site region size for computing the descriptors over multiple enzyme variants remains untested. Here, we conducted convergence tests for dynamics-derived and electrostatic descriptors on 18 Kemp eliminase variants across six active-site regions with various boundary distances to the substrate. The tested descriptors include the root-mean-square deviation of the active-site region, the solvent accessible surface area ratio between the substrate and active site, and the projection of the electric field (EF) on the breaking C-H bond. All descriptors were evaluated using molecular mechanics methods. To understand the effects of electronic structure, the EF was also evaluated using quantum mechanics/molecular mechanics methods. The descriptor values were computed for 18 Kemp eliminase variants. Spearman correlation matrices were used to determine the region size condition under which further expansion of the region boundary does not substantially change the ranking of descriptor values. We observed that protein dynamics-derived descriptors, including RMSDactive_site and SASAratio, converge at a distance cutoff of 5 Å from the substrate. The electrostatic descriptor, EFC-H, converges at 6 Å using molecular mechanics methods with truncated enzyme models and 4 Å using quantum mechanics/molecular mechanics methods with whole enzyme model. This study serves as a future reference to determine descriptors for predictive modeling of enzyme engineering.
Collapse
Affiliation(s)
- Yaoyukun Jiang
- Department of Chemistry, Vanderbilt University, Nashville, TN 37235, United States of America
| | - Sebastian L Stull
- Department of Chemistry, Vanderbilt University, Nashville, TN 37235, United States of America
| | - Qianzhen Shao
- Department of Chemistry, Vanderbilt University, Nashville, TN 37235, United States of America
| | - Zhongyue J Yang
- Department of Chemistry, Vanderbilt University, Nashville, TN 37235, United States of America
- Center for Structural Biology, Vanderbilt University, Nashville, TN 37235, United States of America
- Vanderbilt Institute of Chemical Biology, Vanderbilt University, Nashville, TN 37235, United States of America
- Data Science Institute, Vanderbilt University, Nashville, TN 37235, United States of America
- Department of Chemical and Biomolecular Engineering, Vanderbilt University, Nashville, TN 37235, United States of America
| |
Collapse
|
4
|
Fan M, Cui L, He X, Zou X. Emerging Heterogeneous Supports for Efficient Electrocatalysis. SMALL METHODS 2022; 6:e2200855. [PMID: 36070422 DOI: 10.1002/smtd.202200855] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/03/2022] [Revised: 08/16/2022] [Indexed: 06/15/2023]
Abstract
Electrocatalysis plays a fundamental role in many fields, such as metallurgy, medicine, chemical industry, and energy conversion. Anchoring active electrocatalysts with controllable loading and uniform dispersion onto suitable supports has become an attractive topic. This is because the supports can not only have the potential to improve catalytic activity and stability through the interaction between support and catalytic center, but also can reduce precious metal consumption by improving atomic utilization. Herein, recent theoretical and experimental progresses concerning the development of supports to anchor electrocatalytic materials are first reviewed. Next, their controllable syntheses, characterization techniques, metal-support electronic interactions, and structure-performance relationships are presented. Some representative carbon supports and non-carbonaceous supports, as well as recently reported star supports such as 2D supports, single atom catalysts, and self-supported catalysts are also summarized. In addition, the significant role of support in stabilizing and regulating catalytic active sites is particularly emphasized. Finally, challenges, opportunities, key problems, and further promising solutions for supported catalysts are proposed.
Collapse
Affiliation(s)
- Meihong Fan
- School of Chemistry and Environmental Engineering, Changchun University of Science and Technology, Changchun, 130022, China
- School of Materials Science and Engineering, Changchun University of Science and Technology, Changchun, 130022, China
| | - Lili Cui
- School of Chemistry and Environmental Engineering, Changchun University of Science and Technology, Changchun, 130022, China
| | - Xingquan He
- School of Chemistry and Environmental Engineering, Changchun University of Science and Technology, Changchun, 130022, China
- School of Materials Science and Engineering, Changchun University of Science and Technology, Changchun, 130022, China
| | - Xiaoxin Zou
- State Key Laboratory of Inorganic Synthesis and Preparative Chemistry, College of Chemistry, Jilin University, Changchun, 130012, China
| |
Collapse
|
5
|
Zheng Y, Vaissier Welborn V. Tuning the Catalytic Activity of Synthetic Enzyme KE15 with DNA. J Phys Chem B 2022; 126:3407-3413. [PMID: 35483007 DOI: 10.1021/acs.jpcb.2c00765] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Efficiency improvement of synthetic enzymes through scaffold modifications suffers from limitations in terms of effectiveness, cost, and potential devastating consequences for protein structural stability. Here, we propose an alternative to scaffold modification, within electrostatic preorganization theory, where the enzyme's greater environment is designed to support the evolution of the reaction in the active site. We demonstrate the feasibility of such an approach by placing a (polar) DNA fragment in the surroundings of the Kemp eliminase enzyme KE15 (structure from Houk's group) and computing the resulting change in catalytic activity. We find that the introduction of a DNA fragment magnifies the contribution of protein residues to the stabilization of the transition state, estimated from electric field calculations with polarizable molecular dynamics. Our randomly generated test systems reveal a 2.0 kcal/mol reduction in activation energy, suggesting that even more significant catalytic improvements could be made by optimizing DNA size, sequence, and orientation with respect to the enzyme, validating our approach.
Collapse
Affiliation(s)
- Yi Zheng
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24060, United States
| | | |
Collapse
|
6
|
Abstract
QM/MM simulations have become an indispensable tool in many chemical and biochemical investigations. Considering the tremendous degree of success, including recognition by a 2013 Nobel Prize in Chemistry, are there still "burning challenges" in QM/MM methods, especially for biomolecular systems? In this short Perspective, we discuss several issues that we believe greatly impact the robustness and quantitative applicability of QM/MM simulations to many, if not all, biomolecules. We highlight these issues with observations and relevant advances from recent studies in our group and others in the field. Despite such limited scope, we hope the discussions are of general interest and will stimulate additional developments that help push the field forward in meaningful directions.
Collapse
Affiliation(s)
- Qiang Cui
- Departments of Chemistry, Physics, and Biomedical Engineering, Boston University, 590 Commonwealth Avenue, Boston, Massachusetts 02215, United States
| | - Tanmoy Pal
- Department of Chemistry, Boston University, 590 Commonwealth Avenue, Boston, Massachusetts 02215, United States
| | - Luke Xie
- Department of Chemistry, Boston University, 590 Commonwealth Avenue, Boston, Massachusetts 02215, United States
| |
Collapse
|
7
|
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
|
8
|
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
|
9
|
Dral PO, Wu X, Thiel W. Semiempirical Quantum-Chemical Methods with Orthogonalization and Dispersion Corrections. J Chem Theory Comput 2019; 15:1743-1760. [PMID: 30735388 PMCID: PMC6416713 DOI: 10.1021/acs.jctc.8b01265] [Citation(s) in RCA: 39] [Impact Index Per Article: 7.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/18/2018] [Indexed: 12/31/2022]
Abstract
We present two new semiempirical quantum-chemical methods with orthogonalization and dispersion corrections: ODM2 and ODM3 (ODM x). They employ the same electronic structure model as the OM2 and OM3 (OM x) methods, respectively. In addition, they include Grimme's dispersion correction D3 with Becke-Johnson damping and three-body corrections E ABC for Axilrod-Teller-Muto dispersion interactions as integral parts. Heats of formation are determined by adding explicitly computed zero-point vibrational energy and thermal corrections, in contrast to standard MNDO-type and OM x methods. We report ODM x parameters for hydrogen, carbon, nitrogen, oxygen, and fluorine that are optimized with regard to a wide range of carefully chosen state-of-the-art reference data. Extensive benchmarks show that the ODM x methods generally perform better than the available MNDO-type and OM x methods for ground-state and excited-state properties, while they describe noncovalent interactions with similar accuracy as OM x methods with a posteriori dispersion corrections.
Collapse
Affiliation(s)
- Pavlo O. Dral
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, 45470 Mülheim an der Ruhr, Germany
| | - Xin Wu
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, 45470 Mülheim an der Ruhr, Germany
| | - Walter Thiel
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, 45470 Mülheim an der Ruhr, Germany
| |
Collapse
|
10
|
Xue J, Huang X, Zhu Y. Using molecular dynamics simulations to evaluate active designs of cephradine hydrolase by molecular mechanics/Poisson–Boltzmann surface area and molecular mechanics/generalized Born surface area methods. RSC Adv 2019; 9:13868-13877. [PMID: 35519543 PMCID: PMC9064048 DOI: 10.1039/c9ra02406a] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/31/2019] [Accepted: 04/30/2019] [Indexed: 11/30/2022] Open
Abstract
The poor predictive accuracy of current computational enzyme design methods has led to low success rates of producing highly active variants that target non-natural substrates. In this report, a quantitative assessment approach based on molecular dynamics (MD) simulations was developed to eliminate false-positive enzyme designs at the computational stage. Taking cephradine hydrolase as an example, the apparent Michaelis binding constant (Km) and catalytic efficiency (kcat/Km) of designed variants were correlated with binding free energies and activation energy barriers, respectively, as calculated by molecular mechanics/Poisson–Boltzmann surface area (MM/PBSA) and molecular mechanics/generalized Born surface area (MM/GBSA) methods with explicit water considered based on general MD simulation protocols. The correlation results showed that both the MM/GBSA and MM/PBSA methods with a protein dielectric constant (εp = 4) could rank the variants well based on the predicted binding free energies between enzyme and the substrate. Furthermore, the activation energy barriers calculated by the MM/PBSA method with an εp = 24 correlated well with kcat/Km. Thus, false-positive variants obtained by the enzyme design program PRODA were eliminated prior to experimentation. Therefore, MD simulation-based quantitative assessment of designed variants greatly enhanced the predictive accuracy of computational enzyme design tools and should facilitate the construction of artificial enzymes with high catalytic activities toward non-natural substrates. A quantitative assessment method for computational enzyme design was developed to rank the active designs of cephradine hydrolase based on molecular dynamics simulation.![]()
Collapse
Affiliation(s)
- Jing Xue
- Department of Chemical Engineering
- Tsinghua University
- Beijing 100084
- China
| | - Xiaoqiang Huang
- Department of Chemical Engineering
- Tsinghua University
- Beijing 100084
- China
| | - Yushan Zhu
- Department of Chemical Engineering
- Tsinghua University
- Beijing 100084
- China
- MOE Key Lab for Industrial Biocatalysis
| |
Collapse
|
11
|
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
|
12
|
Hong NS, Petrović D, Lee R, Gryn'ova G, Purg M, Saunders J, Bauer P, Carr PD, Lin CY, Mabbitt PD, Zhang W, Altamore T, Easton C, Coote ML, Kamerlin SCL, Jackson CJ. The evolution of multiple active site configurations in a designed enzyme. Nat Commun 2018; 9:3900. [PMID: 30254369 PMCID: PMC6156567 DOI: 10.1038/s41467-018-06305-y] [Citation(s) in RCA: 56] [Impact Index Per Article: 9.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/04/2017] [Accepted: 08/21/2018] [Indexed: 12/11/2022] Open
Abstract
Developments in computational chemistry, bioinformatics, and laboratory evolution have facilitated the de novo design and catalytic optimization of enzymes. Besides creating useful catalysts, the generation and iterative improvement of designed enzymes can provide valuable insight into the interplay between the many phenomena that have been suggested to contribute to catalysis. In this work, we follow changes in conformational sampling, electrostatic preorganization, and quantum tunneling along the evolutionary trajectory of a designed Kemp eliminase. We observe that in the Kemp Eliminase KE07, instability of the designed active site leads to the emergence of two additional active site configurations. Evolutionary conformational selection then gradually stabilizes the most efficient configuration, leading to an improved enzyme. This work exemplifies the link between conformational plasticity and evolvability and demonstrates that residues remote from the active sites of enzymes play crucial roles in controlling and shaping the active site for efficient catalysis.
Collapse
Affiliation(s)
- Nan-Sook Hong
- Research School of Chemistry, Australian National University, Canberra, ACT, 2601, Australia
| | - Dušan Petrović
- Department of Chemistry, BMC, Uppsala University, Box 576, 751 23, Uppsala, Sweden
| | - Richmond Lee
- Research School of Chemistry, Australian National University, Canberra, ACT, 2601, Australia
| | - Ganna Gryn'ova
- Research School of Chemistry, Australian National University, Canberra, ACT, 2601, Australia.,Institut des Sciences et Ingénierie Chimiques, École Polytechnique Fédérale de Lausanne, 1015, Lausanne, Switzerland
| | - Miha Purg
- Department of Chemistry, BMC, Uppsala University, Box 576, 751 23, Uppsala, Sweden
| | - Jake Saunders
- Research School of Chemistry, Australian National University, Canberra, ACT, 2601, Australia
| | - Paul Bauer
- Department of Chemistry, BMC, Uppsala University, Box 576, 751 23, Uppsala, Sweden
| | - Paul D Carr
- Research School of Chemistry, Australian National University, Canberra, ACT, 2601, Australia
| | - Ching-Yeh Lin
- Research School of Chemistry, Australian National University, Canberra, ACT, 2601, Australia
| | - Peter D Mabbitt
- Research School of Chemistry, Australian National University, Canberra, ACT, 2601, Australia
| | - William Zhang
- Research School of Chemistry, Australian National University, Canberra, ACT, 2601, Australia
| | - Timothy Altamore
- Research School of Chemistry, Australian National University, Canberra, ACT, 2601, Australia
| | - Chris Easton
- Research School of Chemistry, Australian National University, Canberra, ACT, 2601, Australia
| | - Michelle L Coote
- Research School of Chemistry, Australian National University, Canberra, ACT, 2601, Australia
| | - Shina C L Kamerlin
- Department of Chemistry, BMC, Uppsala University, Box 576, 751 23, Uppsala, Sweden.
| | - Colin J Jackson
- Research School of Chemistry, Australian National University, Canberra, ACT, 2601, Australia.
| |
Collapse
|
13
|
Husch T, Reiher M. Comprehensive Analysis of the Neglect of Diatomic Differential Overlap Approximation. J Chem Theory Comput 2018; 14:5169-5179. [DOI: 10.1021/acs.jctc.8b00601] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022]
Affiliation(s)
- Tamara Husch
- Laboratorium für Physikalische Chemie, ETH Zürich, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| | - Markus Reiher
- Laboratorium für Physikalische Chemie, ETH Zürich, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| |
Collapse
|
14
|
Risso VA, Sanchez-Ruiz JM, Ozkan SB. Biotechnological and protein-engineering implications of ancestral protein resurrection. Curr Opin Struct Biol 2018; 51:106-115. [PMID: 29660672 DOI: 10.1016/j.sbi.2018.02.007] [Citation(s) in RCA: 26] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/21/2018] [Revised: 02/18/2018] [Accepted: 02/20/2018] [Indexed: 10/17/2022]
Abstract
Approximations to the sequences of ancestral proteins can be derived from the sequences of their modern descendants. Proteins encoded by such reconstructed sequences can be prepared in the laboratory and subjected to experimental scrutiny. These 'resurrected' ancestral proteins often display remarkable properties, reflecting ancestral adaptations to intra-cellular and extra-cellular environments that differed from the environments hosting modern/extant proteins. Recent experimental and computational work has specifically discussed high stability, substrate and catalytic promiscuity, conformational flexibility/diversity and altered patterns of interaction with other sub-cellular components. In this review, we discuss these remarkable properties as well as recent attempts to explore their biotechnological and protein-engineering potential.
Collapse
Affiliation(s)
- 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.
| | - S Banu Ozkan
- Department of Physics and Center for Biological Physics, Arizona State University, Tempe, AZ 85281, United States.
| |
Collapse
|
15
|
Chen X, Schwartz SD. Directed Evolution as a Probe of Rate Promoting Vibrations Introduced via Mutational Change. Biochemistry 2018; 57:3289-3298. [PMID: 29553716 DOI: 10.1021/acs.biochem.8b00185] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
In this article, we study with transition path sampling and reaction coordinate analysis how directed evolution in the Kemp eliminase family of artificial enzymes makes differential use of rapid rate promoting vibrations as a component of their chemical mechanism. Even though this family was initially created by placing the expected active site in a fixed protein matrix, we find a shift from largely static to more dynamic active sites that make use of donor-acceptor compression as the evolutionary process proceeds. We see that this introduction of dynamics significantly shifts the order of processes in the reaction. We also suggest that the lack of "design for dynamics" may help explain the relatively low proficiency of such designed enzymes.
Collapse
Affiliation(s)
- Xi Chen
- 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
|
16
|
Garcia-Borràs M, Houk KN, Jiménez-Osés G. Computational Design of Protein Function. COMPUTATIONAL TOOLS FOR CHEMICAL BIOLOGY 2017. [DOI: 10.1039/9781788010139-00087] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/28/2023]
Abstract
The computational design of enzymes is a tremendous challenge for both chemistry and biochemistry. The ability to design stable and functional biocatalysts that could operate under different conditions to perform chemical reactions without precedent in nature, allowing the large-scale production of chemicals à la carte, would revolutionise both synthetic, pharmacologic and materials chemistry. Despite the great advances achieved, this highly multidisciplinary area of research is still in its infancy. This chapter describes the ‘inside-out’ protocol for computational enzyme design and both the achievements and limitations of the current technology are highlighted. Furthermore, molecular dynamics simulations have proved to be invaluable in the enzyme design process, constituting an important tool for discovering elusive catalytically relevant conformations of the engineered or designed enzyme. As a complement to the ‘inside-out’ design protocol, different examples where hybrid QM/MM approaches have been directly applied to discover beneficial mutations in rational computational enzyme design are described.
Collapse
Affiliation(s)
- Marc Garcia-Borràs
- Department of Chemistry and Biochemistry, University of California Los Angeles California CA 90095-1569 USA
| | - Kendall N. Houk
- Department of Chemistry and Biochemistry, University of California Los Angeles California CA 90095-1569 USA
| | - Gonzalo Jiménez-Osés
- Departamento de Química, Centro de Investigación en Síntesis Química Universidad de La Rioja 26006 Logroño La Rioja Spain
| |
Collapse
|
17
|
Paladino A, Marchetti F, Rinaldi S, Colombo G. Protein design: from computer models to artificial intelligence. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2017. [DOI: 10.1002/wcms.1318] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Affiliation(s)
- Antonella Paladino
- Biomolecular Simulations & Computational Chemistry Group; Istituto Istituto di Chimica del Riconoscimento Molecolare, CNR; Milano Italy
| | - Filippo Marchetti
- Biomolecular Simulations & Computational Chemistry Group; Istituto Istituto di Chimica del Riconoscimento Molecolare, CNR; Milano Italy
| | - Silvia Rinaldi
- Biomolecular Simulations & Computational Chemistry Group; Istituto Istituto di Chimica del Riconoscimento Molecolare, CNR; Milano Italy
| | - Giorgio Colombo
- Biomolecular Simulations & Computational Chemistry Group; Istituto Istituto di Chimica del Riconoscimento Molecolare, CNR; Milano Italy
| |
Collapse
|
18
|
Świderek K, Tuñón I, Moliner V, Bertran J. Revealing the Origin of the Efficiency of the De Novo Designed Kemp Eliminase HG-3.17 by Comparison with the Former Developed HG-3. Chemistry 2017; 23:7582-7589. [DOI: 10.1002/chem.201700807] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/20/2017] [Indexed: 11/10/2022]
Affiliation(s)
- Katarzyna Świderek
- Departament de Química Física i Analítica; Universitat Jaume I; 12071 Castellón Spain
- Institute of Applied Radiation Chemistry; Lodz University of Technology; 90-924 Lodz Poland
| | - Iñaki Tuñón
- Departament de Química Física; Universitat de València; 46100 Burjassot Spain
| | - Vicent Moliner
- Departament de Química Física i Analítica; Universitat Jaume I; 12071 Castellón Spain
- Department of Chemistry; University of Bath; Bath BA2 7AY United Kingdom
| | - Joan Bertran
- Departament de Química; Universitat Autònoma de Barcelona; 08193 Bellaterra Spain
| |
Collapse
|
19
|
Jindal G, Ramachandran B, Bora RP, Warshel A. Exploring the Development of Ground-State Destabilization and Transition-State Stabilization in Two Directed Evolution Paths of Kemp Eliminases. ACS Catal 2017; 7:3301-3305. [PMID: 29082065 PMCID: PMC5658032 DOI: 10.1021/acscatal.7b00171] [Citation(s) in RCA: 26] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Computer-aided enzyme design presents a major challenge since in most cases it has not resulted in an impressive catalytic power. The reasons for the problems with computational design include the use of nonquantitative approaches, but they may also reflect other difficulties that are not completely obvious. Thus, it is very useful to try to learn from the trend in directed evolution experiments. Here we explore the nature of the refinement of Kemp eliminases by directed evolution, trying to gain an understanding of related requirements from computational design. The observed trend in the directed evolution refinement of KE07 and HG3 are reproduced, showing that in the case of KE07 the directed evolution leads to ground-state destabilization, whereas in the case of HG3 the directed evolution leads to transition-state stabilization. The nature of the different paths of the directed evolution is examined and discussed. The present study seems to indicate that computer-aided enzyme design may require more than calculations of the effect of single mutations and should be extended to calculations of the effect of simultaneous multiple mutations (that make a few residues preorganized effectively). However, the analysis of two known evolution paths can still be accomplished using the relevant sequences and structures. Thus, by comparing two directed evolution paths of Kemp eliminases we reached the important conclusion that the more effective path leads to transition-state stabilization.
Collapse
Affiliation(s)
| | - Balajee Ramachandran
- Department of Chemistry, University of Southern California, 3620 McClintock Avenue, Los Angeles, California 90089, United States
| | - Ram Prasad Bora
- 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
|
20
|
Bhowmick A, Sharma SC, Head-Gordon T. The Importance of the Scaffold for de Novo Enzymes: A Case Study with Kemp Eliminase. J Am Chem Soc 2017; 139:5793-5800. [DOI: 10.1021/jacs.6b12265] [Citation(s) in RCA: 45] [Impact Index Per Article: 6.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
| | | | - Teresa Head-Gordon
- Chemical
Sciences Division, Lawrence Berkeley National Laboratories, Berkeley, California 94720, United States
| |
Collapse
|
21
|
Dral P, Wu X, Spörkel L, Koslowski A, Thiel W. Semiempirical Quantum-Chemical Orthogonalization-Corrected Methods: Benchmarks for Ground-State Properties. J Chem Theory Comput 2016; 12:1097-120. [PMID: 26771261 PMCID: PMC4785506 DOI: 10.1021/acs.jctc.5b01047] [Citation(s) in RCA: 62] [Impact Index Per Article: 7.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/04/2015] [Indexed: 11/30/2022]
Abstract
The semiempirical orthogonalization-corrected OMx methods (OM1, OM2, and OM3) go beyond the standard MNDO model by including additional interactions in the electronic structure calculation. When augmented with empirical dispersion corrections, the resulting OMx-Dn approaches offer a fast and robust treatment of noncovalent interactions. Here we evaluate the performance of the OMx and OMx-Dn methods for a variety of ground-state properties using a large and diverse collection of benchmark sets from the literature, with a total of 13035 original and derived reference data. Extensive comparisons are made with the results from established semiempirical methods (MNDO, AM1, PM3, PM6, and PM7) that also use the NDDO (neglect of diatomic differential overlap) integral approximation. Statistical evaluations show that the OMx and OMx-Dn methods outperform the other methods for most of the benchmark sets.
Collapse
Affiliation(s)
- Pavlo
O. Dral
- Max-Planck-Institut für
Kohlenforschung, Kaiser-Wilhelm-Platz
1, 45470 Mülheim
an der Ruhr, Germany
| | - Xin Wu
- Max-Planck-Institut für
Kohlenforschung, Kaiser-Wilhelm-Platz
1, 45470 Mülheim
an der Ruhr, Germany
| | - Lasse Spörkel
- Max-Planck-Institut für
Kohlenforschung, Kaiser-Wilhelm-Platz
1, 45470 Mülheim
an der Ruhr, Germany
| | - Axel Koslowski
- Max-Planck-Institut für
Kohlenforschung, Kaiser-Wilhelm-Platz
1, 45470 Mülheim
an der Ruhr, Germany
| | - Walter Thiel
- Max-Planck-Institut für
Kohlenforschung, Kaiser-Wilhelm-Platz
1, 45470 Mülheim
an der Ruhr, Germany
| |
Collapse
|
22
|
Bhowmick A, Sharma SC, Honma H, Head-Gordon T. The role of side chain entropy and mutual information for improving the de novo design of Kemp eliminases KE07 and KE70. Phys Chem Chem Phys 2016; 18:19386-96. [DOI: 10.1039/c6cp03622h] [Citation(s) in RCA: 21] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Side chain entropy and mutual entropy information between residue pairs have been calculated for two de novo designed Kemp eliminase enzymes, KE07 and KE70, and for their most improved versions at the end of laboratory directed evolution (LDE).
Collapse
Affiliation(s)
- Asmit Bhowmick
- Department of Chemical and Biomolecular Engineering
- University of California Berkeley
- Berkeley
- USA
| | - Sudhir C. Sharma
- Department of Chemistry
- University of California Berkeley
- Berkeley
- USA
| | - Hallie Honma
- Department of Bioengineering, University of California Berkeley
- Berkeley
- USA
| | - Teresa Head-Gordon
- Department of Chemical and Biomolecular Engineering
- University of California Berkeley
- Berkeley
- USA
- Department of Chemistry
| |
Collapse
|
23
|
Pottel J, Moitessier N. Single-Point Mutation with a Rotamer Library Toolkit: Toward Protein Engineering. J Chem Inf Model 2015; 55:2657-71. [PMID: 26623941 DOI: 10.1021/acs.jcim.5b00525] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022]
Abstract
Protein engineers have long been hard at work to harness biocatalysts as a natural source of regio-, stereo-, and chemoselectivity in order to carry out chemistry (reactions and/or substrates) not previously achieved with these enzymes. The extreme labor demands and exponential number of mutation combinations have induced computational advances in this domain. The first step in our virtual approach is to predict the correct conformations upon mutation of residues (i.e., rebuilding side chains). For this purpose, we opted for a combination of molecular mechanics and statistical data. In this work, we have developed automated computational tools to extract protein structural information and created conformational libraries for each amino acid dependent on a variable number of parameters (e.g., resolution, flexibility, secondary structure). We have also developed the necessary tool to apply the mutation and optimize the conformation accordingly. For side-chain conformation prediction, we obtained overall average root-mean-square deviations (RMSDs) of 0.91 and 1.01 Å for the 18 flexible natural amino acids within two distinct sets of over 3000 and 1500 side-chain residues, respectively. The commonly used dihedral angle differences were also evaluated and performed worse than the state of the art. These two metrics are also compared. Furthermore, we generated a family-specific library for kinases that produced an average 2% lower RMSD upon side-chain reconstruction and a residue-specific library that yielded a 17% improvement. Ultimately, since our protein engineering outlook involves using our docking software, Fitted/Impacts, we applied our mutation protocol to a benchmarked data set for self- and cross-docking. Our side-chain reconstruction does not hinder our docking software, demonstrating differences in pose prediction accuracy of approximately 2% (RMSD cutoff metric) for a set of over 200 protein/ligand structures. Similarly, when docking to a set of over 100 kinases, side-chain reconstruction (using both general and biased conformation libraries) had minimal detriment to the docking accuracy.
Collapse
Affiliation(s)
- Joshua Pottel
- Department of Chemistry, McGill University , 801 Sherbrooke Street West, Montreal, QC, Canada H3A 0B8
| | - Nicolas Moitessier
- Department of Chemistry, McGill University , 801 Sherbrooke Street West, Montreal, QC, Canada H3A 0B8
| |
Collapse
|
24
|
Świderek K, Tuñón I, Moliner V, Bertran J. Computational strategies for the design of new enzymatic functions. Arch Biochem Biophys 2015; 582:68-79. [PMID: 25797438 PMCID: PMC4554825 DOI: 10.1016/j.abb.2015.03.013] [Citation(s) in RCA: 37] [Impact Index Per Article: 4.1] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/26/2014] [Revised: 03/09/2015] [Accepted: 03/13/2015] [Indexed: 11/28/2022]
Abstract
In this contribution, recent developments in the design of biocatalysts are reviewed with particular emphasis in the de novo strategy. Studies based on three different reactions, Kemp elimination, Diels-Alder and Retro-Aldolase, are used to illustrate different success achieved during the last years. Finally, a section is devoted to the particular case of designed metalloenzymes. As a general conclusion, the interplay between new and more sophisticated engineering protocols and computational methods, based on molecular dynamics simulations with Quantum Mechanics/Molecular Mechanics potentials and fully flexible models, seems to constitute the bed rock for present and future successful design strategies.
Collapse
Affiliation(s)
- K Świderek
- Departament de Química Física, Universitat de València, 46100 Burjasot, Spain; Institute of Applied Radiation Chemistry, Lodz University of Technology, 90-924 Lodz, Poland
| | - I Tuñón
- Departament de Química Física, Universitat de València, 46100 Burjasot, Spain
| | - V Moliner
- Departament de Química Física i Analítica, Universitat Jaume I, 12071 Castellón, Spain
| | - J Bertran
- Departament de Química, Universitat Autònoma de Barcelona, 08193 Bellaterra, Spain.
| |
Collapse
|
25
|
Brunk E, Rothlisberger U. Mixed Quantum Mechanical/Molecular Mechanical Molecular Dynamics Simulations of Biological Systems in Ground and Electronically Excited States. Chem Rev 2015; 115:6217-63. [PMID: 25880693 DOI: 10.1021/cr500628b] [Citation(s) in RCA: 292] [Impact Index Per Article: 32.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Affiliation(s)
- Elizabeth Brunk
- †Laboratory of Computational Chemistry and Biochemistry, Ecole Polytechnique Fédérale de Lausanne, 1015 Lausanne, Switzerland.,‡Joint BioEnergy Institute, Lawrence Berkeley National Laboratory, Emeryville, California 94618, United States
| | - Ursula Rothlisberger
- †Laboratory of Computational Chemistry and Biochemistry, Ecole Polytechnique Fédérale de Lausanne, 1015 Lausanne, Switzerland.,§National Competence Center of Research (NCCR) MARVEL-Materials' Revolution: Computational Design and Discovery of Novel Materials, 1015 Lausanne, Switzerland
| |
Collapse
|
26
|
Sharir-Ivry A, Varatharaj R, Shurki A. Valence bond and enzyme catalysis: a time to break down and a time to build up. Chemistry 2015; 21:7159-69. [PMID: 25808731 DOI: 10.1002/chem.201406236] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/26/2014] [Indexed: 11/12/2022]
Abstract
Understanding enzyme catalysis and developing ability to control of it are two great challenges in biochemistry. A few successful examples of computational-based enzyme design have proved the fantastic potential of computational approaches in this field, however, relatively modest rate enhancements have been reported and the further development of complementary methods is still required. Herein we propose a conceptually simple scheme to identify the specific role that each residue plays in catalysis. The scheme is based on a breakdown of the total catalytic effect into contributions of individual protein residues, which are further decomposed into chemically interpretable components by using valence bond theory. The scheme is shown to shed light on the origin of catalysis in wild-type haloalkane dehalogenase (wt-DhlA) and its mutants. Furthermore, the understanding gained through our scheme is shown to have great potential in facilitating the selection of non-optimal sites for catalysis and suggesting effective mutations to enhance the enzymatic rate.
Collapse
Affiliation(s)
- Avital Sharir-Ivry
- Institute for Drug Design, School of Pharmacy, The Lise Meitner-Minerva Center for Computational Quantum Chemistry, The Hebrew University of Jerusalem, Jerusalem 91120 (Israel), Fax: (+972) 2-675-7076
| | | | | |
Collapse
|
27
|
Świderek K, Tuñón I, Moliner V, Bertran J. Protein Flexibility and Preorganization in the Design of Enzymes. The Kemp Elimination Catalyzed by HG3.17. ACS Catal 2015. [DOI: 10.1021/cs501904w] [Citation(s) in RCA: 25] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Katarzyna Świderek
- Departament
de Química Física, Universitat de València, 46100 Burjasot, Spain
- Institute
of Applied Radiation Chemistry, Lodz University of Technology, 90-924 Lodz, Poland
| | - Iñaki Tuñón
- Departament
de Química Física, Universitat de València, 46100 Burjasot, Spain
| | - Vicent Moliner
- Departament
de Química Física i Analítica, Universitat Jaume I, 12071 Castellón, Spain
| | - Joan Bertran
- Departament
de Química, Universitat Autònoma de Barcelona, 08193 Bellaterra, Spain
| |
Collapse
|
28
|
Bora RP, Mills MJL, Frushicheva MP, Warshel A. On the Challenge of Exploring the Evolutionary Trajectory from Phosphotriesterase to Arylesterase Using Computer Simulations. J Phys Chem B 2015; 119:3434-45. [DOI: 10.1021/jp5124025] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/01/2023]
Affiliation(s)
- Ram Prasad Bora
- Department
of Chemistry, University of Southern California, 3620 McClintock Avenue, Los Angeles, California 90089, United States
| | - Matthew J. L. Mills
- Deconstruction Division, Joint BioEnergy Institute, Emeryville, California, United States & Biomass Science and Conversion Technology Department, Sandia National Laboratories, Livermore, California 94550, United States
| | - Maria P. Frushicheva
- Department
of Chemical Engineering, Massachusetts Institute of Technology, 25 Ames
Street, Cambridge, Massachusetts 02142, United States
| | - Arieh Warshel
- Department
of Chemistry, University of Southern California, 3620 McClintock Avenue, Los Angeles, California 90089, United States
| |
Collapse
|
29
|
Qiao Y, Han K, Zhan CG. Reaction pathways and free energy profiles for cholinesterase-catalyzed hydrolysis of 6-monoacetylmorphine. Org Biomol Chem 2014; 12:2214-27. [PMID: 24595354 DOI: 10.1039/c3ob42464b] [Citation(s) in RCA: 28] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
As the most active metabolite of heroin, 6-monoacetylmorphine (6-MAM) can penetrate into the brain for the rapid onset of heroin effects. The primary enzymes responsible for the metabolism of 6-MAM to the less potent morphine in humans are acetylcholinesterase (AChE) and butyrylcholinesterase (BChE). The detailed reaction pathways for AChE- and BChE-catalyzed hydrolysis of 6-MAM to morphine have been explored, for the first time, in the present study by performing first-principles quantum mechanical/molecular mechanical free energy calculations. It has been demonstrated that the two enzymatic reaction processes follow similar catalytic reaction mechanisms, and the whole catalytic reaction pathway for each enzyme consists of four reaction steps. According to the calculated results, the second reaction step associated with the transition state TS2(a)/TS2(b) should be rate-determining for the AChE/BChE-catalyzed hydrolysis, and the free energy barrier calculated for the AChE-catalyzed hydrolysis (18.3 kcal mol(-1)) is 2.5 kcal mol(-1) lower than that for the BChE-catalyzed hydrolysis (20.8 kcal mol(-1)). The free energy barriers calculated for the AChE- and BChE-catalyzed reactions are in good agreement with the experimentally derived activation free energies (17.5 and 20.7 kcal mol(-1) for the AChE- and BChE-catalyzed reactions, respectively). Further structural analysis reveals that the aromatic residues Phe295 and Phe297 in the acyl pocket of AChE (corresponding to Leu286 and Val288 in BChE) contribute to the lower energy of TS2(a) relative to TS2(b). The obtained structural and mechanistic insights could be valuable for use in future rational design of a novel therapeutic treatment of heroin abuse.
Collapse
Affiliation(s)
- Yan Qiao
- State Key Laboratory of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academy of Science, Zhongshan Road 457, Dalian 116023, P. R. China
| | | | | |
Collapse
|
30
|
Sirin S, Pearlman DA, Sherman W. Physics-based enzyme design: predicting binding affinity and catalytic activity. Proteins 2014; 82:3397-409. [PMID: 25243583 DOI: 10.1002/prot.24694] [Citation(s) in RCA: 34] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/12/2014] [Revised: 09/04/2014] [Accepted: 09/10/2014] [Indexed: 11/11/2022]
Abstract
Computational enzyme design is an emerging field that has yielded promising success stories, but where numerous challenges remain. Accurate methods to rapidly evaluate possible enzyme design variants could provide significant value when combined with experimental efforts by reducing the number of variants needed to be synthesized and speeding the time to reach the desired endpoint of the design. To that end, extending our computational methods to model the fundamental physical-chemical principles that regulate activity in a protocol that is automated and accessible to a broad population of enzyme design researchers is essential. Here, we apply a physics-based implicit solvent MM-GBSA scoring approach to enzyme design and benchmark the computational predictions against experimentally determined activities. Specifically, we evaluate the ability of MM-GBSA to predict changes in affinity for a steroid binder protein, catalytic turnover for a Kemp eliminase, and catalytic activity for α-Gliadin peptidase variants. Using the enzyme design framework developed here, we accurately rank the most experimentally active enzyme variants, suggesting that this approach could provide enrichment of active variants in real-world enzyme design applications.
Collapse
Affiliation(s)
- Sarah Sirin
- Department of Biology, Massachusetts Institute of Technology, Cambridge, Massachusetts, 02140
| | | | | |
Collapse
|
31
|
Singh MK, Chu ZT, Warshel A. Simulating the catalytic effect of a designed mononuclear zinc metalloenzyme that catalyzes the hydrolysis of phosphate triesters. J Phys Chem B 2014; 118:12146-52. [PMID: 25233046 PMCID: PMC4207531 DOI: 10.1021/jp507592g] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
![]()
One
of the greatest challenges in biotechnology and in biochemistry
is the ability to design efficient enzymes. In fact, such an ability
would be one of the most convincing manifestations of a full understanding
of the origin of enzyme catalysis. Despite some progress on this front,
most of the advances have been made by placing the reacting fragments
in the proper places rather than by optimizing the preorganization
of the environment, which is the key factor in enzyme catalysis. A
rational improvement of the preorganization and a consistent assessment
of the effectiveness of different design options require approaches
capable of evaluating reliably the actual catalytic effect. In this
work we examine the ability of the empirical valence bond (EVB) to
reproduce the results of directed evolution improvements of the catalysis
of diethyl 7-hydroxycoumarinyl by a designed mononuclear zinc metalloenzyme.
Encouragingly, our study reproduced the catalytic effect obtained
by directed evolution and offers a good start for further studies
of this system.
Collapse
Affiliation(s)
- Manoj Kumar Singh
- Department of Chemistry, University of Southern California , SGM 418, 3620 McClintock Avenue, Los Angeles, California 90089, United States
| | | | | |
Collapse
|
32
|
Frushicheva MP, Mills MJL, Schopf P, Singh MK, Warshel A. Computer aided enzyme design and catalytic concepts. Curr Opin Chem Biol 2014; 21:56-62. [PMID: 24814389 PMCID: PMC4149935 DOI: 10.1016/j.cbpa.2014.03.022] [Citation(s) in RCA: 71] [Impact Index Per Article: 7.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/12/2014] [Revised: 03/19/2014] [Accepted: 03/26/2014] [Indexed: 11/21/2022]
Abstract
Gaining a deeper understanding of enzyme catalysis is of great practical and fundamental importance. Over the years it has become clear that despite advances made in experimental mutational studies, a quantitative understanding of enzyme catalysis will not be possible without the use of computer modeling approaches. While we believe that electrostatic preorganization is by far the most important catalytic factor, convincing the wider scientific community of this may require the demonstration of effective rational enzyme design. Here we make the point that the main current advances in enzyme design are basically advances in directed evolution and that computer aided enzyme design must involve approaches that can reproduce catalysis in well-defined test cases. Such an approach is provided by the empirical valence bond method.
Collapse
Affiliation(s)
- Maria P. Frushicheva
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, MA 02142, USA
| | - Matthew J. L. Mills
- Department of Chemistry, University of Southern California, Los Angeles, CA 90089, USA
| | - Patrick Schopf
- Department of Chemistry, University of Southern California, Los Angeles, CA 90089, USA
| | - Manoj K. Singh
- Department of Chemistry, University of Southern California, Los Angeles, CA 90089, USA
| | - Arieh Warshel
- Department of Chemistry, University of Southern California, Los Angeles, CA 90089, USA
| |
Collapse
|
33
|
Li Q, Huang X, Zhu Y. Evaluation of active designs of cephalosporin C acylase by molecular dynamics simulation and molecular docking. J Mol Model 2014; 20:2314. [PMID: 24935111 DOI: 10.1007/s00894-014-2314-5] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/26/2013] [Accepted: 05/19/2014] [Indexed: 01/04/2023]
Abstract
Optimization to identify the global minimum energy conformation sequence in in silico enzyme design is computationally non-deterministic polynomial-time (NP)-hard, with the search time growing exponentially as the number of design sites increases. This drawback forces the modeling of protein-ligand systems to adopt discrete amino acid rotamers and ligand conformers, as well as continuum solvent treatment of the environment; however, such compromises produce large numbers of false positives in sequence selection. In this report, cephalosporin acylase, which catalyzes the hydrolytic reaction of cephalosporin C to 7-aminocephalosporanic acid, was used to investigate the dynamic features of active-site-transition-state complex structures using molecular dynamics (MD) simulations to potentially eliminate false positives. The molecular docking between cephalosporin C and wild type acylase N176 and its eight mutants showed that the rate-limiting step in the hydrolytic reaction of cephalosporin C is the acylation process. MD simulations of the active-site-transition-state complex structures of the acylation processes for N176 and its eight mutants showed that the geometrical constraints between catalytic residues and small molecule transition states are always well maintained during the 20 ns simulation for mutants with higher activities, and more hydrogen bonds between binding residues and functional groups of the ligand side chain in the active pocket are formed for mutants with higher activities. The conformations of the ligand transition states were changed greatly after the simulation. This indicates that the hydrogen bond network between the ligand and protein could be improved to enhance the activity of cephalosporin C acylase in subsequent design.
Collapse
Affiliation(s)
- Qing Li
- Department of Chemical Engineering, Tsinghua University, Beijing, 100084, People's .Republic of China
| | | | | |
Collapse
|
34
|
Zayner JP, Antoniou C, French AR, Hause RJ, Sosnick TR. Investigating models of protein function and allostery with a widespread mutational analysis of a light-activated protein. Biophys J 2014; 105:1027-36. [PMID: 23972854 DOI: 10.1016/j.bpj.2013.07.010] [Citation(s) in RCA: 43] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/30/2013] [Revised: 06/17/2013] [Accepted: 07/01/2013] [Indexed: 10/26/2022] Open
Abstract
To investigate the relationship between a protein's sequence and its biophysical properties, we studied the effects of more than 100 mutations in Avena sativa light-oxygen-voltage domain 2, a model protein of the Per-Arnt-Sim family. The A. sativa light-oxygen-voltage domain 2 undergoes a photocycle with a conformational change involving the unfolding of the terminal helices. Whereas selection studies typically search for winners in a large population and fail to characterize many sites, we characterized the biophysical consequences of mutations throughout the protein using NMR, circular dichroism, and ultraviolet/visible spectroscopy. Despite our intention to introduce highly disruptive substitutions, most had modest or no effect on function, and many could even be considered to be more photoactive. Substitutions at evolutionarily conserved sites can have minimal effect, whereas those at nonconserved positions can have large effects, contrary to the view that the effects of mutations, especially at conserved positions, are predictable. Using predictive models, we found that the effects of mutations on biophysical function and allostery reflect a complex mixture of multiple characteristics including location, character, electrostatics, and chemistry.
Collapse
Affiliation(s)
- Josiah P Zayner
- Department of Biochemistry and Molecular Biology, The University of Chicago, Chicago, IL, USA
| | | | | | | | | |
Collapse
|
35
|
Acevedo O, Jorgensen WL. Quantum and Molecular Mechanical (QM/MM) Monte Carlo Techniques for Modeling Condensed-Phase Reactions. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2014; 4:422-435. [PMID: 25431625 DOI: 10.1002/wcms.1180] [Citation(s) in RCA: 29] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/12/2022]
Abstract
A recent review (Acc. Chem. Res. 2010, 43:142-151) examined our use and development of a combined quantum and molecular mechanical (QM/MM) technique for modelling organic and enzymatic reactions. Advances included the PDDG/PM3 semiempirical QM (SQM) method, computation of multi-dimensional potentials of mean force (PMF), incorporation of on-the-fly QM in Monte Carlo simulations, and a polynomial quadrature method for rapidly treating proton-transfer reactions. The current article serves as a follow up on our progress. Highlights include new reactions, alternative SQM methods, a polarizable OPLS force field, and novel solvent environments, e.g., "on water" and room temperature ionic liquids. The methodology is strikingly accurate across a wide range of condensed-phase and antibody-catalyzed reactions including substitution, decarboxylation, elimination, isomerization, and pericyclic classes. Comparisons are made to systems treated with continuum-based solvents and ab initio or density functional theory (DFT) methods. Overall, the QM/MM methodology provides detailed characterization of reaction paths, proper configurational sampling, several advantages over implicit solvent models, and a reasonable computational cost.
Collapse
Affiliation(s)
- Orlando Acevedo
- Department of Chemistry and Biochemistry, Auburn University, Auburn, Alabama 36849
| | - Wiliiam L Jorgensen
- Department of Chemistry, Yale University, 225 Prospect Street, New Haven, Connecticut 06520-8107
| |
Collapse
|
36
|
Armacost K, Acevedo O. Exploring the Aldol Reaction using Catalytic Antibodies and “On Water” Organocatalysts from QM/MM Calculations. J Am Chem Soc 2013; 136:147-56. [DOI: 10.1021/ja405614p] [Citation(s) in RCA: 22] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/28/2023]
Affiliation(s)
- Kira Armacost
- Department
of Chemistry,
Biochemistry, Auburn University, Auburn, Alabama 36849, United States
| | - Orlando Acevedo
- Department
of Chemistry,
Biochemistry, Auburn University, Auburn, Alabama 36849, United States
| |
Collapse
|
37
|
Abstract
Diverse engineering strategies have been developed to create enzymes with novel catalytic activities. Among these, computational approaches hold particular promise. Enzymes have been computationally designed to promote several nonbiological reactions, including a Diels-Alder cycloaddition, proton transfer, multistep retroaldol transformations, and metal-dependent hydrolysis of phosphotriesters. Although their efficiencies (kcat/KM = 0.1-100 M(-1) s(-1)) are typically low compared with those of the best natural enzymes (10(6)-10(8) M(-1) s(-1)), these catalysts are excellent starting points for laboratory evolution. This review surveys recent progress in combining computational and evolutionary approaches to enzyme design, together with insights into enzyme function gained from studies of the engineered catalysts.
Collapse
Affiliation(s)
- Donald Hilvert
- Laboratory of Organic Chemistry, ETH Zürich, 8093 Zürich, Switzerland.
| |
Collapse
|
38
|
Kellie JL, Wilson KA, Wetmore SD. Standard role for a conserved aspartate or more direct involvement in deglycosylation? An ONIOM and MD investigation of adenine-DNA glycosylase. Biochemistry 2013; 52:8753-65. [PMID: 24168684 DOI: 10.1021/bi401310w] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Abstract
8-Oxoguanine (OG) is one of the most frequently occurring forms of DNA damage and is particularly deleterious since it forms a stable Hoogsteen base pair with adenine (A). The repair of an OG:A mispair is initiated by adenine-DNA glycosylase (MutY), which hydrolyzes the sugar-nucleobase bond of the adenine residue before the lesion is processed by other proteins. MutY has been proposed to use a two-part chemical step involving protonation of the adenine nucleobase, followed by SN1 hydrolysis of the glycosidic bond. However, differences between a recent (fluorine recognition complex, denoted as the FLRC) crystal structure and the structure on which most mechanistic conclusions have been based to date (namely, the lesion recognition complex or LRC) raise questions regarding the mechanism used by MutY and the discrete role of various active-site residues. The present work uses both molecular dynamics (MD) and quantum mechanical (ONIOM) models to compare the active-site conformational dynamics in the two crystal structures, which suggests that only the understudied FLRC leads to a catalytically competent reactant. Indeed, all previous computational studies on MutY have been initiated from the LRC structure. Subsequently, for the first time, various mechanisms are examined with detailed ONIOM(M06-2X:PM6) reaction potential energy surfaces (PES) based on the FLRC structure, which significantly extends the mechanistic picture. Specifically, our work reveals that the reaction proceeds through a different route than the commonly accepted mechanism and the catalytic function of various active-site residues (Geobacillus stearothermophilus numbering). Specifically, contrary to proposals based on the LRC, E43 is determined to solely be involved in the initial adenine protonation step and not the deglycosylation reaction as the general base. Additionally, a novel catalytic role is proposed for Y126, whereby this residue plays a significant role in stabilizing the highly charged active site, primarily through interactions with E43. More importantly, D144 is found to explicitly catalyze the nucleobase dissociation step through partial nucleophilic attack. Although this is a more direct role than previously proposed for any other DNA glycosylase, comparison to previous work on other glycosylases justifies the larger contribution in the case of MutY and allows us to propose a unified role for the conserved Asp/Glu in the DNA glycosylases, as well as other enzymes that catalyze nucleotide deglycosylation reactions.
Collapse
Affiliation(s)
- Jennifer L Kellie
- Department of Chemistry and Biochemistry, University of Lethbridge , 4401 University Drive West, Lethbridge, Alberta, Canada T1K 3M4
| | | | | |
Collapse
|
39
|
Wei D, Fang L, Tang M, Zhan CG. Fundamental reaction pathway for peptide metabolism by proteasome: insights from first-principles quantum mechanical/molecular mechanical free energy calculations. J Phys Chem B 2013; 117:13418-34. [PMID: 24111489 DOI: 10.1021/jp405337v] [Citation(s) in RCA: 33] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2023]
Abstract
Proteasome is the major component of the crucial non-lysosomal protein degradation pathway in the cells, but the detailed reaction pathway is unclear. In this study, first-principles quantum mechanical/molecular mechanical free energy calculations have been performed to explore, for the first time, possible reaction pathways for proteasomal proteolysis/hydrolysis of a representative peptide, succinyl-leucyl-leucyl-valyl-tyrosyl-7-amino-4-methylcoumarin (Suc-LLVY-AMC). The computational results reveal that the most favorable reaction pathway consists of six steps. The first is a water-assisted proton transfer within proteasome, activating Thr1-O(γ). The second is a nucleophilic attack on the carbonyl carbon of a Tyr residue of substrate by the negatively charged Thr1-O(γ), followed by the dissociation of the amine AMC (third step). The fourth step is a nucleophilic attack on the carbonyl carbon of the Tyr residue of substrate by a water molecule, accompanied by a proton transfer from the water molecule to Thr1-N(z). Then, Suc-LLVY is dissociated (fifth step), and Thr1 is regenerated via a direct proton transfer from Thr1-N(z) to Thr1-O(γ). According to the calculated energetic results, the overall reaction energy barrier of the proteasomal hydrolysis is associated with the transition state (TS3(b)) for the third step involving a water-assisted proton transfer. The determined most favorable reaction pathway and the rate-determining step have provided a reasonable interpretation of the reported experimental observations concerning the substituent and isotopic effects on the kinetics. The calculated overall free energy barrier of 18.2 kcal/mol is close to the experimentally derived activation free energy of ∼18.3-19.4 kcal/mol, suggesting that the computational results are reasonable.
Collapse
Affiliation(s)
- Donghui Wei
- Department of Chemistry, Zhengzhou University , 75 Daxue Road, Zhengzhou, Henan 450052, China
| | | | | | | |
Collapse
|
40
|
Qiao Y, Han K, Zhan CG. Fundamental reaction pathway and free energy profile for butyrylcholinesterase-catalyzed hydrolysis of heroin. Biochemistry 2013; 52:6467-79. [PMID: 23992153 DOI: 10.1021/bi400709v] [Citation(s) in RCA: 34] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/07/2023]
Abstract
The pharmacological function of heroin requires an activation process that transforms heroin into 6-monoacetylmorphine (6-MAM), which is the most active form. The primary enzyme responsible for this activation process in human plasma is butyrylcholinesterase (BChE). The detailed reaction pathway of the activation process via BChE-catalyzed hydrolysis has been explored computationally, for the first time, in this study via molecular dynamics simulation and first-principles quantum mechanical/molecular mechanical free energy calculations. It has been demonstrated that the whole reaction process includes acylation and deacylation stages. The acylation consists of two reaction steps, i.e., the nucleophilic attack on the carbonyl carbon of the 3-acetyl group of heroin by the hydroxyl oxygen of the Ser198 side chain and the dissociation of 6-MAM. The deacylation also consists of two reaction steps, i.e., the nucleophilic attack on the carbonyl carbon of the acyl-enzyme intermediate by a water molecule and the dissociation of the acetic acid from Ser198. The calculated free energy profile reveals that the second transition state (TS2) should be rate-determining. The structural analysis reveals that the oxyanion hole of BChE plays an important role in the stabilization of rate-determining TS2. The free energy barrier (15.9 ± 0.2 or 16.1 ± 0.2 kcal/mol) calculated for the rate-determining step is in good agreement with the experimentally derived activation free energy (~16.2 kcal/mol), suggesting that the mechanistic insights obtained from this computational study are reliable. The obtained structural and mechanistic insights could be valuable for use in the future rational design of a novel therapeutic treatment of heroin abuse.
Collapse
Affiliation(s)
- Yan Qiao
- State Key Laboratory of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academy of Science , Zhongshan Road 457, Dalian 116023, P. R. China
| | | | | |
Collapse
|
41
|
Duarte F, Amrein BA, Kamerlin SCL. Modeling catalytic promiscuity in the alkaline phosphatase superfamily. Phys Chem Chem Phys 2013; 15:11160-77. [PMID: 23728154 PMCID: PMC3693508 DOI: 10.1039/c3cp51179k] [Citation(s) in RCA: 43] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/18/2013] [Accepted: 05/02/2013] [Indexed: 12/19/2022]
Abstract
In recent years, it has become increasingly clear that promiscuity plays a key role in the evolution of new enzyme function. This finding has helped to elucidate fundamental aspects of molecular evolution. While there has been extensive experimental work on enzyme promiscuity, computational modeling of the chemical details of such promiscuity has traditionally fallen behind the advances in experimental studies, not least due to the nearly prohibitive computational cost involved in examining multiple substrates with multiple potential mechanisms and binding modes in atomic detail with a reasonable degree of accuracy. However, recent advances in both computational methodologies and power have allowed us to reach a stage in the field where we can start to overcome this problem, and molecular simulations can now provide accurate and efficient descriptions of complex biological systems with substantially less computational cost. This has led to significant advances in our understanding of enzyme function and evolution in a broader sense. Here, we will discuss currently available computational approaches that can allow us to probe the underlying molecular basis for enzyme specificity and selectivity, discussing the inherent strengths and weaknesses of each approach. As a case study, we will discuss recent computational work on different members of the alkaline phosphatase superfamily (AP) using a range of different approaches, showing the complementary insights they have provided. We have selected this particular superfamily, as it poses a number of significant challenges for theory, ranging from the complexity of the actual reaction mechanisms involved to the reliable modeling of the catalytic metal centers, as well as the very large system sizes. We will demonstrate that, through current advances in methodologies, computational tools can provide significant insight into the molecular basis for catalytic promiscuity, and, therefore, in turn, the mechanisms of protein functional evolution.
Collapse
Affiliation(s)
- Fernanda Duarte
- Uppsala University, Science for Life Laboratory (SciLifeLab), Cell and Molecular Biology, Uppsala, Sweden. ; ;
| | - Beat Anton Amrein
- Uppsala University, Science for Life Laboratory (SciLifeLab), Cell and Molecular Biology, Uppsala, Sweden. ; ;
| | | |
Collapse
|
42
|
Wijma HJ, Janssen DB. Computational design gains momentum in enzyme catalysis engineering. FEBS J 2013; 280:2948-60. [DOI: 10.1111/febs.12324] [Citation(s) in RCA: 53] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/20/2013] [Revised: 04/19/2013] [Accepted: 04/24/2013] [Indexed: 01/19/2023]
Affiliation(s)
- Hein J. Wijma
- Department of Biochemistry; Groningen Biomolecular Sciences and Biotechnology Institute; University of Groningen; The Netherlands
| | - Dick B. Janssen
- Department of Biochemistry; Groningen Biomolecular Sciences and Biotechnology Institute; University of Groningen; The Netherlands
| |
Collapse
|
43
|
van der Kamp MW, Mulholland AJ. Combined quantum mechanics/molecular mechanics (QM/MM) methods in computational enzymology. Biochemistry 2013; 52:2708-28. [PMID: 23557014 DOI: 10.1021/bi400215w] [Citation(s) in RCA: 395] [Impact Index Per Article: 35.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Abstract
Computational enzymology is a rapidly maturing field that is increasingly integral to understanding mechanisms of enzyme-catalyzed reactions and their practical applications. Combined quantum mechanics/molecular mechanics (QM/MM) methods are important in this field. By treating the reacting species with a quantum mechanical method (i.e., a method that calculates the electronic structure of the active site) and including the enzyme environment with simpler molecular mechanical methods, enzyme reactions can be modeled. Here, we review QM/MM methods and their application to enzyme-catalyzed reactions to investigate fundamental and practical problems in enzymology. A range of QM/MM methods is available, from cheaper and more approximate methods, which can be used for molecular dynamics simulations, to highly accurate electronic structure methods. We discuss how modeling of reactions using such methods can provide detailed insight into enzyme mechanisms and illustrate this by reviewing some recent applications. We outline some practical considerations for such simulations. Further, we highlight applications that show how QM/MM methods can contribute to the practical development and application of enzymology, e.g., in the interpretation and prediction of the effects of mutagenesis and in drug and catalyst design.
Collapse
Affiliation(s)
- Marc W van der Kamp
- Centre for Computational Chemistry, School of Chemistry, University of Bristol, Bristol BS8 1TS, UK.
| | | |
Collapse
|
44
|
Kries H, Blomberg R, Hilvert D. De novo enzymes by computational design. Curr Opin Chem Biol 2013; 17:221-8. [DOI: 10.1016/j.cbpa.2013.02.012] [Citation(s) in RCA: 194] [Impact Index Per Article: 17.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/04/2012] [Revised: 02/08/2013] [Accepted: 02/08/2013] [Indexed: 11/29/2022]
|
45
|
Kiss G, Çelebi-Ölçüm N, Moretti R, Baker D, Houk KN. Computational enzyme design. Angew Chem Int Ed Engl 2013; 52:5700-25. [PMID: 23526810 DOI: 10.1002/anie.201204077] [Citation(s) in RCA: 370] [Impact Index Per Article: 33.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/25/2012] [Indexed: 11/07/2022]
Abstract
Recent developments in computational chemistry and biology have come together in the "inside-out" approach to enzyme engineering. Proteins have been designed to catalyze reactions not previously accelerated in nature. Some of these proteins fold and act as catalysts, but the success rate is still low. The achievements and limitations of the current technology are highlighted and contrasted to other protein engineering techniques. On its own, computational "inside-out" design can lead to the production of catalytically active and selective proteins, but their kinetic performances fall short of natural enzymes. When combined with directed evolution, molecular dynamics simulations, and crowd-sourced structure-prediction approaches, however, computational designs can be significantly improved in terms of binding, turnover, and thermal stability.
Collapse
Affiliation(s)
- Gert Kiss
- Department of Chemistry and Biochemistry, University of California, Los Angeles, 607 Charles E. Young Dr. East, Los Angeles, CA 90095, USA
| | | | | | | | | |
Collapse
|
46
|
Kiss G, Çelebi-Ölçüm N, Moretti R, Baker D, Houk KN. Computerbasiertes Enzymdesign. Angew Chem Int Ed Engl 2013. [DOI: 10.1002/ange.201204077] [Citation(s) in RCA: 43] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
|
47
|
Abrahamson MJ, Wong JW, Bommarius AS. The Evolution of an Amine Dehydrogenase Biocatalyst for the Asymmetric Production of Chiral Amines. Adv Synth Catal 2013. [DOI: 10.1002/adsc.201201030] [Citation(s) in RCA: 110] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022]
|
48
|
Masterson JE, Schwartz SD. Changes in protein architecture and subpicosecond protein dynamics impact the reaction catalyzed by lactate dehydrogenase. J Phys Chem A 2013; 117:7107-13. [PMID: 23441954 DOI: 10.1021/jp400376h] [Citation(s) in RCA: 36] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
We have previously established the importance of a promoting vibration, a subpicosecond protein motion that propagates through a specific axis of residues, in the reaction coordinate of lactate dehydrogenase (LDH). To test the effect that perturbation of this motion would have on the enzymatic reaction, we employ transition path sampling to obtain transition path ensembles for four independent LDH enzymatic systems: the wild type enzyme, a version of the enzyme expressing heavy isotopic substitution, and two enzymes with mutations in the promoting vibration axis. We show that even slight changes in the promoting vibration of LDH result in dramatic changes in enzymatic chemistry. In the "heavy" version of the enzyme, we find that the dampening of the subpicosecond dynamics from heavy isotopic substitution leads to a drastic increase in the time of barrier crossing. Furthermore, we see that mutation of the promoting vibration axis causes a decrease in the variability of transition paths available to the enzymatic reaction. The combined results reveal the importance of the protein architecture of LDH in enzymatic catalysis by establishing how the promoting vibration is finely tuned to facilitate chemistry.
Collapse
Affiliation(s)
- Jean E Masterson
- Department of Chemistry and Biochemistry, University of Arizona, 1306 East University Blvd., Tucson, Arizona 85721, USA
| | | |
Collapse
|
49
|
Sparta M, Shirvanyants D, Ding F, Dokholyan NV, Alexandrova AN. Hybrid dynamics simulation engine for metalloproteins. Biophys J 2013; 103:767-76. [PMID: 22947938 DOI: 10.1016/j.bpj.2012.06.024] [Citation(s) in RCA: 24] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/12/2012] [Revised: 06/08/2012] [Accepted: 06/18/2012] [Indexed: 11/15/2022] Open
Abstract
Quality computational description of metalloproteins is a great challenge due to the vast span of time- and lengthscales characteristic of their existence. We present an efficient new method that allows for robust characterization of metalloproteins. It combines quantum mechanical (QM) description of the metal-containing active site, and extensive dynamics of the protein captured by discrete molecular dynamics (DMD) (QM/DMD). DMD samples the entire protein, including the backbone, and most of the active site, except for the immediate coordination region of the metal. QM operates on the part of the protein of electronic and chemical significance, which may include tens to hundreds of atoms. The breathing quantum-classical boundary provides a continuous mutual feedback between the two machineries. We test QM/DMD using the Fe-containing electron transporter protein, rubredoxin, and its three mutants as a model. QM/DMD can provide a reliable balanced description of metalloproteins' structure, dynamics, and electronic structure in a reasonable amount of time. As an illustration of QM/DMD capabilities, we then predict the structure of the Ca(2+) form of the enzyme catechol O-methyl transferase, which, unlike the native Mg(2+) form, is catalytically inactive. The Mg(2+) site is ochtahedral, but the Ca(2+) is 7-coordinate and features the misalignment of the reacting parts of the system. The change is facilitated by the backbone adjustment. QM/DMD is ideal and fast for providing this level of structural insight.
Collapse
Affiliation(s)
- Manuel Sparta
- Department of Chemistry and Biochemistry, University of California, Los Angeles, California, USA
| | | | | | | | | |
Collapse
|
50
|
Emerging themes in the computational design of novel enzymes and protein-protein interfaces. FEBS Lett 2012; 587:1147-54. [DOI: 10.1016/j.febslet.2012.12.009] [Citation(s) in RCA: 44] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/19/2012] [Revised: 12/10/2012] [Accepted: 12/10/2012] [Indexed: 11/21/2022]
|