Navigating the downhill protein folding regime via structural homologues

Folding Free Energy Surfaces from Differential Scanning Calorimetry

Protein folding/unfolding processes involve a large number of weak, non-covalent interactions and are more appropriately described in terms of the movement of a point representing protein conformation in a plot of internal free energy versus conformational degrees of freedom. While these energy landscapes have an astronomically large number of dimensions, it has been shown that many relevant aspects of protein folding can be understood in terms of their projections onto a few relevant coordinates. Remarkably, such low-dimensional free energy surfaces can be obtained from experimental DSC data using suitable analytical models. Here, we describe the experimental procedures to be used to obtain the high-quality DSC data that are required for free-energy surface analysis.

Lessons about Protein Folding and Binding from Archetypal Folds

The function of proteins as biological nanomachines relies on their ability to fold into complex 3D structures, bind selectively to partners, and undergo conformational changes on cue. The native functional structures, and the rates of interconversion between conformational states (folded-unfolded, bound-free), are all encoded in the physical chemistry of their amino acid sequence. However, despite extensive research over decades, this code has proven difficult to fully crack, in terms of both prediction and understanding the molecular mechanisms at play.Earlier work on single-domain proteins reported a commonality of slow rates (10^-2-10² s^-1) and simple behavior in both kinetic and thermodynamic unfolding experiments, which suggested the process was all-or-none and thereby analogous to a chemical reaction (e.g., A ⇄ B). In the absence of a first-principles pre-exponential factor for protein (un)folding dynamics, the rates could only be interpreted in relative terms, e.g., the changes induced by mutation, and hence, neither the height of nor the entropic contribution to the free energy barriers was known. The rates were also many orders of magnitude too slow for direct atomistic simulations, and the computational focus was on predicting rate changes induced by mutation via coarse grained simulations. However, even the effects of mutation proved to be strikingly homogeneous with all experimental data clustering at ∼1/3 of the free energy perturbation recovered on folding and ∼2/3 on unfolding.The implementation of ultrafast kinetic methods turned the field upside down because they allowed researchers to measure the time scales of elementary (un)folding motions, which set the pre-exponential factor for protein conformational transitions at ∼1 μs. In parallel, we and others set out to investigate the simplest possible protein structures capable of autonomous folding, which we defined as archetypal folds. The rationale was to recapitulate the hierarchical organization of protein structure, starting from the bottom up. The study of fold archetypes ended up opening new research avenues in protein (un)folding, but also making unexpected connections with the folding upon binding of intrinsically disordered proteins and suggesting their functioning as conformational rheostats.This Account describes our work on the kinetic, thermodynamic, mechanistic, and functional analysis of fold archetypes. We first discuss the kinetic studies, emphasizing their impact on our understanding of (un)folding rates, of barrierless (downhill) folding, and as benchmarks for atomistic simulations. We continue with the thermodynamic analysis, introducing the differential scanning calorimetry, multiprobe, and NMR approaches that we developed to dissect their gradual, minimally cooperative (un)folding transitions and to probe the underlying mechanisms with unprecedented detail. The last two sections cover single-molecule analyses and some recent, mostly computational, results on the exploration of possible biological and technological roles for the gradual conformational transitions of fold archetypes.

Slow Folding of a Helical Protein: Large Barriers, Strong Internal Friction, or a Shallow, Bumpy Landscape?

The rate at which a protein molecule folds is determined by opposing energetic and entropic contributions to the free energy that shape the folding landscape. Delineating the extent to which they impact the diffusional barrier-crossing events, including the magnitude of internal friction and barrier height, has largely been a challenging task. In this work, we extract the underlying thermodynamic and dynamic contributions to the folding rate of an unusually slow-folding helical DNA-binding domain, PurR, which shares the characteristics of ultrafast downhill-folding proteins but nonetheless appears to exhibit an apparent two-state equilibrium. We combine equilibrium spectroscopy, temperature-viscosity-dependent kinetics, statistical mechanical modeling, and coarse-grained simulations to show that the conformational behavior of PurR is highly heterogeneous characterized by a large spread in melting temperatures, marginal thermodynamic barriers, and populated partially structured states. PurR appears to be at the threshold of disorder arising from frustrated electrostatics and weak packing that in turn slows down folding due to a shallow, bumpy landscape and not due to large thermodynamic barriers or strong internal friction. Our work highlights how a strong temperature dependence on the pre-exponential could signal a shallow landscape and not necessarily a slow-folding diffusion coefficient, thus determining the folding timescales of even millisecond folding proteins and hints at possible structural origins for the shallow landscape.

