The effect of context on the folding of β-hairpins

Small β-hairpin peptides have been widely used as models for the folding of β-sheets. But how applicable is the folding of such models to β-structure in larger proteins with conventional hydrophobic cores? Here we present multiple unfolding simulations of three such proteins that contain the WW domain double hairpin β-sheet motif: cold shock protein A (CspA), cold shock protein B (CspB) and glucose permease IIA domain. We compare the behavior of the free motif in solution and in the context of proteins of different size and architecture. Both Csp proteins lost contacts between the double-hairpin motif and the protein core as the first step of unfolding and proceeded to unfold with loss of the third β-strand, similar to the isolated WW domain. The glucose permease IIA domain is a larger protein and the contacts between the motif and the core were not lost as quickly. Instead the unfolding pathway of glucose permease IIA followed a different pathway with β1 pulling away from the sheet first. Interestingly, when the double hairpin motif was excised from the glucose permease IIA domain and simulated in isolation in water it unfolded by the same pathway as the WW domain, indicating that it is tertiary interactions with the protein that alter the motif’s unfolding not a sequence dependent effect on its intrinsic unfolding behavior. With respect to the unfolding of the hairpins, there was no consistent order to the loss of hydrogen bonds between the β-strands in the hairpins in any of the systems. Our results show that while the folding behavior of the isolated WW domain is generally consistent with the double hairpin motif’s behavior in the cold shock proteins, it is not the case for the glucose permease IIA domain. So, one must be cautious in extrapolating findings from model systems to larger more complicated proteins where tertiary interactions can overwhelm intrinsic behavior.     

Redesign of a WW domain peptide for selective recognition of single-stranded DNA

A β-sheet miniprotein based on the FBP11 WW1 domain sequence has been redesigned for the molecular recognition of ssDNA. A previous report showed that a β-hairpin peptide dimer, (WKWK)(2), binds ssDNA with low micromolar affinity but with little selectivity over duplex DNA. This report extends those studies to a three-stranded β-sheet miniprotein designed to mimic the OB-fold. The new peptide binds ssDNA with low micromolar affinity and shows about 10-fold selectivity for ssDNA over duplex DNA. The redesigned peptide no longer binds its native ligand, the polyproline helix, confirming that the peptide has been redesigned for the function of binding ssDNA. Structural studies provide evidence that this peptide consists of a well-structured β-hairpin made of strands 2 and 3 with a less structured first strand that provides affinity for ssDNA but does not improve the stability of the full peptide. These studies provide insight into protein-DNA interactions as well as a novel example of protein redesign.     

Effects of mutation, truncation, and temperature on the folding kinetics of a WW domain

The purpose of this work is to show how mutation, truncation, and change of temperature can influence the folding kinetics of a protein. This is accomplished by principal component analysis of molecular-dynamics-generated folding trajectories of the triple β-strand WW domain from formin binding protein 28 (FBP28) (Protein Data Bank ID: 1E0L) and its full-size, and singly- and doubly-truncated mutants at temperatures below and very close to the melting point. The reasons for biphasic folding kinetics [i.e., coexistence of slow (three-state) and fast (two-state) phases], including the involvement of a solvent-exposed hydrophobic cluster and another delocalized hydrophobic core in the folding kinetics, are discussed. New folding pathways are identified in free-energy landscapes determined in terms of principal components for full-size mutants. Three-state folding is found to be a main mechanism for folding the FBP28 WW domain and most of the full-size and truncated mutants. The results from the theoretical analysis are compared to those from experiment. Agreements and discrepancies between the theoretical and experimental results are discussed. Because of its importance in understanding protein kinetics and function, the diffusive mechanism by which the FBP28 WW domain and its full-size and truncated mutants explore their conformational space is examined in terms of the mean-square displacement and principal component analysis eigenvalue spectrum analyses. Subdiffusive behavior is observed for all studied systems.     

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.     

