Detailed folding structures of M-lycotoxin-Hc1a and its mutageneses using 2D HP model

Although the hydrophobic-polar (HP) model was proposed a decade ago, it applies almost to no real-case study because of its intense computation. In this study, a 2D HP model was applied to study the folding structures of M-lycotoxin-Hc1a, an antimicrobial peptide, in order to get full pictures of its numerous folding structures. The normalised hydrophobicity index was used to convert M-lycotoxin-Hc1a and its six mutageneses into HP sequences, and then the 2D HP model was used to compute all the possible folding structures (3²⁴ = 282,429,536,481), and finally the normalised hydrophobicity index was used to distinguish the native state. The results showed that M-lycotoxin-Hc1a had 6 and 138 folding structures at their native state with the minimal energy of ? 13 at pH 2 and pH 7 when glycine served as hydrophobic amino acid. When glycine serves as polar amino acid, M-lycotoxin-Hc1a had 12 and 54 folding structures at their native state with the minimal energy of ? 12 and ? 13 at pH 2 and pH 7, respectively. This study advanced the knowledge on how to apply the HP model to real-life study, and how the mutageneses influenced the folding structures of M-lycotoxin-Hc1a, their native states and minimal energy at different pH levels.     

Analysis on folding of misgurin using two-dimensional HP model

Misgurin is an antimicrobial peptide from the loach, while the hydrophobic-polar (HP) model is a way to study the folding conformations and native states in peptide and protein although several amino acids cannot be classified either hydrophobic or polar. Practically, the HP model requires extremely intensive computations, thus it has yet to be used widely. In this study, we use the two-dimensional HP model to analyze all possible folding conformations and native states of misgurin with conversion of natural amino acids according to the normalized amino acid hydrophobicity index as well as the shortest benchmark HP sequence. The results show that the conversion of misgurin into HP sequence with glycine as hydrophobic amino acid at pH 2 has 1212 folding conformations with the same native state of minimal energy -6; the conversion of glycine as polar amino acid at pH 2 has 13,386 folding conformations with three native states of minimal energy -5; the conversion of glycine as hydrophobic amino acid at pH 7 has 2538 folding conformations with three native states of minimal energy -5; and the conversion of glycine as polar amino acid at pH 7 has 12,852 folding conformations with three native states of minimal energy -4. Those native states can be ranked according to the normalized amino acid hydrophobicity index. The detailed discussions suggest two ways to modify misgurin.     

Alternatively folded states of an immunoglobulin

Well-defined, non-native protein structures of low stability have been increasingly observed as intermediates in protein folding or as equilibrium structures populated under specific solvent conditions. These intermediate structures, frequently referred to as molten globule states, are characterized by the presence of secondary structure, a lack of significant tertiary contacts, increased hydrophobicity and partial specific volume as compared to native structures, and low cooperativity in thermal unfolding. The present study demonstrates that under acidic conditions (pH less than 3) the antibody MAK33 can assume a folded stable conformation. This A-state is characterized by a high degree of secondary structure, increased hydrophobicity, a native-like maximum wavelength of fluorescence emission, and a tendency toward slow aggregation. A prominent feature of this low-pH conformation is the stability against denaturant and thermal unfolding that is manifested in highly cooperative reversible phase transitions indicative of the existence of well-defined tertiary contacts. These thermodynamic results are corroborated by the kinetics of folding from the completely unfolded chain to the alternatively folded state at pH 2. The given data suggest that MAK33 at pH 2 adopts a cooperative structure that differs from the native immunoglobulin fold at pH 7. This alternatively folded state exhibits certain characteristics of the molten globule but differs distinctly from it by its extraordinary structural stability that is characteristic for native protein structures.     

Low folding cooperativity of HP35 revealed by single-molecule force spectroscopy and molecular dynamics simulation