Downhill, Ultrafast and Fast Folding Proteins Revised

Research on the protein folding problem differentiates the protein folding process with respect to the duration of this process. The current structure encoded in sequence dogma seems to be clearly justified, especially in the case of proteins referred to as fast-folding, ultra-fast-folding or downhill. In the present work, an attempt to determine the characteristics of this group of proteins using fuzzy oil drop model is undertaken. According to the fuzzy oil drop model, a protein is a specific micelle composed of bi-polar molecules such as amino acids. Protein folding is regarded as a spherical micelle formation process. The presence of covalent peptide bonds between amino acids eliminates the possibility of free mutual arrangement of neighbors. An example would be the construction of co-micelles composed of more than one type of bipolar molecules. In the case of fast folding proteins, the amino acid sequence represents the optimal bipolarity system to generate a spherical micelle. In order to achieve the native form, it is enough to have an external force field provided by the water environment which directs the folding process towards the generation of a centric hydrophobic core. The influence of the external field can be expressed using the 3D Gaussian function which is a mathematical model of the folding process orientation towards the concentration of hydrophobic residues in the center with polar residues exposed on the surface. The set of proteins under study reveals a hydrophobicity distribution compatible with a 3D Gaussian distribution, taken as representing an idealized micelle-like distribution. The structure of the present hydrophobic core is also discussed in relation to the distribution of hydrophobic residues in a partially unfolded form.

A Disordered Loop Mediates Heterogeneous Unfolding of an Ordered Protein by Altering the Native Ensemble

The high flexibility of long disordered or partially structured loops in folded proteins allows for entropic stabilization of native ensembles. Destabilization of such loops could alter the native ensemble or promote alternate conformations within the native ensemble if the ordered regions themselves are held together weakly. This is particularly true of downhill folding systems that exhibit weak unfolding cooperativity. Here, we combine experimental and computational methods to probe the response of the native ensemble of a helical, downhill folding domain PDD, which harbors an 11-residue partially structured loop, to perturbations. Statistical mechanical modeling points to continuous structural changes on both temperature and mutational perturbations driven by entropic stabilization of partially structured conformations within the native ensemble. Long time-scale simulations of the wild-type protein and two mutants showcase a remarkable conformational switching behavior wherein the parallel helices in the wild-type protein sample an antiparallel orientation in the mutants, with the C-terminal helix and the loop connecting the helices displaying high flexibility, disorder, and non-native interactions. We validate these computational predictions via the anomalous fluorescence of a native tyrosine located at the interface of the helices. Our observations highlight the role of long loops in determining the unfolding mechanisms, sensitivity of the native ensembles to mutational perturbations and provide experimentally testable predictions that can be explored in even two-state folding systems.

Ultrafast folding kinetics of WW domains reveal how the amino acid sequence determines the speed limit to protein folding

Protein (un)folding rates depend on the free-energy barrier separating the native and unfolded states and a prefactor term, which sets the timescale for crossing such barrier or folding speed limit. Because extricating these two factors is usually unfeasible, it has been common to assume a constant prefactor and assign all rate variability to the barrier. However, theory and simulations postulate a protein-specific prefactor that contains key mechanistic information. Here, we exploit the special properties of fast-folding proteins to experimentally resolve the folding rate prefactor and investigate how much it varies among structural homologs. We measure the ultrafast (un)folding kinetics of five natural WW domains using nanosecond laser-induced temperature jumps. All five WW domains fold in microseconds, but with a 10-fold difference between fastest and slowest. Interestingly, they all produce biphasic kinetics in which the slower phase corresponds to reequilibration over the small barrier (<3 RT) and the faster phase to the downhill relaxation of the minor population residing at the barrier top [transition state ensemble (TSE)]. The fast rate recapitulates the 10-fold range, demonstrating that the folding speed limit of even the simplest all-β fold strongly depends on the amino acid sequence. Given this fold's simplicity, the most plausible source for such prefactor differences is the presence of nonnative interactions that stabilize the TSE but need to break up before folding resumes. Our results confirm long-standing theoretical predictions and bring into focus the rate prefactor as an essential element for understanding the mechanisms of folding.

Switching Protein Conformational Substates by Protonation and Mutation

