The ensemble folding kinetics of the FBP28 WW domain revealed by an all-atom Monte Carlo simulation in a knowledge-based potential

In this work, we apply a detailed all‐atom model with a transferable knowledge‐based potential to study the folding kinetics of Formin‐Binding protein, FBP28, which is a canonical three‐stranded β‐sheet WW domain. Replica exchange Monte Carlo simulations starting from random coils find native‐like (Cα RMSD of 2.68 Å) lowest energy structure. We also study the folding kinetics of FBP28 WW domain by performing a large number of ab initio Monte Carlo folding simulations. Using these trajectories, we examine the order of formation of two β‐hairpins, the folding mechanism of each individual β‐hairpin, and transition state ensemble (TSE) of FBP28 WW domain and compare our results with experimental data and previous computational studies. To obtain detailed structural information on the folding dynamics viewed as an ensemble process, we perform a clustering analysis procedure based on graph theory. Further, a rigorous P_fold analysis is used to obtain representative samples of the TSEs showing good quantitative agreement between experimental and simulated Φ values. Our analysis shows that the turn structure between first and second β strands is a partially stable structural motif that gets formed before entering the TSE in FBP28 WW domain and there exist two major pathways for the folding of FBP28 WW domain, which differ in the order and mechanism of hairpin formation. Proteins 2011. © 2011 Wiley‐Liss, Inc.     

Understanding the folding and stability of a designed WW domain protein with replica exchange molecular dynamics simulations

WW domain proteins are usually regarded as simple models for understanding the folding mechanism of β-sheet. CC45 is an artificial protein that is capable of folding into the same structure as WW domain. In this article, the replica exchange molecular dynamics simulations are performed to investigate the folding mechanism of CC45. The analysis of thermal stability shows that β-hairpin 1 is more stable than β-hairpin 2 during the unfolding process. Free energy analysis shows that the unfolding of this protein substantially proceeds through solvating the smaller β-hairpin 2, followed by the unfolding of β-hairpin 1. We further propose the unfolding process of CC45 and the folding mechanism of two β-hairpins. These results are similar to the previous folding studies of formin binding protein 28 (FBP28). Compared with FBP28, it is found that CC45 has more aromatic residues in N-terminal loop, and these residues contact with C-terminal loop to form the outer hydrophobic core, which increases the stability of CC45. Knowledge about the stability and folding behaviour of CC45 may help in understanding the folding mechanisms of the β-sheet and in designing new WW domains.     

Folding, misfolding, and amyloid protofibril formation of WW domain FBP28

We study the folding mechanism of a triple beta-strand WW domain from the Formin binding protein 28 (FBP28) at atomic resolution with explicit water model using replica exchange molecular dynamics computer simulations. Extended sampling over a wide range of temperatures to obtain the free energy, enthalpy, and entropy surfaces as a function of structural coordinates has been performed. Simulations were started from different configurations covering the folded and unfolded states. In the free energy landscape a transition state is identified and its structures and -values are compared with experimental data from a homologous protein, the prolyl-isomerase Pin1 WW domain. A stable intermediate state is found to accumulate during the simulation characterized by the carboxyl-terminal beta-strand 3 having misregistered hydrogen bonds and where the structural heterogeneity is due to nonnative turn II formation. Furthermore, the aggregation behavior of the FBP28 WW domain may be related to one such misfolded structure, which has a much lower free energy of dimer formation than that of the native dimer. Based on the misfolded dimer, aggregation to form protofibril structure is discussed.     

Structural characterization of a new binding motif and a novel binding mode in group 2 WW domains

