Thermodynamics, kinetics, and salt dependence of folding of YopM, a large leucine-rich repeat protein

Small globular proteins have many contacts between residues that are distant in primary sequence. These contacts create a complex network between sequence-distant segments of secondary structure, which may be expected to promote the cooperative folding of globular proteins. Although repeat proteins, which are composed of tandem modular units, lack sequence-distant contacts, several of considerable length have been shown to undergo cooperative two-state folding. To explore the limits of cooperativity in repeat proteins, we have studied the unfolding of YopM, a leucine-rich repeat (LRR) protein of over 400 residues. Despite its large size and modular architecture (15 repeats), YopM equilibrium unfolding is highly cooperative, and shows a very strong dependence on the concentration of urea. In contrast, kinetic studies of YopM folding indicate a mechanism that includes one or more transient intermediates. The urea dependence of the folding and unfolding rates suggests a relatively small transition state ensemble. As with the urea dependence, we have found an extreme dependence of the free energy of unfolding on the concentration of salt. This salt dependence likely results from general screening of a large number of unfavorable columbic interactions in the folded state, rather than from specific cation binding.     

Folding with downhill behavior and low cooperativity of proteins

The downhill folding observed experimentally for a small protein BBL is studied using off-lattice Gō-like model. Our simulations show that the downhill folding has low cooperativity and is barrierless, which is consistent with the experimental findings. As an example of comparison in detail, the two-state folding behavior of proteins, for example, protein CI2, is also simulated. By observing the formation of contacts between the residues for these two proteins, it is found that the physical origin of the downhill folding is due to the deficiency of nonlocal contacts which determine the folding cooperatively. From a statistics on contacts of the native structures of 17 well-studied proteins and the calculation of their cooperativity factors kappa2 based on folding simulations, a strong correlation between the number of nonlocal contacts per residue NN and the factors kappa2 is obtained. Protein BBL with a value of NN = 0.73 has the lowest cooperativity factor kappa2 = 0.34 among all 17 proteins. A crossover around NNc approximately 0.9 could be defined to separate the two-state folders and the downhill folder roughly. A protein would behave downhill folding when its NN = NNc. For proteins with their NN values are about (or slightly larger than) NNc, the folding behaves with low cooperativity and the barriers are small, showing a weak two-state behavior or a downhill-like behavior. Furthermore, simulations on mutants of a two-state folder show that a mutant becomes a downhill folder when its NN is reduced to a value smaller than NNc. These could enable us to identify the downhill folding or the cooperative two-state folding behavior solely from the native structures of proteins.     

On the relation between native geometry and conformational plasticity

In protein folding the term plasticity refers to the number of alternative folding pathways encountered in response to free energy perturbations such as those induced by mutation. Here we explore the relation between folding plasticity and a gross, generic feature of the native geometry, namely, the relative number of local and non-local native contacts. The results from our study, which is based on Monte Carlo simulations of simple lattice proteins, show that folding to a structure that is rich in local contacts is considerably more plastic than folding to a native geometry characterized by having a very large number of long-range contacts (i.e., contacts between amino acids that are separated by more than 12 units of backbone distance). The smaller folding plasticity of native geometries is probably a direct consequence of their higher folding cooperativity that renders the folding reaction more robust against single- and multiple-point mutations.     

Cooperativity and the origins of rapid, single-exponential kinetics in protein folding

The folding of naturally occurring, single-domain proteins is usually well described as a simple, single-exponential process lacking significant trapped states. Here we further explore the hypothesis that the smooth energy landscape this implies, and the rapid kinetics it engenders, arises due to the extraordinary thermodynamic cooperativity of protein folding. Studying Miyazawa-Jernigan lattice polymers, we find that, even under conditions where the folding energy landscape is relatively optimized (designed sequences folding at their temperature of maximum folding rate), the folding of protein-like heteropolymers is accelerated when their thermodynamic cooperativity is enhanced by enhancing the nonadditivity of their energy potentials. At lower temperatures, where kinetic traps presumably play a more significant role in defining folding rates, we observe still greater cooperativity-induced acceleration. Consistent with these observations, we find that the folding kinetics of our computational models more closely approximates single-exponential behavior as their cooperativity approaches optimal levels. These observations suggest that the rapid folding of naturally occurring proteins is, in part, a consequence of their remarkably cooperative folding.     