Protein modules that regulate the availability and conformational status of transcription factors determine the rapidity, duration, and magnitude of cellular response to changing conditions. One such system is the single-gene product Cnu, a four-helix bundle transcription co-repressor, which acts as a molecular thermosensor regulating the expression of virulence genes in enterobacteriaceae through modulation of its native conformational ensemble. Cnu and related genes have also been implicated in pH-dependent expression of virulence genes. We hypothesize that protonation of a conserved buried histidine (H45) in Cnu promotes large electrostatic frustration, thus disturbing the H-NS, a transcription factor, binding face. Spectroscopic and calorimetric methods reveal that H45 exhibits a suppressed p K_a of ∼5.1, the protonation of which switches the conformation to an alternate native ensemble in which the fourth helix is disordered. The population redistribution can also be achieved through a mutation H45V, which does not display any switching behavior at pH values greater than 4. The Wako-Saitô-Muñoz-Eaton (WSME) statistical mechanical model predicts specific differences in the conformations and fluctuations of the fourth and first helices of Cnu determining the observed pH response. We validate these predictions through fluorescence lifetime measurements of a sole tryptophan, highlighting the presence of both native and non-native interactions in the regions adjoining the binding face of Cnu. Our combined experimental-computational study thus shows that Cnu acts both as a thermo- and pH-sensor orchestrated via a subtle but quantifiable balance between the weak packing of a structural element and protonation of a buried histidine that promotes electrostatic frustration.

Protein Folding Cooperativity and Thermodynamic Barriers of the Simplest β-Sheet Fold: A Survey of WW Domains

Theory and experiments have shown that microsecond folding proteins exhibit characteristic thermodynamic properties that reflect the limited cooperativity of folding over marginal barriers (downhill folding). Those studies have mostly focused on proteins with large α-helical contents and small size, which tend to be the fastest folders. A key open question is whether such properties are also present in the fastest all-β proteins. We address this issue by investigating the unfolding thermodynamics of a collection of WW domains as representatives of the simplest β-sheet fold. WW domains are small microsecond folders, although they do not fold as fast as their α-helical counterparts. In previous work on the NEDD4-WW4 domain, we reported deviations from two-state thermodynamics that were less apparent and thus suggestive of an incipient downhill scenario. Here we investigate the unfolding thermodynamics of four other WW domains (NEDD4-WW3, YAP65-WW1(L30K), FBP11-WW1, and FBP11-WW2) by performing all of the thermodynamic tests for downhill folding that have been previously developed on α-helical proteins. This set of five WW domains shares low sequence identity and include examples from two specificity classes, thus providing a comprehensive survey. Thermodynamic analysis of the four new WW domains consistently reveals all of the properties of downhill folding equilibria, which are in all cases more marked than what we found before in NEDD4-WW4. Our results show that fast-folding all-β proteins do share limited cooperativity and gradual unfolding thermodynamics with fast α-helical proteins and suggest that the free energy barrier to folding of natural proteins is mostly determined by size and fold topology and much less by the specific amino acid sequence.

Folding mechanisms steer the amyloid fibril formation propensity of highly homologous proteins

Understanding the molecular determinants of fibrillogenesis by studying the aggregation propensities of high homologous proteins with different folding pathways.

Significant advances in the understanding of the molecular determinants of fibrillogenesis can be expected from comparative studies of the aggregation propensities of proteins with highly homologous structures but different folding pathways. Here, we fully characterize, by means of stopped-flow, T-jump, CD and DSC experiments, the unfolding mechanisms of three highly homologous proteins, zinc binding Ros87 and Ml1_53–149 and zinc-lacking Ml4_52–151. The results indicate that the three proteins significantly differ in terms of stability and (un)folding mechanisms. Particularly, Ros87 and Ml1_53–149 appear to be much more stable to guanidine denaturation and are characterized by folding mechanisms including the presence of an intermediate. On the other hand, metal lacking Ml4_52–151 folds according to a classic two-state model. Successively, we have monitored the capabilities of Ros87, Ml4_52–151 and Ml1_53–149 to form amyloid fibrils under native conditions. Particularly, we show, by CD, fluorescence, DLS, TEM and SEM experiments, that after 168 hours, amyloid formation of Ros87 has started, while Ml1_53–149 has formed only amorphous aggregates and Ml4_52–151 is still monomeric in solution. This study shows how metal binding can influence protein folding pathways and thereby control conformational accessibility to aggregation-prone states, which in turn changes aggregation kinetics, shedding light on the role of metal ions in the development of protein deposition diseases.

