1
|
Gholampour-Faroji N, Hemmat J, Haddad-Mashadrizeh A, Asoodeh A. Characterization, structural, and evolutionary analysis of an extremophilic GH5 endoglucanase from Bacillus sp. G131: Insights from ancestral sequence reconstruction. Int J Biol Macromol 2024:134311. [PMID: 39094869 DOI: 10.1016/j.ijbiomac.2024.134311] [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: 02/04/2024] [Revised: 07/16/2024] [Accepted: 07/28/2024] [Indexed: 08/04/2024]
Abstract
Nature has developed extremozymes that catalyze complex reaction processes in extreme environmental conditions. Accordingly, a combined approach consisting of extremozyme screening, ancestral sequence resurrection (ASR), and molecular dynamic simulation was utilized to construct a developed endoglucanase. The primary experimental and in-silico data led to the prediction of a hypothetical sequence of endoglucanase (EG5-G131) using Bacillus sp. G131 confirmed by amplification and sequencing. EG5-G131 exhibited noticeable stability in a broad-pH range, several detergents, organic solvents, and temperatures up to 80 °C. The molecular weight, Vmax, and Km of the purified endoglucanase were estimated to be 36 kDa, 4.32 μmol/min, and 23.62 mg/ml, respectively. The calculated thermodynamic parameters for EG5-G131 confirmed its intrinsic thermostability. Computational analysis revealed Glu142 and Glu230 as active-site residues of the enzyme. Furthermore, the enzyme remained bound to cellotetraose at 298 K, 333 K, 343 K, and 353 K for 300 ns, consistent with our experimental data. ASR of EG5-G131 led to the introduction of ancestral ANC204 and ANC205, which show similar thermodynamic characteristics with the last Firmicute common ancestor. Finally, truncating loops from the N-terminal of two sequences created two variants with desirable thermal stability, suggesting the evolutionary deciphering of the functional domain of the GH5 family in Bacillus sp. G131.
Collapse
Affiliation(s)
- Nazanin Gholampour-Faroji
- Biotechnology Department, Iranian Research Organization for Science and Technology (IROST), Tehran, Iran
| | - Jafar Hemmat
- Biotechnology Department, Iranian Research Organization for Science and Technology (IROST), Tehran, Iran.
| | - Aliakbar Haddad-Mashadrizeh
- Industrial Biotechnology Research Group, Institute of Biotechnology, Ferdowsi University of Mashhad, Mashhad, Iran.
| | - Ahmad Asoodeh
- Industrial Biotechnology Research Group, Institute of Biotechnology, Ferdowsi University of Mashhad, Mashhad, Iran; Department of Chemistry, Faculty of Science, Ferdowsi University of Mashhad, Mashhad, Iran
| |
Collapse
|
2
|
Chisholm LO, Orlandi KN, Phillips SR, Shavlik MJ, Harms MJ. Ancestral Reconstruction and the Evolution of Protein Energy Landscapes. Annu Rev Biophys 2024; 53:127-146. [PMID: 38134334 PMCID: PMC11192866 DOI: 10.1146/annurev-biophys-030722-125440] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2023]
Abstract
A protein's sequence determines its conformational energy landscape. This, in turn, determines the protein's function. Understanding the evolution of new protein functions therefore requires understanding how mutations alter the protein energy landscape. Ancestral sequence reconstruction (ASR) has proven a valuable tool for tackling this problem. In ASR, one phylogenetically infers the sequences of ancient proteins, allowing characterization of their properties. When coupled to biophysical, biochemical, and functional characterization, ASR can reveal how historical mutations altered the energy landscape of ancient proteins, allowing the evolution of enzyme activity, altered conformations, binding specificity, oligomerization, and many other protein features. In this article, we review how ASR studies have been used to dissect the evolution of energy landscapes. We also discuss ASR studies that reveal how energy landscapes have shaped protein evolution. Finally, we propose that thinking about evolution from the perspective of an energy landscape can improve how we approach and interpret ASR studies.
Collapse
Affiliation(s)
- Lauren O Chisholm
- Department of Chemistry and Biochemistry, University of Oregon, Eugene, Oregon, USA;
- Institute of Molecular Biology, University of Oregon, Eugene, Oregon, USA
| | - Kona N Orlandi
- Institute of Molecular Biology, University of Oregon, Eugene, Oregon, USA
- Department of Biology, University of Oregon, Eugene, Oregon, USA
| | - Sophia R Phillips
- Department of Chemistry and Biochemistry, University of Oregon, Eugene, Oregon, USA;
- Institute of Molecular Biology, University of Oregon, Eugene, Oregon, USA
| | - Michael J Shavlik
- Institute of Molecular Biology, University of Oregon, Eugene, Oregon, USA
- Department of Biology, University of Oregon, Eugene, Oregon, USA
| | - Michael J Harms
- Department of Chemistry and Biochemistry, University of Oregon, Eugene, Oregon, USA;
- Institute of Molecular Biology, University of Oregon, Eugene, Oregon, USA
| |
Collapse
|
3
|
Fang N, Wu L, Duan S, Li J. The Structural and Molecular Mechanisms of Mycobacterium tuberculosis Translational Elongation Factor Proteins. Molecules 2024; 29:2058. [PMID: 38731549 PMCID: PMC11085428 DOI: 10.3390/molecules29092058] [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/04/2024] [Revised: 04/19/2024] [Accepted: 04/24/2024] [Indexed: 05/13/2024] Open
Abstract
Targeting translation factor proteins holds promise for developing innovative anti-tuberculosis drugs. During protein translation, many factors cause ribosomes to stall at messenger RNA (mRNA). To maintain protein homeostasis, bacteria have evolved various ribosome rescue mechanisms, including the predominant trans-translation process, to release stalled ribosomes and remove aberrant mRNAs. The rescue systems require the participation of translation elongation factor proteins (EFs) and are essential for bacterial physiology and reproduction. However, they disappear during eukaryotic evolution, which makes the essential proteins and translation elongation factors promising antimicrobial drug targets. Here, we review the structural and molecular mechanisms of the translation elongation factors EF-Tu, EF-Ts, and EF-G, which play essential roles in the normal translation and ribosome rescue mechanisms of Mycobacterium tuberculosis (Mtb). We also briefly describe the structure-based, computer-assisted study of anti-tuberculosis drugs.
Collapse
Affiliation(s)
- Ning Fang
- State Key Laboratory of Genetic Engineering, School of Life Sciences, Shanghai Engineering Research Center of Industrial Microorganisms, Fudan University, Shanghai 200438, China; (N.F.); (L.W.)
| | - Lingyun Wu
- State Key Laboratory of Genetic Engineering, School of Life Sciences, Shanghai Engineering Research Center of Industrial Microorganisms, Fudan University, Shanghai 200438, China; (N.F.); (L.W.)
| | - Shuyan Duan
- State Key Laboratory of Genetic Engineering, School of Life Sciences, Shanghai Engineering Research Center of Industrial Microorganisms, Fudan University, Shanghai 200438, China; (N.F.); (L.W.)
- College of Food Science and Pharmaceutical Engineering, Zaozhuang University, Zaozhuang 277160, China
| | - Jixi Li
- State Key Laboratory of Genetic Engineering, School of Life Sciences, Shanghai Engineering Research Center of Industrial Microorganisms, Fudan University, Shanghai 200438, China; (N.F.); (L.W.)
| |
Collapse
|
4
|
Sucharski F, Gallo G, Coelho C, Hardy L, Würtele M. Modeling the role of charged residues in thermophilic proteins by rotamer and dynamic cross correlation analysis. J Mol Model 2023; 29:132. [PMID: 37036538 DOI: 10.1007/s00894-023-05490-y] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/05/2022] [Accepted: 02/24/2023] [Indexed: 04/11/2023]
Abstract
Discerning the determinants of protein thermostability is very important both from the theoretical and applied perspective. Different lines of evidence seem to indicate that a dynamical network of salt bridges/charged residues plays a fundamental role in the thermostability of enzymes. In this work, we applied measures of dynamic variance, like the Gini coefficients, Kullback-Leibler (KL) divergence and dynamic cross correlation (DCC) coefficients to compare the behavior of 3 pairs of homologous proteins from the thermophilic bacterium Thermus thermophilus and mesophilic Escherichia coli. Molecular dynamic (MD) simulations of these proteins were performed at 303 K and 363 K. In the characterization of their side chain rotamer distributions, the corresponding Gini coefficients and KL-divergence both revealed significant correlations with temperature. Similarly, a DCC analysis revealed a higher trend to de-correlate the movement of charged residues at higher temperatures in the thermophilic proteins, when compared with their mesophilic homologues. These results highlight the importance of dynamic electrostatic network interactions for the thermostability of enzymes.
Collapse
Affiliation(s)
- Fernanda Sucharski
- Department of Science and Technology, Federal University of São Paulo, Talim 330, São José Dos Campos, São Paulo, 12231-280, Brazil
| | - Gloria Gallo
- Department of Science and Technology, Federal University of São Paulo, Talim 330, São José Dos Campos, São Paulo, 12231-280, Brazil
| | - Camila Coelho
- Department of Science and Technology, Federal University of São Paulo, Talim 330, São José Dos Campos, São Paulo, 12231-280, Brazil
| | - Leon Hardy
- Department of Physics, University of South Florida, Tampa, USA
| | - Martin Würtele
- Department of Science and Technology, Federal University of São Paulo, Talim 330, São José Dos Campos, São Paulo, 12231-280, Brazil.
| |
Collapse
|
5
|
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
|
6
|
Rational Design of a Thermostable 2'-Deoxyribosyltransferase for Nelarabine Production by Prediction of Disulfide Bond Engineering Sites. Int J Mol Sci 2022; 23:ijms231911806. [PMID: 36233108 PMCID: PMC9570332 DOI: 10.3390/ijms231911806] [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: 08/03/2022] [Revised: 09/27/2022] [Accepted: 09/29/2022] [Indexed: 11/24/2022] Open
Abstract
One of the major drawbacks of the industrial implementation of enzymatic processes is the low operational stability of the enzymes under tough industrial conditions. In this respect, the use of thermostable enzymes in the industry is gaining ground during the last decades. Herein, we report a structure-guided approach for the development of novel and thermostable 2′-deoxyribosyltransferases (NDTs) based on the computational design of disulfide bonds on hot spot positions. To this end, a small library of NDT variants from Lactobacillus delbrueckii (LdNDT) with introduced cysteine pairs was created. Among them, LdNDTS104C (100% retained activity) was chosen as the most thermostable variant, displaying a six- and two-fold enhanced long-term stability when stored at 55 °C (t1/255 °C ≈ 24 h) and 60 °C (t1/260 °C ≈ 4 h), respectively. Moreover, the biochemical characterization revealed that LdNDTS104C showed >60% relative activity across a broad range of temperature (30−90 °C) and pH (5−7). Finally, to study the potential application of LdNDTS104C as an industrial catalyst, the enzymatic synthesis of nelarabine was successfully carried out under different substrate conditions (1:1 and 3:1) at different reaction times. Under these experimental conditions, the production of nelarabine was increased up to 2.8-fold (72% conversion) compared with wild-type LdNDT.
Collapse
|
7
|
Fer E, McGrath KM, Guy L, Hockenberry AJ, Kaçar B. Early divergence of translation initiation and elongation factors. Protein Sci 2022; 31:e4393. [PMID: 36250475 PMCID: PMC9601768 DOI: 10.1002/pro.4393] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/11/2022] [Revised: 07/05/2022] [Accepted: 07/11/2022] [Indexed: 11/18/2022]
Abstract
Protein translation is a foundational attribute of all living cells. The translation function carried out by the ribosome critically depends on an assortment of protein interaction partners, collectively referred to as the translation machinery. Various studies suggest that the diversification of the translation machinery occurred prior to the last universal common ancestor, yet it is unclear whether the predecessors of the extant translation machinery factors were functionally distinct from their modern counterparts. Here we reconstructed the shared ancestral trajectory and subsequent evolution of essential translation factor GTPases, elongation factor EF-Tu (aEF-1A/eEF-1A), and initiation factor IF2 (aIF5B/eIF5B). Based upon their similar functions and structural homologies, it has been proposed that EF-Tu and IF2 emerged from an ancient common ancestor. We generated the phylogenetic tree of IF2 and EF-Tu proteins and reconstructed ancestral sequences corresponding to the deepest nodes in their shared evolutionary history, including the last common IF2 and EF-Tu ancestor. By identifying the residue and domain substitutions, as well as structural changes along the phylogenetic history, we developed an evolutionary scenario for the origins, divergence and functional refinement of EF-Tu and IF2 proteins. Our analyses suggest that the common ancestor of IF2 and EF-Tu was an IF2-like GTPase protein. Given the central importance of the translation machinery to all cellular life, its earliest evolutionary constraints and trajectories are key to characterizing the universal constraints and capabilities of cellular evolution.
Collapse
Affiliation(s)
- Evrim Fer
- Department of BacteriologyUniversity of Wisconsin‐MadisonMadisonWisconsinUSA
- Microbiology Doctoral Training ProgramUniversity of Wisconsin‐MadisonMadisonWisconsinUSA
- NASA Center for Early Life and EvolutionUniversity of Wisconsin‐MadisonMadisonWisconsinUSA
| | - Kaitlyn M. McGrath
- Department of BacteriologyUniversity of Wisconsin‐MadisonMadisonWisconsinUSA
- NASA Center for Early Life and EvolutionUniversity of Wisconsin‐MadisonMadisonWisconsinUSA
- Department of Molecular and Cellular BiologyUniversity of ArizonaTucsonArizonaUSA
| | - Lionel Guy
- Department of Medical Biochemistry and Microbiology, Science for Life LaboratoryUppsala UniversityUppsalaSweden
| | - Adam J. Hockenberry
- Department of Integrative BiologyThe University of Texas at AustinAustinTexasUSA
| | - Betül Kaçar
- Department of BacteriologyUniversity of Wisconsin‐MadisonMadisonWisconsinUSA
- NASA Center for Early Life and EvolutionUniversity of Wisconsin‐MadisonMadisonWisconsinUSA
| |
Collapse
|
8
|
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] [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 towards 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 (ASR) 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 ASR and what it can reveal about the determinants of stability in proteins.
Collapse
|
9
|
Rahban M, Zolghadri S, Salehi N, Ahmad F, Haertlé T, Rezaei-Ghaleh N, Sawyer L, Saboury AA. Thermal stability enhancement: Fundamental concepts of protein engineering strategies to manipulate the flexible structure. Int J Biol Macromol 2022; 214:642-654. [DOI: 10.1016/j.ijbiomac.2022.06.154] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/17/2022] [Revised: 06/22/2022] [Accepted: 06/23/2022] [Indexed: 01/28/2023]
|
10
|
Gonzalez NA, Li BA, McCully ME. The stability and dynamics of computationally designed proteins. Protein Eng Des Sel 2022; 35:6529794. [PMID: 35174855 DOI: 10.1093/protein/gzac001] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/14/2020] [Revised: 01/12/2022] [Accepted: 01/13/2022] [Indexed: 12/11/2022] Open
Abstract
Protein stability, dynamics and function are intricately linked. Accordingly, protein designers leverage dynamics in their designs and gain insight to their successes and failures by analyzing their proteins' dynamics. Molecular dynamics (MD) simulations are a powerful computational tool for quantifying both local and global protein dynamics. This review highlights studies where MD simulations were applied to characterize the stability and dynamics of designed proteins and where dynamics were incorporated into computational protein design. First, we discuss the structural basis underlying the extreme stability and thermostability frequently observed in computationally designed proteins. Next, we discuss examples of designed proteins, where dynamics were not explicitly accounted for in the design process, whose coordinated motions or active site dynamics, as observed by MD simulation, enhanced or detracted from their function. Many protein functions depend on sizeable or subtle conformational changes, so we finally discuss the computational design of proteins to perform a specific function that requires consideration of motion by multi-state design.
Collapse
Affiliation(s)
- Natali A Gonzalez
- Department of Biology, Santa Clara University, 500 El Camino Real, Santa Clara, CA 95053, USA
| | - Brigitte A Li
- Department of Biology, Santa Clara University, 500 El Camino Real, Santa Clara, CA 95053, USA
| | - Michelle E McCully
- Department of Biology, Santa Clara University, 500 El Camino Real, Santa Clara, CA 95053, USA
| |
Collapse
|
11
|
Akanuma S, Yamaguchi M, Yamagishi A. Comprehensive mutagenesis to identify amino acid residues contributing to the difference in thermostability between two originally thermostable ancestral proteins. PLoS One 2021; 16:e0258821. [PMID: 34673819 PMCID: PMC8530338 DOI: 10.1371/journal.pone.0258821] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/05/2021] [Accepted: 10/05/2021] [Indexed: 11/19/2022] Open
Abstract
Further improvement of the thermostability of inherently thermostable proteins is an attractive challenge because more thermostable proteins are industrially more useful and serve as better scaffolds for protein engineering. To establish guidelines that can be applied for the rational design of hyperthermostable proteins, we compared the amino acid sequences of two ancestral nucleoside diphosphate kinases, Arc1 and Bac1, reconstructed in our previous study. Although Bac1 is a thermostable protein whose unfolding temperature is around 100°C, Arc1 is much more thermostable with an unfolding temperature of 114°C. However, only 12 out of 139 amino acids are different between the two sequences. In this study, one or a combination of amino acid(s) in Bac1 was/were substituted by a residue(s) found in Arc1 at the same position(s). The best mutant, which contained three amino acid substitutions (S108D, G116A and L120P substitutions), showed an unfolding temperature more than 10°C higher than that of Bac1. Furthermore, a combination of the other nine amino acid substitutions also led to improved thermostability of Bac1, although the effects of individual substitutions were small. Therefore, not only the sum of the contributions of individual amino acids, but also the synergistic effects of multiple amino acids are deeply involved in the stability of a hyperthermostable protein. Such insights will be helpful for future rational design of hyperthermostable proteins.
Collapse
Affiliation(s)
- Satoshi Akanuma
- Faculty of Human Sciences, Waseda University, Tokorozawa, Saitama, Japan
- * E-mail:
| | - Minako Yamaguchi
- Department of Applied Life Sciences, Tokyo University of Pharmacy and Life Sciences, Hachioji, Tokyo, Japan
| | - Akihiko Yamagishi
- Department of Applied Life Sciences, Tokyo University of Pharmacy and Life Sciences, Hachioji, Tokyo, Japan
| |
Collapse
|
12
|
Spence MA, Kaczmarski JA, Saunders JW, Jackson CJ. Ancestral sequence reconstruction for protein engineers. Curr Opin Struct Biol 2021; 69:131-141. [PMID: 34023793 DOI: 10.1016/j.sbi.2021.04.001] [Citation(s) in RCA: 53] [Impact Index Per Article: 17.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/15/2020] [Revised: 03/22/2021] [Accepted: 04/07/2021] [Indexed: 12/11/2022]
Abstract
In addition to its value in the study of molecular evolution, ancestral sequence reconstruction (ASR) has emerged as a useful methodology for engineering proteins with enhanced properties. Proteins generated by ASR often exhibit unique or improved activity, stability, and/or promiscuity, all of which are properties that are valued by protein engineers. Comparison between extant proteins and evolutionary intermediates generated by ASR also allows protein engineers to identify substitutions that have contributed to functional innovation or diversification within protein families. As ASR becomes more widely adopted as a protein engineering approach, it is important to understand the applications, limitations, and recent developments of this technique. This review highlights recent exemplifications of ASR, as well as technical aspects of the reconstruction process that are relevant to protein engineering.
Collapse
Affiliation(s)
- Matthew A Spence
- Research School of Chemistry, Australian National University, Canberra, ACT 2601, Australia
| | - Joe A Kaczmarski
- Research School of Chemistry, Australian National University, Canberra, ACT 2601, Australia
| | - Jake W Saunders
- Research School of Chemistry, Australian National University, Canberra, ACT 2601, Australia
| | - Colin J Jackson
- Research School of Chemistry, Australian National University, Canberra, ACT 2601, Australia; ARC Centre of Excellence for Innovations in Peptide & Protein Science, Research School of Chemistry, Australian National University, Canberra, ACT 2601, Australia; ARC Centre of Excellence for Innovations in Synthetic Biology, Research School of Chemistry, Australian National University, Canberra, ACT 2601, Australia.
| |
Collapse
|
13
|
Modi T, Campitelli P, Kazan IC, Ozkan SB. Protein folding stability and binding interactions through the lens of evolution: a dynamical perspective. Curr Opin Struct Biol 2020; 66:207-215. [PMID: 33388636 DOI: 10.1016/j.sbi.2020.11.007] [Citation(s) in RCA: 16] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/17/2020] [Revised: 11/02/2020] [Accepted: 11/26/2020] [Indexed: 01/06/2023]
Abstract
While the function of a protein depends heavily on its ability to fold into a correct 3D structure, billions of years of evolution have tailored proteins from highly stable objects to flexible molecules as they adapted to environmental changes. Nature maintains the fine balance of protein folding and stability while still evolving towards new function through generations of fine-tuning necessary interactions with other proteins and small molecules. Here we focus on recent computational and experimental studies that shed light onto how evolution molds protein folding and the functional landscape from a conformational dynamics' perspective. Particularly, we explore the importance of dynamic allostery throughout protein evolution and discuss how the protein anisotropic network can give rise to allosteric and epistatic interactions.
Collapse
Affiliation(s)
- Tushar Modi
- Department of Physics and Center for Biological Physics, Arizona State University, Tempe, AZ 85287-1504, USA
| | - Paul Campitelli
- Department of Physics and Center for Biological Physics, Arizona State University, Tempe, AZ 85287-1504, USA
| | - Ismail Can Kazan
- Department of Physics and Center for Biological Physics, Arizona State University, Tempe, AZ 85287-1504, USA
| | - Sefika Banu Ozkan
- Department of Physics and Center for Biological Physics, Arizona State University, Tempe, AZ 85287-1504, USA.
| |
Collapse
|
14
|
Sternke M, Tripp KW, Barrick D. The use of consensus sequence information to engineer stability and activity in proteins. Methods Enzymol 2020; 643:149-179. [PMID: 32896279 DOI: 10.1016/bs.mie.2020.06.001] [Citation(s) in RCA: 19] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Abstract
The goal of protein design is to create proteins that are stable, soluble, and active. Here we focus on one approach to protein design in which sequence information is used to create a "consensus" sequence. Such consensus sequences comprise the most common residue at each position in a multiple sequence alignment (MSA). After describing some general ideas that relate MSA and consensus sequences and presenting a statistical thermodynamic framework that relates consensus and non-consensus sequences to stability, we detail the process of designing a consensus sequence and survey reports of consensus design and characterization from the literature. Many of these consensus proteins retain native biological activities including ligand binding and enzyme activity. Remarkably, in most cases the consensus protein shows significantly higher stability than extant versions of the protein, as measured by thermal or chemical denaturation, consistent with the statistical thermodynamic model. To understand this stability increase, we compare various features of consensus sequences with the extant MSA sequences from which they were derived. Consensus sequences show enrichment in charged residues (most notably glutamate and lysine) and depletion of uncharged polar residues (glutamine, serine, and asparagine). Surprisingly, a survey of stability changes resulting from point substitutions show little correlation with residue frequencies at the corresponding positions within the MSA, suggesting that the high stability of consensus proteins may result from interactions among residue pairs or higher-order clusters. Whatever the source, the large number of reported successes demonstrates that consensus design is a viable route to generating active and in many cases highly stabilized proteins.
Collapse
Affiliation(s)
- Matt Sternke
- T.C. Jenkins Department of Biophysics, Johns Hopkins University, Baltimore, MD, United States; Program in Molecular Biophysics, Johns Hopkins University, Baltimore, MD, United States
| | - Katherine W Tripp
- T.C. Jenkins Department of Biophysics, Johns Hopkins University, Baltimore, MD, United States
| | - Doug Barrick
- T.C. Jenkins Department of Biophysics, Johns Hopkins University, Baltimore, MD, United States.
| |
Collapse
|
15
|
Babkova P, Dunajova Z, Chaloupkova R, Damborsky J, Bednar D, Marek M. Structures of hyperstable ancestral haloalkane dehalogenases show restricted conformational dynamics. Comput Struct Biotechnol J 2020; 18:1497-1508. [PMID: 32637047 PMCID: PMC7327271 DOI: 10.1016/j.csbj.2020.06.021] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/07/2020] [Revised: 06/08/2020] [Accepted: 06/10/2020] [Indexed: 12/30/2022] Open
Abstract
Ancestral sequence reconstruction is a powerful method for inferring ancestors of modern enzymes and for studying structure-function relationships of enzymes. We have previously applied this approach to haloalkane dehalogenases (HLDs) from the subfamily HLD-II and obtained thermodynamically highly stabilized enzymes (ΔT m up to 24 °C), showing improved catalytic properties. Here we combined crystallographic structural analysis and computational molecular dynamics simulations to gain insight into the mechanisms by which ancestral HLDs became more robust enzymes with novel catalytic properties. Reconstructed ancestors exhibited similar structure topology as their descendants with the exception of a few loop deviations. Strikingly, molecular dynamics simulations revealed restricted conformational dynamics of ancestral enzymes, which prefer a single state, in contrast to modern enzymes adopting two different conformational states. The restricted dynamics can potentially be linked to their exceptional stabilization. The study provides molecular insights into protein stabilization due to ancestral sequence reconstruction, which is becoming a widely used approach for obtaining robust protein catalysts.
Collapse
Affiliation(s)
- Petra Babkova
- Loschmidt Laboratories, Department of Experimental Biology and RECETOX, Faculty of Science, Masaryk University, Kamenice 5, Bld. A13, 625 00 Brno, Czech Republic
- International Clinical Research Center, St. Anne's University Hospital Brno, Pekarska 53, 656 91 Brno, Czech Republic
| | - Zuzana Dunajova
- Loschmidt Laboratories, Department of Experimental Biology and RECETOX, Faculty of Science, Masaryk University, Kamenice 5, Bld. A13, 625 00 Brno, Czech Republic
| | - Radka Chaloupkova
- Loschmidt Laboratories, Department of Experimental Biology and RECETOX, Faculty of Science, Masaryk University, Kamenice 5, Bld. A13, 625 00 Brno, Czech Republic
| | - Jiri Damborsky
- Loschmidt Laboratories, Department of Experimental Biology and RECETOX, Faculty of Science, Masaryk University, Kamenice 5, Bld. A13, 625 00 Brno, Czech Republic
- International Clinical Research Center, St. Anne's University Hospital Brno, Pekarska 53, 656 91 Brno, Czech Republic
| | - David Bednar
- Loschmidt Laboratories, Department of Experimental Biology and RECETOX, Faculty of Science, Masaryk University, Kamenice 5, Bld. A13, 625 00 Brno, Czech Republic
- International Clinical Research Center, St. Anne's University Hospital Brno, Pekarska 53, 656 91 Brno, Czech Republic
| | - Martin Marek
- Loschmidt Laboratories, Department of Experimental Biology and RECETOX, Faculty of Science, Masaryk University, Kamenice 5, Bld. A13, 625 00 Brno, Czech Republic
| |
Collapse
|
16
|
Noh D, Bower RL, Hay DL, Zhyvoloup A, Raleigh DP. Analysis of Amylin Consensus Sequences Suggests That Human Amylin Is Not Optimized to Minimize Amyloid Formation and Provides Clues to Factors That Modulate Amyloidogenicity. ACS Chem Biol 2020; 15:1408-1416. [PMID: 32364695 DOI: 10.1021/acschembio.9b01050] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Abstract
The neuropancreatic polypeptide hormone amylin forms pancreatic islet amyloid in type-2 diabetes. Islet amyloid formation contributes to β-cell death in the disease and to the failure of islet transplants, but the features which influence amylin amyloidogenicity are not understood. We constructed an amino acid sequence alignment of 202 sequences of amylin and used the alignment to design consensus sequences of vertebrate amylins, mammalian amylins, and primate amylins. Amylin is highly conserved, but there are differences between human amylin and each consensus sequence, ranging from one to six substitutions. Biophysical analysis shows that all of the consensus sequences form amyloid but do so more slowly than human amylin in vitro. The rate of amyloid formation by the primate consensus sequence is 3- to 4-fold slower than human amylin; the mammalian consensus sequence is approximately 20- to 25-fold slower, and the vertebrate consensus sequence is approximately 6-fold slower. All of the consensus sequences are moderately less toxic than human amylin toward a cultured β-cell line, with the vertebrate consensus sequence displaying the largest reduction in toxicity of 3- to 4-fold. All of the consensus sequences activate a human amylin receptor and exhibit only modest reductions in activity, ranging from 3- to 4-fold as judged by a cAMP production assay. The analysis argues that there is no strong selective evolutionary pressure to avoid the formation of islet amyloid and provides information relevant to the design of less amyloidogenic amylin variants.
Collapse
Affiliation(s)
- Daeun Noh
- Graduate Program in Biochemistry and Structural Biology, Stony Brook University, Stony Brook, New York 11790, United States
| | - Rebekah L. Bower
- School of Biological Sciences and Maurice Wilkins Centre for Molecular Biodiscovery, University of Auckland, Auckland 1142, New Zealand
| | - Debbie L. Hay
- School of Biological Sciences and Maurice Wilkins Centre for Molecular Biodiscovery, University of Auckland, Auckland 1142, New Zealand
| | - Alexander Zhyvoloup
- Institute of Structural and Molecular Biology, University College London, Gower Street, London, WC1E 6BT United Kingdom
| | - Daniel P. Raleigh
- Institute of Structural and Molecular Biology, University College London, Gower Street, London, WC1E 6BT United Kingdom
- Department of Chemistry, Stony Brook University, Stony Brook, New York 11790, United States
- Laufer Center for Quantitative Biology, Stony Brook University, Stony Brook, New York 11790, United States
| |
Collapse
|
17
|
Girodat D, Blanchard SC, Wieden HJ, Sanbonmatsu KY. Elongation Factor Tu Switch I Element is a Gate for Aminoacyl-tRNA Selection. J Mol Biol 2020; 432:3064-3077. [PMID: 32061931 DOI: 10.1016/j.jmb.2020.01.038] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/15/2019] [Revised: 01/15/2020] [Accepted: 01/24/2020] [Indexed: 12/16/2022]
Abstract
Selection of correct aminoacyl (aa)-tRNA at the ribosomal A site is fundamental to maintaining translational fidelity. Aa-tRNA selection is a multistep process facilitated by the guanosine triphosphatase elongation factor (EF)-Tu. EF-Tu delivers aa-tRNA to the ribosomal A site and participates in tRNA selection. The structural mechanism of how EF-Tu is involved in proofreading remains to be fully resolved. Here, we provide evidence that switch I of EF-Tu facilitates EF-Tu's involvement during aa-tRNA selection. Using structure-based and explicit solvent molecular dynamics simulations based on recent cryo-electron microscopy reconstructions, we studied the conformational change of EF-Tu from the guanosine triphosphate to guanine diphosphate conformation during aa-tRNA accommodation. Switch I of EF-Tu rapidly converts from an α-helix into a β-hairpin and moves to interact with the acceptor stem of the aa-tRNA. In doing so, switch I gates the movement of the aa-tRNA during accommodation through steric interactions with the acceptor stem. Pharmacological inhibition of the aa-tRNA accommodation pathway prevents the proper positioning of switch I with the aa-tRNA acceptor stem, suggesting that the observed interactions are specific for cognate aa-tRNA substrates, and thus capable of contributing to the fidelity mechanism.
Collapse
Affiliation(s)
- Dylan Girodat
- Theoretical Biology and Biophysics Group, Theoretical Division, Los Alamos National Laboratory, Los Alamos, NM, 87545, USA
| | - Scott C Blanchard
- Department of Structural Biology, St. Jude Children's Research Hospital, Memphis, TN, USA
| | - Hans-Joachim Wieden
- Department of Chemistry and Biochemistry, Alberta RNA Research and Training Institute, University of Lethbridge, Lethbridge, AB T1K 3M4, Canada
| | - Karissa Y Sanbonmatsu
- Theoretical Biology and Biophysics Group, Theoretical Division, Los Alamos National Laboratory, Los Alamos, NM, 87545, USA; New Mexico Consortium, Los Alamos, NM, 87544.
| |
Collapse
|
18
|
Timr S, Madern D, Sterpone F. Protein thermal stability. PROGRESS IN MOLECULAR BIOLOGY AND TRANSLATIONAL SCIENCE 2020; 170:239-272. [PMID: 32145947 DOI: 10.1016/bs.pmbts.2019.12.007] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Abstract
Proteins, in general, fold to a well-organized three-dimensional structure in order to function. The stability of this functional shape can be perturbed by external environmental conditions, such as temperature. Understanding the molecular factors underlying the resistance of proteins to the thermal stress has important consequences. First of all, it can aid the design of thermostable enzymes able to perform efficient catalysis in the high-temperature regime. Second, it is an essential brick of knowledge required to decipher the evolutionary pathways of life adaptation on Earth. Thanks to the development of atomistic simulations and ad hoc enhanced sampling techniques, it is now possible to investigate this problem in silico, and therefore provide support to experiments. After having described the methodological aspects, the chapter proposes an extended discussion on two problems. First, we focus on thermophilic proteins, a perfect model to address the issue of thermal stability and molecular evolution. Second, we discuss the issue of how protein thermal stability is affected by crowded in vivo-like conditions.
Collapse
Affiliation(s)
- Stepan Timr
- CNRS, Université de Paris, UPR 9080, Laboratoire de Biochimie Théorique, Paris, France; Institut de Biologie Physico-Chimique-Fondation Edmond de Rothschild, PSL Research University, Paris, France
| | | | - Fabio Sterpone
- CNRS, Université de Paris, UPR 9080, Laboratoire de Biochimie Théorique, Paris, France; Institut de Biologie Physico-Chimique-Fondation Edmond de Rothschild, PSL Research University, Paris, France.
| |
Collapse
|
19
|
Thomas A, Cutlan R, Finnigan W, van der Giezen M, Harmer N. Highly thermostable carboxylic acid reductases generated by ancestral sequence reconstruction. Commun Biol 2019; 2:429. [PMID: 31799431 PMCID: PMC6874671 DOI: 10.1038/s42003-019-0677-y] [Citation(s) in RCA: 27] [Impact Index Per Article: 5.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/28/2019] [Accepted: 11/04/2019] [Indexed: 12/19/2022] Open
Abstract
Carboxylic acid reductases (CARs) are biocatalysts of industrial importance. Their properties, especially their poor stability, render them sub-optimal for use in a bioindustrial pipeline. Here, we employed ancestral sequence reconstruction (ASR) - a burgeoning engineering tool that can identify stabilizing but enzymatically neutral mutations throughout a protein. We used a three-algorithm approach to reconstruct functional ancestors of the Mycobacterial and Nocardial CAR1 orthologues. Ancestral CARs (AncCARs) were confirmed to be CAR enzymes with a preference for aromatic carboxylic acids. Ancestors also showed varied tolerances to solvents, pH and in vivo-like salt concentrations. Compared to well-studied extant CARs, AncCARs had a Tm up to 35 °C higher, with half-lives up to nine times longer than the greatest previously observed. Using ancestral reconstruction we have expanded the existing CAR toolbox with three new thermostable CAR enzymes, providing access to the high temperature biosynthesis of aldehydes to drive new applications in biocatalysis.
Collapse
Affiliation(s)
- Adam Thomas
- Living Systems Institute, Stocker Road, Exeter, EX4 4QD UK
- Present Address: Department of Biosciences, Geoffrey Pope Building, Stocker Road, Exeter, EX4 4QD UK
| | - Rhys Cutlan
- Living Systems Institute, Stocker Road, Exeter, EX4 4QD UK
- Present Address: Department of Biosciences, Geoffrey Pope Building, Stocker Road, Exeter, EX4 4QD UK
| | - William Finnigan
- Present Address: Department of Biosciences, Geoffrey Pope Building, Stocker Road, Exeter, EX4 4QD UK
| | - Mark van der Giezen
- Present Address: Department of Biosciences, Geoffrey Pope Building, Stocker Road, Exeter, EX4 4QD UK
- Centre for Organelle Research, University of Stavanger, Richard Johnsens gate 4, Stavanger, 4021 Norway
| | - Nicholas Harmer
- Living Systems Institute, Stocker Road, Exeter, EX4 4QD UK
- Present Address: Department of Biosciences, Geoffrey Pope Building, Stocker Road, Exeter, EX4 4QD UK
| |
Collapse
|
20
|
Abstract
Classically, phenotype is what is observed, and genotype is the genetic makeup. Statistical studies aim to project phenotypic likelihoods of genotypic patterns. The traditional genotype-to-phenotype theory embraces the view that the encoded protein shape together with gene expression level largely determines the resulting phenotypic trait. Here, we point out that the molecular biology revolution at the turn of the century explained that the gene encodes not one but ensembles of conformations, which in turn spell all possible gene-associated phenotypes. The significance of a dynamic ensemble view is in understanding the linkage between genetic change and the gained observable physical or biochemical characteristics. Thus, despite the transformative shift in our understanding of the basis of protein structure and function, the literature still commonly relates to the classical genotype-phenotype paradigm. This is important because an ensemble view clarifies how even seemingly small genetic alterations can lead to pleiotropic traits in adaptive evolution and in disease, why cellular pathways can be modified in monogenic and polygenic traits, and how the environment may tweak protein function.
Collapse
Affiliation(s)
- Ruth Nussinov
- Cancer and Inflammation Program, Leidos Biomedical Research, Inc., Frederick National Laboratory for Cancer Research, National Cancer Institute at Frederick, Frederick, Maryland, United States of America
- Sackler Institute of Molecular Medicine, Department of Human Genetics and Molecular Medicine, Sackler School of Medicine, Tel Aviv University, Tel Aviv, Israel
| | - Chung-Jung Tsai
- Cancer and Inflammation Program, Leidos Biomedical Research, Inc., Frederick National Laboratory for Cancer Research, National Cancer Institute at Frederick, Frederick, Maryland, United States of America
| | - Hyunbum Jang
- Cancer and Inflammation Program, Leidos Biomedical Research, Inc., Frederick National Laboratory for Cancer Research, National Cancer Institute at Frederick, Frederick, Maryland, United States of America
| |
Collapse
|
21
|
Trudeau DL, Tawfik DS. Protein engineers turned evolutionists-the quest for the optimal starting point. Curr Opin Biotechnol 2019; 60:46-52. [PMID: 30611116 DOI: 10.1016/j.copbio.2018.12.002] [Citation(s) in RCA: 73] [Impact Index Per Article: 14.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/17/2018] [Revised: 11/22/2018] [Accepted: 12/03/2018] [Indexed: 12/12/2022]
Abstract
The advent of laboratory directed evolution yielded a fruitful crosstalk between the disciplines of molecular evolution and bio-engineering. Here, we outline recent developments in both disciplines with respect to how one can identify the best starting points for directed evolution, such that highly efficient and robust tailor-made enzymes can be obtained with minimal optimization. Directed evolution studies have highlighted essential features of engineer-able enzymes: highly stable, mutationally robust enzymes with the capacity to accept a broad range of substrates. Robust, evolvable enzymes can be inferred from the natural sequence record. Broad substrate spectrum relates to conformational plasticity and can also be predicted by phylogenetic analyses and/or by computational design. Overall, an increasingly powerful toolkit is becoming available for identifying optimal starting points including network analyses of enzyme superfamilies and other bioinformatics methods.
Collapse
Affiliation(s)
- Devin L Trudeau
- Department of Biomolecular Sciences, Weizmann Institute of Science, 234 Herzl Street, Rehovot 7610001, Israel
| | - Dan S Tawfik
- Department of Biomolecular Sciences, Weizmann Institute of Science, 234 Herzl Street, Rehovot 7610001, Israel.
| |
Collapse
|