1
|
Chen L, Yu K, Ma A, Zhu W, Wang H, Tang X, Tang Y, Li Y, Li J. Enhanced Thermostability of Nattokinase by Computation-Based Rational Redesign of Flexible Regions. JOURNAL OF AGRICULTURAL AND FOOD CHEMISTRY 2024; 72:14241-14254. [PMID: 38864682 DOI: 10.1021/acs.jafc.4c02335] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/13/2024]
Abstract
Nattokinase is a nutrient in healthy food natto that has the function of preventing and treating blood thrombus. However, its low thermostability and fibrinolytic activity limit its application in food and pharmaceuticals. In this study, we used bioinformatics analysis to identify two loops (loop10 and loop12) in the flexible region of nattokinase rAprY. Using this basis, we screened the G131S-S161T variant, which showed a 2.38-fold increase in half-life at 55 °C, and the M3 variant, which showed a 2.01-fold increase in activity, by using a thermostability prediction algorithm. Bioinformatics analysis revealed that the enhanced thermostability of the G131S-S161T variant was due to the increased rigidity and structural shrinkage of the overall structure. Additionally, the increased rigidity of the local region surrounding the active center and its mutated sites helps maintain its normal conformation in high-temperature environments. The increased catalytic activity of the M3 variant may be due to its more efficient substrate binding mechanism. We investigated strategies to improve the thermostability and fibrinolytic activity of nattokinase, and the resulting variants show promise for industrial production and application.
Collapse
Affiliation(s)
- Liangqi Chen
- Institute of Materia Medica, College of Pharmacy, Xinjiang University, Urumqi 830017, China
- Xinjiang Key Laboratory of Biological Resources and Genetic Engineering, College of Life Science and Technology, Xinjiang University, Urumqi 830017, China
| | - Kongfang Yu
- Institute of Materia Medica, College of Pharmacy, Xinjiang University, Urumqi 830017, China
| | - Aixia Ma
- Institute of Materia Medica, College of Pharmacy, Xinjiang University, Urumqi 830017, China
- Xinjiang Key Laboratory of Biological Resources and Genetic Engineering, College of Life Science and Technology, Xinjiang University, Urumqi 830017, China
| | - Wenhui Zhu
- Institute of Materia Medica, College of Pharmacy, Xinjiang University, Urumqi 830017, China
| | - Hong Wang
- Institute of Materia Medica, College of Pharmacy, Xinjiang University, Urumqi 830017, China
| | - Xiyu Tang
- Institute of Materia Medica, College of Pharmacy, Xinjiang University, Urumqi 830017, China
- Xinjiang Key Laboratory of Biological Resources and Genetic Engineering, College of Life Science and Technology, Xinjiang University, Urumqi 830017, China
| | - Yaolei Tang
- Xinjiang Key Laboratory of Biological Resources and Genetic Engineering, College of Life Science and Technology, Xinjiang University, Urumqi 830017, China
- The Third People's Hospital of Xinjiang Uygur Autonomous Region, Urumqi 830000, China
| | - Yuan Li
- Institute of Materia Medica, College of Pharmacy, Xinjiang University, Urumqi 830017, China
- Xinjiang Key Laboratory of Biological Resources and Genetic Engineering, College of Life Science and Technology, Xinjiang University, Urumqi 830017, China
| | - Jinyao Li
- Institute of Materia Medica, College of Pharmacy, Xinjiang University, Urumqi 830017, China
- Xinjiang Key Laboratory of Biological Resources and Genetic Engineering, College of Life Science and Technology, Xinjiang University, Urumqi 830017, China
| |
Collapse
|
2
|
Coricello A, Nardone AJ, Lupia A, Gratteri C, Vos M, Chaptal V, Alcaro S, Zhu W, Takagi Y, Richards NGJ. Cryo-EM and Molecular Dynamics Simulations Reveal Hidden Conformational Dynamics Controlling Ammonia Transport in Human Asparagine Synthetase. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2024:2023.05.16.541009. [PMID: 37292727 PMCID: PMC10245805 DOI: 10.1101/2023.05.16.541009] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/10/2023]
Abstract
How motions in enzymes might be linked to catalytic function is of considerable general interest. Recent advances in X-ray crystallography and cryogenic electron microscopy offer the promise of elucidating functionally relevant motions in proteins that are not easily amenable to study by other biophysical methods. Here we use 3D variability analysis (3DVA) on cryo-EM maps for wild type (WT) human asparagine synthetase (ASNS) and the R142I ASNS variant to identify conformational changes in the Arg-142 side chain, which mediates the formation of a catalytically relevant intramolecular tunnel. Our 3DVA results for WT ASNS are consistent with independent molecular dynamics (MD) simulations on a model generated from the X-ray structure of human ASNS. Moreover, MD simulations of computational models for the ASNS/β-aspartyl-AMP/MgPPi and R142I/β-aspartyl-AMP/MgPPi ternary complexes, suggest that the structural integrity of the tunnel is impaired in the R142I variant when β-aspartyl-AMP is present in the synthetase active site. The kinetic properties of the R142I ASNS variant support the proposed function of Arg-142. These studies illustrate the power of cryo-EM to identify localized motions and dissect the conformational landscape of large proteins. When combined with MD simulations, 3DVA is a powerful approach to understanding how conformational dynamics might regulate function in multi-domain enzymes possessing multiple active sites.
Collapse
Affiliation(s)
- Adriana Coricello
- Dipartimento di Scienze della Salute, Università "Magna Græcia" di Catanzaro, 88100 Catanzaro, Italy
- Present address: Dipartimento di Scienze Biomolecolari, Università degli Studi di Urbino "Carlo Bo", 61029 Urbino, Italy
| | - Alanya J Nardone
- Department of Chemistry & Biochemistry, Florida State University, Tallahassee, FL 32306, USA
| | - Antonio Lupia
- Net4Science Academic Spin-Off, Università "Magna Græcia" di Catanzaro, 88100 Catanzaro, Italy
- Present address: Dipartimento di Scienze della vita e dell'ambiente, Università degli Studi di Cagliari, 09042 Cagliari, Italy
| | - Carmen Gratteri
- Dipartimento di Scienze della Salute, Università "Magna Græcia" di Catanzaro, 88100 Catanzaro, Italy
| | - Matthijn Vos
- NanoImaging Core Facility, Centre de Resources et Recherches Technologiques, Institut Pasteur, 75015 Paris, France
| | - Vincent Chaptal
- Molecular Microbiology and Structural Biochemistry Laboratory, CNRS UMR 5086, University of Lyon, 69367 Lyon, France
| | - Stefano Alcaro
- Dipartimento di Scienze della Salute, Università "Magna Græcia" di Catanzaro, 88100 Catanzaro, Italy
- Net4Science Academic Spin-Off, Università "Magna Græcia" di Catanzaro, 88100 Catanzaro, Italy
| | - Wen Zhu
- Department of Chemistry & Biochemistry, Florida State University, Tallahassee, FL 32306, USA
| | - Yuichiro Takagi
- Department of Biochemistry & Molecular Biology, Indiana University School of Medicine, Indianapolis, IN 46202, USA
| | - Nigel G J Richards
- School of Chemistry, Cardiff University, Park Place, Cardiff CF10 3AT, UK
- Foundation for Applied Molecular Evolution, Alachua, FL 32615, USA
| |
Collapse
|
3
|
Nam K, Shao Y, Major DT, Wolf-Watz M. Perspectives on Computational Enzyme Modeling: From Mechanisms to Design and Drug Development. ACS OMEGA 2024; 9:7393-7412. [PMID: 38405524 PMCID: PMC10883025 DOI: 10.1021/acsomega.3c09084] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 11/14/2023] [Revised: 01/15/2024] [Accepted: 01/19/2024] [Indexed: 02/27/2024]
Abstract
Understanding enzyme mechanisms is essential for unraveling the complex molecular machinery of life. In this review, we survey the field of computational enzymology, highlighting key principles governing enzyme mechanisms and discussing ongoing challenges and promising advances. Over the years, computer simulations have become indispensable in the study of enzyme mechanisms, with the integration of experimental and computational exploration now established as a holistic approach to gain deep insights into enzymatic catalysis. Numerous studies have demonstrated the power of computer simulations in characterizing reaction pathways, transition states, substrate selectivity, product distribution, and dynamic conformational changes for various enzymes. Nevertheless, significant challenges remain in investigating the mechanisms of complex multistep reactions, large-scale conformational changes, and allosteric regulation. Beyond mechanistic studies, computational enzyme modeling has emerged as an essential tool for computer-aided enzyme design and the rational discovery of covalent drugs for targeted therapies. Overall, enzyme design/engineering and covalent drug development can greatly benefit from our understanding of the detailed mechanisms of enzymes, such as protein dynamics, entropy contributions, and allostery, as revealed by computational studies. Such a convergence of different research approaches is expected to continue, creating synergies in enzyme research. This review, by outlining the ever-expanding field of enzyme research, aims to provide guidance for future research directions and facilitate new developments in this important and evolving field.
Collapse
Affiliation(s)
- Kwangho Nam
- Department
of Chemistry and Biochemistry, University
of Texas at Arlington, Arlington, Texas 76019, United States
| | - Yihan Shao
- Department
of Chemistry and Biochemistry, University
of Oklahoma, Norman, Oklahoma 73019-5251, United States
| | - Dan T. Major
- Department
of Chemistry and Institute for Nanotechnology & Advanced Materials, Bar-Ilan University, Ramat-Gan 52900, Israel
| | | |
Collapse
|
4
|
Castelli M, Marchetti F, Osuna S, F. Oliveira AS, Mulholland AJ, Serapian SA, Colombo G. Decrypting Allostery in Membrane-Bound K-Ras4B Using Complementary In Silico Approaches Based on Unbiased Molecular Dynamics Simulations. J Am Chem Soc 2024; 146:901-919. [PMID: 38116743 PMCID: PMC10785808 DOI: 10.1021/jacs.3c11396] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/13/2023] [Revised: 12/06/2023] [Accepted: 12/07/2023] [Indexed: 12/21/2023]
Abstract
Protein functions are dynamically regulated by allostery, which enables conformational communication even between faraway residues, and expresses itself in many forms, akin to different "languages": allosteric control pathways predominating in an unperturbed protein are often unintuitively reshaped whenever biochemical perturbations arise (e.g., mutations). To accurately model allostery, unbiased molecular dynamics (MD) simulations require integration with a reliable method able to, e.g., detect incipient allosteric changes or likely perturbation pathways; this is because allostery can operate at longer time scales than those accessible by plain MD. Such methods are typically applied singularly, but we here argue their joint application─as a "multilingual" approach─could work significantly better. We successfully prove this through unbiased MD simulations (∼100 μs) of the widely studied, allosterically active oncotarget K-Ras4B, solvated and embedded in a phospholipid membrane, from which we decrypt allostery using four showcase "languages": Distance Fluctuation analysis and the Shortest Path Map capture allosteric hotspots at equilibrium; Anisotropic Thermal Diffusion and Dynamical Non-Equilibrium MD simulations assess perturbations upon, respectively, either superheating or hydrolyzing the GTP that oncogenically activates K-Ras4B. Chosen "languages" work synergistically, providing an articulate, mutually coherent, experimentally consistent picture of K-Ras4B allostery, whereby distinct traits emerge at equilibrium and upon GTP cleavage. At equilibrium, combined evidence confirms prominent allosteric communication from the membrane-embedded hypervariable region, through a hub comprising helix α5 and sheet β5, and up to the active site, encompassing allosteric "switches" I and II (marginally), and two proposed pockets. Upon GTP cleavage, allosteric perturbations mostly accumulate on the switches and documented interfaces.
Collapse
Affiliation(s)
- Matteo Castelli
- Department
of Chemistry, University of Pavia, viale T. Taramelli 12, 27100 Pavia, Italy
| | - Filippo Marchetti
- Department
of Chemistry, University of Pavia, viale T. Taramelli 12, 27100 Pavia, Italy
- INSTM, via G. Giusti 9, 50121 Florence, Italy
- E4
Computer Engineering, via Martiri delle libertà 66, 42019 Scandiano (RE), Italy
| | - Sílvia Osuna
- Institut
de Química Computacional i Catàlisi (IQCC) and Departament
de Química, Universitat de Girona, Girona, Catalonia E-17071, Spain
- ICREA, Barcelona, Catalonia E-08010, Spain
| | - A. Sofia F. Oliveira
- Centre for
Computational Chemistry, School of Chemistry, University of Bristol, Bristol BS8 1TS, U.K.
| | - Adrian J. Mulholland
- Centre for
Computational Chemistry, School of Chemistry, University of Bristol, Bristol BS8 1TS, U.K.
| | - Stefano A. Serapian
- Department
of Chemistry, University of Pavia, viale T. Taramelli 12, 27100 Pavia, Italy
| | - Giorgio Colombo
- Department
of Chemistry, University of Pavia, viale T. Taramelli 12, 27100 Pavia, Italy
| |
Collapse
|
5
|
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 PMCID: PMC11211065 DOI: 10.1021/acs.jpclett.3c02444] [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] [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
|
6
|
Hu YT, Hong XZ, Li HM, Yang JK, Shen W, Wang YW, Liu YH. Modifying the amino acids in conformational motion pathway of the α-amylase of Geobacillus stearothermophilus improved its activity and stability. Front Microbiol 2023; 14:1261245. [PMID: 38143856 PMCID: PMC10740195 DOI: 10.3389/fmicb.2023.1261245] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/19/2023] [Accepted: 11/21/2023] [Indexed: 12/26/2023] Open
Abstract
Amino acids along the conformational motion pathway of the enzyme molecule correlated to its flexibility and rigidity. To enhance the enzyme activity and thermal stability, the motion pathway of Geobacillus stearothermophilus α-amylase has been identified and molecularly modified by using the neural relational inference model and deep learning tool. The significant differences in substrate specificity, enzymatic kinetics, optimal temperature, and thermal stability were observed among the mutants with modified amino acids along the pathway. Mutants especially the P44E demonstrated enhanced hydrolytic activity and catalytic efficiency (kcat/KM) than the wild-type enzyme to 95.0% and 93.8% respectively, with the optimum temperature increased to 90°C. This mutation from proline to glutamic acid has increased the number and the radius of the bottleneck of the channels, which might facilitate transporting large starch substrates into the enzyme. The mutation could also optimize the hydrogen bonding network of the catalytic center, and diminish the spatial hindering to the substrate entry and exit from the catalytic center.
Collapse
Affiliation(s)
- Yu-Ting Hu
- Pilot Base of Food Microbial Resources Utilization of Hubei Province, College of Life Science and Technology, Wuhan Polytechnic University, Wuhan, China
| | - Xi-Zhi Hong
- Pilot Base of Food Microbial Resources Utilization of Hubei Province, College of Life Science and Technology, Wuhan Polytechnic University, Wuhan, China
| | - Hui-Min Li
- Pilot Base of Food Microbial Resources Utilization of Hubei Province, College of Life Science and Technology, Wuhan Polytechnic University, Wuhan, China
| | - Jiang-Ke Yang
- Pilot Base of Food Microbial Resources Utilization of Hubei Province, College of Life Science and Technology, Wuhan Polytechnic University, Wuhan, China
| | - Wei Shen
- Pilot Base of Food Microbial Resources Utilization of Hubei Province, College of Life Science and Technology, Wuhan Polytechnic University, Wuhan, China
| | - Ya-Wei Wang
- Pilot Base of Food Microbial Resources Utilization of Hubei Province, College of Life Science and Technology, Wuhan Polytechnic University, Wuhan, China
| | - Yi-Han Liu
- Key Laboratory of Industrial Fermentation Microbiology, Ministry of Education, Tianjin Key Laboratory of Industrial Microbiology, The College of Biotechnology, Tianjin University of Science and Technology, Tianjin, China
| |
Collapse
|
7
|
Xiao D, Lu L, Liang B, Xiong Z, Xu X, Chen WH. Identification of 1,3,4-oxadiazolyl-containing β-carboline derivatives as novel α-glucosidase inhibitors with antidiabetic activity. Eur J Med Chem 2023; 261:115795. [PMID: 37688939 DOI: 10.1016/j.ejmech.2023.115795] [Citation(s) in RCA: 14] [Impact Index Per Article: 14.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/30/2023] [Revised: 08/30/2023] [Accepted: 09/04/2023] [Indexed: 09/11/2023]
Abstract
In this study, we designed and synthesized a novel class of 1,3,4-oxadiazolyl-containing β-carboline derivatives, i.e., compounds f1∼f35 as potential α-glucosidase inhibitors. All the synthesized compounds possessed outstanding α-glucosidase inhibitory activity with the IC50 values in the range of 3.07-15.49 μM, representing that they are 36∼183-fold more active than a positive control, acarbose (IC50 = 564.28 μM). Among them, compound f26 exhibited the highest α-glucosidase inhibitory activity (IC50 = 3.07 μM) and was demonstrated to function as a reversible and noncompetitive inhibitor. Mechanistic studies by means of 3D fluorescence spectra, CD spectra and molecular docking suggested that complexation of compound f26 with α-glucosidase through hydrogen bonds and hydrophobic interactions, led to changes in the conformation and secondary strictures of α-glucosidase and further the inhibition of the enzymatic activity. In vivo results showed that oral administration of compound f26 (50 mg/kg/day) could obviously reduce the levels of fasting blood glucose and improve glucose tolerance and dyslipidemia in diabetic mice. The present findings suggest that compound f26 is exploitable as a potential lead compound for the development of new α-glucosidase inhibitors with antidiabetic activity.
Collapse
Affiliation(s)
- Di Xiao
- School of Biotechnology and Health Sciences, Wuyi University, Jiangmen, 529020, China
| | - Li Lu
- School of Biotechnology and Health Sciences, Wuyi University, Jiangmen, 529020, China
| | - Bingwen Liang
- School of Biotechnology and Health Sciences, Wuyi University, Jiangmen, 529020, China
| | - Zhuang Xiong
- School of Biotechnology and Health Sciences, Wuyi University, Jiangmen, 529020, China
| | - Xuetao Xu
- School of Biotechnology and Health Sciences, Wuyi University, Jiangmen, 529020, China.
| | - Wen-Hua Chen
- School of Biotechnology and Health Sciences, Wuyi University, Jiangmen, 529020, China.
| |
Collapse
|
8
|
Abstract
ConspectusThe quantum chemical cluster approach has been used for modeling enzyme active sites and reaction mechanisms for more than two decades. In this methodology, a relatively small part of the enzyme around the active site is selected as a model, and quantum chemical methods, typically density functional theory, are used to calculate energies and other properties. The surrounding enzyme is modeled using implicit solvation and atom fixing techniques. Over the years, a large number of enzyme mechanisms have been solved using this method. The models have gradually become larger as a result of the faster computers, and new kinds of questions have been addressed. In this Account, we review how the cluster approach can be utilized in the field of biocatalysis. Examples from our recent work are chosen to illustrate various aspects of the methodology. The use of the cluster model to explore substrate binding is discussed first. It is emphasized that a comprehensive search is necessary in order to identify the lowest-energy binding mode(s). It is also argued that the best binding mode might not be the productive one, and the full reactions for a number of enzyme-substrate complexes have therefore to be considered to find the lowest-energy reaction pathway. Next, examples are given of how the cluster approach can help in the elucidation of detailed reaction mechanisms of biocatalytically interesting enzymes, and how this knowledge can be exploited to develop enzymes with new functions or to understand the reasons for lack of activity toward non-natural substrates. The enzymes discussed in this context are phenolic acid decarboxylase and metal-dependent decarboxylases from the amidohydrolase superfamily. Next, the application of the cluster approach in the investigation of enzymatic enantioselectivity is discussed. The reaction of strictosidine synthase is selected as a case study, where the cluster calculations could reproduce and rationalize the selectivities of both the natural and non-natural substrates. Finally, we discuss how the cluster approach can be used to guide the rational design of enzyme variants with improved activity and selectivity. Acyl transferase from Mycobacterium smegmatis serves as an instructive example here, for which the calculations could pinpoint the factors controlling the reaction specificity and enantioselectivity. The cases discussed in this Account highlight thus the value of the cluster approach as a tool in biocatalysis. It complements experiments and other computational techniques in this field and provides insights that can be used to understand existing enzymes and to develop new variants with tailored properties.
Collapse
Affiliation(s)
- Xiang Sheng
- Tianjin Institute of Industrial Biotechnology, Chinese Academy of Sciences, Tianjin 300308, PR China
- National Center of Technology Innovation for Synthetic Biology, Tianjin 300308, PR China
| | - Fahmi Himo
- Department of Organic Chemistry, Arrhenius Laboratory, Stockholm University, SE-10691 Stockholm, Sweden
| |
Collapse
|
9
|
Jiang Y, Yan B, Chen Y, Juarez RJ, Yang ZJ. Molecular Dynamics-Derived Descriptor Informs the Impact of Mutation on the Catalytic Turnover Number in Lactonase Across Substrates. J Phys Chem B 2022; 126:2486-2495. [PMID: 35324218 DOI: 10.1021/acs.jpcb.2c00142] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Molecular dynamics simulations have been extensively employed to reveal the roles of protein dynamics in mediating enzyme catalysis. However, simulation-derived predictive descriptors that inform the impacts of mutations on catalytic turnover numbers remain largely unexplored. In this work, we report the identification of molecular modeling-derived descriptors to predict mutation effect on the turnover number of lactonase SsoPox with both native and non-native substrates. The study consists of 10 enzyme-substrate complexes resulting from a combination of five enzyme variants with two substrates. For each complex, we derived 15 descriptors from molecular dynamics simulations and applied principal component analysis to rank the predictive capability of the descriptors. A top-ranked descriptor was identified, which is the solvent-accessible surface area (SASA) ratio of the substrate to the active site pocket. A uniform volcano-shaped plot was observed in the distribution of experimental activation free energy against the SASA ratio. To achieve efficient lactonase hydrolysis, a non-native substrate-bound enzyme variant needs to involve a similar range of the SASA ratio to the native substrate-bound wild-type enzyme. The descriptor reflects how well the enzyme active site pocket accommodates a substrate for reaction, which has the potential of guiding optimization of enzyme reaction turnover for non-native chemical transformations.
Collapse
Affiliation(s)
- Yaoyukun Jiang
- Department of Chemistry, Vanderbilt University, Nashville, Tennessee 37235, United States
| | - Bailu Yan
- Department of Biostatistics, Vanderbilt University, Nashville, Tennessee 37235, United States
| | - Yu Chen
- School of Chemistry and Chemical Engineering, Nanjing University, Nanjing, Jiangsu 210023, China
| | - Reecan J Juarez
- Chemical and Physical Biology Program, 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
| |
Collapse
|
10
|
Qu G, Bi Y, Liu B, Li J, Han X, Liu W, Jiang Y, Qin Z, Sun Z. Unlocking the Stereoselectivity and Substrate Acceptance of Enzymes: Proline‐Induced Loop Engineering Test. Angew Chem Int Ed Engl 2022. [DOI: 10.1002/ange.202110793] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022]
Affiliation(s)
- Ge Qu
- Tianjin Institute of Industrial Biotechnology Chinese Academy of Sciences Tianjin 300308 China
- National Technology Innovation Center of Synthetic Biology Tianjin 300308 China
| | - Yuexin Bi
- Tianjin Institute of Industrial Biotechnology Chinese Academy of Sciences Tianjin 300308 China
- University of Science and Technology of China Hefei 230027 China
| | - Beibei Liu
- Tianjin Institute of Industrial Biotechnology Chinese Academy of Sciences Tianjin 300308 China
| | - Junkuan Li
- Tianjin Institute of Industrial Biotechnology Chinese Academy of Sciences Tianjin 300308 China
- Department of Chemistry School of Science Tianjin University Tianjin 300072 China
| | - Xu Han
- Tianjin Institute of Industrial Biotechnology Chinese Academy of Sciences Tianjin 300308 China
- National Technology Innovation Center of Synthetic Biology Tianjin 300308 China
| | - Weidong Liu
- Tianjin Institute of Industrial Biotechnology Chinese Academy of Sciences Tianjin 300308 China
- National Technology Innovation Center of Synthetic Biology Tianjin 300308 China
| | - Yingying Jiang
- Tianjin Institute of Industrial Biotechnology Chinese Academy of Sciences Tianjin 300308 China
- University of Chinese Academy of Sciences Beijing 100049 China
| | - Zongmin Qin
- Tianjin Institute of Industrial Biotechnology Chinese Academy of Sciences Tianjin 300308 China
- University of Chinese Academy of Sciences Beijing 100049 China
| | - Zhoutong Sun
- Tianjin Institute of Industrial Biotechnology Chinese Academy of Sciences Tianjin 300308 China
- National Technology Innovation Center of Synthetic Biology Tianjin 300308 China
| |
Collapse
|
11
|
Green biomanufacturing promoted by automatic retrobiosynthesis planning and computational enzyme design. Chin J Chem Eng 2022. [DOI: 10.1016/j.cjche.2021.08.017] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
|
12
|
Qu G, Bi Y, Liu B, Li J, Han X, Liu W, Jiang Y, Qin Z, Sun Z. Unlocking the Stereoselectivity and Substrate Acceptance of Enzymes: Proline-Induced Loop Engineering Test. Angew Chem Int Ed Engl 2021; 61:e202110793. [PMID: 34658118 DOI: 10.1002/anie.202110793] [Citation(s) in RCA: 17] [Impact Index Per Article: 5.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/11/2021] [Revised: 09/25/2021] [Indexed: 12/12/2022]
Abstract
Protein stability and evolvability influence each other. Although protein dynamics play essential roles in various catalytically important properties, their high flexibility and diversity makes it difficult to incorporate such properties into rational engineering. Therefore, how to unlock the potential evolvability in a user-friendly rational design process remains a challenge. In this endeavor, we describe a method for engineering an enantioselective alcohol dehydrogenase. It enables synthetically important substrate acceptance for 4-chlorophenyl pyridine-2-yl ketone, and perfect stereocontrol of both (S)- and (R)-configured products. Thermodynamic analysis unveiled the subtle interaction between enzyme stability and evolvability, while computational studies provided insights into the origin of selectivity and substrate recognition. Preparative-scale synthesis of the (S)-product (73 % yield; >99 % ee) was performed on a gram-scale. This proof-of-principle study demonstrates that interfaced proline residues can be rationally engineered to unlock evolvability and thus provide access to new biocatalysts with highly improved catalytic performance.
Collapse
Affiliation(s)
- Ge Qu
- Tianjin Institute of Industrial Biotechnology, Chinese Academy of Sciences, Tianjin, 300308, China.,National Technology Innovation Center of Synthetic Biology, Tianjin, 300308, China
| | - Yuexin Bi
- Tianjin Institute of Industrial Biotechnology, Chinese Academy of Sciences, Tianjin, 300308, China.,University of Science and Technology of China, Hefei, 230027, China
| | - Beibei Liu
- Tianjin Institute of Industrial Biotechnology, Chinese Academy of Sciences, Tianjin, 300308, China
| | - Junkuan Li
- Tianjin Institute of Industrial Biotechnology, Chinese Academy of Sciences, Tianjin, 300308, China.,Department of Chemistry, School of Science, Tianjin University, Tianjin, 300072, China
| | - Xu Han
- Tianjin Institute of Industrial Biotechnology, Chinese Academy of Sciences, Tianjin, 300308, China.,National Technology Innovation Center of Synthetic Biology, Tianjin, 300308, China
| | - Weidong Liu
- Tianjin Institute of Industrial Biotechnology, Chinese Academy of Sciences, Tianjin, 300308, China.,National Technology Innovation Center of Synthetic Biology, Tianjin, 300308, China
| | - Yingying Jiang
- Tianjin Institute of Industrial Biotechnology, Chinese Academy of Sciences, Tianjin, 300308, China.,University of Chinese Academy of Sciences, Beijing, 100049, China
| | - Zongmin Qin
- Tianjin Institute of Industrial Biotechnology, Chinese Academy of Sciences, Tianjin, 300308, China.,University of Chinese Academy of Sciences, Beijing, 100049, China
| | - Zhoutong Sun
- Tianjin Institute of Industrial Biotechnology, Chinese Academy of Sciences, Tianjin, 300308, China.,National Technology Innovation Center of Synthetic Biology, Tianjin, 300308, China
| |
Collapse
|
13
|
Yan B, Ran X, Jiang Y, Torrence SK, Yuan L, Shao Q, Yang ZJ. Rate-Perturbing Single Amino Acid Mutation for Hydrolases: A Statistical Profiling. J Phys Chem B 2021; 125:10682-10691. [PMID: 34524819 DOI: 10.1021/acs.jpcb.1c05901] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Hydrolases are a critical component for modern chemical, pharmaceutical, and environmental sciences. Identifying mutations that enhance catalytic efficiency presents a roadblock to design and to discover new hydrolases for broad academic and industrial uses. Here, we report the statistical profiling for rate-perturbing mutant hydrolases with a single amino acid substitution. We constructed an integrated structure-kinetics database for hydrolases, IntEnzyDB, which contains 3907 kcats, 4175 KMs, and 2715 Protein Data Bank IDs. IntEnzyDB adopts a relational architecture with a flattened data structure, enabling facile and efficient access to clean and tabulated data for machine learning uses. We conducted statistical analyses on how single amino acids mutations influence the turnover number (i.e., kcat) and efficiency (i.e., kcat/KM), with a particular emphasis on profiling the features for rate-enhancing mutations. The results show that mutation to bulky nonpolar residues with a hydrocarbon chain involves a higher likelihood for rate acceleration than to other types of residues. Linear regression models reveal geometric descriptors of substrate and mutation residues that mediate rate-perturbing outcomes for hydrolases with bulky nonpolar mutations. On the basis of the analyses of the structure-kinetics relationship, we observe that the propensity for rate enhancement is independent of protein sizes. In addition, we observe that distal mutations (i.e., >10 Å from the active site) in hydrolases are significantly more prone to induce efficiency neutrality and avoid efficiency deletion but involve similar propensity for rate enhancement. The studies reveal the statistical features for identifying rate-enhancing mutations in hydrolases, which will potentially guide hydrolase discovery in biocatalysis.
Collapse
Affiliation(s)
- Bailu Yan
- Department of Chemistry, Vanderbilt University, Nashville, Tennessee 37235, United States.,Department of Biostatistics, Vanderbilt University, Nashville, Tennessee 37203, United States
| | - Xinchun Ran
- Department of Chemistry, Vanderbilt University, Nashville, Tennessee 37235, United States
| | - Yaoyukun Jiang
- Department of Chemistry, Vanderbilt University, Nashville, Tennessee 37235, United States
| | - Sarah K Torrence
- Data Science Institute, Vanderbilt University, Nashville, Tennessee 37235, United States
| | - Li Yuan
- Data Science Institute, Vanderbilt University, Nashville, Tennessee 37235, United States
| | - Qianzhen Shao
- 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
| |
Collapse
|
14
|
Aledavood E, Forte A, Estarellas C, Javier Luque F. Structural basis of the selective activation of enzyme isoforms: Allosteric response to activators of β1- and β2-containing AMPK complexes. Comput Struct Biotechnol J 2021; 19:3394-3406. [PMID: 34194666 PMCID: PMC8217686 DOI: 10.1016/j.csbj.2021.05.056] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/03/2021] [Revised: 05/30/2021] [Accepted: 05/30/2021] [Indexed: 12/21/2022] Open
Abstract
AMP-activated protein kinase (AMPK) is a key energy sensor regulating the cell metabolism in response to energy supply and demand. The evolutionary adaptation of AMPK to different tissues is accomplished through the expression of distinct isoforms that can form up to 12 complexes, which exhibit notable differences in the sensitivity to allosteric activators. To shed light into the molecular determinants of the allosteric regulation of this energy sensor, we have examined the structural and dynamical properties of β1- and β2-containing AMPK complexes formed with small molecule activators A-769662 and SC4, and dissected the mechanical response leading to active-like enzyme conformations through the analysis of interaction networks between structural domains. The results reveal the mechanical sensitivity of the α2β1 complex, in contrast with a larger resilience of the α2β2 species, especially regarding modulation by A-769662. Furthermore, binding of activators to α2β1 consistently promotes the pre-organization of the ATP-binding site, favoring the adoption of activated states of the enzyme. These findings are discussed in light of the changes in the residue content of β-subunit isoforms, particularly regarding the β1Asn111 → β2Asp111 substitution as a key factor in modulating the mechanical sensitivity of β1- and β2-containing AMPK complexes. Our studies pave the way for the design of activators tailored for improving the therapeutic treatment of tissue-specific metabolic disorders.
Collapse
Affiliation(s)
| | - Alessia Forte
- Department of Nutrition, Food Science and Gastronomy, Faculty of Pharmacy and Food Sciences, Institute of Biomedicine (IBUB) and Institute of Theoretical and Computational Chemistry (IQTCUB), University of Barcelona, Av. Prat de la Riba 171, Santa Coloma de Gramenet 08921, Spain
| | | | | |
Collapse
|
15
|
Sequeiros-Borja CE, Surpeta B, Brezovsky J. Recent advances in user-friendly computational tools to engineer protein function. Brief Bioinform 2021; 22:bbaa150. [PMID: 32743637 PMCID: PMC8138880 DOI: 10.1093/bib/bbaa150] [Citation(s) in RCA: 29] [Impact Index Per Article: 9.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/24/2020] [Revised: 06/03/2020] [Accepted: 06/16/2020] [Indexed: 12/14/2022] Open
Abstract
Progress in technology and algorithms throughout the past decade has transformed the field of protein design and engineering. Computational approaches have become well-engrained in the processes of tailoring proteins for various biotechnological applications. Many tools and methods are developed and upgraded each year to satisfy the increasing demands and challenges of protein engineering. To help protein engineers and bioinformaticians navigate this emerging wave of dedicated software, we have critically evaluated recent additions to the toolbox regarding their application for semi-rational and rational protein engineering. These newly developed tools identify and prioritize hotspots and analyze the effects of mutations for a variety of properties, comprising ligand binding, protein-protein and protein-nucleic acid interactions, and electrostatic potential. We also discuss notable progress to target elusive protein dynamics and associated properties like ligand-transport processes and allosteric communication. Finally, we discuss several challenges these tools face and provide our perspectives on the further development of readily applicable methods to guide protein engineering efforts.
Collapse
Affiliation(s)
- Carlos Eduardo Sequeiros-Borja
- Laboratory of Biomolecular Interactions and Transport, Department of Gene Expression, Institute of Molecular Biology and Biotechnology, Faculty of Biology, Adam Mickiewicz University and the International Institute of Molecular and Cell Biology in Warsaw, Warsaw, Poland
| | - Bartłomiej Surpeta
- Laboratory of Biomolecular Interactions and Transport, Department of Gene Expression, Institute of Molecular Biology and Biotechnology, Faculty of Biology, Adam Mickiewicz University and the International Institute of Molecular and Cell Biology in Warsaw, Warsaw, Poland
| | - Jan Brezovsky
- Laboratory of Biomolecular Interactions and Transport, Department of Gene Expression, Institute of Molecular Biology and Biotechnology, Faculty of Biology, Adam Mickiewicz University and the International Institute of Molecular and Cell Biology in Warsaw
| |
Collapse
|
16
|
Pervasive cooperative mutational effects on multiple catalytic enzyme traits emerge via long-range conformational dynamics. Nat Commun 2021; 12:1621. [PMID: 33712579 PMCID: PMC7955134 DOI: 10.1038/s41467-021-21833-w] [Citation(s) in RCA: 60] [Impact Index Per Article: 20.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/25/2020] [Accepted: 01/29/2021] [Indexed: 12/11/2022] Open
Abstract
Multidimensional fitness landscapes provide insights into the molecular basis of laboratory and natural evolution. To date, such efforts usually focus on limited protein families and a single enzyme trait, with little concern about the relationship between protein epistasis and conformational dynamics. Here, we report a multiparametric fitness landscape for a cytochrome P450 monooxygenase that was engineered for the regio- and stereoselective hydroxylation of a steroid. We develop a computational program to automatically quantify non-additive effects among all possible mutational pathways, finding pervasive cooperative signs and magnitude epistasis on multiple catalytic traits. By using quantum mechanics and molecular dynamics simulations, we show that these effects are modulated by long-range interactions in loops, helices and β-strands that gate the substrate access channel allowing for optimal catalysis. Our work highlights the importance of conformational dynamics on epistasis in an enzyme involved in secondary metabolism and offers insights for engineering P450s. Connecting conformational dynamics and epistasis has so far been limited to a few proteins and a single fitness trait. Here, the authors provide evidence of positive epistasis on multiple catalytic traits in the evolution and dynamics of engineered cytochrome P450 monooxygenase, offering insights for in silico protein design.
Collapse
|
17
|
Karamitros CS, Murray K, Sugiyama Y, Kumada Y, Johnson KA, Georgiou G, D'Arcy S, Stone EM. Conformational Dynamics Contribute to Substrate Selectivity and Catalysis in Human Kynureninase. ACS Chem Biol 2020; 15:3159-3166. [PMID: 33275413 PMCID: PMC11104311 DOI: 10.1021/acschembio.0c00676] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/04/2023]
Abstract
Kynureninases (KYNases) are enzymes that play a key role in tryptophan catabolism through the degradation of intermediate kynurenine and 3'-hydroxy-kynurenine metabolites (KYN and OH-KYN, respectively). Bacterial KYNases exhibit high catalytic efficiency toward KYN and moderate activity toward OH-KYN, whereas animal KYNases are highly selective for OH-KYN, exhibiting only minimal activity toward the smaller KYN substrate. These differences reflect divergent pathways for KYN and OH-KYN utilization in the respective kingdoms. We examined the Homo sapiens and Pseudomonas fluorescens KYNases (HsKYNase and PfKYNase respectively) using pre-steady-state and hydrogen-deuterium exchange mass spectrometry (HDX-MS) methodologies. We discovered that the activity of HsKYNase critically depends on formation of hydrogen bonds with the hydroxyl group of OH-KYN to stabilize the entire active site and allow productive substrate turnover. With the preferred OH-KYN substrate, stabilization is observed at the substrate-binding site and the region surrounding the PLP cofactor. With the nonpreferred KYN substrate, less stabilization occurs, revealing a direct correlation with activity. This correlation holds true for PfKYNases; however there is only a modest stabilization at the substrate-binding site, suggesting that substrate discrimination is simply achieved by steric hindrance. We speculate that eukaryotic KYNases use dynamic mobility as a mechanism of substrate specificity to commit OH-KYN to nicotinamide synthesis and avoid futile hydrolysis of KYN. These findings have important ramifications for the engineering of HsKynase with high KYN activity as required for clinical applications in cancer immunotherapy. Our study shows how homologous enzymes with conserved active sites can use dynamics to discriminate between two highly similar substrates.
Collapse
Affiliation(s)
- Christos S Karamitros
- Department of Chemical Engineering, The University of Texas at Austin, Austin, Texas 78712, United States
| | - Kyle Murray
- Department of Chemistry and Biochemistry, The University of Texas at Dallas, Richardson, Texas 75080, United States
| | - Yusuke Sugiyama
- Department of Molecular Chemistry and Engineering, Graduate School of Science and Technology, Kyoto Institute of Technology, Kyoto, Japan
| | - Yoichi Kumada
- Department of Molecular Chemistry and Engineering, Graduate School of Science and Technology, Kyoto Institute of Technology, Kyoto, Japan
| | - Kenneth A Johnson
- Department of Molecular Biosciences, The University of Texas at Austin, Austin, Texas 78712, United States
| | - George Georgiou
- Department of Chemical Engineering, The University of Texas at Austin, Austin, Texas 78712, United States
- Department of Molecular Biosciences, The University of Texas at Austin, Austin, Texas 78712, United States
- Institute for Cellular and Molecular Biology, The University of Texas at Austin, Austin, Texas 78712, United States
- Department of Oncology, University of Texas Dell Medical School, LiveSTRONG Cancer Institutes, Austin, Texas 78712, United States
| | - Sheena D'Arcy
- Department of Chemistry and Biochemistry, The University of Texas at Dallas, Richardson, Texas 75080, United States
| | - Everett M Stone
- Department of Molecular Biosciences, The University of Texas at Austin, Austin, Texas 78712, United States
- Department of Oncology, University of Texas Dell Medical School, LiveSTRONG Cancer Institutes, Austin, Texas 78712, United States
| |
Collapse
|
18
|
|
19
|
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: 57] [Impact Index Per Article: 14.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
|
20
|
De Raffele D, Martí S, Moliner V. Understanding the Directed Evolution of De Novo Retro-Aldolases from QM/MM Studies. ACS Catal 2020. [DOI: 10.1021/acscatal.0c01165] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022]
Affiliation(s)
- Daria De Raffele
- Departament de Química Física i Analítica; Universitat Jaume I., 12071 Castellón, Spain
| | - Sergio Martí
- 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
|
21
|
Surpeta B, Sequeiros-Borja CE, Brezovsky J. Dynamics, a Powerful Component of Current and Future in Silico Approaches for Protein Design and Engineering. Int J Mol Sci 2020; 21:E2713. [PMID: 32295283 PMCID: PMC7215530 DOI: 10.3390/ijms21082713] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/18/2020] [Revised: 04/10/2020] [Accepted: 04/12/2020] [Indexed: 12/13/2022] Open
Abstract
Computational prediction has become an indispensable aid in the processes of engineering and designing proteins for various biotechnological applications. With the tremendous progress in more powerful computer hardware and more efficient algorithms, some of in silico tools and methods have started to apply the more realistic description of proteins as their conformational ensembles, making protein dynamics an integral part of their prediction workflows. To help protein engineers to harness benefits of considering dynamics in their designs, we surveyed new tools developed for analyses of conformational ensembles in order to select engineering hotspots and design mutations. Next, we discussed the collective evolution towards more flexible protein design methods, including ensemble-based approaches, knowledge-assisted methods, and provable algorithms. Finally, we highlighted apparent challenges that current approaches are facing and provided our perspectives on their further development.
Collapse
Affiliation(s)
- Bartłomiej Surpeta
- Laboratory of Biomolecular Interactions and Transport, Department of Gene Expression, Institute of Molecular Biology and Biotechnology, Faculty of Biology, Adam Mickiewicz University, Uniwersytetu Poznanskiego 6, 61-614 Poznan, Poland; (B.S.); (C.E.S.-B.)
- International Institute of Molecular and Cell Biology in Warsaw, Ks Trojdena 4, 02-109 Warsaw, Poland
| | - Carlos Eduardo Sequeiros-Borja
- Laboratory of Biomolecular Interactions and Transport, Department of Gene Expression, Institute of Molecular Biology and Biotechnology, Faculty of Biology, Adam Mickiewicz University, Uniwersytetu Poznanskiego 6, 61-614 Poznan, Poland; (B.S.); (C.E.S.-B.)
- International Institute of Molecular and Cell Biology in Warsaw, Ks Trojdena 4, 02-109 Warsaw, Poland
| | - Jan Brezovsky
- Laboratory of Biomolecular Interactions and Transport, Department of Gene Expression, Institute of Molecular Biology and Biotechnology, Faculty of Biology, Adam Mickiewicz University, Uniwersytetu Poznanskiego 6, 61-614 Poznan, Poland; (B.S.); (C.E.S.-B.)
- International Institute of Molecular and Cell Biology in Warsaw, Ks Trojdena 4, 02-109 Warsaw, Poland
| |
Collapse
|
22
|
Levina EO, Khrenova MG, Astakhov AA, Tsirelson VG. Revealing electronic features governing hydrolysis of cephalosporins in the active site of the L1 metallo-β-lactamase. RSC Adv 2020; 10:8664-8676. [PMID: 35496524 PMCID: PMC9050041 DOI: 10.1039/c9ra10649a] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/17/2019] [Accepted: 02/14/2020] [Indexed: 11/21/2022] Open
Abstract
The QM/MM simulations followed by electron density feature analysis are carried out to deepen the understanding of the reaction mechanism of cephalosporin hydrolysis in the active site of the L1 metallo-β-lactamase. The differences in reactivity of ten similar cephalosporin compounds are explained by using an extended set of bonding descriptors. The limiting step of the reaction is characterized by the proton transfer to the nitrogen atom of the cephalosporin thiazine ring accompanied with formation of the C4
Created by potrace 1.16, written by Peter Selinger 2001-2019
]]>
C3 double bond in its N–C4–C3 fragment. The temporary N⋯H–Ow hydrogen bond, which is formed in the transition state of the limiting step of the reaction was recognized as a key atomic interaction governing the reactivity of various cephalosporins. Non-local real-space bonding descriptors show that different extent of localization of electron lone pair at N atom in the transition state affect the reactivity of compounds: smaller electron localization is typical for the less reactive species. In particular, the Fermi hole analysis shows how exchange electron correlation in the N⋯H–Ow fragment control electron lone pair localization. Delocalization tensor, linear response kernel and source function indicate that features of electron delocalization in the N–C4–C3 fragment of cephalosporins in the transition state complexes determine the differences in C4–C3 bond for substrates with high and low rate constants. The C4–C3 bond of the N–C4–C3 fragment at the transition state is similar to that of the preceding intermediate for the less reactive species and resembles the features of the enzyme–product complex for more reactive compounds. The power and limitations of the descriptors applied for solving the problem are discussed and the generality of approach is stressed. Combination of QM/MM and modern bonding descriptors explains different reactivity of cephalosporins in the active site of the L1 metallo-β-lactamase.![]()
Collapse
Affiliation(s)
- Elena O. Levina
- Federal Research Centre “Fundamentals of Biotechnology” of the Russian Academy of Sciences
- Moscow
- Russia
- Moscow Institute of Physics and Technology
- Dolgoprudny
| | - Maria G. Khrenova
- Federal Research Centre “Fundamentals of Biotechnology” of the Russian Academy of Sciences
- Moscow
- Russia
- Lomonosov Moscow State University
- Moscow
| | - Andrey A. Astakhov
- Federal Research Centre “Fundamentals of Biotechnology” of the Russian Academy of Sciences
- Moscow
- Russia
- Joint Institute for Nuclear Research
- Dubna
| | - Vladimir G. Tsirelson
- Federal Research Centre “Fundamentals of Biotechnology” of the Russian Academy of Sciences
- Moscow
- Russia
- Mendeleev University of Chemical Technology of Russia
- Moscow
| |
Collapse
|
23
|
Santos YLDL, Chew-Fajardo YL, Brault G, Doucet N. Dissecting the evolvability landscape of the CalB active site toward aromatic substrates. Sci Rep 2019; 9:15588. [PMID: 31666622 PMCID: PMC6821916 DOI: 10.1038/s41598-019-51940-0] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/18/2019] [Accepted: 10/07/2019] [Indexed: 01/17/2023] Open
Abstract
A key event in the directed evolution of enzymes is the systematic use of mutagenesis and selection, a process that can give rise to mutant libraries containing millions of protein variants. To this day, the functional analysis and identification of active variants among such high numbers of mutational possibilities is not a trivial task. Here, we describe a combinatorial semi-rational approach to partly overcome this challenge and help design smaller and smarter mutant libraries. By adapting a liquid medium transesterification assay in organic solvent conditions with a combination of virtual docking, iterative saturation mutagenesis, and residue interaction network (RIN) analysis, we engineered lipase B from P. antarctica (CalB) to improve enzyme recognition and activity against the bulky aromatic substrates and flavoring agents methyl cinnamate and methyl salicylate. Substrate-imprinted docking was used to target active-site positions involved in enzyme-substrate and enzyme-product complexes, in addition to identifying 'hot spots' most likely to yield active variants. This iterative semi-rational design strategy allowed selection of CalB variants exhibiting increased activity in just two rounds of site-saturation mutagenesis. Beneficial replacements were observed by screening only 0.308% of the theoretical library size, illustrating how semi-rational approaches with targeted diversity can quickly facilitate the discovery of improved activity variants relevant to a number of biotechnological applications.
Collapse
Affiliation(s)
- Yossef López de Los Santos
- Centre Armand-Frappier Santé Biotechnologie, Institut National de la Recherche Scientifique (INRS), Université du Québec, 531 Boulevard des Prairies, Laval, QC, H7V 1B7, Canada
| | - Ying Lian Chew-Fajardo
- Centre Armand-Frappier Santé Biotechnologie, Institut National de la Recherche Scientifique (INRS), Université du Québec, 531 Boulevard des Prairies, Laval, QC, H7V 1B7, Canada
| | - Guillaume Brault
- Centre Armand-Frappier Santé Biotechnologie, Institut National de la Recherche Scientifique (INRS), Université du Québec, 531 Boulevard des Prairies, Laval, QC, H7V 1B7, Canada
| | - Nicolas Doucet
- Centre Armand-Frappier Santé Biotechnologie, Institut National de la Recherche Scientifique (INRS), Université du Québec, 531 Boulevard des Prairies, Laval, QC, H7V 1B7, Canada.
- PROTEO, the Québec Network for Research on Protein Function, Engineering, and Applications, 1045 Avenue de la Médecine, Université Laval, Quebec City, QC, G1V 0A6, Canada.
| |
Collapse
|
24
|
Schafer JW, Zoi I, Antoniou D, Schwartz SD. Optimization of the Turnover in Artificial Enzymes via Directed Evolution Results in the Coupling of Protein Dynamics to Chemistry. J Am Chem Soc 2019; 141:10431-10439. [PMID: 31199129 DOI: 10.1021/jacs.9b04515] [Citation(s) in RCA: 24] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The design of artificial enzymes is an emerging field of research. Although progress has been made, the catalytic proficiency of many designed enzymes is low compared to natural enzymes. Nevertheless, recently Hilvert et al. ( Nat. Chem. 2017, 9, 50-56) created a series of five artificial retro-aldolase enzymes via directed evolution, with the final variant exhibiting a rate comparable to the naturally occurring enzyme fructose 1,6 bisphosphate aldolase. We present a study of this system in atomistic detail that elucidates the effects of mutational changes on the chemical step. Transition path sampling is used to create ensembles of reactive trajectories, and committor analysis is used to identify the stochastic separatrix of each ensemble. The application of committor distribution analysis to constrained trajectories allows the identification of changes in important protein motions coupled to reaction across the generated series of the artificial retro-aldolases. We observed two different reaction mechanisms and analyzed the role of the residues participating in the reaction coordinate of each enzyme. However, only in the most evolved variant we identified a fast motion that promotes catalysis, suggesting that this rate promoting vibration was introduced during directed evolution. This study provides further evidence that protein dynamics must be taken into account in designing efficient artificial enzymes.
Collapse
Affiliation(s)
- Joseph W Schafer
- Department of Chemistry and Biochemistry , University of Arizona , Tucson , Arizona 85721 , United States
| | - Ioanna Zoi
- Department of Chemistry and Biochemistry , University of Arizona , Tucson , Arizona 85721 , United States
| | - Dimitri Antoniou
- Department of Chemistry and Biochemistry , University of Arizona , Tucson , Arizona 85721 , United States
| | - Steven D Schwartz
- Department of Chemistry and Biochemistry , University of Arizona , Tucson , Arizona 85721 , United States
| |
Collapse
|
25
|
Curado‐Carballada C, Feixas F, Osuna S. Molecular Dynamics Simulations on
Aspergillus niger
Monoamine Oxidase: Conformational Dynamics and Inter‐monomer Communication Essential for Its Efficient Catalysis. Adv Synth Catal 2019. [DOI: 10.1002/adsc.201900158] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Affiliation(s)
- Christian Curado‐Carballada
- CompBioLab group, Institut de Química Computacional i Catàlisi (IQCC)Departament de Química Carrer Maria Aurèlia Capmany 69 17003 Girona, Catalonia Spain
| | - Ferran Feixas
- CompBioLab group, Institut de Química Computacional i Catàlisi (IQCC)Departament de Química Carrer Maria Aurèlia Capmany 69 17003 Girona, Catalonia Spain
| | - Sílvia Osuna
- CompBioLab group, Institut de Química Computacional i Catàlisi (IQCC)Departament de Química Carrer Maria Aurèlia Capmany 69 17003 Girona, Catalonia Spain
- ICREAPg. Lluís Companys 23 08010 Barcelona Spain
| |
Collapse
|
26
|
Machine learning-assisted directed protein evolution with combinatorial libraries. Proc Natl Acad Sci U S A 2019; 116:8852-8858. [PMID: 30979809 DOI: 10.1073/pnas.1901979116] [Citation(s) in RCA: 279] [Impact Index Per Article: 55.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
To reduce experimental effort associated with directed protein evolution and to explore the sequence space encoded by mutating multiple positions simultaneously, we incorporate machine learning into the directed evolution workflow. Combinatorial sequence space can be quite expensive to sample experimentally, but machine-learning models trained on tested variants provide a fast method for testing sequence space computationally. We validated this approach on a large published empirical fitness landscape for human GB1 binding protein, demonstrating that machine learning-guided directed evolution finds variants with higher fitness than those found by other directed evolution approaches. We then provide an example application in evolving an enzyme to produce each of the two possible product enantiomers (i.e., stereodivergence) of a new-to-nature carbene Si-H insertion reaction. The approach predicted libraries enriched in functional enzymes and fixed seven mutations in two rounds of evolution to identify variants for selective catalysis with 93% and 79% ee (enantiomeric excess). By greatly increasing throughput with in silico modeling, machine learning enhances the quality and diversity of sequence solutions for a protein engineering problem.
Collapse
|
27
|
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]
|
28
|
Voth GA, Yeager M. Editorial overview: COSB biophysical and computational methods. Curr Opin Struct Biol 2018; 52:vi-vii. [PMID: 30554602 DOI: 10.1016/j.sbi.2018.11.011] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Gregory A Voth
- The University of Chicago, Department of Chemistry, 5735 S. Ellis Ave, Chicago, IL 60637, USA
| | - Mark Yeager
- University of Virginia School of Medicine, Department of Molecular Physiology and Biological Physics, Department of Medicine, Division of Cardiology, Sheridan G. Snyder Translational Research Building, 480 Ray C. Hunt Drive, Charlottesville, VA 22908, USA
| |
Collapse
|