1
|
Mu ZC, Tan YL, Liu J, Zhang BG, Shi YZ. Computational Modeling of DNA 3D Structures: From Dynamics and Mechanics to Folding. Molecules 2023; 28:4833. [PMID: 37375388 DOI: 10.3390/molecules28124833] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/09/2023] [Revised: 06/11/2023] [Accepted: 06/14/2023] [Indexed: 06/29/2023] Open
Abstract
DNA carries the genetic information required for the synthesis of RNA and proteins and plays an important role in many processes of biological development. Understanding the three-dimensional (3D) structures and dynamics of DNA is crucial for understanding their biological functions and guiding the development of novel materials. In this review, we discuss the recent advancements in computer methods for studying DNA 3D structures. This includes molecular dynamics simulations to analyze DNA dynamics, flexibility, and ion binding. We also explore various coarse-grained models used for DNA structure prediction or folding, along with fragment assembly methods for constructing DNA 3D structures. Furthermore, we also discuss the advantages and disadvantages of these methods and highlight their differences.
Collapse
Affiliation(s)
- Zi-Chun Mu
- Research Center of Nonlinear Science, School of Mathematical & Physical Sciences, Wuhan Textile University, Wuhan 430073, China
- School of Computer Science and Artificial Intelligence, Wuhan Textile University, Wuhan 430073, China
| | - Ya-Lan Tan
- Research Center of Nonlinear Science, School of Mathematical & Physical Sciences, Wuhan Textile University, Wuhan 430073, China
| | - Jie Liu
- Research Center of Nonlinear Science, School of Mathematical & Physical Sciences, Wuhan Textile University, Wuhan 430073, China
| | - Ben-Gong Zhang
- Research Center of Nonlinear Science, School of Mathematical & Physical Sciences, Wuhan Textile University, Wuhan 430073, China
| | - Ya-Zhou Shi
- Research Center of Nonlinear Science, School of Mathematical & Physical Sciences, Wuhan Textile University, Wuhan 430073, China
| |
Collapse
|
2
|
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
|
3
|
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
|
4
|
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
|
5
|
Cheng Q, Joung I, Lee J, Kuwajima K, Lee J. Exploring the Folding Mechanism of Small Proteins GB1 and LB1. J Chem Theory Comput 2019; 15:3432-3449. [DOI: 10.1021/acs.jctc.8b01163] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Qianyi Cheng
- Department of Chemistry, University of Memphis, Memphis, Tennessee 38152, United States
- School of Computational Sciences, Korea Institute for Advanced Study, Seoul 02455, South Korea
| | - InSuk Joung
- Department of Chemistry, Kangwon National University, Chuncheon 24341, South Korea
- School of Computational Sciences, Korea Institute for Advanced Study, Seoul 02455, South Korea
| | - Juyong Lee
- Department of Chemistry, Kangwon National University, Chuncheon 24341, South Korea
| | - Kunihiro Kuwajima
- Department of Physics, University of Tokyo, Tokyo 113-0033, Japan
- School of Computational Sciences, Korea Institute for Advanced Study, Seoul 02455, South Korea
| | - Jooyoung Lee
- Center for In Silico Protein Science, Korea Institute for Advanced Study, Seoul 02455, South Korea
- School of Computational Sciences, Korea Institute for Advanced Study, Seoul 02455, South Korea
| |
Collapse
|
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
|
Faísca PF. Knotted proteins: A tangled tale of Structural Biology. Comput Struct Biotechnol J 2015; 13:459-68. [PMID: 26380658 PMCID: PMC4556803 DOI: 10.1016/j.csbj.2015.08.003] [Citation(s) in RCA: 58] [Impact Index Per Article: 6.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/28/2015] [Revised: 07/31/2015] [Accepted: 08/07/2015] [Indexed: 01/19/2023] Open
Abstract
Knotted proteins have their native structures arranged in the form of an open knot. In the last ten years researchers have been making significant efforts to reveal their folding mechanism and understand which functional advantage(s) knots convey to their carriers. Molecular simulations have been playing a fundamental role in this endeavor, and early computational predictions about the knotting mechanism have just been confirmed in wet lab experiments. Here we review a collection of simulation results that allow outlining the current status of the field of knotted proteins, and discuss directions for future research.
Collapse
|
9
|
Liu L, Wade RC, Heermann DW. A multiscale approach to simulating the conformational properties of unbound multi-C₂H₂ zinc finger proteins. Proteins 2015; 83:1604-15. [PMID: 26062035 DOI: 10.1002/prot.24845] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/16/2015] [Revised: 05/18/2015] [Accepted: 06/03/2015] [Indexed: 12/25/2022]
Abstract
The conformational properties of unbound multi-Cys2 His2 (mC2H2) zinc finger proteins, in which zinc finger domains are connected by flexible linkers, are studied by a multiscale approach. Three methods on different length scales are utilized. First, atomic detail molecular dynamics simulations of one zinc finger and its adjacent flexible linker confirmed that the zinc finger is more rigid than the flexible linker. Second, the end-to-end distance distributions of mC2H2 zinc finger proteins are computed using an efficient atomistic pivoting algorithm, which only takes excluded volume interactions into consideration. The end-to-end distance distribution gradually changes its profile, from left-tailed to right-tailed, as the number of zinc fingers increases. This is explained by using a worm-like chain model. For proteins of a few zinc fingers, an effective bending constraint favors an extended conformation. Only for proteins containing more than nine zinc fingers, is a somewhat compacted conformation preferred. Third, a mesoscale model is modified to study both the local and the global conformational properties of multi-C2H2 zinc finger proteins. Simulations of the CCCTC-binding factor (CTCF), an important mC2H2 zinc finger protein for genome spatial organization, are presented.
Collapse
Affiliation(s)
- Lei Liu
- Institute for Theoretical Physics Heidelberg University, Heidelberg, Germany.,Interdisciplinary Center for Scientific Computing Heidelberg University, Heidelberg, Germany
| | - Rebecca C Wade
- Interdisciplinary Center for Scientific Computing Heidelberg University, Heidelberg, Germany.,Molecular and Cellular Modeling Group, Heidelberg Institute for Theoretical Studies, Heidelberg, Germany.,Center for Molecular Biology (ZMBH), DKFZ-ZMBH Alliance, Heidelberg University, Heidelberg, Germany
| | - Dieter W Heermann
- Institute for Theoretical Physics Heidelberg University, Heidelberg, Germany.,Interdisciplinary Center for Scientific Computing Heidelberg University, Heidelberg, Germany.,Shanghai Institute of Biological Sciences Chinese Academy of Science, Shanghai, People's Republic of China
| |
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
|
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]
|
12
|
Sinner C, Lutz B, John S, Reinartz I, Verma A, Schug A. Simulating Biomolecular Folding and Function by Native-Structure-Based/Go-Type Models. Isr J Chem 2014. [DOI: 10.1002/ijch.201400012] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
|
13
|
Terakawa T, Takada S. RESPAC: Method to Determine Partial Charges in Coarse-Grained Protein Model and Its Application to DNA-Binding Proteins. J Chem Theory Comput 2014; 10:711-21. [DOI: 10.1021/ct4007162] [Citation(s) in RCA: 34] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Tsuyoshi Terakawa
- Department of Biophysics,
Graduate School of Science, Kyoto University, Kitashirakawa Owake-cho, Sakyo, Kyoto, 606-8501, Japan
| | - Shoji Takada
- Department of Biophysics,
Graduate School of Science, Kyoto University, Kitashirakawa Owake-cho, Sakyo, Kyoto, 606-8501, Japan
| |
Collapse
|
14
|
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
|
15
|
Lin C, Culik RM, Gai F. Using VIPT-jump to distinguish between different folding mechanisms: application to BBL and a Trpzip. J Am Chem Soc 2013; 135:7668-73. [PMID: 23642153 PMCID: PMC3706100 DOI: 10.1021/ja401473m] [Citation(s) in RCA: 21] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Protein folding involves a large number of sequential molecular steps or conformational substates. Thus, experimental characterization of the underlying folding energy landscape for any given protein is difficult. Herein, we present a new method that can be used to determine the major characteristics of the folding energy landscape in question, e.g., to distinguish between activated and barrierless downhill folding scenarios. This method is based on the idea that the conformational relaxation kinetics of different folding mechanisms at a given final condition will show different dependences on the initial condition. We show, using both simulation and experiment, that it is possible to differentiate between disparate kinetic folding models by comparing temperature jump (T-jump) relaxation traces obtained with a fixed final temperature and varied initial temperatures, which effectively varies the initial potential (VIP) of the system of interest. We apply this method (hereafter refer to as VIPT-jump) to two model systems, tryptophan zipper (Trpzip)-2c and BBL, and our results show that BBL exhibits characteristics of barrierless downhill folding, whereas Trpzip-2c folding encounters a free energy barrier. In addition, using the T-jump data of BBL we are able to provide, via Langevin dynamics simulations, a realistic estimate of its conformational diffusion coefficient.
Collapse
Affiliation(s)
- Chun–Wei Lin
- Department of Chemistry, University of Pennsylvania, Philadelphia, PA 19104
| | - Robert M. Culik
- Department of Biochemistry and Biophysics, University of Pennsylvania, Philadelphia, PA 19104
| | - Feng Gai
- Department of Chemistry, University of Pennsylvania, Philadelphia, PA 19104
| |
Collapse
|
16
|
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
|
17
|
Enciso M, Rey A. Simple model for the simulation of peptide folding and aggregation with different sequences. J Chem Phys 2012; 136:215103. [DOI: 10.1063/1.4725883] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
18
|
Noel JK, Whitford PC, Onuchic JN. The shadow map: a general contact definition for capturing the dynamics of biomolecular folding and function. J Phys Chem B 2012; 116:8692-702. [PMID: 22536820 DOI: 10.1021/jp300852d] [Citation(s) in RCA: 164] [Impact Index Per Article: 13.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Abstract
Structure-based models (SBMs) are simplified models of the biomolecular dynamics that arise from funneled energy landscapes. We recently introduced an all-atom SBM that explicitly represents the atomic geometry of a biomolecule. While this initial study showed the robustness of the all-atom SBM Hamiltonian to changes in many of the energetic parameters, an important aspect, which has not been explored previously, is the definition of native interactions. In this study, we propose a general definition for generating atomically grained contact maps called "Shadow". The Shadow algorithm initially considers all atoms within a cutoff distance and then, controlled by a screening parameter, discards the occluded contacts. We show that this choice of contact map is not only well behaved for protein folding, since it produces consistently cooperative folding behavior in SBMs but also desirable for exploring the dynamics of macromolecular assemblies since, it distributes energy similarly between RNAs and proteins despite their disparate internal packing. All-atom structure-based models employing Shadow contact maps provide a general framework for exploring the geometrical features of biomolecules, especially the connections between folding and function.
Collapse
Affiliation(s)
- Jeffrey K Noel
- Center for Theoretical Biological Physics and Department of Physics, Rice University, Houston, Texas 77005, United States
| | | | | |
Collapse
|
19
|
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
|
20
|
Chen Y, Ding J. Roles of non-native hydrogen-bonding interaction in helix-coil transition of a single polypeptide as revealed by comparison between Gō-like and non-Gō models. Proteins 2010; 78:2090-100. [PMID: 20455265 DOI: 10.1002/prot.22724] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022]
Abstract
To explore the role of non-native interactions in the helix-coil transition, a detailed comparison between a Gō-like model and a non-Gō model has been performed via lattice Monte Carlo simulations. Only native hydrogen bonding interactions occur in the Gō-like model, and the non-native ones with sequence interval more than 4 is also included into the non-Gō model. Some significant differences between the results from those two models have been found. The non-native hydrogen bonds were found most populated at temperature around the helix-coil transition. The rearrangement of non-native hydrogen bonds into native ones in the formation of alpha-helix leads to the increase of susceptibility of chain conformation, and even two peaks of susceptibility of radius of gyration versus temperature exist in the case of non-Gō model for a non-short peptide, while just a single peak exists in the case of Gō model for a single polypeptide chain with various chain lengths. The non-native hydrogen bonds have complicated the temperature-dependence of Zimm-Bragg nucleation constant. The increase of relative probability of non-native hydrogen bonding for long polypeptide chains leads to non-monotonous chain length effect on the transition temperature.
Collapse
Affiliation(s)
- Yantao Chen
- Key Laboratory of Molecular Engineering of Polymers of Ministry of Education, Department of Macromolecular Science, Laboratory of Advanced Materials, Fudan University, Shanghai 200433, China
| | | |
Collapse
|
21
|
Enciso M, Rey A. A refined hydrogen bond potential for flexible protein models. J Chem Phys 2010; 132:235102. [DOI: 10.1063/1.3436723] [Citation(s) in RCA: 23] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
22
|
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
|
23
|
Arc-repressor dimerization on DNA: folding rate enhancement by colocalization. Biophys J 2009; 96:4212-20. [PMID: 19450491 DOI: 10.1016/j.bpj.2009.01.057] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/25/2008] [Revised: 01/13/2009] [Accepted: 01/15/2009] [Indexed: 02/04/2023] Open
Abstract
Multimeric proteins are ubiquitous in many cellular processes that require high levels of regulation. Eukaryotic gene expression is often regulated by a mechanism of combinatorial control that involves the binding of dimeric transcription factors to DNA together with the coordinated activity of additional proteins. In this study, we investigated the dimerization of the Arc-repressor on DNA with the aim of achieving microscopic insight into the possible advantages of interacting with DNA as a complex rather than as a monomeric single-domain protein. We used a computational coarse-grained model in which the protein dynamics was governed by native interactions and protein-DNA interactions were dictated by electrostatic forces. Inspired by previous experimental work that showed an enhanced refolding rate for the Arc-repressor in the presence of DNA and other polyanions, we focused on the mechanism and kinetics of the assembly of Arc monomers in the presence of single- (ssDNA) and double-stranded DNA (dsDNA) molecules in a low-salt concentration environment. The electrostatic interactions that attract the protein to the dsDNA were shown to be fundamental in colocalizing the unfolded Arc chains and in accelerating refolding. Arc monomers bind the dsDNA efficiently and nonspecifically, and search for each other via one-dimensional diffusion. The fastest folding of Arc is observed for DNA of 30 bp. Longer DNA is significantly less efficient in accelerating the Arc refolding rate, since the two subunits search distinct regions of the one-dimensional DNA and are therefore much less colocalized. The probability that the two unfolded chains will meet on 200 bp DNA is similar to that in the bulk. The colocalization of Arc subunits on ssDNA results in much faster folding compared to that obtained on dsDNA of the same length. Differences in the rate of Arc refolding, cooperativity, and the structure of its transition state ensemble introduced by ssDNA and dsDNA molecules demonstrate the important role of colocalization in biological self-assembly processes.
Collapse
|
24
|
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]
|
25
|
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]
|
26
|
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
|
27
|
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
|
28
|
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
|
29
|
Kim J, Keyes T. Influence of Go-Like Interactions on Global Shapes of Energy Landscapes in β-Barrel Forming Model Proteins: Inherent Structure Analysis and Statistical Temperature Molecular Dynamics Simulation. J Phys Chem B 2007; 112:954-66. [DOI: 10.1021/jp072872u] [Citation(s) in RCA: 14] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Jaegil Kim
- Department of Chemistry, Boston University, Boston, Massachusetts 02215
| | - Thomas Keyes
- Department of Chemistry, Boston University, Boston, Massachusetts 02215
| |
Collapse
|
30
|
Prieto L, Rey A. Influence of the native topology on the folding barrier for small proteins. J Chem Phys 2007; 127:175101. [DOI: 10.1063/1.2780154] [Citation(s) in RCA: 27] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
31
|
Kmiecik S, Kolinski A. Folding pathway of the b1 domain of protein G explored by multiscale modeling. Biophys J 2007; 94:726-36. [PMID: 17890394 PMCID: PMC2186257 DOI: 10.1529/biophysj.107.116095] [Citation(s) in RCA: 84] [Impact Index Per Article: 4.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
The understanding of the folding mechanisms of single-domain proteins is an essential step in the understanding of protein folding in general. Recently, we developed a mesoscopic CA-CB side-chain protein model, which was successfully applied in protein structure prediction, studies of protein thermodynamics, and modeling of protein complexes. In this research, this model is employed in a detailed characterization of the folding process of a simple globular protein, the B1 domain of IgG-binding protein G (GB1). There is a vast body of experimental facts and theoretical findings for this protein. Performing unbiased, ab initio simulations, we demonstrated that the GB1 folding proceeds via the formation of an extended folding nucleus, followed by slow structure fine-tuning. Remarkably, a subset of native interactions drives the folding from the very beginning. The emerging comprehensive picture of GB1 folding perfectly matches and extends the previous experimental and theoretical studies.
Collapse
Affiliation(s)
| | - Andrzej Kolinski
- Address reprint requests to Andrzej Kolinski, Faculty of Chemistry, University of Warsaw, L. Pasteura 1, 02-093 Warsaw, Poland. Tel.: 48-022-8220211 ext. 320; Fax: 48-022 820221.
| |
Collapse
|
32
|
Prieto L, Rey A. Influence of the chain stiffness on the thermodynamics of a Gō-type model for protein folding. J Chem Phys 2007; 126:165103. [PMID: 17477635 DOI: 10.1063/1.2727465] [Citation(s) in RCA: 12] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The relative importance of local and long range interactions in the characteristics of the protein folding process has long been a matter of controversy. Computer simulations based on Gō-type models have been widely used to study this topic, but without much agreement on which type of interactions is more relevant for the foldability of a protein. In this work, the authors also employ a topology-based potential and simulation model to analyze the influence of local and long range interactions on the thermodynamics of the folding transition. The former are mainly used to control the degree of flexibility (or stiffness) of the chain, mostly appreciable in the unfolded (noncompact) state. Our results show the different effects that local and nonlocal interactions have on the entropy and the energy of the system. This implies that a balance between both types of interactions is required, so that a free energy barrier exists between the native and the denatured states. The variations in the contribution of both types of interactions have also a direct effect on the stability of the chain conformations, including the possible appearance of thermodynamic folding intermediates.
Collapse
Affiliation(s)
- Lidia Prieto
- Departamento de Química Física I, Facultad de Ciencias Químicas, Universidad Complutense, E-28040 Madrid, Spain
| | | |
Collapse
|
33
|
Bruscolini P, Pelizzola A, Zamparo M. Downhill versus two-state protein folding in a statistical mechanical model. J Chem Phys 2007; 126:215103. [PMID: 17567220 DOI: 10.1063/1.2738473] [Citation(s) in RCA: 28] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The authors address the problem of downhill protein folding in the framework of a simple statistical mechanical model, which allows an exact solution for the equilibrium and a semianalytical treatment of the kinetics. Focusing on protein 1BBL, a candidate for downhill folding behavior, and comparing it to the WW domain of protein PIN1, a two-state folder of comparable size, the authors show that there are qualitative differences in both the equilibrium and kinetic properties of the two molecules. However, the barrierless scenario which would be expected if 1BBL were a true downhill folder is observed only at low enough temperature.
Collapse
Affiliation(s)
- Pierpaolo Bruscolini
- Instituto BIFI, Universidad de Zaragoza, Calle Corona de Aragón 42, 50009 Zaragoza, Spain.
| | | | | |
Collapse
|