Mechanism of formation of the C‐terminal β‐hairpin of the B3 domain of the immunoglobulin‐binding protein G from Streptococcus. IV. Implication for the mechanism of folding of the parent protein

A 34‐residue α/β peptide [IG(28–61)], derived from the C‐terminal part of the B3 domain of the immunoglobulin binding protein G from Streptoccocus, was studied using CD and NMR spectroscopy at various temperatures and by differential scanning calorimetry. It was found that the C‐terminal part (a 16‐residue‐long fragment) of this peptide, which corresponds to the sequence of the β‐hairpin in the native structure, forms structure similar to the β‐hairpin only at T = 313 K, and the structure is stabilized by non‐native long‐range hydrophobic interactions (Val47–Val59). On the other hand, the N‐terminal part of IG(28–61), which corresponds to the middle α‐helix in the native structure, is unstructured at low temperature (283 K) and forms an α‐helix‐like structure at 305 K, and only one helical turn is observed at 313 K. At all temperatures at which NMR experiments were performed (283, 305, and 313 K), we do not observe any long‐range connectivities which would have supported packing between the C‐terminal (β‐hairpin) and the N‐terminal (α‐helix) parts of the sequence. Such interactions are absent, in contrast to the folding pathway of the B domain of protein G, proposed recently by Kmiecik and Kolinski (Biophys J 2008, 94, 726–736), based on Monte‐Carlo dynamics studies. Alternative folding mechanisms are proposed and discussed. © 2010 Wiley Periodicals, Inc. Biopolymers 93: 469–480, 2010. This article was originally published online as an accepted preprint. The “Published Online” date corresponds to the preprint version. You can request a copy of the preprint by emailing the Biopolymers editorial office at biopolymers@wiley.com     

Commitment and nucleation in the protein G transition state

An accurate characterization of the transition state ensemble (TSE) is central to furthering our understanding of the protein folding reaction. We have extensively tested a recently reported method for studying a protein's TSE, utilizing phi-value data from protein engineering experiments and computational studies as restraints in all-atom Monte Carlo (MC) simulations. The validity of interpreting experimental phi-values as the fraction of native contacts made by a residue in the TSE was explored, revealing that this definition is unable to uniquely specify a TSE. The identification of protein G's second hairpin, in both pre and post-transition conformations demonstrates that high experimental phi-values do not guarantee a residue's importance in the TSE. An analysis of simulations based on structures restrained by experimental phi-values is necessary to yield this result, which is not obvious from a simplistic interpretation of individual phi-values. The TSE that we obtain corresponds to a single, specific nucleation event, characterized by six residues common to all three observed, convergent folding pathways. The same specific nucleus was independently identified from computational and experimental data, and "Conservation of Conservation" analysis in the protein G fold. When associated strictly with complete nucleus formation and concomitant chain collapse, folding is a well-defined two state event. Once the nucleus has formed, the folding reaction enters a slow relaxation process associated with side-chain packing and small, local backbone rearrangements. A detailed analysis of phi-values and their relationship to the transition state ensemble allows us to construct a unified theoretical model of protein G folding.     

Thermodynamics and kinetics of the hairpin ribozyme from atomistic folding/unfolding simulations

We report a set of atomistic folding/unfolding simulations for the hairpin ribozyme using a Monte Carlo algorithm. The hairpin ribozyme folds in solution and catalyzes self-cleavage or ligation via a specific two-domain structure. The minimal active ribozyme has been studied extensively, showing stabilization of the active structure by cations and dynamic motion of the active structure. Here, we introduce a simple model of tertiary-structure formation that leads to a phase diagram for the RNA as a function of temperature and tertiary-structure strength. We then employ this model to capture many folding/unfolding events and to examine the transition-state ensemble (TSE) of the RNA during folding to its active “docked” conformation. The TSE is compact but with few tertiary interactions formed, in agreement with single-molecule dynamics experiments. To compare with experimental kinetic parameters, we introduce a novel method to benchmark Monte Carlo kinetic parameters to docking/undocking rates collected over many single molecular trajectories. We find that topology alone, as encoded in a biased potential that discriminates between secondary and tertiary interactions, is sufficient to predict the thermodynamic behavior and kinetic folding pathway of the hairpin ribozyme. This method should be useful in predicting folding transition states for many natural or man-made RNA tertiary structures.     

