Context-dependent contributions of backbone hydrogen bonding to beta-sheet folding energetics.

Backbone hydrogen bonds (H-bonds) are prominent features of protein structures; however, their role in protein folding remains controversial because they cannot be selectively perturbed by traditional methods of protein mutagenesis. Here we have assessed the contribution of backbone H-bonds to the folding kinetics and thermodynamics of the PIN WW domain, a small beta-sheet protein, by individually replacing its backbone amides with esters. Amide-to-ester mutations site-specifically perturb backbone H-bonds in two ways: a H-bond donor is eliminated by replacing an amide NH with an ester oxygen, and a H-bond acceptor is weakened by replacing an amide carbonyl with an ester carbonyl. We perturbed the 11 backbone H-bonds of the PIN WW domain by synthesizing 19 amide-to-ester mutants. Thermodynamic studies on these variants show that the protein is most destabilized when H-bonds that are enveloped by a hydrophobic cluster are perturbed. Kinetic studies indicate that native-like secondary structure forms in one of the protein's loops in the folding transition state, but the backbone is less ordered elsewhere in the sequence. Collectively, our results provide an unusually detailed picture of the folding of a beta-sheet protein.

Beta-sheet folding mechanisms from perturbation energetics.

Amide backbone and sidechain mutagenesis data can be used in combination with kinetic and thermodynamic measurements to understand the energetic contributions of backbone hydrogen bonding and the hydrophobic effect to the acquisition of beta-sheet structure. For example, it has been revealed that loop 1 of the WW domain forms in the transition state, consistent with the emerging theme that reverse turn formation is rate limiting in beta-sheet folding. A distinct subset of WW domain residues principally influences thermodynamic stability by forming hydrogen bonds and hydrophobic interactions that stabilize the native state. Energetic data and sequence mining reveal that only a small subset of the molecular information contained in sequences or observed in high-resolution structures is required to generate folded functional beta-sheets, consistent with evolutionary robustness.

Influence of hPin1 WW N-terminal domain boundaries on function, protein stability, and folding.

An N-terminally truncated and cooperatively folded version (residues 6-39) of the human Pin1 WW domain (hPin1 WW hereafter) has served as an excellent model system for understanding triple-stranded beta-sheet folding energetics. Here we report that the negatively charged N-terminal sequence (Met1-Ala-Asp-Glu-Glu5) previously deleted, and which is not conserved in highly homologous WW domain family members from yeast or certain fungi, significantly increases the stability of hPin1 WW (approximately 4 kJ mol(-1) at 65 degrees C), in the context of the 1-39 sequence based on equilibrium measurements. N-terminal truncations and mutations in conjunction with a double mutant cycle analysis and a recently published high-resolution X-ray structure of the hPin1 cis/trans-isomerase suggest that the increase in stability is due to an energetically favorable ionic interaction between the negatively charged side chains in the N terminus of full-length hPin1 WW and the positively charged epsilon-ammonium group of residue Lys13 in beta-strand 1. Our data therefore suggest that the ionic interaction between Lys13 and the charged N terminus is the optimal solution for enhanced stability without compromising function, as ascertained by ligand binding studies. Kinetic laser temperature-jump relaxation studies reveal that this stabilizing interaction has not formed to a significant extent in the folding transition state at near physiological temperature, suggesting a differential contribution of the negatively charged N-terminal sequence to protein stability and folding rate. As neither the N-terminal sequence nor Lys13 are highly conserved among WW domains, our data further suggest that caution must be exercised when selecting domain boundaries for WW domains for structural, functional, or thermodynamic studies.

High-Resolution Mapping of the Folding Transition State of a WW Domain.

Fast-folding WW domains are among the best-characterized systems for comparing experiments and simulations of protein folding. Recent microsecond-resolution experiments and long duration (totaling milliseconds) single-trajectory modeling have shown that even mechanistic changes in folding kinetics due to mutation can now be analyzed. Thus, a comprehensive set of experimental data would be helpful to benchmark the predictions made by simulations. Here, we use T-jump relaxation in conjunction with protein engineering and report mutational Φ-values (Φ(M)) as indicators for folding transition-state structure of 65 side chain, 7 backbone hydrogen bond, and 6 deletion and /or insertion mutants within loop 1 of the 34-residue hPin1 WW domain. Forty-five cross-validated consensus mutants could be identified that provide structural constraints for transition-state structure within all substructures of the WW domain fold (hydrophobic core, loop 1, loop 2, ß-sheet). We probe the robustness of the two hydrophobic clusters in the folding transition state, discuss how local backbone disorder in the native-state can lead to non-classical Φ(M)-values (Φ(M) > 1) in the rate-determining loop 1 substructure, and conclusively identify mutations and positions along the sequence that perturb the folding mechanism from loop 1-limited toward loop 2-limited folding.

Engineering a beta-sheet protein toward the folding speed limit.