When fast is better: protein folding fundamentals and mechanisms from ultrafast approaches

Protein folding research stalled for decades because conventional experiments indicated that proteins fold slowly and in single strokes, whereas theory predicted a complex interplay between dynamics and energetics resulting in myriad microscopic pathways. Ultrafast kinetic methods turned the field upside down by providing the means to probe fundamental aspects of folding, test theoretical predictions and benchmark simulations. Accordingly, experimentalists could measure the timescales for all relevant folding motions, determine the folding speed limit and confirm that folding barriers are entropic bottlenecks. Moreover, a catalogue of proteins that fold extremely fast (microseconds) could be identified. Such fast-folding proteins cross shallow free energy barriers or fold downhill, and thus unfold with minimal co-operativity (gradually). A new generation of thermodynamic methods has exploited this property to map folding landscapes, interaction networks and mechanisms at nearly atomic resolution. In parallel, modern molecular dynamics simulations have finally reached the timescales required to watch fast-folding proteins fold and unfold in silico. All of these findings have buttressed the fundamentals of protein folding predicted by theory, and are now offering the first glimpses at the underlying mechanisms. Fast folding appears to also have functional implications as recent results connect downhill folding with intrinsically disordered proteins, their complex binding modes and ability to moonlight. These connections suggest that the coupling between downhill (un)folding and binding enables such protein domains to operate analogically as conformational rheostats.

Interplay between the folding mechanism and binding modes in folding coupled to binding processes

Proteins that fold upon binding to their partners exhibit complex binding behavior such as induced-fit. Using molecular simulations we find that the ability to bind via induced-fit requires the disordered protein to fold downhill. Thus induced-fit appears as a manifestation of a conformational rheostat.

Toward a quantitative description of microscopic pathway heterogeneity in protein folding

Experimentally consistent statistical modeling of protein folding thermodynamics reveals unprecedented complexity with numerous parallel folding routes in five different proteins.

Activation Barrier-Limited Folding and Conformational Sampling of a Dynamic Protein Domain

Folding reaction mechanisms of globular protein domains have been extensively studied by both experiment and simulation and found to be highly concerted chemical reactions in which numerous noncovalent bonds form in an apparent two-state fashion. However, less is known regarding intrinsically disordered proteins because their folding can usually be studied only in conjunction with binding to a ligand. We have investigated by kinetics the folding mechanism of such a disordered protein domain, the nuclear coactivator-binding domain (NCBD) from CREB-binding protein. While a previous computational study suggested that NCBD folds without an activation free energy barrier, our experimental data demonstrate that NCBD, despite its highly dynamic structure, displays relatively slow folding (∼10 ms at 277 K) consistent with a barrier-limited process. Furthermore, the folding kinetics corroborate previous nuclear magnetic resonance data showing that NCBD exists in two folded conformations and one more denatured conformation at equilibrium and, thus, that the folding mechanism is a three-state mechanism. The refolding kinetics is limited by unfolding of the less populated folded conformation, suggesting that the major route for interconversion between the two folded states is via the denatured state. Because the two folded conformations have been suggested to bind distinct ligands, our results have mechanistic implications for conformational sampling in protein-protein interactions.

Imprints of function on the folding landscape: functional role for an intermediate in a conserved eukaryotic binding protein

In the computational characterization of single domain protein folding, the effective free energies of numerous microstates are projected onto few collective degrees of freedom that in turn serve as well-defined reaction coordinates. In this regard, one-dimensional (1D) free energy profiles are widely used mainly for their simplicity. Since folding and functional landscapes are interlinked, how well can these reduced representations capture the structural and dynamic features of functional states while being simultaneously consistent with experimental observables? We investigate this issue by characterizing the folding of the four-helix bundle bovine acyl-CoA binding protein (bACBP), which exhibits complex equilibrium and kinetic behaviours, employing an Ising-like statistical mechanical model and molecular simulations. We show that the features of the 1D free energy profile are sufficient to quantitatively reproduce multiple experimental observations including millisecond chevron-like kinetics and temperature dependence, a microsecond fast phase, barrier heights, unfolded state movements, the intermediate structure and average ϕ-values. Importantly, we find that the structural features of the native-like intermediate (partial disorder in helix 1) are intricately linked to a unique interplay between packing and electrostatics in this domain. By comparison with available experimental data, we propose that this intermediate determines the promiscuous functional behaviour of bACBP that exhibits broad substrate specificity. Our results present evidence to the possibility of employing the statistical mechanical model and the resulting 1D free energy profile to not just understand folding mechanisms but to even extract features of functionally relevant states and their energetic origins.