Slow formation of aggregation‐resistant β‐sheet folding intermediates

Protein folding has been studied extensively for decades, yet our ability to predict how proteins reach their native state from a mechanistic perspective is still rudimentary at best, limiting our understanding of folding‐related processes in vivo and our ability to manipulate proteins in vitro. Here, we investigate the in vitro refolding mechanism of a large β‐helix protein, pertactin, which has an extended, elongated shape. At 55 kDa, this single domain, all‐β‐sheet protein allows detailed analysis of the formation of β‐sheet structure in larger proteins. Using a combination of fluorescence and far‐UV circular dichroism spectroscopy, we show that the pertactin β‐helix refolds remarkably slowly, with multiexponential kinetics. Surprisingly, despite the slow refolding rates, large size, and β‐sheet‐rich topology, pertactin refolding is reversible and not complicated by off‐pathway aggregation. The slow pertactin refolding rate is not limited by proline isomerization, and 30% of secondary structure formation occurs within the rate‐limiting step. Furthermore, site‐specific labeling experiments indicate that the β‐helix refolds in a multistep but concerted process involving the entire protein, rather than via initial formation of the stable core substructure observed in equilibrium titrations. Hence pertactin provides a valuable system for studying the refolding properties of larger, β‐sheet‐rich proteins, and raises intriguing questions regarding the prevention of aggregation during the prolonged population of partially folded, β‐sheet‐rich refolding intermediates. Proteins 2010. © 2009 Wiley‐Liss, Inc.     

The role of the turn in beta-hairpin formation during WW domain folding

The folding of WW domains is rate limited by formation of a beta-hairpin comprising residues from strands 1 and 2. Residues in the turn of this hairpin have reported Phi-values for folding close to 1 and have been proposed to nucleate folding. High Phi-values do not necessarily imply that the energetics of formation are a driving force for initiating folding. We demonstrate by NMR studies and molecular dynamics simulations that the first turn of the hYAP, FBP28, and PIN1 WW domains is structurally dynamic and solvent exposed in the native and folding transition states. It is, therefore, unlikely that the formation of the beta-turn per se provides the energetic driving force for hairpin folding. It is more likely that the turn acts as an easily formed hinge that facilitates the formation of the hairpin; it is a nucleus as defined by the nucleation-condensation mechanism whereby a diffuse nucleus is stabilized by associated interactions.     

β‐Hairpin folding and stability: molecular dynamics simulations of designed peptides in aqueous solution

The structural properties of a 10‐residue and a 15‐residue peptide in aqueous solution were investigated by molecular dynamics simulation. The two designed peptides, SYINSDGTWT and SESYINSDGTWTVTE, had been studied previously by NMR at 278 K and the resulting model structures were classified as 3:5 β‐hairpins with a type I + G1 β‐bulge turn. In simulations at 278 K, starting from the NMR model structure, the 3:5 β‐hairpin conformers proved to be stable over the time period evaluated (30 ns). Starting from an extended conformation, simulations of the decapeptide at 278 K, 323 K and 353 K were also performed to study folding. Over the relatively short time scales explored (30 ns at 278 K and 323 K, 56 ns at 353 K), folding to the 3:5 β‐hairpin could only be observed at 353 K. At this temperature, the collapse to β‐hairpin‐like conformations is very fast. The conformational space accessible to the peptide is entirely dominated by loop structures with different degrees of β‐hairpin character. The transitions between different types of ordered loops and β‐hairpins occur through two unstructured loop conformations stabilized by a single side‐chain interaction between Tyr2 and Trp9, which facilitates the changes of the hydrogen‐bond register. In agreement with previous experimental results, β‐hairpin formation is initially driven by the bending propensity of the turn segment. Nevertheless, the fine organization of the turn region appears to be a late event in the folding process. Copyright © 2004 European Peptide Society and John Wiley & Sons, Ltd.     

Mechanism of formation of the C‐terminal β‐hairpin of the B3 domain of the immunoglobulin binding protein G from Streptococcus. III. Dynamics of long‐range hydrophobic interactions

