1
|
Freiberger MI, Ruiz-Serra V, Pontes C, Romero-Durana M, Galaz-Davison P, Ramírez-Sarmiento CA, Schuster CD, Marti MA, Wolynes PG, Ferreiro DU, Parra RG, Valencia A. Local energetic frustration conservation in protein families and superfamilies. Nat Commun 2023; 14:8379. [PMID: 38104123 PMCID: PMC10725452 DOI: 10.1038/s41467-023-43801-2] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/23/2023] [Accepted: 11/20/2023] [Indexed: 12/19/2023] Open
Abstract
Energetic local frustration offers a biophysical perspective to interpret the effects of sequence variability on protein families. Here we present a methodology to analyze local frustration patterns within protein families and superfamilies that allows us to uncover constraints related to stability and function, and identify differential frustration patterns in families with a common ancestry. We analyze these signals in very well studied protein families such as PDZ, SH3, ɑ and β globins and RAS families. Recent advances in protein structure prediction make it possible to analyze a vast majority of the protein space. An automatic and unsupervised proteome-wide analysis on the SARS-CoV-2 virus demonstrates the potential of our approach to enhance our understanding of the natural phenotypic diversity of protein families beyond single protein instances. We apply our method to modify biophysical properties of natural proteins based on their family properties, as well as perform unsupervised analysis of large datasets to shed light on the physicochemical signatures of poorly characterized proteins such as the ones belonging to emergent pathogens.
Collapse
Affiliation(s)
- Maria I Freiberger
- Laboratorio de Fisiología de Proteínas, Departamento de Química Biológica - IQUIBICEN/CONICET, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Buenos Aires, C1428EGA, Argentina
| | - Victoria Ruiz-Serra
- Computational Biology Group, Life Sciences Department, Barcelona Supercomputing Center, Barcelona, Spain
| | - Camila Pontes
- Computational Biology Group, Life Sciences Department, Barcelona Supercomputing Center, Barcelona, Spain
| | - Miguel Romero-Durana
- Computational Biology Group, Life Sciences Department, Barcelona Supercomputing Center, Barcelona, Spain
| | - Pablo Galaz-Davison
- Institute for Biological and Medical Engineering, Schools of Engineering, Medicine, and Biological Sciences, Pontificia Universidad Católica de Chile, Santiago, 7820436, Chile
- ANID - Millennium Science Initiative Program - Millennium Institute for Integrative Biology (iBio), Santiago, 8331150, Chile
| | - Cesar A Ramírez-Sarmiento
- Institute for Biological and Medical Engineering, Schools of Engineering, Medicine, and Biological Sciences, Pontificia Universidad Católica de Chile, Santiago, 7820436, Chile
- ANID - Millennium Science Initiative Program - Millennium Institute for Integrative Biology (iBio), Santiago, 8331150, Chile
| | - Claudio D Schuster
- Laboratorio de Bioinformática, Departamento de Química Biológica - IQUIBICEN/CONICET, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, C1428EGA, Buenos Aires, Argentina
| | - Marcelo A Marti
- Laboratorio de Bioinformática, Departamento de Química Biológica - IQUIBICEN/CONICET, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, C1428EGA, Buenos Aires, Argentina
| | - Peter G Wolynes
- Center for Theoretical Biological Physics and Department of Chemistry, Rice University, Houston, TX, 77005, USA
| | - Diego U Ferreiro
- Laboratorio de Fisiología de Proteínas, Departamento de Química Biológica - IQUIBICEN/CONICET, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Buenos Aires, C1428EGA, Argentina
| | - R Gonzalo Parra
- Computational Biology Group, Life Sciences Department, Barcelona Supercomputing Center, Barcelona, Spain.
| | - Alfonso Valencia
- Computational Biology Group, Life Sciences Department, Barcelona Supercomputing Center, Barcelona, Spain
- Catalan Institution for Research and Advanced Studies (ICREA), Barcelona, Spain
| |
Collapse
|
2
|
Oda M. Analysis of the Structural Dynamics of Proteins in the Ligand-Unbound and -Bound States by Diffracted X-ray Tracking. Int J Mol Sci 2023; 24:13717. [PMID: 37762021 PMCID: PMC10531450 DOI: 10.3390/ijms241813717] [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: 08/17/2023] [Revised: 09/03/2023] [Accepted: 09/05/2023] [Indexed: 09/29/2023] Open
Abstract
Although many protein structures have been determined at atomic resolution, the majority of them are static and represent only the most stable or averaged structures in solution. When a protein binds to its ligand, it usually undergoes fluctuation and changes its conformation. One attractive method for obtaining an accurate view of proteins in solution, which is required for applications such as the rational design of proteins and structure-based drug design, is diffracted X-ray tracking (DXT). DXT can detect the protein structural dynamics on a timeline via gold nanocrystals attached to the protein. Here, the structure dynamics of single-chain Fv antibodies, helix bundle-forming de novo designed proteins, and DNA-binding proteins in both ligand-unbound and ligand-bound states were analyzed using the DXT method. The resultant mean square angular displacements (MSD) curves in both the tilting and twisting directions clearly demonstrated that structural fluctuations were suppressed upon ligand binding, and the binding energies determined using the angular diffusion coefficients from the MSD agreed well with the binding thermodynamics determined using isothermal titration calorimetry. In addition, the size of gold nanocrystals is discussed, which is one of the technical concerns of DXT.
Collapse
Affiliation(s)
- Masayuki Oda
- Graduate School of Life and Environmental Sciences, Kyoto Prefectural University, 1-5 Hangi-cho, Shimogamo, Sakyo-ku, Kyoto 606-8522, Japan
| |
Collapse
|
3
|
Banerjee A, Gosavi S. Potential Self-Peptide Inhibitors of the SARS-CoV-2 Main Protease. J Phys Chem B 2023; 127:855-865. [PMID: 36689738 PMCID: PMC9883841 DOI: 10.1021/acs.jpcb.2c05917] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/17/2022] [Revised: 11/23/2022] [Indexed: 01/24/2023]
Abstract
The SARS-CoV-2 main protease (Mpro) plays an essential role in viral replication, cleaving viral polyproteins into functional proteins. This makes Mpro an important drug target. Mpro consists of an N-terminal catalytic domain and a C-terminal α-helical domain (MproC). Previous studies have shown that peptides derived from a given protein sequence (self-peptides) can affect the folding and, in turn, the function of that protein. Since the SARS-CoV-1 MproC is known to stabilize its Mpro and regulate its function, we hypothesized that SARS-CoV-2 MproC-derived self-peptides may modulate the folding and the function of SARS-CoV-2 Mpro. To test this, we studied the folding of MproC in the presence of various self-peptides using coarse-grained structure-based models and molecular dynamics simulations. In these simulations of MproC and one self-peptide, we found that two self-peptides, the α1-helix and the loop between α4 and α5 (loop4), could replace the equivalent native sequences in the MproC structure. Replacement of either sequence in full-length Mpro should, in principle, be able to perturb Mpro function albeit through different mechanisms. Some general principles for the rational design of self-peptide inhibitors emerge: The simulations show that prefolded self-peptides are more likely to replace native sequences than those which do not possess structure. Additionally, the α1-helix self-peptide is kinetically stable and once inserted rarely exchanges with the native α1-helix, while the loop4 self-peptide is easily replaced by the native loop4, making it less useful for modulating function. In summary, a prefolded α1-derived peptide should be able to inhibit SARS-CoV-2 Mpro function.
Collapse
Affiliation(s)
- Arkadeep Banerjee
- Simons Centre for the Study
of Living Machines, National Centre for Biological Sciences, Tata Institute of Fundamental Research, Bengaluru 560065, India
| | - Shachi Gosavi
- Simons Centre for the Study
of Living Machines, National Centre for Biological Sciences, Tata Institute of Fundamental Research, Bengaluru 560065, India
| |
Collapse
|
4
|
Gonet M, Krowarsch D, Schubert J, Tabiś A, Bania J. Stability and Resistance to Proteolysis of Enterotoxins SEC and SEL Produced by Staphylococcus epidermidis and Staphylococcus aureus. Foodborne Pathog Dis 2023; 20:32-37. [PMID: 36622956 DOI: 10.1089/fpd.2022.0059] [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: 01/10/2023] Open
Abstract
The only staphylococcal enterotoxins produced by Staphylococcus epidermidis include SECepi and SELepi, whereas Staphylococcus aureus produces orthologous SECs and SEL having different sequences. We compared S. epidermidis and S. aureus SECs and SELs in terms of resistance to proteolysis and both, thermal and chemical stability. We show that SECepi and SELepi produced by S. epidermidis have similar resistance to proteolysis if compared with their respective orthologues produced by S. aureus. Studied S. epidermidis and S. aureus SEC variants incubated with pepsin at pH 2.0 were found to be more resistant to proteolysis than SELs. SELs turned out to be more resistant than SECs to proteolysis with trypsin at pH 8.0. SECepi was found to be more resistant to thermal denaturation if compared with its S. aureus orthologues. The S. epidermidis and S. aureus SEC variants were found to have higher thermal stability than SELs. Our data indicate that, due to their high stability, the enterotoxins SECepi and SELepi produced in food by S. epidermidis may pose a food safety risk comparable with that posed by S. aureus enterotoxins.
Collapse
Affiliation(s)
- Magdalena Gonet
- Department of Food Hygiene and Consumer Health Protection, the Wroclaw University of Environmental and Life Sciences, Wroclaw, Poland
| | - Daniel Krowarsch
- Department of Protein Biotechnology, Faculty of Biotechnology, University of Wroclaw, Wroclaw, Poland
| | | | - Aleksandra Tabiś
- Department of Food Hygiene and Consumer Health Protection, the Wroclaw University of Environmental and Life Sciences, Wroclaw, Poland
| | - Jacek Bania
- Department of Food Hygiene and Consumer Health Protection, the Wroclaw University of Environmental and Life Sciences, Wroclaw, Poland
| |
Collapse
|
5
|
Yadahalli S, Jayanthi LP, Gosavi S. A Method for Assessing the Robustness of Protein Structures by Randomizing Packing Interactions. Front Mol Biosci 2022; 9:849272. [PMID: 35832734 PMCID: PMC9271847 DOI: 10.3389/fmolb.2022.849272] [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/05/2022] [Accepted: 04/27/2022] [Indexed: 12/02/2022] Open
Abstract
Many single-domain proteins are not only stable and water-soluble, but they also populate few to no intermediates during folding. This reduces interactions between partially folded proteins, misfolding, and aggregation, and makes the proteins tractable in biotechnological applications. Natural proteins fold thus, not necessarily only because their structures are well-suited for folding, but because their sequences optimize packing and fit their structures well. In contrast, folding experiments on the de novo designed Top7 suggest that it populates several intermediates. Additionally, in de novo protein design, where sequences are designed for natural and new non-natural structures, tens of sequences still need to be tested before success is achieved. Both these issues may be caused by the specific scaffolds used in design, i.e., some protein scaffolds may be more tolerant to packing perturbations and varied sequences. Here, we report a computational method for assessing the response of protein structures to packing perturbations. We then benchmark this method using designed proteins and find that it can identify scaffolds whose folding gets disrupted upon perturbing packing, leading to the population of intermediates. The method can also isolate regions of both natural and designed scaffolds that are sensitive to such perturbations and identify contacts which when present can rescue folding. Overall, this method can be used to identify protein scaffolds that are more amenable to whole protein design as well as to identify protein regions which are sensitive to perturbations and where further mutations should be avoided during protein engineering.
Collapse
|
6
|
Du J, Yin H, Lu Y, Lu T, Chen T. Effects of Surface Tethering on the Thermodynamics and Kinetics of Frustrated Protein Folding. J Phys Chem B 2022; 126:4776-4786. [PMID: 35731862 DOI: 10.1021/acs.jpcb.2c01982] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The interaction between the protein and surface plays an important role in biology and biotechnology. To understand how surface tethering influences the folding behavior of frustrated proteins, in this work, we systematically study the thermodynamics and folding kinetics of the bacterial immunity protein Im7 and Fyn SH3 domain tethered to a surface using Langevin dynamics simulations. Upon surface tethering, the stabilization often results from the entropic effect, whereas the destabilization is usually caused by either an energetic or entropic effect. For the Fyn SH3 domain with a two-state folding manner, the influence of nonnative interactions on thermodynamic stability is not significant, while nonnative interactions can weaken the effect of surface tethering on the change in the folding rate. By contrast, for the frustrated protein Im7, depending on where the protein is tethered, the surface tethering can promote or suppress misfolding by modulating specific nonnative contacts, thereby altering the folding rate and folding mechanism. Because surface tethering can change the intrachain diffusivity of unfolding, the kinetic stability cannot be well captured by the thermodynamic stability at some tether points. This study should be helpful in general to understand how surface tethering affects the folding energy landscape of frustrated proteins.
Collapse
Affiliation(s)
- Jiang Du
- Key Laboratory of Synthetic and Natural Functional Molecule of the Ministry of Education, College of Chemistry and Materials Science, Northwest University, Xi'an 710127, P. R. China
| | - Hongmei Yin
- Key Laboratory of Synthetic and Natural Functional Molecule of the Ministry of Education, College of Chemistry and Materials Science, Northwest University, Xi'an 710127, P. R. China
| | - Yanfang Lu
- Key Laboratory of Synthetic and Natural Functional Molecule of the Ministry of Education, College of Chemistry and Materials Science, Northwest University, Xi'an 710127, P. R. China
| | - Teng Lu
- Computer Network Information Center of the Chinese Academy of Sciences, Beijing 100083, P. R. China
| | - Tao Chen
- Key Laboratory of Synthetic and Natural Functional Molecule of the Ministry of Education, College of Chemistry and Materials Science, Northwest University, Xi'an 710127, P. R. China.,Key Laboratory of Polymer Processing Engineering (South China University of Technology), Ministry of Education, Guangzhou 510641, P. R. China
| |
Collapse
|
7
|
Senga A, Numoto N, Yamashita M, Iida A, Ito N, Kawai F, Oda M. Multiple structural states of Ca2+-regulated PET hydrolase, Cut190, and its correlation with activity and stability. J Biochem 2021; 169:207-213. [PMID: 32882044 DOI: 10.1093/jb/mvaa102] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/09/2020] [Accepted: 08/18/2020] [Indexed: 12/14/2022] Open
Abstract
An enzyme, Cut190, from a thermophilic isolate, Saccharomonospora viridis AHK190 could depolymerize polyethylene terephthalate (PET). The catalytic activity and stability of Cut190 and its S226P/R228S mutant, Cut190*, are regulated by Ca2+ binding. We previously determined the crystal structures of the inactive mutant of Cut190*, Cut190*S176A, in complex with metal ions, Ca2+ and Zn2+, and substrates, monoethyl succinate and monoethyl adipate. In this study, we determined the crystal structures of another mutant of Cut190*, Cut190**, in which the three C-terminal residues of Cut190* are deleted, and the inactive mutant, Cut190**S176A, in complex with metal ions. In addition to the previously observed closed, open and engaged forms, we determined the ejecting form, which would allow the product to irreversibly dissociate, followed by proceeding to the next cycle of reaction. These multiple forms would be stable or sub-stable states of Cut190, regulated by Ca2+ binding, and would be closely correlated with the enzyme function. Upon the deletion of the C-terminal residues, we found that the thermal stability increased while retaining the activity. The increased stability could be applied for the protein engineering of Cut190 for PET depolymerization as it requires the reaction above the glass transition temperature of PET.
Collapse
Affiliation(s)
- Akane Senga
- Graduate School of Life and Environmental Sciences, Kyoto Prefectural University, 1-5 Hangi-cho, Shimogamo, Sakyo-ku, Kyoto, Kyoto 606-8522, Japan
| | - Nobutaka Numoto
- Medical Research Institute, Tokyo Medical and Dental University (TMDU), 1-5-45 Yushima Bunkyo-ku, Tokyo 113-8510, Japan
| | - Mitsuaki Yamashita
- Faculty of Agriculture, Kindai University, 3327-204 Nakamachi, Nara, Nara 631-8505, Japan
| | - Akira Iida
- Faculty of Agriculture, Kindai University, 3327-204 Nakamachi, Nara, Nara 631-8505, Japan
| | - Nobutoshi Ito
- Medical Research Institute, Tokyo Medical and Dental University (TMDU), 1-5-45 Yushima Bunkyo-ku, Tokyo 113-8510, Japan
| | | | - Masayuki Oda
- Graduate School of Life and Environmental Sciences, Kyoto Prefectural University, 1-5 Hangi-cho, Shimogamo, Sakyo-ku, Kyoto, Kyoto 606-8522, Japan
| |
Collapse
|
8
|
Shapiro DM, Ney M, Eghtesadi SA, Chilkoti A. Protein Phase Separation Arising from Intrinsic Disorder: First-Principles to Bespoke Applications. J Phys Chem B 2021; 125:6740-6759. [PMID: 34143622 DOI: 10.1021/acs.jpcb.1c01146] [Citation(s) in RCA: 32] [Impact Index Per Article: 10.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The phase separation of biomolecules has become the focus of intense research in the past decade, with a growing body of research implicating this phenomenon in essentially all biological functions, including but not limited to homeostasis, stress responses, gene regulation, cell differentiation, and disease. Excellent reviews have been published previously on the underlying physical basis of liquid-liquid phase separation (LLPS) of biological molecules (Nat. Phys. 2015, 11, 899-904) and LLPS as it occurs natively in physiology and disease (Science 2017, 357, eaaf4382; Biochemistry 2018, 57, 2479-2487; Chem. Rev. 2014, 114, 6844-6879). Here, we review how the theoretical physical basis of LLPS has been used to better understand the behavior of biomolecules that undergo LLPS in natural systems and how this understanding has also led to the development of novel synthetic systems that exhibit biomolecular phase separation, and technologies that exploit these phenomena. In part 1 of this Review, we explore the theory behind the phase separation of biomolecules and synthetic macromolecules and introduce a few notable phase-separating biomolecules. In part 2, we cover experimental and computational methods used to study phase-separating proteins and how these techniques have uncovered the mechanisms underlying phase separation in physiology and disease. Finally, in part 3, we cover the development and applications of engineered phase-separating polypeptides, ranging from control of their self-assembly to create defined supramolecular architectures to reprogramming biological processes using engineered IDPs that exhibit LLPS.
Collapse
Affiliation(s)
- Daniel Mark Shapiro
- Department of Biomedical Engineering, Duke University, Durham, North Carolina 27708, United States
| | - Max Ney
- Department of Biomedical Engineering, Duke University, Durham, North Carolina 27708, United States
| | - Seyed Ali Eghtesadi
- Department of Biomedical Engineering, Duke University, Durham, North Carolina 27708, United States
| | - Ashutosh Chilkoti
- Department of Biomedical Engineering, Duke University, Durham, North Carolina 27708, United States
| |
Collapse
|
9
|
Staś M, Broda MA, Siodłak D. Thiazole-amino acids: influence of thiazole ring on conformational properties of amino acid residues. Amino Acids 2021; 53:673-686. [PMID: 33837859 PMCID: PMC8128816 DOI: 10.1007/s00726-021-02974-0] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/05/2020] [Accepted: 03/29/2021] [Indexed: 12/29/2022]
Abstract
Post-translational modified thiazole-amino acid (Xaa-Tzl) residues have been found in macrocyclic peptides (e.g., thiopeptides and cyanobactins), which mostly inhibit protein synthesis in Gram + bacteria. Conformational study of the series of model compounds containing this structural motif with alanine, dehydroalanine, dehydrobutyrine and dehydrophenylalanine were performed using DFT method in various environments. The solid-state crystal structure conformations of thiazole-amino acid residues retrieved from the Cambridge Structural Database were also analysed. The studied structural units tend to adopt the unique semi-extended β2 conformation; which is stabilised mainly by N-H⋯NTzl hydrogen bond, and for dehydroamino acids also by π-electron conjugation. The conformational preferences of amino acids with a thiazole ring were compared with oxazole analogues and the role of the sulfur atom in stabilising the conformations of studied peptides was discussed.
Collapse
Affiliation(s)
- Monika Staś
- Faculty of Chemistry, University of Opole, 45-052, Opole, Poland.
- Institute of Organic Chemistry and Biochemistry of the Czech Academy of Science, Flemingovo Náměstí 2, 166 10, Praha 6, Czech Republic.
| | | | - Dawid Siodłak
- Faculty of Chemistry, University of Opole, 45-052, Opole, Poland.
| |
Collapse
|
10
|
Bastida A, Zúñiga J, Requena A, Miguel B, Cerezo J. On the Role of Entropy in the Stabilization of α-Helices. J Chem Inf Model 2020; 60:6523-6531. [PMID: 33280379 DOI: 10.1021/acs.jcim.0c01177] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Protein folding evolves by exploring the conformational space with a subtle balance between enthalpy and entropy changes which eventually leads to a decrease of free energy upon reaching the folded structure. A complete understanding of this process requires, therefore, a deep insight into both contributions to free energy. In this work, we clarify the role of entropy in favoring the stabilization of folded structures in polyalanine peptides with up to 12 residues. We use a novel method referred to as K2V that allows us to obtain the potential-energy landscapes in terms of residue conformations extracted from molecular dynamics simulations at conformational equilibrium and yields folding thermodynamic magnitudes, which are in agreement with the experimental data available. Our results demonstrate that the folded structures of the larger polyalanine chains are stabilized with respect to the folded structures of the shorter chains by both an energetic contribution coming from the formation of the intramolecular hydrogen bonds and an entropic contribution coming from an increase of the entropy of the solvent with approximate weights of 60 and 40%, respectively, thus unveiling a key piece in the puzzle of protein folding. In addition, the ability of the K2V method to provide the enthalpic and entropic contributions for individual residues along the peptide chain makes it clear that the energetic and entropic stabilizations are basically governed by the nearest neighbor residue conformations, with the folding propensity being rationalized in terms of triads of residues.
Collapse
Affiliation(s)
- Adolfo Bastida
- Departamento de Química Física, Universidad de Murcia, 30100 Murcia, Spain
| | - José Zúñiga
- Departamento de Química Física, Universidad de Murcia, 30100 Murcia, Spain
| | - Alberto Requena
- Departamento de Química Física, Universidad de Murcia, 30100 Murcia, Spain
| | - Beatriz Miguel
- Departamento de Ingeniería Química y Ambiental, Universidad Politécnica de Cartagena, 30203 Cartagena, Spain
| | - Javier Cerezo
- Departamento de Química, Universidad Autónoma de Madrid, 28049 Madrid, Spain
| |
Collapse
|
11
|
Perumalla DS, Govind G, Anjukandi P. Folding‐Unfolding Dynamics of pH‐Assisted Structures of S‐Peptide. ChemistrySelect 2020. [DOI: 10.1002/slct.202000360] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022]
Affiliation(s)
| | - Gokul Govind
- Department of Chemistry Central University of Tamil Nadu Tiruvarur India
| | - Padmesh Anjukandi
- Department of Chemistry Indian Institute of Technology Palakkad India
| |
Collapse
|
12
|
Gershenson A, Gosavi S, Faccioli P, Wintrode PL. Successes and challenges in simulating the folding of large proteins. J Biol Chem 2020; 295:15-33. [PMID: 31712314 PMCID: PMC6952611 DOI: 10.1074/jbc.rev119.006794] [Citation(s) in RCA: 28] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022] Open
Abstract
Computational simulations of protein folding can be used to interpret experimental folding results, to design new folding experiments, and to test the effects of mutations and small molecules on folding. However, whereas major experimental and computational progress has been made in understanding how small proteins fold, research on larger, multidomain proteins, which comprise the majority of proteins, is less advanced. Specifically, large proteins often fold via long-lived partially folded intermediates, whose structures, potentially toxic oligomerization, and interactions with cellular chaperones remain poorly understood. Molecular dynamics based folding simulations that rely on knowledge of the native structure can provide critical, detailed information on folding free energy landscapes, intermediates, and pathways. Further, increases in computational power and methodological advances have made folding simulations of large proteins practical and valuable. Here, using serpins that inhibit proteases as an example, we review native-centric methods for simulating the folding of large proteins. These synergistic approaches range from Gō and related structure-based models that can predict the effects of the native structure on folding to all-atom-based methods that include side-chain chemistry and can predict how disease-associated mutations may impact folding. The application of these computational approaches to serpins and other large proteins highlights the successes and limitations of current computational methods and underscores how computational results can be used to inform experiments. These powerful simulation approaches in combination with experiments can provide unique insights into how large proteins fold and misfold, expanding our ability to predict and manipulate protein folding.
Collapse
Affiliation(s)
- Anne Gershenson
- Department of Biochemistry and Molecular Biology, University of Massachusetts, Amherst, Massachusetts 01003; Molecular and Cellular Biology Graduate Program, University of Massachusetts, Amherst, Massachusetts 01003.
| | - Shachi Gosavi
- Simons Centre for the Study of Living Machines, National Centre for Biological Sciences, Tata Institute of Fundamental Research, Bangalore-560065, India.
| | - Pietro Faccioli
- Dipartimento di Fisica, Universitá degli Studi di Trento, 38122 Povo (Trento), Italy; Trento Institute for Fundamental Physics and Applications, 38123 Povo (Trento), Italy.
| | - Patrick L Wintrode
- Department of Pharmaceutical Sciences, University of Maryland School of Pharmacy, Baltimore, Maryland 21201.
| |
Collapse
|
13
|
Abstract
In nature, DNA molecules carry the hereditary information. But DNA has physical and chemical properties that make it attractive for uses beyond heredity. In this Review, we discuss the potential of DNA for creating machines that are both encoded by and built from DNA molecules. We review the main methods of DNA nanostructure assembly, describe recent advances in building increasingly complex molecular structures and discuss strategies for creating machine-like nanostructures that can be actuated and move. We highlight opportunities for applications of custom DNA nanostructures as scientific tools to address challenges across biology, chemistry and engineering.
Collapse
|
14
|
Bastida A, Zúñiga J, Requena A, Cerezo J. Energetic Self-Folding Mechanism in α-Helices. J Phys Chem B 2019; 123:8186-8194. [PMID: 31498638 DOI: 10.1021/acs.jpcb.9b05860] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/24/2023]
Abstract
A novel energetic route driving the folding of a polyalanine peptide from an extended conformation to its α-helix native conformation is described, supported by a new method to compute mean potential energy surfaces accurately in terms of the dihedral angles of the peptide chain from extensive molecular dynamics simulations. The energetic self-folding (ESF) route arises specifically from the balance between the intrinsic propensity of alanine residues toward the αR conformation and two, opposite, effects: the destabilizing interaction with neighbor residues and the stabilizing formation of native hydrogen bonds, with the latter being dominant for large peptide lengths. The ESF mechanism provides simple but robust support to the nucleation-elongation or zipper models and offers a quantitative energetic funnel picture of the folding process. The mechanism is validated by the reasonable agreement between the computed folding energies and the experimental values.
Collapse
Affiliation(s)
- Adolfo Bastida
- Departamento de Química Física , Universidad de Murcia , 30100 Murcia , Spain
| | - José Zúñiga
- Departamento de Química Física , Universidad de Murcia , 30100 Murcia , Spain
| | - Alberto Requena
- Departamento de Química Física , Universidad de Murcia , 30100 Murcia , Spain
| | - Javier Cerezo
- Departamento de Química Física , Universidad de Murcia , 30100 Murcia , Spain
| |
Collapse
|
15
|
Bayat M, Karami L, Gourabi H, Ahmad F, Dormiani K, Nasr Esfahani MH, Saboury AA. Stabilizing osmolytes' effects on the structure, stability and function of tc-tenecteplase: A one peptide bond digested form of tenecteplase. Int J Biol Macromol 2019; 130:863-877. [PMID: 30849467 DOI: 10.1016/j.ijbiomac.2019.03.035] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/20/2018] [Revised: 02/20/2019] [Accepted: 03/04/2019] [Indexed: 01/10/2023]
Abstract
Organic osmolytes, as major cellular compounds, cause protein stabilization in the native form. In the present study, the possible chaperone effects of the three naturally occurring osmolytes on the two-chain form of tenecteplase (tc-TNK), a recombinant, genetically engineered mutant tissue plasminogen activator, have been explored by using circular dichroism, steady-state fluorescence, UV-Visible spectroscopy, and in silico experiments. The tc-TNK is derived from the one-chain protein upon disruption of one peptide bond. Thermal denaturation experiments showed a slightly more stabilizing effect of the three co-solvents on the single-chain TNK (sc-TNK) in comparison to that on tc-TNK. Unlike single-chain tenecteplase, the two-chain form undergoes reversible denaturation which is somehow perturbed in some cases as the result of the presence of osmolytes. Very minor changes in the secondary structure and the tertiary structure were observed. The molecular dynamics simulations and comparative structural analysis of catalytic domain of the protein in the single-chain and two-chain forms in pure water, mannitol/water, trehalose/water, and sucrose/water showed that while the stabilizing effect of the three osmolytes on tc-TNK might be induced by preferential accumulation of these molecules around the nonpolar and aromatic residues, that is to say, fewer water-hydrophobic residues' interactions in tc-TNK, sc-TNK is stabilized by preferential exclusion effect.
Collapse
Affiliation(s)
- Mahdieh Bayat
- Department of Genetics at Reproductive Biomedicine Research Center, Royan Institute, Tehran, Iran
| | - Leila Karami
- Department of Cell and Molecular Biology, Faculty of Biological Sciences, Kharazmi University, Tehran, Iran
| | - Hamid Gourabi
- Department of Genetics at Reproductive Biomedicine Research Center, Royan Institute, Tehran, Iran
| | - Faizan Ahmad
- Centre for Interdisciplinary Research in Basic Sciences, Jamia Millia Islamia, New Delhi 110025, India
| | - Kianoush Dormiani
- Department of Molecular Biotechnology, Cell Science Research Center, Royan Institute for Biotechnology, ACECR, Isfahan, Iran
| | - Mohammad H Nasr Esfahani
- Department of Molecular Biotechnology, Cell Science Research Center, Royan Institute for Biotechnology, ACECR, Isfahan, Iran.
| | - Ali A Saboury
- Institute of Biochemistry and Biophysics, University of Tehran, Tehran, Iran.
| |
Collapse
|
16
|
Studying ribosome dynamics with simplified models. Methods 2019; 162-163:128-140. [DOI: 10.1016/j.ymeth.2019.03.023] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/31/2019] [Revised: 03/22/2019] [Accepted: 03/23/2019] [Indexed: 12/24/2022] Open
|
17
|
Levi M, Bandarkar P, Yang H, Wang A, Mohanty U, Noel JK, Whitford PC. Using SMOG 2 to Simulate Complex Biomolecular Assemblies. Methods Mol Biol 2019; 2022:129-151. [PMID: 31396902 DOI: 10.1007/978-1-4939-9608-7_6] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Abstract
Over the last 20 years, the application of structure-based (Gō-like) models has ranged from protein folding with coarse-grained models to all-atom representations of large-scale molecular assemblies. While there are many variants that may be employed, the common feature of these models is that some (or all) of the stabilizing energetic interactions are defined based on the knowledge of a particular experimentally obtained conformation. With the generality of this approach, there was a need for a versatile computational platform for designing and implementing this class of models. To this end, the SMOG 2 software package provides an easy-to-use interface, where the user has full control of the model parameters. This software allows the user to edit XML-formatted files in order to provide definitions of new structure-based models. SMOG 2 reads these "template" files and maps the interactions onto specific structures, which are provided in PDB format. The force field files produced by SMOG 2 may then be used to perform simulations with a variety of popular molecular dynamics suites. In this chapter, we describe some of the key features of the SMOG 2 package, while providing examples and strategies for applying these techniques to complex (often large-scale) molecular assemblies, such as the ribosome.
Collapse
Affiliation(s)
- Mariana Levi
- Department of Physics, Northeastern University, Boston, MA, USA
| | | | - Huan Yang
- Department of Physics, Northeastern University, Boston, MA, USA
| | - Ailun Wang
- Department of Chemistry, Boston College, Chestnut Hill, MA, USA
| | - Udayan Mohanty
- Department of Chemistry, Boston College, Chestnut Hill, MA, USA
| | - Jeffrey K Noel
- Max Delbrueck Center for Molecular Medicine, Kristallographie, Berlin, Germany
| | - Paul C Whitford
- Department of Physics, Northeastern University, Boston, MA, USA.
| |
Collapse
|
18
|
Terse VL, Gosavi S. The Sensitivity of Computational Protein Folding to Contact Map Perturbations: The Case of Ubiquitin Folding and Function. J Phys Chem B 2018; 122:11497-11507. [PMID: 30234303 DOI: 10.1021/acs.jpcb.8b07409] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Abstract
Ubiquitin is a small model protein, commonly used in protein folding experiments and simulations. We simulated ubiquitin using a well-tested structure-based model coarse-grained to a Cα level (Cα-SBM) and found that the simulated folding route did not agree with the experimentally observed one. Simulating the Cα-SBM with a cutoff contact map, instead of a screened contact map, switched the folding route with the new route matching the experimental route. Thus, the simulated folding of ubiquitin is sensitive to contact map definition. The screened contact map, which is used in folding simulations because it captures protein folding cooperativity, removes contacts in which the atoms in contact are occluded by a third atom and is less sensitive to the value of the cutoff distance in well-packed regions of the protein. In sparsely packed regions, the larger cutoff distance creates bridging contacts between atoms which are separated by voids. Such contacts do not seem to affect the folding of most proteins, including those of the ubiquitin fold. However, the surface of ubiquitin has several protruding functional side chains which naturally create bridging contacts. Together, our results show that subtle structural features of a protein that may not be apparent by mere observation can be identified by comparing folding simulations of SBMs in which these features are differently encoded. When such structural features are preserved for functional reasons, differences in computational folding can be leveraged to identify functional features. Notably, such features are accessible to a gradation of SBMs even in commonly studied proteins such as ubiquitin.
Collapse
Affiliation(s)
- Vishram L Terse
- Simons Centre for the Study of Living Machines , National Centre for Biological Sciences , Tata Institute of Fundamental Research, Bangalore 560065 , India
| | - Shachi Gosavi
- Simons Centre for the Study of Living Machines , National Centre for Biological Sciences , Tata Institute of Fundamental Research, Bangalore 560065 , India
| |
Collapse
|
19
|
Mascarenhas NM, Terse VL, Gosavi S. Intrinsic Disorder in a Well-Folded Globular Protein. J Phys Chem B 2018; 122:1876-1884. [PMID: 29304275 DOI: 10.1021/acs.jpcb.7b12546] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Abstract
The folded structure of the heterodimeric sweet protein monellin mimics single-chain proteins with topology β1-α1-β2-β3-β4-β5 (chain A: β3-β4-β5; chain B: β1-α1-β2). Furthermore, like naturally occurring single-chain proteins of a similar size, monellin folds cooperatively with no detectable intermediates. However, the two monellin chains, A and B, are marginally structured in isolation and fold only upon binding to each other. Thus, monellin presents a unique opportunity to understand the design of intrinsically disordered proteins that fold upon binding. Here, we study the folding of a single-chain variant of monellin (scMn) using simulations of an all heavy-atom structure-based model. These simulations can explain mechanistic details derived from scMn experiments performed using several different structural probes. scMn folds cooperatively in our structure-based simulations, as is also seen in experiments. We find that structure formation near the transition-state ensemble of scMn is not uniformly distributed but is localized to a hairpin-like structure which contains one strand from each chain (β2, β3). Thus, the sequence and the underlying energetics of heterodimeric monellin promote the early formation of the interchain interface (β2-β3). By studying computational scMn mutants whose "interchain" interactions are deleted, we infer that this energy distribution allows the two protein chains to remain largely disordered when this interface is not folded. From these results, we suggest that cutting the protein backbone of a globular protein between residues which lie within its folding nucleus may be one way to construct two disordered fragments which fold upon binding.
Collapse
Affiliation(s)
| | - Vishram L Terse
- Simons Centre for the Study of Living Machines, National Centre for Biological Sciences, Tata Institute of Fundamental Research , Bangalore 560065, India
| | - Shachi Gosavi
- Simons Centre for the Study of Living Machines, National Centre for Biological Sciences, Tata Institute of Fundamental Research , Bangalore 560065, India
| |
Collapse
|
20
|
On the folding of a structurally complex protein to its metastable active state. Proc Natl Acad Sci U S A 2018; 115:1998-2003. [PMID: 29343647 DOI: 10.1073/pnas.1708173115] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023] Open
Abstract
For successful protease inhibition, the reactive center loop (RCL) of the two-domain serine protease inhibitor, α1-antitrypsin (α1-AT), needs to remain exposed in a metastable active conformation. The α1-AT RCL is sequestered in a β-sheet in the stable latent conformation. Thus, to be functional, α1-AT must always fold to a metastable conformation while avoiding folding to a stable conformation. We explore the structural basis of this choice using folding simulations of coarse-grained structure-based models of the two α1-AT conformations. Our simulations capture the key features of folding experiments performed on both conformations. The simulations also show that the free energy barrier to fold to the latent conformation is much larger than the barrier to fold to the active conformation. An entropically stabilized on-pathway intermediate lowers the barrier for folding to the active conformation. In this intermediate, the RCL is in an exposed configuration, and only one of the two α1-AT domains is folded. In contrast, early conversion of the RCL into a β-strand increases the coupling between the two α1-AT domains in the transition state and creates a larger barrier for folding to the latent conformation. Thus, unlike what happens in several proteins, where separate regions promote folding and function, the structure of the RCL, formed early during folding, determines both the conformational and the functional fate of α1-AT. Further, the short 12-residue RCL modulates the free energy barrier and the folding cooperativity of the large 370-residue α1-AT. Finally, we suggest experiments to test the predicted folding mechanism for the latent state.
Collapse
|
21
|
Huynh L, Neale C, Pomès R, Chan HS. Molecular recognition and packing frustration in a helical protein. PLoS Comput Biol 2017; 13:e1005909. [PMID: 29261665 PMCID: PMC5757960 DOI: 10.1371/journal.pcbi.1005909] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/17/2017] [Revised: 01/08/2018] [Accepted: 11/28/2017] [Indexed: 01/25/2023] Open
Abstract
Biomolecular recognition entails attractive forces for the functional native states and discrimination against potential nonnative interactions that favor alternate stable configurations. The challenge posed by the competition of nonnative stabilization against native-centric forces is conceptualized as frustration. Experiment indicates that frustration is often minimal in evolved biological systems although nonnative possibilities are intuitively abundant. Much of the physical basis of minimal frustration in protein folding thus remains to be elucidated. Here we make progress by studying the colicin immunity protein Im9. To assess the energetic favorability of nonnative versus native interactions, we compute free energies of association of various combinations of the four helices in Im9 (referred to as H1, H2, H3, and H4) by extensive explicit-water molecular dynamics simulations (total simulated time > 300 μs), focusing primarily on the pairs with the largest native contact surfaces, H1-H2 and H1-H4. Frustration is detected in H1-H2 packing in that a nonnative packing orientation is significantly stabilized relative to native, whereas such a prominent nonnative effect is not observed for H1-H4 packing. However, in contrast to the favored nonnative H1-H2 packing in isolation, the native H1-H2 packing orientation is stabilized by H3 and loop residues surrounding H4. Taken together, these results showcase the contextual nature of molecular recognition, and suggest further that nonnative effects in H1-H2 packing may be largely avoided by the experimentally inferred Im9 folding transition state with native packing most developed at the H1-H4 rather than the H1-H2 interface. Biomolecules need to recognize one another with high specificity: promoting “native” functional intermolecular binding events while avoiding detrimental “nonnative” bound configurations; i.e., “frustration”—the tendency for nonnative interactions—has to be minimized. Folding of globular proteins entails a similar discrimination. To gain physical insight, we computed the binding affinities of helical structures of the protein Im9 in various native or nonnative configurations by atomic simulations, discovering that partial packing of the Im9 core is frustrated. This frustration is overcome when the entire core of the protein is assembled, consistent with experiment indicating no significant kinetic trapping in Im9 folding. Our systematic analysis thus reveals a subtle, contextual aspect of biomolecular recognition and provides a general approach to characterize folding frustration.
Collapse
Affiliation(s)
- Loan Huynh
- Department of Biochemistry, University of Toronto, Toronto, Ontario, Canada
| | - Chris Neale
- Department of Physics, Applied Physics, and Astronomy, Rensselaer Polytechnic Institute, Troy, New York, United States of America
| | - Régis Pomès
- Department of Biochemistry, University of Toronto, Toronto, Ontario, Canada
- Molecular Medicine, The Hospital for Sick Children, Toronto, Ontario, Canada
- * E-mail: (HSC); (RP)
| | - Hue Sun Chan
- Department of Biochemistry, University of Toronto, Toronto, Ontario, Canada
- Department of Molecular Genetics, University of Toronto, Toronto, Ontario, Canada
- * E-mail: (HSC); (RP)
| |
Collapse
|
22
|
Mazurenko S, Kunka A, Beerens K, Johnson CM, Damborsky J, Prokop Z. Exploration of Protein Unfolding by Modelling Calorimetry Data from Reheating. Sci Rep 2017; 7:16321. [PMID: 29176711 PMCID: PMC5701188 DOI: 10.1038/s41598-017-16360-y] [Citation(s) in RCA: 31] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/03/2017] [Accepted: 11/10/2017] [Indexed: 01/14/2023] Open
Abstract
Studies of protein unfolding mechanisms are critical for understanding protein functions inside cells, de novo protein design as well as defining the role of protein misfolding in neurodegenerative disorders. Calorimetry has proven indispensable in this regard for recording full energetic profiles of protein unfolding and permitting data fitting based on unfolding pathway models. While both kinetic and thermodynamic protein stability are analysed by varying scan rates and reheating, the latter is rarely used in curve-fitting, leading to a significant loss of information from experiments. To extract this information, we propose fitting both first and second scans simultaneously. Four most common single-peak transition models are considered: (i) fully reversible, (ii) fully irreversible, (iii) partially reversible transitions, and (iv) general three-state models. The method is validated using calorimetry data for chicken egg lysozyme, mutated Protein A, three wild-types of haloalkane dehalogenases, and a mutant stabilized by protein engineering. We show that modelling of reheating increases the precision of determination of unfolding mechanisms, free energies, temperatures, and heat capacity differences. Moreover, this modelling indicates whether alternative refolding pathways might occur upon cooling. The Matlab-based data fitting software tool and its user guide are provided as a supplement.
Collapse
Affiliation(s)
- Stanislav Mazurenko
- Loschmidt Laboratories, Department of Experimental Biology and Research Centre for Toxic Compounds in the Environment RECETOX, Faculty of Science, Masaryk University, Kamenice 5/A13, 625 00, Brno, Czech Republic
| | - Antonin Kunka
- Loschmidt Laboratories, Department of Experimental Biology and Research Centre for Toxic Compounds in the Environment RECETOX, Faculty of Science, Masaryk University, Kamenice 5/A13, 625 00, Brno, Czech Republic
- International Clinical Research Center, St. Anne's University Hospital, Pekarska 53, 656 91, Brno, Czech Republic
| | - Koen Beerens
- Loschmidt Laboratories, Department of Experimental Biology and Research Centre for Toxic Compounds in the Environment RECETOX, Faculty of Science, Masaryk University, Kamenice 5/A13, 625 00, Brno, Czech Republic
| | - Christopher M Johnson
- Biophysics Facilities, MRC Laboratory of Molecular Biology, Francis Crick Avenue, Cambridge Biomedical Campus, Cambridge, CB2 0QH, UK
| | - Jiri Damborsky
- Loschmidt Laboratories, Department of Experimental Biology and Research Centre for Toxic Compounds in the Environment RECETOX, Faculty of Science, Masaryk University, Kamenice 5/A13, 625 00, Brno, Czech Republic.
- International Clinical Research Center, St. Anne's University Hospital, Pekarska 53, 656 91, Brno, Czech Republic.
| | - Zbynek Prokop
- Loschmidt Laboratories, Department of Experimental Biology and Research Centre for Toxic Compounds in the Environment RECETOX, Faculty of Science, Masaryk University, Kamenice 5/A13, 625 00, Brno, Czech Republic.
- International Clinical Research Center, St. Anne's University Hospital, Pekarska 53, 656 91, Brno, Czech Republic.
| |
Collapse
|
23
|
Ferreiro DU, Komives EA, Wolynes PG. Frustration, function and folding. Curr Opin Struct Biol 2017; 48:68-73. [PMID: 29101782 DOI: 10.1016/j.sbi.2017.09.006] [Citation(s) in RCA: 60] [Impact Index Per Article: 8.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/23/2017] [Revised: 09/25/2017] [Accepted: 09/27/2017] [Indexed: 01/08/2023]
Abstract
Natural protein molecules are exceptional polymers. Encoded in apparently random strings of amino-acids, these objects perform clear physical tasks that are rare to find by simple chance. Accurate folding, specific binding, powerful catalysis, are examples of basic chemical activities that the great majority of polypeptides do not display, and are thought to be the outcome of the natural history of proteins. Function, a concept genuine to Biology, is at the core of evolution and often conflicts with the physical constraints. Locating the frustration between discrepant goals in a recurrent system leads to fundamental insights about the chances and necessities that shape the encoding of biological information.
Collapse
Affiliation(s)
- Diego U Ferreiro
- Protein Physiology Lab, FCEyN-Universidad de Buenos Aires, IQUIBICEN/CONICET, Intendente Güiraldes 2160, Ciudad Universitaria, C1428EGA Buenos Aires, Argentina
| | - Elizabeth A Komives
- Department of Chemistry and Biochemistry, University of California, San Diego, 9500 Gilman Drive, La Jolla, CA 92092-0378, USA
| | - Peter G Wolynes
- Center for Theoretical Biological Physics, Rice University, Houston, TX, USA; Department of Chemistry, Rice University, Houston, TX, USA; Department of Physics, Rice University, Houston, TX, USA; Department of Biosciences, Rice University, Houston, TX, USA
| |
Collapse
|
24
|
Goyal S, Chattopadhyay A, Kasavajhala K, Priyakumar UD. Role of Urea–Aromatic Stacking Interactions in Stabilizing the Aromatic Residues of the Protein in Urea-Induced Denatured State. J Am Chem Soc 2017; 139:14931-14946. [DOI: 10.1021/jacs.7b05463] [Citation(s) in RCA: 37] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/07/2023]
Affiliation(s)
- Siddharth Goyal
- Center for Computational
Natural Sciences and Bioinformatics, International Institute of Information Technology, Hyderabad 500 032, India
| | - Aditya Chattopadhyay
- Center for Computational
Natural Sciences and Bioinformatics, International Institute of Information Technology, Hyderabad 500 032, India
| | - Koushik Kasavajhala
- Center for Computational
Natural Sciences and Bioinformatics, International Institute of Information Technology, Hyderabad 500 032, India
| | - U. Deva Priyakumar
- Center for Computational
Natural Sciences and Bioinformatics, International Institute of Information Technology, Hyderabad 500 032, India
| |
Collapse
|
25
|
|
26
|
Narayan A, Campos LA, Bhatia S, Fushman D, Naganathan AN. Graded Structural Polymorphism in a Bacterial Thermosensor Protein. J Am Chem Soc 2017; 139:792-802. [DOI: 10.1021/jacs.6b10608] [Citation(s) in RCA: 22] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Affiliation(s)
- Abhishek Narayan
- Department of Biotechnology, Bhupat & Jyoti Mehta School of Biosciences, Indian Institute of Technology Madras (IITM), Chennai 600036, India
| | - Luis A. Campos
- National Biotechnology Center, Consejo Superior
de Investigaciones Científicas, Darwin 3, Campus de Cantoblanco, 28049 Madrid, Spain
| | - Sandhya Bhatia
- National Centre for Biological Sciences (NCBS), Tata Institute of Fundamental Research, Bangalore 560065, India
| | - David Fushman
- Department
of Chemistry and Biochemistry, Center for Biomolecular Structure and
Organization, University of Maryland, College Park, Maryland 20742, United States
| | - Athi N. Naganathan
- Department of Biotechnology, Bhupat & Jyoti Mehta School of Biosciences, Indian Institute of Technology Madras (IITM), Chennai 600036, India
| |
Collapse
|
27
|
Yadahalli S, Gosavi S. Packing energetics determine the folding routes of the RNase-H proteins. Phys Chem Chem Phys 2017; 19:9164-9173. [DOI: 10.1039/c6cp08940b] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Abstract
The RNase-H proteins show a diverse range of folding routes with structurally distinct folding nuclei.
Collapse
Affiliation(s)
- Shilpa Yadahalli
- National Centre for Biological Sciences
- Tata Institute of Fundamental Research
- Bangalore-560065
- India
- Manipal University
| | - Shachi Gosavi
- National Centre for Biological Sciences
- Tata Institute of Fundamental Research
- Bangalore-560065
- India
| |
Collapse
|
28
|
Are protein hubs faster folders? Exploration based on Escherichia coli proteome. Amino Acids 2016; 48:2747-2753. [DOI: 10.1007/s00726-016-2309-x] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/01/2016] [Accepted: 08/05/2016] [Indexed: 10/21/2022]
|
29
|
Malhotra P, Udgaonkar JB. How cooperative are protein folding and unfolding transitions? Protein Sci 2016; 25:1924-1941. [PMID: 27522064 PMCID: PMC5079258 DOI: 10.1002/pro.3015] [Citation(s) in RCA: 55] [Impact Index Per Article: 6.9] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/26/2016] [Revised: 08/09/2016] [Accepted: 08/09/2016] [Indexed: 11/12/2022]
Abstract
A thermodynamically and kinetically simple picture of protein folding envisages only two states, native (N) and unfolded (U), separated by a single activation free energy barrier, and interconverting by cooperative two-state transitions. The folding/unfolding transitions of many proteins occur, however, in multiple discrete steps associated with the formation of intermediates, which is indicative of reduced cooperativity. Furthermore, much advancement in experimental and computational approaches has demonstrated entirely non-cooperative (gradual) transitions via a continuum of states and a multitude of small energetic barriers between the N and U states of some proteins. These findings have been instrumental towards providing a structural rationale for cooperative versus noncooperative transitions, based on the coupling between interaction networks in proteins. The cooperativity inherent in a folding/unfolding reaction appears to be context dependent, and can be tuned via experimental conditions which change the stabilities of N and U. The evolution of cooperativity in protein folding transitions is linked closely to the evolution of function as well as the aggregation propensity of the protein. A large activation energy barrier in a fully cooperative transition can provide the kinetic control required to prevent the accumulation of partially unfolded forms, which may promote aggregation. Nevertheless, increasing evidence for barrier-less "downhill" folding, as well as for continuous "uphill" unfolding transitions, indicate that gradual non-cooperative processes may be ubiquitous features on the free energy landscape of protein folding.
Collapse
Affiliation(s)
- Pooja Malhotra
- National Centre for Biological Sciences, Tata Institute of Fundamental Research, Bengaluru, 560065, India
| | - Jayant B Udgaonkar
- National Centre for Biological Sciences, Tata Institute of Fundamental Research, Bengaluru, 560065, India.
| |
Collapse
|
30
|
Mascarenhas NM, Gosavi S. Protein Domain-Swapping Can Be a Consequence of Functional Residues. J Phys Chem B 2016; 120:6929-38. [PMID: 27331242 DOI: 10.1021/acs.jpcb.6b03968] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Abstract
Monomer topology has been implicated in domain-swapping, a potential first step on the route to disease-causing protein aggregation. Despite having the same topology (β1-α1-β2-β3-β4-β5), the cysteine protease inhibitor stefin-B domain swaps more readily than a single-chain variant of the heterodimeric sweet protein monellin (scMn). Here, we computationally study the folding of stefin-B and scMn in order to understand the molecular basis for the difference in their domain-swapping propensities. In agreement with experiments, our structure-based simulations show that scMn folds cooperatively without the population of an intermediate while stefin-B populates an equilibrium intermediate state. Since the simulation intermediate has only one domain structured (β3-β4-β5), it can directly lead to domain-swapping. Using computational variants of stefin-B, we show that the population of this intermediate is caused by regions of stefin-B that have been implicated in protease inhibition. We also find that the protease-binding regions are located on two structural elements and localized in space. In contrast, the residues that contribute to the sweetness of monellin are not localized to a few structural elements but are distributed over the protein fold. We conclude that the distributed functional residues of monellin do not induce large local perturbations in the protein structure, eliminating the formation of folding intermediates and in turn domain-swapping. On the other hand, the localized protease-binding regions of stefin-B promote the formation of a folding intermediate which can lead to domain-swapping. Thus, domain-swapping can be a direct consequence of the constraints that function imposes on the protein structure.
Collapse
Affiliation(s)
- Nahren Manuel Mascarenhas
- Simons Centre for the Study of Living Machines, National Centre for Biological Sciences, Tata Institute of Fundamental Research , Bangalore 560065, India
| | - Shachi Gosavi
- Simons Centre for the Study of Living Machines, National Centre for Biological Sciences, Tata Institute of Fundamental Research , Bangalore 560065, India
| |
Collapse
|
31
|
Using natural sequences and modularity to design common and novel protein topologies. Curr Opin Struct Biol 2016; 38:26-36. [PMID: 27270240 DOI: 10.1016/j.sbi.2016.05.007] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/01/2016] [Revised: 05/13/2016] [Accepted: 05/18/2016] [Indexed: 02/07/2023]
Abstract
Protein design is still a challenging undertaking, often requiring multiple attempts or iterations for success. Typically, the source of failure is unclear, and scoring metrics appear similar between successful and failed cases. Nevertheless, the use of sequence statistics, modularity and symmetry from natural proteins, combined with computational design both at the coarse-grained and atomistic levels is propelling a new wave of design efforts to success. Here we highlight recent examples of design, showing how the wealth of natural protein sequence and topology data may be leveraged to reduce the search space and increase the likelihood of achieving desired outcomes.
Collapse
|