1
|
Ooka K, Arai M. Accurate prediction of protein folding mechanisms by simple structure-based statistical mechanical models. Nat Commun 2023; 14:6338. [PMID: 37857633 PMCID: PMC10587348 DOI: 10.1038/s41467-023-41664-1] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/24/2022] [Accepted: 09/10/2023] [Indexed: 10/21/2023] Open
Abstract
Recent breakthroughs in highly accurate protein structure prediction using deep neural networks have made considerable progress in solving the structure prediction component of the 'protein folding problem'. However, predicting detailed mechanisms of how proteins fold into specific native structures remains challenging, especially for multidomain proteins constituting most of the proteomes. Here, we develop a simple structure-based statistical mechanical model that introduces nonlocal interactions driving the folding of multidomain proteins. Our model successfully predicts protein folding processes consistent with experiments, without the limitations of protein size and shape. Furthermore, slight modifications of the model allow prediction of disulfide-oxidative and disulfide-intact protein folding. These predictions depict details of the folding processes beyond reproducing experimental results and provide a rationale for the folding mechanisms. Thus, our physics-based models enable accurate prediction of protein folding mechanisms with low computational complexity, paving the way for solving the folding process component of the 'protein folding problem'.
Collapse
Affiliation(s)
- Koji Ooka
- Department of Physics, Graduate School of Science, The University of Tokyo, 3-8-1 Komaba, Meguro, Tokyo, 153-8902, Japan
- Komaba Organization for Educational Excellence, College of Arts and Sciences, The University of Tokyo, 3-8-1 Komaba, Meguro, Tokyo, 153-8902, Japan
| | - Munehito Arai
- Department of Physics, Graduate School of Science, The University of Tokyo, 3-8-1 Komaba, Meguro, Tokyo, 153-8902, Japan.
- Komaba Organization for Educational Excellence, College of Arts and Sciences, The University of Tokyo, 3-8-1 Komaba, Meguro, Tokyo, 153-8902, Japan.
- Department of Life Sciences, Graduate School of Arts and Sciences, The University of Tokyo, 3-8-1 Komaba, Meguro, Tokyo, 153-8902, Japan.
| |
Collapse
|
2
|
Ooka K, Liu R, Arai M. The Wako-Saitô-Muñoz-Eaton Model for Predicting Protein Folding and Dynamics. Molecules 2022; 27:molecules27144460. [PMID: 35889332 PMCID: PMC9319528 DOI: 10.3390/molecules27144460] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/20/2022] [Revised: 07/07/2022] [Accepted: 07/08/2022] [Indexed: 11/16/2022] Open
Abstract
Despite the recent advances in the prediction of protein structures by deep neutral networks, the elucidation of protein-folding mechanisms remains challenging. A promising theory for describing protein folding is a coarse-grained statistical mechanical model called the Wako-Saitô-Muñoz-Eaton (WSME) model. The model can calculate the free-energy landscapes of proteins based on a three-dimensional structure with low computational complexity, thereby providing a comprehensive understanding of the folding pathways and the structure and stability of the intermediates and transition states involved in the folding reaction. In this review, we summarize previous and recent studies on protein folding and dynamics performed using the WSME model and discuss future challenges and prospects. The WSME model successfully predicted the folding mechanisms of small single-domain proteins and the effects of amino-acid substitutions on protein stability and folding in a manner that was consistent with experimental results. Furthermore, extended versions of the WSME model were applied to predict the folding mechanisms of multi-domain proteins and the conformational changes associated with protein function. Thus, the WSME model may contribute significantly to solving the protein-folding problem and is expected to be useful for predicting protein folding, stability, and dynamics in basic research and in industrial and medical applications.
Collapse
Affiliation(s)
- Koji Ooka
- Department of Physics, Graduate School of Science, The University of Tokyo, 3-8-1 Komaba, Meguro, Tokyo 153-8902, Japan;
- Komaba Organization for Educational Excellence, College of Arts and Sciences, The University of Tokyo, 3-8-1 Komaba, Meguro, Tokyo 153-8902, Japan
| | - Runjing Liu
- Department of Life Sciences, Graduate School of Arts and Sciences, The University of Tokyo, 3-8-1 Komaba, Meguro, Tokyo 153-8902, Japan;
| | - Munehito Arai
- Department of Physics, Graduate School of Science, The University of Tokyo, 3-8-1 Komaba, Meguro, Tokyo 153-8902, Japan;
- Department of Life Sciences, Graduate School of Arts and Sciences, The University of Tokyo, 3-8-1 Komaba, Meguro, Tokyo 153-8902, Japan;
- Correspondence:
| |
Collapse
|
3
|
Naganathan AN, Kannan A. A hierarchy of coupling free energies underlie the thermodynamic and functional architecture of protein structures. Curr Res Struct Biol 2021; 3:257-267. [PMID: 34704074 PMCID: PMC8526763 DOI: 10.1016/j.crstbi.2021.09.003] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/07/2021] [Revised: 09/08/2021] [Accepted: 09/30/2021] [Indexed: 12/22/2022] Open
Abstract
Protein sequences and structures evolve by satisfying varied physical and biochemical constraints. This multi-level selection is enabled not just by the patterning of amino acids on the sequence, but also via coupling between residues in the native structure. Here, we employ an energetically detailed statistical mechanical model with millions of microstates to extract such long-range structural correlations, i.e. thermodynamic coupling free energies, from a diverse family of protein structures. We find that despite the intricate and anisotropic distribution of coupling patterns, the majority of residues (>70%) are only marginally coupled contributing to functional motions and catalysis. Physical origins of ‘sectors’, determinants of native ensemble heterogeneity in extant, ancient and designed proteins, and the basis for allostery emerge naturally from coupling free energies. The statistical framework highlights how evolutionary selection and optimization occur at the level of global interaction network for a given protein fold impacting folding, function, and allosteric outputs. Evolution of protein structures occurs at the level of global interaction network. More than 70% of the protein residues are weakly or marginally coupled. Functional ‘sector’ regions are a manifestation of marginal coupling. Coupling indices vary across the entire proteins in extant-ancient and natural-designed pairs. The proposed methodology can be used to understand allostery and epistasis.
Collapse
Affiliation(s)
- Athi N Naganathan
- Department of Biotechnology, Bhupat & Jyoti Mehta School of Biosciences, Indian Institute of Technology Madras, Chennai, 600036, India
| | - Adithi Kannan
- Department of Biotechnology, Bhupat & Jyoti Mehta School of Biosciences, Indian Institute of Technology Madras, Chennai, 600036, India
| |
Collapse
|
4
|
Gopi S, Naganathan AN. Non-specific DNA-driven quinary interactions promote structural transitions in proteins. Phys Chem Chem Phys 2020; 22:12671-12677. [PMID: 32458879 DOI: 10.1039/d0cp01758b] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Abstract
The nature and distribution of charged residues on the surface of proteins play a vital role in determining the binding affinity, selectivity and kinetics of association to ligands. When it comes to DNA-binding domains (DBDs), these functional features manifest as anisotropic distribution of positively charged residues on the protein surface driven by the requirement to bind DNA, a highly negatively charged polymer. In this work, we compare the thermodynamic behavior of nine different proteins belonging to three families - LacR, engrailed and Brk - some of which are disordered in solution in the absence of DNA. Combining detailed electrostatic calculations and statistical mechanical modeling of folding landscapes at different distances and relative orientations with respect to DNA, we show that non-specific electrostatic interactions between the protein and DNA can promote structural transitions in DBDs. Such quinary interactions that are strictly agnostic to the DNA sequence induce varied behaviors including folding of disordered domains, partial unfolding of ordered proteins and (de-)population of intermediate states. Our work highlights that the folding landscape of proteins can be tuned as a function of distance from DNA and hints at possible reasons for DBDs exhibiting complex kinetic-thermodynamic behaviors in the absence of DNA.
Collapse
Affiliation(s)
- Soundhararajan Gopi
- Department of Biotechnology, Bhupat & Jyoti Mehta School of Biosciences, Indian Institute of Technology Madras, Chennai 600036, India.
| | - Athi N Naganathan
- Department of Biotechnology, Bhupat & Jyoti Mehta School of Biosciences, Indian Institute of Technology Madras, Chennai 600036, India.
| |
Collapse
|
5
|
Abstract
This review discusses Gō models broadly used in biomolecular simulations. I start with a brief description of the original lattice model study by Nobuhiro Gō. Then, the theory of protein folding behind Gō model, free energy approaches, and off-lattice Gō models are reviewed. I also mention a stringent test for the assumption in Gō models given from all-atom molecular dynamics simulations. Subsequently, I move to application of Gō models to protein dynamical functions. Various extension of Gō models is also reviewed. Finally, some publicly available tools to use Gō models are listed.
Collapse
Affiliation(s)
- Shoji Takada
- Department of Biophysics, Graduate School of Science, Kyoto University, Kyoto 606-8502, Japan
| |
Collapse
|
6
|
Gopi S, Aranganathan A, Naganathan AN. Thermodynamics and folding landscapes of large proteins from a statistical mechanical model. Curr Res Struct Biol 2019; 1:6-12. [PMID: 34235463 PMCID: PMC8244504 DOI: 10.1016/j.crstbi.2019.10.002] [Citation(s) in RCA: 16] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/15/2019] [Revised: 10/10/2019] [Accepted: 10/13/2019] [Indexed: 01/01/2023] Open
Abstract
Statistical mechanical models that afford an intermediate resolution between macroscopic chemical models and all-atom simulations have been successful in capturing folding behaviors of many small single-domain proteins. However, the applicability of one such successful approach, the Wako-Saitô-Muñoz-Eaton (WSME) model, is limited by the size of the protein as the number of conformations grows exponentially with protein length. In this work, we surmount this size limitation by introducing a novel approximation that treats stretches of 3 or 4 residues as blocks, thus reducing the phase space by nearly three orders of magnitude. The performance of the 'bWSME' model is validated by comparing the predictions for a globular enzyme (RNase H) and a repeat protein (IκBα), against experimental observables and the model without block approximation. Finally, as a proof of concept, we predict the free-energy surface of the 370-residue, multi-domain maltose binding protein and identify an intermediate in good agreement with single-molecule force-spectroscopy measurements. The bWSME model can thus be employed as a quantitative predictive tool to explore the conformational landscapes of large proteins, extract the structural features of putative intermediates, identify parallel folding paths, and thus aid in the interpretation of both ensemble and single-molecule experiments.
Collapse
Affiliation(s)
- Soundhararajan Gopi
- Department of Biotechnology, Bhupat & Jyoti Mehta School of Biosciences, Indian Institute of Technology Madras, Chennai 600036, India
| | - Akashnathan Aranganathan
- Department of Biotechnology, Bhupat & Jyoti Mehta School of Biosciences, Indian Institute of Technology Madras, Chennai 600036, India
| | - Athi N Naganathan
- Department of Biotechnology, Bhupat & Jyoti Mehta School of Biosciences, Indian Institute of Technology Madras, Chennai 600036, India
| |
Collapse
|
7
|
|
8
|
Narayan A, Naganathan AN. Switching Protein Conformational Substates by Protonation and Mutation. J Phys Chem B 2018; 122:11039-11047. [PMID: 30048131 DOI: 10.1021/acs.jpcb.8b05108] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Protein modules that regulate the availability and conformational status of transcription factors determine the rapidity, duration, and magnitude of cellular response to changing conditions. One such system is the single-gene product Cnu, a four-helix bundle transcription co-repressor, which acts as a molecular thermosensor regulating the expression of virulence genes in enterobacteriaceae through modulation of its native conformational ensemble. Cnu and related genes have also been implicated in pH-dependent expression of virulence genes. We hypothesize that protonation of a conserved buried histidine (H45) in Cnu promotes large electrostatic frustration, thus disturbing the H-NS, a transcription factor, binding face. Spectroscopic and calorimetric methods reveal that H45 exhibits a suppressed p Ka of ∼5.1, the protonation of which switches the conformation to an alternate native ensemble in which the fourth helix is disordered. The population redistribution can also be achieved through a mutation H45V, which does not display any switching behavior at pH values greater than 4. The Wako-Saitô-Muñoz-Eaton (WSME) statistical mechanical model predicts specific differences in the conformations and fluctuations of the fourth and first helices of Cnu determining the observed pH response. We validate these predictions through fluorescence lifetime measurements of a sole tryptophan, highlighting the presence of both native and non-native interactions in the regions adjoining the binding face of Cnu. Our combined experimental-computational study thus shows that Cnu acts both as a thermo- and pH-sensor orchestrated via a subtle but quantifiable balance between the weak packing of a structural element and protonation of a buried histidine that promotes electrostatic frustration.
Collapse
Affiliation(s)
- Abhishek Narayan
- Department of Biotechnology, Bhupat & Jyoti Mehta School of Biosciences , Indian Institute of Technology Madras , Chennai 600036 , India
| | - Athi N Naganathan
- Department of Biotechnology, Bhupat & Jyoti Mehta School of Biosciences , Indian Institute of Technology Madras , Chennai 600036 , India
| |
Collapse
|
9
|
Grassein P, Delarue P, Scheraga HA, Maisuradze GG, Senet P. Statistical Model To Decipher Protein Folding/Unfolding at a Local Scale. J Phys Chem B 2018; 122:3540-3549. [PMID: 29446945 DOI: 10.1021/acs.jpcb.7b10733] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Protein folding/unfolding can be analyzed experimentally at a local scale by monitoring the physical properties of local probes as a function of the temperature, for example, the distance between fluorophores or the values of chemical shifts of backbone atoms. Here, the analytical Lifson-Roig model for the helix-coil transition is modified to analyze local thermal unfolding of the fast-folder W protein of bacteriophage lambda (gpW) simulated by all-atom molecular dynamics (MD) simulations in explicit solvent at 15 different temperatures. The protein structure is described by the coarse-grained dihedral angles (γ) and bond angles (θ) built between successive Cα-Cα virtual bonds. Each (γ,θ) pair serves as a local probe of protein unfolding. Local native/non-native states are defined for each pair of (γ,θ) angles by analyzing the free-energy landscapes Δ G(γ,θ) computed from MD trajectories. The three local elementary equilibrium constants of the model are extracted for each (γ,θ) pair along the sequence from MD simulations, and the model predictions are compared to the MD data. Using only the local equilibrium constants as an input, we show that the local denaturation curves computed from the model partition function fit their MD simulated counterparts in a satisfying manner without any adjustment. In the model and MD simulations, gpW unfolds gradually between 320 and 340 K, with an average native percentage decreasing from 0.8 (320 K) to 0.2 (340 K). In the prism of the model, there is no stable structure at the local scale in this 20 K unfolding temperature range. The enthalpy change upon local unfolding computed from the model and from MD trajectories suggests that the unfolded state between 320 and 340 K corresponds to a dynamical equilibrium between a large ensemble of constantly changing structures. The present results confirm the downhill unfolding of gpW, which does not obey a two-state global folding/unfolding model, and shed light on the interpretation of local denaturation curves.
Collapse
Affiliation(s)
- Paul Grassein
- Laboratoire Interdisciplinaire Carnot de Bourgogne , UMR 6303 CNRS-Univ. de Bourgogne Franche-Comté , 9 Av. A. Savary, BP 47 870 , F-21078 Dijon Cedex , France
| | - Patrice Delarue
- Laboratoire Interdisciplinaire Carnot de Bourgogne , UMR 6303 CNRS-Univ. de Bourgogne Franche-Comté , 9 Av. A. Savary, BP 47 870 , F-21078 Dijon Cedex , France
| | - Harold A Scheraga
- Baker Laboratory of Chemistry and Chemical Biology , Cornell University , Ithaca , New York 14853-1301 , United States
| | - Gia G Maisuradze
- Baker Laboratory of Chemistry and Chemical Biology , Cornell University , Ithaca , New York 14853-1301 , United States
| | - Patrick Senet
- Laboratoire Interdisciplinaire Carnot de Bourgogne , UMR 6303 CNRS-Univ. de Bourgogne Franche-Comté , 9 Av. A. Savary, BP 47 870 , F-21078 Dijon Cedex , France.,Baker Laboratory of Chemistry and Chemical Biology , Cornell University , Ithaca , New York 14853-1301 , United States
| |
Collapse
|
10
|
Gopi S, Singh A, Suresh S, Paul S, Ranu S, Naganathan AN. Toward a quantitative description of microscopic pathway heterogeneity in protein folding. Phys Chem Chem Phys 2017; 19:20891-20903. [DOI: 10.1039/c7cp03011h] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Experimentally consistent statistical modeling of protein folding thermodynamics reveals unprecedented complexity with numerous parallel folding routes in five different proteins.
Collapse
Affiliation(s)
- Soundhararajan Gopi
- Department of Biotechnology
- Bhupat & Jyoti Mehta School of Biosciences
- Indian Institute of Technology Madras
- Chennai 600036
- India
| | - Animesh Singh
- Department of Computer Science and Engineering
- Indian Institute of Technology Madras
- Chennai 600036
- India
| | | | - Suvadip Paul
- Department of Computer Science and Engineering
- Indian Institute of Technology Madras
- Chennai 600036
- India
| | - Sayan Ranu
- Department of Computer Science and Engineering
- Indian Institute of Technology Madras
- Chennai 600036
- India
| | - Athi N. Naganathan
- Department of Biotechnology
- Bhupat & Jyoti Mehta School of Biosciences
- Indian Institute of Technology Madras
- Chennai 600036
- India
| |
Collapse
|
11
|
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
|
12
|
Gopi S, Rajasekaran N, Singh A, Ranu S, Naganathan AN. Energetic and topological determinants of a phosphorylation-induced disorder-to-order protein conformational switch. Phys Chem Chem Phys 2016; 17:27264-9. [PMID: 26421497 DOI: 10.1039/c5cp04765j] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
We show that the phosphorylation of 4E-BP2 acts as a triggering event to shape its folding-function landscape that is delicately balanced between conflicting favorable energetics and intrinsically unfavorable topological connectivity. We further provide first evidence that the fitness landscapes of proteins at the threshold of disorder can differ considerably from ordered domains.
Collapse
Affiliation(s)
- Soundhararajan Gopi
- Department of Biotechnology, Bhupat & Jyoti Mehta School of Biosciences, Indian Institute of Technology Madras, Chennai 600036, India.
| | | | | | | | | |
Collapse
|
13
|
Rajasekaran N, Gopi S, Narayan A, Naganathan AN. Quantifying Protein Disorder through Measures of Excess Conformational Entropy. J Phys Chem B 2016; 120:4341-50. [PMID: 27111521 DOI: 10.1021/acs.jpcb.6b00658] [Citation(s) in RCA: 29] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Intrinsically disordered proteins (IDPs) and proteins with a large degree of disorder are abundant in the proteomes of eukaryotes and viruses, and play a vital role in cellular homeostasis and disease. One fundamental question that has been raised on IDPs is the process by which they offset the entropic penalty involved in transitioning from a heterogeneous ensemble of conformations to a much smaller collection of binding-competent states. However, this has been a difficult problem to address, as the effective entropic cost of fixing residues in a folded-like conformation from disordered amino acid neighborhoods is itself not known. Moreover, there are several examples where the sequence complexity of disordered regions is as high as well-folded regions. Disorder in such cases therefore arises from excess conformational entropy determined entirely by correlated sequence effects, an entropic code that is yet to be identified. Here, we explore these issues by exploiting the order-disorder transitions of a helix in Pbx-Homeodomain together with a dual entropy statistical mechanical model to estimate the magnitude and sign of the excess conformational entropy of residues in disordered regions. We find that a mere 2.1-fold increase in the number of allowed conformations per residue (∼0.7kBT favoring the unfolded state) relative to a well-folded sequence, or ∼2(N) additional conformations for a N-residue sequence, is sufficient to promote disorder under physiological conditions. We show that this estimate is quite robust and helps in rationalizing the thermodynamic signatures of disordered regions in important regulatory proteins, modeling the conformational folding-binding landscapes of IDPs, quantifying the stability effects characteristic of disordered protein loops and their subtle roles in determining the partitioning of folding flux in ordered domains. In effect, the dual entropy model we propose provides a statistical thermodynamic basis for the relative conformational propensities of amino acids in folded and disordered environments in proteins. Our work thus lays the foundation for understanding and quantifying protein disorder through measures of excess conformational entropy.
Collapse
Affiliation(s)
- Nandakumar Rajasekaran
- Department of Biotechnology, Bhupat & Jyoti Mehta School of Biosciences, Indian Institute of Technology Madras , Chennai 600036, India
| | - Soundhararajan Gopi
- Department of Biotechnology, Bhupat & Jyoti Mehta School of Biosciences, Indian Institute of Technology Madras , Chennai 600036, India
| | - Abhishek Narayan
- Department of Biotechnology, Bhupat & Jyoti Mehta School of Biosciences, Indian Institute of Technology Madras , Chennai 600036, India
| | - Athi N Naganathan
- Department of Biotechnology, Bhupat & Jyoti Mehta School of Biosciences, Indian Institute of Technology Madras , Chennai 600036, India
| |
Collapse
|
14
|
Munshi S, Naganathan AN. Imprints of function on the folding landscape: functional role for an intermediate in a conserved eukaryotic binding protein. Phys Chem Chem Phys 2016; 17:11042-52. [PMID: 25824585 DOI: 10.1039/c4cp06102k] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/15/2023]
Abstract
In the computational characterization of single domain protein folding, the effective free energies of numerous microstates are projected onto few collective degrees of freedom that in turn serve as well-defined reaction coordinates. In this regard, one-dimensional (1D) free energy profiles are widely used mainly for their simplicity. Since folding and functional landscapes are interlinked, how well can these reduced representations capture the structural and dynamic features of functional states while being simultaneously consistent with experimental observables? We investigate this issue by characterizing the folding of the four-helix bundle bovine acyl-CoA binding protein (bACBP), which exhibits complex equilibrium and kinetic behaviours, employing an Ising-like statistical mechanical model and molecular simulations. We show that the features of the 1D free energy profile are sufficient to quantitatively reproduce multiple experimental observations including millisecond chevron-like kinetics and temperature dependence, a microsecond fast phase, barrier heights, unfolded state movements, the intermediate structure and average ϕ-values. Importantly, we find that the structural features of the native-like intermediate (partial disorder in helix 1) are intricately linked to a unique interplay between packing and electrostatics in this domain. By comparison with available experimental data, we propose that this intermediate determines the promiscuous functional behaviour of bACBP that exhibits broad substrate specificity. Our results present evidence to the possibility of employing the statistical mechanical model and the resulting 1D free energy profile to not just understand folding mechanisms but to even extract features of functionally relevant states and their energetic origins.
Collapse
Affiliation(s)
- Sneha Munshi
- Department of Biotechnology, Bhupat & Jyoti Mehta School of Biosciences, Indian Institute of Technology Madras, Chennai 600036, India.
| | | |
Collapse
|
15
|
van Dijk E, Varilly P, Knowles TPJ, Frenkel D, Abeln S. Consistent Treatment of Hydrophobicity in Protein Lattice Models Accounts for Cold Denaturation. PHYSICAL REVIEW LETTERS 2016; 116:078101. [PMID: 26943560 DOI: 10.1103/physrevlett.116.078101] [Citation(s) in RCA: 24] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/24/2015] [Indexed: 05/04/2023]
Abstract
The hydrophobic effect stabilizes the native structure of proteins by minimizing the unfavorable interactions between hydrophobic residues and water through the formation of a hydrophobic core. Here, we include the entropic and enthalpic contributions of the hydrophobic effect explicitly in an implicit solvent model. This allows us to capture two important effects: a length-scale dependence and a temperature dependence for the solvation of a hydrophobic particle. This consistent treatment of the hydrophobic effect explains cold denaturation and heat capacity measurements of solvated proteins.
Collapse
Affiliation(s)
- Erik van Dijk
- Department of Chemistry, University of Cambridge, Cambridge CB2 1EW, United Kingdom
- Centre for Integrative Bioinformatics (IBIVU), Vrije Universiteit, De Boelelaan 1081A, 1081 HV Amsterdam, Netherlands
| | - Patrick Varilly
- Department of Chemistry, University of Cambridge, Cambridge CB2 1EW, United Kingdom
| | - Tuomas P J Knowles
- Department of Chemistry, University of Cambridge, Cambridge CB2 1EW, United Kingdom
| | - Daan Frenkel
- Department of Chemistry, University of Cambridge, Cambridge CB2 1EW, United Kingdom
| | - Sanne Abeln
- Centre for Integrative Bioinformatics (IBIVU), Vrije Universiteit, De Boelelaan 1081A, 1081 HV Amsterdam, Netherlands
| |
Collapse
|
16
|
Naganathan AN, De Sancho D. Bridging Experiments and Native-Centric Simulations of a Downhill Folding Protein. J Phys Chem B 2015; 119:14925-33. [DOI: 10.1021/acs.jpcb.5b09568] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Affiliation(s)
- Athi N. Naganathan
- Department of Biotechnology, Bhupat & Jyoti Mehta School of Biosciences, Indian Institute of Technology Madras, Chennai 600036, India
| | - David De Sancho
- CIC nanoGUNE, Tolosa Hiribidea,
76, E-20018 Donostia-San
Sebastián, Spain
- IKERBASQUE, Basque Foundation for Science, María Díaz de Haro 3, 48013 Bilbao, Spain
| |
Collapse
|
17
|
Hutton RD, Wilkinson J, Faccin M, Sivertsson EM, Pelizzola A, Lowe AR, Bruscolini P, Itzhaki LS. Mapping the Topography of a Protein Energy Landscape. J Am Chem Soc 2015; 137:14610-25. [PMID: 26561984 DOI: 10.1021/jacs.5b07370] [Citation(s) in RCA: 23] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Abstract
Protein energy landscapes are highly complex, yet the vast majority of states within them tend to be invisible to experimentalists. Here, using site-directed mutagenesis and exploiting the simplicity of tandem-repeat protein structures, we delineate a network of these states and the routes between them. We show that our target, gankyrin, a 226-residue 7-ankyrin-repeat protein, can access two alternative (un)folding pathways. We resolve intermediates as well as transition states, constituting a comprehensive series of snapshots that map early and late stages of the two pathways and show both to be polarized such that the repeat array progressively unravels from one end of the molecule or the other. Strikingly, we find that the protein folds via one pathway but unfolds via a different one. The origins of this behavior can be rationalized using the numerical results of a simple statistical mechanics model that allows us to visualize the equilibrium behavior as well as single-molecule folding/unfolding trajectories, thereby filling in the gaps that are not accessible to direct experimental observation. Our study highlights the complexity of repeat-protein folding arising from their symmetrical structures; at the same time, however, this structural simplicity enables us to dissect the complexity and thereby map the precise topography of the energy landscape in full breadth and remarkable detail. That we can recapitulate the key features of the folding mechanism by computational analysis of the native structure alone will help toward the ultimate goal of designed amino-acid sequences with made-to-measure folding mechanisms-the Holy Grail of protein folding.
Collapse
Affiliation(s)
- Richard D Hutton
- Hutchison/MRC Research Centre , Hills Road, Cambridge CB2 0XZ, U.K
| | - James Wilkinson
- Hutchison/MRC Research Centre , Hills Road, Cambridge CB2 0XZ, U.K
| | - Mauro Faccin
- ICTEAM, Université Catholique de Lovain , Euler Building 4, Avenue Lemaître, B-1348 Louvain-la-Neuve, Belgium
| | - Elin M Sivertsson
- Department of Pharmacology, University of Cambridge , Tennis Court Road, Cambridge CB2 1PD, U.K
| | - Alessandro Pelizzola
- Dipartimento di Scienza Applicata e Tecnologia, CNISM, and Center for Computational Studies, Politecnico di Torino , Corso Duca degli Abruzzi 24, I-10129 Torino, Italy.,INFN, Sezione di Torino , via Pietro Giuria 1, I-10125 Torino, Italy.,Human Genetics Foundation (HuGeF) , Via Nizza 52, I-10126 Torino, Italy
| | - Alan R Lowe
- Institute of Structural and Molecular Biology and London Centre for Nanotechnology, University College London and Birkbeck College , London WC1E 7HX, U.K
| | - Pierpaolo Bruscolini
- Departamento de Física Teórica and Instituto de Biocomputacíon y Física de Sistemas Complejos (BIFI), Universidad de Zaragoza , c/Mariano Esquillor s/n, 50018 Zaragoza, Spain
| | - Laura S Itzhaki
- Department of Pharmacology, University of Cambridge , Tennis Court Road, Cambridge CB2 1PD, U.K
| |
Collapse
|
18
|
Naganathan AN, Sanchez-Ruiz JM, Munshi S, Suresh S. Are Protein Folding Intermediates the Evolutionary Consequence of Functional Constraints? J Phys Chem B 2015; 119:1323-33. [DOI: 10.1021/jp510342m] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/03/2023]
Affiliation(s)
- Athi N. Naganathan
- Department of Biotechnology, Bhupat & Jyoti Mehta School of Biosciences, Indian Institute of Technology Madras, Chennai 600036, India
| | - Jose M. Sanchez-Ruiz
- Departamento de Quimica Fisica,
Facultad de Ciencias, Universidad de Granada, 18071 Granada, Spain
| | - Sneha Munshi
- Department of Biotechnology, Bhupat & Jyoti Mehta School of Biosciences, Indian Institute of Technology Madras, Chennai 600036, India
| | | |
Collapse
|
19
|
Rollins GC, Dill KA. General mechanism of two-state protein folding kinetics. J Am Chem Soc 2014; 136:11420-7. [PMID: 25056406 DOI: 10.1021/ja5049434] [Citation(s) in RCA: 50] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Abstract
We describe here a general model of the kinetic mechanism of protein folding. In the Foldon Funnel Model, proteins fold in units of secondary structures, which form sequentially along the folding pathway, stabilized by tertiary interactions. The model predicts that the free energy landscape has a volcano shape, rather than a simple funnel, that folding is two-state (single-exponential) when secondary structures are intrinsically unstable, and that each structure along the folding path is a transition state for the previous structure. It shows how sequential pathways are consistent with multiple stochastic routes on funnel landscapes, and it gives good agreement with the 9 order of magnitude dependence of folding rates on protein size for a set of 93 proteins, at the same time it is consistent with the near independence of folding equilibrium constant on size. This model gives estimates of folding rates of proteomes, leading to a median folding time in Escherichia coli of about 5 s.
Collapse
Affiliation(s)
- Geoffrey C Rollins
- Department of Biochemistry and Biophysics, University of California , San Francisco, California 94143, United States
| | | |
Collapse
|
20
|
Naganathan AN, Muñoz V. Thermodynamics of Downhill Folding: Multi-Probe Analysis of PDD, a Protein that Folds Over a Marginal Free Energy Barrier. J Phys Chem B 2014; 118:8982-94. [DOI: 10.1021/jp504261g] [Citation(s) in RCA: 32] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Affiliation(s)
- Athi N. Naganathan
- Department
of Chemistry and Biochemistry, University of Maryland, College Park, Maryland 20742, United States
- Department of Biotechnology, Bhupat & Jyoti Mehta School of Biosciences, Indian Institute of Technology Madras, Chennai 600036, India
| | - Victor Muñoz
- Department
of Chemistry and Biochemistry, University of Maryland, College Park, Maryland 20742, United States
- Centro Nacional
de Biotecnología, Consejo Superior de Investigaciones Científicas, 28049 Madrid, Spain
| |
Collapse
|
21
|
Narayan A, Naganathan AN. Evidence for the sequential folding mechanism in RNase H from an ensemble-based model. J Phys Chem B 2014; 118:5050-8. [PMID: 24762044 DOI: 10.1021/jp500934f] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
The number of distinct protein folding pathways starting from an unfolded ensemble, and hence, the folding mechanism is an intricate function of protein size, sequence complexity, and stability conditions. This has traditionally been a contentious issue particularly because of the ensemble nature of conventional experiments that can mask the complexity of the underlying folding landscape. Recent hydrogen-exchange experiments combined with mass spectrometry (HX-MS) reveal that the folding of RNase H proceeds in a hierarchical fashion with distinct intermediates and following a single discrete path. In our current work, we provide computational evidence for this unique folding mechanism by employing a structure-based statistical mechanical model. Upon calibrating the energetic terms of the model with equilibrium measurements, we predict multiple intermediate states in the folding of RNase H that closely resemble experimental observations. Remarkably, a simplified landscape representation adequately captures the folding complexity and predicts the possibility of a well-defined sequence of folding events. We supplement the statistical model study with both explicit solvent molecular simulations of the helical units and electrostatic calculations to provide structural and energetic insights into the early and late stages of RNase H folding that hint at the frustrated nature of its folding landscape.
Collapse
Affiliation(s)
- Abhishek Narayan
- Department of Biotechnology, Bhupat & Jyoti Mehta School of Biosciences, Indian Institute of Technology Madras , Chennai 600036, India
| | | |
Collapse
|
22
|
A disorder-induced domino-like destabilization mechanism governs the folding and functional dynamics of the repeat protein IκBα. PLoS Comput Biol 2013; 9:e1003403. [PMID: 24367251 PMCID: PMC3868533 DOI: 10.1371/journal.pcbi.1003403] [Citation(s) in RCA: 23] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/18/2013] [Accepted: 11/07/2013] [Indexed: 11/19/2022] Open
Abstract
The stability of the repeat protein IκBα, a transcriptional inhibitor in mammalian cells, is critical in the functioning of the NF-κB signaling module implicated in an array of cellular processes, including cell growth, disease, immunity and apoptosis. Structurally, IκBα is complex, with both ordered and disordered regions, thus posing a challenge to the available computational protocols to model its conformational behavior. Here, we introduce a simple procedure to model disorder in systems that undergo binding-induced folding that involves modulation of the contact map guided by equilibrium experimental observables in combination with an Ising-like Wako-Saitô-Muñoz-Eaton model. This one-step procedure alone is able to reproduce a variety of experimental observables, including ensemble thermodynamics (scanning calorimetry, pre-transitions, m-values) and kinetics (roll-over in chevron plot, intermediates and their identity), and is consistent with hydrogen-deuterium exchange measurements. We further capture the intricate distance-dynamics between the domains as measured by single-molecule FRET by combining the model predictions with simple polymer physics arguments. Our results reveal a unique mechanism at work in IκBα folding, wherein disorder in one domain initiates a domino-like effect partially destabilizing neighboring domains, thus highlighting the effect of symmetry-breaking at the level of primary sequences. The offshoot is a multi-state and a dynamic conformational landscape that is populated by increasingly partially folded ensembles upon destabilization. Our results provide, in a straightforward fashion, a rationale to the promiscuous binding and short intracellular half-life of IκBα evolutionarily engineered into it through repeats with variable stabilities and expand the functional repertoire of disordered regions in proteins.
Collapse
|
23
|
Comparing a simple theoretical model for protein folding with all-atom molecular dynamics simulations. Proc Natl Acad Sci U S A 2013; 110:17880-5. [PMID: 24128764 DOI: 10.1073/pnas.1317105110] [Citation(s) in RCA: 83] [Impact Index Per Article: 7.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
Advances in computing have enabled microsecond all-atom molecular dynamics trajectories of protein folding that can be used to compare with and test critical assumptions of theoretical models. We show that recent simulations by the Shaw group (10, 11, 14, 15) are consistent with a key assumption of an Ising-like theoretical model that native structure grows in only a few regions of the amino acid sequence as folding progresses. The distribution of mechanisms predicted by simulating the master equation of this native-centric model for the benchmark villin subdomain, with only two adjustable thermodynamic parameters and one temperature-dependent kinetic parameter, is remarkably similar to the distribution in the molecular dynamics trajectories.
Collapse
|
24
|
Tsai M, Yuan J, Yamaki M, Lin C, Lin SH. Molecular Dynamics Insight into the Diverse Thermodynamic Behavior of a Beta‐Hairpin Peptide. J CHIN CHEM SOC-TAIP 2013. [DOI: 10.1002/jccs.201300173] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Affiliation(s)
- Min‐Yeh Tsai
- Department of Applied Chemistry, National Chiao Tung University, Hsinchu, Taiwan 30010, ROC
- Department of Chemistry, National Taiwan University, Taipei, Taiwan 10617, ROC
| | - Jian‐Min Yuan
- Department of Physics, Drexel University, Philadelphia, Pennsylvania 19104, USA
| | - Masahiro Yamaki
- Department of Applied Chemistry, National Chiao Tung University, Hsinchu, Taiwan 30010, ROC
| | - Chih‐Kai Lin
- Department of Applied Chemistry, National Chiao Tung University, Hsinchu, Taiwan 30010, ROC
| | - Sheng Hsien Lin
- Department of Applied Chemistry, National Chiao Tung University, Hsinchu, Taiwan 30010, ROC
| |
Collapse
|
25
|
Naganathan AN. A Rapid, Ensemble and Free Energy Based Method for Engineering Protein Stabilities. J Phys Chem B 2013; 117:4956-64. [PMID: 23541220 DOI: 10.1021/jp401588x] [Citation(s) in RCA: 25] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Athi N. Naganathan
- Department of Biotechnology, Indian Institute of Technology Madras, Chennai 600036, India
| |
Collapse
|
26
|
Naganathan AN. Predictions from an Ising-like Statistical Mechanical Model on the Dynamic and Thermodynamic Effects of Protein Surface Electrostatics. J Chem Theory Comput 2012; 8:4646-56. [PMID: 26605620 DOI: 10.1021/ct300676w] [Citation(s) in RCA: 48] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Abstract
Charged residues on the surface of a protein are known hot-spots for post-translational modification, protein/ligand-binding, and tuning conformational stabilities. Recent experimental evidence points to the fact that surface electrostatics can also modulate thermodynamic barriers and hence folding mechanisms. To probe for this behavior across different proteins, we develop a novel version of the Wako-Saitô-Muñoz-Eaton (WSME) model in which we include an electrostatic potential term in the energy function while simplifying the treatment of solvation free energy. Both of the energy terms are obtained by quantitatively fitting the model to differential scanning calorimetry (DSC) experiments that carry critical information on the protein partition function. We characterize four sets of structural/functional homologues (HEWL/BLA, CspB, engrailed, α-spectrin) either by fitting the experimental data of a single domain in the homologous set and predicting the conformational behavior of the rest with the same set of parameters or by performing semiblind predictions. The model with the added electrostatic term is able to successfully reproduce the order of thermodynamic stabilities and relaxation rates of most of the homologues. In parallel, we predict diverse conformational features including a wide range of thermodynamic barriers (∼9-40 kJ/mol), broad native ensembles in helical proteins, structured unfolded states and intermediates, rugged folding landscapes, and further provide an independent protein-specific estimate of the folding speed limit at 298 K (1/(7-300 μs)). Our results are evidence that protein surface electrostatics can be tailored to not only engineer stabilities but also folding mechanisms and the ruggedness of the underlying landscape.
Collapse
Affiliation(s)
- Athi N Naganathan
- Department of Biotechnology, Indian Institute of Technology Madras , Chennai-600036, India
| |
Collapse
|
27
|
Tsai MY, Yuan JM, Teranishi Y, Lin SH. Thermodynamics of protein folding using a modified Wako-Saitô-Muñoz-Eaton model. J Biol Phys 2012; 38:543-71. [PMID: 24615219 PMCID: PMC3473134 DOI: 10.1007/s10867-012-9271-y] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/26/2011] [Accepted: 05/07/2012] [Indexed: 10/28/2022] Open
Abstract
Herein, we propose a modified version of the Wako-Saitô-Muñoz-Eaton (WSME) model. The proposed model introduces an empirical temperature parameter for the hypothetical structural units (i.e., foldons) in proteins to include site-dependent thermodynamic behavior. The thermodynamics for both our proposed model and the original WSME model were investigated. For a system with beta-hairpin topology, a mathematical treatment (contact-pair treatment) to facilitate the calculation of its partition function was developed. The results show that the proposed model provides better insight into the site-dependent thermodynamic behavior of the system, compared with the original WSME model. From this site-dependent point of view, the relationship between probe-dependent experimental results and model's thermodynamic predictions can be explained. The model allows for suggesting a general principle to identify foldon behavior. We also find that the backbone hydrogen bonds may play a role of structural constraints in modulating the cooperative system. Thus, our study may contribute to the understanding of the fundamental principles for the thermodynamics of protein folding.
Collapse
Affiliation(s)
- Min-Yeh Tsai
- National Chiao Tung University, 1001 Ta Hsuen Road, Hsinchu, Taiwan, Republic of China,
| | | | | | | |
Collapse
|
28
|
Aioanei D, Brucale M, Tessari I, Bubacco L, Samorì B. Worm-like Ising model for protein mechanical unfolding under the effect of osmolytes. Biophys J 2012; 102:342-50. [PMID: 22339871 PMCID: PMC3260691 DOI: 10.1016/j.bpj.2011.12.007] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/23/2011] [Revised: 10/12/2011] [Accepted: 12/06/2011] [Indexed: 10/14/2022] Open
Abstract
We show via single-molecule mechanical unfolding experiments that the osmolyte glycerol stabilizes the native state of the human cardiac I27 titin module against unfolding without shifting its unfolding transition state on the mechanical reaction coordinate. Taken together with similar findings on the immunoglobulin-binding domain of streptococcal protein G (GB1), these experimental results suggest that osmolytes act on proteins through a common mechanism that does not entail a shift of their unfolding transition state. We investigate the above common mechanism via an Ising-like model for protein mechanical unfolding that adds worm-like-chain behavior to a recent generalization of the Wako-Saitô-Muñoz-Eaton model with support for group-transfer free energies. The thermodynamics of the model are exactly solvable, while protein kinetics under mechanical tension can be simulated via Monte Carlo algorithms. Notably, our force-clamp and velocity-clamp simulations exhibit no shift in the position of the unfolding transition state of GB1 and I27 under the effect of various osmolytes. The excellent agreement between experiment and simulation strongly suggests that osmolytes do not assume a structural role at the mechanical unfolding transition state of proteins, acting instead by adjusting the solvent quality for the protein chain analyte.
Collapse
Affiliation(s)
- Daniel Aioanei
- Department of Biochemistry G. Moruzzi, University of Bologna, Bologna, Italy.
| | | | | | | | | |
Collapse
|
29
|
Sborgi L, Verma A, Muñoz V, de Alba E. Revisiting the NMR structure of the ultrafast downhill folding protein gpW from bacteriophage λ. PLoS One 2011; 6:e26409. [PMID: 22087227 PMCID: PMC3208555 DOI: 10.1371/journal.pone.0026409] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/20/2011] [Accepted: 09/26/2011] [Indexed: 11/18/2022] Open
Abstract
GpW is a 68-residue protein from bacteriophage λ that participates in virus head morphogenesis. Previous NMR studies revealed a novel α+β fold for this protein. Recent experiments have shown that gpW folds in microseconds by crossing a marginal free energy barrier (i.e., downhill folding). These features make gpW a highly desirable target for further experimental and computational folding studies. As a step in that direction, we have re-determined the high-resolution structure of gpW by multidimensional NMR on a construct that eliminates the purification tags and unstructured C-terminal tail present in the prior study. In contrast to the previous work, we have obtained a full manual assignment and calculated the structure using only unambiguous distance restraints. This new structure confirms the α+β topology, but reveals important differences in tertiary packing. Namely, the two α-helices are rotated along their main axis to form a leucine zipper. The β-hairpin is orthogonal to the helical interface rather than parallel, displaying most tertiary contacts through strand 1. There also are differences in secondary structure: longer and less curved helices and a hairpin that now shows the typical right-hand twist. Molecular dynamics simulations starting from both gpW structures, and calculations with CS-Rosetta, all converge to our gpW structure. This confirms that the original structure has strange tertiary packing and strained secondary structure. A comparison of NMR datasets suggests that the problems were mainly caused by incomplete chemical shift assignments, mistakes in NOE assignment and the inclusion of ambiguous distance restraints during the automated procedure used in the original study. The new gpW corrects these problems, providing the appropriate structural reference for future work. Furthermore, our results are a cautionary tale against the inclusion of ambiguous experimental information in the determination of protein structures.
Collapse
Affiliation(s)
- Lorenzo Sborgi
- Centro de Investigaciones Biológicas, Consejo Superior de Investigaciones Científicas, Madrid, Spain
| | - Abhinav Verma
- Centro de Investigaciones Biológicas, Consejo Superior de Investigaciones Científicas, Madrid, Spain
| | - Victor Muñoz
- Centro de Investigaciones Biológicas, Consejo Superior de Investigaciones Científicas, Madrid, Spain
- Department of Chemistry and Biochemistry, University of Maryland, College Park, Maryland, United States of America
- * E-mail: (EdA); (VM)
| | - Eva de Alba
- Centro de Investigaciones Biológicas, Consejo Superior de Investigaciones Científicas, Madrid, Spain
- * E-mail: (EdA); (VM)
| |
Collapse
|
30
|
Enciso M, Rey A. Improvement of structure-based potentials for protein folding by native and nonnative hydrogen bonds. Biophys J 2011; 101:1474-82. [PMID: 21943429 PMCID: PMC3177075 DOI: 10.1016/j.bpj.2011.08.017] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/13/2011] [Revised: 07/26/2011] [Accepted: 08/11/2011] [Indexed: 12/29/2022] Open
Abstract
Pure Gō models (where every native interaction equally stabilizes the folded state) have widely proved their convenience in the computational investigation of protein folding. However, a chemistry-based description of the real interactions also provides a desirable tune in the analysis of the folding process, and thus some hybrid Gō potentials that combine both aspects have been proposed. Among all the noncovalent interactions that contribute to protein folding, hydrogen bonds are the only ones with a partial covalent character. This feature makes them directional and, thus, more difficult to model as part of the coarse-grained descriptions that are typically employed in Gō models. Thanks to a simplified but rigorous representation of backbone hydrogen bonds that we have recently proposed, we present in this article a combined potential (Gō + backbone hydrogen bond) to study the thermodynamics of protein folding in the frame of very simple simulation models. We show that the explicit inclusion of hydrogen bonds leads to a systematic improvement in the description of protein folding. We discuss a representative set of examples (from two-state folders to downhill proteins, with different types of native structures) that reveal a relevant agreement with experimental data.
Collapse
Affiliation(s)
| | - Antonio Rey
- Departamento de Química Física I, Facultad de Ciencias Químicas, Universidad Complutense, Madrid, Spain
| |
Collapse
|
31
|
Naganathan AN, Perez-Jimenez R, Muñoz V, Sanchez-Ruiz JM. Estimation of protein folding free energy barriers from calorimetric data by multi-model Bayesian analysis. Phys Chem Chem Phys 2011; 13:17064-76. [PMID: 21769353 DOI: 10.1039/c1cp20156e] [Citation(s) in RCA: 30] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The realization that folding free energy barriers can be small enough to result in significant population of the species at the barrier top has sprouted in several methods to estimate folding barriers from equilibrium experiments. Some of these approaches are based on fitting the experimental thermogram measured by differential scanning calorimetry (DSC) to a one-dimensional representation of the folding free-energy surface (FES). Different physical models have been used to represent the FES: (1) a Landau quartic polynomial as a function of the total enthalpy, which acts as an order parameter; (2) the projection onto a structural order parameter (i.e. number of native residues or native contacts) of the free energy of all the conformations generated by Ising-like statistical mechanical models; and (3) mean-field models that define conformational entropy and stabilization energy as functions of a continuous local order parameter. The fundamental question that emerges is how can we obtain robust, model-independent estimates of the thermodynamic folding barrier from the analysis of DSC experiments. Here we address this issue by comparing the performance of various FES models in interpreting the thermogram of a protein with a marginal folding barrier. We chose the small α-helical protein PDD, which folds-unfolds in microseconds crossing a free energy barrier previously estimated as ~1 RT. The fits of the PDD thermogram to the various models and assumptions produce FES with a consistently small free energy barrier separating the folded and unfolded ensembles. However, the fits vary in quality as well as in the estimated barrier. Applying Bayesian probabilistic analysis we rank the fit performance using a statistically rigorous criterion that leads to a global estimate of the folding barrier and its precision, which for PDD is 1.3 ± 0.4 kJ mol(-1). This result confirms that PDD folds over a minor barrier consistent with the downhill folding regime. We have further validated the multi-model Bayesian approach through the analysis of two additional protein systems: gpW, a midsize single-domain with α + β topology that also folds in microseconds and has been previously catalogued as a downhill folder, and α-spectrin SH3, a domain of similar size but with a β-barrel fold, slow-folding kinetics and two-state-like thermodynamics. From a general viewpoint, the Bayesian analysis developed here results in a statistically robust, virtually model-independent, method to estimate the thermodynamic free-energy barriers to protein folding from DSC thermograms. Our method appears to be sufficiently accurate to consistently detect small differences in the barrier height, and thus opens up the possibility of characterizing experimentally the changes in thermodynamic folding barriers induced by single-point mutations on proteins within the downhill regime.
Collapse
Affiliation(s)
- Athi N Naganathan
- Centro de Investigaciones Biológicas, Consejo Superior de Investigaciones Científicas (CSIC), Ramiro de Maeztu 9, Madrid 28040, Spain
| | | | | | | |
Collapse
|