A 20‐residue peptide, IG(42–61), derived from the C‐terminal β‐hairpin of the B3 domain of the immunoglobulin binding protein G from Streptoccocus was studied using circular dichroism, nuclear magnetic resonance (NMR) spectroscopy at various temperatures and by differential scanning calorimetry (DSC). Unlike other related peptides studied so far, this peptide displays two heat capacity peaks in DSC measurements (at a scanning rate of 1.5 deg/min at a peptide concentration of 0.07 mM), which suggests a three‐state folding/unfolding process. The results from DSC and NMR measurements suggest the formation of a dynamic network of hydrophobic interactions stabilizing the structure, which resembles a β‐hairpin shape over a wide range of temperatures (283–313 K). Our results show that IG (42–61) possesses a well‐organized three‐dimensional structure stabilized by long‐range hydrophobic interactions (Tyr50 ··· Phe57 and Trp48 ··· Val59) at T = 283 K and (Trp48 ··· Val59) at 305 and 313 K. The mechanism of β‐hairpin folding and unfolding, as well as the influence of peptide length on its conformational properties, are also discussed. Proteins 2010. © 2009 Wiley‐Liss, Inc.     

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.     

β‐Hairpin stabilization in a 28‐residue peptide derived from the β‐subunit sequence of human chorionic gonadotropin hormone

The β‐subunit of the human chorionic gonadotropin (hCG) hormone, which is believed to be related to certain types of cancer, contains three hairpin‐like fragments. To investigate the role of β‐hairpin formation in the early stages of the hCGβ folding, a 28‐residue peptide with the sequence RDVRFESIRLPGSPRGVNPVVSYAVALS, corresponding to the H3‐β hairpin fragment (residues 60–87) of the hCGβ subunit, was studied under various conditions using three optical spectroscopic methods: Fourier transform ir spectroscopy, electronic CD, and vibrational CD. Environmental conditions are critical factors for formation of secondary structure in this peptide. TFE : H₂O mixed solvents induced helical formation. Formation of β‐structure in this peptide, which may be related to the native β‐hairpin formation in the intact hormone, was found to be induced only under conditions such as high concentration, high temperature, and the presence of nonmicellar sodium dodecyl sulfate concentrations. These findings support a protein folding mechanism for the hCGβ subunit in which an initial hydrophobic collapse, which increases intermolecular interactions in hCGβ, is needed to induce the H3‐β hairpin formation. © 1999 John Wiley & Sons, Inc. Biopoly 50: 413–423, 1999     

Dependence of folding dynamics and structural stability on the location of a hydrophobic pair in beta-hairpins

We study the dependence of folding time, nucleation site, and stability of a model beta-hairpin on the location of a cross-strand hydrophobic pair, using a coarse-grained off-lattice model with the aid of Monte Carlo simulations. Our simulations have produced 6500 independent folding trajectories dynamically, forming the basis for extensive statistical analysis. Four folding pathways, zipping-out, middle-out, zipping-in, and reptation, have been closely monitored and discussed in all seven sequences studied. A hydrophobic pair placed near the beta-turn or in the middle section effectively speed up folding; a hydrophobic pair placed close to the terminal ends or next to the beta-turn encourages stability of the entire chain.     

Mechanism of formation of the C‐terminal β‐hairpin of the B3 domain of the immunoglobulin binding protein G from Streptococcus. II. Interplay of local backbone conformational dynamics and long‐range hydrophobic interactions in hairpin formation

Two peptides, corresponding to the turn region of the C‐terminal β‐hairpin of the B3 domain of the immunoglobulin binding protein G from Streptococcus, consisting of residues 51–56 [IG(51–56)] and 50–57 [IG(50–57)], respectively, were studied by circular dichroism and NMR spectroscopy at various temperatures and by differential scanning calorimetry. Our results show that the part of the sequence corresponding to the β‐turn in the native structure (DDATKT) of the B3 domain forms bent conformations similar to those observed in the native protein. The formation of a turn is observed for both peptides in a broad range of temperatures (T = 283–323 K), which confirms the conclusion drawn from our previous studies of longer sequences from the C‐terminal β‐hairpin of the B3 domain of the immunoglobulin binding protein G (16, 14, and 12 residues), that the DDATKT sequence forms a nucleation site for formation of the β‐hairpin structure of peptides corresponding to the C‐terminal part of all the B domains of the immunoglobulin binding protein G. We also show and discuss the role of long‐range hydrophobic interactions as well as local conformational properties of polypeptide chains in the mechanism of formation of the β‐hairpin structure. Proteins 2009. © 2009 Wiley‐Liss, Inc.     

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.     

