1
|
Especial JNC, Faísca PFN. Effects of sequence-dependent non-native interactions in equilibrium and kinetic folding properties of knotted proteins. J Chem Phys 2023; 159:065101. [PMID: 37551809 DOI: 10.1063/5.0160886] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/06/2023] [Accepted: 07/24/2023] [Indexed: 08/09/2023] Open
Abstract
Determining the role of non-native interactions in folding dynamics, kinetics, and mechanisms is a classic problem in protein folding. More recently, this question has witnessed a renewed interest in light of the hypothesis that knotted proteins require the assistance of non-native interactions to fold efficiently. Here, we conduct extensive equilibrium and kinetic Monte Carlo simulations of a simple off-lattice C-alpha model to explore the role of non-native interactions in the thermodynamics and kinetics of three proteins embedding a trefoil knot in their native structure. We find that equilibrium knotted conformations are stabilized by non-native interactions that are non-local, and proximal to native ones, thus enhancing them. Additionally, non-native interactions increase the knotting frequency at high temperatures, and in partially folded conformations below the transition temperatures. Although non-native interactions clearly enhance the efficiency of transition from an unfolded conformation to a partially folded knotted one, they are not required to efficiently fold a knotted protein. Indeed, a native-centric interaction potential drives the most efficient folding transition, provided that the simulation temperature is well below the transition temperature of the considered model system.
Collapse
Affiliation(s)
- João N C Especial
- Departamento de Física, Faculdade de Ciências, Ed. C8, Universidade de Lisboa, Campo Grande, Lisboa, Portugal
- BioISI-Biosystems and Integrative Sciences Institute, Faculdade de Ciências, Universidade de Lisboa, Campo Grande, Lisboa, Portugal
| | - Patrícia F N Faísca
- Departamento de Física, Faculdade de Ciências, Ed. C8, Universidade de Lisboa, Campo Grande, Lisboa, Portugal
- BioISI-Biosystems and Integrative Sciences Institute, Faculdade de Ciências, Universidade de Lisboa, Campo Grande, Lisboa, Portugal
| |
Collapse
|
2
|
A Note on the Effects of Linear Topology Preservation in Monte Carlo Simulations of Knotted Proteins. Int J Mol Sci 2022; 23:ijms232213871. [PMID: 36430350 PMCID: PMC9695063 DOI: 10.3390/ijms232213871] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/17/2022] [Revised: 11/07/2022] [Accepted: 11/08/2022] [Indexed: 11/12/2022] Open
Abstract
Monte Carlo simulations are a powerful technique and are widely used in different fields. When applied to complex molecular systems with long chains, such as those in synthetic polymers and proteins, they have the advantage of providing a fast and computationally efficient way to sample equilibrium ensembles and calculate thermodynamic and structural properties under desired conditions. Conformational Monte Carlo techniques employ a move set to perform the transitions in the simulation Markov chain. While accepted conformations must preserve the sequential bonding of the protein chain model and excluded volume among its units, the moves themselves may take the chain across itself. We call this a break in linear topology preservation. In this manuscript, we show, using simple protein models, that there is no difference in equilibrium properties calculated with a move set that preserves linear topology and one that does not. However, for complex structures, such as those of deeply knotted proteins, the preservation of linear topology provides correct equilibrium results but only after long relaxation. In any case, to analyze folding pathways, knotting mechanisms and folding kinetics, the preservation of linear topology may be an unavoidable requirement.
Collapse
|
3
|
Fernández Del Río B, Rey A. Behavior of Proteins under Pressure from Experimental Pressure-Dependent Structures. J Phys Chem B 2021; 125:6179-6191. [PMID: 34100621 PMCID: PMC8478274 DOI: 10.1021/acs.jpcb.1c03313] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Structure-based models are coarse-grained representations of the interactions responsible for the protein folding process. In their simplest form, they use only the native contact map of a given protein to predict the main features of its folding process by computer simulation. Given their limitations, these models are frequently complemented with sequence-dependent contributions or additional information. Specifically, to analyze the effect of pressure on the folding/unfolding transition, special forms of these interaction potentials are employed, which may a priori determine the outcome of the simulations. In this work, we have tried to keep the original simplicity of structure-based models. Therefore, we have used folded structures that have been experimentally determined at different pressures to define native contact maps and thus interactions dependent on pressure. Despite the apparently tiny structural differences induced by pressure, our simulation results provide different thermodynamic and kinetic behaviors, which roughly correspond to experimental observations (when there is a possible comparison) of two proteins used as benchmarks, hen egg-white lysozyme and dihydrofolate reductase. Therefore, this work shows the feasibility of using experimental native structures at different pressures to analyze the global effects of this physical property on the protein folding process.
Collapse
Affiliation(s)
- Beatriz Fernández Del Río
- Departamento de Química Física, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, E-28040 Madrid, Spain
| | - Antonio Rey
- Departamento de Química Física, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, E-28040 Madrid, Spain
| |
Collapse
|
4
|
Perego C, Potestio R. Computational methods in the study of self-entangled proteins: a critical appraisal. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2019; 31:443001. [PMID: 31269476 DOI: 10.1088/1361-648x/ab2f19] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/09/2023]
Abstract
The existence of self-entangled proteins, the native structure of which features a complex topology, unveils puzzling, and thus fascinating, aspects of protein biology and evolution. The discovery that a polypeptide chain can encode the capability to self-entangle in an efficient and reproducible way during folding, has raised many questions, regarding the possible function of these knots, their conservation along evolution, and their role in the folding paradigm. Understanding the function and origin of these entanglements would lead to deep implications in protein science, and this has stimulated the scientific community to investigate self-entangled proteins for decades by now. In this endeavour, advanced experimental techniques are more and more supported by computational approaches, that can provide theoretical guidelines for the interpretation of experimental results, and for the effective design of new experiments. In this review we provide an introduction to the computational study of self-entangled proteins, focusing in particular on the methodological developments related to this research field. A comprehensive collection of techniques is gathered, ranging from knot theory algorithms, that allow detection and classification of protein topology, to Monte Carlo or molecular dynamics strategies, that constitute crucial instruments for investigating thermodynamics and kinetics of this class of proteins.
Collapse
Affiliation(s)
- Claudio Perego
- Max Panck Institute for Polymer Research, Ackermannweg 10, Mainz 55128, Germany
| | | |
Collapse
|
5
|
Especial J, Nunes A, Rey A, Faísca PF. Hydrophobic confinement modulates thermal stability and assists knotting in the folding of tangled proteins. Phys Chem Chem Phys 2019; 21:11764-11775. [PMID: 31114834 DOI: 10.1039/c9cp01701a] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/27/2023]
Abstract
There is growing support for the idea that the in vivo folding process of knotted proteins is assisted by chaperonins, but the mechanism of chaperonin assisted folding remains elusive. Here, we conduct extensive Monte Carlo simulations of lattice and off-lattice models to explore the effects of confinement and hydrophobic intermolecular interactions with the chaperonin cage in the folding and knotting processes. We find that moderate to high protein-cavity interactions (which are likely to be established in the beginning of the chaperonin working cycle) cause an energetic destabilization of the protein that overcomes the entropic stabilization driven by excluded volume, and leads to a decrease of the melting temperature relative to bulk conditions. Moreover, mild-to-moderate hydrophobic interactions with the cavity (which would be established later in the cycle) lead to a significant enhancement of knotting probability in relation to bulk conditions while simultaneously moderating the effect of steric confinement in the enhancement of thermal stability. Our results thus indicate that the chaperonin may be able to assist knotting without simultaneously thermally stabilizing potential misfolded states to a point that would hamper productive folding thus compromising its functional role.
Collapse
Affiliation(s)
- João Especial
- Departamento de Física, Universidade de Lisboa, Campo Grande, Ed. C8, 1749-016 Lisboa, Portugal.
| | | | | | | |
Collapse
|
6
|
Rubio AM, Rey A. Design of a structure-based model for protein folding from flexible conformations. Phys Chem Chem Phys 2019; 21:6544-6552. [DOI: 10.1039/c9cp00168a] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/04/2023]
Abstract
We introduce a coarse-grained, structure-based model for protein folding that considers the flexibility of the native state in the definition of the model interactions.
Collapse
Affiliation(s)
- Ana M. Rubio
- Departamento de Química Física
- Facultad de Ciencias Químicas
- Universidad Complutense de Madrid
- E-28040 Madrid
- Spain
| | - Antonio Rey
- Departamento de Química Física
- Facultad de Ciencias Químicas
- Universidad Complutense de Madrid
- E-28040 Madrid
- Spain
| |
Collapse
|
7
|
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
|
8
|
Exploring the Sequence-based Prediction of Folding Initiation Sites in Proteins. Sci Rep 2017; 7:8826. [PMID: 28821744 PMCID: PMC5562875 DOI: 10.1038/s41598-017-08366-3] [Citation(s) in RCA: 26] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/28/2017] [Accepted: 07/10/2017] [Indexed: 11/23/2022] Open
Abstract
Protein folding is a complex process that can lead to disease when it fails. Especially poorly understood are the very early stages of protein folding, which are likely defined by intrinsic local interactions between amino acids close to each other in the protein sequence. We here present EFoldMine, a method that predicts, from the primary amino acid sequence of a protein, which amino acids are likely involved in early folding events. The method is based on early folding data from hydrogen deuterium exchange (HDX) data from NMR pulsed labelling experiments, and uses backbone and sidechain dynamics as well as secondary structure propensities as features. The EFoldMine predictions give insights into the folding process, as illustrated by a qualitative comparison with independent experimental observations. Furthermore, on a quantitative proteome scale, the predicted early folding residues tend to become the residues that interact the most in the folded structure, and they are often residues that display evolutionary covariation. The connection of the EFoldMine predictions with both folding pathway data and the folded protein structure suggests that the initial statistical behavior of the protein chain with respect to local structure formation has a lasting effect on its subsequent states.
Collapse
|
9
|
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
|
10
|
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
|
11
|
Matsuoka M, Sugita M, Kikuchi T. Implication of the cause of differences in 3D structures of proteins with high sequence identity based on analyses of amino acid sequences and 3D structures. BMC Res Notes 2014; 7:654. [PMID: 25231773 PMCID: PMC4180342 DOI: 10.1186/1756-0500-7-654] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/29/2014] [Accepted: 09/05/2014] [Indexed: 11/10/2022] Open
Abstract
BACKGROUND Proteins that share a high sequence homology while exhibiting drastically different 3D structures are investigated in this study. Recently, artificial proteins related to the sequences of the GA and IgG binding GB domains of human serum albumin have been designed. These artificial proteins, referred to as GA and GB, share 98% amino acid sequence identity but exhibit different 3D structures, namely, a 3α bundle versus a 4β + α structure. Discriminating between their 3D structures based on their amino acid sequences is a very difficult problem. In the present work, in addition to using bioinformatics techniques, an analysis based on inter-residue average distance statistics is used to address this problem. RESULTS It was hard to distinguish which structure a given sequence would take only with the results of ordinary analyses like BLAST and conservation analyses. However, in addition to these analyses, with the analysis based on the inter-residue average distance statistics and our sequence tendency analysis, we could infer which part would play an important role in its structural formation. CONCLUSIONS The results suggest possible determinants of the different 3D structures for sequences with high sequence identity. The possibility of discriminating between the 3D structures based on the given sequences is also discussed.
Collapse
Affiliation(s)
| | | | - Takeshi Kikuchi
- Department of Bioinformatics, College of Life Sciences, Ritsumeikan University, 1-1-1 Nojihigashi, Kusatsu, Shiga, Japan.
| |
Collapse
|
12
|
González-Charro V, Rey A. Intermediates in the folding equilibrium of repeat proteins from the TPR family. EUROPEAN BIOPHYSICS JOURNAL: EBJ 2014; 43:433-43. [DOI: 10.1007/s00249-014-0975-8] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/23/2014] [Revised: 06/20/2014] [Accepted: 07/03/2014] [Indexed: 11/29/2022]
|
13
|
Enciso M, Rey A. Sketching protein aggregation with a physics-based toy model. J Chem Phys 2013; 139:115101. [DOI: 10.1063/1.4820793] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
14
|
Sugita M, Kikuchi T. Incorporating into a Cα Go model the effects of geometrical restriction on Cα atoms caused by side chain orientations. Proteins 2013; 81:1434-45. [DOI: 10.1002/prot.24294] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/18/2012] [Revised: 03/13/2013] [Accepted: 03/17/2013] [Indexed: 11/10/2022]
Affiliation(s)
- Masatake Sugita
- Department of Bioinformatics; College of Life Sciences; Ritsumeikan University; Kusatsu; Shiga; Japan
| | - Takeshi Kikuchi
- Department of Bioinformatics; College of Life Sciences; Ritsumeikan University; Kusatsu; Shiga; Japan
| |
Collapse
|
15
|
Perezzan R, Rey A. Simulating protein unfolding under pressure with a coarse-grained model. J Chem Phys 2012; 137:185102. [DOI: 10.1063/1.4765057] [Citation(s) in RCA: 11] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
16
|
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
|
17
|
Chan HS, Zhang Z, Wallin S, Liu Z. Cooperativity, local-nonlocal coupling, and nonnative interactions: principles of protein folding from coarse-grained models. Annu Rev Phys Chem 2011; 62:301-26. [PMID: 21453060 DOI: 10.1146/annurev-physchem-032210-103405] [Citation(s) in RCA: 172] [Impact Index Per Article: 13.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Abstract
Coarse-grained, self-contained polymer models are powerful tools in the study of protein folding. They are also essential to assess predictions from less rigorous theoretical approaches that lack an explicit-chain representation. Here we review advances in coarse-grained modeling of cooperative protein folding, noting in particular that the Levinthal paradox was raised in response to the experimental discovery of two-state-like folding in the late 1960s, rather than to the problem of conformational search per se. Comparisons between theory and experiment indicate a prominent role of desolvation barriers in cooperative folding, which likely emerges generally from a coupling between local conformational preferences and nonlocal packing interactions. Many of these principles have been elucidated by native-centric models, wherein nonnative interactions may be treated perturbatively. We discuss these developments as well as recent applications of coarse-grained chain modeling to knotted proteins and to intrinsically disordered proteins.
Collapse
Affiliation(s)
- Hue Sun Chan
- Department of Biochemistry, University of Toronto, Toronto, Ontario M5S 1A8, Canada.
| | | | | | | |
Collapse
|
18
|
Su JG, Chen WZ, Wang CX. Role of electrostatic interactions for the stability and folding behavior of cold shock protein. Proteins 2010; 78:2157-69. [DOI: 10.1002/prot.22730] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/04/2023]
|
19
|
Wu L, Li WF, Liu F, Zhang J, Wang J, Wang W. Understanding protein folding cooperativity based on topological consideration. J Chem Phys 2009; 131:065105. [PMID: 19691415 DOI: 10.1063/1.3200952] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The folding cooperativity is an important issue of protein folding dynamics. Since the native topology plays a significant role in determining the folding behavior of proteins, we believe that it also has close relationship with the folding cooperativity. In the present work, we perform simulations on proteins Naf-BBL, QNND-BBL, CI2, and SH3 with the Gō model and compare their different folding behaviors. By analyzing the weak cooperative folding of protein Naf-BBL in detail, we found that the folding of Naf-BBL shows relatively weak thermodynamic coupling between residues, and such weak coupling is found mainly between the nonlocal native contacts. This finding complements our understandings on the source of barrierless folding of Naf-BBL and promotes us to analyze the topological origins of the poor thermodynamic coupling of Naf-BBL. Then, we further extend our analysis to other two-state and multistate proteins. Based on the considerations of the thermodynamic coupling and kinetic coupling, we conclude that the fraction of scattered native contacts, the difference in loop entropy of contacts, and the long range relative contact order are the major topological factors that influence the folding cooperativity. The combination of these three tertiary structural features shows significant correlations with the folding types of proteins. Moreover, we also discuss the topological factors related to downhill folding. Finally, the generic role of tertiary structure in determining the folding cooperativity is summarized.
Collapse
Affiliation(s)
- L Wu
- Department of Physics and National Laboratory of Solid State Microstructure, Nanjing University, Nanjing 210093, China
| | | | | | | | | | | |
Collapse
|
20
|
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
|
21
|
Larriva M, Prieto L, Bruscolini P, Rey A. A simple simulation model can reproduce the thermodynamic folding intermediate of apoflavodoxin. Proteins 2009; 78:73-82. [DOI: 10.1002/prot.22521] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
|
22
|
Rey-Stolle MF, Enciso M, Rey A. Topology-based models and NMR structures in protein folding simulations. J Comput Chem 2009; 30:1212-9. [DOI: 10.1002/jcc.21149] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/21/2022]
|
23
|
Prieto L, Rey A. Topology-based potentials and the study of the competition between protein folding and aggregation. J Chem Phys 2009; 130:115101. [DOI: 10.1063/1.3089708] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
24
|
Hills RD, Brooks CL. Insights from coarse-grained Gō models for protein folding and dynamics. Int J Mol Sci 2009; 10:889-905. [PMID: 19399227 PMCID: PMC2672008 DOI: 10.3390/ijms10030889] [Citation(s) in RCA: 204] [Impact Index Per Article: 13.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/22/2009] [Revised: 02/23/2009] [Accepted: 02/26/2009] [Indexed: 12/17/2022] Open
Abstract
Exploring the landscape of large scale conformational changes such as protein folding at atomistic detail poses a considerable computational challenge. Coarse-grained representations of the peptide chain have therefore been developed and over the last decade have proved extremely valuable. These include topology-based Gō models, which constitute a smooth and funnel-like approximation to the folding landscape. We review the many variations of the Gō model that have been employed to yield insight into folding mechanisms. Their success has been interpreted as a consequence of the dominant role of the native topology in folding. The role of local contact density in determining protein dynamics is also discussed and is used to explain the ability of Gō-like models to capture sequence effects in folding and elucidate conformational transitions.
Collapse
Affiliation(s)
- Ronald D. Hills
- Department of Molecular Biology and Kellogg School of Science and Technology, The Scripps Research Institute, 10550 N. Torrey Pines Rd. TPC6 La Jolla, CA 92037, USA
| | - Charles L. Brooks
- Department of Molecular Biology and Kellogg School of Science and Technology, The Scripps Research Institute, 10550 N. Torrey Pines Rd. TPC6 La Jolla, CA 92037, USA
- Department of Chemistry and Biophysics Program, University of Michigan, 930 N. University Ave, Ann Arbor, MI 48109, USA
- Author to whom correspondence should be addressed; E-Mail:
; Tel. +1-734-647-6682; Fax: +1-734-647-1604
| |
Collapse
|
25
|
Li P, Oliva FY, Naganathan AN, Muñoz V. Dynamics of one-state downhill protein folding. Proc Natl Acad Sci U S A 2009; 106:103-8. [PMID: 19118204 PMCID: PMC2629219 DOI: 10.1073/pnas.0802986106] [Citation(s) in RCA: 86] [Impact Index Per Article: 5.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/11/2008] [Indexed: 11/18/2022] Open
Abstract
The small helical protein BBL has been shown to fold and unfold in the absence of a free energy barrier according to a battery of quantitative criteria in equilibrium experiments, including probe-dependent equilibrium unfolding, complex coupling between denaturing agents, characteristic DSC thermogram, gradual melting of secondary structure, and heterogeneous atom-by-atom unfolding behaviors spanning the entire unfolding process. Here, we present the results of nanosecond T-jump experiments probing backbone structure by IR and end-to-end distance by FRET. The folding dynamics observed with these two probes are both exponential with common relaxation times but have large differences in amplitude following their probe-dependent equilibrium unfolding. The quantitative analysis of amplitude and relaxation time data for both probes shows that BBL folding dynamics are fully consistent with the one-state folding scenario and incompatible with alternative models involving one or several barrier crossing events. At 333 K, the relaxation time for BBL is 1.3 micros, in agreement with previous folding speed limit estimates. However, late folding events at room temperature are an order of magnitude slower (20 micros), indicating a relatively rough underlying energy landscape. Our results in BBL expose the dynamic features of one-state folding and chart the intrinsic time-scales for conformational motions along the folding process. Interestingly, the simple self-averaging folding dynamics of BBL are the exact dynamic properties required in molecular rheostats, thus supporting a biological role for one-state folding.
Collapse
Affiliation(s)
- Peng Li
- Department of Chemistry and Biochemistry, and Center for Biomolecular Structure and Organization, University of Maryland, College Park, MD 20742; and
| | - Fabiana Y. Oliva
- Department of Chemistry and Biochemistry, and Center for Biomolecular Structure and Organization, University of Maryland, College Park, MD 20742; and
| | - Athi N. Naganathan
- Department of Chemistry and Biochemistry, and Center for Biomolecular Structure and Organization, University of Maryland, College Park, MD 20742; and
- Department of Protein Science, Centro de Investigaciones Biológicas, Consejo Superior de Investigaciones Científicas, Ramiro de Maeztu 9, Madrid 28040, Spain
| | - Victor Muñoz
- Department of Chemistry and Biochemistry, and Center for Biomolecular Structure and Organization, University of Maryland, College Park, MD 20742; and
- Department of Protein Science, Centro de Investigaciones Biológicas, Consejo Superior de Investigaciones Científicas, Ramiro de Maeztu 9, Madrid 28040, Spain
| |
Collapse
|
26
|
Prieto L, Rey A. Simulations of the protein folding process using topology-based models depend on the experimental structure. J Chem Phys 2008; 129:115101. [DOI: 10.1063/1.2977744] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
27
|
Badasyan A, Liu Z, Chan HS. Probing possible downhill folding: native contact topology likely places a significant constraint on the folding cooperativity of proteins with approximately 40 residues. J Mol Biol 2008; 384:512-30. [PMID: 18823994 DOI: 10.1016/j.jmb.2008.09.023] [Citation(s) in RCA: 35] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/16/2008] [Revised: 09/06/2008] [Accepted: 09/10/2008] [Indexed: 10/21/2022]
Abstract
Experiments point to appreciable variations in folding cooperativity among natural proteins with approximately 40 residues, indicating that the behaviors of these proteins are valuable for delineating the contributing factors to cooperative folding. To explore the role of native topology in a protein's propensity to fold cooperatively and how native topology might constrain the degree of cooperativity achievable by a given set of physical interactions, we compared folding/unfolding kinetics simulated using three classes of native-centric C(alpha) chain models with different interaction schemes. The approach was applied to two homologous 45-residue fragments from the peripheral subunit-binding domain family and a 39-residue fragment of the N-terminal domain of ribosomal protein L9. Free-energy profiles as functions of native contact number were computed to assess the heights of thermodynamic barriers to folding. In addition, chevron plots of folding/unfolding rates were constructed as functions of native stability to facilitate comparison with available experimental data. Although common Gō-like models with pairwise Lennard-Jones-type interactions generally fold less cooperatively than real proteins, the rank ordering of cooperativity predicted by these models is consistent with experiment for the proteins investigated, showing increasing folding cooperativity with increasing nonlocality of a protein's native contacts. Models that account for water-expulsion (desolvation) barriers and models with many-body (nonadditive) interactions generally entail higher degrees of folding cooperativity indicated by more linear model chevron plots, but the rank ordering of cooperativity remains unchanged. A robust, experimentally valid rank ordering of model folding cooperativity independent of the multiple native-centric interaction schemes tested here argues that native topology places significant constraints on how cooperatively a protein can fold.
Collapse
Affiliation(s)
- Artem Badasyan
- Department of Biochemistry, Faculty of Medicine, University of Toronto, Toronto, Ontario, Canada M5S 1A8
| | | | | |
Collapse
|