Exploiting the downhill folding regime via experiment

Protein folding rate evolution upon mutations

Despite the spectacular success of cutting-edge protein fold prediction methods, many critical questions remain unanswered, including why proteins can reach their native state in a biologically reasonable time. A satisfactory answer to this simple question could shed light on the slowest folding rate of proteins as well as how mutations-amino-acid substitutions and/or post-translational modifications-might affect it. Preliminary results indicate that (i) Anfinsen's dogma validity ensures that proteins reach their native state on a reasonable timescale regardless of their sequence or length, and (ii) it is feasible to determine the evolution of protein folding rates without accounting for epistasis effects or the mutational trajectories between the starting and target sequences. These results have direct implications for evolutionary biology because they lay the groundwork for a better understanding of why, and to what extent, mutations-a crucial element of evolution and a factor influencing it-affect protein evolvability. Furthermore, they may spur significant progress in our efforts to solve crucial structural biology problems, such as how a sequence encodes its folding.

The Molten Globule State of a Globular Protein in a Cell Is More or Less Frequent Case Rather than an Exception

Quite a long time ago, Oleg B. Ptitsyn put forward a hypothesis about the possible functional significance of the molten globule (MG) state for the functioning of proteins. MG is an intermediate between the unfolded and the native state of a protein. Its experimental detection and investigation in a cell are extremely difficult. In the last decades, intensive studies have demonstrated that the MG-like state of some globular proteins arises from either their modifications or interactions with protein partners or other cell components. This review summarizes such reports. In many cases, MG was evidenced to be functionally important. Thus, the MG state is quite common for functional cellular proteins. This supports Ptitsyn’s hypothesis that some globular proteins may switch between two active states, rigid (N) and soft (MG), to work in solution or interact with partners.

Predicting and Simulating Mutational Effects on Protein Folding Kinetics

Mutational perturbations of protein structures, i.e., phi-value analysis, are commonly employed to probe the extent of involvement of a particular residue in the rate-determining step(s) of folding. This generally involves the measurement of folding thermodynamic parameters and kinetic rate constants for the wild-type and mutant proteins. While computational approaches have been reasonably successful in understanding and predicting the effect of mutations on folding thermodynamics, it has been challenging to explore the same on kinetics due to confounding structural, energetic, and dynamic factors. Accordingly, the frequent observation of fractional phi-values (mean of ~0.3) has resisted a precise and consistent interpretation. Here, we describe how to construct, parameterize, and employ a simple one-dimensional free energy surface model that is grounded in the basic tenets of the energy landscape theory to predict and simulate the effect of mutations on folding kinetics. As a proof of principle, we simulate one-dimensional free energy profiles of 806 mutations from 24 different proteins employing just the experimental destabilization as input, reproduce the relative unfolding activation free energies with a correlation of 0.91, and show that the mean phi-value of 0.3 essentially corresponds to the extent of stabilization energy gained at the barrier top while folding.

The change of conditions does not affect Ros87 downhill folding mechanism

Downhill folding has been defined as a unique thermodynamic process involving a conformations ensemble that progressively loses structure with the decrease of protein stability. Downhill folders are estimated to be rather rare in nature as they miss an energetically substantial folding barrier that can protect against aggregation and proteolysis. We have previously demonstrated that the prokaryotic zinc finger protein Ros87 shows a bipartite folding/unfolding process in which a metal binding intermediate converts to the native structure through a delicate barrier-less downhill transition. Significant variation in folding scenarios can be detected within protein families with high sequence identity and very similar folds and for the same sequence by varying conditions. For this reason, we here show, by means of DSC, CD and NMR, that also in different pH and ionic strength conditions Ros87 retains its partly downhill folding scenario demonstrating that, at least in metallo-proteins, the downhill mechanism can be found under a much wider range of conditions and coupled to other different transitions. We also show that mutations of Ros87 zinc coordination sphere produces a different folding scenario demonstrating that the organization of the metal ion core is determinant in the folding process of this family of proteins.

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.

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.

Modulation of allosteric coupling by mutations: from protein dynamics and packing to altered native ensembles and function

A large body of work has gone into understanding the effect of mutations on protein structure and function. Conventional treatments have involved quantifying the change in stability, activity and relaxation rates of the mutants with respect to the wild-type protein. However, it is now becoming increasingly apparent that mutational perturbations consistently modulate the packing and dynamics of a significant fraction of protein residues, even those that are located >10–15 Å from the mutated site. Such long-range modulation of protein features can distinctly tune protein stability and the native conformational ensemble contributing to allosteric modulation of function. In this review, I summarize a series of experimental and computational observations that highlight the incredibly pliable nature of proteins and their response to mutational perturbations manifested via the intra-protein interaction network. I highlight how an intimate understanding of mutational effects could pave the way for integrating stability, folding, cooperativity and even allostery within a single physical framework.