Some small proteins, such as HP35, fold at submicrosecond timescale with low folding cooperativity. Although these proteins have been extensively investigated, still relatively little is known about their folding mechanism. Here, using single-molecule force spectroscopy and steered molecule dynamics simulation, we study the unfolding of HP35 under external force. Our results show that HP35 unfolds at extremely low forces without a well-defined unfolding transition state. Subsequently, we probe the structure of unfolded HP35 using the persistence length obtained in the force spectroscopy. We found that the persistence length of unfolded HP35 is around 0.72 nm, >40% longer than typical unstructured proteins, suggesting that there are a significant amount of residual secondary structures in the unfolded HP35. Molecular dynamics simulation further confirmed this finding and revealed that many native contacts are preserved in HP35, even its two ends have been extended up to 8 nm. Our results therefore suggest that retaining a significant amount of secondary structures in the unfolded state of HP35 may be an efficient way to reduce the entropic cost for the formation of tertiary structure and increase the folding speed, although the folding cooperativity is compromised. Moreover, we anticipate that the methods we used in this work can be extended to the study of other proteins with complex folding behaviors and even intrinsically disordered ones.     

Conformational states of annexin VI in solution induced by acidic pH

Acidic pH-induced folding of annexin (Anx)VI in solution was investigated in order to study the mechanism of formation of ion channels by the protein in membranes. Using 2-(p-toluidino)naphthalene-6-sulfonic acid as a hydrophobic probe, it was demonstrated that AnxVI exerts a large change in hydrophobicity at acidic pH. Moreover, circular dichroism spectra indicated that the native state of AnxVI changes at acidic pH towards a state characterized by a significant loss of alpha-helix content and appearance of new beta-structures. These changes are reversible upon an increase of pH. It is postulated that the structural folding of AnxVI could explain how a soluble protein may undergo transition into a molecule able to penetrate the membrane hydrophobic region. The physiological significance of these observations is discussed.     

Routes are trees: the parsing perspective on protein folding

An important puzzle in structural biology is the question of how proteins are able to fold so quickly into their unique native structures. There is much evidence that protein folding is hierarchic. In that case, folding routes are not linear, but have a tree structure. Trees are commonly used to represent the grammatical structure of natural language sentences, and chart parsing algorithms efficiently search the space of all possible trees for a given input string. Here we show that one such method, the CKY algorithm, can be useful both for providing novel insight into the physical protein folding process, and for computational protein structure prediction. As proof of concept, we apply this algorithm to the HP lattice model of proteins. Our algorithm identifies all direct folding route trees to the native state and allows us to construct a simple model of the folding process. Despite its simplicity, our model provides an account for the fact that folding rates depend only on the topology of the native state but not on sequence composition.     

Folding network of villin headpiece subdomain

Protein folding is a complex multidimensional process that is difficult to illustrate by the traditional analyses based on one- or two-dimensional profiles. Analyses based on transition networks have become an alternative approach that has the potential to reveal detailed features of protein folding dynamics. However, due to the lack of successful reversible folding of proteins from conventional molecular-dynamics simulations, this approach has rarely been utilized. Here, we analyzed the folding network from several 10 μs conventional molecular-dynamics reversible folding trajectories of villin headpiece subdomain (HP35). The folding network revealed more complexity than the traditional two-dimensional map and demonstrated a variety of conformations in the unfolded state, intermediate states, and the native state. Of note, deep enthalpic traps at the unfolded state were observed on the folding landscape. Furthermore, in contrast to the clear separation of the native state and the primary intermediate state shown on the two-dimensional map, the two states were mingled on the folding network, and prevalent interstate transitions were observed between these two states. A more complete picture of the folding mechanism of HP35 emerged when the traditional and network analyses were considered together.     

Crystal structure of a pH-stabilized mutant of villin headpiece