Formin homology 1 (FH1), is a long proline-rich region of formins, shown to bind to five WW containing proteins named formin binding proteins (FBPs). FH1 has several potential binding regions but only the PPLPx motif and its interaction with FBP11WW1 has been characterized structurally. To detect whether additional motifs exist in FH1, we synthesized five peptides and investigated their interaction with FBP28WW2, FBP11WW1 and FBP11WW2 domains. Peptides of sequence PTPPPLPP (positive control), PPPLIPPPP and PPLIPPPP (new motifs) interact with the domains with micromolar affinity. We observed that FBP28WW2 and FBP11WW2 behave differently from FBP11WW1 in terms of motif selection and affinity, since they prefer a doubly interrupted proline stretch of sequence PPLIPP. We determined the NMR structure of three complexes involving the FBP28WW2 domain and the three ligands. Depending on the peptide under study, the domain interacts with two proline residues accommodated in either the XP or the XP2 groove. This difference represents a one-turn displacement of the domain along the ligand sequence. To understand what drives this behavior, we performed further structural studies with the FBP11WW1 and a mutant of FBP28WW2 mimicking the XP2 groove of FBP11WW1. Our observations suggest that the nature of the XP2 groove and the balance of flexibility/rigidity around loop 1 of the domain contribute to the selection of the final ligand positioning in fully independent domains. Additionally, we analyzed the binding of a double WW domain region, FBP11WW1-2, to a long stretch of FH1 using fluorescence spectroscopy and NMR titrations. With this we show that the presence of two consecutive WW domains may also influence the selection of the binding mode, particularly if both domains can interact with consecutive motifs in the ligand. Our results represent the first observation of protein-ligand recognition where a pair of WW and two consecutive motifs in a ligand participate simultaneously.     

(Un)Folding Mechanisms of the FBP28 WW Domain in Explicit Solvent Revealed by Multiple Rare Event Simulation Methods

We report a numerical study of the (un)folding routes of the truncated FBP28 WW domain at ambient conditions using a combination of four advanced rare event molecular simulation techniques. We explore the free energy landscape of the native state, the unfolded state, and possible intermediates, with replica exchange molecular dynamics. Subsequent application of bias-exchange metadynamics yields three tentative unfolding pathways at room temperature. Using these paths to initiate a transition path sampling simulation reveals the existence of two major folding routes, differing in the formation order of the two main hairpins, and in hydrophobic side-chain interactions. Having established that the hairpin strand separation distances can act as reasonable reaction coordinates, we employ metadynamics to compute the unfolding barriers and find that the barrier with the lowest free energy corresponds with the most likely pathway found by transition path sampling. The unfolding barrier at 300 K is ∼17 k_BT ≈ 42 kJ/mol, in agreement with the experimental unfolding rate constant. This work shows that combining several powerful simulation techniques provides a more complete understanding of the kinetic mechanism of protein folding.     

The folding mechanism of a beta-sheet: the WW domain

The folding thermodynamics and kinetics of the Pin WW domain, a three-stranded antiparallel beta-sheet, have been characterized extensively. Folding and activation free energies were determined as a function of temperature for 16 mutants, which sample all strands and turns of the molecule. The mutational phi value (Phi(m)) diagram is a smooth function of sequence, indicating a prevalence of local interactions in the transition state (TS). At 37 degrees C, the diagram has a single pronounced maximum at turn 1: the rate-limiting step during folding is the formation of loop 1. In contrast, key residues for thermodynamic stability are located in the strand hydrophobic clusters, indicating that factors contributing to protein stability and folding kinetics are not correlated. The location of the TS along the entropic reaction coordinate Phi(T), obtained by temperature-tuning the kinetics, reveals that sufficiently destabilizing mutants in loop 2 or in the Leu7-Trp11-Tyr24-Pro37 hydrophobic cluster can cause a switch to a late TS. Phi(m) analysis is usually applied "perturbatively" (methyl truncation), but with Phi(T) to quantitatively assess TS shifts along a reaction coordinate, more severe mutations can be used to probe regions of the free energy surface beyond the TS.     

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.     

FBP WW domains and the Abl SH3 domain bind to a specific class of proline-rich ligands.

WW domains are conserved protein motifs of 38-40 amino acids found in a broad spectrum of proteins. They mediate protein-protein interactions by binding proline-rich modules in ligands. A 10 amino acid proline-rich portion of the morphogenic protein, formin, is bound in vitro by both the WW domain of the formin-binding protein 11 (FBP11) and the SH3 domain of Abl. To explore whether the FBP11 WW domain and Abl SH3 domain bind to similar ligands, we screened a mouse limb bud expression library for putative ligands of the FBP11 WW domain. In so doing, we identified eight ligands (WBP3 through WBP10), each of which contains a proline-rich region or regions. Peptide sequence comparisons of the ligands revealed a conserved motif of 10 amino acids that acts as a modular sequence binding the FBP11 WW domain, but not the WW domain of the putative signal transducing factor, hYAP65. Interestingly, the consensus ligand for the FBP11 WW domain contains residues that are also required for binding by the Abl SH3 domain. These findings support the notion that the FBP11 WW domain and the Abl SH3 domain can compete for the same proline-rich ligands and suggest that at least two subclasses of WW domains exist, namely those that bind a PPLP motif, and those that bind a PPXY motif.     