Extracting the Hidden Distributions Underlying the Mean Transition State Structures in Protein Folding

The inherent conflict between noncovalent interactions and the large conformational entropy of the polypeptide chain forces folding reactions and their mechanisms to deviate significantly from chemical reactions. Accordingly, measures of structure in the transition state ensemble (TSE) are strongly influenced by the underlying distributions of microscopic folding pathways that are challenging to discern experimentally. Here, we present a detailed analysis of 150,000 folding transition paths of five proteins at three different thermodynamic conditions from an experimentally consistent statistical mechanical model. We find that the underlying TSE structural distributions are rarely unimodal, and the average experimental measures arise from complex underlying distributions. Unfolding pathways also exhibit subtle differences from folding counterparts due to a combination of Hammond behavior and native-state movements. Local interactions and topological complexity, to a lesser extent, are found to determine pathway heterogeneity, underscoring the importance of the balance between local and nonlocal energetics in protein folding.

Roles of conformational disorder and downhill folding in modulating protein-DNA recognition

Transcription factors are thought to efficiently search for their target DNA site via a combination of conventional 3D diffusion and 1D diffusion along the DNA molecule mediated by non-specific electrostatic interactions. This process requires the DNA-binding protein to quickly exchange between a search competent and a target recognition mode, but little is known as to how these two binding modes are encoded in the conformational properties of the protein. Here, we investigate this issue on the engrailed homeodomain (EngHD), a DNA-binding domain that folds ultrafast and exhibits a complex conformational behavior consistent with the downhill folding scenario. We explore the interplay between folding and DNA recognition using a coarse-grained computational model that allows us to manipulate the folding properties of the protein and monitor its non-specific and specific binding to DNA. We find that conformational disorder increases the search efficiency of EngHD by promoting a fast gliding search mode in addition to sliding. When gliding, EngHD remains loosely bound to DNA moving linearly along its length. A partially disordered EngHD also binds more dynamically to the target site, reducing the half-life of the specific complex via a spring-loaded mechanism. These findings apply to all conditions leading to partial disorder. However, we also find that at physiologically relevant temperatures EngHD is well folded and can only obtain the conformational flexibility required to accelerate 1D diffusion when it folds/unfolds within the downhill scenario (crossing a marginal free energy barrier). In addition, the conformational flexibility of native downhill EngHD enables its fast reconfiguration to lock into the specific binding site upon arrival, thereby affording finer control of the on- and off-rates of the specific complex. Our results provide key mechanistic insights into how DNA-binding domains optimize specific DNA recognition through the control of their conformational dynamics and folding mechanism.

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.

Evidence for the sequential folding mechanism in RNase H from an ensemble-based model

The number of distinct protein folding pathways starting from an unfolded ensemble, and hence, the folding mechanism is an intricate function of protein size, sequence complexity, and stability conditions. This has traditionally been a contentious issue particularly because of the ensemble nature of conventional experiments that can mask the complexity of the underlying folding landscape. Recent hydrogen-exchange experiments combined with mass spectrometry (HX-MS) reveal that the folding of RNase H proceeds in a hierarchical fashion with distinct intermediates and following a single discrete path. In our current work, we provide computational evidence for this unique folding mechanism by employing a structure-based statistical mechanical model. Upon calibrating the energetic terms of the model with equilibrium measurements, we predict multiple intermediate states in the folding of RNase H that closely resemble experimental observations. Remarkably, a simplified landscape representation adequately captures the folding complexity and predicts the possibility of a well-defined sequence of folding events. We supplement the statistical model study with both explicit solvent molecular simulations of the helical units and electrostatic calculations to provide structural and energetic insights into the early and late stages of RNase H folding that hint at the frustrated nature of its folding landscape.

How quickly can a β-hairpin fold from its transition state?