Villin-type headpiece domains are compact F-actin-binding motifs that have been used extensively as a model system to investigate protein folding by both experimental and computational methods. Villin headpiece (HP67) harbors a highly helical, thermostable, and autonomously folding subdomain in the C terminus (HP35), and because of this feature, HP67 is usually considered to be composed of a N- and C-terminal subdomain. Unlike the C-terminal subdomain, the N-terminal subdomain consists mainly of loops and turns, and the folding is dependent upon the presence of the C-terminal subdomain. The pH sensitivity of this subdomain is thought to arise from, at least partially, protonation of H41 buried in the hydrophobic core. Substitution of this histidine with tyrosine, another permissive residue at this position for naturally occurring sequences, increases not only the pH stability of HP67 but also the thermal stability and the cooperativity of thermal unfolding over a wide pH range (0.9-7.5). The crystal structures of wild-type HP67 and the H41Y mutant, determined under the same conditions, indicate that the H41Y substitution causes only localized rearrangement around the mutated residue. The F-actin-binding motif remains essentially the same after the mutation, accounting for the negligible effect of the mutation on F-actin affinity. The hydrogen bond formed between the imidazole ring of H41 and the backbone carbonyl of E14 of HP67 is eliminated by the H41Y mutation, which renders the extreme N terminus of H41Y more mobile; the hydrogen bond formed between the imidazole ring of H41 and the backbone nitrogen of D34 is replaced with that between the hydroxyl group of Y41 and the backbone nitrogen of D34 after the H41Y substitution. The increased hydrophobicity of tyrosine compensates for the loss of hydrogen bonds in the extreme N terminus and accounts for the increased stability and cooperativity of the H41Y mutant.     

Rough energy landscapes in protein folding: dimeric E. coli Trp repressor folds through three parallel channels

The folding mechanism of the dimeric Escherichia coli Trp repressor (TR) is a kinetically complex process that involves three distinguishable stages of development. Following the formation of a partially folded, monomeric ensemble of species, within 5 ms, folding to the native dimer is controlled by three kinetic phases. The rate-limiting step in each phase is either a non-proline isomerization reaction or a dimerization reaction, depending on the final denaturant concentration. Two approaches have been employed to test the previously proposed folding mechanism of TR through three parallel channels: (1) unfolding double-jump experiments demonstrate that all three folding channels lead directly to native dimer; and (2) the differential stabilization of the transition state for the final step in folding and the native dimer, by the addition of salt, shows that all three channels involve isomerization of a dimeric species. A refined model for the folding of Trp repressor is presented, in which all three channels involve a rapid dimerization reaction between partially folded monomers followed by the isomerization of the dimeric intermediates to yield native dimer. The ensemble of partially folded monomers can be captured at equilibrium by low pH; one-dimensional proton NMR spectra at pH 2.5 demonstrate that monomers exist in two distinct, slowly interconverting conformations. These data provide a potential structural explanation for the three-channel folding mechanism of TR: random association of two different monomeric forms, which are distinguished by alternative packing modes of the core dimerization domain and the DNA-binding, helix-turn-helix, domain. One, perhaps both, of these packing modes contains non-native contacts.     

Native state EX2 and EX1 hydrogen exchange of Escherichia coli CspA, a small beta-sheet protein

Escherichia coli CspA is a small all-beta-sheet protein that folds fast (tau = 4 ms) via an apparent two-state mechanism. Our previous studies have shown that a large aromatic cluster on the surface of the protein participates in the rate-limiting step of folding and thus may be part of the folding nucleus of this protein. To obtain a more detailed picture of molecular events at the peptide backbone during unfolding and folding of CspA, we used native state hydrogen exchange and nuclear magnetic resonance spectroscopy (NMR). The experiments with native CspA were performed over a range of pH values from low pH, where exchange is governed by a rapid equilibrium before chemical exchange (EX2 exchange), to high pH, where exchange is dictated by the rate of unfolding (EX1 exchange). Rates of folding and unfolding were determined for 11 residues. The distribution of rates of folding within the structure of CspA suggests that hairpin turns, including one near the aromatic cluster, may nucleate the folding of CspA.     

Structure-approximating inverse protein folding problem in the 2D HP model.

The inverse protein folding problem is that of designing an amino acid sequence which has a particular native protein fold. This problem arises in drug design where a particular structure is necessary to ensure proper protein-protein interactions. In this paper, we show that in the 2D HP model of Dill it is possible to solve this problem for a broad class of structures. These structures can be used to closely approximate any given structure. One of the most important properties of a good protein (in drug design) is its stability--the aptitude not to fold simultaneously into other structures. We show that for a number of basic structures, our sequences have a unique fold.     