Repeat-protein folding: new insights into origins of cooperativity, stability, and topology

Although our understanding of globular protein folding continues to advance, the irregular tertiary structures and high cooperativity of globular proteins complicates energetic dissection. Recently, proteins with regular, repetitive tertiary structures have been identified that sidestep limitations imposed by globular protein architecture. Here we review recent studies of repeat-protein folding. These studies uniquely advance our understanding of both the energetics and kinetics of protein folding. Equilibrium studies provide detailed maps of local stabilities, access to energy landscapes, insights into cooperativity, determination of nearest-neighbor interaction parameters using statistical thermodynamics, relationships between consensus sequences and repeat-protein stability. Kinetic studies provide insight into the influence of short-range topology on folding rates, the degree to which folding proceeds by parallel (versus localized) pathways, and the factors that select among multiple potential pathways. The recent application of force spectroscopy to repeat-protein unfolding is providing a unique route to test and extend many of these findings.     

The contribution of entropy, enthalpy, and hydrophobic desolvation to cooperativity in repeat-protein folding

Cooperativity is a defining feature of protein folding, but its thermodynamic and structural origins are not completely understood. By constructing consensus ankyrin repeat protein arrays that have nearly identical sequences, we quantify cooperativity by resolving stability into intrinsic and interfacial components. Heteronuclear NMR and CD spectroscopy show that these constructs adopt ankyrin repeat structures. Applying a one-dimensional Ising model to a series of constructs chosen to maximize information content in unfolding transitions, we quantify stabilities of the terminal capping repeats, and resolve the effects of denaturant into intrinsic and interfacial components. Reversible thermal denaturation resolves interfacial and intrinsic free energies into enthalpic, entropic, and heat capacity terms. Intrinsic folding is entropically disfavored, whereas interfacial interaction is entropically favored and attends a decrease in heat capacity. These results suggest that helix formation and backbone ordering occurs upon intrinsic folding, whereas hydrophobic desolvation occurs upon interfacial interaction, contributing to cooperativity.     

Inter-residue interactions in protein folding and stability

During the process of protein folding, the amino acid residues along the polypeptide chain interact with each other in a cooperative manner to form the stable native structure. The knowledge about inter-residue interactions in protein structures is very helpful to understand the mechanism of protein folding and stability. In this review, we introduce the classification of inter-residue interactions into short, medium and long range based on a simple geometric approach. The features of these interactions in different structural classes of globular and membrane proteins, and in various folds have been delineated. The development of contact potentials and the application of inter-residue contacts for predicting the structural class and secondary structures of globular proteins, solvent accessibility, fold recognition and ab initio tertiary structure prediction have been evaluated. Further, the relationship between inter-residue contacts and protein-folding rates has been highlighted. Moreover, the importance of inter-residue interactions in protein-folding kinetics and for understanding the stability of proteins has been discussed. In essence, the information gained from the studies on inter-residue interactions provides valuable insights for understanding protein folding and de novo protein design.     

Amino acid replacements can selectively affect the interaction energy of autonomous folding units in the alpha subunit of tryptophan synthase.

Amino acid replacements were made at the interface between two autonomous folding units in the alpha subunit of tryptophan synthase from Salmonella typhimurium to test their mutual interaction energy. The results of equilibrium studies of the urea-induced unfolding reaction of the wild-type and mutant proteins in which phenylalanine 22 is replaced by leucine, isoleucine, and valine can be understood in terms of a selective decrease in the interaction energy between the two folding units; the intrinsic stability of each folding unit is not significantly altered. Kinetic studies of the rate-limiting step in unfolding show that the interaction energy appears in the transition state preceding the native conformation. Comparisons of the individual effects of these nonpolar side chains show that both hydrophobic and steric effects play important roles in the interaction energy between the folding units. The implication of these results is that the high cooperativity observed in the folding of many globular proteins may be reduced by appropriate amino acid replacements.     

Influence of medium- and long-range interactions in different folding types of globular proteins

