Folding specificity induced by loop stiffness

Collapse kinetics and chevron plots from simulations of denaturant-dependent folding of globular proteins

Quantitative description of how proteins fold under experimental conditions remains a challenging problem. Experiments often use urea and guanidinium chloride to study folding whereas the natural variable in simulations is temperature. To bridge the gap, we use the molecular transfer model that combines measured denaturant-dependent transfer free energies for the peptide group and amino acid residues, and a coarse-grained C(α)-side chain model for polypeptide chains to simulate the folding of src SH(3) domain. Stability of the native state decreases linearly as [C] (the concentration of guanidinium chloride) increases with the slope, m, that is in excellent agreement with experiments. Remarkably, the calculated folding rate at [C] = 0 is only 16-fold larger than the measured value. Most importantly ln k(obs) (k(obs) is the sum of folding and unfolding rates) as a function of [C] has the characteristic V (chevron) shape. In every folding trajectory, the times for reaching the native state, interactions stabilizing all the substructures, and global collapse coincide. The value of (m(f) is the slope of the folding arm of the chevron plot) is identical to the fraction of buried solvent accessible surface area in the structures of the transition state ensemble. In the dominant transition state, which does not vary significantly at low [C], the core of the protein and certain loops are structured. Besides solving the long-standing problem of computing the chevron plot, our work lays the foundation for incorporating denaturant effects in a physically transparent manner either in all-atom or coarse-grained simulations.

Constraining local structure can speed up folding by promoting structural polarization of the folding pathway

The pathway which proteins take to fold can be influenced from the earliest events of structure formation. In this light, it was both predicted and confirmed that increasing the stiffness of a beta hairpin turn decreased the size of the transition state ensemble (TSE), while increasing the folding rate. Thus, there appears to be a relationship between conformationally restricting the TSE and increasing the folding rate, at least for beta hairpin turns. In this study, we hypothesize that the enormous sampling necessary to fold even two-state folding proteins in silico could be reduced if local structure constraints were used to restrict structural heterogeneity by polarizing folding pathways or forcing folding into preferred routes. Using a Gō model, we fold Chymotrypsin Inhibitor 2 (CI-2) and the src SH3 domain after constraining local sequence windows to their native structure by rigid body dynamics (RBD). Trajectories were monitored for any changes to the folding pathway and differences in the kinetics compared with unconstrained simulations. Constraining local structure decreases folding time two-fold for 41% of src SH3 windows and 45% of CI-2 windows. For both proteins, folding times are never significantly increased after constraining any window. Structural polarization of the folding pathway appears to explain these rate increases. Folding rate enhancements are consistent with the goal to reduce sampling time necessary to reach native structures during folding simulations. As anticipated, not all constrained windows showed an equal decrease in folding time. We conclude by analyzing these differences and explain why RBD may be the preferred way to constrain structure.

Theoretical perspectives on protein folding

Understanding how monomeric proteins fold under in vitro conditions is crucial to describing their functions in the cellular context. Significant advances in theory and experiments have resulted in a conceptual framework for describing the folding mechanisms of globular proteins. The sizes of proteins in the denatured and folded states, cooperativity of the folding transition, dispersions in the melting temperatures at the residue level, and timescales of folding are, to a large extent, determined by N, the number of residues. The intricate details of folding as a function of denaturant concentration can be predicted by using a novel coarse-grained molecular transfer model. By watching one molecule fold at a time, using single-molecule methods, investigators have established the validity of the theoretically anticipated heterogeneity in the folding routes and the N-dependent timescales for the three stages in the approach to the native state. Despite the successes of theory, of which only a few examples are documented here, we conclude that much remains to be done to solve the protein folding problem in the broadest sense.

Roles of beta-turns in protein folding: from peptide models to protein engineering

Reverse turns are a major class of protein secondary structure; they represent sites of chain reversal and thus sites where the globular character of a protein is created. It has been speculated for many years that turns may nucleate the formation of structure in protein folding, as their propensity to occur will favor the approximation of their flanking regions and their general tendency to be hydrophilic will favor their disposition at the solvent-accessible surface. Reverse turns are local features, and it is therefore not surprising that their structural properties have been extensively studied using peptide models. In this article, we review research on peptide models of turns to test the hypothesis that the propensities of turns to form in short peptides will relate to the roles of corresponding sequences in protein folding. Turns with significant stability as isolated entities should actively promote the folding of a protein, and by contrast, turn sequences that merely allow the chain to adopt conformations required for chain reversal are predicted to be passive in the folding mechanism. We discuss results of protein engineering studies of the roles of turn residues in folding mechanisms. Factors that correlate with the importance of turns in folding indeed include their intrinsic stability, as well as their topological context and their participation in hydrophobic networks within the protein's structure.

The amyloid stretch hypothesis: recruiting proteins toward the dark side

A detailed understanding of the molecular events underlying the conversion and self-association of normally soluble proteins into amyloid fibrils is fundamental to the identification of therapeutic strategies to prevent or cure amyloid-related disorders. Recent investigations indicate that amyloid fibril formation is not just a general property of the polypeptide backbone depending on external factors, but that it is strongly modulated by amino acid side chains. Here, we propose and address the validation of the premise that the amyloidogenicity of a protein is indeed localized in short protein stretches (amyloid stretch hypothesis). We demonstrate that the conversion of a soluble nonamyloidogenic protein into an amyloidogenic prone molecule can be triggered by a nondestabilizing six-residue amyloidogenic insertion in a particular structural environment. Interestingly enough, although the inserted amyloid sequences clearly cause the process, the protease-resistant core of the fiber also includes short adjacent sequences from the otherwise soluble globular domain. Thus, short amyloid stretches accessible for intermolecular interactions trigger the self-assembly reaction and pull the rest of the protein into the fibrillar aggregate. The reliable identification of such amyloidogenic stretches in proteins opens the possibility of using them as targets for the inhibition of the amyloid fibril formation process.

Nucleation and the transition state of the SH3 domain

We present a verified computational model of the SH3 domain transition state (TS) ensemble. This model was built for three separate SH3 domains using experimental phi-values as structural constraints in all-atom protein folding simulations. While averaging over all conformations incorrectly considers non-TS conformations as transition states, quantifying structures as pre-TS, TS, and post-TS by measurement of their transmission coefficient ("probability to fold", or p(fold)) allows for rigorous conclusions regarding the structure of the folding nucleus and a full mechanistic analysis of the folding process. Through analysis of the TS, we observe a highly polarized nucleus in which many residues are solvent-exposed. Mechanistic analysis suggests the hydrophobic core forms largely after an early nucleation step. SH3 presents an ideal system for studying the nucleation-condensation mechanism and highlights the synergistic relationship between experiment and simulation in the study of protein folding.