Recent experimental studies have shown that alpha-helical proteins can approach the folding "speed limit", where folding switches from an activated to a downhill process in free energy. beta-sheet proteins are generally thought to fold more slowly than helix bundles. However, based on studies of hairpins, folding should still be able to approach the microsecond time scale. Here we demonstrate how the hPin1 WW domain, a triple-stranded beta-sheet protein with a sharp thermodynamic melting transition, can be engineered toward the folding "speed limit" without a significant loss in thermal denaturation cooperativity.

Understanding the mechanism of beta-sheet folding from a chemical and biological perspective.

Perturbing the structure of the Pin1 WW domain, a 34-residue protein comprised of three beta-strands and two intervening loops has provided significant insight into the structural and energetic basis of beta-sheet folding. We will review our current perspective on how structure acquisition is influenced by the sequence, which determines local conformational propensities and mediates the hydrophobic effect, hydrogen bonding, and analogous intramolecular interactions. We have utilized both traditional site-directed mutagenesis and backbone mutagenesis approaches to alter the primary structure of this beta-sheet protein. Traditional site-directed mutagenesis experiments are excellent for altering side-chain structure, whereas amide-to-ester backbone mutagenesis experiments modify backbone-backbone hydrogen bonding capacity. The transition state structure associated with the folding of the Pin1 WW domain features a partially H-bonded, near-native reverse turn secondary structure in loop 1 that has little influence on thermodynamic stability. The thermodynamic stability of the Pin1 WW domain is largely determined by the formation of a small hydrophobic core and by the formation of desolvated backbone-backbone H-bonds enveloped by this hydrophobic core. Loop 1 engineering to the consensus five-residue beta-bulge-turn found in most WW domains or a four-residue beta-turn found in most beta-hairpins accelerates folding substantially relative to the six-residue turn found in the wild type Pin1 WW domain. Furthermore, the more efficient five- and four-residue reverse turns now contribute to the stability of the three-stranded beta-sheet. These insights have allowed the design of Pin1 WW domains that fold at rates that approach the theoretical speed limit of folding.

Structure-function-folding relationship in a WW domain.

Protein folding barriers result from a combination of factors including unavoidable energetic frustration from nonnative interactions, natural variation and selection of the amino acid sequence for function, and/or selection pressure against aggregation. The rate-limiting step for human Pin1 WW domain folding is the formation of the loop 1 substructure. The native conformation of this six-residue loop positions side chains that are important for mediating protein-protein interactions through the binding of Pro-rich sequences. Replacement of the wild-type loop 1 primary structure by shorter sequences with a high propensity to fold into a type-I' beta-turn conformation or the statistically preferred type-I G1 bulge conformation accelerates WW domain folding by almost an order of magnitude and increases thermodynamic stability. However, loop engineering to optimize folding energetics has a significant downside: it effectively eliminates WW domain function according to ligand-binding studies. The energetic contribution of loop 1 to ligand binding appears to have evolved at the expense of fast folding and additional protein stability. Thus, the two-state barrier exhibited by the wild-type human Pin1 WW domain principally results from functional requirements, rather than from physical constraints inherent to even the most efficient loop formation process.

Absolute comparison of simulated and experimental protein-folding dynamics.

Protein folding is difficult to simulate with classical molecular dynamics. Secondary structure motifs such as alpha-helices and beta-hairpins can form in 0.1-10 micros (ref. 1), whereas small proteins have been shown to fold completely in tens of microseconds. The longest folding simulation to date is a single 1- micro s simulation of the villin headpiece; however, such single runs may miss many features of the folding process as it is a heterogeneous reaction involving an ensemble of transition states. Here, we have used a distributed computing implementation to produce tens of thousands of 5-20-ns trajectories (700 micros) to simulate mutants of the designed mini-protein BBA5. The fast relaxation dynamics these predict were compared with the results of laser temperature-jump experiments. Our computational predictions are in excellent agreement with the experimentally determined mean folding times and equilibrium constants. The rapid folding of BBA5 is due to the swift formation of secondary structure. The convergence of experimentally and computationally accessible timescales will allow the comparison of absolute quantities characterizing in vitro and in silico (computed) protein folding.

Tuning the free-energy landscape of a WW domain by temperature, mutation, and truncation.

The equilibrium unfolding of the Formin binding protein 28 (FBP) WW domain, a stable three-stranded beta-sheet protein, can be described as reversible apparent two-state folding. Kinetics studied by laser temperature jump reveal a third state at temperatures below the midpoint of unfolding. The FBP free-energy surface can be tuned between three-state and two-state kinetics by changing the temperature, by truncation of the C terminus, or by selected point mutations. FBP WW domain is the smallest three-state folder studied to date and the only one that can be freely tuned between three-state and apparent two-state folding by several methods (temperature, truncation, and mutation). Its small size (28-37 residues), the availability of a quantitative reaction coordinate (phi(T)), the fast folding time scale (10s of micros), and the tunability of the folding routes by small temperature or sequence changes make this system the ideal prototype for studying more subtle features of the folding free-energy landscape by simulations or analytical theory.