Collapse Kinetics and Chevron Plots from Simulations of Denaturant-Dependent Folding of Globular Proteins

Total Page:16

File Type:pdf, Size:1020Kb

Collapse Kinetics and Chevron Plots from Simulations of Denaturant-Dependent Folding of Globular Proteins Collapse kinetics and chevron plots from simulations of denaturant-dependent folding of globular proteins Zhenxing Liua, Govardhan Reddya, Edward P. O’Briena, and D. Thirumalaia,b,1 aBiophysics Program, Institute for Physical Science and Technology, University of Maryland, College Park, MD 20742; and bDepartment of Chemistry, University of Maryland, College Park, MD 20742 Edited* by William A. Eaton, National Institutes of Health-National Institute of Diabetes and Digestive and Kidney Diseases, Bethesda, MD, and approved March 8, 2011 (received for review December 24, 2010) Quantitative description of how proteins fold under experimental (folding arm) followed by an increase (unfolding arm) at high conditions remains a challenging problem. Experiments often use ½C. In previous computational studies, chevron plots were gen- urea and guanidinium chloride to study folding whereas the nat- erated by mapping temperature or certain interaction parameters ural variable in simulations is temperature. To bridge the gap, we in the force-field to ½C (16). Although these studies have pro- use the molecular transfer model that combines measured denatur- vided insights into the role of nonnative interactions they fall ant-dependent transfer free energies for the peptide group and short of directly probing folding of proteins in the presence of amino acid residues, and a coarse-grained Cα-side chain model denaturants. Here, we solve the long-standing problem of how for polypeptide chains to simulate the folding of src SH3 domain. to calculate chevron plots, and in the process describe the entire Stability of the native state decreases linearly as ½C (the concentra- folding reaction of a small protein in near quantitative detail. tion of guanidinium chloride) increases with the slope, m, that is in Molecular simulations of coarse-grained off-lattice models (14, excellent agreement with experiments. Remarkably, the calculated 24, 25), are ideally suited to make (semi)quantitative predictions folding rate at ½C¼0 is only 16-fold larger than the measured of the folding reaction. A combination of statistical mechanical k k value. Most importantly ln obs ( obs is the sum of folding and approach and simulations using simple off-lattice models has unfolding rates) as a function of ½C has the characteristic V (chev- provided insights into protein folding mechanisms. In particular, ron) shape. In every folding trajectory, the times for reaching the these approaches have shown how linear free energy relations for native state, interactions stabilizing all the substructures, and the seemingly complex folding reactions emerge naturally (26). BIOPHYSICS AND mf COMPUTATIONAL BIOLOGY global collapse coincide. The value of m (mf is the slope of the fold- More importantly, the theoretical ideas rationalize the robustness ing arm of the chevron plot) is identical to the fraction of buried of folding in terms of topology alone as illustrated in the simula- solvent accessible surface area in the structures of the transition tions of structure-based models of small globular proteins (27). state ensemble. In the dominant transition state, which does not Of particular relevance to the present study is the demonstration vary significantly at low ½C, the core of the protein and certain that solvation forces determine the origin of barriers in SH3 loops are structured. Besides solving the long-standing problem domains, which helped rationalize the effect of specific mutations of computing the chevron plot, our work lays the foundation for on the folding rates (28). Here, we extend the molecular transfer CHEMISTRY incorporating denaturant effects in a physically transparent man- model (MTM) (29–31) to map in quantitative detail the thermo- ner either in all-atom or coarse-grained simulations. dynamics and to a lesser extent kinetics of folding of src SH3 do- Δ main as a function of GdmCl and urea. We show that GNUð½CÞ kinetic cooperativity ∣ self-organized polymer model ∣ pathway diversity ∣ decreases linearly as ½C, with the slopes (m) being in quantitative protein denaturation agreement with experiments. Simulations of the ½C-dependent ln kobs has the characteristic chevron shape, and the values of nderstanding how proteins fold in quantitative molecular the slopes of the folding (mf ) and unfolding (mu) arms are in ex- Udetail can provide insights into protein aggregation and dy- cellent agreement with experiments. There is a single dominant namics of formation of multisubunit complexes. As a result there transition state whose structure shows that the core β-sheets are have been intense efforts in deciphering the folding mechanisms well packed and the stiff distal loop is ordered. Contacts involving of proteins using a variety of experiments (1–6), theories (7–13), the hydrophobic core and the other secondary structural ele- and simulations (14–16). Despite these advances, a number of ments, overall collapse, and formation of the entire native struc- issues such as the denaturant-dependent characteristics of the ture form nearly simultaneously at zero and low ½C. Our work unfolded states and the link between protein collapse and fold- solves the long-standing problem of describing the entire dena- ing remain unclear (17–21). Single molecule experiments have turant-dependent folding of a protein, including the creation of provided clear evidence for folding pathway diversity, and have the chevron plot, albeit at a coarse-grained level. shown that polypeptide chains undergo almost continuous col- lapse as the denaturant concentration (½C) is decreased (1, Results 19–21). These and other experiments (5, 22) raise the need for Structure of src SH3. To establish the utility of MTM we computed computational models whose predictions can be directly com- several experimentally measurable thermodynamic and kinetic pared to experiments, which often use denaturants [guanidinium properties of the 56-residue src SH3 domain (Fig. 1A) (29) (see chloride (GdmCl) and urea] to initiate folding and unfolding. Methods for details). The core of SH3 domains, whose folding has Global thermodynamic and kinetic behavior of small proteins been extensively investigated using experiments (32–35) and si- often exhibit two-state behavior (23), which implies that at all values of ½C only the population of native (N) and the unfolded Author contributions: Z.L., G.R., E.P.O, and D.T. designed research; Z.L., G.R., and D.T. (U) states are detectable. Experimentally two-state behavior is performed research; D.T. contributed new reagents/analytic tools; Z.L., G.R., and D.T. Δ characterized by (i) the free energy of stability, GNUð½CÞ analyzed data; and Z.L., G.R., and D.T. wrote the paper. − ð¼ GN ð½CÞ GU ð½CÞÞ,ofN with respect to U, varies linearly The authors declare no conflict of interest. with ½C, and (ii) the logarithm of the relaxation rate kobs, which *This Direct Submission article had a prearranged editor. is the sum of folding and unfolding rates, is V-shaped (“chevron 1To whom correspondence should be addressed. E-mail: [email protected]. ” plot ) when plotted as a function of ½C. For two-state systems This article contains supporting information online at www.pnas.org/lookup/suppl/ chevron plots show a linear decrease in ln kobs as ½C increases doi:10.1073/pnas.1019500108/-/DCSupplemental. www.pnas.org/cgi/doi/10.1073/pnas.1019500108 PNAS Early Edition ∣ 1of6 Downloaded by guest on September 29, 2021 mulations (28, 36–38), consists of two well-packed three stranded Δ − f NBA ð½CÞ f NBAð½CÞ leading to GNUð½CÞ ¼ kBTs lnð1− Þ. The line- β f NBA ð½CÞ -sheets that are connected by RT, n-src, and short distal loops Δ Δ 0 ar fit (Fig. 1C), GNUð½CÞ ¼ GNUð½ Þ þ m½C, yields (Fig. 1A). Δ 0 −3 4 ∕ 1 34 ∕ GNUð½ Þ ¼ . kcal mol and m ¼ . kcal mol·M. We also computed ΔG C using free energy profiles as a function of χ Denaturant-Dependent Folding Thermodynamics. Given the approx- NUð½ Þ (Eq. S16), and obtained m ¼ 1.47 kcal∕mol·M. The experimen- imations in our model (see Methods and SI Text), it is difficult 1 5 ∕ to get absolute free energy differences between the folded and tally inferred m value using circular dichroism is . kcal mol·M 1 6 ∕ unfolded states from simulations. To make comparisons with whereas m ¼ . kcal mol·M using fluorescence measurements. The values from our simulation for m ð1.34–1.47Þ kcal∕mol·M experiments, we choose a simulation temperature, Ts, at which the calculated free energy of stability of the native state (N) with are in excellent agreement with experimental estimates, which Δ − further shows that MTM can predict the global thermodynamic respect to the unfolded state (U), GNUðTsÞ (GN ðTsÞ GU ðTsÞ), 295 Δ properties accurately. and the measured free energy at TEð¼ KÞ GNUðTEÞ coin- Δ Δ 0 Although the calculated value of ΔG 0 is in good agree- cide. The use of GNUðTsÞ¼ GNUðTEÞ (in water with ½C¼ ) NUð½ Þ to fix Ts is equivalent to choosing a overall reference energy ment with the experimental result, it is not as accurate as m-value Δ 295 scale in the simulations. For src SH3, GNUðTE ¼ KÞ¼ because of the known difficulties in obtaining accurate results at −4 1 ∕ 0 339 . kcal mol at ½C¼ (34), which results in Ts ¼ K. low GdmCl concentrations (40). To ensure that MTM can predict Δ 0 Besides the choice of Ts no other parameter is adjusted to fit GNUð½ Þ accurately we calculated f NBAð½CÞ for urea (green 339 Δ any other aspect of src SH3 folding. With Ts ¼ K fixed, we squares in Fig. 1B) from which we obtained GNUð½CÞ (Fig. 1C). Δ 0 −4 2 ∕ calculated the dependence of fraction of molecules in the native The linear fit yields GNUð½ Þ ¼ . kcal mol, which is in S18 basin of attraction (39), f NBAð½CÞ,on½C using Eq.
Recommended publications
  • State Folding of Globular Proteins
    Preprints (www.preprints.org) | NOT PEER-REVIEWED | Posted: 20 February 2020 doi:10.20944/preprints202002.0290.v1 Peer-reviewed version available at Biomolecules 2020, 10, 407; doi:10.3390/biom10030407 Review The Molten Globule, and Two-State vs. Non-Two- State Folding of Globular Proteins Kunihiro Kuwajima,* Department of Physics, School of Science, the University of Tokyo, 7-3-1 Hongo, Bunkyo-ku, Tokyo 113-0033, Japan, and School of Computational Sciences, Korea Institute for Advanced Study (KIAS), Seoul 02455, Korea ; [email protected] * Correspondence: [email protected] Abstract: From experimental studies of protein folding, it is now clear that there are two types of folding behavior, i.e., two-state folding and non-two-state folding, and understanding the relationships between these apparently different folding behaviors is essential for fully elucidating the molecular mechanisms of protein folding. This article describes how the presence of the two types of folding behavior has been confirmed experimentally, and discusses the relationships between the two-state and the non-two-state folding reactions, on the basis of available data on the correlations of the folding rate constant with various structure-based properties, which are determined primarily by the backbone topology of proteins. Finally, a two-stage hierarchical model is proposed as a general mechanism of protein folding. In this model, protein folding occurs in a hierarchical manner, reflecting the hierarchy of the native three-dimensional structure, as embodied in the case of non-two-state folding with an accumulation of the molten globule state as a folding intermediate. The two-state folding is thus merely a simplified version of the hierarchical folding caused either by an alteration in the rate-limiting step of folding or by destabilization of the intermediate.
    [Show full text]
  • Kinetic Folding Mechanism of Erythropoietin
    Biophysical Journal Volume 96 May 2009 4221–4230 4221 Kinetic Folding Mechanism of Erythropoietin Douglas D. Banks,* Joanna L. Scavezze, and Christine C. Siska Department of Analytical and Formulation Sciences, Amgen, Seattle, Washington 98119-3105 ABSTRACT This report describes what to our knowledge is the first kinetic folding studies of erythropoietin, a glycosylated four-helical bundle cytokine responsible for the regulation of red blood cell production. Kinetic responses for folding and unfolding reactions initiated by manual mixing were monitored by far-ultraviolet circular dichroism and fluorescence spectroscopy, and folding reactions initiated by stopped-flow mixing were monitored by fluorescence. The urea concentration dependence of the observed kinetics were best described by a three-state model with a transiently populated intermediate species that is on-pathway and obligatory. This folding scheme was further supported by the excellent agreement between the free energy of unfolding and m-value calculated from the microscopic rate constants derived from this model and these parameters determined from separate equilibrium unfolding experiments. Compared to the kinetics of other members of the four-helical bundle cytokine family, erythro- poietin folding and unfolding reactions were slower and less susceptible to aggregation. We tentatively attribute these slower rates and protection from association events to the large amount of carbohydrate attached to erythropoietin at four sites. INTRODUCTION Erythropoietin (EPO) is a 30 kD hormone produced EPO is an attractive protein-folding model due to its inter- primarily in the kidneys and is responsible for regulating mediate size and high degree of folding reversibility. A red blood cell production by stimulating late erythroid wealth of information has already been learned studying precursor cells (1).
    [Show full text]
  • Protein Folding: New Methods Unveil Rate-Limiting Structures
    UNIVERSITY OF CHICAGO PROTEIN FOLDING: NEW METHODS UNVEIL RATE-LIMITING STRUCTURES A DISSERTATION SUBMITTED TO THE FACULTY OF THE DIVISION OF THE BIOLOGICAL SCIENCES AND THE PRITZKER SCHOOL OF MEDICINE IN CANDIDACY FOR THE DEGREE OF DOCTOR OF PHILOSOPHY DEPARTMENT OF BIOCHEMISTRY AND MOLECULAR BIOLOGY BY BRYAN ANDREW KRANTZ CHICAGO, ILLINOIS AUGUST 2002 Copyright © 2002 by Bryan Andrew Krantz All rights reserved ii TO MY WIFE, KRIS, AND MY FAMILY iii Table of Contents Page List of Figures viii List of Tables x List of Abbreviations xi Acknowledgements xii Abstract xiii Chapter 1.0 Introduction to Protein Folding 1 1.1 Beginnings: Anfinsen and Levinthal 1 1.1.1 Intermediates populated by proline isomerization 2 1.1.2 Intermediates populated by disulfide bond formation 4 1.1.3 H/D labeling folding intermediates 7 1.1.4 Kinetically two-state folding 8 1.2 Present thesis within historical context 11 1.3 The time scales 12 1.3.1 From microseconds to decades 12 1.3.2 A diffusive process 13 1.3.3 Misfolding and chaperones 13 1.4 The energies 14 1.4.1 The hydrophobic effect 14 1.4.2 Conformational entropy 15 1.4.3 Hydrogen bonds 15 1.4.4 Electrostatics 16 1.4.5 Van der Waals 19 1.5 Surface burial 19 1.6 Modeling protein folding 20 1.6.1 Diffusion-collision and hydrophobic collapse 21 1.6.2 Search nucleation 22 1.6.3 Are there intermediates? 23 1.6.4 Pathways or landscapes? 27 2.0 The Initial Barrier Hypothesis in Protein Folding 30 2.1 Abstract 30 2.2 Introduction 30 2.3 Early intermediates 41 2.3.1 Cytochrome c 41 2.3.2 Barnase 46 2.3.3 Ubiquitin
    [Show full text]
  • The Inverted Chevron Plot Measured by NMR Relaxation Reveals a Native-Like Unfolding Intermediate in Acyl-Coa Binding Protein
    The inverted chevron plot measured by NMR relaxation reveals a native-like unfolding intermediate in acyl-CoA binding protein Kaare Teilum*, Flemming M. Poulsen†, and Mikael Akke*‡ *Department of Biophysical Chemistry, Lund University, P.O. Box 124, SE-221 00 Lund, Sweden; and †Institute of Molecular Biology and Physiology, University of Copenhagen, Øster Farimagsgade 2A, DK-1353 Copenhagen, Denmark Edited by Alan R. Fersht, University of Cambridge, Cambridge, United Kingdom, and approved March 14, 2006 (received for review October 18, 2005) The folding kinetics of bovine acyl-CoA binding protein was stud- equilibrium after rapid perturbation of the system. For example, ied by 15N relaxation dispersion measurements under equilibrium the response to a rapid change in denaturant concentration can conditions. Relaxation dispersion profiles were measured at sev- be monitored by using stopped-flow methods. In this fashion, eral concentrations of guanidine hydrochloride (GuHCl). The un- folding (unfolding) rates can be measured sensitively only at folding rate constant (ku) was determined under conditions favor- denaturant concentrations below (above) the midpoint of the ing folding, for which the folding rate constant (kf) dominates the folding transition, where relaxation to the new equilibrium is relaxation in stopped-flow kinetic measurements. Conversely, kf dominated by the folding (unfolding) rate constant. Hence, long was determined under conditions favoring unfolding, for which ku extrapolations are usually necessary to extract the unfolding dominates stopped-flow data. The rates determined by NMR there- rates under physiological conditions at zero denaturant concen- fore complement those from stopped-flow kinetics and define an tration. Protein unfolding under conditions favoring the native ‘‘inverted chevron’’ plot.
    [Show full text]
  • Global Analysis of Protein Stability by Temperature and Chemical Denaturation
    bioRxiv preprint doi: https://doi.org/10.1101/2020.04.20.049429; this version posted April 20, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY-NC-ND 4.0 International license. Global analysis of protein stability by temperature and chemical denaturation Louise Hamborg, Emma Wenzel Horsted, Kristoffer Enøe Johansson, Martin Willemoës, Kresten Lindorff-Larsen and Kaare Teilum* Structural Biology and NMR laboratory and the Linderstrøm-Lang Centre for Protein Science, Department of Biology, University of Copenhagen, Ole Maaloes Vej 5, 2200 Copenhagen N, Denmark * Corresponding author Phone: +45 35 32 20 29 Postal Address: Ole Maaloes Vej 5, 2200 Copenhagen N, Denmark Email: [email protected] 1 bioRxiv preprint doi: https://doi.org/10.1101/2020.04.20.049429; this version posted April 20, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY-NC-ND 4.0 International license. Abstract The stability of a protein is a fundamental property that determines under which conditions, the protein is functional. Equilibrium unfolding with denaturants requires preparation of several samples and only provides the free energy of folding when performed at a single temperature. The typical sample requirement is around 0.5 – 1 mg of protein. If the stability of many proteins or protein variants needs to be determined, substantial protein production may be needed.
    [Show full text]
  • Identification and Characterization of Protein Folding Intermediates Stefano Gianni, Ylva Ivarsson, Per Jemth, Maurizio Brunori, Carlo Travaglini-Allocatelli
    Identification and characterization of protein folding intermediates Stefano Gianni, Ylva Ivarsson, Per Jemth, Maurizio Brunori, Carlo Travaglini-Allocatelli To cite this version: Stefano Gianni, Ylva Ivarsson, Per Jemth, Maurizio Brunori, Carlo Travaglini-Allocatelli. Identifica- tion and characterization of protein folding intermediates. Biophysical Chemistry, Elsevier, 2007, 128 (2-3), pp.105. 10.1016/j.bpc.2007.04.008. hal-00501662 HAL Id: hal-00501662 https://hal.archives-ouvertes.fr/hal-00501662 Submitted on 12 Jul 2010 HAL is a multi-disciplinary open access L’archive ouverte pluridisciplinaire HAL, est archive for the deposit and dissemination of sci- destinée au dépôt et à la diffusion de documents entific research documents, whether they are pub- scientifiques de niveau recherche, publiés ou non, lished or not. The documents may come from émanant des établissements d’enseignement et de teaching and research institutions in France or recherche français ou étrangers, des laboratoires abroad, or from public or private research centers. publics ou privés. ÔØ ÅÒÙ×Ö ÔØ Identification and characterization of protein folding intermediates Stefano Gianni, Ylva Ivarsson, Per Jemth, Maurizio Brunori, Carlo Travaglini- Allocatelli PII: S0301-4622(07)00087-7 DOI: doi: 10.1016/j.bpc.2007.04.008 Reference: BIOCHE 4952 To appear in: Biophysical Chemistry Received date: 29 March 2007 Revised date: 16 April 2007 Accepted date: 16 April 2007 Please cite this article as: Stefano Gianni, Ylva Ivarsson, Per Jemth, Maurizio Brunori, Carlo Travaglini-Allocatelli, Identification and characterization of protein folding inter- mediates, Biophysical Chemistry (2007), doi: 10.1016/j.bpc.2007.04.008 This is a PDF file of an unedited manuscript that has been accepted for publication.
    [Show full text]
  • GB1 Is Not a Two-State Folder: Identification and Characterization of an On-Pathway Intermediate
    GB1 is not a two-state folder: identification and characterization of an on-pathway intermediate. Angela Morrone, Rajanish Giri, Rudesh D Toofanny, Carlo Travaglini-Allocatelli, Maurizio Brunori, Valerie Daggett, Stefano Gianni To cite this version: Angela Morrone, Rajanish Giri, Rudesh D Toofanny, Carlo Travaglini-Allocatelli, Maurizio Brunori, et al.. GB1 is not a two-state folder: identification and characterization of an on-pathway intermedi- ate.. Biophysical Journal, Biophysical Society, 2011, 101 (8), pp.2053-60. 10.1016/j.bpj.2011.09.013. pasteur-00982098 HAL Id: pasteur-00982098 https://hal-riip.archives-ouvertes.fr/pasteur-00982098 Submitted on 23 Apr 2014 HAL is a multi-disciplinary open access L’archive ouverte pluridisciplinaire HAL, est archive for the deposit and dissemination of sci- destinée au dépôt et à la diffusion de documents entific research documents, whether they are pub- scientifiques de niveau recherche, publiés ou non, lished or not. The documents may come from émanant des établissements d’enseignement et de teaching and research institutions in France or recherche français ou étrangers, des laboratoires abroad, or from public or private research centers. publics ou privés. Biophysical Journal Volume 101 October 2011 2053–2060 2053 GB1 Is Not a Two-State Folder: Identification and Characterization of an On-Pathway Intermediate Angela Morrone,†6 Rajanish Giri,†6 Rudesh D. Toofanny,‡ Carlo Travaglini-Allocatelli,† Maurizio Brunori,† Valerie Daggett,‡* and Stefano Gianni†* † Istituto Pasteur-Fondazione Cenci Bolognetti and Istituto di Biologia e Patologia Molecolari del CNR, Dipartimento di Scienze Biochimiche ‡ ‘‘A. Rossi Fanelli’’, Universita` di Roma ‘‘La Sapienza’’, Rome, Italy; and Department of Bioengineering, University of Washington, Seattle, Washington ABSTRACT The folding pathway of the small a/b protein GB1 has been extensively studied during the past two decades using both theoretical and experimental approaches.
    [Show full text]
  • Martin Kurnik
    PROTEIN FOLDING WITHOUT LOOPS AND CHARGES Martin Kurnik Protein folding without loops and charges Martin Kurnik ©Martin Kurnik, Stockholm 2012 ISBN 978-91-7447-545-6, pp. 1-70 Printed in Sweden by US-AB, Stockholm 2012 Distributor: Department of Biochemistry and Biophysics, Stockholm University Till M M M List of Publications Primary Publications I Danielsson J, Kurnik M, Lang L, Oliveberg M (2011) Cutting off functional loops from homodimeric enzyme superoxide dismutase 1 (SOD1) leaves monomeric β-barrels. J Biol Chem 286:33070–33083. II Kurnik M, Saraboji, K, Logan, DT, Oliveberg M (2012) Curved un- folding kinetics of superoxide dismutase 1 originate from the central β-barrel. Manuscript. III Lang L, Kurnik M, Danielsson J, Oliveberg M (2012) Fibrillation precursor of superoxide dismutase 1 revealed by gradual tuning of the protein-folding equilibrium. Proc Natl Acad Sci USA. Jul 13 Epub ahead of print. IV Kurnik M, Hedberg L, Danielsson J, Oliveberg M (2012) Folding without charges. Proc Natl Acad Sci USA 109:5705–5710. Additional Publications Haglund E*, Kurnik M*, Oliveberg M (2012) Probing entropy by loop elongation. Manuscript. Björnerås J, Kurnik M, Oliveberg M, Gräslund A, Mäler L, Dan- ielsson J (2012) Direct detection of neuropeptide Dynorphin A bind- ing to the second extracellular loop of the κ-opioid receptor using a soluble protein scaffold and NMR relaxation. Manuscript. * these authors contributed equally Abstract Going down the folding funnel, proteins may sample a wide variety of con- formations, some being outright detrimental to the organism. Yet, the vast majority of polypeptide molecules avoid such pitfalls. Not only do they reach the native minimum of the energy landscape; they do so via blazingly fast, biased, routes.
    [Show full text]
  • The Transition State for Integral Membrane Protein Folding
    The transition state for integral membrane protein folding Paul Curnow1 and Paula J. Booth1 Department of Biochemistry, School of Medical Sciences, University of Bristol, University Walk, Bristol BS8 1TD, United Kingdom Edited by Alan Fersht, University of Cambridge, Cambridge, United Kingdom, and approved November 21, 2008 (received for review July 18, 2008) Biology relies on the precise self-assembly of its molecular com- water-soluble proteins. In particular, there are few quantitative ponents. Generic principles of protein folding have emerged from studies of the kinetics and thermodynamics of membrane protein extensive studies on small, water-soluble proteins, but it is unclear folding. Yet, understanding these folding parameters gives valu- how these ideas are translated into more complex situations. In able information on a major problem: how to fold and stabilize particular, the one-third of cellular proteins that reside in biological membrane proteins for structural and functional studies. membranes will not fold like water-soluble proteins because mem- A ⌽ value is a measure of the change in activation energy, brane proteins need to expose, not hide, their hydrophobic sur- relative to the change in the overall free energy of folding, faces. Here, we apply the powerful protein engineering method of induced by the directed mutation of a single amino acid. ⌽ values ⌽-value analysis to investigate the folding transition state of the are expected to fall between one and zero (these extreme values alpha-helical membrane protein, bacteriorhodopsin, from a par- mean the change in the free energy of the transition state is the tially unfolded state. Our results imply that much of helix B of the same as either that in the folded or unfolded state, respectively).
    [Show full text]
  • Solvation Effects and Driving Forces for Protein Thermodynamic and Kinetic Cooperativity: How Adequate Is Native-Centric Topological Modeling?
    November 29, 2002 Solvation Effects and Driving Forces for Protein Thermodynamic And Kinetic Cooperativity: How Adequate Is Native-Centric Topological Modeling? H¨useyin KAYA and Hue Sun CHAN† Protein Engineering Network of Centres of Excellence (PENCE), Department of Biochemistry, and Department of Medical Genetics & Microbiology, Faculty of Medicine, University of Toronto, Toronto, Ontario M5S 1A8, Canada Running title: Continuum G¯oModel Chevron Plots Key words: calorimetric cooperativity / single-exponential kinetics / unfolding / chevron plot / desolvation barrier / continuum G¯omodels / heat capacity † Corresponding author. E-mail address of Hue Sun CHAN: [email protected] arXiv:cond-mat/0212105v2 [cond-mat.stat-mech] 12 Dec 2002 Tel: (416)978-2697; Fax: (416)978-8548 Mailing address: Department of Biochemistry, University of Toronto, Medical Sciences Building – 5th Fl., 1 King’s College Circle, Toronto, Ontario M5S 1A8, Canada. 1 Summary What energetic and solvation effects underlie the remarkable two-state thermodynamics and folding/unfolding kinetics of small single-domain proteins? To address this question, we investigate the folding and unfolding of a hierarchy of continuum Langevin dynamics models of chymotrypsin inhibitor 2. We find that residue-based additive G¯o-like contact energies, although native-centric, are by themselves insufficient for proteinlike calorimet- ric two-state cooperativity. Further native biases by local conformational preferences are necessary for proteinlike thermodynamics. Kinetically, however, even models with both contact and local native-centric energies do not produce simple two-state chevron plots. Thus a model protein’s thermodynamic cooperativity is not sufficient for sim- ple two-state kinetics. The models tested appear to have increasing internal friction with increasing native stability, leading to chevron rollovers that typify kinetics that are commonly referred to as non-two-state.
    [Show full text]
  • Probing Possible Downhill Folding: Native Contact Topology Likely Places a Significant Constraint on the Folding Cooperativity of Proteins with ∼40 Residues
    doi:10.1016/j.jmb.2008.09.023 J. Mol. Biol. (2008) 384, 512–530 Available online at www.sciencedirect.com Probing Possible Downhill Folding: Native Contact Topology Likely Places a Significant Constraint on the Folding Cooperativity of Proteins with ∼40 Residues Artem Badasyan, Zhirong Liu and Hue Sun Chan⁎ Department of Biochemistry Experiments point to appreciable variations in folding cooperativity among and Department of Molecular natural proteins with approximately 40 residues, indicating that the Genetics, Faculty of Medicine, behaviors of these proteins are valuable for delineating the contributing University of Toronto, Toronto, factors to cooperative folding. To explore the role of native topology in a Ontario, Canada M5S 1A8 protein's propensity to fold cooperatively and how native topology might constrain the degree of cooperativity achievable by a given set of physical Received 16 July 2008; interactions, we compared folding/unfolding kinetics simulated using received in revised form α three classes of native–centric C chain models with different interaction 6 September 2008; schemes. The approach was applied to two homologous 45-residue accepted 10 September 2008 fragments from the peripheral subunit-binding domain family and a 39- Available online residue fragment of the N-terminal domain of ribosomal protein L9. Free- 17 September 2008 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.
    [Show full text]
  • Criteria for Downhill Protein Folding: Calorimetry, Chevron Plot, Kinetic
    PROTEINS: Structure, Function, and Bioinformatics 65:373–391 (2006) Criteria for Downhill Protein Folding: Calorimetry, Chevron Plot, Kinetic Relaxation, and Single-Molecule Radius of Gyration in Chain Models With Subdued Degrees of Cooperativity Michael Knott and Hue Sun Chan* Department of Biochemistry, and of Medical Genetics and Microbiology, Protein Engineering Network of Centres of Excellence, Faculty of Medicine, University of Toronto, Toronto, Ontario M5S 1A8, Canada ABSTRACT Recent investigations of possible INTRODUCTION downhill folding of small proteins such as BBL Polymer physics is important to the study of proteins have focused on the thermodynamics of non-two- in at least two respects. First, it is a logical starting state, ‘‘barrierless’’ folding/denaturation transitions. point for an analysis aimed at providing a physico-chem- Downhill folding is noncooperative and thermody- ical basis for the structure and stability of proteins.1–6 namically ‘‘one-state,’’ a phenomenon underpinned Second, natural proteins are a subset of all possible by a unimodal conformational distribution over polypeptides. As polymers, polypeptides can, in princi- chain properties such as enthalpy, hydrophobic ex- ple, possess a broad range of properties shared with posure, and conformational dimension. In contrast, other chain molecules. By contrasting the properties of corresponding distributions for cooperative two- natural proteins with the general behavior of polymers, state folding are bimodal with well-separated popu- we can clarify the biological
    [Show full text]