Transition states in protein folding kinetics: modeling phi-values of small beta-sheet proteins

Small single-domain proteins often exhibit only a single free-energy barrier, or transition state, between the denatured and the native state. The folding kinetics of these proteins is usually explored via mutational analysis. A central question is which structural information on the transition state can be derived from the mutational data. In this article, we model and structurally interpret mutational Φ-values for two small β-sheet proteins, the PIN and the FBP WW domains. The native structure of these WW domains comprises two β-hairpins that form a three-stranded β-sheet. In our model, we assume that the transition state consists of two conformations in which either one of the hairpins is formed. Such a transition state has been recently observed in molecular dynamics folding-unfolding simulations of a small designed three-stranded β-sheet protein. We obtain good agreement with the experimental data 1), by splitting up the mutation-induced free-energy changes into terms for the two hairpins and for the small hydrophobic core of the proteins; and 2), by fitting a single parameter, the relative degree to which hairpins 1 and 2 are formed in the transition state. The model helps us to understand how mutations affect the folding kinetics of WW domains, and captures also negative Φ-values that have been difficult to interpret.     

Principal component analysis for protein folding dynamics

Protein folding is considered here by studying the dynamics of the folding of the triple β-strand WW domain from the Formin-binding protein 28. Starting from the unfolded state and ending either in the native or nonnative conformational states, trajectories are generated with the coarse-grained united residue (UNRES) force field. The effectiveness of principal components analysis (PCA), an already established mathematical technique for finding global, correlated motions in atomic simulations of proteins, is evaluated here for coarse-grained trajectories. The problems related to PCA and their solutions are discussed. The folding and nonfolding of proteins are examined with free-energy landscapes. Detailed analyses of many folding and nonfolding trajectories at different temperatures show that PCA is very efficient for characterizing the general folding and nonfolding features of proteins. It is shown that the first principal component captures and describes in detail the dynamics of a system. Anomalous diffusion in the folding/nonfolding dynamics is examined by the mean-square displacement (MSD) and the fractional diffusion and fractional kinetic equations. The collisionless (or ballistic) behavior of a polypeptide undergoing Brownian motion along the first few principal components is accounted for.     

Mapping the transition state of the WW domain beta-sheet

The folding kinetics of a three-stranded antiparallel beta-sheet (WW domain) have been measured by temperature jump relaxation. Folding and activation free energies were determined as a function of temperature for both the wild-type and the mutant domain, W39F, which modifies the beta(2)-beta(3) hydrophobic interface. The folding rate decreases at higher temperatures as a result of the increase in the activation free energy for folding. Phi-Values were obtained for thermal perturbations allowing the primary features of the folding free energy surface to be determined. The results of this analysis indicate a significant shift from an "early" (Phi(T)=0. 4) to a "late" (Phi(T)=0.8) transition state with increasing temperature. The temperature-dependent Phi-value analysis of the wild-type WW domain and of its more stable W39F hydrophobic cluster mutant reveals little participation of residue 39 in the transition state at lower temperature. As the temperature is raised, hydrophobic interactions at the beta(2)-beta(3) interface gain importance in the transition state and the barrier height of the wild-type, which contains the larger tryptophan residue, increases more slowly than the barrier height of the mutant.     

Increasing protein stability using a rational approach combining sequence homology and structural alignment: Stabilizing the WW domain

This study shows that a combination of sequence homology and structural information can be used to increase the stability of the WW domain by 2.5 kcal mol(-1) and increase the T(m) by 28 degrees C. Previous homology-based protein design efforts typically investigate positions with low sequence identity, whereas this study focuses on semi-conserved core residues and proximal residues, exploring their role(s) in mediating stabilizing interactions on the basis of structural considerations. The A20R and L30Y mutations allow increased hydrophobic interactions because of complimentary surfaces and an electrostatic interaction with a third residue adjacent to the ligand-binding hydrophobic cluster, increasing stability significantly beyond what additivity would predict for the single mutations. The D34T mutation situated in a pi-turn possibly disengages Asn31, allowing it to make up to three hydrogen bonds with the backbone in strand 1 and loop 2. The synergistic mutations A20R/L30Y in combination with the remotely located mutation D34T add together to create a hYap WW domain that is significantly more stable than any of the protein structures on which the design was based (Pin and FBP28 WW domains).     