Sequence, structure, and cooperativity in folding of elementary protein structural motifs

Residue-level unfolding of two helix-turn-helix proteins--one naturally occurring and one de novo designed--is reconstructed from multiple sets of site-specific (13)C isotopically edited infrared (IR) and circular dichroism (CD) data using Ising-like statistical-mechanical models. Several model variants are parameterized to test the importance of sequence-specific interactions (approximated by Miyazawa-Jernigan statistical potentials), local structural flexibility (derived from the ensemble of NMR structures), interhelical hydrogen bonds, and native contacts separated by intervening disordered regions (through the Wako-Saitô-Muñoz-Eaton scheme, which disallows such configurations). The models are optimized by directly simulating experimental observables: CD ellipticity at 222 nm for model proteins and their fragments and (13)C-amide I' bands for multiple isotopologues of each protein. We find that data can be quantitatively reproduced by the model that allows two interacting segments flanking a disordered loop (double sequence approximation) and incorporates flexibility in the native contact maps, but neither sequence-specific interactions nor hydrogen bonds are required. The near-identical free energy profiles as a function of the global order parameter are consistent with expected similar folding kinetics for nearly identical structures. However, the predicted folding mechanism for the two motifs is different, reflecting the order of local stability. We introduce free energy profiles for "experimental" reaction coordinates--namely, the degree of local folding as sensed by site-specific (13)C-edited IR, which highlight folding heterogeneity and contrast its overall, average description with the detailed, local picture.

A disorder-induced domino-like destabilization mechanism governs the folding and functional dynamics of the repeat protein IκBα

The stability of the repeat protein IκBα, a transcriptional inhibitor in mammalian cells, is critical in the functioning of the NF-κB signaling module implicated in an array of cellular processes, including cell growth, disease, immunity and apoptosis. Structurally, IκBα is complex, with both ordered and disordered regions, thus posing a challenge to the available computational protocols to model its conformational behavior. Here, we introduce a simple procedure to model disorder in systems that undergo binding-induced folding that involves modulation of the contact map guided by equilibrium experimental observables in combination with an Ising-like Wako-Saitô-Muñoz-Eaton model. This one-step procedure alone is able to reproduce a variety of experimental observables, including ensemble thermodynamics (scanning calorimetry, pre-transitions, m-values) and kinetics (roll-over in chevron plot, intermediates and their identity), and is consistent with hydrogen-deuterium exchange measurements. We further capture the intricate distance-dynamics between the domains as measured by single-molecule FRET by combining the model predictions with simple polymer physics arguments. Our results reveal a unique mechanism at work in IκBα folding, wherein disorder in one domain initiates a domino-like effect partially destabilizing neighboring domains, thus highlighting the effect of symmetry-breaking at the level of primary sequences. The offshoot is a multi-state and a dynamic conformational landscape that is populated by increasingly partially folded ensembles upon destabilization. Our results provide, in a straightforward fashion, a rationale to the promiscuous binding and short intracellular half-life of IκBα evolutionarily engineered into it through repeats with variable stabilities and expand the functional repertoire of disordered regions in proteins.

Slow proton transfer coupled to unfolding explains the puzzling results of single-molecule experiments on BBL, a paradigmatic downhill folding protein

