1
|
Chang L, Perez A. Deciphering the Folding Mechanism of Proteins G and L and Their Mutants. J Am Chem Soc 2022; 144:14668-14677. [PMID: 35930769 DOI: 10.1021/jacs.2c04488] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/24/2023]
Abstract
Much of our understanding of folding mechanisms comes from interpretations of experimental ϕ and ψ value analysis, relating the differences in stability of the transition state ensemble (TSE) and folded state. We introduce a unified approach combining simulations and Bayesian inference to provide atomistic detail for the folding mechanism of proteins G and L and their mutants. Proteins G and L fold to similar topologies despite low sequence similarity, but differ in their folding pathways. A fast folding redesign of protein G, NuG2, switches folding pathways and folds through a similar pathway with protein L. A redesign of protein L also leads to faster folding, respecting the original folding pathway. Our Bayesian inference approach starts from the same prior on all systems and correctly identifies the folding mechanism for each of the four proteins, a success of the force field and sampling strategy. The approach is computationally efficient and correctly identifies the TSE and intermediate structures along the folding pathway in good agreement with experiments. We complement our findings by using two orthogonal approaches that differ in computational cost and interpretability. Adaptive sampling MD combined with the Markov state model provides a kinetic model that confirms the more complex folding mechanism of protein G and its mutant. Finally, a novel fragment decomposition approach using AlphaFold identifies preferences for secondary structure element combinations that follow the order of events observed in the folding pathways.
Collapse
Affiliation(s)
- Liwei Chang
- Department of Chemistry, University of Florida, Gainesville, Florida 32611, United States.,Quantum Theory Project, University of Florida, Gainesville, Florida 32611, United States
| | - Alberto Perez
- Department of Chemistry, University of Florida, Gainesville, Florida 32611, United States.,Quantum Theory Project, University of Florida, Gainesville, Florida 32611, United States
| |
Collapse
|
2
|
Baxa MC, Sosnick TR. Engineered Metal-Binding Sites to Probe Protein Folding Transition States: Psi Analysis. Methods Mol Biol 2022; 2376:31-63. [PMID: 34845602 DOI: 10.1007/978-1-0716-1716-8_2] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/13/2023]
Abstract
The formation of the transition state ensemble (TSE) represents the rate-limiting step in protein folding. The TSE is the least populated state on the pathway, and its characterization remains a challenge. Properties of the TSE can be inferred from the effects on folding and unfolding rates for various perturbations. A difficulty remains on how to translate these kinetic effects to structural properties of the TSE. Several factors can obscure the translation of point mutations in the frequently used method, "mutational Phi analysis." We take a complementary approach in "Psi analysis," employing rationally inserted metal binding sites designed to probe pairwise contacts in the TSE. These contacts can be confidently identified and used to construct structural models of the TSE. The method has been applied to multiple proteins and consistently produces a considerably more structured and native-like TSE than Phi analysis. This difference has significant implications to our understanding of protein folding mechanisms. Here we describe the application of the method and discuss how it can be used to study other conformational transitions such as binding.
Collapse
Affiliation(s)
- Michael C Baxa
- Department of Biochemistry and Molecular Biology, University of Chicago, Chicago, IL, USA
| | - Tobin R Sosnick
- Department of Biochemistry and Molecular Biology, University of Chicago, Chicago, IL, USA.
| |
Collapse
|
3
|
A Glimpse into the Structural Properties of the Intermediate and Transition State in the Folding of Bromodomain 2 Domain 2 by Φ Value Analysis. Int J Mol Sci 2021; 22:ijms22115953. [PMID: 34073056 PMCID: PMC8199192 DOI: 10.3390/ijms22115953] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/03/2021] [Revised: 05/21/2021] [Accepted: 05/27/2021] [Indexed: 11/23/2022] Open
Abstract
Bromodomains (BRDs) are small protein interaction modules of about 110 amino acids that selectively recognize acetylated lysine in histones and other proteins. These domains have been identified in a variety of multi-domain proteins involved in transcriptional regulation or chromatin remodeling in eukaryotic cells. BRD inhibition is considered an attractive therapeutic approach in epigenetic disorders, particularly in oncology. Here, we present a Φ value analysis to investigate the folding pathway of the second domain of BRD2 (BRD2(2)). Using an extensive mutational analysis based on 25 site-directed mutants, we provide structural information on both the intermediate and late transition state of BRD2(2). The data reveal that the C-terminal region represents part of the initial folding nucleus, while the N-terminal region of the domain consolidates its structure only later in the folding process. Furthermore, only a small number of native-like interactions have been identified, suggesting the presence of a non-compact, partially folded state with scarce native-like characteristics. Taken together, these results indicate that, in BRD2(2), a hierarchical mechanism of protein folding can be described with non-native interactions that play a significant role in folding.
Collapse
|
4
|
Bedford JT, Poutsma J, Diawara N, Greene LH. The nature of persistent interactions in two model β-grasp proteins reveals the advantage of symmetry in stability. J Comput Chem 2021; 42:600-607. [PMID: 33534913 DOI: 10.1002/jcc.26477] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/13/2020] [Revised: 10/27/2020] [Accepted: 11/22/2020] [Indexed: 01/25/2023]
Abstract
Two proteins within the β-grasp superfamily, the B1-domain of protein G and the small archaeal modifier protein 1, were investigated to elucidate the key determinants of structural stability at the level of individual interactions. These symmetrical proteins both contain two β-hairpins which form a sheet flanked by a central α-helix. They were subjected to high temperature molecular dynamics simulations and the detailed behavior of each long-range interaction was characterized. The results revealed that in GB1 the most stable region was the C-terminal hairpin and in SAMP1 it was the opposite, the N-terminal hairpin. Experimental results for GB1 support this finding. In conclusion, it appears that the difference in the location and number of hydrophobic interactions dictate the differential stability which is accommodated due to structural symmetry of the β-grasp fold. Thus, the hairpins are interchangeable and in nature this lends itself to adaptability and flexibility.
Collapse
|
5
|
Validation of DBFOLD: An efficient algorithm for computing folding pathways of complex proteins. PLoS Comput Biol 2020; 16:e1008323. [PMID: 33196646 PMCID: PMC7704049 DOI: 10.1371/journal.pcbi.1008323] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/26/2020] [Revised: 11/30/2020] [Accepted: 10/17/2020] [Indexed: 11/19/2022] Open
Abstract
Atomistic simulations can provide valuable, experimentally-verifiable insights into protein folding mechanisms, but existing ab initio simulation methods are restricted to only the smallest proteins due to severe computational speed limits. The folding of larger proteins has been studied using native-centric potential functions, but such models omit the potentially crucial role of non-native interactions. Here, we present an algorithm, entitled DBFOLD, which can predict folding pathways for a wide range of proteins while accounting for the effects of non-native contacts. In addition, DBFOLD can predict the relative rates of different transitions within a protein’s folding pathway. To accomplish this, rather than directly simulating folding, our method combines equilibrium Monte-Carlo simulations, which deploy enhanced sampling, with unfolding simulations at high temperatures. We show that under certain conditions, trajectories from these two types of simulations can be jointly analyzed to compute unknown folding rates from detailed balance. This requires inferring free energies from the equilibrium simulations, and extrapolating transition rates from the unfolding simulations to lower, physiologically-reasonable temperatures at which the native state is marginally stable. As a proof of principle, we show that our method can accurately predict folding pathways and Monte-Carlo rates for the well-characterized Streptococcal protein G. We then show that our method significantly reduces the amount of computation time required to compute the folding pathways of large, misfolding-prone proteins that lie beyond the reach of existing direct simulation. Our algorithm, which is available online, can generate detailed atomistic models of protein folding mechanisms while shedding light on the role of non-native intermediates which may crucially affect organismal fitness and are frequently implicated in disease. Many proteins must adopt a specific structure in order to function. Computational simulations have been used to shed light on the mechanisms of protein folding, but unfortunately, realistic simulations can typically only be run for small proteins, due to severe limits in computational speed. Here, we present a method to solve this problem, whereby instead of directly simulating folding from an unfolded state, we run simulations that allow for computation of equilibrium folding free energies, alongside high temperature simulations to compute unfolding rates. From these quantities, folding rates can be computed using detailed balance. Importantly, our method can account for the effects of nonnative contacts which transiently form during folding and must be broken prior to adoption of the native state. Such contacts, which are often excluded from simple models of folding, may crucially affect real protein folding pathways and are often observed in folding intermediates implicated in disease.
Collapse
|
6
|
Becerra D, Butyaev A, Waldispühl J. Fast and flexible coarse-grained prediction of protein folding routes using ensemble modeling and evolutionary sequence variation. Bioinformatics 2020; 36:1420-1428. [PMID: 31584628 DOI: 10.1093/bioinformatics/btz743] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/02/2019] [Revised: 09/22/2019] [Accepted: 09/28/2019] [Indexed: 11/15/2022] Open
Abstract
MOTIVATION Protein folding is a dynamic process through which polypeptide chains reach their native 3D structures. Although the importance of this mechanism is widely acknowledged, very few high-throughput computational methods have been developed to study it. RESULTS In this paper, we report a computational platform named P3Fold that combines statistical and evolutionary information for predicting and analyzing protein folding routes. P3Fold uses coarse-grained modeling and efficient combinatorial schemes to predict residue contacts and evaluate the folding routes of a protein sequence within minutes or hours. To facilitate access to this technology, we devise graphical representations and implement an interactive web interface that allows end-users to leverage P3Fold predictions. Finally, we use P3Fold to conduct large and short scale experiments on the human proteome that reveal the broad conservation and variations of structural intermediates within protein families. AVAILABILITY AND IMPLEMENTATION A Web server of P3Fold is freely available at http://csb.cs.mcgill.ca/P3Fold. SUPPLEMENTARY INFORMATION Supplementary data are available at Bioinformatics online.
Collapse
Affiliation(s)
- David Becerra
- School of Computer Science, McGill University, Montréal, QC H3A 0E9, Canada
| | - Alexander Butyaev
- School of Computer Science, McGill University, Montréal, QC H3A 0E9, Canada
| | - Jérôme Waldispühl
- School of Computer Science, McGill University, Montréal, QC H3A 0E9, Canada
| |
Collapse
|
7
|
Clark PL, Plaxco KW, Sosnick TR. Water as a Good Solvent for Unfolded Proteins: Folding and Collapse are Fundamentally Different. J Mol Biol 2020; 432:2882-2889. [PMID: 32044346 DOI: 10.1016/j.jmb.2020.01.031] [Citation(s) in RCA: 16] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/22/2019] [Revised: 01/28/2020] [Accepted: 01/29/2020] [Indexed: 12/30/2022]
Abstract
The argument that the hydrophobic effect is the primary effect driving the folding of globular proteins is nearly universally accepted (including by the authors). But does this view also imply that water is a "poor" solvent for the unfolded states of these same proteins? Here we argue that the answer is "no," that is, folding to a well-packed, extensively hydrogen-bonded native structure differs fundamentally from the nonspecific chain collapse that defines a poor solvent. Thus, the observation that a protein folds in water does not necessitate that water is a poor solvent for its unfolded state. Indeed, chain-solvent interactions that are marginally more favorable than nonspecific intrachain interactions are beneficial to protein function because they destabilize deleterious misfolded conformations and inter-chain interactions.
Collapse
Affiliation(s)
- Patricia L Clark
- Department of Chemistry and Biochemistry, University of Notre Dame, Notre Dame, IN, 46556, USA.
| | - Kevin W Plaxco
- Department of Chemistry and Biochemistry, University of California, Santa Barbara, CA, 93106, USA.
| | - Tobin R Sosnick
- Department of Biochemistry and Molecular Biology, Institute for Biophysical Dynamics, Pritzker School of Molecular Engineering, University of Chicago, Chicago, IL, 60637, USA.
| |
Collapse
|
8
|
Cheng Q, Joung I, Lee J, Kuwajima K, Lee J. Exploring the Folding Mechanism of Small Proteins GB1 and LB1. J Chem Theory Comput 2019; 15:3432-3449. [DOI: 10.1021/acs.jctc.8b01163] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Qianyi Cheng
- Department of Chemistry, University of Memphis, Memphis, Tennessee 38152, United States
- School of Computational Sciences, Korea Institute for Advanced Study, Seoul 02455, South Korea
| | - InSuk Joung
- Department of Chemistry, Kangwon National University, Chuncheon 24341, South Korea
- School of Computational Sciences, Korea Institute for Advanced Study, Seoul 02455, South Korea
| | - Juyong Lee
- Department of Chemistry, Kangwon National University, Chuncheon 24341, South Korea
| | - Kunihiro Kuwajima
- Department of Physics, University of Tokyo, Tokyo 113-0033, Japan
- School of Computational Sciences, Korea Institute for Advanced Study, Seoul 02455, South Korea
| | - Jooyoung Lee
- Center for In Silico Protein Science, Korea Institute for Advanced Study, Seoul 02455, South Korea
- School of Computational Sciences, Korea Institute for Advanced Study, Seoul 02455, South Korea
| |
Collapse
|
9
|
Zerze GH, Zheng W, Best RB, Mittal J. Evolution of All-Atom Protein Force Fields to Improve Local and Global Properties. J Phys Chem Lett 2019; 10:2227-2234. [PMID: 30990694 PMCID: PMC7507668 DOI: 10.1021/acs.jpclett.9b00850] [Citation(s) in RCA: 50] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/06/2023]
Abstract
Experimental studies on intrinsically disordered and unfolded proteins have shown that in isolation they typically have low populations of secondary structure and exhibit distance scalings suggesting that they are at near-theta-solvent conditions. Until recently, however, all-atom force fields failed to reproduce these fundamental properties of intrinsically disordered proteins (IDPs). Recent improvements by refining against ensemble-averaged experimental observables for polypeptides in aqueous solution have addressed deficiencies including secondary structure bias, global conformational properties, and thermodynamic parameters of biophysical reactions such as folding and collapse. To date, studies utilizing these improved all-atom force fields have mostly been limited to a small set of unfolded or disordered proteins. Here, we present data generated for a diverse library of unfolded or disordered proteins using three progressively improved generations of Amber03 force fields, and we explore how global and local properties are affected by each successive change in the force field. We find that the most recent force field refinements significantly improve the agreement of the global properties such as radii of gyration and end-to-end distances with experimental estimates. However, these global properties are largely independent of the local secondary structure propensity. This result stresses the need to validate force fields with reference to a combination of experimental data providing information about both local and global structure formation.
Collapse
Affiliation(s)
- Gül H Zerze
- Department of Chemical and Biomolecular Engineering , Lehigh University , Bethlehem , Pennsylvania 18015 , United States
| | - Wenwei Zheng
- College of Integrative Sciences and Arts , Arizona State University , Mesa , Arizona 85212 , United States
| | - Robert B Best
- Laboratory of Chemical Physics, National Institute of Diabetes and Digestive and Kidney Diseases , National Institutes of Health , Bethesda , Maryland 20892 , United States
| | - Jeetain Mittal
- Department of Chemical and Biomolecular Engineering , Lehigh University , Bethlehem , Pennsylvania 18015 , United States
| |
Collapse
|
10
|
Abstract
The folding simulations of three ββα-motifs and β-barrel structured proteins (NTL9, NuG2b, and CspA) were performed to determine the important roles of native and nonnative contacts in protein folding.
Collapse
Affiliation(s)
- Qiang Shao
- Drug Discovery and Design Center
- CAS Key Laboratory of Receptor Research
- Shanghai Institute of Materia Medica
- Chinese Academy of Sciences
- Shanghai
| | - Weiliang Zhu
- Drug Discovery and Design Center
- CAS Key Laboratory of Receptor Research
- Shanghai Institute of Materia Medica
- Chinese Academy of Sciences
- Shanghai
| |
Collapse
|
11
|
Ge B, Jiang X, Chen Y, Sun T, Yang Q, Huang F. Kinetic and thermodynamic studies reveal chemokine homologues CC11 and CC24 with an almost identical tertiary structure have different folding pathways. BMC BIOPHYSICS 2017; 10:7. [PMID: 28919974 PMCID: PMC5596964 DOI: 10.1186/s13628-017-0039-4] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 03/31/2017] [Accepted: 09/06/2017] [Indexed: 11/10/2022]
Abstract
BACKGROUND Proteins with low sequence identity but almost identical tertiary structure and function have been valuable to uncover the relationship between sequence, tertiary structure, folding mechanism and functions. Two homologous chemokines, CCL11 and CCL24, with low sequence identity but similar tertiary structure and function, provide an excellent model system for respective studies. RESULTS The kinetics and thermodynamics of the two homologous chemokines were systematically characterized. Despite their similar tertiary structures, CCL11 and CCL24 show different thermodynamic stability in guanidine hydrochloride titration, with D50% = 2.20 M and 4.96 M, respectively. The kinetics curves clearly show two phases in the folding/unfolding processes of both CCL11 and CCL24, which suggests the existence of an intermediate state in their folding/unfolding processes. The folding pathway of both CCL11 and CCL24 could be well described using a folding model with an on-pathway folding intermediate. However, the folding kinetics and stability of the intermediate state of CCL11 and CCL24 are obviously different. CONCLUSION Our results suggest homologous proteins with low sequence identity can display almost identical tertiary structure, but very different folding mechanisms, which applies to homologues in the chemokine protein family, extending the general applicability of the above observation.
Collapse
Affiliation(s)
- Baosheng Ge
- Center for Bioengineering and Biotechnology, China University of Petroleum (East China), Qingdao, 266580 People's Republic of China
| | - Xiaoyong Jiang
- Center for Bioengineering and Biotechnology, China University of Petroleum (East China), Qingdao, 266580 People's Republic of China
| | - Yao Chen
- Center for Bioengineering and Biotechnology, China University of Petroleum (East China), Qingdao, 266580 People's Republic of China
| | - Tingting Sun
- Center for Bioengineering and Biotechnology, China University of Petroleum (East China), Qingdao, 266580 People's Republic of China
| | - Qiuxia Yang
- Center for Bioengineering and Biotechnology, China University of Petroleum (East China), Qingdao, 266580 People's Republic of China
| | - Fang Huang
- Center for Bioengineering and Biotechnology, China University of Petroleum (East China), Qingdao, 266580 People's Republic of China
| |
Collapse
|
12
|
Jacobs WM, Shakhnovich EI. Structure-Based Prediction of Protein-Folding Transition Paths. Biophys J 2017; 111:925-36. [PMID: 27602721 PMCID: PMC5018131 DOI: 10.1016/j.bpj.2016.06.031] [Citation(s) in RCA: 21] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/22/2016] [Revised: 06/08/2016] [Accepted: 06/27/2016] [Indexed: 12/24/2022] Open
Abstract
We propose a general theory to describe the distribution of protein-folding transition paths. We show that transition paths follow a predictable sequence of high-free-energy transient states that are separated by free-energy barriers. Each transient state corresponds to the assembly of one or more discrete, cooperative units, which are determined directly from the native structure. We show that the transition state on a folding pathway is reached when a small number of critical contacts are formed between a specific set of substructures, after which folding proceeds downhill in free energy. This approach suggests a natural resolution for distinguishing parallel folding pathways and provides a simple means to predict the rate-limiting step in a folding reaction. Our theory identifies a common folding mechanism for proteins with diverse native structures and establishes general principles for the self-assembly of polymers with specific interactions.
Collapse
Affiliation(s)
- William M Jacobs
- Department of Chemistry and Chemical Biology, Harvard University, Cambridge, Massachusetts
| | - Eugene I Shakhnovich
- Department of Chemistry and Chemical Biology, Harvard University, Cambridge, Massachusetts.
| |
Collapse
|
13
|
Shao Q, Shi J, Zhu W. Determining Protein Folding Pathway and Associated Energetics through Partitioned Integrated-Tempering-Sampling Simulation. J Chem Theory Comput 2017; 13:1229-1243. [DOI: 10.1021/acs.jctc.6b00967] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/13/2023]
Affiliation(s)
- Qiang Shao
- Drug
Discovery and Design Center, CAS Key Laboratory of Receptor Research,
Shanghai Institute of Materia Medica, Chinese Academy of Sciences, 555 Zuchongzhi Road, Shanghai 201203, China
| | - Jiye Shi
- UCB Biopharma
SPRL, Chemin du Foriest, 1420 Braine-l’Alleud, Belgium
| | - Weiliang Zhu
- Drug
Discovery and Design Center, CAS Key Laboratory of Receptor Research,
Shanghai Institute of Materia Medica, Chinese Academy of Sciences, 555 Zuchongzhi Road, Shanghai 201203, China
| |
Collapse
|
14
|
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
|
15
|
Collins JC, Bedford JT, Greene LH. Elucidating the Key Determinants of Structure, Folding, and Stability for the ( 4β+ α ) Conformation of the B1 Domain of Protein G Using Bioinformatics Approaches. IEEE Trans Nanobioscience 2016; 15:140-7. [PMID: 27071185 DOI: 10.1109/tnb.2016.2546247] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/10/2022]
Abstract
The B1 domain of protein G (GB1) is a small, 56 amino acid bacterial immunoglobulin-binding protein with a 4β+ α fold. Architecturally, it is composed of a two-layer sandwich consisting of a four-stranded β -sheet that packs against an α -helix. Using several bioinformatics approaches, we investigated which residues may be key determinants of this fold. We identified nine structurally conserved amino acids using a conservation analysis and propose they are critical to forming and stabilizing the fold. The nine conserved residues form a predominantly hydrophobic nucleus within the core of GB1. A network analysis of all the long-range interactions in the structure of GB1 in concert with a betweenness centrality analysis revealed the relative significance of each conserved amino acid residue based on the number and location of the interactions. This bioinformatics analysis provides an important foundation for the design and interpretation of both computational and experimental work which may be helpful in solving the protein folding problem.
Collapse
|
16
|
Cooperative folding near the downhill limit determined with amino acid resolution by hydrogen exchange. Proc Natl Acad Sci U S A 2016; 113:4747-52. [PMID: 27078098 DOI: 10.1073/pnas.1522500113] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022] Open
Abstract
The relationship between folding cooperativity and downhill, or barrier-free, folding of proteins under highly stabilizing conditions remains an unresolved topic, especially for proteins such as λ-repressor that fold on the microsecond timescale. Under aqueous conditions where downhill folding is most likely to occur, we measure the stability of multiple H bonds, using hydrogen exchange (HX) in a λYA variant that is suggested to be an incipient downhill folder having an extrapolated folding rate constant of 2 × 10(5) s(-1) and a stability of 7.4 kcal·mol(-1) at 298 K. At least one H bond on each of the three largest helices (α1, α3, and α4) breaks during a common unfolding event that reflects global denaturation. The use of HX enables us to both examine folding under highly stabilizing, native-like conditions and probe the pretransition state region for stable species without the need to initiate the folding reaction. The equivalence of the stability determined at zero and high denaturant indicates that any residual denatured state structure minimally affects the stability even under native conditions. Using our ψ analysis method along with mutational ϕ analysis, we find that the three aforementioned helices are all present in the folding transition state. Hence, the free energy surface has a sufficiently high barrier separating the denatured and native states that folding appears cooperative even under extremely stable and fast folding conditions.
Collapse
|
17
|
Microscopic interpretation of folding ϕ-values using the transition path ensemble. Proc Natl Acad Sci U S A 2016; 113:3263-8. [PMID: 26957599 DOI: 10.1073/pnas.1520864113] [Citation(s) in RCA: 44] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
All-atom molecular dynamics simulations now allow us to create movies of proteins folding and unfolding. However, it is difficult to assess the accuracy of the folding mechanisms observed because experiments cannot yet directly resolve events occurring along the transition paths between unfolded and folded states. Protein folding ϕ-values provide residue-resolved information about folding mechanisms by comparing the effects of mutations on folding rates and stability, but determining ϕ-values by separately simulating mutant proteins would be computationally demanding and prone to large statistical errors. Here we use transition path theory to develop a method for computing ϕ-values directly from the transition path ensemble, without the need for additional simulations. This path-based approach uses the full transition path information available from equilibrium folding and unfolding trajectories, or from transition path sampling, and does not require identification of folding transition states. Applying our approach to a set of simulations of 10 small proteins by Shaw and coworkers [Lindorff-Larsen K, Piana S, Dror RO, Shaw DE (2011) Science 334(6055):517-520; Piana S, Lindorff-Larsen K, Shaw DE (2011) Biophys J100(9):L47-L49; and Piana S, Lindorff-Larsen K, Shaw DE (2013) Proc Natl Acad Sci USA 110(15):5915-5920], we find good agreement with experiments in most cases where data are available. We can further resolve the contributions to fractional ϕ-values coming from partial contact formation versus transition path heterogeneity. Although in some cases, there is substantial heterogeneity of folding mechanism, in others, such as Ubiquitin, the mechanism is strongly conserved.
Collapse
|
18
|
Maity H, Reddy G. Folding of Protein L with Implications for Collapse in the Denatured State Ensemble. J Am Chem Soc 2016; 138:2609-16. [PMID: 26835789 DOI: 10.1021/jacs.5b11300] [Citation(s) in RCA: 38] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
A fundamental question in protein folding is whether the coil to globule collapse transition occurs during the initial stages of folding (burst phase) or simultaneously with the protein folding transition. Single molecule fluorescence resonance energy transfer (FRET) and small-angle X-ray scattering (SAXS) experiments disagree on whether Protein L collapse transition occurs during the burst phase of folding. We study Protein L folding using a coarse-grained model and molecular dynamics simulations. The collapse transition in Protein L is found to be concomitant with the folding transition. In the burst phase of folding, we find that FRET experiments overestimate radius of gyration, Rg, of the protein due to the application of Gaussian polymer chain end-to-end distribution to extract Rg from the FRET efficiency. FRET experiments estimate ≈6 Å decrease in Rg when the actual decrease is ≈3 Å on guanidinium chloride denaturant dilution from 7.5 to 1 M, thereby suggesting pronounced compaction in the protein dimensions in the burst phase. The ≈3 Å decrease is close to the statistical uncertainties of the Rg data measured from SAXS experiments, which suggest no compaction, leading to a disagreement with the FRET experiments. The transition-state ensemble (TSE) structures in Protein L folding are globular and extensive in agreement with the Ψ-analysis experiments. The results support the hypothesis that the TSE of single domain proteins depends on protein topology and is not stabilized by local interactions alone.
Collapse
Affiliation(s)
- Hiranmay Maity
- Solid State and Structural Chemistry Unit, Indian Institute of Science , Bangalore, Karnataka 560012, India
| | - Govardhan Reddy
- Solid State and Structural Chemistry Unit, Indian Institute of Science , Bangalore, Karnataka 560012, India
| |
Collapse
|
19
|
Where the complex things are: single molecule and ensemble spectroscopic investigations of protein folding dynamics. Curr Opin Struct Biol 2015; 36:1-9. [PMID: 26687767 DOI: 10.1016/j.sbi.2015.11.006] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/18/2015] [Accepted: 11/10/2015] [Indexed: 01/11/2023]
Abstract
Progress in our understanding of the simple folding dynamics of small proteins and the complex dynamics of large proteins is reviewed. Recent characterizations of the folding transition path of small proteins revealed a simple dynamics explainable by the native centric model. In contrast, the accumulated data showed the substates containing residual structures in the unfolded state and partially populated intermediates, causing complexity in the early folding dynamics of small proteins. The size of the unfolded proteins in the absence of denaturants is likely expanded but still controversial. The steady progress in the observation of folding of large proteins has clarified the rapid formation of long-range contacts that seem inconsistent with the native centric model, suggesting that the folding strategy of large proteins is distinct from that of small proteins.
Collapse
|