1
|
Fersht AR. From covalent transition states in chemistry to noncovalent in biology: from β- to Φ-value analysis of protein folding. Q Rev Biophys 2024; 57:e4. [PMID: 38597675 DOI: 10.1017/s0033583523000045] [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] [Indexed: 04/11/2024]
Abstract
Solving the mechanism of a chemical reaction requires determining the structures of all the ground states on the pathway and the elusive transition states linking them. 2024 is the centenary of Brønsted's landmark paper that introduced the β-value and structure-activity studies as the only experimental means to infer the structures of transition states. It involves making systematic small changes in the covalent structure of the reactants and analysing changes in activation and equilibrium-free energies. Protein engineering was introduced for an analogous procedure, Φ-value analysis, to analyse the noncovalent interactions in proteins central to biological chemistry. The methodology was developed first by analysing noncovalent interactions in transition states in enzyme catalysis. The mature procedure was then applied to study transition states in the pathway of protein folding - 'part (b) of the protein folding problem'. This review describes the development of Φ-value analysis of transition states and compares and contrasts the interpretation of β- and Φ-values and their limitations. Φ-analysis afforded the first description of transition states in protein folding at the level of individual residues. It revealed the nucleation-condensation folding mechanism of protein domains with the transition state as an expanded, distorted native structure, containing little fully formed secondary structure but many weak tertiary interactions. A spectrum of transition states with various degrees of structural polarisation was then uncovered that spanned from nucleation-condensation to the framework mechanism of fully formed secondary structure. Φ-analysis revealed how movement of the expanded transition state on an energy landscape accommodates the transition from framework to nucleation-condensation mechanisms with a malleability of structure as a unifying feature of folding mechanisms. Such movement follows the rubric of analysis of classical covalent chemical mechanisms that began with Brønsted. Φ-values are used to benchmark computer simulation, and Φ and simulation combine to describe folding pathways at atomic resolution.
Collapse
Affiliation(s)
- Alan R Fersht
- MRC Laboratory of Molecular Biology, Cambridge, UK
- Yusuf Hamied Department of Chemistry, University of Cambridge, Cambridge, UK
- Gonville and Caius College, University of Cambridge, Cambridge, UK
| |
Collapse
|
2
|
Das D, Ainavarapu SRK. Circular permutation at azurin's active site slows down its folding. J Biol Inorg Chem 2023; 28:737-749. [PMID: 37957357 DOI: 10.1007/s00775-023-02023-z] [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: 05/08/2023] [Accepted: 09/26/2023] [Indexed: 11/15/2023]
Abstract
Circular permutation (CP) is a technique by which the primary sequence of a protein is rearranged to create new termini. The connectivity of the protein is altered but the overall protein structure generally remains unperturbed. Understanding the effect of CP can help design robust proteins for numerous applications such as in genetic engineering, optoelectronics, and improving catalytic activity. Studies on different protein topologies showed that CP usually affects protein stability as well as unfolding rates. Though a significant number of proteins contain metals or other cofactors, reports of metalloprotein CPs are rare. Thus, we chose a bacterial metalloprotein, azurin, and its CP within the metal-binding site (cpF114). We studied the stabilities, folding, and unfolding rates of apo- and Zn2+-bound CP azurin using fluorescence and circular dichroism. The introduced CP had destabilizing effects on the protein. Also, the folding of the Zn2+-CP protein was much slower than that of the Zn2+-WT or apo-protein. We compared this study to our previously reported azurin-cpN42, where we had observed an equilibrium and kinetic intermediate. cpF114 exhibits an apparent two-state equilibrium unfolding but has an off-pathway kinetic intermediate. Our study hinted at CP as a method to modify the energy landscape of proteins to alter their folding pathways. WT azurin, being a faster folder, may have evolved to optimize the folding rate of metal-bound protein compared to its CPs, albeit all of them have the same structure and function. Our study underscores that protein sequence and protein termini positions are crucial for metalloproteins. TOC Figure. (Top) Zn2+-azurin WT structure (PDB code: 1E67) and 2-D topology diagram of Zn2+-cpF114 azurin. (Bottom) Cartoon diagram representing folding (red arrows) and unfolding (blue arrows) of apo- and Zn2+- WT and cpF114 azurins. The width of the arrows represents the rate of the corresponding processes.
Collapse
Affiliation(s)
- Debanjana Das
- Department of Chemical Sciences, Tata Institute of Fundamental Research, Dr. Homi Bhabha Road, Colaba, Mumbai, 400005, India
| | - Sri Rama Koti Ainavarapu
- Department of Chemical Sciences, Tata Institute of Fundamental Research, Dr. Homi Bhabha Road, Colaba, Mumbai, 400005, India.
| |
Collapse
|
3
|
Das D, Yadav P, Mitra S, Ainavarapu SRK. Metal-binding and circular permutation-dependent thermodynamic and kinetic stability of azurin. Proteins 2023; 91:634-648. [PMID: 36511110 DOI: 10.1002/prot.26454] [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: 09/03/2022] [Revised: 11/18/2022] [Accepted: 11/25/2022] [Indexed: 12/15/2022]
Abstract
Native topology is known to determine the folding kinetics and the energy landscape of proteins. Furthermore, the circular permutation (CP) of proteins alters the order of the secondary structure connectivity while retaining the three-dimensional structure, making it an elegant and powerful approach to altering native topology. Previous studies elucidated the influence of CP in proteins with different folds such as Greek key β-barrel, β-sandwich, β-α-β, and all α-Greek key. CP mainly affects the protein stability and unfolding kinetics, while folding kinetics remains mostly unaltered. However, the effect of CP on metalloproteins is yet to be elaborately studied. The active site of metalloproteins poses an additional complexity in studying protein folding. Here, we investigate a CP variant (cpN42) of azurin-in both metal-free and metal-bound (holo) forms. As observed earlier in other proteins, apo-forms of wild-type (WT) and cpN42 fold with similar rates. In contrast, zinc-binding accelerates the folding of WT but decelerates the folding of cpN42. On zinc-binding, the spontaneous folding rate of WT increases by >250 times that of cpN42, which is unprecedented and the highest for any CP to date. On the other hand, zinc-binding reduces the spontaneous unfolding rate of cpN42 by ~100 times, making the WT and CP azurins unfold at similar rates. Our study demonstrates metal binding as a novel way to modulate the unfolding and folding rates of CPs compared to their WT counterparts. We hope our study increases the understanding of the effect of CP on the folding mechanism and energy landscape of metalloproteins.
Collapse
Affiliation(s)
- Debanjana Das
- Department of Chemical Sciences, Dr. Homi Bhabha Road, Colaba, Tata Institute of Fundamental Research, Mumbai, India
| | - Priya Yadav
- Department of Chemical Sciences, Dr. Homi Bhabha Road, Colaba, Tata Institute of Fundamental Research, Mumbai, India
| | - Soumyajit Mitra
- Department of Chemical Sciences, Dr. Homi Bhabha Road, Colaba, Tata Institute of Fundamental Research, Mumbai, India
| | - Sri Rama Koti Ainavarapu
- Department of Chemical Sciences, Dr. Homi Bhabha Road, Colaba, Tata Institute of Fundamental Research, Mumbai, India
| |
Collapse
|
4
|
Zegarra FC, Homouz D, Wittung-Stafshede P, Cheung MS. The Zero-Order Loop in Apoazurin Modulates Folding Mechanism In Silico. J Phys Chem B 2021; 125:3501-3509. [PMID: 33818090 DOI: 10.1021/acs.jpcb.1c00219] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/02/2023]
Abstract
Pseudomonas aeruginosa apoazurin (apo, without the copper cofactor) has a single disulfide bond between residues 3 and 26 and unfolds in a two-state reaction in vitro. The disulfide bond covalently connects the N-termini of β-strands 1 and 3; thereby, it creates a zero-order loop or a "cinch" that restricts conformational space. Covalent loops and threaded topologies are emerging as important structural elements in folded proteins and may be important for function. In order to understand the role of a zero-order loop in the folding process of a protein, here we used coarse-grained molecular dynamics (CGMD) simulations in silico to compare two variants of apoazurin: one named "loop" which contained the disulfide, and another named "open" in which the disulfide bond between residues 3 and 26 was removed. CGMD simulations were performed to probe the stability and unfolding pathway of the two apoazurin variants at different urea concentrations and temperatures. Our results show that the covalent loop plays a prominent role in the unfolding mechanism of apoazurin; its removal alters both the folding-transition state and the unfolded-state ensemble of conformations. We propose that modulation of azurin's folding landscape by the disulfide bridge may be related to both copper capturing and redox sensing.
Collapse
Affiliation(s)
- Fabio C Zegarra
- Laboratorio de Sistemas Inteligentes, EPIME, Universidad Nacional Tecnológica de Lima Sur (UNTELS), Lima 15834, Peru.,Department of Physics, University of Houston, Houston, Texas 77204, United States.,Center for Theoretical Biological Physics, Rice University, Houston, Texas 77005, United States
| | - Dirar Homouz
- Department of Physics, Khalifa University of Science and Technology, P.O. Box 127788, Abu Dhabi, 127788 United Arab Emirates.,Department of Physics, University of Houston, Houston, Texas 77204, United States.,Center for Theoretical Biological Physics, Rice University, Houston, Texas 77005, United States
| | - Pernilla Wittung-Stafshede
- Department of Biology and Biological Engineering, Chalmers University of Technology, 412 96 Gothenburg, Sweden
| | - Margaret S Cheung
- Pacific Northwest National Laboratory, Research Science Center, Seattle, Washington 98109, United States.,Department of Physics, University of Washington, Seattle, Washington 98195, United States.,Department of Physics, University of Houston, Houston, Texas 77204, United States.,Center for Theoretical Biological Physics, Rice University, Houston, Texas 77005, United States
| |
Collapse
|
5
|
Abstract
Copper is a redox-active transition metal ion required for the function of many essential human proteins. For biosynthesis of proteins coordinating copper, the metal may bind before, during or after folding of the polypeptide. If the metal binds to unfolded or partially folded structures of the protein, such coordination may modulate the folding reaction. The molecular understanding of how copper is incorporated into proteins requires descriptions of chemical, thermodynamic, kinetic and structural parameters involved in the formation of protein-metal complexes. Because free copper ions are toxic, living systems have elaborate copper-transport systems that include particular proteins that facilitate efficient and specific delivery of copper ions to target proteins. Therefore, these pathways become an integral part of copper protein folding in vivo. This review summarizes biophysical-molecular in vitro work assessing the role of copper in folding and stability of copper-binding proteins as well as protein-protein copper exchange reactions between human copper transport proteins. We also describe some recent findings about the participation of copper ions and copper proteins in protein misfolding and aggregation reactions in vitro.
Collapse
|
6
|
Selection originating from protein stability/foldability: Relationships between protein folding free energy, sequence ensemble, and fitness. J Theor Biol 2017; 433:21-38. [DOI: 10.1016/j.jtbi.2017.08.018] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/27/2017] [Revised: 07/27/2017] [Accepted: 08/21/2017] [Indexed: 11/19/2022]
|
7
|
Wako H, Abe H. Characterization of protein folding by a Φ-value calculation with a statistical-mechanical model. Biophys Physicobiol 2016; 13:263-279. [PMID: 28409079 PMCID: PMC5221509 DOI: 10.2142/biophysico.13.0_263] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/01/2016] [Accepted: 09/20/2016] [Indexed: 12/01/2022] Open
Abstract
The Φ-value analysis approach provides information about transition-state structures along the folding pathway of a protein by measuring the effects of an amino acid mutation on folding kinetics. Here we compared the theoretically calculated Φ values of 27 proteins with their experimentally observed Φ values; the theoretical values were calculated using a simple statistical-mechanical model of protein folding. The theoretically calculated Φ values reflected the corresponding experimentally observed Φ values with reasonable accuracy for many of the proteins, but not for all. The correlation between the theoretically calculated and experimentally observed Φ values strongly depends on whether the protein-folding mechanism assumed in the model holds true in real proteins. In other words, the correlation coefficient can be expected to illuminate the folding mechanisms of proteins, providing the answer to the question of which model more accurately describes protein folding: the framework model or the nucleation-condensation model. In addition, we tried to characterize protein folding with respect to various properties of each protein apart from the size and fold class, such as the free-energy profile, contact-order profile, and sensitivity to the parameters used in the Φ-value calculation. The results showed that any one of these properties alone was not enough to explain protein folding, although each one played a significant role in it. We have confirmed the importance of characterizing protein folding from various perspectives. Our findings have also highlighted that protein folding is highly variable and unique across different proteins, and this should be considered while pursuing a unified theory of protein folding.
Collapse
Affiliation(s)
- Hiroshi Wako
- School of Social Sciences, Waseda University, Shinjuku, Tokyo 169-8050, Japan
| | - Haruo Abe
- Department of Electrical Engineering, Nishinippon Institute of Technology, Miyako, Fukuoka 800-0394, Japan
| |
Collapse
|
8
|
Affiliation(s)
- Pernilla Wittung-Stafshede
- Department of Biology and Biological Engineering; Chalmers University of Technology; 41296 Gothenburg Sweden
| |
Collapse
|
9
|
Contini A, Tiana G. A many-body term improves the accuracy of effective potentials based on protein coevolutionary data. J Chem Phys 2016; 143:025103. [PMID: 26178131 DOI: 10.1063/1.4926665] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
The study of correlated mutations in alignments of homologous proteins proved to be successful not only in the prediction of their native conformation but also in the development of a two-body effective potential between pairs of amino acids. In the present work, we extend the effective potential, introducing a many-body term based on the same theoretical framework, making use of a principle of maximum entropy. The extended potential performs better than the two-body one in predicting the energetic effect of 308 mutations in 14 proteins (including membrane proteins). The average value of the parameters of the many-body term correlates with the degree of hydrophobicity of the corresponding residues, suggesting that this term partly reflects the effect of the solvent.
Collapse
Affiliation(s)
- A Contini
- Department of Physics, Università degli Studi di Milano, via Celoria 16, 20133 Milano, Italy
| | - G Tiana
- Department of Physics, Università degli Studi di Milano, and INFN, via Celoria 16, 20133 Milano, Italy
| |
Collapse
|
10
|
Lui S, Tiana G. The network of stabilizing contacts in proteins studied by coevolutionary data. J Chem Phys 2014; 139:155103. [PMID: 24160546 DOI: 10.1063/1.4826096] [Citation(s) in RCA: 28] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The primary structure of proteins, that is their sequence, represents one of the most abundant sets of experimental data concerning biomolecules. The study of correlations in families of co-evolving proteins by means of an inverse Ising-model approach allows to obtain information on their native conformation. Following up on a recent development along this line, we optimize the algorithm to calculate effective energies between the residues, validating the approach both back-calculating interaction energies in a model system, and predicting the free energies associated to mutations in real systems. Making use of these effective energies, we study the network of interactions which stabilizes the native conformation of some well-studied proteins, showing that it displays different properties than the associated contact network.
Collapse
Affiliation(s)
- Sara Lui
- Department of Physics, Università degli Studi di Milano, via Celoria 16, 20133 Milano, Italy
| | | |
Collapse
|
11
|
Christiansen A, Wittung-Stafshede P. Quantification of excluded volume effects on the folding landscape of Pseudomonas aeruginosa apoazurin in vitro. Biophys J 2014; 105:1689-99. [PMID: 24094410 DOI: 10.1016/j.bpj.2013.08.038] [Citation(s) in RCA: 51] [Impact Index Per Article: 5.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/20/2013] [Revised: 08/09/2013] [Accepted: 08/15/2013] [Indexed: 11/15/2022] Open
Abstract
Proteins fold and function inside cells that are crowded with macromolecules. Here, we address the role of the resulting excluded volume effects by in vitro spectroscopic studies of Pseudomonas aeruginosa apoazurin stability (thermal and chemical perturbations) and folding kinetics (chemical perturbation) as a function of increasing levels of crowding agents dextran (sizes 20, 40, and 70 kDa) and Ficoll 70. We find that excluded volume theory derived by Minton quantitatively captures the experimental effects when crowding agents are modeled as arrays of rods. This finding demonstrates that synthetic crowding agents are useful for studies of excluded volume effects. Moreover, thermal and chemical perturbations result in free energy effects by the presence of crowding agents that are identical, which shows that the unfolded state is energetically the same regardless of method of unfolding. This also underscores the two-state approximation for apoazurin's unfolding reaction and suggests that thermal and chemical unfolding experiments can be used in an interchangeable way. Finally, we observe increased folding speed and invariant unfolding speed for apoazurin in the presence of macromolecular crowding agents, a result that points to unfolded-state perturbations. Although the absolute magnitude of excluded volume effects on apoazurin is only on the order of 1-3 kJ/mol, differences of this scale may be biologically significant.
Collapse
|
12
|
Tobin PH, Wilson CJ. Examining Photoinduced Energy Transfer in Pseudomonas aeruginosa Azurin. J Am Chem Soc 2014; 136:1793-802. [DOI: 10.1021/ja412308r] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/21/2023]
Affiliation(s)
- Peter H. Tobin
- Department
of Chemical and Environmental Engineering, ‡Department of Biomedical
Engineering, and §Department of Molecular Biophysics and Biochemistry, Yale University, New Haven, Connecticut 06520-8286, United States
| | - Corey J. Wilson
- Department
of Chemical and Environmental Engineering, ‡Department of Biomedical
Engineering, and §Department of Molecular Biophysics and Biochemistry, Yale University, New Haven, Connecticut 06520-8286, United States
| |
Collapse
|
13
|
Warren JJ, Gray HB, Winkler JR, Kozak JJ. Euclidean perspective on the unfolding of azurin: angular correlations. Mol Phys 2013. [DOI: 10.1080/00268976.2013.787153] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/27/2022]
|
14
|
Role of metal in folding and stability of copper proteins in vitro. BIOCHIMICA ET BIOPHYSICA ACTA-MOLECULAR CELL RESEARCH 2012; 1823:1594-603. [PMID: 22306006 DOI: 10.1016/j.bbamcr.2012.01.013] [Citation(s) in RCA: 61] [Impact Index Per Article: 5.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/26/2011] [Revised: 01/09/2012] [Accepted: 01/18/2012] [Indexed: 01/04/2023]
Abstract
Metal coordination is required for function of many proteins. For biosynthesis of proteins coordinating a metal, the question arises if the metal binds before, during or after folding of the polypeptide. Moreover, when the metal is bound to the protein, how does its coordination affect biophysical properties such as stability and dynamics? Understanding how metals are utilized by proteins in cells on a molecular level requires accurate descriptions of the thermodynamic and kinetic parameters involved in protein-metal complexes. Copper is one of the essential transition metals found in the active sites of many key proteins. To avoid toxicity of free copper ions, living systems have developed elaborate copper-transport systems that involve dedicated proteins that facilitate efficient and specific delivery of copper to target proteins. This review describes in vitro and in silico biophysical work assessing the role of copper in folding and stability of copper-binding proteins. Examples of proteins discussed are: a blue-copper protein (Pseudomonas aeruginosa azurin), members of copper-transport systems (bacterial CopZ, human Atox1 and ATP7B domains) and multi-copper ferroxidases (yeast Fet3p and human ceruloplasmin). The consequences of interactions between copper proteins and platinum-complexes are also discussed. This article is part of a Special Issue entitled: Cell Biology of Metals.
Collapse
|
15
|
Investigation of an anomalously accelerating substitution in the folding of a prototypical two-state protein. J Mol Biol 2010; 403:446-58. [PMID: 20816985 DOI: 10.1016/j.jmb.2010.08.049] [Citation(s) in RCA: 13] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/30/2010] [Revised: 08/26/2010] [Accepted: 08/27/2010] [Indexed: 11/22/2022]
Abstract
The folding rates of two-state single-domain proteins are generally resistant to small-scale changes in amino acid sequence. For example, having surveyed here over 700 single-residue substitutions in 24 well-characterized two-state proteins, we find that the majority (55%) of these substitutions affect folding rates by less than a factor of 2, and that only 9% affect folding rates by more than a factor of 8. Among those substitutions that significantly affect folding rates, we find that accelerating substitutions are an order of magnitude less common than those that decelerate the process. One of the most extreme outliers in this data set, an arginine-to-phenylalanine substitution at position 48 (R48F) of chymotrypsin inhibitor 2 (CI2), accelerates the protein's folding rate by a factor of 36 relative to that of the wild-type protein and is the most accelerating substitution reported to date in a two-state protein. In order to better understand the origins of this anomalous behavior, we have characterized the kinetics of multiple additional substitutions at this position. We find that substitutions at position 48 in CI2 fall into two distinct classes. The first, comprising residues that ablate the charge of the wild-type arginine but retain the hydrophobicity of its alkane chain, accelerate folding by at least 10-fold. The second class, comprising all other residues, produces folding rates within a factor of two of the wild-type rate. A significant positive correlation between hydrophobicity and folding rate across all of the residues we have characterized at this position suggests that the hydrophobic methylene units of the wild-type arginine play a significant role in stabilizing the folding transition state. Likewise, studies of the pH dependence of the histidine substitution indicate a strong correlation between folding rate and charge state. Thus, mutations that ablate the arginine's positive charge while retaining the hydrophobic contacts of its methylene units tend to dramatically accelerate folding. Previous studies have suggested that arginine 48 plays an important functional role in CI2, which may explain why it is highly conserved despite the anomalously large deceleration it produces in the folding of this protein.
Collapse
|
16
|
Christiansen A, Wang Q, Samiotakis A, Cheung MS, Wittung-Stafshede P. Factors Defining Effects of Macromolecular Crowding on Protein Stability: An in Vitro/in Silico Case Study Using Cytochrome c. Biochemistry 2010; 49:6519-30. [DOI: 10.1021/bi100578x] [Citation(s) in RCA: 122] [Impact Index Per Article: 8.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Alexander Christiansen
- Department of Chemistry, Chemical Biological Center, Umeå University, 901 87 Umeå, Sweden
| | - Qian Wang
- Department of Physics, University of Houston, Houston, Texas 77204
| | | | | | | |
Collapse
|
17
|
Hills RD, Kathuria SV, Wallace LA, Day IJ, Brooks CL, Matthews CR. Topological frustration in beta alpha-repeat proteins: sequence diversity modulates the conserved folding mechanisms of alpha/beta/alpha sandwich proteins. J Mol Biol 2010; 398:332-50. [PMID: 20226790 DOI: 10.1016/j.jmb.2010.03.001] [Citation(s) in RCA: 27] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/09/2009] [Revised: 02/27/2010] [Accepted: 03/03/2010] [Indexed: 10/19/2022]
Abstract
The thermodynamic hypothesis of Anfinsen postulates that structures and stabilities of globular proteins are determined by their amino acid sequences. Chain topology, however, is known to influence the folding reaction, in that motifs with a preponderance of local interactions typically fold more rapidly than those with a larger fraction of nonlocal interactions. Together, the topology and sequence can modulate the energy landscape and influence the rate at which the protein folds to the native conformation. To explore the relationship of sequence and topology in the folding of beta alpha-repeat proteins, which are dominated by local interactions, we performed a combined experimental and simulation analysis on two members of the flavodoxin-like, alpha/beta/alpha sandwich fold. Spo0F and the N-terminal receiver domain of NtrC (NT-NtrC) have similar topologies but low sequence identity, enabling a test of the effects of sequence on folding. Experimental results demonstrated that both response-regulator proteins fold via parallel channels through highly structured submillisecond intermediates before accessing their cis prolyl peptide bond-containing native conformations. Global analysis of the experimental results preferentially places these intermediates off the productive folding pathway. Sequence-sensitive Gō-model simulations conclude that frustration in the folding in Spo0F, corresponding to the appearance of the off-pathway intermediate, reflects competition for intra-subdomain van der Waals contacts between its N- and C-terminal subdomains. The extent of transient, premature structure appears to correlate with the number of isoleucine, leucine, and valine (ILV) side chains that form a large sequence-local cluster involving the central beta-sheet and helices alpha2, alpha 3, and alpha 4. The failure to detect the off-pathway species in the simulations of NT-NtrC may reflect the reduced number of ILV side chains in its corresponding hydrophobic cluster. The location of the hydrophobic clusters in the structure may also be related to the differing functional properties of these response regulators. Comparison with the results of previous experimental and simulation analyses on the homologous CheY argues that prematurely folded unproductive intermediates are a common property of the beta alpha-repeat motif.
Collapse
Affiliation(s)
- Ronald D Hills
- Department of Molecular Biology and Kellogg School of Science and Technology, The Scripps Research Institute, 10550 North Torrey Pines Road TPC6, La Jolla, CA 92037, USA
| | | | | | | | | | | |
Collapse
|
18
|
Different Folding Pathways Taken by Highly Homologous Proteins, Goat α-Lactalbumin and Canine Milk Lysozyme. J Mol Biol 2010; 396:1361-78. [DOI: 10.1016/j.jmb.2010.01.021] [Citation(s) in RCA: 23] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/23/2009] [Revised: 01/10/2010] [Accepted: 01/11/2010] [Indexed: 11/19/2022]
|
19
|
Stagg L, Samiotakis A, Homouz D, Cheung MS, Wittung-Stafshede P. Residue-specific analysis of frustration in the folding landscape of repeat beta/alpha protein apoflavodoxin. J Mol Biol 2009; 396:75-89. [PMID: 19913555 DOI: 10.1016/j.jmb.2009.11.008] [Citation(s) in RCA: 12] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/24/2009] [Revised: 11/04/2009] [Accepted: 11/05/2009] [Indexed: 11/17/2022]
Abstract
Flavodoxin adopts the common repeat beta/alpha topology and folds in a complex kinetic reaction with intermediates. To better understand this reaction, we analyzed a set of Desulfovibrio desulfuricans apoflavodoxin variants with point mutations in most secondary structure elements by in vitro and in silico methods. By equilibrium unfolding experiments, we first revealed how different secondary structure elements contribute to overall protein resistance to heat and urea. Next, using stopped-flow mixing coupled with far-UV circular dichroism, we probed how individual residues affect the amount of structure formed in the experimentally detected burst-phase intermediate. Together with in silico folding route analysis of the same point-mutated variants and computation of growth in nucleation size during early folding, computer simulations suggested the presence of two competing folding nuclei at opposite sides of the central beta-strand 3 (i.e., at beta-strands 1 and 4), which cause early topological frustration (i.e., misfolding) in the folding landscape. Particularly, the extent of heterogeneity in folding nuclei growth correlates with the in vitro burst-phase circular dichroism amplitude. In addition, phi-value analysis (in vitro and in silico) of the overall folding barrier to apoflavodoxin's native state revealed that native-like interactions in most of the beta-strands must form in transition state. Our study reveals that an imbalanced competition between the two sides of apoflavodoxin's central beta-sheet directs initial misfolding, while proper alignment on both sides of beta-strand 3 is necessary for productive folding.
Collapse
Affiliation(s)
- Loren Stagg
- Department of Biochemistry and Cell Biology, Rice University, Houston, TX 77251, USA
| | | | | | | | | |
Collapse
|
20
|
Li H, Zhao Y, Zhen X, Yang B. Fluorescence spectra study the perturbations of CopC native fold by 2-p-toluidinynaphthalene-6-sulfonate. SPECTROCHIMICA ACTA. PART A, MOLECULAR AND BIOMOLECULAR SPECTROSCOPY 2009; 72:56-60. [PMID: 19019726 DOI: 10.1016/j.saa.2008.07.025] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/11/2008] [Revised: 07/04/2008] [Accepted: 07/17/2008] [Indexed: 05/27/2023]
Abstract
2-p-Toluidinynaphthalene-6-sulfonate (TNS) was discovered to perturb native fold of CopC protein and to induce loss of biological activity to some extent which was dependent on TNS concentration. Hydrophobic and electrostatic interactions were revealed to account for the perturbation by comparison with some analogy. TNS, with far low concentration of 10(-5) to 10(-4)M, is presented as a denaturant. So TNS should be deliberated in detecting macromolecular conformation change as single evidence at higher concentration.
Collapse
Affiliation(s)
- Huiqing Li
- Institute of Molecular Science, Key Laboratory of Chemical Biology and Molecular Engineering of Ministry of Education, Shanxi University, Taiyuan, Shanxi 030006, PR China
| | | | | | | |
Collapse
|
21
|
In vitro unfolding of yeast multicopper oxidase Fet3p variants reveals unique role of each metal site. Proc Natl Acad Sci U S A 2008; 105:19258-63. [PMID: 19033465 DOI: 10.1073/pnas.0806431105] [Citation(s) in RCA: 30] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
Fet3p from Saccharomyces cerevisiae is a multicopper oxidase (MCO) that contains 3 cupredoxin-like beta-barrel domains and 4 copper ions located in 3 distinct metal sites (T1 in domain 3, T2, and the binuclear T3 at the interface between domains 1 and 3). To better understand how protein structure and stability is defined by cofactor coordination in MCO proteins, we assessed thermal unfolding of apo and metallated forms of Fet3p by using spectroscopic and calorimetric methods in vitro (pH 7). We find that unfolding reactions of apo and different holo forms of Fet3p are irreversible reactions that depend on the scan rate. The domains in apo-Fet3p unfold sequentially [thermal midpoint (T(m)) of 45 degrees C, 62 degrees C, and 72 degrees C; 1 K/min]. Addition of T3 imposes strain in the apo structure that results in coupled domain unfolding and low stability (T(m) of 50 degrees C; 1 K/min). Further inclusion of T2 (i.e., only T1 absent) increases overall stability by approximately 5 degrees C but unfolding remains coupled in 1 step. Introduction of T1, producing fully-loaded holo-Fet3p (or in the absence of T2), results in stabilization of domain 3, which uncouples unfolding of the domains; unfolding of domain 2 occurs first along with Cu-site perturbations (T(m) 50-55 degrees C; 1 K/min), followed by unfolding of domains 1 and 3 ( approximately 65-70 degrees C; 1 K/min). Our results suggest that there is a metal-induced tradeoff between overall protein stability and metal coordination in members of the MCO family.
Collapse
|
22
|
Weikl TR. Loop-closure principles in protein folding. Arch Biochem Biophys 2008; 469:67-75. [PMID: 17662688 DOI: 10.1016/j.abb.2007.06.018] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/28/2007] [Revised: 06/20/2007] [Accepted: 06/22/2007] [Indexed: 10/23/2022]
Abstract
Simple theoretical concepts and models have been helpful to understand the folding rates and routes of single-domain proteins. As reviewed in this article, a physical principle that appears to underly these models is loop closure.
Collapse
Affiliation(s)
- Thomas R Weikl
- Max Planck Institute of Colloids and Interfaces, Department of Theory and Bio-Systems, 14424 Potsdam, Germany.
| |
Collapse
|
23
|
Weikl TR. Transition states in protein folding kinetics: modeling phi-values of small beta-sheet proteins. Biophys J 2007; 94:929-37. [PMID: 17905840 PMCID: PMC2186242 DOI: 10.1529/biophysj.107.109868] [Citation(s) in RCA: 14] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
Small single-domain proteins often exhibit only a single free-energy barrier, or transition state, between the denatured and the native state. The folding kinetics of these proteins is usually explored via mutational analysis. A central question is which structural information on the transition state can be derived from the mutational data. In this article, we model and structurally interpret mutational Phi-values for two small beta-sheet proteins, the PIN and the FBP WW domains. The native structure of these WW domains comprises two beta-hairpins that form a three-stranded beta-sheet. In our model, we assume that the transition state consists of two conformations in which either one of the hairpins is formed. Such a transition state has been recently observed in molecular dynamics folding-unfolding simulations of a small designed three-stranded beta-sheet protein. We obtain good agreement with the experimental data 1), by splitting up the mutation-induced free-energy changes into terms for the two hairpins and for the small hydrophobic core of the proteins; and 2), by fitting a single parameter, the relative degree to which hairpins 1 and 2 are formed in the transition state. The model helps us to understand how mutations affect the folding kinetics of WW domains, and captures also negative Phi-values that have been difficult to interpret.
Collapse
Affiliation(s)
- Thomas R Weikl
- Max Planck Institute of Colloids and Interfaces, Department of Theory and Bio-Systems, Potsdam, Germany.
| |
Collapse
|
24
|
Minimizing frustration by folding in an aqueous environment. Arch Biochem Biophys 2007; 469:118-31. [PMID: 17719000 DOI: 10.1016/j.abb.2007.07.007] [Citation(s) in RCA: 11] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/19/2007] [Accepted: 07/08/2007] [Indexed: 10/23/2022]
Abstract
Although life as we know it evolved in an aqueous medium, the properties of water are not completely understood. In this review, we focus on the role of water in guiding protein folding and stability. Specifically, we discuss the mechanisms of protein folding in an aqueous environment, the effects of water on the folding energy landscape as well as the transition state ensemble, and interactions of water with the folded state. We show that water cannot be viewed as a passive solvent, but rather, plays a very active role in the life of a protein.
Collapse
|
25
|
Zong C, Wilson CJ, Shen T, Wittung-Stafshede P, Mayo SL, Wolynes PG. Establishing the entatic state in folding metallated Pseudomonas aeruginosa azurin. Proc Natl Acad Sci U S A 2007; 104:3159-64. [PMID: 17301232 PMCID: PMC1805512 DOI: 10.1073/pnas.0611149104] [Citation(s) in RCA: 23] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
Understanding how the folding of proteins establishes their functional characteristics at the molecular level challenges both theorists and experimentalists. The simplest test beds for confronting this issue are provided by electron transfer proteins. The environment provided by the folded protein to the cofactor tunes the metal's electron transport capabilities as envisioned in the entatic hypothesis. To see how the entatic state is achieved one must study how the folding landscape affects and in turn is affected by the metal. Here, we develop a coarse-grained functional to explicitly model how the coordination of the metal (which results in a so-called entatic or rack-induced state) modifies the folding of the metallated Pseudomonas aeruginosa azurin. Our free-energy functional-based approach directly yields the proper nonlinear extra-thermodynamic free energy relationships for the kinetics of folding the wild type and several point-mutated variants of the metallated protein. The results agree quite well with corresponding laboratory experiments. Moreover, our modified free-energy functional provides a sufficient level of detail to explicitly model how the geometric entatic state of the metal modifies the dynamic folding nucleus of azurin.
Collapse
Affiliation(s)
- Chenghang Zong
- *Department of Chemistry and Biochemistry
- Center for Theoretical Biological Physics, University of California at San Diego, La Jolla, CA 92093
- To whom correspondence may be addressed. E-mail: , , or
| | - Corey J. Wilson
- Division of Biology, California Institute of Technology, Mail Code 147-75, Pasadena, CA 91125; and
- To whom correspondence may be addressed. E-mail: , , or
| | - Tongye Shen
- *Department of Chemistry and Biochemistry
- Center for Theoretical Biological Physics, University of California at San Diego, La Jolla, CA 92093
| | - Pernilla Wittung-Stafshede
- Departments of Biochemistry and Cell Biology and Chemistry, Rice University, 6100 Main Street, Houston, TX 77251
| | - Steven L. Mayo
- Division of Biology, California Institute of Technology, Mail Code 147-75, Pasadena, CA 91125; and
| | - Peter G. Wolynes
- *Department of Chemistry and Biochemistry
- Center for Theoretical Biological Physics, University of California at San Diego, La Jolla, CA 92093
- To whom correspondence may be addressed. E-mail: , , or
| |
Collapse
|
26
|
Lindberg MO, Oliveberg M. Malleability of protein folding pathways: a simple reason for complex behaviour. Curr Opin Struct Biol 2007; 17:21-9. [PMID: 17251003 DOI: 10.1016/j.sbi.2007.01.008] [Citation(s) in RCA: 92] [Impact Index Per Article: 5.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/01/2006] [Revised: 12/13/2006] [Accepted: 01/12/2007] [Indexed: 10/23/2022]
Abstract
Although the structures of native proteins are generally unique, the pathways by which they form are often free to vary. Some proteins fold by a multitude of different pathways, whereas others seem restricted to only one choice. An explanation for this variation in folding behaviour has recently emerged from studies of transition state changes: the number of accessible pathways is linked to the number of nucleation motifs contained within the native topology. We refer to these nucleation motifs as 'foldons', as they approach the size of an independent cooperative unit. Thus, with respect to pathway malleability and the composition of the folding funnel, proteins can be seen as modular assemblies of competing foldons. For the split beta-alpha-beta fold, these foldons are two-strand-helix motifs coupled by spatial overlap.
Collapse
Affiliation(s)
- Magnus O Lindberg
- Department of Biochemistry and Biophysics, Arrhenius Laboratories of Natural Sciences, Stockholm University, S-106 91 Stockholm, Sweden
| | | |
Collapse
|
27
|
Weikl TR, Dill KA. Transition-states in protein folding kinetics: the structural interpretation of Phi values. J Mol Biol 2006; 365:1578-86. [PMID: 17141267 DOI: 10.1016/j.jmb.2006.10.082] [Citation(s) in RCA: 25] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/28/2006] [Revised: 10/07/2006] [Accepted: 10/25/2006] [Indexed: 10/24/2022]
Abstract
Phi values are experimental measures of the effects of mutations on the folding kinetics of a protein. A central question is what structural information Phi values give about the transition-state of folding. Traditionally, a Phi value is interpreted as representing the "nativeness" of a mutated residue in the transition-state. However, this interpretation is often problematic. We present here a better structural interpretation of Phi values for mutations within a given helix. Our interpretation is based on a simple physical model that distinguishes between secondary and tertiary free energy contributions of helical residues. From a linear fit of the model to experimental data, we obtain two structural parameters: the extent of helix formation in the transition-state, and the nativeness of tertiary interactions in the transition-state. We apply the model to all proteins with well-characterized helices for which more than 10 Phi values are available: protein A, CI2, and protein L. The model is simple to apply to experimental data, captures nonclassical Phi values <0 or >1 in these helices, and explains how different mutations at a given site can lead to different Phi values.
Collapse
Affiliation(s)
- Thomas R Weikl
- Max Planck Institute of Colloids and Interfaces, Theory Department, 14424 Potsdam, Germany.
| | | |
Collapse
|
28
|
Chen M, Wilson CJ, Wu Y, Wittung-Stafshede P, Ma J. Correlation between Protein Stability Cores and Protein Folding Kinetics: A Case Study on Pseudomonas aeruginosa Apo-Azurin. Structure 2006; 14:1401-10. [PMID: 16962971 DOI: 10.1016/j.str.2006.07.007] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/08/2006] [Revised: 06/30/2006] [Accepted: 07/13/2006] [Indexed: 10/24/2022]
Abstract
This paper reports a combined computational and experimental study of the correlation between protein stability cores and folding kinetics. An empirical potential function was developed, and it was used for analyzing interaction energies among secondary structure elements. Studies on a beta sandwich protein, Pseudomonas aeruginosa azurin, showed that the computationally identified substructure with the strongest interactions in the native state is identical to the "interlocked pair" of beta strands, an invariant motif found in most sandwich-like proteins. Moreover, previous and new in vitro folding results revealed that the identified substructure harbors most residues that form native-like interactions in the folding transition state. These observations demonstrate that the potential function is effective in revealing the relative strength of interactions among various protein parts; they also strengthen the suggestion that the most stable regions in native proteins favor stable interactions early during folding.
Collapse
Affiliation(s)
- Mingzhi Chen
- Graduate Program of Structural and Computational Biology and Molecular Biophysics, Baylor College of Medicine, Houston, Texas 77030, USA
| | | | | | | | | |
Collapse
|
29
|
Wilson CJ, Apiyo D, Wittung-Stafshede P. Solvation of the folding-transition state in Pseudomonas aeruginosa azurin is modulated by metal: Solvation of azurin's folding nucleus. Protein Sci 2006; 15:843-52. [PMID: 16522792 PMCID: PMC2242485 DOI: 10.1110/ps.051838206] [Citation(s) in RCA: 11] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/24/2022]
Abstract
The role of water in protein folding, specifically its presence or not in the transition-state structure, is an unsolved question. There are two common classes of folding-transition states: diffuse transition states, in which almost all side chains have similar, rather low phi (phi) values, and polarized transition states, which instead display distinct substructures with very high phi-values. Apo-and zinc-forms of Pseudomonas aeruginosa azurin both fold in two-state equilibrium and kinetic reactions; while the apo-form exhibits a polarized transition state, the zinc form entails a diffuse, moving transition state. To examine the presence of water in these two types of folding-transition states, we probed the equilibrium and kinetic consequences of replacing core valines with isosteric threonines at six positions in azurin. In contrast to regular hydrophobic-to-alanine phi-value analysis, valine-to-threonine mutations do not disrupt the core packing but stabilize the unfolded state and can be used to assess the degree of solvation in the folding-transition state upon combination with regular phi-values. We find that the transition state for folding of apo-azurin appears completely dry, while that for zinc-azurin involves partially formed interactions that engage water molecules. This distinct difference between the apo-and holo-folding nuclei can be rationalized in terms of the shape of the free-energy barrier.
Collapse
Affiliation(s)
- Corey J Wilson
- Department of Biochemistry and Cell Biology, Rice University, Houston, Texas 77251, USA
| | | | | |
Collapse
|