(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.     

Synthesis and β-sheet propensity of constrained N-amino peptides

The stabilization of β-sheet secondary structure through peptide backbone modification represents an attractive approach to protein mimicry. Here, we present strategies toward stable β-hairpin folds based on peptide strand N-amination. Novel pyrazolidinone and tetrahydropyridazinone dipeptide constraints were introduced via on-resin Mitsunobu cyclization between α-hydrazino acid residues and a serine or homoserine side chain. Acyclic and cyclic N-amino peptide building blocks were then evaluated for their effect on β-hairpin stability in water using a GB1-derived model system. Our results demonstrate the strong β-sheet stabilizing effect of the peptide N-amino substituent, and provide useful insights into the impact of covalent dipeptide constraint on β-sheet folding.     

Stability of the beta-sheet of the WW domain: A molecular dynamics simulation study.

The WW domain consists of approximately 40 residues, has no disulfide bridges, and forms a three-stranded antiparallel beta-sheet that is monomeric in solution. It thus provides a model system for studying beta-sheet stability in native proteins. We performed molecular dynamics simulations of two WW domains, YAP65 and FBP28, with very different stability characteristics, in order to explore the initial unfolding of the beta-sheet. The less stable YAP domain is much more sensitive to simulation conditions than the FBP domain. Under standard simulation conditions in water (with or without charge-balancing counterions) at 300 K, the beta-sheet of the YAP WW domain disintegrated at early stages of the simulations. Disintegration commenced with the breakage of a hydrogen bond between the second and third strands of the beta-sheet due to an anticorrelated transition of the Tyr-28 psi and Phe-29 phi angles. Electrostatic interactions play a role in this event, and the YAP WW domain structure is more stable when simulated with a complete explicit model of the surrounding ionic strength. Other factors affecting stability of the beta-sheet are side-chain packing, the conformational entropy of the flexible chain termini, and the binding of cognate peptide.     

An Overlapping Region between the Two Terminal Folding Units of the Outer Surface Protein A (OspA) Controls Its Folding Behavior

Although many naturally occurring proteins consist of multiple domains, most studies on protein folding to date deal with single-domain proteins or isolated domains of multi-domain proteins. Studies of multi-domain protein folding are required for further advancing our understanding of protein folding mechanisms. Borrelia outer surface protein A (OspA) is a β-rich two-domain protein, in which two globular domains are connected by a rigid and stable single-layer β-sheet. Thus, OspA is particularly suited as a model system for studying the interplays of domains in protein folding. Here, we studied the equilibria and kinetics of the urea-induced folding–unfolding reactions of OspA probed with tryptophan fluorescence and ultraviolet circular dichroism. Global analysis of the experimental data revealed compelling lines of evidence for accumulation of an on-pathway intermediate during kinetic refolding and for the identity between the kinetic intermediate and a previously described equilibrium unfolding intermediate. The results suggest that the intermediate has the fully native structure in the N-terminal domain and the single layer β-sheet, with the C-terminal domain still unfolded. The observation of the productive on-pathway folding intermediate clearly indicates substantial interactions between the two domains mediated by the single-layer β-sheet. We propose that a rigid and stable intervening region between two domains creates an overlap between two folding units and can energetically couple their folding reactions.     

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).     

Contributions of aromatic pairs to the folding and stability of long-lived human γD-crystallin

Human γD-crystallin (HγD-Crys) is a highly stable protein that remains folded in the eye lens for the majority of an individual's lifetime. HγD-Crys exhibits two homologous crystallin domains, each containing two Greek key motifs with eight β-strands. Six aromatic pairs (four Tyr/Tyr, one Tyr/Phe and one Phe/Phe) are present in the β-hairpin sequences of the Greek keys. Ultraviolet damage to the aromatic residues in lens crystallins may contribute to the genesis of cataract. Mutant proteins with these aromatic residues substituted with alanines were constructed and expressed in E. coli. All mutant proteins except F115A and F117A had lower thermal stability than the WT protein. In equilibrium experiments in guanidine hydrochloride (GuHCl), all mutant proteins had lower thermodynamic stability than the WT protein. N-terminal domain (N-td) substitutions shifted the N-td transition to lower GuHCl concentration, but the C-terminal domain (C-td) transition remained unaffected. C-td substitutions led to a more cooperative unfolding/refolding process, with both the N-td and C-td transitions shifted to lower GuHCl concentration. The aromatic pairs conserved for each Greek key motif (Greek key pairs) had larger contributions to both thermal stability and thermodynamic stability than the other pairs. Aromatic-aromatic interaction was estimated as 1.5-2.0 kcal/mol. In kinetic experiments, N-td substitutions accelerated the early phase of unfolding, while C-td substitutions accelerated the late phase, suggesting independent domain unfolding. Only substitutions of the second Greek key pair of each crystallin domain slowed refolding. The second Greek keys may provide nucleation sites during the folding of the double-Greek-key crystallin domains.     