A battery of thermodynamic, kinetic, and structural approaches has indicated that the small α-helical protein BBL folds-unfolds via the one-state downhill scenario. Yet, single-molecule fluorescence spectroscopy offers a more conflicting view. Single-molecule experiments at pH 6 show a unique half-unfolded conformational ensemble at mid denaturation, whereas other experiments performed at higher pH show a bimodal distribution, as expected for two-state folding. Here we use thermodynamic and laser T-jump kinetic experiments combined with theoretical modeling to investigate the pH dependence of BBL stability, folding kinetics and mechanism within the pH 6-11 range. We find that BBL unfolding is tightly coupled to the protonation of one of its residues with an apparent pKa of ~ 7. Therefore, in chemical denaturation experiments around neutral pH BBL unfolds gradually, and also converts in binary fashion to the protonated species. Moreover, under the single-molecule experimental conditions (denaturant midpoint and 279 K), we observe that proton transfer is much slower than the ~ 15 microseconds folding-unfolding kinetics of BBL. The relaxation kinetics is distinctly biphasic, and the overall relaxation time (i.e. 0.2-0.5 ms) becomes controlled by the proton transfer step. We then show that a simple theoretical model of protein folding coupled to proton transfer explains quantitatively all these results as well as the two sets of single-molecule experiments, including their more puzzling features. Interestingly, this analysis suggests that BBL unfolds following a one-state downhill folding mechanism at all conditions. Accordingly, the source of the bimodal distributions observed during denaturation at pH 7-8 is the splitting of the unique conformational ensemble of BBL onto two slowly inter-converting protonation species. Both, the unprotonated and protonated species unfold gradually (one-state downhill), but they exhibit different degree of unfolding at any given condition because the native structure is less stable for the protonated form.

Predictions from an Ising-like Statistical Mechanical Model on the Dynamic and Thermodynamic Effects of Protein Surface Electrostatics

Charged residues on the surface of a protein are known hot-spots for post-translational modification, protein/ligand-binding, and tuning conformational stabilities. Recent experimental evidence points to the fact that surface electrostatics can also modulate thermodynamic barriers and hence folding mechanisms. To probe for this behavior across different proteins, we develop a novel version of the Wako-Saitô-Muñoz-Eaton (WSME) model in which we include an electrostatic potential term in the energy function while simplifying the treatment of solvation free energy. Both of the energy terms are obtained by quantitatively fitting the model to differential scanning calorimetry (DSC) experiments that carry critical information on the protein partition function. We characterize four sets of structural/functional homologues (HEWL/BLA, CspB, engrailed, α-spectrin) either by fitting the experimental data of a single domain in the homologous set and predicting the conformational behavior of the rest with the same set of parameters or by performing semiblind predictions. The model with the added electrostatic term is able to successfully reproduce the order of thermodynamic stabilities and relaxation rates of most of the homologues. In parallel, we predict diverse conformational features including a wide range of thermodynamic barriers (∼9-40 kJ/mol), broad native ensembles in helical proteins, structured unfolded states and intermediates, rugged folding landscapes, and further provide an independent protein-specific estimate of the folding speed limit at 298 K (1/(7-300 μs)). Our results are evidence that protein surface electrostatics can be tailored to not only engineer stabilities but also folding mechanisms and the ruggedness of the underlying landscape.

Exploring one-state downhill protein folding in single molecules

A one-state downhill protein folding process is barrierless at all conditions, resulting in gradual melting of native structure that permits resolving folding mechanisms step-by-step at atomic resolution. Experimental studies of one-state downhill folding have typically focused on the thermal denaturation of proteins that fold near the speed limit (ca. 10(6) s(-1)) at their unfolding temperature, thus being several orders of magnitude too fast for current single-molecule methods, such as single-molecule FRET. An important open question is whether one-state downhill folding kinetics can be slowed down to make them accessible to single-molecule approaches without turning the protein into a conventional activated folder. Here we address this question on the small helical protein BBL, a paradigm of one-state downhill thermal (un)folding. We decreased 200-fold the BBL folding-unfolding rate by combining chemical denaturation and low temperature, and carried out free-diffusion single-molecule FRET experiments with 50-μs resolution and maximal photoprotection using a recently developed Trolox-cysteamine cocktail. These experiments revealed a single conformational ensemble at all denaturing conditions. The chemical unfolding of BBL was then manifested by the gradual change of this unique ensemble, which shifts from high to low FRET efficiency and becomes broader at increasing denaturant. Furthermore, using detailed quantitative analysis, we could rule out the possibility that the BBL single-molecule data are produced by partly overlapping folded and unfolded peaks. Thus, our results demonstrate the one-state downhill folding regime at the single-molecule level and highlight that this folding scenario is not necessarily associated with ultrafast kinetics.

Estimation of protein folding free energy barriers from calorimetric data by multi-model Bayesian analysis