Recognition of protein fold from amino acid sequence is a challenging task. The structure and stability of proteins from different fold are mainly dictated by inter-residue interactions. In our earlier work, we have successfully used the medium- and long-range contacts for predicting the protein folding rates, discriminating globular and membrane proteins and for distinguishing protein structural classes. In this work, we analyze the role of inter-residue interactions in commonly occurring folds of globular proteins in order to understand their folding mechanisms. In the medium-range contacts, the globin fold and four-helical bundle proteins have more contacts than that of DNA-RNA fold although they all belong to all-alpha class. In long-range contacts, only the ribonuclease fold prefers 4-10 range and the other folding types prefer the range 21-30 in alpha/beta class proteins. Further, the preferred residues and residue pairs influenced by these different folds are discussed. The information about the preference of medium- and long-range contacts exhibited by the 20 amino acid residues can be effectively used to predict the folding type of each protein.     

A simple lattice model that exhibits a protein-like cooperative all-or-none folding transition

In a recent paper (D. Gront et al., Journal of Chemical Physics, Vol. 115, pp. 1569, 2001) we applied a simple combination of the Replica Exchange Monte Carlo and the Histogram methods in the computational studies of a simplified protein lattice model containing hydrophobic and polar units and sequence-dependent local stiffness. A well-defined, relatively complex Greek-key topology, ground (native) conformations was found; however, the cooperativity of the folding transition was very low. Here we describe a modified minimal model of the same Greek-key motif for which the folding transition is very cooperative and has all the features of the "all-or-none" transition typical of real globular proteins. It is demonstrated that the all-or-none transition arises from the interplay between local stiffness and properly defined tertiary interactions. The tertiary interactions are directional, mimicking the packing preferences seen in proteins. The model properties are compared with other minimal protein-like models, and we argue that the model presented here captures essential physics of protein folding (structurally well-defined protein-like native conformation and cooperative all-or-none folding transition).     

Role of medium--and long-range interactions in discriminating globular and membrane proteins.

The analysis of inter-residue interactions in protein structures provides considerable insight to understand their folding and stability. We have previously analyzed the role of medium- and long-range interactions in the folding of globular proteins. In this work, we study the distinct role of such interactions in the three-dimensional structures of membrane proteins. We observed a higher number of long-range contacts in the termini of transmembrane helical (TMH) segments, implying their role in the stabilization of helix-helix interactions. The transmembrane strand (TMS) proteins are having appreciably higher long-range contacts than that in all-beta class of globular proteins, indicating closer packing of the strands in TMS proteins. The residues in membrane spanning segments of TMH proteins have 1.3 times higher medium-range contacts than long-range contacts whereas that of TMS proteins have 14 times higher long-range contacts than medium-range contacts. Residue-wise analysis indicates that in TMH proteins, the residues Cys, Glu, Gly, Pro, Gln, Ser and Tyr have higher long-range contacts than medium-range contacts in contrast with all-alpha class of globular proteins. The charged residue pairs have higher medium-range contacts in all-alpha proteins, whereas hydrophobic residue pairs are dominant in TMH proteins. The information on the preference of residue pairs to form medium-range contacts has been successfully used to discriminate the TMH proteins from all-alpha proteins. The statistical significance of the results obtained from the present study has been verified using randomized structures of TMH and TMS protein templates.     

Transition states for folding of circular-permuted proteins

To explore the role of entropy and chain connectivity in protein folding, a particularly interesting scheme, namely, the circular permutation, has been used. Recently, experimental observations showed that there are large differences in the folding mechanisms between the wild-type proteins and their circular permutants. These differences are strongly related to the change in the intrachain connectivity. Some results obtained by molecular dynamics simulations also showed a good agreement with the experimental findings. Here, we use a topology-based free-energy functional method to study the role of the chain connectivity in folding by comparing features of transition states of the wild-type proteins with those of their circular permutants. We concentrate our study on 3 small globular proteins, namely, the alpha-spectrin SH3 domain (SH3), the chymotrypsin inhibitor 2 (CI2), and the ribosomal protein S6, and obtain exciting results that are consistent with the available experimental and simulation results. A heterogeneity of the interaction energies between contacts for protein CI2 and for protein S6 is also introduced, which characterizes the strong interactions between contacts with long loops, as speculated from experiments for protein S6. The comparison between the folding nucleus of the wild-type proteins and those of their circular permutants indicates that chain connectivity affects remarkably the shapes of the energy profiles and thus the folding mechanism. Further comparisons between our theoretical calculated phi(th) values and the experimental observed phi(exp) values for the 3 proteins and their permutants show that our results are in good agreement with experimental ones and that correlations between them are high. These indicate that the free-energy functional method really provides a way to analyze the folding behavior of the circular-permuted proteins and therefore the folding mechanism of the wild-type proteins.     

