1
|
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
|
2
|
Soler MA, Rey A, Faísca PFN. Steric confinement and enhanced local flexibility assist knotting in simple models of protein folding. Phys Chem Chem Phys 2018; 18:26391-26403. [PMID: 27722468 DOI: 10.1039/c6cp05086g] [Citation(s) in RCA: 22] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Abstract
The chaperonin complex GroEL-GroES is able to accelerate the folding process of knotted proteins considerably. However, the folding mechanism inside the chaperonin cage is elusive. Here we use a combination of lattice and off-lattice Monte Carlo simulations of simple Gō models to study the effect of physical confinement and local flexibility on the folding process of protein model systems embedding a trefoil knot in their native structure. This study predicts that steric confinement plays a specific role in the folding of knotted proteins by increasing the knotting probability for very high degrees of confinement. This effect is observed for protein MJ0366 even above the melting temperature for confinement sizes compatible with the size of the GroEL/GroES chaperonin cage. An enhanced local flexibility produces the same qualitative effects on the folding process. In particular, we observe that knotting probability increases up to 40% in the transition state of protein MJ0366 when flexibility is enhanced. This is underlined by a structural change in the transition state, which becomes devoid of helical content. No relation between the knotting mechanism and flexibility was found in the context of the off-lattice model adopted in this work.
Collapse
Affiliation(s)
- Miguel A Soler
- Dipartimento di Scienze Mediche e Biologiche, Universita' di Udine, Piazzale Kolbe 4, 33100 Udine, Italy
| | - Antonio Rey
- Departamento de Química Física I, Facultad de Ciencias Químicas, Universidad Complutense, 28040 Madrid, Spain.
| | - Patrícia F N Faísca
- Departamento de Física and BioISI - Biosystems and Integrative Sciences Institute, Faculdade de Ciências, Universidade de Lisboa, Campo Grande, Ed. C8, 1749-016 Lisboa, Portugal.
| |
Collapse
|
3
|
Affiliation(s)
- A. Subha Mahadevi
- Centre for Molecular Modelling, CSIR-Indian Institute of Chemical Technology, Tarnaka, Hyderabad, India 500607
| | - G. Narahari Sastry
- Centre for Molecular Modelling, CSIR-Indian Institute of Chemical Technology, Tarnaka, Hyderabad, India 500607
| |
Collapse
|
4
|
Krobath H, Rey A, Faísca PFN. How determinant is N-terminal to C-terminal coupling for protein folding? Phys Chem Chem Phys 2015; 17:3512-24. [DOI: 10.1039/c4cp05178e] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The existence of native interactions between the protein termini is a major determinant of the free energy barrier in a two-state folding transition being therefore a critical modulator of protein folding cooperativity.
Collapse
Affiliation(s)
- Heinrich Krobath
- Centro de Física da Matéria Condensada and Departamento de Física
- Faculdade de Ciências da Universidade de Lisboa
- Portugal
| | - Antonio Rey
- Departamento de Química Física I
- Facultad de Ciencias Químicas
- Universidad Complutense
- Madrid
- Spain
| | - Patrícia F. N. Faísca
- Centro de Física da Matéria Condensada and Departamento de Física
- Faculdade de Ciências da Universidade de Lisboa
- Portugal
| |
Collapse
|
5
|
Lane TJ, Schwantes CR, Beauchamp KA, Pande VS. Probing the origins of two-state folding. J Chem Phys 2014; 139:145104. [PMID: 24116650 DOI: 10.1063/1.4823502] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/26/2023] Open
Abstract
Many protein systems fold in a two-state manner. Random models, however, rarely display two-state kinetics and thus such behavior should not be accepted as a default. While theories for the prevalence of two-state kinetics have been presented, none sufficiently explain the breadth of experimental observations. A model, making minimal assumptions, is introduced that suggests two-state behavior is likely for any system with an overwhelmingly populated native state. We show two-state folding is a natural consequence of such two-state thermodynamics, and is strengthened by increasing the population of the native state. Further, the model exhibits hub-like behavior, with slow interconversions between unfolded states. Despite this, the unfolded state equilibrates quickly relative to the folding time. This apparent paradox is readily understood through this model. Finally, our results compare favorable with measurements of folding rates as a function of chain length and Keq, providing new insight into these relations.
Collapse
Affiliation(s)
- Thomas J Lane
- Department of Chemistry, Stanford University, Stanford, California 94305, USA
| | | | | | | |
Collapse
|
6
|
Maher B, Albrecht AA, Loomes M, Yang XS, Steinhöfel K. A firefly-inspired method for protein structure prediction in lattice models. Biomolecules 2014; 4:56-75. [PMID: 24970205 PMCID: PMC4030990 DOI: 10.3390/biom4010056] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/01/2013] [Revised: 12/17/2013] [Accepted: 12/27/2013] [Indexed: 02/05/2023] Open
Abstract
We introduce a Firefly-inspired algorithmic approach for protein structure prediction over two different lattice models in three-dimensional space. In particular, we consider three-dimensional cubic and three-dimensional face-centred-cubic (FCC) lattices. The underlying energy models are the Hydrophobic-Polar (H-P) model, the Miyazawa–Jernigan (M-J) model and a related matrix model. The implementation of our approach is tested on ten H-P benchmark problems of a length of 48 and ten M-J benchmark problems of a length ranging from 48 until 61. The key complexity parameter we investigate is the total number of objective function evaluations required to achieve the optimum energy values for the H-P model or competitive results in comparison to published values for the M-J model. For H-P instances and cubic lattices, where data for comparison are available, we obtain an average speed-up over eight instances of 2.1, leaving out two extreme values (otherwise, 8.8). For six M-J instances, data for comparison are available for cubic lattices and runs with a population size of 100, where, a priori, the minimum free energy is a termination criterion. The average speed-up over four instances is 1.2 (leaving out two extreme values, otherwise 1.1), which is achieved for a population size of only eight instances. The present study is a test case with initial results for ad hoc parameter settings, with the aim of justifying future research on larger instances within lattice model settings, eventually leading to the ultimate goal of implementations for off-lattice models.
Collapse
Affiliation(s)
- Brian Maher
- Department of Informatics, King's College London, Strand, London WC2R 2LS, UK.
| | - Andreas A Albrecht
- School of Science and Technology, Middlesex University, The Burroughs, London, NW4 4BT, UK.
| | - Martin Loomes
- School of Science and Technology, Middlesex University, The Burroughs, London, NW4 4BT, UK.
| | - Xin-She Yang
- School of Science and Technology, Middlesex University, The Burroughs, London, NW4 4BT, UK.
| | - Kathleen Steinhöfel
- Department of Informatics, King's College London, Strand, London WC2R 2LS, UK.
| |
Collapse
|
7
|
Krobath H, Faísca PFN. Interplay between native topology and non-native interactions in the folding of tethered proteins. Phys Biol 2013; 10:016002. [DOI: 10.1088/1478-3975/10/1/016002] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
|
8
|
Estácio SG, Fernandes CS, Krobath H, Faísca PFN, Shakhnovich EI. Robustness of atomistic Gō models in predicting native-like folding intermediates. J Chem Phys 2012; 137:085102. [DOI: 10.1063/1.4747492] [Citation(s) in RCA: 26] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
9
|
|
10
|
Travasso RDM, Faísca PFN, Rey A. The protein folding transition state: Insights from kinetics and thermodynamics. J Chem Phys 2010; 133:125102. [DOI: 10.1063/1.3485286] [Citation(s) in RCA: 12] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
11
|
Vorov OK, Livesay DR, Jacobs DJ. Helix/coil nucleation: a local response to global demands. Biophys J 2010; 97:3000-9. [PMID: 19948130 DOI: 10.1016/j.bpj.2009.09.013] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/23/2009] [Revised: 09/04/2009] [Accepted: 09/09/2009] [Indexed: 11/18/2022] Open
Abstract
A complete description of protein structure and function must include a proper treatment of mechanisms that lead to cooperativity. The helix/coil transition serves as a simple example of a cooperative folding process, commonly described by a nucleation-propagation mechanism. The prevalent view is that coil structure must first form a short segment of helix in a localized region despite paying a free energy cost (nucleation). Afterward, helical structure propagates outward from the nucleation site. Both processes entail enthalpy-entropy compensation that derives from the loss in conformational entropy on helix formation with concomitant gain in favorable interactions. Nucleation-propagation models inherently assume that cooperativity arises from a sequential series of local events. An alternative distance constraint model asserts there is a direct link between available degrees of freedom and cooperativity through the nonadditivity in conformational entropy. That is, helix nucleation is a concerted manifestation of rigidity propagating through atomic structure. The link between network rigidity and nonadditivity of conformational entropy is shown in this study by solving the distance constraint model using a simple global constraint counting approximation. Cooperativity arises from competition between excess and deficiency in available degrees of freedom in the coil and helix states respectively.
Collapse
Affiliation(s)
- Oleg K Vorov
- Department of Physics, University of North Carolina at Charlotte, Charlotte, North Carolina, USA
| | | | | |
Collapse
|
12
|
Faísca PFN, Travasso RDM, Charters T, Nunes A, Cieplak M. The folding of knotted proteins: insights from lattice simulations. Phys Biol 2010; 7:16009. [PMID: 20130340 DOI: 10.1088/1478-3975/7/1/016009] [Citation(s) in RCA: 42] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022]
Abstract
We carry out systematic Monte Carlo simulations of Gō lattice proteins to investigate and compare the folding processes of two model proteins whose native structures differ from each other due to the presence of a trefoil knot located near the terminus of one of the protein chains. We show that the folding time of the knotted fold is larger than that of the unknotted protein and that this difference in folding time is particularly striking in the temperature region below the optimal folding temperature. Both proteins display similar folding transition temperatures, which is indicative of similar thermal stabilities. By using the folding probability reaction coordinate as an estimator of folding progression we have found out that the formation of the knot is mainly a late folding event in our shallow knot system.
Collapse
Affiliation(s)
- Patrícia F N Faísca
- Centro de Física da Matéria Condensada, Universidade de Lisboa, Av. Prof. Gama Pinto 2, 1649-003 Lisboa, Portugal.
| | | | | | | | | |
Collapse
|
13
|
Abstract
Empirical or knowledge-based potentials have many applications in structural biology such as the prediction of protein structure, protein-protein, and protein-ligand interactions and in the evaluation of stability for mutant proteins, the assessment of errors in experimentally solved structures, and the design of new proteins. Here, we describe a simple procedure to derive and use pairwise distance-dependent potentials that rely on the definition of effective atomic interactions, which attempt to capture interactions that are more likely to be physically relevant. Based on a difficult benchmark test composed of proteins with different secondary structure composition and representing many different folds, we show that the use of effective atomic interactions significantly improves the performance of potentials at discriminating between native and near-native conformations. We also found that, in agreement with previous reports, the potentials derived from the observed effective atomic interactions in native protein structures contain a larger amount of mutual information. A detailed analysis of the effective energy functions shows that atom connectivity effects, which mostly arise when deriving the potential by the incorporation of those indirect atomic interactions occurring beyond the first atomic shell, are clearly filtered out. The shape of the energy functions for direct atomic interactions representing hydrogen bonding and disulfide and salt bridges formation is almost unaffected when effective interactions are taken into account. On the contrary, the shape of the energy functions for indirect atom interactions (i.e., those describing the interaction between two atoms bound to a direct interacting pair) is clearly different when effective interactions are considered. Effective energy functions for indirect interacting atom pairs are not influenced by the shape or the energy minimum observed for the corresponding direct interacting atom pair. Our results suggest that the dependency between the signals in different energy functions is a key aspect that need to be addressed when empirical energy functions are derived and used, and also highlight the importance of additivity assumptions in the use of potential energy functions.
Collapse
Affiliation(s)
- Evandro Ferrada
- Departamento de Genética Molecular y Microbiología, Facultad de Ciencias Biológicas, Pontificia Universidad Católica de Chile, Alameda 340, Santiago, Chile
| | | |
Collapse
|
14
|
Faísca PFN. The nucleation mechanism of protein folding: a survey of computer simulation studies. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2009; 21:373102. [PMID: 21832332 DOI: 10.1088/0953-8984/21/37/373102] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/31/2023]
Abstract
The nucleation mechanism of protein folding, originally proposed by Baldwin in the early 1970s, was firstly observed by Shakhnovich and co-workers two decades later in the context of Monte Carlo simulations of a simple lattice model. At about the same time the extensive use of φ-value analysis provided the first experimental evidence that the folding of Chymotrypsin-inhibitor 2, a small single-domain protein, which folds with two-state kinetics, is also driven by a nucleation mechanism. Since then, the nucleation mechanism is generally considered the most common form of folding mechanism amongst two-state proteins. However, recent experimental data has put forward the idea that this may not necessarily be so, since the accuracy of the experimentally determined φ values, which are used to identify the critical (i.e. nucleating) residues, is typically poor. Here, we provide a survey of in silico results on the nucleation mechanism, ranging from simple lattice Monte Carlo to more sophisticated off-lattice molecular dynamics simulations, and discuss them in light of experimental data.
Collapse
Affiliation(s)
- Patrícia F N Faísca
- Centro de Física Teórica e Computacional, Universidade de Lisboa, Avenida Professor Gama Pinto 2, 1649-003 Lisboa, Portugal
| |
Collapse
|
15
|
Kapsokalivas L, Gan X, Albrecht AA, Steinhöfel K. Population-based local search for protein folding simulation in the MJ energy model and cubic lattices. Comput Biol Chem 2009; 33:283-94. [PMID: 19647489 DOI: 10.1016/j.compbiolchem.2009.06.006] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/06/2009] [Accepted: 06/17/2009] [Indexed: 10/20/2022]
Abstract
We present experimental results on benchmark problems in 3D cubic lattice structures with the Miyazawa-Jernigan energy function for two local search procedures that utilise the pull-move set: (i) population-based local search (PLS) that traverses the energy landscape with greedy steps towards (potential) local minima followed by upward steps up to a certain level of the objective function; (ii) simulated annealing with a logarithmic cooling schedule (LSA). The parameter settings for PLS are derived from short LSA-runs executed in pre-processing and the procedure utilises tabu lists generated for each member of the population. In terms of the total number of energy function evaluations both methods perform equally well, however, PLS has the potential of being parallelised with an expected speed-up in the region of the population size. Furthermore, both methods require a significant smaller number of function evaluations when compared to Monte Carlo simulations with kink-jump moves.
Collapse
Affiliation(s)
- L Kapsokalivas
- King's College London, Department of Computer Science, London WC2R 2LS, England, United Kingdom
| | | | | | | |
Collapse
|
16
|
Barrick D. What have we learned from the studies of two-state folders, and what are the unanswered questions about two-state protein folding? Phys Biol 2009; 6:015001. [PMID: 19208936 DOI: 10.1088/1478-3975/6/1/015001] [Citation(s) in RCA: 39] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
Abstract
Small proteins with globular structures often fold by simple all-or-none mechanisms, both in an equilibrium and a kinetic sense, despite the very large number of partly folded conformations available. This type of 'two-state' folding will be discussed in terms of experimental tests, underlying molecular mechanisms, and limits to two-state behavior. Factors that appear to be important for two-state folding include topology (sequence distance of contacts in the native structure), molecular cooperativity and local energy distribution. Because their local stability distributions and cooperativities can be dissected and analyzed separately from topological features, recent studies of the folding of symmetric proteins will be discussed as a means to better understand the origins of two-state folding.
Collapse
Affiliation(s)
- Doug Barrick
- T C Department of Biophysics, The Johns Hopkins University, 3400 N Charles St, Baltimore, MD 21218, USA.
| |
Collapse
|
17
|
Sattin BD, Zhao W, Travers K, Chu S, Herschlag D. Direct measurement of tertiary contact cooperativity in RNA folding. J Am Chem Soc 2008; 130:6085-7. [PMID: 18429611 DOI: 10.1021/ja800919q] [Citation(s) in RCA: 56] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
All structured biological macromolecules must overcome the thermodynamic folding problem to populate a unique functional state among a vast ensemble of unfolded and alternate conformations. The exploration of cooperativity in protein folding has helped reveal and distinguish the underlying mechanistic solutions to this folding problem. Analogous dissections of RNA tertiary stability remain elusive, however, despite the central biological importance of folded RNA molecules and the potential to reveal fundamental properties of structured macromolecules via comparisons of protein and RNA folding. We report a direct quantitative measure of tertiary contact cooperativity in a folded RNA. We precisely measured the stability of an independently folding P4-P6 domain from the Tetrahymena thermophila group I intron by single molecule fluorescence resonance energy transfer (smFRET). Using wild-type and mutant RNAs, we found that cooperativity between the two tertiary contacts enhances P4-P6 stability by 3.2 +/- 0.2 kcal/mol.
Collapse
Affiliation(s)
- Bernie D Sattin
- Department of Biochemistry, Stanford University, Stanford, California 94305, USA
| | | | | | | | | |
Collapse
|
18
|
Travasso RDM, Telo da Gama MM, Faísca PFN. Pathways to folding, nucleation events, and native geometry. J Chem Phys 2007; 127:145106. [DOI: 10.1063/1.2777150] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
19
|
Street TO, Bradley CM, Barrick D. Predicting coupling limits from an experimentally determined energy landscape. Proc Natl Acad Sci U S A 2007; 104:4907-12. [PMID: 17360387 PMCID: PMC1829238 DOI: 10.1073/pnas.0608756104] [Citation(s) in RCA: 33] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
Repeat proteins are composed of tandem structural modules in which close contacts do not extend beyond adjacent repeats. Despite the local nature of these close contacts, repeat proteins often unfold as a single, highly coupled unit. Previous studies on the Notch ankyrin domain suggest that this lack of equilibrium unfolding intermediates results both from stabilizing interfaces between each repeat and from a roughly uniform distribution of stability across the folding energy landscape. To investigate this idea, we have generated 15 variants of the Notch ankyrin domain with single and multiple destabilizing substitutions that make the energy landscape uneven. By applying a free energy additivity analysis to these variants, we quantified the destabilization threshold over which repeats 6 and 7 decouple from repeats 1-5. The free energy coupling limit suggested by this additivity analysis ( approximately 4 kcal/mol) is also reflected in m-value analysis and in differences among equilibrium unfolding transitions as monitored by CD versus fluorescence for all 15 variants. All of these observations are quantitatively predicted by analyzing the response of the experimentally determined energy landscape to increasing unevenness. These results highlight the importance of a uniform distribution of local stability in achieving cooperative unfolding.
Collapse
Affiliation(s)
- Timothy O. Street
- T. C. Jenkins Department of Biophysics, The Johns Hopkins University, 3400 North Charles Street, Baltimore, MD 21218
| | - Christina M. Bradley
- T. C. Jenkins Department of Biophysics, The Johns Hopkins University, 3400 North Charles Street, Baltimore, MD 21218
| | - Doug Barrick
- T. C. Jenkins Department of Biophysics, The Johns Hopkins University, 3400 North Charles Street, Baltimore, MD 21218
- To whom correspondence should be addressed. E-mail:
| |
Collapse
|