Phi-analysis at the experimental limits: mechanism of beta-hairpin formation

The 37-residue Formin-binding protein, FBP28, is a canonical three-stranded beta-sheet WW domain. Because of its small size, it is so insensitive to chemical denaturation that it is barely possible to determine accurately a denaturation curve, as the transition spans 0-7 M guanidinium hydrochloride (GdmCl). It is also only marginally stable, with a free energy of denaturation of just 2.3 kcal/mol at 10 degrees Celsius so only small changes in energy upon mutation can be tolerated. But these properties and relaxation times for folding of 25 micros-400 micros conspire to allow the rapid acquisition of accurate and reproducible kinetic data for Phi-analysis using classical temperature-jump methods. The transition state for folding is highly polarized with some regions having Phi-values of 0 and others 1, as readily seen in chevron plots, with Phi-values of 0 having the refolding arms overlaying and those of 1 the unfolding arms superimposable. Good agreement is seen with transition state structures identified from independent molecular dynamics (MD) simulations at 60, 75, and 100 degrees Celsius, which allows us to explore further the details of the folding and unfolding pathway of FBP28. The first beta-turn is near native-like in the transition state for folding (experimental) and unfolding (MD and experiment). The simulations show that there are transient contacts between the aromatic side-chains of the beta-strands in the denatured state and that these interactions provide the driving force for folding of the first beta-hairpin of this three-stranded sheet. Only after the backbone hydrogen bonds are formed between beta1 and beta2 does a hydrogen bond form to stabilize the intervening turn, or the first beta-turn.     

Solution structure and binding specificity of FBP11/HYPA WW domain as Group-II/III

The Group-II/III WW domains bind Pro-rich sequences, the most frequent protein motif found in eucaryotic genomes. We have proposed that the Group-II and -III WW domains be merged into a larger group because the members of each group have relatively wide specificity and bind to the common ligands [Kato et al., J Biol Chem 2004;279:31833-31841]. We have also proposed that Group-II/III has a common surface patch, the XP2 groove, to bind the ligands. The first WW domain of FBP11/HYPA is one of the Group-II/III WW domains. The solution structure of the 26 residue-long converged region exhibits an antiparallel triple stranded beta-sheet with a small hydrophobic core. The WW domain of FBP11/HYPA has both XP and XP2 grooves on its surface. Ligand titration by 1H-15N HSQC NMR spectra revealed that the WW domain of FBP11/HYPA binds all the peptides with the PL, PP, and PR motifs. The profile patterns of chemical shift perturbation were quite similar among the spectra titrated with all three ligands. In addition, the titration significantly shifts the signals of the residues that compose the XP2 groove. All these findings suggest the functional importance of the XP2 groove and group definition of Group-II/III of the WW domains.     

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.     

Solution structures of the YAP65 WW domain and the variant L30 K in complex with the peptides GTPPPPYTVG, N-(n-octyl)-GPPPY and PLPPY and the application of peptide libraries reveal a minimal binding epitope.

The single mutation L30 K in the Hu-Yap65 WW domain increased the stability of the complex with the peptide GTPPPPYTVG (K(d)=40(+/-5) microM). Here we report the refined solution structure of this complex by NMR spectroscopy and further derived structure-activity relationships by using ligand peptide libraries with truncated sequences and a substitution analysis that yielded acetyl-PPPPY as the smallest high-affinity binding peptide (K(d)=60 microM). The structures of two new complexes with weaker binding ligands chosen based on these results (N-(n-octyl)-GPPPYNH(2) and Ac-PLPPY) comprising the wild-type WW domain of Hu-Yap65 were determined. Comparison of the structures of the three complexes were useful for identifying the molecular basis of high-affinity: hydrophobic and specific interactions between the side-chains of Y28 and W39 and P5' and P4', respectively, and hydrogen bonds between T37 (donnor) and P5' (acceptor) and between W39 (donnor) and T2' (acceptor) stabilize the complex.The structure of the complex L30 K Hu-Yap65 WW domain/GTPPPPYTVG is compared to the published crystal structure of the dystrophin WW domain bound to a segment of the beta-dystroglycan protein and to the solution structure of the first Nedd4 WW domain and its prolin-rich ligand, suggesting that WW sequences bind proline-rich peptides in an evolutionary conserved fashion. The position equivalent to T22 in the Hu-Yap65 WW domain sequence is seen as responsible for differentiation in the binding mode among the WW domains of group I.     

