1
|
Zhang S, Zheng R, Long J, Zhu Y, Tan T. Computational design of carboxylase for the synthesis of 4-hydroxyisophthalic acid from p-hydroxybenzoic acid by fixing CO 2. JOURNAL OF ENVIRONMENTAL MANAGEMENT 2024; 366:121703. [PMID: 38996602 DOI: 10.1016/j.jenvman.2024.121703] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/28/2024] [Revised: 06/11/2024] [Accepted: 07/02/2024] [Indexed: 07/14/2024]
Abstract
Carbon dioxide (CO2) emissions constitute the primary contribution to global climate change. Synthetic CO2 fixation represents an exceptionally appealing and sustainable method for carbon neutralization. Unlike the limitations of chemical catalysis, biological CO2 fixation displays high selectivity and the ability to operate under mild conditions. The superfamily of amidohydrolases has demonstrated the ability to synthesize a range of aromatic monocarboxylic acids. However, there is a scarcity of reported carboxylases capable of synthesizing aromatic dicarboxylic acids. Among these, 4-hydroxyisophthalic acid holds significant potential for applications across various fields, yet no enzyme has been reported for its synthesis. In this study, we developed for the first time that exhibits starting activity in fixing CO2 to synthesize 4-hydroxyisophthalic acid. Furthermore, we have devised a computational strategy that effectively enhances the catalytic activity of this enzyme. A focused library comprising only 13 variants was generated. Experimental validation confirmed a threefold improvement in the carboxylation activity of the optimal variant (L47M). The computational enzyme design strategy proposed in this paper demonstrates broad applicability in developing carboxylases for synthesizing other aromatic dicarboxylic acids. This lays the groundwork for leveraging biocatalysis in industrial synthesis for CO2 fixation.
Collapse
Affiliation(s)
- Shiding Zhang
- College of Life Science and Technology, Beijing University of Chemical Technology, Beijing, 100029, China
| | - Ruonan Zheng
- College of Life Science and Technology, Beijing University of Chemical Technology, Beijing, 100029, China
| | - Jianyu Long
- College of Life Science and Technology, Beijing University of Chemical Technology, Beijing, 100029, China
| | - Yushan Zhu
- College of Life Science and Technology, Beijing University of Chemical Technology, Beijing, 100029, China; National Energy R&D Center for Biorefinery, Beijing University of Chemical Technology, Beijing, 100029, China.
| | - Tianwei Tan
- College of Life Science and Technology, Beijing University of Chemical Technology, Beijing, 100029, China; National Energy R&D Center for Biorefinery, Beijing University of Chemical Technology, Beijing, 100029, China.
| |
Collapse
|
2
|
Rakotoharisoa RV, Seifinoferest B, Zarifi N, Miller JDM, Rodriguez JM, Thompson MC, Chica RA. Design of Efficient Artificial Enzymes Using Crystallographically Enhanced Conformational Sampling. J Am Chem Soc 2024; 146:10001-10013. [PMID: 38532610 DOI: 10.1021/jacs.4c00677] [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: 03/28/2024]
Abstract
The ability to create efficient artificial enzymes for any chemical reaction is of great interest. Here, we describe a computational design method for increasing the catalytic efficiency of de novo enzymes by several orders of magnitude without relying on directed evolution and high-throughput screening. Using structural ensembles generated from dynamics-based refinement against X-ray diffraction data collected from crystals of Kemp eliminases HG3 (kcat/KM 125 M-1 s-1) and KE70 (kcat/KM 57 M-1 s-1), we design from each enzyme ≤10 sequences predicted to catalyze this reaction more efficiently. The most active designs display kcat/KM values improved by 100-250-fold, comparable to mutants obtained after screening thousands of variants in multiple rounds of directed evolution. Crystal structures show excellent agreement with computational models, with catalytic contacts present as designed and transition-state root-mean-square deviations of ≤0.65 Å. Our work shows how ensemble-based design can generate efficient artificial enzymes by exploiting the true conformational ensemble to design improved active sites.
Collapse
Affiliation(s)
- Rojo V Rakotoharisoa
- Department of Chemistry and Biomolecular Sciences, University of Ottawa, Ottawa, Ontario K1N 6N5, Canada
- Center for Catalysis Research and Innovation, University of Ottawa, Ottawa, Ontario K1N 6N5, Canada
| | - Behnoush Seifinoferest
- Department of Chemistry and Biochemistry, University of California Merced, Merced, California 95343, United States
| | - Niayesh Zarifi
- Department of Chemistry and Biomolecular Sciences, University of Ottawa, Ottawa, Ontario K1N 6N5, Canada
- Center for Catalysis Research and Innovation, University of Ottawa, Ottawa, Ontario K1N 6N5, Canada
| | - Jack D M Miller
- Department of Chemistry and Biomolecular Sciences, University of Ottawa, Ottawa, Ontario K1N 6N5, Canada
- Center for Catalysis Research and Innovation, University of Ottawa, Ottawa, Ontario K1N 6N5, Canada
| | - Joshua M Rodriguez
- Department of Chemistry and Biochemistry, University of California Merced, Merced, California 95343, United States
| | - Michael C Thompson
- Department of Chemistry and Biochemistry, University of California Merced, Merced, California 95343, United States
| | - Roberto A Chica
- Department of Chemistry and Biomolecular Sciences, University of Ottawa, Ottawa, Ontario K1N 6N5, Canada
- Center for Catalysis Research and Innovation, University of Ottawa, Ottawa, Ontario K1N 6N5, Canada
| |
Collapse
|
3
|
Lu L, Gou X, Tan SK, Mann SI, Yang H, Zhong X, Gazgalis D, Valdiviezo J, Jo H, Wu Y, Diolaiti ME, Ashworth A, Polizzi NF, DeGrado WF. De novo design of drug-binding proteins with predictable binding energy and specificity. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2023:2023.12.23.573178. [PMID: 38187746 PMCID: PMC10769398 DOI: 10.1101/2023.12.23.573178] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/09/2024]
Abstract
The de novo design of small-molecule-binding proteins has seen exciting recent progress; however, the ability to achieve exquisite affinity for binding small molecules while tuning specificity has not yet been demonstrated directly from computation. Here, we develop a computational procedure that results in the highest affinity binders to date with predetermined relative affinities, targeting a series of PARP1 inhibitors. Two of four designed proteins bound with affinities ranging from < 5 nM to low μM, in a predictable manner. X-ray crystal structures confirmed the accuracy of the designed protein-drug interactions. Molecular dynamics simulations informed the role of water in binding. Binding free-energy calculations performed directly on the designed models are in excellent agreement with the experimentally measured affinities, suggesting that the de novo design of small-molecule-binding proteins with tuned interaction energies is now feasible entirely from computation. We expect these methods to open many opportunities in biomedicine, including rapid sensor development, antidote design, and drug delivery vehicles.
Collapse
Affiliation(s)
- Lei Lu
- Department of Pharmaceutical Chemistry, University of California, San Francisco, San Francisco, CA 94158, USA
| | - Xuxu Gou
- Helen Diller Family Comprehensive Cancer Center, University of California, San Francisco, CA 94158, USA
| | - Sophia K Tan
- Department of Pharmaceutical Chemistry, University of California, San Francisco, San Francisco, CA 94158, USA
| | - Samuel I. Mann
- Department of Pharmaceutical Chemistry, University of California, San Francisco, San Francisco, CA 94158, USA
| | - Hyunjun Yang
- Department of Pharmaceutical Chemistry, University of California, San Francisco, San Francisco, CA 94158, USA
| | | | - Dimitrios Gazgalis
- Department of Cellular and Molecular Pharmacology, University of California, San Francisco, CA 94158, USA
- Dana Farber Cancer Institute, Harvard Medical School, Boston, MA 02215, USA
| | - Jesús Valdiviezo
- Dana Farber Cancer Institute, Harvard Medical School, Boston, MA 02215, USA
| | - Hyunil Jo
- Department of Pharmaceutical Chemistry, University of California, San Francisco, San Francisco, CA 94158, USA
| | - Yibing Wu
- Department of Pharmaceutical Chemistry, University of California, San Francisco, San Francisco, CA 94158, USA
| | - Morgan E. Diolaiti
- Helen Diller Family Comprehensive Cancer Center, University of California, San Francisco, CA 94158, USA
| | - Alan Ashworth
- Helen Diller Family Comprehensive Cancer Center, University of California, San Francisco, CA 94158, USA
| | | | - William F. DeGrado
- Department of Pharmaceutical Chemistry, University of California, San Francisco, San Francisco, CA 94158, USA
| |
Collapse
|
4
|
Zhang J, Wang H, Luo Z, Yang Z, Zhang Z, Wang P, Li M, Zhang Y, Feng Y, Lu D, Zhu Y. Computational design of highly efficient thermostable MHET hydrolases and dual enzyme system for PET recycling. Commun Biol 2023; 6:1135. [PMID: 37945666 PMCID: PMC10636135 DOI: 10.1038/s42003-023-05523-5] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/09/2023] [Accepted: 10/30/2023] [Indexed: 11/12/2023] Open
Abstract
Recently developed enzymes for the depolymerization of polyethylene terephthalate (PET) such as FAST-PETase and LCC-ICCG are inhibited by the intermediate PET product mono(2-hydroxyethyl) terephthalate (MHET). Consequently, the conversion of PET enzymatically into its constituent monomers terephthalic acid (TPA) and ethylene glycol (EG) is inefficient. In this study, a protein scaffold (1TQH) corresponding to a thermophilic carboxylesterase (Est30) was selected from the structural database and redesigned in silico. Among designs, a double variant KL-MHETase (I171K/G130L) with a similar protein melting temperature (67.58 °C) to that of the PET hydrolase FAST-PETase (67.80 °C) exhibited a 67-fold higher activity for MHET hydrolysis than FAST-PETase. A fused dual enzyme system comprising KL-MHETase and FAST-PETase exhibited a 2.6-fold faster PET depolymerization rate than FAST-PETase alone. Synergy increased the yield of TPA by 1.64 fold, and its purity in the released aromatic products reached 99.5%. In large reaction systems with 100 g/L substrate concentrations, the dual enzyme system KL36F achieved over 90% PET depolymerization into monomers, demonstrating its potential applicability in the industrial recycling of PET plastics. Therefore, a dual enzyme system can greatly reduce the reaction and separation cost for sustainable enzymatic PET recycling.
Collapse
Affiliation(s)
- Jun Zhang
- College of Life Science and Technology, Beijing University of Chemical Technology, Beijing, 100029, China
- Department of Chemical Engineering, Tsinghua University, Beijing, 100084, China
| | - Hongzhao Wang
- College of Life Science and Technology, Beijing University of Chemical Technology, Beijing, 100029, China
| | - Zhaorong Luo
- Beijing Advanced Innovation Center for Soft Matter Science and Engineering, State Key Laboratory of Chemical Resource Engineering, Beijing University of Chemical Technology, Beijing, 100029, China
| | - Zhenwu Yang
- College of Life Science and Technology, Beijing University of Chemical Technology, Beijing, 100029, China
| | - Zixuan Zhang
- College of Life Science and Technology, Beijing University of Chemical Technology, Beijing, 100029, China
| | - Pengyu Wang
- Department of Chemical Engineering, Tsinghua University, Beijing, 100084, China
| | - Mengyu Li
- College of Life Science and Technology, Beijing University of Chemical Technology, Beijing, 100029, China
| | - Yi Zhang
- Beijing Advanced Innovation Center for Soft Matter Science and Engineering, State Key Laboratory of Chemical Resource Engineering, Beijing University of Chemical Technology, Beijing, 100029, China
| | - Yue Feng
- Beijing Advanced Innovation Center for Soft Matter Science and Engineering, State Key Laboratory of Chemical Resource Engineering, Beijing University of Chemical Technology, Beijing, 100029, China
| | - Diannan Lu
- Department of Chemical Engineering, Tsinghua University, Beijing, 100084, China.
| | - Yushan Zhu
- College of Life Science and Technology, Beijing University of Chemical Technology, Beijing, 100029, China.
- National Energy R&D Center for Biorefinery, Beijing University of Chemical Technology, Beijing, 100029, China.
| |
Collapse
|
5
|
Rakotoharisoa RV, Seifinoferest B, Zarifi N, Miller JD, Rodriguez JM, Thompson MC, Chica RA. Design of efficient artificial enzymes using crystallographically-enhanced conformational sampling. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2023:2023.11.01.564846. [PMID: 37961474 PMCID: PMC10635043 DOI: 10.1101/2023.11.01.564846] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/15/2023]
Abstract
The ability to create efficient artificial enzymes for any chemical reaction is of great interest. Here, we describe a computational design method for increasing catalytic efficiency of de novo enzymes to a level comparable to their natural counterparts without relying on directed evolution. Using structural ensembles generated from dynamics-based refinement against X-ray diffraction data collected from crystals of Kemp eliminases HG3 (kcat/KM 125 M-1 s-1) and KE70 (kcat/KM 57 M-1 s-1), we design from each enzyme ≤10 sequences predicted to catalyze this reaction more efficiently. The most active designs display kcat/KM values improved by 100-250-fold, comparable to mutants obtained after screening thousands of variants in multiple rounds of directed evolution. Crystal structures show excellent agreement with computational models. Our work shows how computational design can generate efficient artificial enzymes by exploiting the true conformational ensemble to more effectively stabilize the transition state.
Collapse
Affiliation(s)
- Rojo V. Rakotoharisoa
- Department of Chemistry and Biomolecular Sciences, University of Ottawa, Ottawa, Ontario, Canada, K1N 6N5
- Center for Catalysis Research and Innovation, University of Ottawa, Ottawa, Ontario, Canada, K1N 6N5
| | - Behnoush Seifinoferest
- Department of Chemistry and Biochemistry, University of California Merced, Merced, California 95343, United States
| | - Niayesh Zarifi
- Department of Chemistry and Biomolecular Sciences, University of Ottawa, Ottawa, Ontario, Canada, K1N 6N5
- Center for Catalysis Research and Innovation, University of Ottawa, Ottawa, Ontario, Canada, K1N 6N5
| | - Jack D.M. Miller
- Department of Chemistry and Biomolecular Sciences, University of Ottawa, Ottawa, Ontario, Canada, K1N 6N5
- Center for Catalysis Research and Innovation, University of Ottawa, Ottawa, Ontario, Canada, K1N 6N5
| | - Joshua M. Rodriguez
- Department of Chemistry and Biochemistry, University of California Merced, Merced, California 95343, United States
| | - Michael C. Thompson
- Department of Chemistry and Biochemistry, University of California Merced, Merced, California 95343, United States
| | - Roberto A. Chica
- Department of Chemistry and Biomolecular Sciences, University of Ottawa, Ottawa, Ontario, Canada, K1N 6N5
- Center for Catalysis Research and Innovation, University of Ottawa, Ottawa, Ontario, Canada, K1N 6N5
| |
Collapse
|
6
|
Kratochvil HT, Watkins LC, Mravic M, Thomaston JL, Nicoludis JM, Somberg NH, Liu L, Hong M, Voth GA, DeGrado WF. Transient water wires mediate selective proton transport in designed channel proteins. Nat Chem 2023; 15:1012-1021. [PMID: 37308712 PMCID: PMC10475958 DOI: 10.1038/s41557-023-01210-4] [Citation(s) in RCA: 11] [Impact Index Per Article: 11.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/19/2022] [Accepted: 04/19/2023] [Indexed: 06/14/2023]
Abstract
Selective proton transport through proteins is essential for forming and using proton gradients in cells. Protons are conducted along hydrogen-bonded 'wires' of water molecules and polar side chains, which, somewhat surprisingly, are often interrupted by dry apolar stretches in the conduction pathways, inferred from static protein structures. Here we hypothesize that protons are conducted through such dry spots by forming transient water wires, often highly correlated with the presence of the excess protons in the water wire. To test this hypothesis, we performed molecular dynamics simulations to design transmembrane channels with stable water pockets interspersed by apolar segments capable of forming flickering water wires. The minimalist designed channels conduct protons at rates similar to viral proton channels, and they are at least 106-fold more selective for H+ over Na+. These studies inform the mechanisms of biological proton conduction and the principles for engineering proton-conductive materials.
Collapse
Affiliation(s)
- Huong T Kratochvil
- Department of Pharmaceutical Chemistry, University of California-San Francisco, San Francisco, CA, USA.
- Department of Chemistry, University of North Carolina at Chapel Hill, Chapel Hill, NC, USA.
| | - Laura C Watkins
- Department of Chemistry, Chicago Center for Theoretical Chemistry, Institute for Biophysical Dynamics and James Franck Institute, The University of Chicago, Chicago, IL, USA
- Kemper Insurance, Chicago, IL, USA
| | - Marco Mravic
- Department of Pharmaceutical Chemistry, University of California-San Francisco, San Francisco, CA, USA
- Department of Integrative Structural and Computational Biology Scripps Research Institute, La Jolla, CA, USA
| | - Jessica L Thomaston
- Department of Pharmaceutical Chemistry, University of California-San Francisco, San Francisco, CA, USA
| | - John M Nicoludis
- Department of Pharmaceutical Chemistry, University of California-San Francisco, San Francisco, CA, USA
- Genentech, San Francisco, CA, USA
| | - Noah H Somberg
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, MA, USA
| | | | - Mei Hong
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, MA, USA
| | - Gregory A Voth
- Department of Chemistry, Chicago Center for Theoretical Chemistry, Institute for Biophysical Dynamics and James Franck Institute, The University of Chicago, Chicago, IL, USA.
| | - William F DeGrado
- Department of Pharmaceutical Chemistry, University of California-San Francisco, San Francisco, CA, USA.
| |
Collapse
|
7
|
Wang P, Zhang J, Zhang S, Lu D, Zhu Y. Using High-Throughput Molecular Dynamics Simulation to Enhance the Computational Design of Kemp Elimination Enzymes. J Chem Inf Model 2023; 63:1323-1337. [PMID: 36782360 DOI: 10.1021/acs.jcim.3c00002] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/15/2023]
Abstract
Computational enzyme design has been successfully applied to identify new alternatives to natural enzymes for the biosynthesis of important compounds. However, the moderate catalytic activities of de novo designed enzymes indicate that the modeling accuracy of current computational enzyme design methods should be improved. Here, high-throughput molecular dynamics simulations were used to enhance computational enzyme design, thus allowing the identification of variants with higher activities in silico. Different time schemes of high-throughput molecular dynamics simulations were tested to identify the catalytic features of evolved Kemp eliminases. The 20 × 1 ns molecular dynamics simulation scheme was sufficiently accurate and computationally viable to screen the computationally designed massive variants of Kemp elimination enzymes. The developed hybrid computational strategy was used to redesign the most active Kemp eliminase, HG3.17, and five variants were generated and experimentally confirmed to afford higher catalytic efficiencies than that of HG3.17, with one double variant (D52Q/A53S) exhibiting a 55% increase. The hybrid computational enzyme design strategy is general and computationally economical, with which we anticipate the efficient creation of practical enzymes for industrial biocatalysis.
Collapse
Affiliation(s)
- Pengyu Wang
- College of Life Science and Technology, Beijing University of Chemical Technology, Beijing 100029, China.,Department of Chemical Engineering, Tsinghua University, Beijing 100084, China
| | - Jun Zhang
- Department of Chemical Engineering, Tsinghua University, Beijing 100084, China
| | - Shengyu Zhang
- College of Life Science and Technology, Beijing University of Chemical Technology, Beijing 100029, China
| | - Diannan Lu
- Department of Chemical Engineering, Tsinghua University, Beijing 100084, China
| | - Yushan Zhu
- College of Life Science and Technology, Beijing University of Chemical Technology, Beijing 100029, China
| |
Collapse
|
8
|
Zhang S, Zhang J, Luo W, Wang P, Zhu Y. A preorganization oriented computational method for de novo design of Kemp elimination enzymes. Enzyme Microb Technol 2022; 160:110093. [DOI: 10.1016/j.enzmictec.2022.110093] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/07/2022] [Revised: 06/13/2022] [Accepted: 06/30/2022] [Indexed: 11/26/2022]
|
9
|
Kretschmer S, Kortemme T. Advances in the Computational Design of Small-Molecule-Controlled Protein-Based Circuits for Synthetic Biology. PROCEEDINGS OF THE IEEE. INSTITUTE OF ELECTRICAL AND ELECTRONICS ENGINEERS 2022; 110:659-674. [PMID: 36531560 PMCID: PMC9754107 DOI: 10.1109/jproc.2022.3157898] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/16/2023]
Abstract
Synthetic biology approaches living systems with an engineering perspective and promises to deliver solutions to global challenges in healthcare and sustainability. A critical component is the design of biomolecular circuits with programmable input-output behaviors. Such circuits typically rely on a sensor module that recognizes molecular inputs, which is coupled to a functional output via protein-level circuits or regulating the expression of a target gene. While gene expression outputs can be customized relatively easily by exchanging the target genes, sensing new inputs is a major limitation. There is a limited repertoire of sensors found in nature, and there are often difficulties with interfacing them with engineered circuits. Computational protein design could be a key enabling technology to address these challenges, as it allows for the engineering of modular and tunable sensors that can be tailored to the circuit's application. In this article, we review recent computational approaches to design protein-based sensors for small-molecule inputs with particular focus on those based on the widely used Rosetta software suite. Furthermore, we review mechanisms that have been harnessed to couple ligand inputs to functional outputs. Based on recent literature, we illustrate how the combination of protein design and synthetic biology enables new sensors for diverse applications ranging from biomedicine to metabolic engineering. We conclude with a perspective on how strategies to address frontiers in protein design and cellular circuit design may enable the next generation of sense-response networks, which may increasingly be assembled from de novo components to display diverse and engineerable input-output behaviors.
Collapse
Affiliation(s)
- Simon Kretschmer
- Department of Bioengineering and Therapeutic Sciences, University of California at San Francisco, San Francisco, CA 94158 USA, and affiliated with the California Quantitative Biosciences Institute (QBI) at UCSF, San Francisco, CA 94158 USA
| | - Tanja Kortemme
- Department of Bioengineering and Therapeutic Sciences, University of California at San Francisco, San Francisco, CA 94158 USA, and affiliated with the California Quantitative Biosciences Institute (QBI) at UCSF, San Francisco, CA 94158 USA
| |
Collapse
|
10
|
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]
|
11
|
Pan X, Kortemme T. De novo protein fold families expand the designable ligand binding site space. PLoS Comput Biol 2021; 17:e1009620. [PMID: 34807909 PMCID: PMC8648124 DOI: 10.1371/journal.pcbi.1009620] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/13/2021] [Revised: 12/06/2021] [Accepted: 11/08/2021] [Indexed: 11/19/2022] Open
Abstract
A major challenge in designing proteins de novo to bind user-defined ligands with high affinity is finding backbones structures into which a new binding site geometry can be engineered with high precision. Recent advances in methods to generate protein fold families de novo have expanded the space of accessible protein structures, but it is not clear to what extend de novo proteins with diverse geometries also expand the space of designable ligand binding functions. We constructed a library of 25,806 high-quality ligand binding sites and developed a fast protocol to place (“match”) these binding sites into both naturally occurring and de novo protein families with two fold topologies: Rossman and NTF2. Each matching step involves engineering new binding site residues into each protein “scaffold”, which is distinct from the problem of comparing already existing binding pockets. 5,896 and 7,475 binding sites could be matched to the Rossmann and NTF2 fold families, respectively. De novo designed Rossman and NTF2 protein families can support 1,791 and 678 binding sites that cannot be matched to naturally existing structures with the same topologies, respectively. While the number of protein residues in ligand binding sites is the major determinant of matching success, ligand size and primary sequence separation of binding site residues also play important roles. The number of matched binding sites are power law functions of the number of members in a fold family. Our results suggest that de novo sampling of geometric variations on diverse fold topologies can significantly expand the space of designable ligand binding sites for a wealth of possible new protein functions. De novo design of proteins that can bind to novel and highly diverse user-defined small molecule ligands could have broad biomedical and synthetic biology applications. Because ligand binding site geometries need to be accommodated by protein backbone scaffolds at high accuracy, the diversity of scaffolds is a major limitation for designing new ligand binding functions. Advances in computational protein structure design methods have significantly increased the number of accessible stable scaffold structures. Understanding how many new ligand binding sites can be designed into the de novo scaffolds is important for engineering novel ligand binding proteins. To answer this question, we constructed a large library of ligand binding sites from the Protein Data Bank (PDB). We tested the number of ligand binding sites that can be designed into de novo scaffolds and naturally existing scaffolds with the same fold topologies. The results showed that de novo scaffolds significantly expanded the potential ligand binding space of their respective fold topologies. We also identified factors that affect difficulties of binding site accommodation, as well as the relationship between the number of scaffolds and the accessible ligand binding site space. We believe our findings will benefit future method development and applications of ligand binding protein design.
Collapse
Affiliation(s)
- Xingjie Pan
- Department of Bioengineering and Therapeutic Sciences, University of California San Francisco, San Francisco, California, United States of America
- UC Berkeley–UCSF Graduate Program in Bioengineering, University of California San Francisco, San Francisco, California, United States of America
- * E-mail: (XP); (TK)
| | - Tanja Kortemme
- Department of Bioengineering and Therapeutic Sciences, University of California San Francisco, San Francisco, California, United States of America
- UC Berkeley–UCSF Graduate Program in Bioengineering, University of California San Francisco, San Francisco, California, United States of America
- Quantitative Biosciences Institute (QBI), University of California San Francisco, San Francisco, California, United States of America
- Chan Zuckerberg Biohub, San Francisco, California, United States of America
- * E-mail: (XP); (TK)
| |
Collapse
|
12
|
Wang P, Zhang S, Zhang J, Zhu Y. Computational design of penicillin acylase variants with improved kinetic selectivity for the enzymatic synthesis of cefazolin. Biochem Eng J 2021. [DOI: 10.1016/j.bej.2021.108149] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/04/2023]
|
13
|
Zhang J, Pang Q, Wang Q, Qi Q, Wang Q. Modular tuning engineering and versatile applications of genetically encoded biosensors. Crit Rev Biotechnol 2021; 42:1010-1027. [PMID: 34615431 DOI: 10.1080/07388551.2021.1982858] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/20/2022]
Abstract
Genetically encoded biosensors have a diverse range of detectable signals and potential applications in many fields, including metabolism control and high-throughput screening. Their ability to be used in situ with minimal interference to the bioprocess of interest could revolutionize synthetic biology and microbial cell factories. The performance and functions of these biosensors have been extensively studied and have been rapidly improved. We review here current biosensor tuning strategies and attempt to unravel how to obtain ideal biosensor functions through experimental adjustments. Strategies for expanding the biosensor input signals that increases the number of detectable compounds have also been summarized. Finally, different output signals and their practical requirements for biotechnology and biomedical applications and environmental safety concerns have been analyzed. This in-depth review of the responses and regulation mechanisms of genetically encoded biosensors will assist to improve their design and optimization in various application scenarios.
Collapse
Affiliation(s)
- Jian Zhang
- National Glycoengineering Research Center, State Key Laboratory of Microbial Technology, Shandong University, Qingdao, P. R. China
| | - Qingxiao Pang
- National Glycoengineering Research Center, State Key Laboratory of Microbial Technology, Shandong University, Qingdao, P. R. China
| | - Qi Wang
- National Glycoengineering Research Center, State Key Laboratory of Microbial Technology, Shandong University, Qingdao, P. R. China
| | - Qingsheng Qi
- National Glycoengineering Research Center, State Key Laboratory of Microbial Technology, Shandong University, Qingdao, P. R. China.,CAS Key Lab of Biobased Materials, Qingdao Institute of Bioenergy and Bioprocess Technology, Chinese Academy of Sciences, Qingdao, P. R. China
| | - Qian Wang
- National Glycoengineering Research Center, State Key Laboratory of Microbial Technology, Shandong University, Qingdao, P. R. China.,CAS Key Lab of Biobased Materials, Qingdao Institute of Bioenergy and Bioprocess Technology, Chinese Academy of Sciences, Qingdao, P. R. China
| |
Collapse
|
14
|
Zhao T, Yong X, Zhao Z, Dolce V, Li Y, Curcio R. Research status of Bacillus phytase. 3 Biotech 2021; 11:415. [PMID: 34485008 DOI: 10.1007/s13205-021-02964-9] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/28/2021] [Accepted: 07/27/2021] [Indexed: 10/20/2022] Open
Abstract
Phytic acid is abundant in seeds, roots and stems of plants, it acts as an anti-nutrient in food and feed industry, since it affects the absorption of nutrients by humans and monogastric animals. Furthermore, phosphorus produced through its decomposition by microorganisms can cause environmental pollution. Phytase degrades phytic acid generating precursors of inositol that can be used in clinical practice; in addition, phytase treatment can minimize the anti-nutritional effect of phytic acid. The use of phytase synthesized from Bacillus is more advantageous due to its high activity. Additionally, its good heat resistance under neutral conditions greatly fills the gap of commercial utilization of acid phytase. In this review, we summarize the latest research results on Bacillus phytase, including its physiological and biochemical characteristics, molecular structure information, calcium effects on its catalytic activity and stability, its catalytic mechanism and molecular modification.
Collapse
Affiliation(s)
- Ting Zhao
- College of Life Science and Technology, Xinjiang University, Urumqi, China
- National Facility for Protein Science in Shanghai, Zhangjiang Lab, Shanghai Advanced Research Institute, Chinese Academy of Science, Shanghai, China
| | - Xihao Yong
- College of Life Science and Technology, Xinjiang University, Urumqi, China
- Faculty of Bioengineering, Sichuan University of Science and Engineering, Zigong, People's Republic of China
| | - Ziming Zhao
- Faculty of Bioengineering, Sichuan University of Science and Engineering, Zigong, People's Republic of China
| | - Vincenza Dolce
- Department of Pharmacy, Health and Nutritional Sciences, University of Calabria, Rende, Italy
| | - Yuan Li
- College of Life Science and Technology, Xinjiang University, Urumqi, China
- Department of Bioscience, Biotechnology and Biopharmaceutics, University of Bari, Bari, Italy
| | - Rosita Curcio
- Department of Pharmacy, Health and Nutritional Sciences, University of Calabria, Rende, Italy
| |
Collapse
|
15
|
Pagar AD, Patil MD, Flood DT, Yoo TH, Dawson PE, Yun H. Recent Advances in Biocatalysis with Chemical Modification and Expanded Amino Acid Alphabet. Chem Rev 2021; 121:6173-6245. [PMID: 33886302 DOI: 10.1021/acs.chemrev.0c01201] [Citation(s) in RCA: 41] [Impact Index Per Article: 13.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Abstract
The two main strategies for enzyme engineering, directed evolution and rational design, have found widespread applications in improving the intrinsic activities of proteins. Although numerous advances have been achieved using these ground-breaking methods, the limited chemical diversity of the biopolymers, restricted to the 20 canonical amino acids, hampers creation of novel enzymes that Nature has never made thus far. To address this, much research has been devoted to expanding the protein sequence space via chemical modifications and/or incorporation of noncanonical amino acids (ncAAs). This review provides a balanced discussion and critical evaluation of the applications, recent advances, and technical breakthroughs in biocatalysis for three approaches: (i) chemical modification of cAAs, (ii) incorporation of ncAAs, and (iii) chemical modification of incorporated ncAAs. Furthermore, the applications of these approaches and the result on the functional properties and mechanistic study of the enzymes are extensively reviewed. We also discuss the design of artificial enzymes and directed evolution strategies for enzymes with ncAAs incorporated. Finally, we discuss the current challenges and future perspectives for biocatalysis using the expanded amino acid alphabet.
Collapse
Affiliation(s)
- Amol D Pagar
- Department of Systems Biotechnology, Konkuk University, 120 Neungdong-ro, Gwangjin-gu, Seoul 05029, Korea
| | - Mahesh D Patil
- Department of Systems Biotechnology, Konkuk University, 120 Neungdong-ro, Gwangjin-gu, Seoul 05029, Korea
| | - Dillon T Flood
- Department of Chemistry, The Scripps Research Institute, 10550 North Torrey Pines Road, La Jolla, California 92037, United States
| | - Tae Hyeon Yoo
- Department of Molecular Science and Technology, Ajou University, 206 World cup-ro, Yeongtong-gu, Suwon 16499, Korea
| | - Philip E Dawson
- Department of Chemistry, The Scripps Research Institute, 10550 North Torrey Pines Road, La Jolla, California 92037, United States
| | - Hyungdon Yun
- Department of Systems Biotechnology, Konkuk University, 120 Neungdong-ro, Gwangjin-gu, Seoul 05029, Korea
| |
Collapse
|
16
|
Ghattas W, Mahy JP, Réglier M, Simaan AJ. Artificial Enzymes for Diels-Alder Reactions. Chembiochem 2020; 22:443-459. [PMID: 32852088 DOI: 10.1002/cbic.202000316] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/15/2020] [Revised: 08/17/2020] [Indexed: 12/13/2022]
Abstract
The Diels-Alder (DA) reaction is a cycloaddition of a conjugated diene and an alkene (dienophile) leading to the formation of a cyclohexene derivative through a concerted mechanism. As DA reactions generally proceed with a high degree of regio- and stereoselectivity, they are widely used in synthetic organic chemistry. Considering eco-conscious public and governmental movements, efforts are now directed towards the development of synthetic processes that meet environmental concerns. Artificial enzymes, which can be developed to catalyze abiotic reactions, appear to be important synthetic tools in the synthetic biology field. This review describes the different strategies used to develop protein-based artificial enzymes for DA reactions, including for in cellulo approaches.
Collapse
Affiliation(s)
- Wadih Ghattas
- Institut de Chimie Moléculaire et des Matériaux d'Orsay (ICMMO), UMR 8182 CNRS, Université Paris Sud, Université Paris-Saclay, Orsay, 91405 Cedex 8, France
| | - Jean-Pierre Mahy
- Institut de Chimie Moléculaire et des Matériaux d'Orsay (ICMMO), UMR 8182 CNRS, Université Paris Sud, Université Paris-Saclay, Orsay, 91405 Cedex 8, France
| | - Marius Réglier
- Aix Marseille Univ, CNRS, Centrale Marseille, iSm2, Avenue Escadrille Normandie Niemen, Service 342, Marseille, 13397, France
| | - A Jalila Simaan
- Aix Marseille Univ, CNRS, Centrale Marseille, iSm2, Avenue Escadrille Normandie Niemen, Service 342, Marseille, 13397, France
| |
Collapse
|
17
|
Lucas JE, Kortemme T. New computational protein design methods for de novo small molecule binding sites. PLoS Comput Biol 2020; 16:e1008178. [PMID: 33017412 PMCID: PMC7575090 DOI: 10.1371/journal.pcbi.1008178] [Citation(s) in RCA: 16] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/30/2020] [Revised: 10/20/2020] [Accepted: 07/22/2020] [Indexed: 11/19/2022] Open
Abstract
Protein binding to small molecules is fundamental to many biological processes, yet it remains challenging to predictively design this functionality de novo. Current state-of-the-art computational design methods typically rely on existing small molecule binding sites or protein scaffolds with existing shape complementarity for a target ligand. Here we introduce new methods that utilize pools of discrete contacts between protein side chains and defined small molecule ligand substructures (ligand fragments) observed in the Protein Data Bank. We use the Rosetta Molecular Modeling Suite to recombine protein side chains in these contact pools to generate hundreds of thousands of energetically favorable binding sites for a target ligand. These composite binding sites are built into existing scaffold proteins matching the intended binding site geometry with high accuracy. In addition, we apply pools of side chain rotamers interacting with the target ligand to augment Rosetta's conventional design machinery and improve key metrics known to be predictive of design success. We demonstrate that our method reliably builds diverse binding sites into different scaffold proteins for a variety of target molecules. Our generalizable de novo ligand binding site design method provides a foundation for versatile design of protein to interface previously unattainable molecules for applications in medical diagnostics and synthetic biology.
Collapse
Affiliation(s)
- James E. Lucas
- UC Berkeley–UCSF Graduate Program in Bioengineering, University of California San Francisco, San Francisco, CA, United States of America
- Department of Bioengineering and Therapeutic Sciences, University of California San Francisco, San Francisco, CA, United States of America
| | - Tanja Kortemme
- UC Berkeley–UCSF Graduate Program in Bioengineering, University of California San Francisco, San Francisco, CA, United States of America
- Department of Bioengineering and Therapeutic Sciences, University of California San Francisco, San Francisco, CA, United States of America
| |
Collapse
|
18
|
Ensemble-based enzyme design can recapitulate the effects of laboratory directed evolution in silico. Nat Commun 2020; 11:4808. [PMID: 32968058 PMCID: PMC7511930 DOI: 10.1038/s41467-020-18619-x] [Citation(s) in RCA: 57] [Impact Index Per Article: 14.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/27/2020] [Accepted: 08/25/2020] [Indexed: 01/30/2023] Open
Abstract
The creation of artificial enzymes is a key objective of computational protein design. Although de novo enzymes have been successfully designed, these exhibit low catalytic efficiencies, requiring directed evolution to improve activity. Here, we use room-temperature X-ray crystallography to study changes in the conformational ensemble during evolution of the designed Kemp eliminase HG3 (kcat/KM 146 M−1s−1). We observe that catalytic residues are increasingly rigidified, the active site becomes better pre-organized, and its entrance is widened. Based on these observations, we engineer HG4, an efficient biocatalyst (kcat/KM 103,000 M−1s−1) containing key first and second-shell mutations found during evolution. HG4 structures reveal that its active site is pre-organized and rigidified for efficient catalysis. Our results show how directed evolution circumvents challenges inherent to enzyme design by shifting conformational ensembles to favor catalytically-productive sub-states, and suggest improvements to the design methodology that incorporate ensemble modeling of crystallographic data. Kemp eliminases are artificial enzymes that catalyze the concerted deprotonation and ring-opening of benzisoxazoles. Here, the authors use room-temperature X-ray crystallography to investigate changes to the conformational ensemble of the Kemp eliminase HG3 along a directed evolutionary trajectory, and develop an experimentally guided, ensemble-based computational enzyme design procedure.
Collapse
|
19
|
Polizzi NF, DeGrado WF. A defined structural unit enables de novo design of small-molecule-binding proteins. Science 2020; 369:1227-1233. [PMID: 32883865 PMCID: PMC7526616 DOI: 10.1126/science.abb8330] [Citation(s) in RCA: 52] [Impact Index Per Article: 13.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/21/2020] [Accepted: 06/29/2020] [Indexed: 12/21/2022]
Abstract
The de novo design of proteins that bind highly functionalized small molecules represents a great challenge. To enable computational design of binders, we developed a unit of protein structure-a van der Mer (vdM)-that maps the backbone of each amino acid to statistically preferred positions of interacting chemical groups. Using vdMs, we designed six de novo proteins to bind the drug apixaban; two bound with low and submicromolar affinity. X-ray crystallography and mutagenesis confirmed a structure with a precisely designed cavity that forms favorable interactions in the drug-protein complex. vdMs may enable design of functional proteins for applications in sensing, medicine, and catalysis.
Collapse
Affiliation(s)
- Nicholas F Polizzi
- Department of Pharmaceutical Chemistry, Cardiovascular Research Institute, University of California, San Francisco, San Francisco, CA 94158, USA.
| | - William F DeGrado
- Department of Pharmaceutical Chemistry, Cardiovascular Research Institute, University of California, San Francisco, San Francisco, CA 94158, USA.
| |
Collapse
|
20
|
Peacock A. Can proteins be truly designed sans function? Science 2020; 369:1166-1167. [DOI: 10.1126/science.abd4791] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/02/2022]
Abstract
A new unit of local protein structure can aid in the de novo design of ligand-binding proteins
Collapse
Affiliation(s)
- Anna Peacock
- School of Chemistry, University of Birmingham, Birmingham, UK
| |
Collapse
|
21
|
|
22
|
Xue J, Wang P, Kuang J, Zhu Y. Computational design of new enzymes for hydrolysis and synthesis of third-generation cephalosporin antibiotics. Enzyme Microb Technol 2020; 140:109649. [PMID: 32912699 DOI: 10.1016/j.enzmictec.2020.109649] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/14/2020] [Revised: 08/07/2020] [Accepted: 08/08/2020] [Indexed: 10/23/2022]
Abstract
Engineering active sites in inert scaffolds to catalyze chemical transformations with unnatural substrates is still a great challenge for enzyme catalysis. In this research, a p-nitrobenzyl esterase from Bacillus subtilis was identified from the structural database, and a double mutant E115A/E188A was designed to afford catalytic activities toward the hydrolysis of ceftizoxime. A quadruple mutant E115A/E188A/L362S/I270A with enhanced catalytic efficiency was created to catalyze the condensation reaction of ethyl-2-methoxy-amino-2-(2-aminothiazole-4-yl) acetate with 7-amino-3-nor-cephalosporanic acid to produce ceftizoxime in a fully aqueous medium. The catalytic efficiencies of the computationally designed mutants E115A/E188A/L362S/I270A and E115A/Y118 K/E188 V/I270A/L362S can be taken as starting points to further improve their properties towards the practical application in designing more ecology-friendly production of third-generation cephalosporins.
Collapse
Affiliation(s)
- Jing Xue
- Department of Chemical Engineering, Tsinghua University, Beijing, 100084, China
| | - Pengyu Wang
- Department of Chemical Engineering, Tsinghua University, Beijing, 100084, China
| | - Jianyong Kuang
- Department of Chemical Engineering, Tsinghua University, Beijing, 100084, China
| | - Yushan Zhu
- Department of Chemical Engineering, Tsinghua University, Beijing, 100084, China; MOE Key Lab of Industrial Biocatalysis, Tsinghua University, Beijing, 100084, China.
| |
Collapse
|
23
|
Weitzner BD, Kipnis Y, Daniel AG, Hilvert D, Baker D. A computational method for design of connected catalytic networks in proteins. Protein Sci 2020; 28:2036-2041. [PMID: 31642127 PMCID: PMC6863703 DOI: 10.1002/pro.3757] [Citation(s) in RCA: 22] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/27/2019] [Revised: 10/21/2019] [Accepted: 10/21/2019] [Indexed: 02/05/2023]
Abstract
Computational design of new active sites has generally proceeded by geometrically defining interactions between the reaction transition state(s) and surrounding side‐chain functional groups which maximize transition‐state stabilization, and then searching for sites in protein scaffolds where the specified side‐chain–transition‐state interactions can be realized. A limitation of this approach is that the interactions between the side chains themselves are not constrained. An extensive connected hydrogen bond network involving the catalytic residues was observed in a designed retroaldolase following directed evolution. Such connected networks could increase catalytic activity by preorganizing active site residues in catalytically competent orientations, and enabling concerted interactions between side chains during catalysis, for example, proton shuffling. We developed a method for designing active sites in which the catalytic side chains, in addition to making interactions with the transition state, are also involved in extensive hydrogen bond networks. Because of the added constraint of hydrogen‐bond connectivity between the catalytic side chains, to find solutions, a wider range of interactions between these side chains and the transition state must be considered. Our new method starts from a ChemDraw‐like two‐dimensional representation of the transition state with hydrogen‐bond donors, acceptors, and covalent interaction sites indicated, and all placements of side‐chain functional groups that make the indicated interactions with the transition state, and are fully connected in a single hydrogen‐bond network are systematically enumerated. The RosettaMatch method can then be used to identify realizations of these fully‐connected active sites in protein scaffolds. The method generates many fully‐connected active site solutions for a set of model reactions that are promising starting points for the design of fully‐preorganized enzyme catalysts.
Collapse
Affiliation(s)
- Brian D Weitzner
- Department of Biochemistry, University of Washington, Seattle, Washington.,Institute for Protein Design, University of Washington, Seattle, Washington
| | - Yakov Kipnis
- Department of Biochemistry, University of Washington, Seattle, Washington.,Institute for Protein Design, University of Washington, Seattle, Washington.,Howard Hughes Medical Institute, University of Washington, Seattle, Washington
| | - A Gerard Daniel
- Department of Biochemistry, University of Washington, Seattle, Washington.,Institute for Protein Design, University of Washington, Seattle, Washington
| | - Donald Hilvert
- Laboratory of Organic Chemistry, ETH Zurich, Zurich, Switzerland
| | - David Baker
- Department of Biochemistry, University of Washington, Seattle, Washington.,Institute for Protein Design, University of Washington, Seattle, Washington.,Howard Hughes Medical Institute, University of Washington, Seattle, Washington
| |
Collapse
|
24
|
Ali M, Ishqi HM, Husain Q. Enzyme engineering: Reshaping the biocatalytic functions. Biotechnol Bioeng 2020; 117:1877-1894. [DOI: 10.1002/bit.27329] [Citation(s) in RCA: 40] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/18/2019] [Revised: 01/13/2020] [Accepted: 03/09/2020] [Indexed: 12/19/2022]
Affiliation(s)
- Misha Ali
- Department of Biochemistry, Faculty of Life SciencesAligarh Muslim University Aligarh Uttar Pradesh India
| | | | - Qayyum Husain
- Department of Biochemistry, Faculty of Life SciencesAligarh Muslim University Aligarh Uttar Pradesh India
| |
Collapse
|
25
|
Hayes RL, Vilseck JZ, Brooks CL. Approaching protein design with multisite λ dynamics: Accurate and scalable mutational folding free energies in T4 lysozyme. Protein Sci 2019; 27:1910-1922. [PMID: 30175503 DOI: 10.1002/pro.3500] [Citation(s) in RCA: 23] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/17/2018] [Revised: 08/06/2018] [Accepted: 08/15/2018] [Indexed: 12/14/2022]
Abstract
The estimation of changes in free energy upon mutation is central to the problem of protein design. Modern protein design methods have had remarkable success over a wide range of design targets, but are reaching their limits in ligand binding and enzyme design due to insufficient accuracy in mutational free energies. Alchemical free energy calculations have the potential to supplement modern design methods through more accurate molecular dynamics based prediction of free energy changes, but suffer from high computational cost. Multisite λ dynamics (MSλD) is a particularly efficient and scalable free energy method with potential to explore combinatorially large sequence spaces inaccessible with other free energy methods. This work aims to quantify the accuracy of MSλD and demonstrate its scalability. We apply MSλD to the classic problem of calculating folding free energies in T4 lysozyme, a system with a wealth of experimental measurements. Single site mutants considering 32 mutations show remarkable agreement with experiment with a Pearson correlation of 0.914 and mean unsigned error of 1.19 kcal/mol. Multisite mutants in systems with up to five concurrent mutations spanning 240 different sequences show comparable agreement with experiment. These results demonstrate the promise of MSλD in exploring large sequence spaces for protein design.
Collapse
Affiliation(s)
- Ryan L Hayes
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan, 48109
| | - Jonah Z Vilseck
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan, 48109
| | - Charles L Brooks
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan, 48109.,Biophysics Program, University of Michigan, Ann Arbor, Michigan, 48109
| |
Collapse
|
26
|
Learning the pattern of epistasis linking genotype and phenotype in a protein. Nat Commun 2019; 10:4213. [PMID: 31527666 PMCID: PMC6746860 DOI: 10.1038/s41467-019-12130-8] [Citation(s) in RCA: 73] [Impact Index Per Article: 14.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/21/2018] [Accepted: 08/19/2019] [Indexed: 02/02/2023] Open
Abstract
Understanding the pattern of epistasis-the non-independence of mutations-is critical for relating genotype and phenotype. However, the combinatorial complexity of potential epistatic interactions has severely limited the analysis of this problem. Using new mutational approaches, we report a comprehensive experimental study of all 213 mutants that link two phenotypically distinct variants of the Entacmaea quadricolor fluorescent protein-an opportunity to examine epistasis up to the 13th order. The data show the existence of many high-order epistatic interactions between mutations, but also reveal extraordinary sparsity, enabling novel experimental and computational strategies for learning the relevant epistasis. We demonstrate that such information, in turn, can be used to accurately predict phenotypes in practical situations where the number of measurements is limited. Finally, we show how the observed epistasis shapes the solution space of single-mutation trajectories between the parental fluorescent proteins, informative about the protein's evolutionary potential. This work provides conceptual and experimental strategies to profoundly characterize epistasis in a protein, relevant to both natural and laboratory evolution.
Collapse
|
27
|
Abstract
The application of native enzymes may not be economical owing to the stability factor. A smaller protein molecule may be less susceptible to external stresses. Haloalkane dehalogenases (HLDs) that act on toxic haloalkanes may be incorporated as bioreceptors to detect haloalkane contaminants. Therefore, this study aims to develop mini proteins of HLD as an alternative bioreceptor which was able to withstand extreme conditions. Initially, the mini proteins were designed through computer modeling. Based on the results, five designed mini proteins were deemed to be viable stable mini proteins. They were then validated through experimental study. The smallest mini protein (model 5) was chosen for subsequent analysis as it was expressed in soluble form. No dehalogenase activity was detected, thus the specific binding interaction of between 1,3-dibromopropane with mini protein was investigated using isothermal titration calorimetry. Higher binding affinity between 1,3-dibromopropane and mini protein was obtained than the native. Thermal stability study with circular dichroism had proven that the mini protein possessed two times higher Tm value at 83.73 °C than the native at 43.97 °C. In conclusion, a stable mini protein was successfully designed and may be used as bioreceptors in the haloalkane sensor that is suitable for industrial application.
Collapse
|
28
|
Gutte B, Klauser S. Design of catalytic polypeptides and proteins. Protein Eng Des Sel 2018; 31:457-470. [PMID: 31241746 DOI: 10.1093/protein/gzz009] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/14/2019] [Indexed: 11/13/2022] Open
Abstract
The first part of this review article lists examples of complete, empirical de novo design that made important contributions to the development of the field and initiated challenging projects. The second part of this article deals with computational design of novel enzymes in native protein scaffolds; active designs were refined through random and site-directed mutagenesis producing artificial enzymes with nearly native enzyme- like activities against a number of non-natural substrates. Combining aspects of de novo design and biological evolution of nature's enzymes has started and will accelerate the development of novel enzyme activities.
Collapse
Affiliation(s)
- B Gutte
- Department of Biochemistry, University of Zürich, Winterthurerstrasse 190, Zürich, Switzerland
| | - S Klauser
- Department of Biochemistry, University of Zürich, Winterthurerstrasse 190, Zürich, Switzerland
| |
Collapse
|
29
|
Abstract
Directed evolution is a powerful technique for generating tailor-made enzymes for a wide range of biocatalytic applications. Following the principles of natural evolution, iterative cycles of mutagenesis and screening or selection are applied to modify protein properties, enhance catalytic activities, or develop completely new protein catalysts for non-natural chemical transformations. This review briefly surveys the experimental methods used to generate genetic diversity and screen or select for improved enzyme variants. Emphasis is placed on a key challenge, namely how to generate novel catalytic activities that expand the scope of natural reactions. Two particularly effective strategies, exploiting catalytic promiscuity and rational design, are illustrated by representative examples of successfully evolved enzymes. Opportunities for extending these approaches to more complex biocatalytic systems are also considered.
Collapse
Affiliation(s)
- Cathleen Zeymer
- Laboratory of Organic Chemistry, ETH Zürich, CH-8093 Zürich, Switzerland;,
| | - Donald Hilvert
- Laboratory of Organic Chemistry, ETH Zürich, CH-8093 Zürich, Switzerland;,
| |
Collapse
|
30
|
Gainza-Cirauqui P, Correia BE. Computational protein design-the next generation tool to expand synthetic biology applications. Curr Opin Biotechnol 2018; 52:145-152. [PMID: 29729544 DOI: 10.1016/j.copbio.2018.04.001] [Citation(s) in RCA: 26] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/17/2017] [Revised: 04/03/2018] [Accepted: 04/04/2018] [Indexed: 11/25/2022]
Abstract
One powerful approach to engineer synthetic biology pathways is the assembly of proteins sourced from one or more natural organisms. However, synthetic pathways often require custom functions or biophysical properties not displayed by natural proteins, limitations that could be overcome through modern protein engineering techniques. Structure-based computational protein design is a powerful tool to engineer new functional capabilities in proteins, and it is beginning to have a profound impact in synthetic biology. Here, we review efforts to increase the capabilities of synthetic biology using computational protein design. We focus primarily on computationally designed proteins not only validated in vitro, but also shown to modulate different activities in living cells. Efforts made to validate computational designs in cells can illustrate both the challenges and opportunities in the intersection of protein design and synthetic biology. We also highlight protein design approaches, which although not validated as conveyors of new cellular function in situ, may have rapid and innovative applications in synthetic biology. We foresee that in the near-future, computational protein design will vastly expand the functional capabilities of synthetic cells.
Collapse
Affiliation(s)
- Pablo Gainza-Cirauqui
- Institute of Bioengineering, École Polytechnique Fédérale de Lausanne, Lausanne CH-1015, Switzerland; Swiss Institute of Bioinformatics (SIB), Lausanne CH-1015, Switzerland
| | - Bruno Emanuel Correia
- Institute of Bioengineering, École Polytechnique Fédérale de Lausanne, Lausanne CH-1015, Switzerland; Swiss Institute of Bioinformatics (SIB), Lausanne CH-1015, Switzerland.
| |
Collapse
|
31
|
Enabling tools for high-throughput detection of metabolites: Metabolic engineering and directed evolution applications. Biotechnol Adv 2017; 35:950-970. [PMID: 28723577 DOI: 10.1016/j.biotechadv.2017.07.005] [Citation(s) in RCA: 79] [Impact Index Per Article: 11.3] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/16/2017] [Revised: 06/07/2017] [Accepted: 07/11/2017] [Indexed: 12/21/2022]
Abstract
Within the Design-Build-Test Cycle for strain engineering, rapid product detection and selection strategies remain challenging and limit overall throughput. Here we summarize a wide variety of modalities that transduce chemical concentrations into easily measured absorbance, luminescence, and fluorescence signals. Specifically, we cover protein-based biosensors (including transcription factors), nucleic acid-based biosensors, coupled enzyme reactions, bioorthogonal chemistry, and fluorescent and chromogenic dyes and substrates as modalities for detection. We focus on the use of these methods for strain engineering and enzyme discovery and conclude with remarks on the current and future state of biosensor development for application in the metabolic engineering field.
Collapse
|
32
|
Computational design of ligand-binding proteins. Curr Opin Struct Biol 2016; 45:67-73. [PMID: 27951448 DOI: 10.1016/j.sbi.2016.11.021] [Citation(s) in RCA: 27] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/16/2016] [Revised: 11/24/2016] [Accepted: 11/25/2016] [Indexed: 02/07/2023]
Abstract
Custom-designed ligand-binding proteins with novel functions hold the potential for numerous applications. In recent years, the developments of computational methods together with high-throughput experimental screening techniques have led to the generation of novel, high-affinity ligand-binding proteins for given ligands. In addition, naturally occurring ligand-binding proteins have been computationally designed to recognize new ligands while keeping their original biological functions at the same time. Furthermore, metalloproteins have been successfully designed for novel functions and applications. Though much has been learned in these successful design cases, advances in our understanding of protein dynamics and functions related to ligand binding and development of novel computational strategies are necessary to further increase the success rate of computational protein-ligand binding design.
Collapse
|
33
|
De Paepe B, Peters G, Coussement P, Maertens J, De Mey M. Tailor-made transcriptional biosensors for optimizing microbial cell factories. J Ind Microbiol Biotechnol 2016; 44:623-645. [PMID: 27837353 DOI: 10.1007/s10295-016-1862-3] [Citation(s) in RCA: 68] [Impact Index Per Article: 8.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/25/2016] [Accepted: 10/30/2016] [Indexed: 12/24/2022]
Abstract
Monitoring cellular behavior and eventually properly adapting cellular processes is key to handle the enormous complexity of today's metabolic engineering questions. Hence, transcriptional biosensors bear the potential to augment and accelerate current metabolic engineering strategies, catalyzing vital advances in industrial biotechnology. The development of such transcriptional biosensors typically starts with exploring nature's richness. Hence, in a first part, the transcriptional biosensor architecture and the various modi operandi are briefly discussed, as well as experimental and computational methods and relevant ontologies to search for natural transcription factors and their corresponding binding sites. In the second part of this review, various engineering approaches are reviewed to tune the main characteristics of these (natural) transcriptional biosensors, i.e., the response curve and ligand specificity, in view of specific industrial biotechnology applications, which is illustrated using success stories of transcriptional biosensor engineering.
Collapse
Affiliation(s)
- Brecht De Paepe
- Department of Biochemical and Microbial Technology, Ghent University, Coupure Links 653, 9000, Ghent, Belgium
| | - Gert Peters
- Department of Biochemical and Microbial Technology, Ghent University, Coupure Links 653, 9000, Ghent, Belgium
| | - Pieter Coussement
- Department of Biochemical and Microbial Technology, Ghent University, Coupure Links 653, 9000, Ghent, Belgium
| | - Jo Maertens
- Department of Biochemical and Microbial Technology, Ghent University, Coupure Links 653, 9000, Ghent, Belgium
| | - Marjan De Mey
- Department of Biochemical and Microbial Technology, Ghent University, Coupure Links 653, 9000, Ghent, Belgium.
| |
Collapse
|
34
|
Obexer R, Pott M, Zeymer C, Griffiths AD, Hilvert D. Efficient laboratory evolution of computationally designed enzymes with low starting activities using fluorescence-activated droplet sorting. Protein Eng Des Sel 2016; 29:355-66. [DOI: 10.1093/protein/gzw032] [Citation(s) in RCA: 47] [Impact Index Per Article: 5.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/19/2016] [Accepted: 06/21/2016] [Indexed: 11/14/2022] Open
|
35
|
Huang X, Xue J, Lin M, Zhu Y. Use of an Improved Matching Algorithm to Select Scaffolds for Enzyme Design Based on a Complex Active Site Model. PLoS One 2016; 11:e0156559. [PMID: 27243223 PMCID: PMC4887040 DOI: 10.1371/journal.pone.0156559] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/21/2016] [Accepted: 05/16/2016] [Indexed: 11/24/2022] Open
Abstract
Active site preorganization helps native enzymes electrostatically stabilize the transition state better than the ground state for their primary substrates and achieve significant rate enhancement. In this report, we hypothesize that a complex active site model for active site preorganization modeling should help to create preorganized active site design and afford higher starting activities towards target reactions. Our matching algorithm ProdaMatch was improved by invoking effective pruning strategies and the native active sites for ten scaffolds in a benchmark test set were reproduced. The root-mean squared deviations between the matched transition states and those in the crystal structures were < 1.0 Å for the ten scaffolds, and the repacking calculation results showed that 91% of the hydrogen bonds within the active sites are recovered, indicating that the active sites can be preorganized based on the predicted positions of transition states. The application of the complex active site model for de novo enzyme design was evaluated by scaffold selection using a classic catalytic triad motif for the hydrolysis of p-nitrophenyl acetate. Eighty scaffolds were identified from a scaffold library with 1,491 proteins and four scaffolds were native esterase. Furthermore, enzyme design for complicated substrates was investigated for the hydrolysis of cephalexin using scaffold selection based on two different catalytic motifs. Only three scaffolds were identified from the scaffold library by virtue of the classic catalytic triad-based motif. In contrast, 40 scaffolds were identified using a more flexible, but still preorganized catalytic motif, where one scaffold corresponded to the α-amino acid ester hydrolase that catalyzes the hydrolysis and synthesis of cephalexin. Thus, the complex active site modeling approach for de novo enzyme design with the aid of the improved ProdaMatch program is a promising approach for the creation of active sites with high catalytic efficiencies towards target reactions.
Collapse
Affiliation(s)
- Xiaoqiang Huang
- Department of Chemical Engineering, Tsinghua University, Beijing 100084, P.R. China
| | - Jing Xue
- Department of Chemical Engineering, Tsinghua University, Beijing 100084, P.R. China
| | - Min Lin
- Department of Chemical Engineering, Tsinghua University, Beijing 100084, P.R. China
| | - Yushan Zhu
- Department of Chemical Engineering, Tsinghua University, Beijing 100084, P.R. China
- * E-mail:
| |
Collapse
|
36
|
de los Santos ELC, Meyerowitz JT, Mayo SL, Murray RM. Engineering Transcriptional Regulator Effector Specificity Using Computational Design and In Vitro Rapid Prototyping: Developing a Vanillin Sensor. ACS Synth Biol 2016; 5:287-95. [PMID: 26262913 DOI: 10.1021/acssynbio.5b00090] [Citation(s) in RCA: 59] [Impact Index Per Article: 7.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
The pursuit of circuits and metabolic pathways of increasing complexity and robustness in synthetic biology will require engineering new regulatory tools. Feedback control based on relevant molecules, including toxic intermediates and environmental signals, would enable genetic circuits to react appropriately to changing conditions. In this work, variants of qacR, a tetR family repressor, were generated by computational protein design and screened in a cell-free transcription-translation (TX-TL) system for responsiveness to a new targeted effector. The modified repressors target vanillin, a growth-inhibiting small molecule found in lignocellulosic hydrolysates and other industrial processes. Promising candidates from the in vitro screen were further characterized in vitro and in vivo in a gene circuit. The screen yielded two qacR mutants that respond to vanillin both in vitro and in vivo. While the mutants exhibit some toxicity to cells, presumably due to off-target effects, they are prime starting points for directed evolution toward vanillin sensors with the specifications required for use in a dynamic control loop. We believe this process, a combination of the generation of variants coupled with in vitro screening, can serve as a framework for designing new sensors for other target compounds.
Collapse
Affiliation(s)
- Emmanuel L. C. de los Santos
- Division of Biology and Biological
Engineering, California Institute of Technology, Pasadena, California 91125, United States
| | - Joseph T. Meyerowitz
- Division of Biology and Biological
Engineering, California Institute of Technology, Pasadena, California 91125, United States
| | - Stephen L. Mayo
- Division of Biology and Biological
Engineering, California Institute of Technology, Pasadena, California 91125, United States
| | - Richard M. Murray
- Division of Biology and Biological
Engineering, California Institute of Technology, Pasadena, California 91125, United States
| |
Collapse
|
37
|
A fast loop-closure algorithm to accelerate residue matching in computational enzyme design. J Mol Model 2016; 22:49. [DOI: 10.1007/s00894-016-2915-2] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/02/2015] [Accepted: 01/11/2016] [Indexed: 01/04/2023]
|
38
|
Świderek K, Moliner V. Computational Studies of Candida Antarctica Lipase B to Test Its Capability as a Starting Point To Redesign New Diels-Alderases. J Phys Chem B 2015; 120:2053-70. [PMID: 26624234 DOI: 10.1021/acs.jpcb.5b10527] [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/28/2022]
Abstract
The design of new biocatalysts is a target that is receiving increasing attention. One of the most popular reactions in this regard is the Diels-Alder cycloaddition because of its applications in organic synthesis and the absence of efficient natural enzymes that catalyze it. In this paper, the possibilities of using the highly promiscuous Candida Antarctica lipase B as a protein scaffold to redesign a Diels-Alderase has been explored by means of theoretical quantum mechanics/molecular mechanics (QM/MM) molecular dynamics simulations. Free energy surfaces have been computed for two reactions in the wild-type and in several mutants with hybrid AM1/MM potentials with corrections at M06-2X/MM level. The study of the counterpart reactions in solution has allowed performing comparative analysis that render interesting conclusions. Since the dienophile anchors very well in the oxyanion hole of all tested protein variants, the slight electronic changes from reactant complex to the transition state suggest that mutations should be focused in favoring the formation of reactive conformations of a reactant complex that, in turn, would reduce the energy barrier.
Collapse
Affiliation(s)
- Katarzyna Świderek
- Institute of Applied Radiation Chemistry, Lodz University of Technology , 90-924 Lodz, Poland.,Departamento de Química Física y Analítica, Universitat Jaume I , 12071 Castellón, Spain
| | - Vicent Moliner
- Departamento de Química Física y Analítica, Universitat Jaume I , 12071 Castellón, Spain
| |
Collapse
|
39
|
Druart K, Palmai Z, Omarjee E, Simonson T. Protein:Ligand binding free energies: A stringent test for computational protein design. J Comput Chem 2015; 37:404-15. [PMID: 26503829 DOI: 10.1002/jcc.24230] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/01/2015] [Revised: 10/01/2015] [Accepted: 10/02/2015] [Indexed: 01/29/2023]
Abstract
A computational protein design method is extended to allow Monte Carlo simulations where two ligands are titrated into a protein binding pocket, yielding binding free energy differences. These provide a stringent test of the physical model, including the energy surface and sidechain rotamer definition. As a test, we consider tyrosyl-tRNA synthetase (TyrRS), which has been extensively redesigned experimentally. We consider its specificity for its substrate l-tyrosine (l-Tyr), compared to the analogs d-Tyr, p-acetyl-, and p-azido-phenylalanine (ac-Phe, az-Phe). We simulate l- and d-Tyr binding to TyrRS and six mutants, and compare the structures and binding free energies to a more rigorous "MD/GBSA" procedure: molecular dynamics with explicit solvent for structures and a Generalized Born + Surface Area model for binding free energies. Next, we consider l-Tyr, ac- and az-Phe binding to six other TyrRS variants. The titration results are sensitive to the precise rotamer definition, which involves a short energy minimization for each sidechain pair to help relax bad contacts induced by the discrete rotamer set. However, when designed mutant structures are rescored with a standard GBSA energy model, results agree well with the more rigorous MD/GBSA. As a third test, we redesign three amino acid positions in the substrate coordination sphere, with either l-Tyr or d-Tyr as the ligand. For two, we obtain good agreement with experiment, recovering the wildtype residue when l-Tyr is the ligand and a d-Tyr specific mutant when d-Tyr is the ligand. For the third, we recover His with either ligand, instead of wildtype Gln.
Collapse
Affiliation(s)
- Karen Druart
- Laboratoire De Biochimie (UMR CNRS 7654), Department of Biology, Ecole Polytechnique, Palaiseau, France
| | - Zoltan Palmai
- Laboratoire De Biochimie (UMR CNRS 7654), Department of Biology, Ecole Polytechnique, Palaiseau, France
| | - Eyaz Omarjee
- Laboratoire De Biochimie (UMR CNRS 7654), Department of Biology, Ecole Polytechnique, Palaiseau, France
| | - Thomas Simonson
- Laboratoire De Biochimie (UMR CNRS 7654), Department of Biology, Ecole Polytechnique, Palaiseau, France
| |
Collapse
|
40
|
Ollikainen N, de Jong RM, Kortemme T. Coupling Protein Side-Chain and Backbone Flexibility Improves the Re-design of Protein-Ligand Specificity. PLoS Comput Biol 2015; 11:e1004335. [PMID: 26397464 PMCID: PMC4580623 DOI: 10.1371/journal.pcbi.1004335] [Citation(s) in RCA: 62] [Impact Index Per Article: 6.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/10/2015] [Accepted: 05/10/2015] [Indexed: 11/25/2022] Open
Abstract
Interactions between small molecules and proteins play critical roles in regulating and facilitating diverse biological functions, yet our ability to accurately re-engineer the specificity of these interactions using computational approaches has been limited. One main difficulty, in addition to inaccuracies in energy functions, is the exquisite sensitivity of protein–ligand interactions to subtle conformational changes, coupled with the computational problem of sampling the large conformational search space of degrees of freedom of ligands, amino acid side chains, and the protein backbone. Here, we describe two benchmarks for evaluating the accuracy of computational approaches for re-engineering protein-ligand interactions: (i) prediction of enzyme specificity altering mutations and (ii) prediction of sequence tolerance in ligand binding sites. After finding that current state-of-the-art “fixed backbone” design methods perform poorly on these tests, we develop a new “coupled moves” design method in the program Rosetta that couples changes to protein sequence with alterations in both protein side-chain and protein backbone conformations, and allows for changes in ligand rigid-body and torsion degrees of freedom. We show significantly increased accuracy in both predicting ligand specificity altering mutations and binding site sequences. These methodological improvements should be useful for many applications of protein – ligand design. The approach also provides insights into the role of subtle conformational adjustments that enable functional changes not only in engineering applications but also in natural protein evolution. Designing new protein–ligand interactions has tremendous potential for engineering sensitive biosensors for diagnostics or new enzymes useful in biotechnology, but these applications are extremely challenging, both because of inaccuracies of the energy functions used in modeling and design, and because protein active and binding sites are highly sensitive to subtle changes in structure. Here we describe a new method that addresses the second problem and couples changes in the structure of the protein backbone and of the amino acid side chains, the amino acid sequence, and the conformation of the ligand and its orientation in the binding site. We show that our method improvements significantly increase the accuracy of designing protein–ligand interactions compared to current state-of-the-art design methods. We assess these improvements in two important tests: the first predicts mutations that change ligand-binding preferences in enzymes, and the second predicts protein sequences that bind a given ligand. In these tests, subtle conformational changes made in our model are essential to recapitulate both the results from engineering experiments and the sequence diversity occurring in natural protein families. These results therefore shed light on the mechanisms of how new protein functions might have emerged and can be engineered in the laboratory.
Collapse
Affiliation(s)
- Noah Ollikainen
- Graduate Program in Bioinformatics, University of California San Francisco, San Francisco, California, United States of America
| | - René M. de Jong
- DSM Biotechnology Center, Alexander Fleminglaan 1, Delft, The Netherlands
| | - Tanja Kortemme
- Graduate Program in Bioinformatics, University of California San Francisco, San Francisco, California, United States of America
- California Institute for Quantitative Biosciences (QB3), University of California San Francisco, San Francisco, California, United States of America
- Department of Bioengineering and Therapeutic Science, University of California San Francisco, San Francisco, California, United States of America
- * E-mail:
| |
Collapse
|
41
|
Ś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
|
42
|
Tian Y, Huang X, Zhu Y. Computational design of enzyme-ligand binding using a combined energy function and deterministic sequence optimization algorithm. J Mol Model 2015; 21:191. [PMID: 26162695 DOI: 10.1007/s00894-015-2742-x] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/27/2015] [Accepted: 06/24/2015] [Indexed: 01/06/2023]
Abstract
Enzyme amino-acid sequences at ligand-binding interfaces are evolutionarily optimized for reactions, and the natural conformation of an enzyme-ligand complex must have a low free energy relative to alternative conformations in native-like or non-native sequences. Based on this assumption, a combined energy function was developed for enzyme design and then evaluated by recapitulating native enzyme sequences at ligand-binding interfaces for 10 enzyme-ligand complexes. In this energy function, the electrostatic interaction between polar or charged atoms at buried interfaces is described by an explicitly orientation-dependent hydrogen-bonding potential and a pairwise-decomposable generalized Born model based on the general side chain in the protein design framework. The energy function is augmented with a pairwise surface-area based hydrophobic contribution for nonpolar atom burial. Using this function, on average, 78% of the amino acids at ligand-binding sites were predicted correctly in the minimum-energy sequences, whereas 84% were predicted correctly in the most-similar sequences, which were selected from the top 20 sequences for each enzyme-ligand complex. Hydrogen bonds at the enzyme-ligand binding interfaces in the 10 complexes were usually recovered with the correct geometries. The binding energies calculated using the combined energy function helped to discriminate the active sequences from a pool of alternative sequences that were generated by repeatedly solving a series of mixed-integer linear programming problems for sequence selection with increasing integer cuts.
Collapse
Affiliation(s)
- Ye Tian
- Department of Chemical Engineering, Tsinghua University, Beijing, 100084, People's Republic of China
| | | | | |
Collapse
|
43
|
Mahapa A, Mandal S, Biswas A, Jana B, Polley S, Sau S, Sau K. Chemical and thermal unfolding of a global staphylococcal virulence regulator with a flexible C-terminal end. PLoS One 2015; 10:e0122168. [PMID: 25822635 PMCID: PMC4379015 DOI: 10.1371/journal.pone.0122168] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/20/2014] [Accepted: 02/07/2015] [Indexed: 11/19/2022] Open
Abstract
SarA, a Staphylococcus aureus-specific dimeric protein, modulates the expression of numerous proteins including various virulence factors. Interestingly, S. aureus synthesizes multiple SarA paralogs seemingly for optimizing the expression of its virulence factors. To understand the domain structure/flexibility and the folding/unfolding mechanism of the SarA protein family, we have studied a recombinant SarA (designated rSarA) using various in vitro probes. Limited proteolysis of rSarA and the subsequent analysis of the resulting protein fragments suggested it to be a single-domain protein with a long, flexible C-terminal end. rSarA was unfolded by different mechanisms in the presence of different chemical and physical denaturants. While urea-induced unfolding of rSarA occurred successively via the formation of a dimeric and a monomeric intermediate, GdnCl-induced unfolding of this protein proceeded through the production of two dimeric intermediates. The surface hydrophobicity and the structures of the intermediates were not identical and also differed significantly from those of native rSarA. Of the intermediates, the GdnCl-generated intermediates not only possessed a molten globule-like structure but also exhibited resistance to dissociation during their unfolding. Compared to the native rSarA, the intermediate that was originated at lower GdnCl concentration carried a compact shape, whereas, other intermediates owned a swelled shape. The chemical-induced unfolding, unlike thermal unfolding of rSarA, was completely reversible in nature.
Collapse
Affiliation(s)
- Avisek Mahapa
- Department of Biotechnology, Haldia Institute of Technology, Haldia, West Bengal, India
| | - Sukhendu Mandal
- Department of Biochemistry, Bose Institute, Kolkata, West Bengal, India
| | - Anindya Biswas
- Department of Biochemistry, Bose Institute, Kolkata, West Bengal, India
| | - Biswanath Jana
- Department of Biochemistry, Bose Institute, Kolkata, West Bengal, India
| | - Soumitra Polley
- Department of Biochemistry, Bose Institute, Kolkata, West Bengal, India
| | - Subrata Sau
- Department of Biochemistry, Bose Institute, Kolkata, West Bengal, India
- * E-mail: (SS); (KS)
| | - Keya Sau
- Department of Biotechnology, Haldia Institute of Technology, Haldia, West Bengal, India
- * E-mail: (SS); (KS)
| |
Collapse
|
44
|
Muñoz Robles V, Ortega-Carrasco E, Alonso-Cotchico L, Rodriguez-Guerra J, Lledós A, Maréchal JD. Toward the Computational Design of Artificial Metalloenzymes: From Protein–Ligand Docking to Multiscale Approaches. ACS Catal 2015. [DOI: 10.1021/acscatal.5b00010] [Citation(s) in RCA: 44] [Impact Index Per Article: 4.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Victor Muñoz Robles
- Departament de Química, Universitat Autònoma de Barcelona, 08193 Cerdanyola
del Vallès, Barcelona, Spain
| | - Elisabeth Ortega-Carrasco
- Departament de Química, Universitat Autònoma de Barcelona, 08193 Cerdanyola
del Vallès, Barcelona, Spain
| | - Lur Alonso-Cotchico
- Departament de Química, Universitat Autònoma de Barcelona, 08193 Cerdanyola
del Vallès, Barcelona, Spain
| | - Jaime Rodriguez-Guerra
- Departament de Química, Universitat Autònoma de Barcelona, 08193 Cerdanyola
del Vallès, Barcelona, Spain
| | - Agustí Lledós
- Departament de Química, Universitat Autònoma de Barcelona, 08193 Cerdanyola
del Vallès, Barcelona, Spain
| | - Jean-Didier Maréchal
- Departament de Química, Universitat Autònoma de Barcelona, 08193 Cerdanyola
del Vallès, Barcelona, Spain
| |
Collapse
|
45
|
Zanghellini A. de novo computational enzyme design. Curr Opin Biotechnol 2014; 29:132-8. [DOI: 10.1016/j.copbio.2014.03.002] [Citation(s) in RCA: 58] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/26/2014] [Accepted: 03/10/2014] [Indexed: 11/25/2022]
|
46
|
Subramaniam S, Senes A. Backbone dependency further improves side chain prediction efficiency in the Energy-based Conformer Library (bEBL). Proteins 2014; 82:3177-87. [PMID: 25212195 DOI: 10.1002/prot.24685] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/02/2014] [Revised: 08/21/2014] [Accepted: 09/03/2014] [Indexed: 12/11/2022]
Abstract
Side chain optimization is an integral component of many protein modeling applications. In these applications, the conformational freedom of the side chains is often explored using libraries of discrete, frequently occurring conformations. Because side chain optimization can pose a computationally intensive combinatorial problem, the nature of these conformer libraries is important for ensuring efficiency and accuracy in side chain prediction. We have previously developed an innovative method to create a conformer library with enhanced performance. The Energy-based Library (EBL) was obtained by analyzing the energetic interactions between conformers and a large number of natural protein environments from crystal structures. This process guided the selection of conformers with the highest propensity to fit into spaces that should accommodate a side chain. Because the method requires a large crystallographic data-set, the EBL was created in a backbone-independent fashion. However, it is well established that side chain conformation is strongly dependent on the local backbone geometry, and that backbone-dependent libraries are more efficient in side chain optimization. Here we present the backbone-dependent EBL (bEBL), whose conformers are independently sorted for each populated region of Ramachandran space. The resulting library closely mirrors the local backbone-dependent distribution of side chain conformation. Compared to the EBL, we demonstrate that the bEBL uses fewer conformers to produce similar side chain prediction outcomes, thus further improving performance with respect to the already efficient backbone-independent version of the library.
Collapse
Affiliation(s)
- Sabareesh Subramaniam
- Department of Biochemistry, University of Wisconsin-Madison, Madison, Wisconsin, 53706
| | | |
Collapse
|
47
|
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
|
48
|
Howell S, Inampudi K, Bean D, Wilson C. Understanding Thermal Adaptation of Enzymes through the Multistate Rational Design and Stability Prediction of 100 Adenylate Kinases. Structure 2014; 22:218-29. [DOI: 10.1016/j.str.2013.10.019] [Citation(s) in RCA: 23] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/22/2013] [Revised: 10/21/2013] [Accepted: 10/26/2013] [Indexed: 10/25/2022]
|
49
|
An accurate binding interaction model in de novo computational protein design of interactions: If you build it, they will bind. J Struct Biol 2014; 185:136-46. [DOI: 10.1016/j.jsb.2013.03.012] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/22/2012] [Revised: 03/15/2013] [Accepted: 03/21/2013] [Indexed: 01/07/2023]
|
50
|
Feldmeier K, Höcker B. Computational protein design of ligand binding and catalysis. Curr Opin Chem Biol 2013; 17:929-33. [DOI: 10.1016/j.cbpa.2013.10.002] [Citation(s) in RCA: 40] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
|