Structure of FBP11 WW1-PL ligand complex reveals the mechanism of proline-rich ligand recognition by group II/III WW domains

FBP11/HYPA is a mammalian homologue of yeast splicing factor Prp40. The first WW domain of FBP11/HYPA (FBP11 WW1) is essential for preventing severe neurological diseases such as Huntington disease and Rett syndrome and strongly resembles the WW domain of FCA, the essential regulator for flowering time control. We have solved the structure of FBP11 WW1 and a Pro-Pro-Leu-Pro ligand complex, and demonstrated the binding mechanism with mutational analysis using surface plasmon resonance. The overall structure of FBP11 WW1 in the complex form is quite similar to the structures of WW domains from Group I and IV in complexes. In addition, conformation of FBP11 WW1 does not change much upon ligand binding. The binding orientation of the ligand against FBP11 WW1 is the same as that of the Group IV WW domain-ligand complex, but opposite to that of the Group I complex. The ligand interacts with two grooves formed by surface aromatic residues. The Pro and Leu residues in the ligand interact with the grooves and the Loop I region of FBP11 WW1, respectively, which are necessary interactions for binding the ligand. Interestingly, the two aromatic grooves recognize the Pro residues in entirely different manners, which allows FBP11 WW1 to recognize shorter sequences than the SH3 domain. Combined with homology models of other WW domains, the present report shows the detailed mechanism of ligand binding by Group II/III WW domains, and provides information useful in designing drugs to treat neurodegenerative diseases.     

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.     

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.     

Computational design and experimental testing of the fastest-folding β-sheet protein

One of the most important and elusive goals of molecular biology is the formulation of a detailed, atomic-level understanding of the process of protein folding. Fast-folding proteins with low free-energy barriers have proved to be particularly productive objects of investigation in this context, but the design of fast-folding proteins was previously driven largely by experiment. Dramatic advances in the attainable length of molecular dynamics simulations have allowed us to characterize in atomic-level detail the folding mechanism of the fast-folding all-β WW domain FiP35. In the work reported here, we applied the biophysical insights gained from these studies to computationally design an even faster-folding variant of FiP35 containing only naturally occurring amino acids. The increased stability and high folding rate predicted by our simulations were subsequently validated by temperature-jump experiments. The experimentally measured folding time was 4.3 μs at 80 °C—about three times faster than the fastest previously known protein with β-sheet content and in good agreement with our prediction. These results provide a compelling demonstration of the potential utility of very long molecular dynamics simulations in redesigning proteins well beyond their evolved stability and folding speed.     

Dual energy landscape: the functional state of the β-barrel outer membrane protein G molds its unfolding energy landscape

We applied dynamic single-molecule force spectroscopy to quantify the parameters (free energy of activation and distance of the transition state from the folded state) characterizing the energy barriers in the unfolding energy landscape of the outer membrane protein G (OmpG) from Escherichia coli. The pH-dependent functional switching of OmpG directs the protein along different regions on the unfolding energy landscape. The two functional states of OmpG take the same unfolding pathway during the sequential unfolding of β-hairpins I-IV. After the initial unfolding events, the unfolding pathways diverge. In the open state, the unfolding of β-hairpin V in one step precedes the unfolding of β-hairpin VI. In the closed state, β-hairpin V and β-strand S11 with a part of extracellular loop L6 unfold cooperatively, and subsequently β-strand S12 unfolds with the remaining loop L6. These two unfolding pathways in the open and closed states join again in the last unfolding step of β-hairpin VII. Also, the conformational change from the open to the closed state witnesses a rigidified extracellular gating loop L6. Thus, a change in the conformational state of OmpG not only bifurcates its unfolding pathways but also tunes its mechanical properties for optimum function.     

