1
|
Tripathy RK, Khandave PY, Bzdrenga J, Nachon F, Brazzolotto X, Pande AH. Role of paraoxonase 1 in organophosphate G-series nerve agent poisoning and future therapeutic strategies. Arch Toxicol 2024:10.1007/s00204-024-03884-2. [PMID: 39356346 DOI: 10.1007/s00204-024-03884-2] [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: 07/12/2024] [Accepted: 09/24/2024] [Indexed: 10/03/2024]
Abstract
Chemical warfare nerve agents (CWNA) are neurotoxic chemicals unethically used as agents of mass destruction by terrorist outfits and during war. The available antidote against CWNA-mediated toxicity is not sufficiently effective and possesses several limitations. As a countermeasure, paraoxonase 1 (PON1), a catalytic bioscavenger, is being developed as a prophylactic treatment. However, the catalytic activity and substrate specificity of human PON1 are insufficient to be used as a potential antidote. Several laboratories have made different approaches to enhance the CWNA hydrolytic activity against various nerve agents. This review explores the holistic view of PON1 as a potential prophylactic agent against G-series CWNA poisoning, from its initial development to recent advancements and limitations. Apart from this, the review also provides an overview of all available PON1 variants that could be used as a potential prophylactic agent and discusses several possible ways to counteract immunogenicity.
Collapse
Affiliation(s)
- Rajan K Tripathy
- Department of Biotechnology, National Institute of Pharmaceutical Education and Research (NIPER), S.A.S. Nagar, Sector 67, Mohali, 160062, Punjab, India
| | - Prakash Y Khandave
- Department of Biotechnology, National Institute of Pharmaceutical Education and Research (NIPER), S.A.S. Nagar, Sector 67, Mohali, 160062, Punjab, India
| | - Janek Bzdrenga
- Département de Toxicologie et Risques Chimiques, Institut de Recherche Biomédicale des Armées, Place du Général Valérie André, 91220, Brétigny-Sur-Orge, France
| | - Florian Nachon
- Département de Toxicologie et Risques Chimiques, Institut de Recherche Biomédicale des Armées, Place du Général Valérie André, 91220, Brétigny-Sur-Orge, France
| | - Xavier Brazzolotto
- Département de Toxicologie et Risques Chimiques, Institut de Recherche Biomédicale des Armées, Place du Général Valérie André, 91220, Brétigny-Sur-Orge, France
| | - Abhay H Pande
- Department of Biotechnology, National Institute of Pharmaceutical Education and Research (NIPER), S.A.S. Nagar, Sector 67, Mohali, 160062, Punjab, India.
| |
Collapse
|
2
|
Lipsh-Sokolik R, Fleishman SJ. Addressing epistasis in the design of protein function. Proc Natl Acad Sci U S A 2024; 121:e2314999121. [PMID: 39133844 PMCID: PMC11348311 DOI: 10.1073/pnas.2314999121] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/29/2024] Open
Abstract
Mutations in protein active sites can dramatically improve function. The active site, however, is densely packed and extremely sensitive to mutations. Therefore, some mutations may only be tolerated in combination with others in a phenomenon known as epistasis. Epistasis reduces the likelihood of obtaining improved functional variants and dramatically slows natural and lab evolutionary processes. Research has shed light on the molecular origins of epistasis and its role in shaping evolutionary trajectories and outcomes. In addition, sequence- and AI-based strategies that infer epistatic relationships from mutational patterns in natural or experimental evolution data have been used to design functional protein variants. In recent years, combinations of such approaches and atomistic design calculations have successfully predicted highly functional combinatorial mutations in active sites. These were used to design thousands of functional active-site variants, demonstrating that, while our understanding of epistasis remains incomplete, some of the determinants that are critical for accurate design are now sufficiently understood. We conclude that the space of active-site variants that has been explored by evolution may be expanded dramatically to enhance natural activities or discover new ones. Furthermore, design opens the way to systematically exploring sequence and structure space and mutational impacts on function, deepening our understanding and control over protein activity.
Collapse
Affiliation(s)
- Rosalie Lipsh-Sokolik
- Department of Biomolecular Sciences, Weizmann Institute of Science, Rehovot 7610001, Israel
| | - Sarel J Fleishman
- Department of Biomolecular Sciences, Weizmann Institute of Science, Rehovot 7610001, Israel
| |
Collapse
|
3
|
Chen A, Zhang XD, Đelmaš AĐ, Weitz DA, Milcic K. Systems and Methods for Continuous Evolution of Enzymes. Chemistry 2024; 30:e202400880. [PMID: 38780896 DOI: 10.1002/chem.202400880] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/01/2024] [Revised: 05/21/2024] [Accepted: 05/23/2024] [Indexed: 05/25/2024]
Abstract
Directed evolution generates novel biomolecules with desired functions by iteratively diversifying the genetic sequence of wildtype biomolecules, relaying the genetic information to the molecule with function, and selecting the variants that progresses towards the properties of interest. While traditional directed evolution consumes significant labor and time for each step, continuous evolution seeks to automate all steps so directed evolution can proceed with minimum human intervention and dramatically shortened time. A major application of continuous evolution is the generation of novel enzymes, which catalyze reactions under conditions that are not favorable to their wildtype counterparts, or on altered substrates. The challenge to continuously evolve enzymes lies in automating sufficient, unbiased gene diversification, providing selection for a wide array of reaction types, and linking the genetic information to the phenotypic function. Over years of development, continuous evolution has accumulated versatile strategies to address these challenges, enabling its use as a general tool for enzyme engineering. As the capability of continuous evolution continues to expand, its impact will increase across various industries. In this review, we summarize the working mechanisms of recently developed continuous evolution strategies, discuss examples of their applications focusing on enzyme evolution, and point out their limitations and future directions.
Collapse
Affiliation(s)
- Anqi Chen
- School of Engineering and Applied Sciences, Harvard University, Cambridge, Massachusetts, 02138, USA E-mail: Dr David A. Weitz: E-mail: Dr. Karla Milcic
| | - Xinge Diana Zhang
- School of Engineering and Applied Sciences, Harvard University, Cambridge, Massachusetts, 02138, USA E-mail: Dr David A. Weitz: E-mail: Dr. Karla Milcic
| | | | - David A Weitz
- School of Engineering and Applied Sciences, Harvard University, Cambridge, Massachusetts, 02138, USA E-mail: Dr David A. Weitz: E-mail: Dr. Karla Milcic
- Wyss Institute for Biologically Inspired Engineering, Harvard University, 3 Blackfan Circle, Boston, MA, 02115, USA
- Department of Physics, Harvard University, Cambridge, MA, 02138, USA
| | - Karla Milcic
- School of Engineering and Applied Sciences, Harvard University, Cambridge, Massachusetts, 02138, USA E-mail: Dr David A. Weitz: E-mail: Dr. Karla Milcic
- University of Belgrade-Faculty of Chemistry, Studentski trg 12-16, 11000, Belgrade, Serbia
| |
Collapse
|
4
|
Wardman JF, Withers SG. Carbohydrate-active enzyme (CAZyme) discovery and engineering via (Ultra)high-throughput screening. RSC Chem Biol 2024; 5:595-616. [PMID: 38966674 PMCID: PMC11221537 DOI: 10.1039/d4cb00024b] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/23/2024] [Accepted: 05/16/2024] [Indexed: 07/06/2024] Open
Abstract
Carbohydrate-active enzymes (CAZymes) constitute a diverse set of enzymes that catalyze the assembly, degradation, and modification of carbohydrates. These enzymes have been fashioned into potent, selective catalysts by millennia of evolution, and yet are also highly adaptable and readily evolved in the laboratory. To identify and engineer CAZymes for different purposes, (ultra)high-throughput screening campaigns have been frequently utilized with great success. This review provides an overview of the different approaches taken in screening for CAZymes and how mechanistic understandings of CAZymes can enable new approaches to screening. Within, we also cover how cutting-edge techniques such as microfluidics, advances in computational approaches and synthetic biology, as well as novel assay designs are leading the field towards more informative and effective screening approaches.
Collapse
Affiliation(s)
- Jacob F Wardman
- Department of Biochemistry and Molecular Biology, University of British Columbia Vancouver BC V6T 1Z3 Canada
- Michael Smith Laboratories, University of British Columbia Vancouver BC V6T 1Z4 Canada
| | - Stephen G Withers
- Department of Biochemistry and Molecular Biology, University of British Columbia Vancouver BC V6T 1Z3 Canada
- Michael Smith Laboratories, University of British Columbia Vancouver BC V6T 1Z4 Canada
- Department of Chemistry, University of British Columbia Vancouver BC V6T 1Z1 Canada
| |
Collapse
|
5
|
Daffern N, Johansson KE, Baumer ZT, Robertson NR, Woojuh J, Bedewitz MA, Davis Z, Wheeldon I, Cutler SR, Lindorff-Larsen K, Whitehead TA. GMMA Can Stabilize Proteins Across Different Functional Constraints. J Mol Biol 2024; 436:168586. [PMID: 38663544 DOI: 10.1016/j.jmb.2024.168586] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/30/2024] [Revised: 04/16/2024] [Accepted: 04/17/2024] [Indexed: 05/06/2024]
Abstract
Stabilizing proteins without otherwise hampering their function is a central task in protein engineering and design. PYR1 is a plant hormone receptor that has been engineered to bind diverse small molecule ligands. We sought a set of generalized mutations that would provide stability without affecting functionality for PYR1 variants with diverse ligand-binding capabilities. To do this we used a global multi-mutant analysis (GMMA) approach, which can identify substitutions that have stabilizing effects and do not lower function. GMMA has the added benefit of finding substitutions that are stabilizing in different sequence contexts and we hypothesized that applying GMMA to PYR1 with different functionalities would identify this set of generalized mutations. Indeed, conducting FACS and deep sequencing of libraries for PYR1 variants with two different functionalities and applying a GMMA analysis identified 5 substitutions that, when inserted into four PYR1 variants that each bind a unique ligand, provided an increase of 2-6 °C in thermal inactivation temperature and no decrease in functionality.
Collapse
Affiliation(s)
- Nicolas Daffern
- Department of Chemical and Biological Engineering, University of Colorado Boulder, Boulder, CO 80305, USA
| | - Kristoffer E Johansson
- Linderstrøm-Lang Centre for Protein Science, Department of Biology, University of Copenhagen, Copenhagen, Denmark
| | - Zachary T Baumer
- Department of Chemical and Biological Engineering, University of Colorado Boulder, Boulder, CO 80305, USA
| | | | - Janty Woojuh
- Department of Botany and Plant Sciences, University of California, Riverside, USA
| | - Matthew A Bedewitz
- Department of Chemical and Biological Engineering, University of Colorado Boulder, Boulder, CO 80305, USA
| | - Zoë Davis
- Department of Chemical and Biological Engineering, University of Colorado Boulder, Boulder, CO 80305, USA
| | - Ian Wheeldon
- Department of Chemical and Environmental Engineering, University of California, Riverside, USA; Institute for Integrative Genome Biology, University of California, Riverside, Riverside, CA, USA
| | - Sean R Cutler
- Department of Botany and Plant Sciences, University of California, Riverside, USA; Institute for Integrative Genome Biology, University of California, Riverside, Riverside, CA, USA; Center for Plant Cell Biology, University of California, Riverside, Riverside, CA, USA
| | - Kresten Lindorff-Larsen
- Linderstrøm-Lang Centre for Protein Science, Department of Biology, University of Copenhagen, Copenhagen, Denmark.
| | - Timothy A Whitehead
- Department of Chemical and Biological Engineering, University of Colorado Boulder, Boulder, CO 80305, USA.
| |
Collapse
|
6
|
Carrera-Pacheco SE, Mueller A, Puente-Pineda JA, Zúñiga-Miranda J, Guamán LP. Designing cytochrome P450 enzymes for use in cancer gene therapy. Front Bioeng Biotechnol 2024; 12:1405466. [PMID: 38860140 PMCID: PMC11164052 DOI: 10.3389/fbioe.2024.1405466] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/22/2024] [Accepted: 04/30/2024] [Indexed: 06/12/2024] Open
Abstract
Cancer is a significant global socioeconomic burden, as millions of new cases and deaths occur annually. In 2020, almost 10 million cancer deaths were recorded worldwide. Advancements in cancer gene therapy have revolutionized the landscape of cancer treatment. An approach with promising potential for cancer gene therapy is introducing genes to cancer cells that encode for chemotherapy prodrug metabolizing enzymes, such as Cytochrome P450 (CYP) enzymes, which can contribute to the effective elimination of cancer cells. This can be achieved through gene-directed enzyme prodrug therapy (GDEPT). CYP enzymes can be genetically engineered to improve anticancer prodrug conversion to its active metabolites and to minimize chemotherapy side effects by reducing the prodrug dosage. Rational design, directed evolution, and phylogenetic methods are some approaches to developing tailored CYP enzymes for cancer therapy. Here, we provide a compilation of genetic modifications performed on CYP enzymes aiming to build highly efficient therapeutic genes capable of bio-activating different chemotherapeutic prodrugs. Additionally, this review summarizes promising preclinical and clinical trials highlighting engineered CYP enzymes' potential in GDEPT. Finally, the challenges, limitations, and future directions of using CYP enzymes for GDEPT in cancer gene therapy are discussed.
Collapse
Affiliation(s)
- Saskya E. Carrera-Pacheco
- Centro de Investigación Biomédica (CENBIO), Facultad de Ciencias de la Salud Eugenio Espejo, Universidad UTE, Quito, Ecuador
| | | | | | | | | |
Collapse
|
7
|
Sennett MA, Theobald DL. Extant Sequence Reconstruction: The Accuracy of Ancestral Sequence Reconstructions Evaluated by Extant Sequence Cross-Validation. J Mol Evol 2024; 92:181-206. [PMID: 38502220 PMCID: PMC10978691 DOI: 10.1007/s00239-024-10162-3] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/12/2023] [Accepted: 02/20/2024] [Indexed: 03/21/2024]
Abstract
Ancestral sequence reconstruction (ASR) is a phylogenetic method widely used to analyze the properties of ancient biomolecules and to elucidate mechanisms of molecular evolution. Despite its increasingly widespread application, the accuracy of ASR is currently unknown, as it is generally impossible to compare resurrected proteins to the true ancestors. Which evolutionary models are best for ASR? How accurate are the resulting inferences? Here we answer these questions using a cross-validation method to reconstruct each extant sequence in an alignment with ASR methodology, a method we term "extant sequence reconstruction" (ESR). We thus can evaluate the accuracy of ASR methodology by comparing ESR reconstructions to the corresponding known true sequences. We find that a common measure of the quality of a reconstructed sequence, the average probability, is indeed a good estimate of the fraction of correct amino acids when the evolutionary model is accurate or overparameterized. However, the average probability is a poor measure for comparing reconstructions from different models, because, surprisingly, a more accurate phylogenetic model often results in reconstructions with lower probability. While better (more predictive) models may produce reconstructions with lower sequence identity to the true sequences, better models nevertheless produce reconstructions that are more biophysically similar to true ancestors. In addition, we find that a large fraction of sequences sampled from the reconstruction distribution may have fewer errors than the single most probable (SMP) sequence reconstruction, despite the fact that the SMP has the lowest expected error of all possible sequences. Our results emphasize the importance of model selection for ASR and the usefulness of sampling sequence reconstructions for analyzing ancestral protein properties. ESR is a powerful method for validating the evolutionary models used for ASR and can be applied in practice to any phylogenetic analysis of real biological sequences. Most significantly, ESR uses ASR methodology to provide a general method by which the biophysical properties of resurrected proteins can be compared to the properties of the true protein.
Collapse
Affiliation(s)
- Michael A Sennett
- Department of Biochemistry, Brandeis University, Waltham, MA, 02453, USA
| | - Douglas L Theobald
- Department of Biochemistry, Brandeis University, Waltham, MA, 02453, USA.
| |
Collapse
|
8
|
Jian X, Li C, Feng X. Strategies for modulating transglycosylation activity, substrate specificity, and product polymerization degree of engineered transglycosylases. Crit Rev Biotechnol 2023; 43:1284-1298. [PMID: 36154438 DOI: 10.1080/07388551.2022.2105687] [Citation(s) in RCA: 7] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/23/2021] [Accepted: 06/21/2022] [Indexed: 01/18/2023]
Abstract
Glycosides are widely used in many fields due to their favorable biological activity. The traditional plant extractions and chemical methods for glycosides production are limited by environmentally unfriendly, laborious protecting group strategies and low yields. Alternatively, enzymatic glycosylation has drawn special attention due to its mild reaction conditions, high catalytic efficiency, and specific stereo-/regioselectivity. Glycosyltransferases (GTs) and retaining glycoside hydrolases (rGHs) are two major enzymes for the formation of glycosidic linkages. Therein GTs generally use nucleotide phosphate activated donors. In contrast, GHs can use broader simple and affordable glycosyl donors, showing great potential in industrial applications. However, most rGHs mainly show hydrolysis activity and only a few rGHs, namely non-Leloir transglycosylases (TGs), innately present strong transglycosylation activities. To address this problem, various strategies have recently been developed to successfully tailor rGHs to alleviate their hydrolysis activity and obtain the engineered TGs. This review summarizes the current modification strategies in TGs engineering, with a special focus on transglycosylation activity enhancement, substrate specificity modulation, and product polymerization degree distribution, which provides a reference for exploiting the transglycosylation potentials of rGHs.
Collapse
Affiliation(s)
- Xing Jian
- Key Laboratory of Medical Molecule Science and Pharmaceutics Engineering, Ministry of Industry and Information Technology, Institute of Biochemical Engineering, School of Chemistry and Chemical Engineering, Beijing Institute of Technology, Beijing, China
| | - Chun Li
- Key Laboratory of Medical Molecule Science and Pharmaceutics Engineering, Ministry of Industry and Information Technology, Institute of Biochemical Engineering, School of Chemistry and Chemical Engineering, Beijing Institute of Technology, Beijing, China
- Key Lab for Industrial Biocatalysis, Ministry of Education, Department of Chemical Engineering, Tsinghua University, Beijing, China
- Center for Synthetic & Systems Biology, Tsinghua University, Beijing, China
| | - Xudong Feng
- Key Laboratory of Medical Molecule Science and Pharmaceutics Engineering, Ministry of Industry and Information Technology, Institute of Biochemical Engineering, School of Chemistry and Chemical Engineering, Beijing Institute of Technology, Beijing, China
| |
Collapse
|
9
|
Buda K, Miton CM, Fan XC, Tokuriki N. Molecular determinants of protein evolvability. Trends Biochem Sci 2023; 48:751-760. [PMID: 37330341 DOI: 10.1016/j.tibs.2023.05.009] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/27/2023] [Revised: 05/18/2023] [Accepted: 05/23/2023] [Indexed: 06/19/2023]
Abstract
The plethora of biological functions that sustain life is rooted in the remarkable evolvability of proteins. An emerging view highlights the importance of a protein's initial state in dictating evolutionary success. A deeper comprehension of the mechanisms that govern the evolvability of these initial states can provide invaluable insights into protein evolution. In this review, we describe several molecular determinants of protein evolvability, unveiled by experimental evolution and ancestral sequence reconstruction studies. We further discuss how genetic variation and epistasis can promote or constrain functional innovation and suggest putative underlying mechanisms. By establishing a clear framework for these determinants, we provide potential indicators enabling the forecast of suitable evolutionary starting points and delineate molecular mechanisms in need of deeper exploration.
Collapse
Affiliation(s)
- Karol Buda
- Michael Smith Laboratories, University of British Columbia, Vancouver, Canada
| | - Charlotte M Miton
- Michael Smith Laboratories, University of British Columbia, Vancouver, Canada
| | - Xingyu Cara Fan
- Michael Smith Laboratories, University of British Columbia, Vancouver, Canada
| | - Nobuhiko Tokuriki
- Michael Smith Laboratories, University of British Columbia, Vancouver, Canada.
| |
Collapse
|
10
|
Liechty ET, Hren A, Kramer L, Donovan G, Friedman AJ, Shirts MR, Fox JM. Analysis of neutral mutational drift in an allosteric enzyme. Protein Sci 2023; 32:e4719. [PMID: 37402140 PMCID: PMC10364584 DOI: 10.1002/pro.4719] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/10/2023] [Revised: 06/26/2023] [Accepted: 06/28/2023] [Indexed: 07/05/2023]
Abstract
Neutral mutational drift is an important source of biological diversity that remains underexploited in fundamental studies of protein biophysics. This study uses a synthetic transcriptional circuit to study neutral drift in protein tyrosine phosphatase 1B (PTP1B), a mammalian signaling enzyme for which conformational changes are rate limiting. Kinetic assays of purified mutants indicate that catalytic activity, rather than thermodynamic stability, guides enrichment under neutral drift, where neutral or mildly activating mutations can mitigate the effects of deleterious ones. In general, mutants show a moderate activity-stability tradeoff, an indication that minor improvements in the activity of PTP1B do not require concomitant losses in its stability. Multiplexed sequencing of large mutant pools suggests that substitutions at allosterically influential sites are purged under biological selection, which enriches for mutations located outside of the active site. Findings indicate that the positional dependence of neutral mutations within drifting populations can reveal the presence of allosteric networks and illustrate an approach for using synthetic transcriptional systems to explore these mutations in regulatory enzymes.
Collapse
Affiliation(s)
- Evan T. Liechty
- Department of Chemical and Biological EngineeringUniversity of ColoradoBoulderColoradoUSA
| | - Andrew Hren
- Department of Chemical and Biological EngineeringUniversity of ColoradoBoulderColoradoUSA
| | - Levi Kramer
- Department of Chemical and Biological EngineeringUniversity of ColoradoBoulderColoradoUSA
| | - Gregory Donovan
- Department of Chemical and Biological EngineeringUniversity of ColoradoBoulderColoradoUSA
| | - Anika J. Friedman
- Department of Chemical and Biological EngineeringUniversity of ColoradoBoulderColoradoUSA
| | - Michael R. Shirts
- Department of Chemical and Biological EngineeringUniversity of ColoradoBoulderColoradoUSA
| | - Jerome M. Fox
- Department of Chemical and Biological EngineeringUniversity of ColoradoBoulderColoradoUSA
| |
Collapse
|
11
|
Seo D, Koh B, Eom GE, Kim HW, Kim S. A dual gene-specific mutator system installs all transition mutations at similar frequencies in vivo. Nucleic Acids Res 2023; 51:e59. [PMID: 37070179 PMCID: PMC10250238 DOI: 10.1093/nar/gkad266] [Citation(s) in RCA: 5] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/24/2022] [Accepted: 03/31/2023] [Indexed: 04/19/2023] Open
Abstract
Targeted in vivo hypermutation accelerates directed evolution of proteins through concurrent DNA diversification and selection. Although systems employing a fusion protein of a nucleobase deaminase and T7 RNA polymerase present gene-specific targeting, their mutational spectra have been limited to exclusive or dominant C:G→T:A mutations. Here we describe eMutaT7transition, a new gene-specific hypermutation system, that installs all transition mutations (C:G→T:A and A:T→G:C) at comparable frequencies. By using two mutator proteins in which two efficient deaminases, PmCDA1 and TadA-8e, are separately fused to T7 RNA polymerase, we obtained similar numbers of C:G→T:A and A:T→G:C substitutions at a sufficiently high frequency (∼6.7 substitutions in 1.3 kb gene during 80-h in vivo mutagenesis). Through eMutaT7transition-mediated TEM-1 evolution for antibiotic resistance, we generated many mutations found in clinical isolates. Overall, with a high mutation frequency and wider mutational spectrum, eMutaT7transition is a potential first-line method for gene-specific in vivo hypermutation.
Collapse
Affiliation(s)
- Daeje Seo
- Department of Chemistry, Seoul National University, 1 Gwanak-ro, Gwanak-gu, Seoul 08826, Republic of Korea
| | - Bonghyun Koh
- Department of Chemistry, Seoul National University, 1 Gwanak-ro, Gwanak-gu, Seoul 08826, Republic of Korea
| | - Ga-eul Eom
- Department of Chemistry, Seoul National University, 1 Gwanak-ro, Gwanak-gu, Seoul 08826, Republic of Korea
| | - Hye Won Kim
- Department of Chemistry, Seoul National University, 1 Gwanak-ro, Gwanak-gu, Seoul 08826, Republic of Korea
| | - Seokhee Kim
- Department of Chemistry, Seoul National University, 1 Gwanak-ro, Gwanak-gu, Seoul 08826, Republic of Korea
| |
Collapse
|
12
|
Sanejouand YH. On the Unknown Proteins of Eukaryotic Proteomes. J Mol Evol 2023:10.1007/s00239-023-10116-1. [PMID: 37219573 DOI: 10.1007/s00239-023-10116-1] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/30/2022] [Accepted: 05/07/2023] [Indexed: 05/24/2023]
Abstract
To study unknown proteins on a large scale, a reference system has been set up for the three better studied eukaryotic kingdoms, built with 36 proteomes as taxonomically diverse as possible. Proteins from 362 other eukaryotic proteomes with no known homologue in this set were then analyzed, focusing noteworthy on singletons, that is, on such proteins with no known homologue in their own proteome. Consistently, for a given species, no more than 12% of the singletons thus found are known at the protein level, according to Uniprot. In addition, since they rely on the information found in the alignment of homologous sequences, predictions of AlphaFold2 for their tridimensional structure are poor. In the case of metazoan species, the number of singletons rarely exceeds 1000 for the species the closest to the reference system (divergence times below 75 Myr). Interestingly, in the cases of viridiplantae and fungi, larger amounts of singletons are found for such species, as if the timescale on which singletons are added to proteomes were different in metazoa and in other eukaryotic kingdoms. In order to confirm this phenomenon, further studies of proteomes closer to those of the reference system are, however, needed.
Collapse
Affiliation(s)
- Yves-Henri Sanejouand
- US2B, UMR 6286 of CNRS, Nantes University, rue de la Houssinière, 44322, Nantes, France.
| |
Collapse
|
13
|
Gantz M, Neun S, Medcalf EJ, van Vliet LD, Hollfelder F. Ultrahigh-Throughput Enzyme Engineering and Discovery in In Vitro Compartments. Chem Rev 2023; 123:5571-5611. [PMID: 37126602 PMCID: PMC10176489 DOI: 10.1021/acs.chemrev.2c00910] [Citation(s) in RCA: 12] [Impact Index Per Article: 12.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/30/2022] [Indexed: 05/03/2023]
Abstract
Novel and improved biocatalysts are increasingly sourced from libraries via experimental screening. The success of such campaigns is crucially dependent on the number of candidates tested. Water-in-oil emulsion droplets can replace the classical test tube, to provide in vitro compartments as an alternative screening format, containing genotype and phenotype and enabling a readout of function. The scale-down to micrometer droplet diameters and picoliter volumes brings about a >107-fold volume reduction compared to 96-well-plate screening. Droplets made in automated microfluidic devices can be integrated into modular workflows to set up multistep screening protocols involving various detection modes to sort >107 variants a day with kHz frequencies. The repertoire of assays available for droplet screening covers all seven enzyme commission (EC) number classes, setting the stage for widespread use of droplet microfluidics in everyday biochemical experiments. We review the practicalities of adapting droplet screening for enzyme discovery and for detailed kinetic characterization. These new ways of working will not just accelerate discovery experiments currently limited by screening capacity but profoundly change the paradigms we can probe. By interfacing the results of ultrahigh-throughput droplet screening with next-generation sequencing and deep learning, strategies for directed evolution can be implemented, examined, and evaluated.
Collapse
Affiliation(s)
| | | | | | | | - Florian Hollfelder
- Department of Biochemistry, University of Cambridge, 80 Tennis Court Rd, Cambridge CB2 1GA, U.K.
| |
Collapse
|
14
|
Santos-Moreno J, Tasiudi E, Kusumawardhani H, Stelling J, Schaerli Y. Robustness and innovation in synthetic genotype networks. Nat Commun 2023; 14:2454. [PMID: 37117168 PMCID: PMC10147661 DOI: 10.1038/s41467-023-38033-3] [Citation(s) in RCA: 5] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/11/2022] [Accepted: 04/13/2023] [Indexed: 04/30/2023] Open
Abstract
Genotype networks are sets of genotypes connected by small mutational changes that share the same phenotype. They facilitate evolutionary innovation by enabling the exploration of different neighborhoods in genotype space. Genotype networks, first suggested by theoretical models, have been empirically confirmed for proteins and RNAs. Comparative studies also support their existence for gene regulatory networks (GRNs), but direct experimental evidence is lacking. Here, we report the construction of three interconnected genotype networks of synthetic GRNs producing three distinct phenotypes in Escherichia coli. Our synthetic GRNs contain three nodes regulating each other by CRISPR interference and governing the expression of fluorescent reporters. The genotype networks, composed of over twenty different synthetic GRNs, provide robustness in face of mutations while enabling transitions to innovative phenotypes. Through realistic mathematical modeling, we quantify robustness and evolvability for the complete genotype-phenotype map and link these features mechanistically to GRN motifs. Our work thereby exemplifies how GRN evolution along genotype networks might be driving evolutionary innovation.
Collapse
Affiliation(s)
- Javier Santos-Moreno
- Department of Fundamental Microbiology, University of Lausanne, Biophore Building, 1015, Lausanne, Switzerland
- Department of Medicine and Life Sciences, Pompeu Fabra University, 00803, Barcelona, Spain
| | - Eve Tasiudi
- Department of Biosystems Science and Engineering, ETH Zurich and SIB Swiss Institute of Bioinformatics, Basel, Switzerland
| | - Hadiastri Kusumawardhani
- Department of Fundamental Microbiology, University of Lausanne, Biophore Building, 1015, Lausanne, Switzerland
| | - Joerg Stelling
- Department of Biosystems Science and Engineering, ETH Zurich and SIB Swiss Institute of Bioinformatics, Basel, Switzerland.
| | - Yolanda Schaerli
- Department of Fundamental Microbiology, University of Lausanne, Biophore Building, 1015, Lausanne, Switzerland.
| |
Collapse
|
15
|
Chen J, Huang Z, Shi T, Ni D, Zhu Y, Xu W, Zhang W, Mu W. Engineering D-allulose 3-epimerase from Clostridium cellulolyticum for improved thermostability using directed evolution facilitated by a nonenzymatic colorimetric screening assay. FOOD BIOSCI 2023. [DOI: 10.1016/j.fbio.2023.102607] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/03/2023]
|
16
|
Inferring protein fitness landscapes from laboratory evolution experiments. PLoS Comput Biol 2023; 19:e1010956. [PMID: 36857380 PMCID: PMC10010530 DOI: 10.1371/journal.pcbi.1010956] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/05/2022] [Revised: 03/13/2023] [Accepted: 02/16/2023] [Indexed: 03/02/2023] Open
Abstract
Directed laboratory evolution applies iterative rounds of mutation and selection to explore the protein fitness landscape and provides rich information regarding the underlying relationships between protein sequence, structure, and function. Laboratory evolution data consist of protein sequences sampled from evolving populations over multiple generations and this data type does not fit into established supervised and unsupervised machine learning approaches. We develop a statistical learning framework that models the evolutionary process and can infer the protein fitness landscape from multiple snapshots along an evolutionary trajectory. We apply our modeling approach to dihydrofolate reductase (DHFR) laboratory evolution data and the resulting landscape parameters capture important aspects of DHFR structure and function. We use the resulting model to understand the structure of the fitness landscape and find numerous examples of epistasis but an overall global peak that is evolutionarily accessible from most starting sequences. Finally, we use the model to perform an in silico extrapolation of the DHFR laboratory evolution trajectory and computationally design proteins from future evolutionary rounds.
Collapse
|
17
|
Clifton BE, Kozome D, Laurino P. Efficient Exploration of Sequence Space by Sequence-Guided Protein Engineering and Design. Biochemistry 2023; 62:210-220. [PMID: 35245020 DOI: 10.1021/acs.biochem.1c00757] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/02/2023]
Abstract
The rapid growth of sequence databases over the past two decades means that protein engineers faced with optimizing a protein for any given task will often have immediate access to a vast number of related protein sequences. These sequences encode information about the evolutionary history of the protein and the underlying sequence requirements to produce folded, stable, and functional protein variants. Methods that can take advantage of this information are an increasingly important part of the protein engineering tool kit. In this Perspective, we discuss the utility of sequence data in protein engineering and design, focusing on recent advances in three main areas: the use of ancestral sequence reconstruction as an engineering tool to generate thermostable and multifunctional proteins, the use of sequence data to guide engineering of multipoint mutants by structure-based computational protein design, and the use of unlabeled sequence data for unsupervised and semisupervised machine learning, allowing the generation of diverse and functional protein sequences in unexplored regions of sequence space. Altogether, these methods enable the rapid exploration of sequence space within regions enriched with functional proteins and therefore have great potential for accelerating the engineering of stable, functional, and diverse proteins for industrial and biomedical applications.
Collapse
Affiliation(s)
- Ben E Clifton
- Protein Engineering and Evolution Unit, Okinawa Institute of Science and Technology, 1919-1 Tancha, Onna, Okinawa 904-0495, Japan
| | - Dan Kozome
- Protein Engineering and Evolution Unit, Okinawa Institute of Science and Technology, 1919-1 Tancha, Onna, Okinawa 904-0495, Japan
| | - Paola Laurino
- Protein Engineering and Evolution Unit, Okinawa Institute of Science and Technology, 1919-1 Tancha, Onna, Okinawa 904-0495, Japan
| |
Collapse
|
18
|
Kinateder T, Drexler L, Straub K, Merkl R, Sterner R. Experimental and computational analysis of the ancestry of an evolutionary young enzyme from histidine biosynthesis. Protein Sci 2023; 32:e4536. [PMID: 36502290 PMCID: PMC9798254 DOI: 10.1002/pro.4536] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/13/2022] [Revised: 12/01/2022] [Accepted: 12/03/2022] [Indexed: 12/14/2022]
Abstract
The conservation of fold and chemistry of the enzymes associated with histidine biosynthesis suggests that this pathway evolved prior to the diversification of Bacteria, Archaea, and Eukaryotes. The only exception is the histidinol phosphate phosphatase (HolPase). So far, non-homologous HolPases that possess distinct folds and belong to three different protein superfamilies have been identified in various phylogenetic clades. However, their evolution has remained unknown to date. Here, we analyzed the evolutionary history of the HolPase from γ-Proteobacteria (HisB-N). It has been argued that HisB-N and its closest homologue d-glycero-d-manno-heptose-1,7-bisphosphate 7-phosphatase (GmhB) have emerged from the same promiscuous ancestral phosphatase. GmhB variants catalyze the hydrolysis of the anomeric d-glycero-d-manno-heptose-1,7-bisphosphate (αHBP or βHBP) with a strong preference for one anomer (αGmhB or βGmhB). We found that HisB-N from Escherichia coli shows promiscuous activity for βHBP but not αHBP, while βGmhB from Crassaminicella sp. shows promiscuous activity for HolP. Accordingly, a combined phylogenetic tree of αGmhBs, βGmhBs, and HisB-N sequences revealed that HisB-Ns form a compact subcluster derived from βGmhBs. Ancestral sequence reconstruction and in vitro analysis revealed a promiscuous HolPase activity in the resurrected enzymes prior to functional divergence of the successors. The following increase in catalytic efficiency of the HolP turnover is reflected in the shape and electrostatics of the active site predicted by AlphaFold. An analysis of the phylogenetic tree led to a revised evolutionary model that proposes the horizontal gene transfer of a promiscuous βGmhB from δ- to γ-Proteobacteria where it evolved to the modern HisB-N.
Collapse
Affiliation(s)
- Thomas Kinateder
- Institute of Biophysics and Physical Biochemistry and Regensburg Center for Biochemistry, University of RegensburgRegensburgGermany
| | - Lukas Drexler
- Institute of Biophysics and Physical Biochemistry and Regensburg Center for Biochemistry, University of RegensburgRegensburgGermany
| | - Kristina Straub
- Institute of Biophysics and Physical Biochemistry and Regensburg Center for Biochemistry, University of RegensburgRegensburgGermany
| | - Rainer Merkl
- Institute of Biophysics and Physical Biochemistry and Regensburg Center for Biochemistry, University of RegensburgRegensburgGermany
| | - Reinhard Sterner
- Institute of Biophysics and Physical Biochemistry and Regensburg Center for Biochemistry, University of RegensburgRegensburgGermany
| |
Collapse
|
19
|
King E, Maxel S, Zhang Y, Kenney KC, Cui Y, Luu E, Siegel JB, Weiss GA, Luo R, Li H. Orthogonal glycolytic pathway enables directed evolution of noncanonical cofactor oxidase. Nat Commun 2022; 13:7282. [PMID: 36435948 PMCID: PMC9701214 DOI: 10.1038/s41467-022-35021-x] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/11/2022] [Accepted: 11/15/2022] [Indexed: 11/28/2022] Open
Abstract
Noncanonical cofactor biomimetics (NCBs) such as nicotinamide mononucleotide (NMN+) provide enhanced scalability for biomanufacturing. However, engineering enzymes to accept NCBs is difficult. Here, we establish a growth selection platform to evolve enzymes to utilize NMN+-based reducing power. This is based on an orthogonal, NMN+-dependent glycolytic pathway in Escherichia coli which can be coupled to any reciprocal enzyme to recycle the ensuing reduced NMN+. With a throughput of >106 variants per iteration, the growth selection discovers a Lactobacillus pentosus NADH oxidase variant with ~10-fold increase in NMNH catalytic efficiency and enhanced activity for other NCBs. Molecular modeling and experimental validation suggest that instead of directly contacting NCBs, the mutations optimize the enzyme's global conformational dynamics to resemble the WT with the native cofactor bound. Restoring the enzyme's access to catalytically competent conformation states via deep navigation of protein sequence space with high-throughput evolution provides a universal route to engineer NCB-dependent enzymes.
Collapse
Affiliation(s)
- Edward King
- Department of Molecular Biology and Biochemistry, University of California Irvine, Irvine, CA, 92697, USA
| | - Sarah Maxel
- Department Chemical and Biomolecular Engineering University of California Irvine, Irvine, CA, 92697, USA
| | - Yulai Zhang
- Department Chemical and Biomolecular Engineering University of California Irvine, Irvine, CA, 92697, USA
| | - Karissa C Kenney
- Department of Chemistry, University of California Irvine, Irvine, CA, 92697, USA
| | - Youtian Cui
- Genome Center, University of California Davis, Davis, CA, 95616, USA
| | - Emma Luu
- Genome Center, University of California Davis, Davis, CA, 95616, USA
| | - Justin B Siegel
- Genome Center, University of California Davis, Davis, CA, 95616, USA
- Department of Chemistry, Molecular Medicine University of California, Davis, Davis, CA, USA
- Department of Biochemistry and Molecular Medicine University of California, Davis, Davis, CA, USA
| | - Gregory A Weiss
- Department of Molecular Biology and Biochemistry, University of California Irvine, Irvine, CA, 92697, USA
- Department of Chemistry, University of California Irvine, Irvine, CA, 92697, USA
- Department of Pharmaceutical Sciences, University of California Irvine, Irvine, CA, 92697, USA
| | - Ray Luo
- Department of Molecular Biology and Biochemistry, University of California Irvine, Irvine, CA, 92697, USA
- Department Chemical and Biomolecular Engineering University of California Irvine, Irvine, CA, 92697, USA
- Department Materials Science and Engineering, University of California Irvine, Irvine, CA, 92697, USA
- Department of Biomedical Engineering, University of California Irvine, Irvine, CA, 92697, USA
| | - Han Li
- Department Chemical and Biomolecular Engineering University of California Irvine, Irvine, CA, 92697, USA.
- Department of Biomedical Engineering, University of California Irvine, Irvine, CA, 92697, USA.
| |
Collapse
|
20
|
Wittmund M, Cadet F, Davari MD. Learning Epistasis and Residue Coevolution Patterns: Current Trends and Future Perspectives for Advancing Enzyme Engineering. ACS Catal 2022. [DOI: 10.1021/acscatal.2c01426] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022]
Affiliation(s)
- Marcel Wittmund
- Department of Bioorganic Chemistry, Leibniz Institute of Plant Biochemistry, Weinberg 3, 06120 Halle, Germany
| | - Frederic Cadet
- Laboratory of Excellence LABEX GR, DSIMB, Inserm UMR S1134, University of Paris city & University of Reunion, Paris 75014, France
| | - Mehdi D. Davari
- Department of Bioorganic Chemistry, Leibniz Institute of Plant Biochemistry, Weinberg 3, 06120 Halle, Germany
| |
Collapse
|
21
|
Stability and expression of SARS-CoV-2 spike-protein mutations. Mol Cell Biochem 2022; 478:1269-1280. [PMID: 36302994 PMCID: PMC9612610 DOI: 10.1007/s11010-022-04588-w] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/21/2022] [Accepted: 10/12/2022] [Indexed: 12/02/2022]
Abstract
Protein fold stability likely plays a role in SARS-CoV-2 S-protein evolution, together with ACE2 binding and antibody evasion. While few thermodynamic stability data are available for S-protein mutants, many systematic experimental data exist for their expression. In this paper, we explore whether such expression levels relate to the thermodynamic stability of the mutants. We studied mutation-induced SARS-CoV-2 S-protein fold stability, as computed by three very distinct methods and eight different protein structures to account for method- and structure-dependencies. For all methods and structures used (24 comparisons), computed stability changes correlate significantly (99% confidence level) with experimental yeast expression from the literature, such that higher expression is associated with relatively higher fold stability. Also significant, albeit weaker, correlations were seen between stability and ACE2 binding effects. The effect of thermodynamic fold stability may be direct or a correlate of amino acid or site properties, notably the solvent exposure of the site. Correlation between computed stability and experimental expression and ACE2 binding suggests that functional properties of the SARS-CoV-2 S-protein mutant space are largely determined by a few simple features, due to underlying correlations. Our study lends promise to the development of computational tools that may ideally aid in understanding and predicting SARS-CoV-2 S-protein evolution.
Collapse
|
22
|
Thomson RES, Carrera-Pacheco SE, Gillam EMJ. Engineering functional thermostable proteins using ancestral sequence reconstruction. J Biol Chem 2022; 298:102435. [PMID: 36041629 PMCID: PMC9525910 DOI: 10.1016/j.jbc.2022.102435] [Citation(s) in RCA: 15] [Impact Index Per Article: 7.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/31/2022] [Revised: 08/23/2022] [Accepted: 08/24/2022] [Indexed: 11/20/2022] Open
Abstract
Natural proteins are often only slightly more stable in the native state than the denatured state, and an increase in environmental temperature can easily shift the balance toward unfolding. Therefore, the engineering of proteins to improve protein stability is an area of intensive research. Thermostable proteins are required to withstand industrial process conditions, for increased shelf-life of protein therapeutics, for developing robust 'biobricks' for synthetic biology applications, and for research purposes (e.g., structure determination). In addition, thermostability buffers the often destabilizing effects of mutations introduced to improve other properties. Rational design approaches to engineering thermostability require structural information, but even with advanced computational methods, it is challenging to predict or parameterize all the relevant structural factors with sufficient precision to anticipate the results of a given mutation. Directed evolution is an alternative when structures are unavailable but requires extensive screening of mutant libraries. Recently, however, bioinspired approaches based on phylogenetic analyses have shown great promise. Leveraging the rapid expansion in sequence data and bioinformatic tools, ancestral sequence reconstruction can generate highly stable folds for novel applications in industrial chemistry, medicine, and synthetic biology. This review provides an overview of the factors important for successful inference of thermostable proteins by ancestral sequence reconstruction and what it can reveal about the determinants of stability in proteins.
Collapse
Affiliation(s)
- Raine E S Thomson
- School of Chemistry and Molecular Biosciences, The University of Queensland, Brisbane, Australia
| | - Saskya E Carrera-Pacheco
- Centro de Investigación Biomédica (CENBIO), Facultad de Ciencias de la Salud Eugenio Espejo, Universidad UTE, Quito, Ecuador
| | - Elizabeth M J Gillam
- School of Chemistry and Molecular Biosciences, The University of Queensland, Brisbane, Australia.
| |
Collapse
|
23
|
Structural heterogeneity and precision of implications drawn from cryo-electron microscopy structures: SARS-CoV-2 spike-protein mutations as a test case. EUROPEAN BIOPHYSICS JOURNAL 2022; 51:555-568. [PMID: 36167828 PMCID: PMC9514682 DOI: 10.1007/s00249-022-01619-8] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 08/31/2022] [Accepted: 09/19/2022] [Indexed: 11/18/2022]
Abstract
Protein structures may be used to draw functional implications at the residue level, but how sensitive are these implications to the exact structure used? Calculation of the effects of SARS-CoV-2 S-protein mutations based on experimental cryo-electron microscopy structures have been abundant during the pandemic. To understand the precision of such estimates, we studied three distinct methods to estimate stability changes for all possible mutations in 23 different S-protein structures (3.69 million ΔΔG values in total) and explored how random and systematic errors can be remedied by structure-averaged mutation group comparisons. We show that computational estimates have low precision, due to method and structure heterogeneity making results for single mutations uninformative. However, structure-averaged differences in mean effects for groups of substitutions can yield significant results. Illustrating this protocol, functionally important natural mutations, despite individual variations, average to a smaller stability impact compared to other possible mutations, independent of conformational state (open, closed). In summary, we document substantial issues with precision in structure-based protein modeling and recommend sensitivity tests to quantify these effects, but also suggest partial solutions to the problem in the form of structure-averaged “ensemble” estimates for groups of residues when multiple structures are available.
Collapse
|
24
|
A Simulation Analysis and Screening of Deleterious Nonsynonymous Single Nucleotide Polymorphisms (nsSNPs) in Sheep LEP Gene. BIOMED RESEARCH INTERNATIONAL 2022; 2022:7736485. [PMID: 35978633 PMCID: PMC9377880 DOI: 10.1155/2022/7736485] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 05/24/2022] [Accepted: 07/22/2022] [Indexed: 11/17/2022]
Abstract
Leptin is a polypeptide hormone produced in the adipose tissue and governs many processes in the body. Recently, polymorphisms in the LEP gene revealed a significant change in body weight regulation, energy balance, food intake, and reproductive hormone secretion. This study considers its crucial role in the regulation of the economically important traits of sheep. Several computational tools, including SIFT, Predict SNP2, SNAP2, and PROVEAN, have been used to screen out the deleterious nsSNPs. Following the screening of 11 nsSNPs in the sheep genome, 5 nsSNPs, T86M (C → T), D98N (G → A), N136T (A → C), R142Q (G → A), and P157Q (C → A), were predicted to have a significant deleterious effect on the LEP protein function, leading to phenotypic difference. The analysis of proteins’ stability change due to amino acid substitution using the I-stable, SDM, and DynaMut consistently confirmed that three nsSNPs (T86M (C → T), D98N (G → A), and P157Q (C → A)) increased protein stability. It is suggested that these three nsSNPs may enhance the evolvability of LEP protein, which is vital for the evolutionary adaptation of sheep. Our findings demonstrate that the five nsSNPs reported in this study might be responsible for sheep’s structural and functional modifications of LEP protein. This is the first comprehensive report on the sheep LEP gene. It narrow downs the candidate nsSNPs for in vitro experiments to facilitate the development of reliable molecular markers for associated traits.
Collapse
|
25
|
Chattopadhyay G, Bhowmick J, Manjunath K, Ahmed S, Goyal P, Varadarajan R. Mechanistic insights into global suppressors of protein folding defects. PLoS Genet 2022; 18:e1010334. [PMID: 36037221 PMCID: PMC9491731 DOI: 10.1371/journal.pgen.1010334] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/16/2022] [Revised: 09/09/2022] [Accepted: 07/11/2022] [Indexed: 01/14/2023] Open
Abstract
Most amino acid substitutions in a protein either lead to partial loss-of-function or are near neutral. Several studies have shown the existence of second-site mutations that can rescue defects caused by diverse loss-of-function mutations. Such global suppressor mutations are key drivers of protein evolution. However, the mechanisms responsible for such suppression remain poorly understood. To address this, we characterized multiple suppressor mutations both in isolation and in combination with inactive mutants. We examined six global suppressors of the bacterial toxin CcdB, the known M182T global suppressor of TEM-1 β-lactamase, the N239Y global suppressor of p53-DBD and three suppressors of the SARS-CoV-2 spike Receptor Binding Domain. When coupled to inactive mutants, they promote increased in-vivo solubilities as well as regain-of-function phenotypes. In the case of CcdB, where novel suppressors were isolated, we determined the crystal structures of three such suppressors to obtain insight into the specific molecular interactions responsible for the observed effects. While most individual suppressors result in small stability enhancements relative to wildtype, which can be combined to yield significant stability increments, thermodynamic stabilisation is neither necessary nor sufficient for suppressor action. Instead, in diverse systems, we observe that individual global suppressors greatly enhance the foldability of buried site mutants, primarily through increase in refolding rate parameters measured in vitro. In the crowded intracellular environment, mutations that slow down folding likely facilitate off-pathway aggregation. We suggest that suppressor mutations that accelerate refolding can counteract this, enhancing the yield of properly folded, functional protein in vivo.
Collapse
Affiliation(s)
| | - Jayantika Bhowmick
- Molecular Biophysics Unit, Indian Institute of Science, Bangalore,
India
| | - Kavyashree Manjunath
- Centre for Chemical Biology and Therapeutics, Institute For Stem Cell
Science and Regenerative Medicine, Bangalore, India
| | - Shahbaz Ahmed
- Molecular Biophysics Unit, Indian Institute of Science, Bangalore,
India
| | - Parveen Goyal
- Institute for Stem Cell Science and Regenerative Medicine, Bangalore,
India
| | | |
Collapse
|
26
|
Low protein expression enhances phenotypic evolvability by intensifying selection on folding stability. Nat Ecol Evol 2022; 6:1155-1164. [PMID: 35798838 PMCID: PMC7613228 DOI: 10.1038/s41559-022-01797-w] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/17/2021] [Accepted: 05/19/2022] [Indexed: 01/09/2023]
Abstract
Protein abundance affects the evolution of protein genotypes, but we do not know how it affects the evolution of protein phenotypes. Here we investigate the role of protein abundance in the evolvability of green fluorescent protein (GFP) towards the novel phenotype of cyan fluorescence. We evolve GFP in E. coli through multiple cycles of mutation and selection and show that low GFP expression facilitates the evolution of cyan fluorescence. A computational model whose predictions we test experimentally helps explain why: lowly expressed proteins are under stronger selection for proper folding, which facilitates their evolvability on short evolutionary time scales. The reason is that high fluorescence can be achieved by either few proteins that fold well or by many proteins that fold less well. In other words, we observe a synergy between a protein's scarcity and its stability. Because many proteins meet the essential requirements for this scarcity-stability synergy, it may be a widespread mechanism by which low expression helps proteins evolve new phenotypes and functions.
Collapse
|
27
|
Matsumura I, Patrick WM. Dan Tawfik's Lessons for Protein Engineers about Enzymes Adapting to New Substrates. Biochemistry 2022; 62:158-162. [PMID: 35820168 PMCID: PMC9851151 DOI: 10.1021/acs.biochem.2c00230] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/02/2023]
Abstract
Natural evolution has been creating new complex systems for billions of years. The process is spontaneous and requires neither intelligence nor moral purpose but is nevertheless difficult to understand. The late Dan Tawfik spent years studying enzymes as they adapted to recognize new substrates. Much of his work focused on gaining fundamental insights, so the practical utility of his experiments may not be obvious even to accomplished protein engineers. Here we focus on two questions fundamental to any directed evolution experiment. Which proteins are the best starting points for such experiments? Which trait(s) of the chosen parental protein should be evolved to achieve the desired outcome? We summarize Tawfik's contributions to our understanding of these problems, to honor his memory and encourage those unfamiliar with his ideas to read his publications.
Collapse
Affiliation(s)
- Ichiro Matsumura
- O.
Wayne Rollins Research Center, 1510 Clifton Road NE, Room 4001, Atlanta, Georgia 30322, United States,E-mail:
| | - Wayne M. Patrick
- Centre
for Biodiscovery, School of Biological Sciences, Victoria University of Wellington, Wellington 6012, New Zealand,E-mail:
| |
Collapse
|
28
|
Genetic and Phenotypic Study of the Pectobacterium versatile Beta-Lactamase, the Enzyme Most Similar to the Plasmid-Encoded TEM-1. Appl Environ Microbiol 2022; 88:e0022022. [PMID: 35575550 DOI: 10.1128/aem.00220-22] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/20/2022] Open
Abstract
Genus Pectobacterium bacteria include important agricultural pathogens. Pectobacterium versatile isolates contain a chromosome-borne beta-lactamase, PEC-1. This enzyme is the closest relative of TEM-1, a plasmid-borne beta-lactamase widespread in the Enterobacterales. We performed bioinformatics and phenotypic analyses to investigate the genetic and phenotypic features of PEC-1 and its frequency and ability to spread within genus Pectobacterium. We also compared the characteristics of PEC-1 and TEM-1 and evaluated the likelihood of transfer. We found that blaPEC-1 was present principally in a small number of genetic environments in P. versatile. Identical blaPEC-1 genetic environments were present in closely related species, consistent with the high frequency of genetic exchange within the genus Pectobacterium. Despite the similarities between PEC-1 and TEM-1, their genetic environments displayed no significant identity, suggesting an absence of recent transfer. Phenotypic analyses on clonal constructs revealed similar hydrolysis spectra. Our results suggest that P. versatile is the main reservoir of PEC-1, which seems to transfer to closely related species. The genetic distance between PEC-1 and TEM-1, and the lack of conserved elements in their genetic environments, suggest that any transfer that may have occurred must have taken place well before the antibiotic era. IMPORTANCE This study aimed to compare the chromosomal beta-lactamase from Pectobacterium versatile, PEC-1, with the well-known and globally distributed TEM-1 in terms of genetic and functional properties. Despite the similarities between the enzymes, we obtained no definitive proof of gene transfer for the emergence of blaPEC-1 from blaTEM-1. Indeed, given the limited degree of sequence identity and the absence of a common genetic environment, it seems unlikely that any transfer of this gene has occurred recently. However, although blaPEC-1 was found mostly in one specific clade of the P. versatile species, certain isolates from other closely related species, such as Pectobacterium brasiliense and Pectobacterium polaris, may also carry this gene inserted into common genetic environments. This observation suggests that genetic exchanges are frequent, accounting for the diffusion of blaPEC-1 between isolates from different Pectobacterium species and, potentially, to exogenous mobile genetic elements.
Collapse
|
29
|
Local and Global Protein Interactions Contribute to Residue Entrenchment in Beta-Lactamase TEM-1. Antibiotics (Basel) 2022; 11:antibiotics11050652. [PMID: 35625296 PMCID: PMC9137480 DOI: 10.3390/antibiotics11050652] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/12/2022] [Revised: 04/29/2022] [Accepted: 05/05/2022] [Indexed: 11/22/2022] Open
Abstract
Due to their rapid evolution and their impact on healthcare, beta-lactamases, protein degrading beta-lactam antibiotics, are used as generic models of protein evolution. Therefore, we investigated the mutation effects in two distant beta-lactamases, TEM-1 and CTX-M-15. Interestingly, we found a site with a complex pattern of genetic interactions. Mutation G251W in TEM-1 inactivates the protein’s function, just as the reciprocal mutation, W251G, does in CTX-M-15. The phylogenetic analysis revealed that mutation G has been entrenched in TEM-1’s background: while rarely observed throughout the phylogeny, it is essential in TEM-1. Using a rescue experiment, in the TEM-1 G251W mutant, we identified sites that alleviate the deviation from G to W. While few of these mutations could potentially involve local interactions, most of them were found on distant residues in the 3D structure. Many well-known mutations that have an impact on protein stability, such as M182T, were recovered. Our results therefore suggest that entrenchment of an amino acid may rely on diffuse interactions among multiple sites, with a major impact on protein stability.
Collapse
|
30
|
Eliasof M, Boesen T, Haber E, Keasar C, Treister E. Mimetic Neural Networks: A Unified Framework for Protein Design and Folding. FRONTIERS IN BIOINFORMATICS 2022; 2:715006. [DOI: 10.3389/fbinf.2022.715006] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/26/2021] [Accepted: 03/29/2022] [Indexed: 11/13/2022] Open
Abstract
Recent advancements in machine learning techniques for protein structure prediction motivate better results in its inverse problem–protein design. In this work we introduce a new graph mimetic neural network, MimNet, and show that it is possible to build a reversible architecture that solves the structure and design problems in tandem, allowing to improve protein backbone design when the structure is better estimated. We use the ProteinNet data set and show that the state of the art results in protein design can be met and even improved, given recent architectures for protein folding.
Collapse
|
31
|
Abstract
The spike protein (S-protein) of SARS-CoV-2, the protein that enables the virus to infect human cells, is the basis for many vaccines and a hotspot of concerning virus evolution. Here, we discuss the outstanding progress in structural characterization of the S-protein and how these structures facilitate analysis of virus function and evolution. We emphasize the differences in reported structures and that analysis of structure-function relationships is sensitive to the structure used. We show that the average residue solvent exposure in nearly complete structures is a good descriptor of open vs closed conformation states. Because of structural heterogeneity of functionally important surface-exposed residues, we recommend using averages of a group of high-quality protein structures rather than a single structure before reaching conclusions on specific structure-function relationships. To illustrate these points, we analyze some significant chemical tendencies of prominent S-protein mutations in the context of the available structures. In the discussion of new variants, we emphasize the selectivity of binding to ACE2 vs prominent antibodies rather than simply the antibody escape or ACE2 affinity separately. We note that larger chemical changes, in particular increased electrostatic charge or side-chain volume of exposed surface residues, are recurring in mutations of concern, plausibly related to adaptation to the negative surface potential of human ACE2. We also find indications that the fixated mutations of the S-protein in the main variants are less destabilizing than would be expected on average, possibly pointing toward a selection pressure on the S-protein. The richness of available structures for all of these situations provides an enormously valuable basis for future research into these structure-function relationships.
Collapse
Affiliation(s)
- Rukmankesh Mehra
- Department of Chemistry, Indian Institute
of Technology Bhilai, Sejbahar, Raipur 492015, Chhattisgarh,
India
| | - Kasper P. Kepp
- DTU Chemistry, Technical University of
Denmark, Building 206, 2800 Kongens Lyngby,
Denmark
| |
Collapse
|
32
|
Jackson C, Toth-Petroczy A, Kolodny R, Hollfelder F, Fuxreiter M, Caroline Lynn Kamerlin S, Tokuriki N. Adventures on the routes of protein evolution — in memoriam Dan Salah Tawfik (1955 - 2021). J Mol Biol 2022; 434:167462. [DOI: 10.1016/j.jmb.2022.167462] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/06/2022] [Accepted: 01/17/2022] [Indexed: 12/21/2022]
|
33
|
Schneider S, Kozuch J, Boxer SG. The Interplay of Electrostatics and Chemical Positioning in the Evolution of Antibiotic Resistance in TEM β-Lactamases. ACS CENTRAL SCIENCE 2021; 7:1996-2008. [PMID: 34963893 PMCID: PMC8704030 DOI: 10.1021/acscentsci.1c00880] [Citation(s) in RCA: 17] [Impact Index Per Article: 5.7] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/21/2021] [Indexed: 05/25/2023]
Abstract
The interplay of enzyme active site electrostatics and chemical positioning is important for understanding the origin(s) of enzyme catalysis and the design of novel catalysts. We reconstruct the evolutionary trajectory of TEM-1 β-lactamase to TEM-52 toward extended-spectrum activity to better understand the emergence of antibiotic resistance and to provide insights into the structure-function paradigm and noncovalent interactions involved in catalysis. Utilizing a detailed kinetic analysis and the vibrational Stark effect, we quantify the changes in rates and electric fields in the Michaelis and acyl-enzyme complexes for penicillin G and cefotaxime to ascertain the evolutionary role of electric fields to modulate function. These data are combined with MD simulations to interpret and quantify the substrate-dependent structural changes during evolution. We observe that this evolutionary trajectory utilizes a large preorganized electric field and substrate-dependent chemical positioning to facilitate catalysis. This governs the evolvability, substrate promiscuity, and protein fitness landscape in TEM β-lactamase antibiotic resistance.
Collapse
Affiliation(s)
| | | | - Steven G. Boxer
- Chemistry Department, Stanford University, Stanford, California 94305, United States
| |
Collapse
|
34
|
Buchholz PCF, van Loo B, Eenink BDG, Bornberg-Bauer E, Pleiss J. Ancestral sequences of a large promiscuous enzyme family correspond to bridges in sequence space in a network representation. J R Soc Interface 2021; 18:20210389. [PMID: 34727710 DOI: 10.1098/rsif.2021.0389] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022] Open
Abstract
Evolutionary relationships of protein families can be characterized either by networks or by trees. Whereas trees allow for hierarchical grouping and reconstruction of the most likely ancestral sequences, networks lack a time axis but allow for thresholds of pairwise sequence identity to be chosen and, therefore, the clustering of family members with presumably more similar functions. Here, we use the large family of arylsulfatases and phosphonate monoester hydrolases to investigate similarities, strengths and weaknesses in tree and network representations. For varying thresholds of pairwise sequence identity, values of betweenness centrality and clustering coefficients were derived for nodes of the reconstructed ancestors to measure the propensity to act as a bridge in a network. Based on these properties, ancestral protein sequences emerge as bridges in protein sequence networks. Interestingly, many ancestral protein sequences appear close to extant sequences. Therefore, reconstructed ancestor sequences might also be interpreted as yet-to-be-identified homologues. The concept of ancestor reconstruction is compared to consensus sequences, too. It was found that hub sequences in a network, e.g. reconstructed ancestral sequences that are connected to many neighbouring sequences, share closer similarity with derived consensus sequences. Therefore, some reconstructed ancestor sequences can also be interpreted as consensus sequences.
Collapse
Affiliation(s)
- Patrick C F Buchholz
- Institute of Biochemistry and Technical Biochemistry, University of Stuttgart, Allmandring 31, Stuttgart 70569, Germany
| | - Bert van Loo
- Department of Applied Sciences, Northumbria University, Newcastle-upon-Tyne NE1 8ST, UK.,Institute for Evolution and Biodiversity, University of Münster, Hüfferstraße 1, Münster 48149, Germany
| | - Bernard D G Eenink
- Institute for Evolution and Biodiversity, University of Münster, Hüfferstraße 1, Münster 48149, Germany
| | - Erich Bornberg-Bauer
- Institute for Evolution and Biodiversity, University of Münster, Hüfferstraße 1, Münster 48149, Germany.,Department of Protein Evolution, Max Planck Institute for Developmental Biology, Max-Planck-Ring 5, Tübingen 72076, Germany
| | - Jürgen Pleiss
- Institute of Biochemistry and Technical Biochemistry, University of Stuttgart, Allmandring 31, Stuttgart 70569, Germany
| |
Collapse
|
35
|
Jensen ED, Ambri F, Bendtsen MB, Javanpour AA, Liu CC, Jensen MK, Keasling JD. Integrating continuous hypermutation with high-throughput screening for optimization of cis,cis-muconic acid production in yeast. Microb Biotechnol 2021; 14:2617-2626. [PMID: 33645919 PMCID: PMC8601171 DOI: 10.1111/1751-7915.13774] [Citation(s) in RCA: 17] [Impact Index Per Article: 5.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/10/2020] [Revised: 01/31/2021] [Accepted: 02/02/2021] [Indexed: 12/16/2022] Open
Abstract
Directed evolution is a powerful method to optimize proteins and metabolic reactions towards user-defined goals. It usually involves subjecting genes or pathways to iterative rounds of mutagenesis, selection and amplification. While powerful, systematic searches through large sequence-spaces is a labour-intensive task, and can be further limited by a priori knowledge about the optimal initial search space, and/or limits in terms of screening throughput. Here, we demonstrate an integrated directed evolution workflow for metabolic pathway enzymes that continuously generate enzyme variants using the recently developed orthogonal replication system, OrthoRep and screens for optimal performance in high-throughput using a transcription factor-based biosensor. We demonstrate the strengths of this workflow by evolving a rate-limiting enzymatic reaction of the biosynthetic pathway for cis,cis-muconic acid (CCM), a precursor used for bioplastic and coatings, in Saccharomyces cerevisiae. After two weeks of simply iterating between passaging of cells to generate variant enzymes via OrthoRep and high-throughput sorting of best-performing variants using a transcription factor-based biosensor for CCM, we ultimately identified variant enzymes improving CCM titers > 13-fold compared with reference enzymes. Taken together, the combination of synthetic biology tools as adopted in this study is an efficient approach to debottleneck repetitive workflows associated with directed evolution of metabolic enzymes.
Collapse
Affiliation(s)
- Emil D. Jensen
- Novo Nordisk Foundation Center for BiosustainabilityTechnical University of DenmarkKgs. LyngbyDenmark
| | - Francesca Ambri
- Novo Nordisk Foundation Center for BiosustainabilityTechnical University of DenmarkKgs. LyngbyDenmark
| | - Marie B. Bendtsen
- Novo Nordisk Foundation Center for BiosustainabilityTechnical University of DenmarkKgs. LyngbyDenmark
| | - Alex A. Javanpour
- Department of Biomedical EngineeringUniversity of California, IrvineIrvineCA92697USA
| | - Chang C. Liu
- Department of Biomedical EngineeringUniversity of California, IrvineIrvineCA92697USA
- Department of ChemistryUniversity of California, IrvineIrvineCA92697USA
- Department of Molecular Biology and BiochemistryUniversity of California, IrvineIrvineCA92697USA
| | - Michael K. Jensen
- Novo Nordisk Foundation Center for BiosustainabilityTechnical University of DenmarkKgs. LyngbyDenmark
| | - Jay D. Keasling
- Novo Nordisk Foundation Center for BiosustainabilityTechnical University of DenmarkKgs. LyngbyDenmark
- Joint BioEnergy InstituteEmeryvilleCAUSA
- Biological Systems and Engineering DivisionLawrence Berkeley National LaboratoryBerkeleyCAUSA
- Department of Chemical and Biomolecular EngineeringDepartment of BioengineeringUniversity of CaliforniaBerkeleyCAUSA
- Center for Synthetic BiochemistryInstitute for Synthetic BiologyShenzhen Institutes of Advanced TechnologiesShenzhenChina
| |
Collapse
|
36
|
Zheng J, Guo N, Wagner A. Selection enhances protein evolvability by increasing mutational robustness and foldability. Science 2020; 370:370/6521/eabb5962. [DOI: 10.1126/science.abb5962] [Citation(s) in RCA: 39] [Impact Index Per Article: 9.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/09/2020] [Accepted: 09/25/2020] [Indexed: 01/14/2023]
Abstract
Natural selection can promote or hinder a population’s evolvability—the ability to evolve new and adaptive phenotypes—but the underlying mechanisms are poorly understood. To examine how the strength of selection affects evolvability, we subjected populations of yellow fluorescent protein to directed evolution under different selection regimes and then evolved them toward the new phenotype of green fluorescence. Populations under strong selection for the yellow phenotype evolved the green phenotype most rapidly. They did so by accumulating mutations that increase both robustness to mutations and foldability. Under weak selection, neofunctionalizing mutations rose to higher frequency at first, but more frequent deleterious mutations undermined their eventual success. Our experiments show how selection can enhance evolvability by enhancing robustness and create the conditions necessary for evolutionary success.
Collapse
Affiliation(s)
- Jia Zheng
- Department of Evolutionary Biology and Environmental Studies, University of Zurich, Zurich, Switzerland
- Swiss Institute of Bioinformatics, Quartier Sorge-Batiment Genopode, Lausanne, Switzerland
| | - Ning Guo
- Zwirnereistrasse 11, Wallisellen, Zurich, Switzerland
| | - Andreas Wagner
- Department of Evolutionary Biology and Environmental Studies, University of Zurich, Zurich, Switzerland
- Swiss Institute of Bioinformatics, Quartier Sorge-Batiment Genopode, Lausanne, Switzerland
- The Santa Fe Institute, Santa Fe, NM, USA
| |
Collapse
|
37
|
Rix G, Watkins-Dulaney EJ, Almhjell PJ, Boville CE, Arnold FH, Liu CC. Scalable continuous evolution for the generation of diverse enzyme variants encompassing promiscuous activities. Nat Commun 2020; 11:5644. [PMID: 33159067 PMCID: PMC7648111 DOI: 10.1038/s41467-020-19539-6] [Citation(s) in RCA: 44] [Impact Index Per Article: 11.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/24/2020] [Accepted: 10/19/2020] [Indexed: 01/11/2023] Open
Abstract
Enzyme orthologs sharing identical primary functions can have different promiscuous activities. While it is possible to mine this natural diversity to obtain useful biocatalysts, generating comparably rich ortholog diversity is difficult, as it is the product of deep evolutionary processes occurring in a multitude of separate species and populations. Here, we take a first step in recapitulating the depth and scale of natural ortholog evolution on laboratory timescales. Using a continuous directed evolution platform called OrthoRep, we rapidly evolve the Thermotoga maritima tryptophan synthase β-subunit (TmTrpB) through multi-mutation pathways in many independent replicates, selecting only on TmTrpB's primary activity of synthesizing L-tryptophan from indole and L-serine. We find that the resulting sequence-diverse TmTrpB variants span a range of substrate profiles useful in industrial biocatalysis and suggest that the depth and scale of evolution that OrthoRep affords will be generally valuable in enzyme engineering and the evolution of biomolecular functions.
Collapse
Affiliation(s)
- Gordon Rix
- Department of Molecular Biology and Biochemistry, University of California, Irvine, CA, USA
| | - Ella J Watkins-Dulaney
- Division of Biology and Biological Engineering, California Institute of Technology, Pasadena, CA, USA
| | - Patrick J Almhjell
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, CA, USA
| | - Christina E Boville
- Division of Biology and Biological Engineering, California Institute of Technology, Pasadena, CA, USA
- Aralez Bio, Emeryville, CA, USA
| | - Frances H Arnold
- Division of Biology and Biological Engineering, California Institute of Technology, Pasadena, CA, USA
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, CA, USA
| | - Chang C Liu
- Department of Molecular Biology and Biochemistry, University of California, Irvine, CA, USA.
- Department of Biomedical Engineering, University of California, Irvine, CA, USA.
- Department of Chemistry, University of California, Irvine, CA, USA.
| |
Collapse
|
38
|
Kurahashi R, Tanaka SI, Takano K. Highly active enzymes produced by directed evolution with stability-based selection. Enzyme Microb Technol 2020; 140:109626. [DOI: 10.1016/j.enzmictec.2020.109626] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/09/2020] [Revised: 06/12/2020] [Accepted: 06/12/2020] [Indexed: 12/22/2022]
|
39
|
Zhong Z, Wong BG, Ravikumar A, Arzumanyan GA, Khalil AS, Liu CC. Automated Continuous Evolution of Proteins in Vivo. ACS Synth Biol 2020; 9:1270-1276. [PMID: 32374988 PMCID: PMC7370864 DOI: 10.1021/acssynbio.0c00135] [Citation(s) in RCA: 31] [Impact Index Per Article: 7.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/20/2023]
Abstract
We present automated continuous evolution (ACE), a platform for the hands-free directed evolution of biomolecules. ACE pairs OrthoRep, a genetic system for continuous targeted mutagenesis of user-selected genes in vivo, with eVOLVER, a scalable and automated continuous culture device for precise, multiparameter regulation of growth conditions. By implementing real-time feedback-controlled tuning of selection stringency with eVOLVER, genes of interest encoded on OrthoRep autonomously traversed multimutation adaptive pathways to reach desired functions, including drug resistance and improved enzyme activity. The durability, scalability, and speed of biomolecular evolution with ACE should be broadly applicable to protein engineering as well as prospective studies on how selection parameters and schedules shape adaptation.
Collapse
Affiliation(s)
- Ziwei Zhong
- Department of Biomedical Engineering, University of California, Irvine, California, USA
| | - Brandon G. Wong
- Biological Design Center, Boston University, Boston, Massachusetts, USA
- Department of Biomedical Engineering, Boston University, Boston, Massachusetts, USA
| | - Arjun Ravikumar
- Department of Biomedical Engineering, University of California, Irvine, California, USA
- Biological Design Center, Boston University, Boston, Massachusetts, USA
- Department of Biomedical Engineering, Boston University, Boston, Massachusetts, USA
| | - Garri A. Arzumanyan
- Department of Biomedical Engineering, University of California, Irvine, California, USA
| | - Ahmad S. Khalil
- Biological Design Center, Boston University, Boston, Massachusetts, USA
- Department of Biomedical Engineering, Boston University, Boston, Massachusetts, USA
- Wyss Institute for Biologically Inspired Engineering, Harvard University, Boston, Massachusetts, USA
| | - Chang C. Liu
- Department of Biomedical Engineering, University of California, Irvine, California, USA
- Department of Chemistry, University of California, Irvine, California, USA
- Department of Molecular Biology & Biochemistry, University of California, Irvine, California, USA
| |
Collapse
|
40
|
Abstract
The distribution of fitness effects of mutation plays a central role in constraining protein evolution. The underlying mechanisms by which mutations lead to fitness effects are typically attributed to changes in protein specific activity or abundance. Here, we reveal the importance of a mutation's collateral fitness effects, which we define as effects that do not derive from changes in the protein's ability to perform its physiological function. We comprehensively measured the collateral fitness effects of missense mutations in the Escherichia coli TEM-1 β-lactamase antibiotic resistance gene using growth competition experiments in the absence of antibiotic. At least 42% of missense mutations in TEM-1 were deleterious, indicating that for some proteins collateral fitness effects occur as frequently as effects on protein activity and abundance. Deleterious mutations caused improper posttranslational processing, incorrect disulfide-bond formation, protein aggregation, changes in gene expression, and pleiotropic effects on cell phenotype. Deleterious collateral fitness effects occurred more frequently in TEM-1 than deleterious effects on antibiotic resistance in environments with low concentrations of the antibiotic. The surprising prevalence of deleterious collateral fitness effects suggests they may play a role in constraining protein evolution, particularly for highly expressed proteins, for proteins under intermittent selection for their physiological function, and for proteins whose contribution to fitness is buffered against deleterious effects on protein activity and protein abundance.
Collapse
|
41
|
Janesch B, Baumann L, Mark A, Thompson N, Rahmani S, Sim L, Withers SG, Wakarchuk WW. Directed evolution of bacterial polysialyltransferases. Glycobiology 2020; 29:588-598. [PMID: 30976781 DOI: 10.1093/glycob/cwz021] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/06/2018] [Revised: 02/15/2019] [Accepted: 03/18/2019] [Indexed: 01/30/2023] Open
Abstract
Polysialyltransferases (polySTs) are glycosyltransferases that synthesize polymers of sialic acid found in vertebrates and some bacterial pathogens. Bacterial polySTs have utility in the modification of therapeutic proteins to improve serum half-life, and the potential for tissue engineering. PolySTs are membrane-associated proteins and as recombinant proteins suffer from inherently low solubility, low expression levels and poor thermal stability. To improve their physicochemical and biochemical properties, we applied a directed evolution approach using a FACS-based ultrahigh-throughput assay as a simple, robust and reliable screening method. We were able to enrich a large mutant library and, in combination with plate-based high-throughput secondary screening, we discovered mutants with increased enzymatic activity and improved stability compared to the wildtype enzyme. This work presents a powerful strategy for the screening of directed evolution libraries of bacterial polySTs to identify better catalysts for in vitro polysialylation of therapeutics.
Collapse
Affiliation(s)
- Bettina Janesch
- Department of Chemistry and Biology, Ryerson University, Toronto, ON, Canada
| | - Lars Baumann
- Departments of Chemistry and Biochemistry, and Michael Smith Laboratory, University of British Columbia, Vancouver, BC, Canada
| | - Alison Mark
- Department of Chemistry and Biology, Ryerson University, Toronto, ON, Canada
| | - Nicole Thompson
- Department of Chemistry and Biology, Ryerson University, Toronto, ON, Canada
| | - Sadia Rahmani
- Department of Chemistry and Biology, Ryerson University, Toronto, ON, Canada
| | - Lyann Sim
- Departments of Chemistry and Biochemistry, and Michael Smith Laboratory, University of British Columbia, Vancouver, BC, Canada
| | - Stephen G Withers
- Departments of Chemistry and Biochemistry, and Michael Smith Laboratory, University of British Columbia, Vancouver, BC, Canada
| | - Warren W Wakarchuk
- Department of Chemistry and Biology, Ryerson University, Toronto, ON, Canada
| |
Collapse
|
42
|
Stiffler MA, Poelwijk FJ, Brock KP, Stein RR, Riesselman A, Teyra J, Sidhu SS, Marks DS, Gauthier NP, Sander C. Protein Structure from Experimental Evolution. Cell Syst 2020; 10:15-24.e5. [DOI: 10.1016/j.cels.2019.11.008] [Citation(s) in RCA: 24] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/13/2019] [Revised: 09/06/2019] [Accepted: 11/20/2019] [Indexed: 11/27/2022]
|
43
|
van Loo B, Heberlein M, Mair P, Zinchenko A, Schüürmann J, Eenink BDG, Holstein JM, Dilkaute C, Jose J, Hollfelder F, Bornberg-Bauer E. High-Throughput, Lysis-Free Screening for Sulfatase Activity Using Escherichia coli Autodisplay in Microdroplets. ACS Synth Biol 2019; 8:2690-2700. [PMID: 31738524 DOI: 10.1021/acssynbio.9b00274] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Directed evolution of enzymes toward improved catalytic performance has become a powerful tool in protein engineering. To be effective, a directed evolution campaign requires the use of high-throughput screening. In this study we describe the development of an ultra high-throughput lysis-free procedure to screen for improved sulfatase activity by combining microdroplet-based single-variant activity sorting with E. coli autodisplay. For the first step in a 4-step screening procedure, we quantitatively screened >105 variants of the homodimeric arylsulfatase from Silicibacter pomeroyi (SpAS1), displayed on the E. coli cell surface, for improved sulfatase activity using fluorescence activated droplet sorting. Compartmentalization of the fluorescent reaction product with living E. coli cells autodisplaying the sulfatase variants ensured the continuous linkage of genotype and phenotype during droplet sorting and allowed for direct recovery by simple regrowth of the sorted cells. The use of autodisplay on living cells simplified and reduced the degree of liquid handling during all steps in the screening procedure to the single event of simply mixing substrate and cells. The percentage of apparent improved variants was enriched >10-fold as a result of droplet sorting. We ultimately identified 25 SpAS1 variants with improved performance toward 4-nitrophenyl sulfate (up to 6.2-fold) and/or fluorescein disulfate (up to 30-fold). In SpAS1 variants with improved performance toward the bulky fluorescein disulfate, many of the beneficial mutations occur in residues that form hydrogen bonds between α-helices in the C-terminal oligomerization region, suggesting a previously unknown role for the dimer interface in shaping the substrate binding site of SpAS1.
Collapse
Affiliation(s)
- Bert van Loo
- Institute for Evolution and Biodiversity, University of Münster, 48149 Münster, Germany
| | - Magdalena Heberlein
- Institute for Evolution and Biodiversity, University of Münster, 48149 Münster, Germany
- Department of Biochemistry, University of Cambridge, Cambridge CB2 1GA, United Kingdom
| | - Philip Mair
- Department of Biochemistry, University of Cambridge, Cambridge CB2 1GA, United Kingdom
| | - Anastasia Zinchenko
- Department of Biochemistry, University of Cambridge, Cambridge CB2 1GA, United Kingdom
| | - Jan Schüürmann
- Institute of Pharmaceutical and Medicinal Chemistry, University of Münster, 48149 Münster, Germany
| | - Bernard D. G. Eenink
- Institute for Evolution and Biodiversity, University of Münster, 48149 Münster, Germany
| | - Josephin M. Holstein
- Department of Biochemistry, University of Cambridge, Cambridge CB2 1GA, United Kingdom
| | - Carina Dilkaute
- Institute of Pharmaceutical and Medicinal Chemistry, University of Münster, 48149 Münster, Germany
| | - Joachim Jose
- Institute of Pharmaceutical and Medicinal Chemistry, University of Münster, 48149 Münster, Germany
| | - Florian Hollfelder
- Department of Biochemistry, University of Cambridge, Cambridge CB2 1GA, United Kingdom
| | - Erich Bornberg-Bauer
- Institute for Evolution and Biodiversity, University of Münster, 48149 Münster, Germany
| |
Collapse
|
44
|
Leveson-Gower RB, Mayer C, Roelfes G. The importance of catalytic promiscuity for enzyme design and evolution. Nat Rev Chem 2019. [DOI: 10.1038/s41570-019-0143-x] [Citation(s) in RCA: 121] [Impact Index Per Article: 24.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
|
45
|
Activity-stability trade-off in random mutant proteins. J Biosci Bioeng 2019; 128:405-409. [DOI: 10.1016/j.jbiosc.2019.03.017] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/19/2019] [Revised: 03/18/2019] [Accepted: 03/27/2019] [Indexed: 12/23/2022]
|
46
|
Tomala K, Zrebiec P, Hartl DL. Limits to Compensatory Mutations: Insights from Temperature-Sensitive Alleles. Mol Biol Evol 2019; 36:1874-1883. [PMID: 31058959 PMCID: PMC6735812 DOI: 10.1093/molbev/msz110] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/13/2022] Open
Abstract
Previous experiments with temperature-sensitive mutants of the yeast enzyme orotidine 5'-phosphate decarboxylase (encoded in gene URA3) yielded the unexpected result that reversion occurs only through exact reversal of the original mutation (Jakubowska A, Korona R. 2009. Lack of evolutionary conservation at positions important for thermal stability in the yeast ODCase protein. Mol Biol Evol. 26(7):1431-1434.). We recreated a set of these mutations in which the codon had two nucleotide substitutions, making exact reversion much less likely. We screened these double mutants for reversion and obtained a number of compensatory mutations occurring at alternative sites in the molecule. None of these compensatory mutations fully restored protein performance. The mechanism of partial compensation is consistent with a model in which protein stabilization is additive, as the same secondary mutations can compensate different primary alternations. The distance between primary and compensatory residues precludes direct interaction between the sites. Instead, most of the compensatory mutants were clustered in proximity to the catalytic center. All of the second-site compensatory substitutions occurred at relatively conserved sites, and the amino acid replacements were to residues found at these sites in a multispecies alignment of the protein. Based on the estimated distribution of changes in Gibbs free energy among a large number of amino acid replacements, we estimate that, for most proteins, the probability that a second-site mutation would have a sufficiently large stabilizing effect to offset a temperature-sensitive mutation in the order of 10-4 or less. Hence compensation is likely to take place only for slightly destabilizing mutations because highly stabilizing mutations are exceeding rare.
Collapse
Affiliation(s)
- Katarzyna Tomala
- Institute of Environmental Sciences, Faculty of Biology, Jagiellonian University, Krakow, Poland
| | - Piotr Zrebiec
- Institute of Environmental Sciences, Faculty of Biology, Jagiellonian University, Krakow, Poland
| | - Daniel L Hartl
- Department of Organismic and Evolutionary Biology, Harvard University, Cambridge, MA
| |
Collapse
|
47
|
Zadorin AS, Rondelez Y. Selection strategies for randomly partitioned genetic replicators. Phys Rev E 2019; 99:062416. [PMID: 31330587 DOI: 10.1103/physreve.99.062416] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/06/2018] [Indexed: 11/07/2022]
Abstract
The amplification cycle of many replicators (natural or artificial) involves the usage of a host compartment, inside of which the replicator expresses phenotypic compounds necessary to carry out its genetic replication. For example, viruses infect cells, where they express their own proteins and replicate. In this process, the host cell boundary limits the diffusion of the viral protein products, thereby ensuring that phenotypic compounds, such as proteins, promote the replication of the genes that encoded them. This role of maintaining spatial colocalization, also called genotype-phenotype linkage, is a critical function of compartments in natural selection. In most cases, however, individual replicating elements do not distribute systematically among the hosts, but are randomly partitioned. Depending on the replicator-to-host ratio, more than one variant may thus occupy some compartments, blurring the genotype-phenotype linkage and affecting the effectiveness of natural selection. We derive selection equations for a variety of such random multiple occupancy situations, in particular considering the effect of replicator population polymorphism and internal replication dynamics. We conclude that the deleterious effect of random multiple occupancy on selection is relatively benign, and may even completely vanish is some specific cases. In addition, given that higher mean occupancy allows larger populations to be channeled through the selection process, and thus provide a better exploration of phenotypic diversity, we show that it may represent a valid strategy in both natural and technological cases.
Collapse
Affiliation(s)
- Anton S Zadorin
- Laboratoire Gulliver, CNRS, ESPCI Paris, PSL Research University, 75005 Paris, France
| | - Yannick Rondelez
- Laboratoire Gulliver, CNRS, ESPCI Paris, PSL Research University, 75005 Paris, France
| |
Collapse
|
48
|
Validation and Stabilization of a Prophage Lysin of Clostridium perfringens by Using Yeast Surface Display and Coevolutionary Models. Appl Environ Microbiol 2019; 85:AEM.00054-19. [PMID: 30850429 DOI: 10.1128/aem.00054-19] [Citation(s) in RCA: 16] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/08/2019] [Accepted: 02/27/2019] [Indexed: 01/15/2023] Open
Abstract
Bacteriophage lysins are compelling antimicrobial proteins whose biotechnological utility and evolvability would be aided by elevated stability. Lysin catalytic domains, which evolved as modular entities distinct from cell wall binding domains, can be classified into one of several families with highly conserved structure and function, many of which contain thousands of annotated homologous sequences. Motivated by the quality of these evolutionary data, the performance of generative protein models incorporating coevolutionary information was analyzed to predict the stability of variants in a collection of 9,749 multimutants across 10 libraries diversified at different regions of a putative lysin from a prophage region of a Clostridium perfringens genome. Protein stability was assessed via a yeast surface display assay with accompanying high-throughput sequencing. Statistical fitness of mutant sequences, derived from second-order Potts models inferred with different levels of sequence homolog information, was predictive of experimental stability with areas under the curve (AUCs) ranging from 0.78 to 0.85. To extract an experimentally derived model of stability, a logistic model with site-wise score contributions was regressed on the collection of multimutants. This achieved a cross-validated classification performance of 0.95. Using this experimentally derived model, 5 designs incorporating 5 or 6 mutations from multiple libraries were constructed. All designs retained enzymatic activity, with 4 of 5 increasing the melting temperature and with the highest-performing design achieving an improvement of +4°C.IMPORTANCE Bacteriophage lysins exhibit high specificity and activity toward host bacteria with which the phage coevolved. These properties of lysins make them attractive for use as antimicrobials. Although there has been significant effort to develop platforms for rapid lysin engineering, there have been numerous shortcomings when pursuing the ultrahigh throughput necessary for the discovery of rare combinations of mutations to improve performance. In addition to validation of a putative lysin and stabilization thereof, the experimental and computational methods presented here offer a new avenue for improving protein stability and are easily scalable to analysis of tens of millions of mutations in single experiments.
Collapse
|
49
|
Zhong Z, Liu CC. Probing pathways of adaptation with continuous evolution. CURRENT OPINION IN SYSTEMS BIOLOGY 2019; 14:18-24. [PMID: 31608311 PMCID: PMC6788780 DOI: 10.1016/j.coisb.2019.02.002] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/11/2023]
Affiliation(s)
- Ziwei Zhong
- Department of Biomedical Engineering, University of California, Irvine, Irvine, CA 92697, USA
| | - Chang C. Liu
- Department of Biomedical Engineering, University of California, Irvine, Irvine, CA 92697, USA
- Department of Chemistry, University of California, Irvine, Irvine, CA 92697, USA
- Department of Molecular Biology and Biochemistry, University of California, Irvine, Irvine, CA 92697, USA
- Lead Contact
| |
Collapse
|
50
|
Daudé D, Vergès A, Cambon E, Emond S, Tranier S, André I, Remaud-Siméon M. Neutral Genetic Drift-Based Engineering of a Sucrose-Utilizing Enzyme toward Glycodiversification. ACS Catal 2018. [DOI: 10.1021/acscatal.8b03609] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- David Daudé
- Laboratoire d’Ingénierie des Systèmes Biologiques et Procédés (LISBP), Université de Toulouse, CNRS, INRA, INSA, 135, avenue de Rangueil, F-31077 cedex 04 Toulouse, France
| | - Alizée Vergès
- Laboratoire d’Ingénierie des Systèmes Biologiques et Procédés (LISBP), Université de Toulouse, CNRS, INRA, INSA, 135, avenue de Rangueil, F-31077 cedex 04 Toulouse, France
| | - Emmanuelle Cambon
- Laboratoire d’Ingénierie des Systèmes Biologiques et Procédés (LISBP), Université de Toulouse, CNRS, INRA, INSA, 135, avenue de Rangueil, F-31077 cedex 04 Toulouse, France
| | - Stéphane Emond
- Laboratoire d’Ingénierie des Systèmes Biologiques et Procédés (LISBP), Université de Toulouse, CNRS, INRA, INSA, 135, avenue de Rangueil, F-31077 cedex 04 Toulouse, France
| | - Samuel Tranier
- Département Biophysique Structurale, Institut de Pharmacologie et de Biologie Structurale, Université de Toulouse, Université Paul Sabatier, CNRS, F-31077 Toulouse, France
| | - Isabelle André
- Laboratoire d’Ingénierie des Systèmes Biologiques et Procédés (LISBP), Université de Toulouse, CNRS, INRA, INSA, 135, avenue de Rangueil, F-31077 cedex 04 Toulouse, France
| | - Magali Remaud-Siméon
- Laboratoire d’Ingénierie des Systèmes Biologiques et Procédés (LISBP), Université de Toulouse, CNRS, INRA, INSA, 135, avenue de Rangueil, F-31077 cedex 04 Toulouse, France
| |
Collapse
|