Exploring the kinetic requirements for enhancement of protein folding rates in the GroEL cavity

The chaperonin system, GroEL and GroES of Escherichia coli enable certain proteins to fold under conditions when spontaneous folding is prohibitively slow as to compete with other non-productive channels such as aggregation. We investigated the plausible mechanisms of GroEL-mediated folding using simple lattice models. In particular, we have investigated protein folding in a confined environment, such as those offered by the GroEL, to decipher whether rate and yield enhancement can occur when the substrate protein is allowed to fold within the cavity of the chaperonins. The GroEL cavity is modeled as a cubic box and a simple bead model is used to represent the substrate chain. We consider three distinct characteristic of the confining environment. First, the cavity is taken to be a passive Anfinsen cage in which the walls merely reduce the available conformation space. We find that at temperatures when the native conformation is stable, the folding rate is retarded in the Anfinsen cage. We then assumed that the interior of the wall is hydrophobic. In this case the folding times exhibit a complex behavior. When the strength of the interaction between the polypeptide chain and the cavity is too strong or too weak we find that the rates of folding are retarded compared to spontaneous folding. There is an optimum range of the interaction strength that enhances the rates. Thus, above this value there is an inverse correlation between the folding rates and the strength of the substrate-cavity interactions. The optimal hydrophobic walls essentially pull the kinetically trapped states which leads to a smoother the energy landscape. It is known that upon addition of ATP and GroES the interior cavity of GroEL offers a hydrophilic-like environment to the substrate protein. In order to mimic this within the context of the dynamic Anfinsen cage model, we allow for changes in the hydrophobicity of the walls of the cavity. The duration for which the walls remain hydrophobic during one cycle of ATP hydrolysis is allowed to vary. These calculations show that frequent cycling of the wall hydrophobicity can dramatically reduce the folding times and increase the yield as well under non-permissive conditions. Examination of the structures of the substrate proteins before and after the change in hydrophobicity indicates that there is global unfolding involved. In addition, it is found that a fraction of the molecules kinetically partition to the native state in accordabce with the iterative annealing mechanism. Thus, frequent "unfoldase" activity of chaperonins leading to global unfolding of the polypeptide chain results in enhancement of the folding rates and yield of the folded protein. We suggest that chaperonin efficiency can be greatly enhanced if the cycling time is reduced. The calculations are used to interpret a few experiments on chaperonin-mediated protein folding.     

Hydrophobic Core Formation and Dehydration in Protein Folding Studied by Generalized-Ensemble Simulations

Despite its small size, chicken villin headpiece subdomain HP36 folds into the native structure with a stable hydrophobic core within several microseconds. How such a small protein keeps up its conformational stability and fast folding in solution is an important issue for understanding molecular mechanisms of protein folding. In this study, we performed multicanonical replica-exchange simulations of HP36 in explicit water, starting from a fully extended conformation. We observed at least five events of HP36 folding into nativelike conformations. The smallest backbone root mean-square deviation from the crystal structure was 1.1 Å. In the nativelike conformations, the stably formed hydrophobic core was fully dehydrated. Statistical analyses of the simulation trajectories show the following sequential events in folding of HP36: 1), Helix 3 is formed at the earliest stage; 2), the backbone and the side chains near the loop between Helices 2 and 3 take nativelike conformations; and 3), the side-chain packing at the hydrophobic core and the dehydration of the core side chains take place simultaneously at the later stage of folding. This sequence suggests that the initial folding nucleus is not necessarily the same as the hydrophobic core, consistent with a recent experimental ϕ-value analysis.     

Acid-induced partly folded conformation resembling a molten globule state of xylanase from an alkalothermophilic Bacillus sp.