Arginine methylation inhibits the binding of proline-rich ligands to Src homology 3, but not WW, domains

Src homology 3 (SH3) and WW domains are known to associate with proline-rich motifs within their respective ligands. Here we demonstrate that the proposed adapter protein for Src kinases, Sam68, is a ligand whose proline-rich motifs interact with the SH3 domains of p59(fyn) and phospholipase Cgamma-1 as well as with the WW domains of FBP30 and FBP21. These proline-rich motifs, in turn, are flanked by RG repeats that represent targets for the type I protein arginine N-methyltransferase. The asymmetrical dimethylation of arginine residues within these RG repeats dramatically reduces the binding of the SH3 domains of p59(fyn) and phospholipase Cgamma-1, but has no effect on their binding to the WW domain of FBP30. These results suggest that protein arginine methylation can selectively modulate certain protein-protein interactions and that mechanisms exist for the irreversible regulation of SH3 domain-mediated interactions.     

The folding pathway of an FF domain: characterization of an on-pathway intermediate state under folding conditions by (15)N, (13)C(alpha) and (13)C-methyl relaxation dispersion and (1)H/(2)H-exchange NMR spectroscopy

The FF domain from the human protein HYPA/FBP11 folds via a low-energy on-pathway intermediate (I). Elucidation of the structure of such folding intermediates and denatured states under conditions that favour folding are difficult tasks. Here, we investigated the millisecond time-scale equilibrium folding transition of the 71-residue four-helix bundle wild-type protein by (15)N, (13)C(alpha) and methyl(13)C Carr-Purcell-Meiboom-Gill (CPMG) NMR relaxation dispersion experiments and by (1)H/(2)H-exchange measurements. The relaxation data for the wild-type protein fitted a simple two-site exchange process between the folded state (F) and I. Destabilization of F in mutants A17G and Q19G allowed the detection of the unfolded state U by (15)N CPMG relaxation dispersion. The dispersion data for these mutants fitted a three-site exchange scheme, U<-->I<-->F, with I populated higher than U. The kinetics and thermodynamics of the folding reaction were obtained via temperature and urea-dependent relaxation dispersion experiments, along with structural information on I from backbone (15)N, (13)C(alpha) and side-chain methyl (13)C chemical shifts, with further information from protection factors for the backbone amide groups from (1)H/(2)H-exchange. Notably, helices H1-H3 are at least partially formed in I, while helix H4 is largely disordered. Chemical shift differences for the methyl (13)C nuclei suggest a paucity of stable, native-like hydrophobic interactions in I. These data are consistent with Phi-analysis of the rate-limiting transition state between I and F. The combination of relaxation dispersion and Phi data can elucidate whole experimental folding pathways.     

Temperature-dependent folding pathways of Pin1 WW domain: an all-atom molecular dynamics simulation of a Gō model

We study the folding thermodynamics and kinetics of the Pin1 WW domain, a three-stranded beta-sheet protein, by using all-atom (except nonpolar hydrogens) discontinuous molecular dynamics simulations at various temperatures with a Gō model. The protein exhibits a two-state folding kinetics near the folding transition temperature. A good agreement between our simulations and the experimental measurements by the Gruebele group has been found, and the simulation sheds new insights into the structure of transition state, which is hard to be straightforwardly captured in experiments. The simulation also reveals that the folding pathways at approximately the transition temperature and at low temperatures are much different, and an intermediate state at a low temperature is predicted. The transition state of this small beta-protein at its folding transition temperature has a well-established hairpin 1 made of beta1 and beta2 strands while its low-temperature kinetic intermediate has a formed hairpin 2 composed of beta2 and beta3 strands. Theoretical results are compared with other simulation results as well as available experimental data. This study confirms that specific side-chain packing in an all-atom Gō model can yield a reasonable prediction of specific folding kinetics for a given protein. Different folding behaviors at different temperatures are interpreted in terms of the interplay of entropy and enthalpy in folding process.