Solvation effects and driving forces for protein thermodynamic and kinetic cooperativity: how adequate is native-centric topological modeling?

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ō-like contact energies, although native-centric, are by themselves insufficient for protein-like calorimetric two-state cooperativity. Further native biases by local conformational preferences are necessary for protein-like 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 simple 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. The free energy profiles of these models are found to be sensitive to the choice of native contacts and the presumed spatial ranges of the contact interactions. Motivated by explicit-water considerations, we explore recent treatments of solvent granularity that incorporate desolvation free energy barriers into effective implicit-solvent intraprotein interactions. This additional feature reduces both folding and unfolding rates vis-à-vis that of the corresponding models without desolvation barriers, but the kinetics remain non-two-state. Taken together, our observations suggest that interaction mechanisms more intricate than simple Gō-like constructs and pairwise additive solvation-like contributions are needed to rationalize some of the most basic generic protein properties. Therefore, as experimental constraints on protein chain models, requiring a consistent account of protein-like thermodynamic and kinetic cooperativity can be more stringent and productive for some applications than simply requiring a model heteropolymer to fold to a target structure.     

Folding thermodynamics and kinetics of the leucine-rich repeat domain of the virulence factor Internalin B

Although the folding of alpha-helical repeat proteins has been well characterized, much less is known about the folding of repeat proteins containing beta-sheets. Here we investigate the folding thermodynamics and kinetics of the leucine-rich repeat (LRR) domain of Internalin B (InlB), an extracellular virulence factor from the bacterium Lysteria monocytogenes. This domain contains seven tandem leucine-rich repeats, of which each contribute a single beta-strand that forms a continuous beta-sheet with neighboring repeats, and an N-terminal alpha-helical capping motif. Despite its modular structure, InlB folds in an equilibrium two-state manner, as reflected by the identical thermodynamic parameters obtained by monitoring its sigmoidal urea-induced unfolding transition by different spectroscopic probes. Although equilibrium two-state folding is common in alpha-helical repeat proteins, to date, InlB is the only beta-sheet-containing repeat protein for which this behavior is observed. Surprisingly, unlike other repeat proteins exhibiting equilibrium two-state folding, InlB also folds by a simple two-state kinetic mechanism lacking intermediates, aside from the effects of prolyl isomerization on the denatured state. However, like other repeat proteins, InlB also folds significantly more slowly than expected from contact order. When plotted against urea, the rate constants for the fast refolding and single unfolding phases constitute a linear chevron that, when fitted with a kinetic two-state model, yields thermodynamic parameters matching those observed for equilibrium folding. Based on these kinetic parameters, the transition state is estimated to comprise 40% of the total surface area buried upon folding, indicating that a large fraction of the native contacts are formed in the rate-limiting step to folding.     

Direct Observation of Parallel Folding Pathways Revealed Using a Symmetric Repeat Protein System

Although progress has been made to determine the native fold of a polypeptide from its primary structure, the diversity of pathways that connect the unfolded and folded states has not been adequately explored. Theoretical and computational studies predict that proteins fold through parallel pathways on funneled energy landscapes, although experimental detection of pathway diversity has been challenging. Here, we exploit the high translational symmetry and the direct length variation afforded by linear repeat proteins to directly detect folding through parallel pathways. By comparing folding rates of consensus ankyrin repeat proteins (CARPs), we find a clear increase in folding rates with increasing size and repeat number, although the size of the transition states (estimated from denaturant sensitivity) remains unchanged. The increase in folding rate with chain length, as opposed to a decrease expected from typical models for globular proteins, is a clear demonstration of parallel pathways. This conclusion is not dependent on extensive curve-fitting or structural perturbation of protein structure. By globally fitting a simple parallel-Ising pathway model, we have directly measured nucleation and propagation rates in protein folding, and have quantified the fluxes along each path, providing a detailed energy landscape for folding. This finding of parallel pathways differs from results from kinetic studies of repeat-proteins composed of sequence-variable repeats, where modest repeat-to-repeat energy variation coalesces folding into a single, dominant channel. Thus, for globular proteins, which have much higher variation in local structure and topology, parallel pathways are expected to be the exception rather than the rule.     