Nonnative protein structures having a compact secondary, but not rigid tertiary structure, have been increasingly observed as intermediate states in protein folding. We have shown for the first time during acid-induced unfolding of xylanase (Xyl II) the presence of a partially structured intermediate form resembling a molten globule state. The conformation and stability of Xyl II at acidic pH was investigated by equilibrium unfolding methods. Using intrinsic fluorescence and CD spectroscopic studies, we have established that Xyl II at pH 1.8 (A-state) retains the helical secondary structure of the native protein at pH 7.0, while the tertiary interactions are much weaker. At variance, from the native species (N-state), Xyl II in the A-state binds 1-anilino-8-sulfonic acid (ANS) indicating a considerable exposure of aromatic side chains. Lower concentration of Gdn HCl are required to unfold the A-state. For denaturation by Gdn HCl, the midpoint of the cooperative unfolding transition measured by fluorescence for the N-state is 3.5 +/- 0.1 M, which is higher than the value (2.2 +/- 0.1 M) observed for the A-state at pH 1.8. This alternatively folded state exhibits certain characteristics of the molten globule but differs distinctly from it by its structural stability that is characteristic for native proteins.     

Extracting contact energies from protein structures: A study using a simplified model

In this study, we exploited an elementary 2-dimensional square lattice model of HP polymers to test the premise of extracting contact energies from protein structures. Given a set of prespecified energies for H–H, H–P, and P–P contacts, all possible sequences of various lengths were exhaustively enumerated to find sequences that have unique lowest-energy conformations. The lowest-energy structures (or native structures) of such (native) sequences were used to extract contact energies using the Miyazawa-Jernigan procedure and here-defined reference state. The relative magnitudes of the original energies were restored reasonably well, but the extracted contact energies were independent of the absolute magnitudes of the initial energies. We turned to a more detailed characterization of the energy landscapes of the native sequences in light of a new theoretical framework on protein folding. Foldability of such sequences imposes two limits on the absolute value of the prespecified energies: a lower bound entailed by the minimum requirement for thermodynamic stability and an upper bound associated with the entrapment of the chain to local minima. We found that these two limits confine the prespecified energy values to a rather narrow range which, surprisingly, also contains the extracted energies in all the cases examined. These results indicate that the quasi-chemical approximation can be used to connect quantitatively the occurrence of various residue–residue contacts in an ensemble of native structures with the energies of the contacts. More importantly, they suggest that the extracted contact energies do contain information on structural stability and can be used to estimate actual structural energetics. This study also encourages the use of structure-derived contact energies in threading. The finding that there is a rather narrow range of energies that are optimal for folding a sequence also cautions the use of arbitrary energy Hamiltonion in minimal folding models. Proteins 31:299–308, 1998. © 1998 Wiley-Liss, Inc.     

Folding pathway of apo-metallothionein induced by Zn2+, Cd2+ and Co2+.

Metal ion binding to the sulfhydryl groups of apometallothionein (apo-MT) causes both the formation of native metal-thiolate clusters and the folding of the polypeptide chain of each domain. Cd2+ and Zn2+ react with apo-MT to form metal-thiolate bonds in reactions that are complete within milliseconds and which are pH-dependent. Dual mixing experiments were conducted that involve the initial reaction of metal ion and apo-MT followed by mixing with 5,5'-N-dithio-bis(2-nitrobenzoate) or EDTA after 26 ms. They showed that structures had formed within the brief reaction period which were resistant to rapid reaction with reagents that interact with sulfhydryl groups or metal ions, respectively. It was concluded that native metallothionein domains had been constituted within this brief period. Apo-MT was also titrated with Co2+ to yield Co(n)-MT (n=1-7). Initially, Co2+ bound to independent, tetrahedral thiolate sites. Spectrophotometric analysis of the titration suggested that the independent Co(II) sites began to coalesce into clusters at n=4 (pH 7.2) or n=5 (pH 8.4). Back titration of free sulfhydryl groups (S) in Co(n)-MT (n=1-7) with iodoacetamide at pH 7.2 confirmed that clustering began at n=4. Upon conversion of these alkylated structures to the corresponding 113Cd2+ species 113Cd NMR spectroscopy established that the location of Co(II) in Co(n)-MT (n=1-3) was non-specific and that at n=4, the only observable structure was Co(II)4S11. The results suggest possible kinetic pathways of folding that are conceptually similar to those hypothesized for other small proteins.     

Folding mechanism of three structurally similar β-sheet proteins