Understanding the structural nature of the free energy bottleneck(s) encountered in protein folding is essential to elucidating the underlying dynamics and mechanism. For this reason, several techniques, including Φ-value analysis, have previously been developed to infer the structural characteristics of such high free-energy or transition states. Herein we propose that one (or few) appropriately placed backbone and/or side chain cross-linkers, such as disulfides, could be used to populate a thermodynamically accessible conformational state that mimics the folding transition state. Specifically, we test this hypothesis on a model β-hairpin, Trpzip4, as its folding mechanism has been extensively studied and is well understood. Our results show that cross-linking the two β-strands near the turn region increases the folding rate by an order of magnitude, to about (500 ns)⁻¹, whereas cross-linking the termini results in a hyperstable β-hairpin that has essentially the same folding rate as the uncross-linked peptide. Taken together, these findings suggest that cross-linking is not only a useful strategy to manipulate folding free energy barriers, as shown in other studies, but also, in some cases, it can be used to stabilize a folding transition state analogue and allow for direct assessment of the folding process on the downhill side of the free energy barrier. The calculated free energy landscape of the cross-linked Trpzip4 also supports this picture. An empirical analysis further suggests, when folding of β-hairpins does not involve a significant free energy barrier, the folding time (τ) follows a power law dependence on the number of hydrogen bonds to be formed (n_H), namely, τ = τ₀n_H^α, with τ₀ = 20 ns and α = 2.3.

Isomerization- and temperature-jump-induced dynamics of a photoswitchable β-hairpin

Conformational changes in proteins and peptides can be initiated by diverse processes. This raises the question how the variation of initiation mechanisms is connected to differences in folding or unfolding processes. In this work structural dynamics of a photoswitchable β-hairpin model peptide were initiated by two different mechanisms: temperature jump (T-jump) and isomerization of a backbone element. In both experiments the structural changes were followed by time-resolved IR spectroscopy in the nanosecond to microsecond range. When the photoisomerization of the azobenzene backbone switch initiated the folding reaction, pronounced absorption changes related to folding into the hairpin structure were found with a time constant of about 16 μs. In the T-jump experiment kinetics with the same time constant were observed. For both initiation processes the reaction dynamics revealed the same strong dependence of the reaction time on temperature. The highly similar transients in the microsecond range show that the peptide dynamics induced by T-jump and isomerization are both determined by the same mechanism and exclude a downhill-folding process. Furthermore, the combination of the two techniques allows a detailed model for folding and unfolding to be presented: The isomerization-induced folding process ends in a transition-state reaction scheme, in which a high energetic barrier of 48 kJ mol(-1) separates unfolded and folded structures.

Revisiting the NMR structure of the ultrafast downhill folding protein gpW from bacteriophage λ

GpW is a 68-residue protein from bacteriophage λ that participates in virus head morphogenesis. Previous NMR studies revealed a novel α+β fold for this protein. Recent experiments have shown that gpW folds in microseconds by crossing a marginal free energy barrier (i.e., downhill folding). These features make gpW a highly desirable target for further experimental and computational folding studies. As a step in that direction, we have re-determined the high-resolution structure of gpW by multidimensional NMR on a construct that eliminates the purification tags and unstructured C-terminal tail present in the prior study. In contrast to the previous work, we have obtained a full manual assignment and calculated the structure using only unambiguous distance restraints. This new structure confirms the α+β topology, but reveals important differences in tertiary packing. Namely, the two α-helices are rotated along their main axis to form a leucine zipper. The β-hairpin is orthogonal to the helical interface rather than parallel, displaying most tertiary contacts through strand 1. There also are differences in secondary structure: longer and less curved helices and a hairpin that now shows the typical right-hand twist. Molecular dynamics simulations starting from both gpW structures, and calculations with CS-Rosetta, all converge to our gpW structure. This confirms that the original structure has strange tertiary packing and strained secondary structure. A comparison of NMR datasets suggests that the problems were mainly caused by incomplete chemical shift assignments, mistakes in NOE assignment and the inclusion of ambiguous distance restraints during the automated procedure used in the original study. The new gpW corrects these problems, providing the appropriate structural reference for future work. Furthermore, our results are a cautionary tale against the inclusion of ambiguous experimental information in the determination of protein structures.

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.

Navigating the downhill protein folding regime via structural homologues

Proteins that fold over free-energy barriers <or= 3RT are classified as downhill folders. This regime is characterized by equilibrium unfolding that is proportionally broader and more complex the lower the folding barrier. Downhill proteins are also expected to fold up in a few microseconds. However, the relationship between rate and equilibrium signatures is affected by other factors such as protein size and folding topology. Here we perform a direct comparison of the kinetics and equilibrium unfolding of two structural homologues: BBL and PDD. BBL folds-unfolds in just approximately 1 micros at 335 K and displays the equilibrium signatures expected for a protein at the bottom of the downhill folding regime. PDD, which has the same 3D structure and size, folds-unfolds approximately 8 times more slowly and, concomitantly, exhibits all the downhill equilibrium signatures to a lesser degree. Our results demonstrate that the equilibrium signatures of downhill folding are proportional to the changes in folding rate once structural and size-scaling effects are factored out. This conclusion has two important implications: (1) it confirms that the quantitative analysis of equilibrium experiments in ultrafast folding proteins does provide direct information about free-energy barriers, a result that is incompatible with the conventional view of protein folding as a highly activated process, and (2) it advocates for equilibrium-kinetic studies of homologous proteins as a powerful tool to navigate the downhill folding regime via comparative analysis. The latter should prove extremely useful for the investigation of sequence, functional, and evolutionary determinants of protein folding barriers.