A Study of the Influence of Charged Residues on β-Hairpin Formation by Nuclear Magnetic Resonance and Molecular Dynamics

Chain reversals are often nucleation sites in protein folding. The β-hairpins of FBP28 WW domain and IgG are stable and have been proved to initiate the folding and are, therefore, suitable for studying the influence of charged residues on β-hairpin conformation. In this paper, we carried out NMR examination of the conformations in solution of two fragments from the FPB28 protein (PDB code: 1E0L) (N-terminal part) namely KTADGKT-NH₂ (1E0L 12–18, D7) and YKTADGKTY-NH₂ (1E0L 11–19, D9), one from the B3 domain of the protein G (PDB code: 1IGD), namely DDATKT-NH₂ (1IGD 51–56) (Dag1), and three variants of Dag1 peptide: DVATKT-NH₂ (Dag2), OVATKT-NH₂ (Dag3) and KVATKT-NH₂ (Dag4), respectively, in which the original charged residue were replaced with non-polar residues or modified charged residues. It was found that both the D7 and D9 peptides form a large fraction bent conformations. However, no hydrophobic contacts between the terminal Tyr residues of D9 occur, which suggests that the presence of a pair of like-charged residues stabilizes chain reversal. Conversely, only the Dag1 and Dag2 peptides exhibit some chain reversal; replacing the second aspartic-acid residue with a valine and the first one with a basic residue results in a nearly extended conformation. These results suggest that basic residues farther away in sequence can result in stabilization of chain reversal owing to screening of the non-polar core. Conversely, smaller distance in sequence prohibits this screening, while the presence oppositely-charged residues can stabilize a turn because of salt-bridge formation.     

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.     

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.     

Special Interaction of Anionic Phosphatidic Acid Promotes High Secondary Structure in Tetrameric Potassium Channel

Anionic phosphatidic acid (PA) has been shown to stabilize and bind stronger than phosphatidylglycerol via electrostatic and hydrogen bond interaction with the positively charged residues of potassium channel KcsA. However, the effects of these lipids on KcsA folding or secondary structure are not clear. In this study, the secondary structure analyses of KcsA potassium channel was carried out using circular dichroism spectroscopy. It was found that PA interaction leads to increases in α-helical and β-sheet content of KcsA protein. In PA, KcsA α-helical structure was further stabilized by classical membrane-active cosolvent trifluoroethanol followed by reduction in the β-sheet content indicating cooperative transformation from the β-sheet to an α-helical structure. The data further uncover the role of anionic PA in KcsA folding and provide mechanism by which strong hydrogen bonds/electrostatic interaction among PA headgroup and basic residues on lipid binding domains may induce high helical structure thereby altering the protein folding and increasing the stability of tetrameric assembly.     

Structural and thermodynamic studies of Bergerac-SH3 chimeras

Bergerac-type chimeras of spectrin SH3 were designed by extending a β-hairpin by eight amino acids so that the extension protruded from the domain body like a “nose” being exposed to the solvent. A calorimetric study of several Bergerac-SH3 variants was carried out over a wide range of pH values and protein concentrations and the three-dimensional structure of one of them, SHH, was determined. X-ray studies confirmed that the nose had a well defined β-structure whilst the chimera formed a stable tetramer within the crystal unit because of four tightly packed noses. In the pH range of 4–7 the heat-induced unfolding of some chimeras was complex and concentration dependent, whilst at pH values below 3.5, low protein concentrations of all the chimeras studied, including SHH, seemed to obey a monomolecular two-state unfolding model. The best set of data was obtained for the SHA variant, the unfolding heat effects of which were systematically higher than those of the WT protein (about 16.4 kJ/mol at 323 K), which may be close to the upper limit of the enthalpy gain due to 10 residue β-hairpin folding. At the same time, the chimeras with high nose stability, which, like SHH, have a hydrophobic (IVY) cluster on their surface, showed a lower apparent unfolding heat effect, much closer to that of the WT protein. The possible reasons for this difference are discussed.