The folding mechanism of cellular retinoic acid binding protein I (CRABP I), cellular retinol binding protein II (CRBP II), and intestinal fatty acid binding protein (IFABP) were investigated to determine if proteins with similar native structures have similar folding mechanisms. These mostly β-sheet proteins have very similar structures, despite having as little as 33% sequence similarity. The reversible urea denaturation of these proteins was characterized at equilibrium by circular dichroism and fluorescence. The data were best fit by a two-state model for each of these proteins, suggesting that no significant population of folding intermediates were present at equilibrium. The native states were of similar stability with free energies (linearly extrapolated to 0 M urea, ΔG) of 6.5, 8.3, and 5.5 kcal/mole for CRABP I, CRBP II, and IFABP, respectively. The kinetics of the folding and unfolding processes for these proteins was monitored by stopped-flow CD and fluorescence. Intermediates were observed during both the folding and unfolding of all of these proteins. However, the overall rates of folding and unfolding differed by nearly three orders of magnitude. Further, the spectroscopic properties of the intermediate states were different for each protein, suggesting that different amounts of secondary and/or tertiary structure were associated with each intermediate state for each protein. These data show that the folding path for proteins in the same structural family can be quite different, and provide evidence for different folding landscapes for these sequences. Proteins 33:107–118, 1998. © 1998 Wiley-Liss, Inc.     

pH-jump-induced folding and unfolding studies of barstar: evidence for multiple folding and unfolding pathways.

Equilibrium and kinetic characterization of the high pH-induced unfolding transition of the small protein barstar have been carried out in the pH range 7-12. A mutant form of barstar, containing a single tryptophan, Trp 53, completely buried in the core of the native protein, has been used. It is shown that the protein undergoes reversible unfolding above pH 10. The pH 12 form (the D form) appears to be as unfolded as the form unfolded by 6 M guanidine hydrochloride (GdnHCl) at pH 7 (the U form): both forms have similar fluorescence and far-UV circular dichroism (CD) signals and have similar sizes, as determined by dynamic light scattering and size-exclusion chromatography. No residual structure is detected in the D form: addition of GdnHCl does not alter its fluorescence and far-UV CD properties. The fluorescence signal of Trp 53 has been used to monitor folding and unfolding kinetics. The kinetics of folding of the D form in the pH range 7-11 are complex and are described by four exponential processes, as are the kinetics of unfolding of the native state (N state) in the pH range 10.5-12. Each kinetic phase of folding decreases in rate with increase in pH from 7 to 10.85, and each kinetic phase of unfolding decreases in rate with decrease in pH from 12 to 10.85. At pH 10.85, the folding and unfolding rates for any particular kinetic phase are identical and minimal. The two slowest phases of folding and unfolding have identical kinetics whether measured by Trp 53 fluorescence or by mean residue ellipticity at 222 nm. Direct determination of the increase in the N state with time of folding at pH 7 and of the D form with time of unfolding at pH 12, by means of double-jump assays, show that between 85 and 95% of protein molecules fold or unfold via fast pathways between the two forms. The remaining 5-15% of protein molecules appear to fold or unfold via slower pathways, on which at least two intermediates accumulate. The mechanism of folding from the high pH-denatured D form is remarkably similar to the mechanism of folding from the urea or GdnHCl-denatured U form.     

Salt effects on hydrophobic‐core formation in folding of a helical miniprotein studied by molecular dynamics simulations

We have investigated effects of salt ions on folding events of a helical miniprotein chicken villin headpiece subdomain HP36. Low concentrations of ions alter electrostatic interactions between charged groups of a protein and can change the populations of conformers. Here, we compare two data sets of folding simulations of HP36 in explicit water solvent with or without ions. For efficient sampling of the conformational space of HP36, the multicanonical replica‐exchange molecular dynamics method was employed. Our analyses suggest that salt alters salt‐bridging nature of the protein at later stages of folding at room temperature. Especially, more nonnative, nonlocal salt bridges are formed at near‐native conformations in pure water. Our analyses also show that such salt‐bridge formation hinders the fully native hydrophobic‐core packing at the final stages of folding. Proteins 2014; 82:933–943. © 2013 Wiley Periodicals, Inc.