Thermodynamics and folding kinetics analysis of the SH3 domain form discrete molecular dynamics

We perform a detailed analysis of the thermodynamics and folding kinetics of the SH3 domain fold with discrete molecular dynamic simulations. We propose a protein model that reproduces some of the experimentally observed thermodynamic and folding kinetic properties of proteins. Specifically, we use our model to study the transition state ensemble of the SH3 fold family of proteins, a set of unstable conformations that fold to the protein native state with probability 1/2. We analyze the participation of each secondary structure element formed at the transition state ensemble. We also identify the folding nucleus of the SH3 fold and test extensively its importance for folding kinetics. We predict that a set of amino acid contacts between the RT-loop and the distal hairpin are the critical folding nucleus of the SH3 fold and propose a hypothesis that explains this result.     

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.     

Mechanism of formation of the C‐terminal β‐hairpin of the B3 domain of the immunoglobulin binding protein G from Streptococcus. I. Importance of hydrophobic interactions in stabilization of β‐hairpin structure

We previously studied a 16‐amino acid‐residue fragment of the C‐terminal β‐hairpin of the B3 domain (residues 46–61), [IG(46–61)] of the immunoglobulin binding protein G from Streptoccocus, and found that hydrophobic interactions and the turn region play an important role in stabilizing the structure. Based on these results, we carried out systematic structural studies of peptides derived from the sequence of IG (46–61) by systematically shortening the peptide by one residue at a time from both the C‐ and the N‐terminus. To determine the structure and stability of two resulting 12‐ and 14‐amino acid‐residue peptides, IG(48–59) and IG(47–60), respectively, we carried out circular dichroism, NMR, and calorimetric studies of these peptides in pure water. Our results show that IG(48–59) possesses organized three‐dimensional structure stabilized by hydrophobic interactions (Tyr50–Phe57 and Trp48–Val59) at T = 283 and 305 K. At T = 313 K, the structure breaks down because of increased chain entropy, but the turn region is preserved in the same position observed for the structure of the whole protein. The breakdown of structure occurs near the melting temperature of this peptide (T_m = 310 K) measured by differential scanning calorimetry (DSC). The melting temperature of IG(47–60) determined by DSC is T_m = 330 K and its structure is similar to that of the native β‐hairpin at all (lower) temperatures examined (283–313 K). Both of these truncated sequences are conserved in all known amino acid sequences of the B domains of the immunoglobulin binding protein G from bacteria. Thus, this study contributes to an understanding of the mechanism of folding of this whole family of proteins, and provides information about the mechanism of formation and stabilization of a β‐hairpin structural element. Proteins 2009. © 2008 Wiley‐Liss, Inc.     

Molecular dynamics simulations of certain mutant peptide models from staphylococcal nuclease reveal that initial hydrophobic collapse associated with turn propensity drives β‐hairpin folding