The effect of electrostatics on the marginal cooperativity of an ultrafast folding protein

Proteins fold up by coordinating the different segments of their polypeptide chain through a network of weak cooperative interactions. Such cooperativity results in unfolding curves that are typically sigmoidal. However, we still do not know what factors modulate folding cooperativity or the minimal amount that ensures folding into specific three-dimensional structures. Here, we address these issues on BBL, a small helical protein that folds in microseconds via a marginally cooperative downhill process (Li, P., Oliva, F. Y., Naganathan, A. N., and Muñoz, V. (2009) Proc. Natl. Acad. Sci. USA. 106, 103-108). Particularly, we explore the effects of salt-induced screening of the electrostatic interactions in BBL at neutral pH and in acid-denatured BBL. Our results show that electrostatic screening stabilizes the native state of the neutral and protonated forms, inducing complete refolding of acid-denatured BBL. Furthermore, without net electrostatic interactions, the unfolding process becomes much less cooperative, as judged by the broadness of the equilibrium unfolding curve and the relaxation rate. Our experiments show that the marginally cooperative unfolding of BBL can still be made twice as broad while the protein retains its ability to fold into the native three-dimensional structure in microseconds. This result demonstrates experimentally that efficient folding does not require cooperativity, confirming predictions from theory and computer simulations and challenging the conventional biochemical paradigm. Furthermore, we conclude that electrostatic interactions are an important factor in determining folding cooperativity. Thus, electrostatic modulation by pH-salt and/or mutagenesis of charged residues emerges as an attractive tool for tuning folding cooperativity.

Insights into protein folding mechanisms from large scale analysis of mutational effects

Protein folding mechanisms are probed experimentally using single-point mutant perturbations. The relative effects on the folding (phi-values) and unfolding (1 - phi) rates are used to infer the detailed structure of the transition-state ensemble (TSE). Here we analyze kinetic data on > 800 mutations carried out for 24 proteins with simple kinetic behavior. We find two surprising results: (i) all mutant effects are described by the equation: DeltaDeltaG(double dagger)(unfold)=0.76DeltaDeltaG(eq) +/- 1.8kJ/mol. Therefore all data are consistent with a single phi-value (0.24) with accuracy comparable to experimental precision, suggesting that the structural information in conventional phi-values is low. (ii) phi-values change with stability, increasing in mean value and spread from native to unfolding conditions, and thus cannot be interpreted without proper normalization. We eliminate stability effects calculating the phi-values at the mutant denaturation midpoints; i.e., conditions of zero stability (phi(0)). We then show that the intrinsic variability is phi(0) = 0.36 +/- 0.11, being somewhat larger for beta-sheet-rich proteins than for alpha-helical proteins. Importantly, we discover that phi(0)-values are proportional to how many of the residues surrounding the mutated site are local in sequence. High phi(0)-values correspond to protein surface sites, which have few nonlocal neighbors, whereas core residues with many tertiary interactions produce the lowest phi(0)-values. These results suggest a general mechanism in which the TSE at zero stability is a broad conformational ensemble stabilized by local interactions and without specific tertiary interactions, reconciling phi-values with many other empirical observations.

Understanding protein folding cooperativity based on topological consideration

The folding cooperativity is an important issue of protein folding dynamics. Since the native topology plays a significant role in determining the folding behavior of proteins, we believe that it also has close relationship with the folding cooperativity. In the present work, we perform simulations on proteins Naf-BBL, QNND-BBL, CI2, and SH3 with the Gō model and compare their different folding behaviors. By analyzing the weak cooperative folding of protein Naf-BBL in detail, we found that the folding of Naf-BBL shows relatively weak thermodynamic coupling between residues, and such weak coupling is found mainly between the nonlocal native contacts. This finding complements our understandings on the source of barrierless folding of Naf-BBL and promotes us to analyze the topological origins of the poor thermodynamic coupling of Naf-BBL. Then, we further extend our analysis to other two-state and multistate proteins. Based on the considerations of the thermodynamic coupling and kinetic coupling, we conclude that the fraction of scattered native contacts, the difference in loop entropy of contacts, and the long range relative contact order are the major topological factors that influence the folding cooperativity. The combination of these three tertiary structural features shows significant correlations with the folding types of proteins. Moreover, we also discuss the topological factors related to downhill folding. Finally, the generic role of tertiary structure in determining the folding cooperativity is summarized.