The realization that folding free energy barriers can be small enough to result in significant population of the species at the barrier top has sprouted in several methods to estimate folding barriers from equilibrium experiments. Some of these approaches are based on fitting the experimental thermogram measured by differential scanning calorimetry (DSC) to a one-dimensional representation of the folding free-energy surface (FES). Different physical models have been used to represent the FES: (1) a Landau quartic polynomial as a function of the total enthalpy, which acts as an order parameter; (2) the projection onto a structural order parameter (i.e. number of native residues or native contacts) of the free energy of all the conformations generated by Ising-like statistical mechanical models; and (3) mean-field models that define conformational entropy and stabilization energy as functions of a continuous local order parameter. The fundamental question that emerges is how can we obtain robust, model-independent estimates of the thermodynamic folding barrier from the analysis of DSC experiments. Here we address this issue by comparing the performance of various FES models in interpreting the thermogram of a protein with a marginal folding barrier. We chose the small α-helical protein PDD, which folds-unfolds in microseconds crossing a free energy barrier previously estimated as ~1 RT. The fits of the PDD thermogram to the various models and assumptions produce FES with a consistently small free energy barrier separating the folded and unfolded ensembles. However, the fits vary in quality as well as in the estimated barrier. Applying Bayesian probabilistic analysis we rank the fit performance using a statistically rigorous criterion that leads to a global estimate of the folding barrier and its precision, which for PDD is 1.3 ± 0.4 kJ mol(-1). This result confirms that PDD folds over a minor barrier consistent with the downhill folding regime. We have further validated the multi-model Bayesian approach through the analysis of two additional protein systems: gpW, a midsize single-domain with α + β topology that also folds in microseconds and has been previously catalogued as a downhill folder, and α-spectrin SH3, a domain of similar size but with a β-barrel fold, slow-folding kinetics and two-state-like thermodynamics. From a general viewpoint, the Bayesian analysis developed here results in a statistically robust, virtually model-independent, method to estimate the thermodynamic free-energy barriers to protein folding from DSC thermograms. Our method appears to be sufficiently accurate to consistently detect small differences in the barrier height, and thus opens up the possibility of characterizing experimentally the changes in thermodynamic folding barriers induced by single-point mutations on proteins within the downhill regime.

The native ensemble and folding of a protein molten-globule: functional consequence of downhill folding

The continually emerging functional significance of intrinsic disorder and conformational flexibility in proteins has challenged the long-standing dogma of a well-defined structure contributing to a specific function. Molten-globular states, a class of proteins with significant secondary-structure but a fluid hydrophobic core, is one such example. They have however been difficult to characterize due to the complexity of experimental data and lack of computational avenues. Here, we dissect the folding mechanism of the α-helical molten-globular protein NCBD from three fundamentally different approaches: statistical-mechanical variable barrier model, C(α)-based Gō-model and explicit water all-atom molecular dynamics (MD) simulations. We find that NCBD displays the characteristics of a one-state globally downhill folder but is significantly destabilized. Using simulation techniques, we generate a highly constrained but a heterogeneous native ensemble of the molten-globule for the first time that is consistent with experimental data including small angle X-ray scattering (SAXS), circular dichroism (CD), and nuclear magnetic resonance (NMR). The resulting native ensemble populates conformations reported in other bound-forms providing direct evidence to the mechanism of conformational selection for binding multiple partners in this domain. Importantly, our simulations reveal a connection between downhill folding and large conformational flexibility in this domain that has been evolutionarily selected and functionally exploited resulting in large binding promiscuity. Finally, the multimodel approach we employ here serves as a powerful methodology to study mechanisms and suggests that the thermodynamic features of molten-globules fall within the array of folding mechanisms available to small single-domain proteins.

Quantitative prediction of protein folding behaviors from a simple statistical model

The statistical nature of the protein folding process requires the use of equally detailed yet simple models that lend themselves to characterize experiments. One such model is the Wako-Saitô-Muñoz-Eaton model, that we extend here to include solvation effects (WSME-S), introduced via empirical terms. We employ the novel version to analyze the folding of two proteins, gpW and SH3, that have similar size and thermodynamic stability but with the former folding 3 orders of magnitude faster than SH3. A quantitative analysis reveals that gpW presents at most marginal barriers, in contrast to SH3 that folds following a simple two-state approximation. We reproduce the observed experimental differences in melting temperature in gpW as seen by different experimental spectroscopic probes and the shape of the rate-temperature plot. In parallel, we predict the folding complexity expected in gpW from the analysis of both the residue-level thermodynamics and kinetics. SH3 serves as a stringent control with neither folding complexity nor dispersion in melting temperatures being observed. The extended model presented here serves as an ideal tool not only to characterize folding data but also to make experimentally testable predictions.