An important nucleation event during the folding of staphylococcal nuclease involves the formation of a β‐hairpin by the sequence ²¹DTVKLMYKGQPMTFR³⁵. Earlier studies show that the turn sequence ‘YKGQP’ has an important role in the folding of this β‐hairpin. To understand the active or passive nature of the turn sequence ‘YKGQP’ in the folding of the aforementioned β‐hairpin sequence, we studied glycine mutant peptides Ac‐²DTVKLMYGGQPMTFR¹⁶‐NMe (K9G:15), Ac‐²DTVKLMYKGGPMTFR¹⁶‐NMe (Q11G:15), Ac‐²DTVKLMYGGGPMTFR¹⁶‐NMe (K9G/Q11G:15), and Ac‐²DTVKLMGGGGGMTFR¹⁶‐NMe (penta‐G:15) by using molecular dynamics simulations, starting with two different unfolded states, polyproline II and extended conformational forms. Further, 5mer mutant turn peptides Ac‐²YGGQP⁶‐NMe (K3G:5), Ac‐²YKGGP⁶‐NMe (Q5G:5), Ac‐²YGGGP⁶‐NMe (K3G/Q5G:5), and Ac‐²GGGGG⁶‐NMe (penta‐G:5) were also studied individually. Our results show that an initial hydrophobic collapse and loop closure occurs in all 15mer mutants, but only K9G:15 mutant forms a stable native‐like β‐hairpin. In the other 15mer mutants, the hydrophobic collapsed state would not proceed to β‐hairpin formation. Of the different simulations performed for the penta‐G:15 mutant, in only one simulation a nonnative β‐hairpin conformation is sampled with highly flexible loop region (⁸GGGGG¹²), which has no specific conformational preference as a 5mer. While the sequence ‘YGGQP’ in the K3G:5 simulation shows relatively higher β‐turn propensity, the presence of this sequence in K9G:15 peptide seems to be driving the β‐hairpin formation. Thus, these results seem to suggest that for the formation of a stable β‐hairpin, the initial hydrophobic collapse is to be assisted by a turn propensity. Initial hydrophobic collapse alone is not sufficient to guide β‐hairpin formation. Copyright © 2013 European Peptide Society and John Wiley & Sons, Ltd.     

Folding thermodynamics of β‐hairpins studied by replica‐exchange molecular dynamics simulations

We study the differences in folding stability of β‐hairpin peptides, including GB1 hairpin and a point mutant GB1 K10G, as well as tryptophan zippers (TrpZips): TrpZip1, TrpZip2, TrpZip3‐1, and TrpZip4. By performing replica‐exchange molecular dynamics simulations with Amber03* force field (a modified version of Amber ff03) in explicit solvent, we observe ab initio folding of all the peptides except TrpZip3‐1, which is experimentally known to be the least stable among the peptides studied here. By calculating the free energies of unfolding of the peptides at room temperature and folding midpoint temperatures for thermal unfolding of peptides, we find that TrpZip4 and GB1 K10G peptides are the most stable β‐hairpins followed by TrpZip1, GB1, and TrpZip2 in the given order. Hence, the proposed K10G mutation of GB1 peptide results in enhanced stability compared to wild‐type GB1. An important goal of our study is to test whether simulations with Amber 03* model can reproduce experimentally predicted folding stability differences between these peptides. While the stabilities of GB1 and TrpZip1 yield close agreement with experiment, TrpZip2 is found to be less stable than predicted by experiment. However, as heterogenous folding of TrpZip2 may yield divergent thermodynamic parameters by different spectroscopic methods, mismatching of results with previous experimental values are not conclusive of model shortcomings. For most of the cases, molecular simulations with Amber03* can successfully reproduce experimentally known differences between the mutated peptides, further highlighting the predictive capabilities of current state‐of‐the‐art all‐atom protein force fields. Proteins 2015; 83:1307–1315. © 2015 Wiley Periodicals, Inc.     

Crystal structure of BinB: A receptor binding component of the binary toxin from Lysinibacillus sphaericus

The binary toxin (Bin), produced by Lysinibacillus sphaericus, is composed of BinA (42 kDa) and BinB (51 kDa) proteins, which are both required for full toxicity against Culex and Anopheles mosquito larvae. Specificity of Bin toxin is determined by the binding of BinB component to a receptor present on the midgut epithelial membranes, while BinA is proposed to be a toxic component. Here, we determined the first crystal structure of the active form of BinB at a resolution of 1.75 Å. BinB possesses two distinct structural domains in its N‐ and C‐termini. The globular N‐terminal domain has a β‐trefoil scaffold which is a highly conserved architecture of some sugar binding proteins or lectins, suggesting a role of this domain in receptor‐binding. The BinB β‐rich C‐terminal domain shares similar three‐dimensional folding with aerolysin type β‐pore forming toxins, despite a low sequence identity. The BinB structure, therefore, is a new member of the aerolysin‐like toxin family, with probably similarities in the cytolytic mechanism that takes place via pore formation. Proteins 2014; 82:2703–2712. © 2014 Wiley Periodicals, Inc.     

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.