Cooperativity,smooth energy landscapes and the origins of topology-dependent protein folding rates

The relative folding rates of simple, single-domain proteins, proteins whose folding energy landscapes are smooth, are highly dispersed and strongly correlated with native-state topology. In contrast, the relative folding rates of small, Gō-potential lattice polymers, which also exhibit smooth energy landscapes, are poorly dispersed and insignificantly correlated with native-state topology. Here, we investigate this discrepancy in light of a recent, quantitative theory of two-state folding kinetics, the topomer search model. This model stipulates that the topology-dependence of two-state folding rates is a direct consequence of the extraordinarily cooperative equilibrium folding of simple proteins. We demonstrate that traditional Gō polymers lack the extreme cooperativity that characterizes the folding of naturally occurring, two-state proteins and confirm that the folding rates of a diverse set of Gō 27-mers are poorly dispersed and effectively uncorrelated with native state topology. Upon modestly increasing the cooperativity of the Gō-potential, however, significantly increased dispersion and strongly topology-dependent kinetics are observed. These results support previous arguments that the cooperative folding of simple, single-domain proteins gives rise to their topology-dependent folding rates. We speculate that this cooperativity, and thus, indirectly, the topology-rate relationship, may have arisen in order to generate the smooth energetic landscapes upon which rapid folding can occur.     

Probing possible downhill folding: native contact topology likely places a significant constraint on the folding cooperativity of proteins with approximately 40 residues

Experiments point to appreciable variations in folding cooperativity among natural proteins with approximately 40 residues, indicating that the behaviors of these proteins are valuable for delineating the contributing factors to cooperative folding. To explore the role of native topology in a protein's propensity to fold cooperatively and how native topology might constrain the degree of cooperativity achievable by a given set of physical interactions, we compared folding/unfolding kinetics simulated using three classes of native-centric C^α chain models with different interaction schemes. The approach was applied to two homologous 45-residue fragments from the peripheral subunit-binding domain family and a 39-residue fragment of the N-terminal domain of ribosomal protein L9. Free-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. Models that account for water-expulsion (desolvation) barriers and models with many-body (nonadditive) interactions generally entail higher degrees of folding cooperativity indicated by more linear model chevron plots, but the rank ordering of cooperativity remains unchanged. A robust, experimentally valid rank ordering of model folding cooperativity independent of the multiple native-centric interaction schemes tested here argues that native topology places significant constraints on how cooperatively a protein can fold.     

Intermediates in the folding equilibrium of repeat proteins from the TPR family

In recent decades, advances in computational methods and experimental biophysical techniques have improved our understanding of protein folding. Although some of these advances have been remarkable, the structural variability of globular proteins usually encountered makes it difficult to extract general features of their folding processes. To overcome this difficulty, experimental and computational studies of the folding of repeat (or modular) proteins are of interest. Because their native structures can be described as linear arrays of the same, repeated, supersecondary structure unit, it is possible to seek  a possibly independent behavior of the different modules without taking into account the intrinsic stability associated with different secondary structure motifs. In this work we have used a Monte Carlo-based simulation to study the folding equilibrium of four repeat proteins belonging to the tetratricopeptide repeat family. Our studies provide new insights into their energy profiles, enabling investigation about the existence of intermediate states and their relative stabilities. We have also performed structural analyses to describe the structure of these intermediates, going through the vast number of conformations obtained from the simulations. In this way, we have tried to identify the regions of each protein in which the modular structure yields a different behavior and, more specifically, regions of the proteins that can stay folded when the rest of the chain has been thermally denatured.     

Folding landscapes of ankyrin repeat proteins: experiments meet theory

Nearly 6% of eukaryotic protein sequences contain ankyrin repeat (AR) domains, which consist of several repeats and often function in binding. AR proteins show highly cooperative folding despite a lack of long-range contacts. Both theory and experiment converge to explain that formation of the interface between elements is more favorable than formation of any individual repeat unit. IkappaBalpha and Notch both undergo partial folding upon binding perhaps influencing the binding free energy. The simple architecture, combined with identification of consensus residues that are important for stability, has enabled systematic perturbation of the energy landscape by single point mutations that affect stability or by addition of consensus repeats. The folding energy landscapes appear highly plastic, with small perturbations re-routing folding pathways.