1
|
Magi Meconi G, Sasselli IR, Bianco V, Onuchic JN, Coluzza I. Key aspects of the past 30 years of protein design. REPORTS ON PROGRESS IN PHYSICS. PHYSICAL SOCIETY (GREAT BRITAIN) 2022; 85:086601. [PMID: 35704983 DOI: 10.1088/1361-6633/ac78ef] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/03/2021] [Accepted: 06/15/2022] [Indexed: 06/15/2023]
Abstract
Proteins are the workhorse of life. They are the building infrastructure of living systems; they are the most efficient molecular machines known, and their enzymatic activity is still unmatched in versatility by any artificial system. Perhaps proteins' most remarkable feature is their modularity. The large amount of information required to specify each protein's function is analogically encoded with an alphabet of just ∼20 letters. The protein folding problem is how to encode all such information in a sequence of 20 letters. In this review, we go through the last 30 years of research to summarize the state of the art and highlight some applications related to fundamental problems of protein evolution.
Collapse
Affiliation(s)
- Giulia Magi Meconi
- Computational Biophysics Lab, Center for Cooperative Research in Biomaterials (CIC biomaGUNE), Basque Research and Technology Alliance (BRTA), Paseo de Miramon 182, 20014, Donostia-San Sebastián, Spain
| | - Ivan R Sasselli
- Computational Biophysics Lab, Center for Cooperative Research in Biomaterials (CIC biomaGUNE), Basque Research and Technology Alliance (BRTA), Paseo de Miramon 182, 20014, Donostia-San Sebastián, Spain
| | | | - Jose N Onuchic
- Center for Theoretical Biological Physics, Department of Physics & Astronomy, Department of Chemistry, Department of Biosciences, Rice University, Houston, TX 77251, United States of America
| | - Ivan Coluzza
- BCMaterials, Basque Center for Materials, Applications and Nanostructures, Bld. Martina Casiano, UPV/EHU Science Park, Barrio Sarriena s/n, 48940 Leioa, Spain
- Basque Foundation for Science, Ikerbasque, 48009, Bilbao, Spain
| |
Collapse
|
2
|
Nerattini F, Figliuzzi M, Cardelli C, Tubiana L, Bianco V, Dellago C, Coluzza I. Identification of Protein Functional Regions. Chemphyschem 2020; 21:335-347. [PMID: 31944517 DOI: 10.1002/cphc.201900898] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/15/2019] [Revised: 11/01/2019] [Indexed: 11/12/2022]
Abstract
Protein sequence stores the information relative to both functionality and stability, thus making it difficult to disentangle the two contributions. However, the identification of critical residues for function and stability has important implications for the mapping of the proteome interactions, as well as for many pharmaceutical applications, e. g. the identification of ligand binding regions for targeted pharmaceutical protein design. In this work, we propose a computational method to identify critical residues for protein functionality and stability and to further categorise them in strictly functional, structural and intermediate. We evaluate single site conservation and use Direct Coupling Analysis (DCA) to identify co-evolved residues both in natural and artificial evolution processes. We reproduce artificial evolution using protein design and base our approach on the hypothesis that artificial evolution in the absence of any functional constraint would exclusively lead to site conservation and co-evolution events of the structural type. Conversely, natural evolution intrinsically embeds both functional and structural information. By comparing the lists of conserved and co-evolved residues, outcomes of the analysis on natural and artificial evolution, we identify the functional residues without the need of any a priori knowledge of the biological role of the analysed protein.
Collapse
Affiliation(s)
- Francesca Nerattini
- Faculty of Physics, University of Vienna, Boltzmanngasse 5, 1090, Vienna, Austria
| | - Matteo Figliuzzi
- Sorbonne Universites, UPMC, Institut de Biologie Paris-Seine, CNRS, Laboratoire de Biologie Computationnelle et Quantitative UMR, 7238, Paris, France
| | - Chiara Cardelli
- Faculty of Physics, University of Vienna, Boltzmanngasse 5, 1090, Vienna, Austria
| | - Luca Tubiana
- Physics Department, Universitá degli studi di Trento, via Sommarive 14, 38123, Trento, IT
| | - Valentino Bianco
- Faculty of Physics, University of Vienna, Boltzmanngasse 5, 1090, Vienna, Austria.,Faculty of Chemistry, Chemical Physics Department, Universidad Complutense de Madrid, Plaza de las Ciencias, Ciudad Universitaria, Madrid, 28040, Spain
| | - Christoph Dellago
- Faculty of Physics, University of Vienna, Boltzmanngasse 5, 1090, Vienna, Austria
| | - Ivan Coluzza
- CIC biomaGUNE, Paseo Miramon 182, 20014 San Sebastian, Spain, and IKERBASQUE, Basque Foundation for Science, 48013, Bilbao, Spain
| |
Collapse
|
3
|
Magyar C, Gromiha MM, Sávoly Z, Simon I. The role of stabilization centers in protein thermal stability. Biochem Biophys Res Commun 2016; 471:57-62. [DOI: 10.1016/j.bbrc.2016.01.181] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/14/2016] [Accepted: 01/30/2016] [Indexed: 11/25/2022]
|
4
|
Balcão VM, Vila MMDC. Structural and functional stabilization of protein entities: state-of-the-art. Adv Drug Deliv Rev 2015; 93:25-41. [PMID: 25312675 DOI: 10.1016/j.addr.2014.10.005] [Citation(s) in RCA: 146] [Impact Index Per Article: 16.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/04/2013] [Revised: 08/03/2014] [Accepted: 10/01/2014] [Indexed: 12/13/2022]
Abstract
Within the context of biomedicine and pharmaceutical sciences, the issue of (therapeutic) protein stabilization assumes particular relevance. Stabilization of protein and protein-like molecules translates into preservation of both structure and functionality during storage and/or targeting, and such stabilization is mostly attained through establishment of a thermodynamic equilibrium with the (micro)environment. The basic thermodynamic principles that govern protein structural transitions and the interactions of the protein molecule with its (micro)environment are, therefore, tackled in a systematic fashion. Highlights are given to the major classes of (bio)therapeutic molecules, viz. enzymes, recombinant proteins, (macro)peptides, (monoclonal) antibodies and bacteriophages. Modification of the microenvironment of the biomolecule via multipoint covalent attachment onto a solid surface followed by hydrophilic polymer co-immobilization, or physical containment within nanocarriers, are some of the (latest) strategies discussed aiming at full structural and functional stabilization of said biomolecules.
Collapse
Affiliation(s)
- Victor M Balcão
- LaBNUS - Biomaterials and Nanotechnology Laboratory, i(bs)(2) - intelligent biosensing and biomolecule stabilization research group, University of Sorocaba, Sorocaba, SP, Brazil; CEB - Centre of Biological Engineering, University of Minho, Braga, Portugal.
| | - Marta M D C Vila
- LaBNUS - Biomaterials and Nanotechnology Laboratory, i(bs)(2) - intelligent biosensing and biomolecule stabilization research group, University of Sorocaba, Sorocaba, SP, Brazil
| |
Collapse
|
5
|
Relationship between Metabolic Fluxes and Sequence-Derived Properties of Enzymes. INTERNATIONAL SCHOLARLY RESEARCH NOTICES 2014; 2014:817102. [PMID: 27437461 PMCID: PMC4897147 DOI: 10.1155/2014/817102] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 04/16/2014] [Accepted: 08/24/2014] [Indexed: 11/17/2022]
Abstract
Metabolic fluxes are key parameters of metabolic pathways being closely related to the kinetic properties of enzymes, thereby could be dependent on. This study examines possible relationships between the metabolic fluxes and the physical-chemical/structural features of enzymes from the yeast Saccharomyces cerevisiae glycolysis pathway. Metabolic fluxes were quantified by the COPASI tool using the kinetic models of Hynne and Teusink at varied concentrations of external glucose. The enzyme sequences were taken from the UniProtKB and the average amino acid (AA) properties were computed using the set of Georgiev's uncorrelated scales that satisfy the VARIMAX criterion and specific AA indices that show the highest correlations with those. Multiple linear regressions (88.41% <R adjusted (2) < 93.32%; P < 0.00001) were found between the values of metabolic fluxes and the selected sets of the average AA properties. The hydrophobicity, α-helicity, and net charge were pointed out as the most influential characteristics of the sequences. The results provide an evidence that metabolic fluxes of the yeast glycolysis pathway are closely related to certain physical-chemical properties of relevant enzymes and support the view on the interdependence of catalytic, binding, and structural AA residues to ensure the efficiency of biocatalysts and, hence, physiologically adequate metabolic processes.
Collapse
|
6
|
Zikmanis P, Kampenusa I. Relationships between kinetic constants and the amino acid composition of enzymes from the yeast Saccharomyces cerevisiae glycolysis pathway. EURASIP JOURNAL ON BIOINFORMATICS & SYSTEMS BIOLOGY 2012; 2012:11. [PMID: 22867018 PMCID: PMC3494524 DOI: 10.1186/1687-4153-2012-11] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 01/11/2012] [Accepted: 07/21/2012] [Indexed: 01/23/2023]
Abstract
The kinetic models of metabolic pathways represent a system of biochemical reactions in terms of metabolic fluxes and enzyme kinetics. Therefore, the apparent differences of metabolic fluxes might reflect distinctive kinetic characteristics, as well as sequence-dependent properties of the employed enzymes. This study aims to examine possible linkages between kinetic constants and the amino acid (AA) composition (AAC) for enzymes from the yeast Saccharomyces cerevisiae glycolytic pathway. The values of Michaelis-Menten constant (KM), turnover number (kcat), and specificity constant (ksp = kcat/KM) were taken from BRENDA (15, 17, and 16 values, respectively) and protein sequences of nine enzymes (HXK, GADH, PGK, PGM, ENO, PK, PDC, TIM, and PYC) from UniProtKB. The AAC and sequence properties were computed by ExPASy/ProtParam tool and data processed by conventional methods of multivariate statistics. Multiple linear regressions were found between the log-values of kcat (3 models, 85.74% < Radj.2 <94.11%, p < 0.00001), KM (1 model, Radj.2 = 96.70%, p < 0.00001), ksp (3 models, 96.15% < Radj.2 < 96.50%, p < 0.00001), and the sets of AA frequencies (four to six for each model) selected from enzyme sequences while assessing the potential multicollinearity between variables. It was also found that the selection of independent variables in multiple regression models may reflect certain advantages for definite AA physicochemical and structural propensities, which could affect the properties of sequences. The results support the view on the actual interdependence of catalytic, binding, and structural residues to ensure the efficiency of biocatalysts, since the kinetic constants of the yeast enzymes appear as closely related to the overall AAC of sequences.
Collapse
Affiliation(s)
- Peteris Zikmanis
- Institute of Microbiology and Biotechnology, University of Latvia, Kronvalda Boulevard 4, Riga, LV-1010, Latvia
| | - Inara Kampenusa
- Institute of Microbiology and Biotechnology, University of Latvia, Kronvalda Boulevard 4, Riga, LV-1010, Latvia
| |
Collapse
|
7
|
Ágoston V, Cemazar M, Kaján L, Pongor S. Graph-representation of oxidative folding pathways. BMC Bioinformatics 2005; 6:19. [PMID: 15676070 PMCID: PMC549202 DOI: 10.1186/1471-2105-6-19] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/20/2004] [Accepted: 01/27/2005] [Indexed: 11/17/2022] Open
Abstract
Background The process of oxidative folding combines the formation of native disulfide bond with conformational folding resulting in the native three-dimensional fold. Oxidative folding pathways can be described in terms of disulfide intermediate species (DIS) which can also be isolated and characterized. Each DIS corresponds to a family of folding states (conformations) that the given DIS can adopt in three dimensions. Results The oxidative folding space can be represented as a network of DIS states interconnected by disulfide interchange reactions that can either create/abolish or rearrange disulfide bridges. We propose a simple 3D representation wherein the states having the same number of disulfide bridges are placed on separate planes. In this representation, the shuffling transitions are within the planes, and the redox edges connect adjacent planes. In a number of experimentally studied cases (bovine pancreatic trypsin inhibitor, insulin-like growth factor and epidermal growth factor), the observed intermediates appear as part of contiguous oxidative folding pathways. Conclusions Such networks can be used to visualize folding pathways in terms of the experimentally observed intermediates. A simple visualization template written for the Tulip package can be obtained from V.A.
Collapse
Affiliation(s)
- Vilmos Ágoston
- Bioinformatics Group, Biological Research Center, Hungarian Academy of Sciences, Temesvári krt. 62, 6726 Szeged, Hungary
| | - Masa Cemazar
- Protein Structure and Bioinformatics Group, International Centre for Genetic Engineering and Biotechnology, Area Science Park, 34012 Trieste, Italy
- Institute for Molecular Bioscience, University of Queensland, St. Lucia 4072, QLD, Australia
| | - László Kaján
- Protein Structure and Bioinformatics Group, International Centre for Genetic Engineering and Biotechnology, Area Science Park, 34012 Trieste, Italy
| | - Sándor Pongor
- Protein Structure and Bioinformatics Group, International Centre for Genetic Engineering and Biotechnology, Area Science Park, 34012 Trieste, Italy
| |
Collapse
|