1
|
Farías-Rico JA, Ruud Selin F, Myronidi I, Frühauf M, von Heijne G. Effects of protein size, thermodynamic stability, and net charge on cotranslational folding on the ribosome. Proc Natl Acad Sci U S A 2018; 115:E9280-E9287. [PMID: 30224455 PMCID: PMC6176590 DOI: 10.1073/pnas.1812756115] [Citation(s) in RCA: 47] [Impact Index Per Article: 7.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/23/2023] Open
Abstract
During the last five decades, studies of protein folding in dilute buffer solutions have produced a rich picture of this complex process. In the cell, however, proteins can start to fold while still attached to the ribosome (cotranslational folding) and it is not yet clear how the ribosome affects the folding of protein domains of different sizes, thermodynamic stabilities, and net charges. Here, by using arrest peptides as force sensors and on-ribosome pulse proteolysis, we provide a comprehensive picture of how the distance from the peptidyl transferase center in the ribosome at which proteins fold correlates with protein size. Moreover, an analysis of a large collection of mutants of the Escherichia coli ribosomal protein S6 shows that the force exerted on the nascent chain by protein folding varies linearly with the thermodynamic stability of the folded state, and that the ribosome environment disfavors folding of domains of high net-negative charge.
Collapse
Affiliation(s)
| | - Frida Ruud Selin
- Department of Biochemistry and Biophysics, Stockholm University, SE-106 91 Stockholm, Sweden
| | - Ioanna Myronidi
- Department of Biochemistry and Biophysics, Stockholm University, SE-106 91 Stockholm, Sweden
| | - Marie Frühauf
- Department of Biochemistry and Biophysics, Stockholm University, SE-106 91 Stockholm, Sweden
| | - Gunnar von Heijne
- Department of Biochemistry and Biophysics, Stockholm University, SE-106 91 Stockholm, Sweden;
- Science for Life Laboratory, Stockholm University, SE-171 21 Solna, Sweden
| |
Collapse
|
2
|
Meuzelaar H, Panman MR, van Dijk CN, Woutersen S. Folding of a Zinc-Finger ββα-Motif Investigated Using Two-Dimensional and Time-Resolved Vibrational Spectroscopy. J Phys Chem B 2016; 120:11151-11158. [PMID: 27723346 DOI: 10.1021/acs.jpcb.6b08883] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Small proteins provide good model systems for studying the fundamental forces that control protein folding. Here, we investigate the folding dynamics of the 28-residue zinc-finger mutant FSD-1, which is designed to form a metal-independent folded ββα-motif, and which provides a testing ground for proteins containing a mixed α/β fold. Although the folding of FSD-1 has been actively studied, the folding mechanism remains largely unclear. In particular, it is unclear in what stage of folding the α-helix is formed. To address this issue we investigate the folding mechanism of FSD-1 using a combination of temperature-dependent UV circular dichroism (UV-CD), Fourier transform infrared (FTIR) spectroscopy, two-dimensional infrared (2D-IR) spectroscopy, and temperature-jump (T-jump) transient-IR spectroscopy. Our UV-CD and FTIR data show different thermal melting transitions, indicating multistate folding behavior. Temperature-dependent 2D-IR spectra indicate that the α-helix is the most stable structural element of FSD-1. To investigate the folding/unfolding re-equilibration dynamics of FSD-1, the conformational changes induced by a nanosecond T-jump are probed with transient-IR and transient dispersed-pump-probe (DPP) IR spectroscopy. We observe biexponential T-jump relaxation kinetics (with time constants of 80 ± 13 ns and 1300 ± 100 ns at 322 K), confirming that the folding involves an intermediate state. The IR and dispersed-pump-probe IR spectra associated with the two kinetic components suggest that the folding of FSD-1 involves early formation of the α-helix, followed by the formation of the β-hairpin and hydrophobic contacts.
Collapse
Affiliation(s)
- Heleen Meuzelaar
- Van 't Hoff Institute for Molecular Sciences, University of Amsterdam , Science Park 904, 1098 XH Amsterdam, The Netherlands
| | - Matthijs R Panman
- Van 't Hoff Institute for Molecular Sciences, University of Amsterdam , Science Park 904, 1098 XH Amsterdam, The Netherlands
| | - Chris N van Dijk
- Van 't Hoff Institute for Molecular Sciences, University of Amsterdam , Science Park 904, 1098 XH Amsterdam, The Netherlands
| | - Sander Woutersen
- Van 't Hoff Institute for Molecular Sciences, University of Amsterdam , Science Park 904, 1098 XH Amsterdam, The Netherlands
| |
Collapse
|
3
|
Luna-Martínez OD, Vidal-Limón A, Villalba-Velázquez MI, Sánchez-Alcalá R, Garduño-Juárez R, Uversky VN, Becerril B. Simple approach for ranking structure determining residues. PeerJ 2016; 4:e2136. [PMID: 27366642 PMCID: PMC4924125 DOI: 10.7717/peerj.2136] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/07/2015] [Accepted: 05/25/2016] [Indexed: 12/30/2022] Open
Abstract
Mutating residues has been a common task in order to study structural properties of the protein of interest. Here, we propose and validate a simple method that allows the identification of structural determinants; i.e., residues essential for preservation of the stability of global structure, regardless of the protein topology. This method evaluates all of the residues in a 3D structure of a given globular protein by ranking them according to their connectivity and movement restrictions without topology constraints. Our results matched up with sequence-based predictors that look up for intrinsically disordered segments, suggesting that protein disorder can also be described with the proposed methodology.
Collapse
Affiliation(s)
- Oscar D Luna-Martínez
- Instituto de Biotecnología, Universidad Nacional Autónoma de México , Cuernavaca , Morelos , Mexico
| | - Abraham Vidal-Limón
- División de Biología Molecular, Instituto Potosino de Investigación Científica y Tecnológica , San Luis Potosí , Mexico
| | | | - Rosalba Sánchez-Alcalá
- Instituto de Biotecnología, Universidad Nacional Autónoma de México , Cuernavaca , Morelos , Mexico
| | - Ramón Garduño-Juárez
- Instituto de Ciencias Físicas, Universidad Nacional Autónoma de México , Cuernavaca , Morelos , Mexico
| | - Vladimir N Uversky
- Department of Molecular Medicine and USF Health Byrd Alzheimer's Research Institute, University of South Florida, Tampa, FL, United States; Institute for Biological Instrumentation, Russian Academy of Sciences, Puschino, Moscow Region, Russia; Laboratory of Structural Dynamics, Stability and Folding of Proteins, Russian Academy of Sciences, St. Petersburg, Russia
| | - Baltazar Becerril
- Instituto de Biotecnología, Universidad Nacional Autónoma de México , Cuernavaca , Morelos , Mexico
| |
Collapse
|
4
|
Sokkar P, Choi SM, Rhee YM. Simple Method for Simulating the Mixture of Atomistic and Coarse-Grained Molecular Systems. J Chem Theory Comput 2013; 9:3728-39. [DOI: 10.1021/ct400091a] [Citation(s) in RCA: 25] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Pandian Sokkar
- Center for
Self-assembly and Complexity, Institute for Basic Science (IBS), Pohang 790-784, Korea
| | - Sun Mi Choi
- Center for
Self-assembly and Complexity, Institute for Basic Science (IBS), Pohang 790-784, Korea
- Department of Chemistry, Pohang University of Science and Technology (POSTECH), Pohang 790-784,
Korea
| | - Young Min Rhee
- Center for
Self-assembly and Complexity, Institute for Basic Science (IBS), Pohang 790-784, Korea
- Department of Chemistry, Pohang University of Science and Technology (POSTECH), Pohang 790-784,
Korea
| |
Collapse
|
5
|
Abstract
We have carried out all-atom action-derived molecular dynamics (ADMD) folding simulations of the full-size FSD-1. FSD-1 is a designed mini-protein of 28 residues containing both α and β secondary structure elements. Multiple folding pathways are found for FSD-1, which is consistent with existing computational studies. Hydrophobic collapse is observed first, and then subsequent folding events proceeds by forming either α-helix or β-hairpin. Concurrent formation of the full tertiary structure and the secondary structure elements of α-helix and β-hairpin is observed. The folding pathway of FSD-1 elucidated by ADMD simulations does not follow the scenario of the framework model. ADMD simulations provide significant insights for the general mechanisms of protein folding and conformational changes.
Collapse
Affiliation(s)
- In-Ho Lee
- Korea Research Institute of Standards and Science, Daejeon 305-340, Korea
| | | | | |
Collapse
|
6
|
TSAI YILEN, CHEN HWUNGWEN, LIN TOPP, WANG WEIZHOU, SUN YINGCHIEH. MOLECULAR DYNAMICS SIMULATION OF FOLDING OF A SHORT HELICAL TOXIN PEPTIDE. JOURNAL OF THEORETICAL & COMPUTATIONAL CHEMISTRY 2011. [DOI: 10.1142/s0219633607002964] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
Abstract
A molecular dynamics simulation of the folding of a short helical toxin peptide was carried out. The simulation gave a folding time of ~10 ns, which is longer than typical time of ~1 ns for the formation of 1–2 helical turns. The simulation demonstrates that a helical peptide with disulfide bonds, which may encounter extra steric hindrance compared with the peptide without disulfide bonds, can fold in nanosecond timescale. An analysis shows that this folding time should correspond to the folding time in weak denaturation condition in experiment. Interactions and factors affecting folding pathways are analyzed and discussed.
Collapse
Affiliation(s)
- YI-LEN TSAI
- Department of Chemistry, National Taiwan Normal University, 88, TingChow Road Section 4, Taipei 116, Taiwan
| | - HWUNG-WEN CHEN
- Department of Chemistry, National Taiwan Normal University, 88, TingChow Road Section 4, Taipei 116, Taiwan
| | - TOPP LIN
- Department of Chemistry, National Taiwan Normal University, 88, TingChow Road Section 4, Taipei 116, Taiwan
| | - WEI-ZHOU WANG
- Department of Chemistry, National Taiwan Normal University, 88, TingChow Road Section 4, Taipei 116, Taiwan
| | - YING-CHIEH SUN
- Department of Chemistry, National Taiwan Normal University, 88, TingChow Road Section 4, Taipei 116, Taiwan
| |
Collapse
|
7
|
Stability and folding behavior analysis of zinc-finger using simple models. Int J Mol Sci 2010; 11:4014-34. [PMID: 21152317 PMCID: PMC2996801 DOI: 10.3390/ijms11104014] [Citation(s) in RCA: 19] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/07/2010] [Revised: 10/01/2010] [Accepted: 10/09/2010] [Indexed: 01/03/2023] Open
Abstract
Zinc-fingers play crucial roles in regulating gene expression and mediating protein-protein interactions. In this article, two different proteins (Sp1f2 and FSD-1) are investigated using the Gaussian network model and anisotropy elastic network model. By using these simple coarse-grained methods, we analyze the structural stabilization and establish the unfolding pathway of the two different proteins, in good agreement with related experimental and molecular dynamics simulation data. From the analysis, it is also found that the folding process of the zinc-finger motif is predominated by several factors. Both the zinc ion and C-terminal loop affect the folding pathway of the zinc-finger motif. Knowledge about the stability and folding behavior of zinc-fingers may help in understanding the folding mechanisms of the zinc-finger motif and in designing new zinc-fingers. Meanwhile, these simple coarse-grained analyses can be used as a general and quick method for mechanistic studies of metalloproteins.
Collapse
|
8
|
Feng JA, Kao J, Marshall GR. A second look at mini-protein stability: analysis of FSD-1 using circular dichroism, differential scanning calorimetry, and simulations. Biophys J 2010; 97:2803-10. [PMID: 19917235 DOI: 10.1016/j.bpj.2009.08.046] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/01/2009] [Revised: 08/24/2009] [Accepted: 08/25/2009] [Indexed: 10/20/2022] Open
Abstract
Mini-proteins that contain <50 amino acids often serve as model systems for studying protein folding because their small size makes long timescale simulations possible. However, not all mini-proteins are created equal. The stability and structure of FSD-1, a 28-residue mini-protein that adopted the betabetaalpha zinc-finger motif independent of zinc binding, was investigated using circular dichroism, differential scanning calorimetry, and replica-exchange molecular dynamics. The broad melting transition of FSD-1, similar to that of a helix-to-coil transition, was observed by using circular dichroism, differential scanning calorimetry, and replica-exchange molecular dynamics. The N-terminal beta-hairpin was found to be flexible. The FSD-1 apparent melting temperature of 41 degrees C may be a reflection of the melting of its alpha-helical segment instead of the entire protein. Thus, despite its attractiveness due to small size and purposefully designed helix, sheet, and turn structures, the status of FSD-1 as a model system for studying protein folding should be reconsidered.
Collapse
Affiliation(s)
- Jianwen A Feng
- Center for Computational Biology, Department of Biochemistry and Molecular Biophysics, Washington University in St. Louis, St. Louis, Missouri, USA.
| | | | | |
Collapse
|
9
|
Kim E, Jang S, Pak Y. All-atom ab initio native structure prediction of a mixed fold (1FME): a comparison of structural and folding characteristics of various beta beta alpha miniproteins. J Chem Phys 2010; 131:195102. [PMID: 19929079 DOI: 10.1063/1.3266510] [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
We performed an all-atom ab initio native structure prediction of 1FME, which is one of the computationally challenging mixed fold beta beta alpha miniproteins, by combining a novel conformational search algorithm (multiplexed Q-replica exchange molecular dynamics scheme) with a well-balanced all-atom force field employing a generalized Born implicit solvation model (param99MOD5/GBSA). The nativelike structure of 1FME was identified from the lowest free energy minimum state and in excellent agreement with the NMR structure. Based on the interpretation of the free energy landscape, the structural properties as well as the folding behaviors of 1FME were compared with other beta beta alpha miniproteins (1FSD, 1PSV, and BBA5) that we have previously studied with the same force field. Our simulation showed that the 28-residue beta beta alpha miniproteins (1FME, 1FSD, and 1PSV) share a common feature of the free energy topography and exhibit the three local minimum states on each computed free energy map, but the 23-residue miniprotein (BBA5) follows a downhill folding with a single minimum state. Also, the structure and stability changes resulting from the two point mutation (Gln1-->Glu1 and Ile7-->Tyr7) of 1FSD were investigated in details for direct comparison with the experiment. The comparison shows that upon mutation, the experimentally observed turn type switch from an irregular turn (1FSD) to type I(') turn (1FME) was well reproduced with the present simulation.
Collapse
Affiliation(s)
- Eunae Kim
- Department of Chemistry, Pusan National University, Busan 609-735, Korea
| | | | | |
Collapse
|
10
|
Lei H, Wang ZX, Wu C, Duan Y. Dual folding pathways of an alpha/beta protein from all-atom ab initio folding simulations. J Chem Phys 2009; 131:165105. [PMID: 19894980 PMCID: PMC2780466 DOI: 10.1063/1.3238567] [Citation(s) in RCA: 21] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/04/2009] [Accepted: 09/07/2009] [Indexed: 11/14/2022] Open
Abstract
Successful ab initio folding of proteins with both alpha-helix and beta-sheet requires a delicate balance among a variety of forces in the simulation model, which may explain that the successful folding of any alpha/beta proteins to within experimental error has yet to be reported. Here we demonstrate that it is an achievable goal to fold alpha/beta proteins with a force field emphasizing the balance between the two major secondary structures. Using our newly developed force field, we conducted extensive ab initio folding simulations on an alpha/beta protein full sequence design (FSD) employing both conventional molecular dynamics and replica exchange molecular dynamics in combination with a generalized-Born solvation model. In these simulations, the folding of FSD to the native state with high population (>64.2%) and high fidelity (C(alpha)-Root Mean Square Deviation of 1.29 A for the most sampled conformation when compared to the experimental structure) was achieved. The folding of FSD was found to follow two pathways. In the major pathway, the folding started from the formation of the helix. In the minor pathway, however, folding of the beta-hairpin started first. Further examination revealed that the helix initiated from the C-terminus and propagated toward the N-terminus. The formation of the hydrophobic contacts coincided with the global folding. Therefore the hydrophobic force does not appear to be the driving force of the folding of this protein.
Collapse
Affiliation(s)
- Hongxing Lei
- Beijing Institute of Genomics, Chinese Academy of Sciences, Beijing 100029, China
| | | | | | | |
Collapse
|
11
|
Chebaro Y, Dong X, Laghaei R, Derreumaux P, Mousseau N. Replica exchange molecular dynamics simulations of coarse-grained proteins in implicit solvent. J Phys Chem B 2009; 113:267-74. [PMID: 19067549 DOI: 10.1021/jp805309e] [Citation(s) in RCA: 66] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Abstract
Current approaches aimed at determining the free energy surface of all-atom medium-size proteins in explicit solvent are slow and are not sufficient to converge to equilibrium properties. To ensure a proper sampling of the configurational space, it is preferable to use reduced representations such as implicit solvent and/or coarse-grained protein models, which are much lighter computationally. Each model must be verified, however, to ensure that it can recover experimental structures and thermodynamics. Here we test the coarse-grained implicit solvent OPEP model with replica exchange molecular dynamics (REMD) on six peptides ranging in length from 10 to 28 residues: two alanine-based peptides, the second beta-hairpin from protein G, the Trp-cage and zinc-finger motif, and a dimer of a coiled coil peptide. We show that REMD-OPEP recovers the proper thermodynamics of the systems studied, with accurate structural description of the beta-hairpin and Trp-cage peptides (within 1-2 A from experiments). The light computational burden of REMD-OPEP, which enables us to generate many hundred nanoseconds at each temperature and fully assess convergence to equilibrium ensemble, opens the door to the determination of the free energy surface of larger proteins and assemblies.
Collapse
Affiliation(s)
- Yassmine Chebaro
- Institut de Biologie Physico-Chimique et Universite Paris 7 Denis Diderot, 75005 Paris, France
| | | | | | | | | |
Collapse
|
12
|
Abstract
Computational protein folding can be classified into pathway and sampling approaches. Here, we use the AMBER simulation package as an example to illustrate the protocols for all-atom molecular simulations of protein folding, including system setup, simulation, and analysis. We introduced two traditional pathway approaches: ab inito folding and high-temperature unfolding. The popular replica exchange method was chosen to represent sampling approaches. Our emphasis is placed on the analysis of the simulation trajectories, and some in-depth discussions are provided for commonly encountered problems.
Collapse
Affiliation(s)
- Hongxing Lei
- Genome Center, University of California, Davis, CA, USA
| | | |
Collapse
|
13
|
Wei CC, Ho MH, Wang WH, Sun YC. Molecular dynamics simulation of folding of a short helical peptide with many charged residues. J Phys Chem B 2007; 109:19980-6. [PMID: 16853583 DOI: 10.1021/jp052349k] [Citation(s) in RCA: 11] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
A molecular dynamics simulation of the folding of conantokin-T (con-T), a short helical peptide with 5 helical turns of 21 amino acids with 10 charged residues, was carried out to examine folding pathways for this peptide and to predict the folding rate. In the 18 trajectories run at 300 K, 16 trajectories folded, with an averaged folding time of approximately 50 ns. Two trajectories did not fold in up to 200 ns simulation. The folded structure in folded trajectories is in good agreement with experimental structure. An analysis of the trajectories showed that, at the beginning of a few nanoseconds, helix formation started from residues 5-9 with assistance of a hydrophobic clustering involving Tyr5, Met8, and Leu9. The peptide formed a U-shape mainly due to charge-charge interactions between charged residues at the N- and C-terminus segments. In the next approximately 10 ns, several nonnative charge-charge interactions were broken and nonnative Gla10-Lys18 (this denotes a salt bridge between Gal10 and Lys18) and/or Gla10-Lys19 interactions appeared more frequently in this folding step and the peptide became a fishhook J-shape. From this structure, the peptide folded to the folded state in 7 of all 16 folded trajectories in approximately 15 ns. Alternatively, in approximately 30 ns, the con-T went to a conformation in an L-shape with 4 helical turns and a kink at the Arg13 and Gla14 segment in the other 9 trajectories. Con-T in the L-shape then required another approximately 15 ns to fold into the folded state. In addition, in overall folding times, the former 7 trajectories folded faster with the total folding times all shorter than 45 ns, while the latter 9 trajectories folded at a time longer than 45 ns, resulting in an average folding time of approximately 50 ns. Two major folding intermediates found in 2 nonfolded trajectories are stabilized by charge clusters of 5 and 6 charged residues, respectively. With inclusion of friction and solvent-solvent interactions, which were ignored in the present GB/SA solvation model, the folding time obtained above should be multiplied by a factor of 1.25-1.7 according to a previous, similar simulation study. This results in a folding time of 65-105 ns, slightly shorter than the folding time of 127 ns for an alanine-based peptide of the same length. This suggests that the energy barrier of folding for this type of peptides with many charged residues is slightly lower than alanine-based helical peptides by less than 1 kcal/mol.
Collapse
Affiliation(s)
- Chung-Cheng Wei
- Department of Chemistry, National Taiwan Normal University, 88, TingChow Road Section 4, Taipei 116, Taiwan
| | | | | | | |
Collapse
|
14
|
Abstract
Using the 28 residue betabetaalpha protein FSD-EY as a target system, we examine correction terms for the ECEPP/3 force field. We find an increased probability of formation of the native state at low temperatures resulting from a reduced propensity to form alpha helices and increased formation of beta sheets. Our analysis of the observed folding events suggests that the C-terminal helix of FSD-EY is much more stable than the N-terminal beta hairpin and forms first. The hydrophobic groups of the helix provide a template which promotes the formation of the beta hairpin that is never observed to form without the helix.
Collapse
Affiliation(s)
- Sandipan Mohanty
- John von Neumann Institut für Computing, Forschungszentrum Jülich, Jülich D-52425, Germany.
| | | |
Collapse
|
15
|
Lei H, Duan Y. Two-stage folding of HP-35 from ab initio simulations. J Mol Biol 2007; 370:196-206. [PMID: 17512537 PMCID: PMC2701201 DOI: 10.1016/j.jmb.2007.04.040] [Citation(s) in RCA: 81] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/27/2006] [Revised: 04/10/2007] [Accepted: 04/13/2007] [Indexed: 11/24/2022]
Abstract
Accurate ab initio simulation of protein folding is a critical step toward elucidation of protein-folding mechanisms. Here, we demonstrate highly accurate folding of the 35 residue villin headpiece subdomain (HP35) by all-atom molecular dynamics simulations using AMBER FF03 and the generalized-Born solvation model. In a set of 20 micros long simulations, the protein folded to the native state in multiple trajectories, with the lowest C(alpha) RMSD being 0.39 A for residues 2-34 (excluding residues 1 and 35). The native state had the highest population among all sampled conformations, and the center of most populated cluster had a C(alpha) RMSD of 1.63 A. Folding of this protein can be described as a two-stage process that followed a well-defined pathway. In the first stage, formation of helices II and III as a folding intermediate constituted the rate-limiting step and was initiated at a folding nucleus around residues Phe17 and Pro21. The folding intermediate further acted as a template that facilitated the folding and docking of helix I in the second stage. Detailed descriptions of the folding kinetics and the roles of key residues are presented.
Collapse
Affiliation(s)
- Hongxing Lei
- Genome Center and Department of Applied Science, University of California Davis, One Shields Avenue, Davis, CA 95616, USA
| | | |
Collapse
|
16
|
Wang WZ, Lin T, Sun YC. Examination of the folding of a short alanine-based helical peptide with salt bridges using molecular dynamics simulation. J Phys Chem B 2007; 111:3508-14. [PMID: 17388513 DOI: 10.1021/jp067637a] [Citation(s) in RCA: 15] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
A molecular dynamics simulation of the folding of a short alanine-based helical peptide of 17 residues with three Glu...Lys (i, i + 4) salt bridge pairs, referred to as the AEK17 peptide, was carried out. The simulation gave an estimated simulation folding time of 2.5 ns, shorter than 12 ns for an alanine-based peptide of 16 residues with three Lys residues only, referred to as the AK16 peptide, simulated previously. After folded, the AEK17 peptide had a helical content of 77%, in excellent agreement with the experimentally determined value of 80%. An examination of the folding pathways of AEK17 indicated that the peptide proceeded via three-turn helix conformations more than the helix-turn-helix conformation in the folding pathways. An analysis of interactions indicated that the formation of hydrogen bonds between Lys residue side chains and backbone carbonyls is a major factor in the abundant conformation of the three-turn helix intermediate. The substitution of three Ala with Glu residues reduces the extent of hydrophobic interaction in alanine-based AK peptides with the result that the breaking of the interactions of Lys epsilon-NH3+(side chain)...C=O(backbone) is a major activation action for the AEK17 to achieve a complete fold, in contrast to the AK16 peptide, in which breaking non-native hydrophobic interaction is the rate-determining step.
Collapse
Affiliation(s)
- Wei-Zhou Wang
- Department of Chemistry, National Taiwan Normal University, 88, TingChow Road Section 4, Taipei 116, Taiwan
| | | | | |
Collapse
|
17
|
Lei H, Wu C, Liu H, Duan Y. Folding free-energy landscape of villin headpiece subdomain from molecular dynamics simulations. Proc Natl Acad Sci U S A 2007; 104:4925-30. [PMID: 17360390 PMCID: PMC1829241 DOI: 10.1073/pnas.0608432104] [Citation(s) in RCA: 186] [Impact Index Per Article: 10.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
High-accuracy ab initio folding has remained an elusive objective despite decades of effort. To explore the folding landscape of villin headpiece subdomain HP35, we conducted two sets of replica exchange molecular dynamics for 200 ns each and three sets of conventional microsecond-long molecular dynamics simulations, using AMBER FF03 force field and a generalized-Born solvation model. The protein folded consistently to the native state; the lowest C(alpha)-rmsd from the x-ray structure was 0.46 A, and the C(alpha)- rmsd of the center of the most populated cluster was 1.78 A at 300 K. ab initio simulations have previously not reached this level. The folding landscape of HP35 can be partitioned into the native, denatured, and two intermediate-state regions. The native state is separated from the major folding intermediate state by a small barrier, whereas a large barrier exists between the major folding intermediate and the denatured states. The melting temperature T(m) = 339 K extracted from the heat-capacity profile was in close agreement with the experimentally derived T(m) = 342 K. A comprehensive picture of the kinetics and thermodynamics of HP35 folding emerges when the results from replica exchange and conventional molecular dynamics simulations are combined.
Collapse
Affiliation(s)
- Hongxing Lei
- Genome Center and Department of Applied Science, University of California, Davis, CA 95616
| | - Chun Wu
- Genome Center and Department of Applied Science, University of California, Davis, CA 95616
| | - Haiguang Liu
- Genome Center and Department of Applied Science, University of California, Davis, CA 95616
| | - Yong Duan
- Genome Center and Department of Applied Science, University of California, Davis, CA 95616
- *To whom correspondence should be addressed. E-mail:
| |
Collapse
|
18
|
Lei H, Dastidar SG, Duan Y. Folding transition-state and denatured-state ensembles of FSD-1 from folding and unfolding simulations. J Phys Chem B 2007; 110:22001-8. [PMID: 17064170 PMCID: PMC3382983 DOI: 10.1021/jp063716a] [Citation(s) in RCA: 12] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Characterization of the folding transition-state ensemble and the denatured-state ensemble is an important step toward a full elucidation of protein folding mechanisms. We report herein an investigation of the free-energy landscape of FSD-1 protein by a total of four sets of folding and unfolding molecular dynamics simulations with explicit solvent. The transition-state ensemble was initially identified from unfolding simulations at 500 K and was verified by simulations at 300 K starting from the ensemble structures. The denatured-state ensemble and the early-stage folding were studied by a combination of unfolding simulations at 500 K and folding simulations at 300 K starting from the extended conformation. A common feature of the transition-state ensemble was the substantial formation of the native secondary structures, including both the alpha-helix and beta-sheet, with partial exposure of the hydrophobic core in the solvent. Both the native and non-native secondary structures were observed in the denatured-state ensemble and early-stage folding, consistent with the smooth experimental melting curve. Interestingly, the contact orders of the transition-state ensemble structures were similar to that of the native structure and were notably lower than those of the compact structures found in early-stage folding, implying that chain and topological entropy might play significant roles in protein folding. Implications for FSD-1 folding mechanisms and the rate-limiting step are discussed. Analyses further revealed interesting non-native interactions in the denatured-state ensemble and early-stage folding and the possibility that destabilization of these interactions could help to enhance the stability and folding rate of the protein.
Collapse
Affiliation(s)
| | | | - Yong Duan
- Corresponding author Telephone: (530) 754-7632, Fax: (530) 754-9648,
| |
Collapse
|
19
|
Li W, Zhang J, Wang W. Understanding the folding and stability of a zinc finger-based full sequence design protein with replica exchange molecular dynamics simulations. Proteins 2007; 67:338-49. [PMID: 17285627 DOI: 10.1002/prot.21312] [Citation(s) in RCA: 22] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/05/2022]
Abstract
Full sequence design protein FSD-1 is a designed protein based on the motif of zinc finger protein. In this work, its folding mechanism and thermal stability are investigated using the replica exchange molecular dynamics model with the water molecules being treated explicitly. The results show that the folding of the FSD-1 is initiated by the hydrophobic collapse, which is accompanied with the formation of the C-terminal alpha-helix. Then the folding proceeds with the formation of the beta-hairpin and the further package of the hydrophobic core. Compared with the beta-hairpin, the alpha-helix has much higher stability. It is also found that the N-capping motif adopted by the FSD-1 contributes to the stability of the alpha-helix dramatically. The hydrophobic contacts made by the side chain of Tyr3 in the native state are essential for the stabilization of the beta-hairpin. It is also found that the folding of the N-terminal beta-hairpin and the C-terminal alpha-helix exhibits weak cooperativity, which is consistent with the experimental data. Meanwhile, the folding pathway is compared between the FSD-1 and the target zinc finger peptide, and the possible role of the zinc ion on the folding pathway of zinc finger is proposed. Proteins 2007. (c) 2007 Wiley-Liss, Inc.
Collapse
Affiliation(s)
- Wenfei Li
- National laboratory of Solid State Microstructure and Department of Physics, Nanjing University, Nanjing 210093, China
| | | | | |
Collapse
|
20
|
|
21
|
Jang S, Kim E, Pak Y. Free energy surfaces of miniproteins with a ββα motif: Replica exchange molecular dynamics simulation with an implicit solvation model. Proteins 2005; 62:663-71. [PMID: 16329109 DOI: 10.1002/prot.20771] [Citation(s) in RCA: 33] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Abstract
Designed miniproteins with a betabetaalpha motif, such as BBA5, 1FSD, and 1PSV can serve as a benchmark set to test the validity of all-atom force fields with computer simulation, because they contain all the basic structural elements in protein folding. Unfortunately, it was found that the standard all-atom force fields with the generalized Born (GB) implicit solvation model tend to produce distorted free energy surfaces for the betabetaalpha proteins, not only because energetically those proteins need to be described by more balanced weights of the alpha- and beta-strands, but also because the GB implicit solvation model suffers from overestimated salt bridge effects. In an attempt to resolve these problems, we have modified one of the standard all-atom force fields in conjunction with the GB model, such that each native state of the betabetaalpha proteins is in its free energy minimum state with reasonable energy barriers separating local minima. With this modified energy model, the free energy contour map in each protein was constructed from the replica exchange molecular dynamics REMD simulation. The resulting free energy surfaces are significantly improved in comparison with previous simulation results and consistent with general views on small protein folding behaviors with realistic topology and energetics of all three proteins.
Collapse
Affiliation(s)
- Soonmin Jang
- Department of Chemistry, Seoul National University, Seoul, Korea
| | | | | |
Collapse
|