Peptide models of protein folding initiation sites. 3. The G-H helical hairpin of myoglobin

From Immunogenic Peptides to Intrinsically Disordered Proteins

It is hard to evaluate the role of individual mentors in the genesis of important ideas. In the case of our realization that proteins do not have to be stably folded to be functional, the influence of Richard Lerner and our collaborative work in the 1980s on the conformations of immunogenic peptides provided a base level of thinking about the nature of polypeptides in water solutions that led us to formulate and develop our ideas on the importance of intrinsic disorder in proteins. This review describes how the insights gained into the behavior of peptides led directly to the realization that proteins were not only capable of being functional while disordered, but also that disorder provided a distinct functional advantage in many important cellular processes.

Thermal unfolding of myoglobin in the Landau-Ginzburg-Wilson approach

The Landau-Ginzburg-Wilson paradigm is applied to model the low-temperature crystallographic Cα backbone structure of sperm whale myoglobin. The Glauber protocol is employed to simulate its response to an increase in ambient temperature. The myoglobin is found to unfold from its native state by a succession of α-helical intermediates, fully in line with the observed folding and unfolding patterns in denaturation experiments. In particular, a molten globule intermediate is identified with experimentally correct attributes. A detailed, experimentally testable contact map is constructed to characterize the specifics of the unfolding pathway, including the formation of long-range interactions. The results reveal how the unfolding process of a protein is driven by the interplay between, and a successive melting of, its modular secondary structure components.

Evidence for a Shared Mechanism in the Formation of Urea-Induced Kinetic and Equilibrium Intermediates of Horse Apomyoglobin from Ultrarapid Mixing Experiments

In this study, the equivalence of the kinetic mechanisms of the formation of urea-induced kinetic folding intermediates and non-native equilibrium states was investigated in apomyoglobin. Despite having similar structural properties, equilibrium and kinetic intermediates accumulate under different conditions and via different mechanisms, and it remains unknown whether their formation involves shared or distinct kinetic mechanisms. To investigate the potential mechanisms of formation, the refolding and unfolding kinetics of horse apomyoglobin were measured by continuous- and stopped-flow fluorescence over a time range from approximately 100 μs to 10 s, along with equilibrium unfolding transitions, as a function of urea concentration at pH 6.0 and 8°C. The formation of a kinetic intermediate was observed over a wider range of urea concentrations (0–2.2 M) than the formation of the native state (0–1.6 M). Additionally, the kinetic intermediate remained populated as the predominant equilibrium state under conditions where the native and unfolded states were unstable (at ~0.7–2 M urea). A continuous shift from the kinetic to the equilibrium intermediate was observed as urea concentrations increased from 0 M to ~2 M, which indicates that these states share a common kinetic folding mechanism. This finding supports the conclusion that these intermediates are equivalent. Our results in turn suggest that the regions of the protein that resist denaturant perturbations form during the earlier stages of folding, which further supports the structural equivalence of transient and equilibrium intermediates. An additional folding intermediate accumulated within ~140 μs of refolding and an unfolding intermediate accumulated in <1 ms of unfolding. Finally, by using quantitative modeling, we showed that a five-state sequential scheme appropriately describes the folding mechanism of horse apomyoglobin.

Hierarchical folding mechanism of apomyoglobin revealed by ultra-fast H/D exchange coupled with 2D NMR

The earliest steps in the folding of proteins are complete on an extremely rapid time scale that is difficult to access experimentally. We have used rapid-mixing quench-flow methods to extend the time resolution of folding studies on apomyoglobin and elucidate the structural and dynamic features of members of the ensemble of intermediate states that are populated on a submillisecond time scale during this process. The picture that emerges is of a continuum of rapidly interconverting states. Even after only 0.4 ms of refolding time a compact state is formed that contains major parts of the A, G, and H helices, which are sufficiently well folded to protect amides from exchange. The B, C, and E helix regions fold more slowly and fluctuate rapidly between open and closed states as they search docking sites on this core; the secondary structure in these regions becomes stabilized as the refolding time is increased from 0.4 to 6 ms. No further stabilization occurs in the A, G, H core at 6 ms of folding time. These studies begin to time-resolve a progression of compact states between the fully unfolded and native folded states and confirm the presence an ensemble of intermediates that interconvert in a hierarchical sequence as the protein searches conformational space on its folding trajectory.

Interactions responsible for secondary structure formation during folding of equine beta-lactoglobulin

Equine beta-lactoglobulin forms a compact intermediate at an acidic pH (A state). It also forms an expanded and helical conformation at low temperatures (C state). The structure of a single disulfide mutant C66A/C160A is similar to the A state in the presence of salts, while it is similar to the C state at low anion concentrations. We have investigated the temperature-dependent change in the secondary structure using circular dichroism and proline scanning mutagenesis. At low anion concentrations, the helical content increased linearly as temperature decreased. In the presence of salts, the A state was cooperatively transformed into the C state at low temperatures. This suggests the importance of hydrophobic interactions for stabilizing the A state. Peptides encompassing native-like and non-native alpha-helices were synthesized to investigate the interactions responsible for helix formation in the A and C states. These did not form stable helices, indicating that not only the helices in the A state but also the helices in the C state are stabilized by long-range interactions. A longer fragment, CHIBL, which encompasses the structured region in the A and C states, showed a helical structure. Proline-substituted mutants of CHIBL showed CD spectral changes similar to the corresponding mutants of the full-length protein in the C state. Therefore, CHIBL has a structure similar to the corresponding region of the full-length protein in the C state. This result indicates that interactions responsible for helix formation in the C state reside in the sequence of CHIBL, and that the sequences outside CHIBL are essential for secondary structure formation in the A state.

The role of hydrophobic interactions in initiation and propagation of protein folding

Globular proteins fold by minimizing the nonpolar surface that is exposed to water, while simultaneously providing hydrogen-bonding interactions for buried backbone groups, usually in the form of secondary structures such as alpha-helices, beta-sheets, and tight turns. A primary thermodynamic driving force for the formation of globular structure is thus the sequestration of nonpolar groups, but the correlation between the parts of proteins that are observed to fold first (termed folding initiation sites) and the "hydrophobicity" (as customarily defined) of the amino acids in these regions has been quite weak. It has previously been noted that many amino acid side chains contain considerable nonpolar sections, even if they also contain polar or charged groups. For example, a lysine side chain contains four methylenes, which may undergo hydrophobic interactions if the charged epsilon-NH(3)(+) group is salt-bridged or hydrogen-bonded. Folding initiation sites might therefore contain not only accepted "hydrophobic" amino acids, but also larger charged side chains. Recent experiments on the folding of mutant apomyoglobins provides corroboration for models based on the hypothesis that folding initiation sites arise from hydrophobic interactions. A near-perfect correlation was observed between the areas of the molecule that are present in the burst-phase kinetic intermediate and both the free energy of formation of hydrophobic initiation sites and the parameter "average area buried upon folding," which pinpoints large side chains, even those containing charged or polar portions. These results provide a putative mechanism for the control of protein-folding initiation and growth by polar/nonpolar sequence propensity alone.

Folding very short peptides using molecular dynamics

Peptides often have conformational preferences. We simulated 133 peptide 8-mer fragments from six different proteins, sampled by replica-exchange molecular dynamics using Amber7 with a GB/SA (generalized-Born/solvent-accessible electrostatic approximation to water) implicit solvent. We found that 85 of the peptides have no preferred structure, while 48 of them converge to a preferred structure. In 85% of the converged cases (41 peptides), the structures found by the simulations bear some resemblance to their native structures, based on a coarse-grained backbone description. In particular, all seven of the β hairpins in the native structures contain a fragment in the turn that is highly structured. In the eight cases where the bioinformatics-based I-sites library picks out native-like structures, the present simulations are largely in agreement. Such physics-based modeling may be useful for identifying early nuclei in folding kinetics and for assisting in protein-structure prediction methods that utilize the assembly of peptide fragments.

To carry out specific biochemical reactions, proteins must adopt precise three-dimensional conformations. During the folding of a protein, the protein picks out the right conformation out of billions of other conformations. It is not yet possible to do this computationally. Picking out the native conformation using physics-based atomically detailed models, sampled by molecular dynamics, is presently beyond the reach of computer methods. How can we speed up computational protein-structure prediction? One idea is that proteins start folding at specific parts of a chain that kink up early in the folding process. If we can identify these kinks, we should be able to speed up protein-structure prediction. Previous studies have identified likely kinks through bioinformatic analysis of existing protein structures. The goal of the authors here is to identify these putative folding initiation sites with a physical model instead. In this study, Ho and Dill show that, by chopping a protein chain into peptide pieces, then simulating the pieces in molecular dynamics, they can identify those peptide fragments that have conformational biases. These peptides identify the kinks in the protein chain.

Steady-state and transient ultraviolet resonance Raman spectrometer for the 193-270 nm spectral region

We describe a state-of-the-art tunable ultraviolet (UV) Raman spectrometer for the 193-270 nm spectral region. This instrument allows for steady-state and transient UV Raman measurements. We utilize a 5 kHz Ti-sapphire continuously tunable laser (approximately 20 ns pulse width) between 193 nm and 240 nm for steady-state measurements. For transient Raman measurements we utilize one Coherent Infinity YAG laser to generate nanosecond infrared (IR) pump laser pulses to generate a temperature jump (T-jump) and a second Coherent Infinity YAG laser that is frequency tripled and Raman shifted into the deep UV (204 nm) for transient UV Raman excitation. Numerous other UV excitation frequencies can be utilized for selective excitation of chromophoric groups for transient Raman measurements. We constructed a subtractive dispersion double monochromator to minimize stray light. We utilize a new charge-coupled device (CCD) camera that responds efficiently to UV light, as opposed to the previous CCD and photodiode detectors, which required intensifiers for detecting UV light. For the T-jump measurements we use a second camera to simultaneously acquire the Raman spectra of the water stretching bands (2500-4000 cm(-1)) whose band-shape and frequency report the sample temperature.

Searching for folding initiation sites of staphylococcal nuclease: A study of N-terminal short fragments

The N-terminal short fragments of staphylococcal nuclease (SNase), SNase20, SNase28, and SNase36, corresponding to the sequence regions, Ala1-Gly20, Ala1-Lys28, and Ala1-Leu36, respectively, as well as an 8-residue peptide (Ala17-Ile18-Asp19-Gly20-Asp21-Thr22-Val23-Lys24) have been synthesized. The conformational states of these fragments were investigated using CD and NMR spectroscopy in aqueous solution and in trifluoroethanol (TFE)-H(2)O mixture. SNase20 containing a sequence corresponding to a bent peptide in native SNase shows a transient population of bend-like conformation around Ala12-Thr13-Leu14 in TFE-H(2)O mixture. The sequence region of Ala17-Thr22 of SNase28 displays a localized propensity for turn-like conformation in both aqueous solution and TFE-H(2)O mixture. The conformational ensemble of SNase36 in aqueous solution includes populated turn-like conformations localized in sequence regions Ala17-Thr22 and Tyr27-Gln30. The analysis suggests that these sequence regions, which form the regular secondary structures in native protein, may serve as the folding nucleation sites of SNase fragments of different chain lengths starting from the N-terminal end. Thus, the formation of bend- and turn-like conformations of these sequence regions may be involved in the early folding events of the SNase polypeptide chain in vitro.

Similarity of force-induced unfolding of apomyoglobin to its chemical-induced unfolding: an atomistic molecular dynamics simulation approach

We have compared force-induced unfolding with traditional unfolding methods using apomyoglobin as a model protein. Using molecular dynamics simulation, we have investigated the structural stability as a function of the degree of mechanical perturbation. Both anisotropic perturbation by stretching two terminal atoms and isotropic perturbation by increasing the radius of gyration of the protein show the same key event of force-induced unfolding. Our primary results show that the native structure of apomyoglobin becomes destabilized against the mechanical perturbation as soon as the interhelical packing between the G and H helices is broken, suggesting that our simulation results share a common feature with the experimental observation that the interhelical contact is more important for the folding of apomyoglobin than the stability of individual helices. This finding is further confirmed by simulating both helix destabilizing and interhelical packing destabilizing mutants.

Chain length dependence of apomyoglobin folding: structural evolution from misfolded sheets to native helices

Very little is known about how protein structure evolves during the polypeptide chain elongation that accompanies cotranslational protein folding. This in vitro model study is aimed at probing how conformational space evolves for purified N-terminal polypeptides of increasing length. These peptides are derived from the sequence of an all-alpha-helical single domain protein, Sperm whale apomyoglobin (apoMb). Even at short chain lengths, ordered structure is found. The nature of this structure is strongly chain length dependent. At relatively short lengths, a predominantly non-native beta-sheet conformation is present, and self-associated amyloid-like species are generated. As chain length increases, alpha-helix progressively takes over, and it replaces the beta-strand. The observed trends correlate with the specific fraction of solvent-accessible nonpolar surface area present at different chain lengths. The C-terminal portion of the chain plays an important role by promoting a large and cooperative overall increase in helical content and by consolidating the monomeric association state of the full-length protein. Thus, a native-like energy landscape develops late during apoMb chain elongation. This effect may provide an important driving force for chain expulsion from the ribosome and promote nearly-posttranslational folding of single domain proteins in the cell. Nature has been able to overcome the above intrinsic misfolding trends by modulating the composition of the intracellular environment. An imbalance or improper functioning by the above modulating factors during translation may play a role in misfolding-driven intracellular disorders.

Primary folding dynamics of sperm whale apomyoglobin: core formation

The structure, thermodynamics, and kinetics of heat-induced unfolding of sperm whale apomyoglobin core formation have been studied. The most rudimentary core is formed at pH(*) 3.0 and up to 60 mM NaCl. Steady state for ultraviolet circular dichroism and fluorescence melting studies indicate that the core in this acid-destabilized state consists of a heterogeneous composition of structures of approximately 26 residues, two-thirds of the number involved for horse heart apomyoglobin under these conditions. Fluorescence temperature-jump relaxation studies show that there is only one process involved in Trp burial. This occurs in 20 micro s for a 7 degrees jump to 52 degrees C, which is close to the limits placed by diffusion on folding reactions. However, infrared temperature jump studies monitoring native helix burial are biexponential with times of 5 micro s and 56 micro s for a similar temperature jump. Both fluorescence and infrared fast phases are energetically favorable but the slow infrared absorbance phase is highly temperature-dependent, indicating a substantial enthalpic barrier for this process. The kinetics are best understood by a multiple-pathway kinetics model. The rapid phases likely represent direct burial of one or both of the Trp residues and parts of the G- and H-helices. We attribute the slow phase to burial and subsequent rearrangement of a misformed core or to a collapse having a high energy barrier wherein both Trps are solvent-exposed.

NMR evidence for different conformations of the bioactive region of rat CCK-8 and CCK-58

Sulfated CCK-58 and CCK-8 have identical bioactive C-terminal primary sequences but distinct C-terminal solution structures and different bioactivities. To examine structural differences in greater detail, rat CCK-58 and -8 were synthesized with isotopic enrichment of C-terminal residues with (15)N at alpha-amino nitrogens. Proton and nitrogen chemical shift assignments of peptide solutions were obtained by homo- and heteronuclear NMR methods. These data show that the tertiary structure ensembles of C-terminal CCK-8 and CCK-58 differ significantly. Thus, distinct solution conformations may explain differences in CCK(A) and CCK(B) receptor interactions of large and small molecular forms of CCK.

Structural determinants outside the PXDLS sequence affect the interaction of adenovirus E1A, C-terminal interacting protein and Drosophila repressors with C-terminal binding protein

C-Terminal binding protein (CtBP) interacts with a highly conserved amino acid motif (PXDLS) at the C terminus of adenovirus early region 1A (AdE1A) protein. This amino acid sequence has recently been demonstrated in the mammalian protein C-terminal interacting protein (CtIP) and a number of Drosophila repressors including Snail, Knirps and Hairy. In the study described here we have examined the structures of synthetic peptides identical to the CtBP binding sites on these proteins using NMR spectroscopy. It has been shown that peptides identical to the CtBP binding site in CtIP and at the N terminus of Snail form a series of beta-turns similar to those seen in AdE1A. The PXDLS motif towards the C terminus of Snail forms an alpha-helix. However, the motifs in Knirps and Hairy did not adopt well-defined structures in TFE/water mixtures as shown by the absence of medium range NOEs and a high proportion of signal overlap. The affinities of peptides for Drosophila and mammalian CtBP were compared using enzyme-linked immunosorbent assay. CtIP, Snail (N-terminal peptide) and Knirps peptides all bind to mammalian CtBP with high affinity (K(i) of 1.04, 1.34 and 0.52 microM, respectively). However, different effects were observed with dCtBP, most notably the affinity for the Snail (N-terminal peptide) and Knirps peptides were markedly reduced (K(i) of 332 and 56 microM, respectively) whilst the Hairy peptide bound much more strongly (K(i) for dCtBP of 6.22 compared to 133 microM for hCtBP). In addition we have shown that peptides containing identical PXDLS motifs but with different N and C terminal sequences have appreciably different affinities for mammalian CtBP and different structures in solution. We conclude that the factors governing the interactions of CtBPs with partner proteins are more complex than simple possession of the PXDLS motif. In particular the overall secondary structures and amino acid side chains in the binding sites of partner proteins are of importance as well as possible global structural effects in both members of the complex. These data are considered evidence for a multiplicity of CtBPs and partner proteins in the cell.

Anatomy of protein structures: visualizing how a one-dimensional protein chain folds into a three-dimensional shape

Here, we depict the anatomy of protein structures in terms of the protein folding process. Via an iterative, top-down dissecting procedure, tertiary structures are spliced down to reveal their anatomy: first, to produce domains (defined by visual three-dimensional inspection criteria); then, hydrophobic folding units (HFU); and, at the end of a multilevel process, a set of building blocks. The resulting anatomy tree organization not only clearly depicts the organization of a one-dimensional polypeptide chain in three-dimensional space but also straightforwardly describes the most likely folding pathway(s). Comparison of the tree with the formation of the hydrophobic folding units through combinatorial assembly of the building blocks illustrates how the chain folds in a sequential or a complex folding pathway. Further, the tree points to the kinetics of the folding, whether the chain is a fast or a slow folder, and the probability of misfolding. Our ability to successfully dissect the protein into an anatomy tree illustrates that protein folding is a hierarchical process and further validates a building blocks protein folding model.

High-pressure denaturation of apomyoglobin

The pressure denaturation of wild type and mutant apomyoglobin (apoMb) was investigated using a high-pressure, high-resolution nuclear magnetic resonance and high-pressure fluorescence techniques. Wild type apoMb is resistant to pressures up to 80 MPa, and denatures to a high-pressure intermediate, I(p), between 80 and 200 MPa. A further increase of pressure to 500 MPa results in denaturation of the intermediate. The two tryptophans, both in the A helix, remain sequestered from solvent in the high-pressure intermediate, which retains some native NOESY cross peaks in the AGH core as well as between F33 and F43. High-pressure fluorescence shows that the tryptophans remain inaccessible to solvent in the I(p) state. Thus the high-pressure intermediate has some structural properties in common with the apoMb I(2) acid intermediate. The resistance of the AGH core to pressures up to 200 MPa provides further evidence that the intrinsic stability of these alpha-helices is responsible for their presence in a number of equilibrium intermediates as well as in the earliest kinetic folding intermediate. Mutations in the AGH core designed to disrupt packing by burying a charge or increasing the size of a hydrophobic residue significantly perturbed the unfolding of native apoMb to the high-pressure intermediate. The F123W and S108L mutants both unfolded at lower pressures, while retaining some resistance to pressures below 50 MPa. The charge burial mutants, A130K and S108K, are not stable at very low pressures and both denature to the intermediate by 100 MPa, half of the pressure required for wild type apoMb. Thus a similar intermediate state is created independent of the method of perturbation, and mutations have similar effects on native state destabilization for both methods of denaturation. These data suggest that equilibrium intermediates that can be formed through different means are likely to resemble a kinetic intermediate.

Peptide analogs from E-cadherin with different calcium-binding affinities

Cadherins are a family of calcium-dependent cell-surface proteins that are fundamental in controlling the development and maintenance of tissues. Motif B of E-cadherin seems to be a crucial calcium-binding site as single point mutations (D134A and D134K) completely inactivate its adhesion activity. We analyzed peptide models corresponding to motif B (amino acids 128-144) as well as selected mutations of this motif. Our NMR studies showed that this motif B sequence is actually an active calcium-binding region, even in the absence of the rest of the cadherin molecule. We found that the binding affinity of this motif is very sensitive to mutations. For example, our peptide P128-144 with the native calcium-binding sequence has an affinity of Kd 0.4 mM, whereas the mutants P128-144/ D134A and P128-144/D134K containing the replacement of Asp134 by Ala and Lys, have Kd values of only 1.5 and 11 mM, respectively. Removing Asp at position 134, which correlates with the loss of adhesion activity, decreases calcium-binding affinity 20-fold. Ala132, along with residues Asp134, Asp136 and Asn143, is involved in calcium binding in solution. We also demonstrated that the calcium-binding affinity can be increased 3-fold when an additional Asp is introduced at position 132. In 50% organic solvent, this binding affinity of peptide P128-144/A132D (17-mer) from E-cadherin is similar to that of peptide P72-100/C73-77-91A (29-mer) from alpha-lactalbumin.

Distinguishing between sequential and nonsequentially folded proteins: implications for folding and misfolding

We describe here an algorithm for distinguishing sequential from nonsequentially folding proteins. Several experiments have recently suggested that most of the proteins that are synthesized in the eukaryotic cell may fold sequentially. This proposed folding mechanism in vivo is particularly advantageous to the organism. In the absence of chaperones, the probability that a sequentially folding protein will misfold is reduced significantly. The problem we address here is devising a procedure that would differentiate between the two types of folding patterns. Footprints of sequential folding may be found in structures where consecutive fragments of the chain interact with each other. In such cases, the folding complexity may be viewed as being lower. On the other hand, higher folding complexity suggests that at least a portion of the polypeptide backbone folds back upon itself to form three-dimensional (3D) interactions with noncontiguous portion(s) of the chain. Hence, we look at the mechanism of folding of the molecule via analysis of its complexity, that is, through the 3D interactions formed by contiguous segments on the polypeptide chain. To computationally splice the structure into consecutively interacting fragments, we either cut it into compact hydrophobic folding units or into a set of hypothetical, transient, highly populated, contiguous fragments ("building blocks" of the structure). In sequential folding, successive building blocks interact with each other from the amino to the carboxy terminus of the polypeptide chain. Consequently, the results of the parsing differentiate between sequentially vs. nonsequentially folded chains. The automated assessment of the folding complexity provides insight into both the likelihood of misfolding and the kinetic folding rate of the given protein. In terms of the funnel free energy landscape theory, a protein that truly follows the mechanism of sequential folding, in principle, encounters smoother free energy barriers. A simple sequentially folded protein should, therefore, be less error prone and fold faster than a protein with a complex folding pattern.

Helix capping in the GCN4 leucine zipper

Capping interactions associated with specific sequences at or near the ends of alpha-helices are important determinants of the stability of protein secondary and tertiary structure. We investigate here the role of the helix-capping motif Ser-X-X-Glu, a sequence that occurs frequently at the N termini of alpha helices in proteins, on the conformation and stability of the GCN4 leucine zipper. The 1.8 A resolution crystal structure of the capped molecule reveals distinct conformations, packing geometries and hydrogen-bonding networks at the amino terminus of the two helices in the leucine zipper dimer. The free energy of helix stabilization associated with the hydrogen-bonding and hydrophobic interactions in this capping structure is -1.2 kcal/mol, evaluated from thermal unfolding experiments. A single cap thus contributes appreciably to stabilizing the terminated helix and thereby the native state. These results suggest that helix capping plays a further role in protein folding, providing a sensitive connector linking alpha-helix formation to the developing tertiary structure of a protein.

Influence of disulfide bonds on the induction of helical conformation in proteins

The effect(s) of TFE (2,2,2-trifluoroethanol) on three different conformational states (native, denatured, and carboxymethylated) of CTX III and RNase A has been examined. Contrary to the general belief, the results of the present study reveal that TFE can induce helical conformation in a protein which has no sequence propensity to form a helix. It is found that the helix induction in TFE is intricately related to the destabilization of the tertiary structural conformation in proteins. More importantly, the disulfide bonds in proteins are found to have significant influence on the TFE-mediated helix induction. The results obtained in this study strongly suggest that information pertaining to the influence of disulfide bonds on helix induction need to be considered to improve the accuracy of secondary structure prediction algorithms.

Fast events in protein folding: the time evolution of primary processes

Most experimental studies on the dynamics of protein folding have been confined to timescales of 1 ms and longer. Yet it is obvious that many phenomena that are obligatory elements of the folding process occur on much faster timescales. For example, it is also now clear that the formation of secondary and tertiary structures can occur on nanosecond and microsecond times, respectively. Although fast events are essential to, and sometimes dominate, the overall folding process, with a few exceptions their experimental study has become possible only recently with the development of appropriate techniques. This review discusses new approaches that are capable of initiating and monitoring the fast events in protein folding with temporal resolution down to picoseconds. The first important results from those techniques, which have been obtained for the folding of some globular proteins and polypeptide models, are also discussed.

Effect of H helix destabilizing mutations on the kinetic and equilibrium folding of apomyoglobin

Acid-denatured apomyoglobin (apoMb) contains residual helical structure in the region of the polypeptide which corresponds to the H helix of the folded protein. In order to elucidate the role of this residual secondary structure in the protein folding process and to determine whether residual structure in the denatured state affects either the overall rate of folding or the rate of formation of a burst phase intermediate, we have examined the equilibrium and kinetic folding behavior of a mutant designed to destabilize residual secondary structure in the H helix region. Both Asn132 and Glu136 were changed to Gly (N132G,E136G) to effect this destabilization. Circular dichroism spectra show that the mutant protein contains less helical structure in the acid-denatured state and in the equilibrium intermediate state at pH 4.2 than does the wild-type protein. The CD spectra of the native states of the two proteins are nearly identical. The refolding kinetics for each of the species were measured by stopped-flow CD in the far-UV region and by NMR quench-flow pulse labeling. Under identical conditions, the CD-detected refolding of wild-type and mutant apomyoglobin from the acid-denatured state or from the urea-denatured state occurs at very similar rates following a burst phase that occurs too rapidly to measure by the stopped-flow technique. The urea dependence of the unfolding and refolding rates is consistent with the presence of at least one obligatory on-pathway intermediate in both wild-type and mutant proteins. The kinetic intermediate of the mutant protein is considerably less stable than that of the wild-type protein. Hydrogen exchange pulse labeling experiments indicate that, in contrast to the wild-type protein, the H helix is not stabilized during the burst phase refolding of the mutant but becomes stabilized during the slower phases. While the wild-type and mutant proteins both form compact intermediates, these differ in the content and location of secondary structure. The rate of folding of the AGH subdomain, which takes place prior to the transition state, is substantially slower for the N132G,E136G mutant protein. A strong propensity for spontaneous formation of helical structure in the H helix region is not a prerequisite for efficient folding nor for formation of equilibrium or kinetic intermediates. These observations suggest that while folding of apomyoglobin proceeds through an obligatory intermediate, the precise structure of this intermediate is not critical and its secondary structure may be altered without substantially affecting either the overall refolding kinetics or the integrity of the final folded state.

Elucidating the folding problem of alpha-helices: local motifs, long-range electrostatics, ionic-strength dependence and prediction of NMR parameters

The information about the conformational behavior of monomeric helical peptides in solution, as well as the alpha-helix stability in proteins, has been previously utilized to derive a database with the energy contributions for various interactions taking place in an alpha-helix: intrinsic helical propensities, side-chain-side-chain interactions, main-chain-main-chain hydrogen bonds, and capping effects. This database was implemented in an algorithm based on the helix/coil transition theory (AGADIR). Here, we have modified this algorithm to include previously described local motifs: hydrophobic staple, Schellman motif and Pro-capping motif, new variants of these, and newly described side-chain-side-chain interactions. Based on recent experimental data we have introduced a position dependence of the helical propensities for some of the 20 amino acid residues. A new electrostatic model that takes into consideration all electrostatic interactions up to 12 residues in distance in the helix and random-coil conformations, as well as the effect of ionic strength, has been implemented. We have synthesized and analyzed several peptides, and used data from peptides already analysed by other groups, to test the validity of our electrostatic model. The modified algorithm predicts, with an overall standard deviation value of 6.6 (maximum helix is 100%), the helical, content of 778 peptides of which 223 correspond to wild-type and modified protein fragments. To improve the prediction potential of the algorithm and to have a direct comparison with nuclear magnetic resonance data, the algorithm now predicts the conformational shift of the CalphaH protons, 13Calpha and 3JalphaN values. We have found that for those peptides correctly predicted from the point of view of circular dichroism, the prediction of the NMR parameters is very good.

Peptide models of local and long-range interactions in the molten globule state of human alpha-lactalbumin

alpha-Lactalbumin, a small calcium-binding protein, forms an equilibrium molten globule state under a variety of conditions. A set of four peptides designed to probe the role of local interactions and the role of potential long-range interactions in stabilizing the molten globule of alpha-lactalbumin has been prepared. The first peptide consists of residues 20 through 36 of human alpha-lactalbumin and includes the entire B-helix. This peptide is unstructured in solution as judged by CD. The second peptide is derived from residues 101 through 120 and contains both the D and 310 helices. When this peptide is crosslinked via the native 28 to 111 disulfide to the B-helix peptide, a dramatic increase in helicity is observed. The crosslinked peptide is monomeric, as judged by analytical ultracentrifugation. The peptide binds 1-anilinonaphthalene-8-sulphonate (ANS) and the fluorescence emission maximum of the construct is consistent with partial solvent exposure of the tryptophan residues. The peptide corresponding to residues 101 to 120 adopts significant non-random structure in aqueous solution at low pH. Two hydrophobic clusters, one involving residues 101 through 104 and the other residues 115 through 119 have been identified and characterized by NMR. The hydrophobic cluster formed by residues 101 through 104 is still present in a smaller peptide containing only residues 101 to 111 of alpha-lactalbumin. The cluster also persists in 6 M urea. A non-native, pH-dependent interaction between the Y103 and H107 side-chains that was previously identified in the acid-denatured molten globule state was examined. This interaction was found to be more prevalent at low pH and may therefore be an example of a local interaction that stabilizes preferentially the acid-induced molten globule state.

Structural determinants present in the C-terminal binding protein binding site of adenovirus early region 1A proteins

The C-terminal binding protein (CtBP) has previously been shown to bind to a highly conserved six-amino acid motif very close to the C terminus of adenovirus early region 1A (Ad E1A) proteins. We have developed an enzyme-linked immunosorbent assay that has facilitated the screening of synthetic peptides identical or similar to the binding site on Ad E1A for their ability to bind CtBP and thus inhibit its interaction with Ad12 E1A. It has been shown that amino acids both C-terminal and N-terminal to the original proposed binding site contribute to the interaction of peptides with CtBP. Single amino acid substitutions across the binding site appreciably alter the Kd of the peptide for CtBP, indicative of a marked reduction in the affinity of the peptide for CtBP. The solution structures of synthetic peptides equivalent to the C termini of both Ad5 and Ad12 E1A and two substituted forms of these have been determined by proton NMR spectroscopy. Both the Ad12 and Ad5 peptides dissolved in trifluoroethanol/water mixtures were found to adopt regular secondary structural conformations seen as a series of beta-turns. An Ad12 peptide bearing a substitution that resulted in only very weak binding to CtBP (Ad12 L258G) was found to be random coil in solution. However, a second mutant (Ad12 V256K), which bound to CtBP rather more strongly (although not as well as the wild type), adopted a conformation similar to that of the wild type. We conclude that secondary structure (beta-turns) and an appropriate series of amino acid side chains are necessary for recognition by CtBP.

UV resonance Raman determination of protein acid denaturation: selective unfolding of helical segments of horse myoglobin

We have used UV resonance Raman spectroscopy to study the acid denaturation of horse heart aquometmyoglobin (Mb) between pH 7.5 and 1.5. Raman spectra excited at 206.5 nm are dominated by amide vibrations, which are analyzed by using a new methodology to quantitatively determine the Mb secondary structure. In contrast, the 229-nm Raman spectra are dominated by the Tyr and Trp Raman bands, which are analyzed to examine changes in Tyr and Trp environments, such as exposure to water, hydrogen bonding, and, for Trp, any alterations of the dihedral angle between the Trp ring and its linkage to the protein backbone. We uniquely determined which Mb alpha-helices melt by combining the amide, Tyr, and Trp Raman spectral information with heme absorption spectral information. We calculate that the Mb alpha-helical composition decreases from approximately 80% at neutral pH to approximately 19% below pH 3.5. The Trp Raman cross sections dramatically decrease at low pH to values which indicate that they are fully exposed to water; this result indicates that the A helix melts. The Tyr Raman bands are pH independent, which indicates that the G and H helices around the Tyr residues do not melt. The dramatic heme absorption acid denaturation changes indicate major alterations of the heme pocket and changes in heme binding. These results indicate that the A, B, C, D, E, and F helices melt in a concerted fashion, while the antiparallel G and H helices only partially melt.

The early folding kinetics of apomyoglobin

The folding pathway of apomyoglobin has been experimentally shown to have early kinetic intermediates involving the A, B, G, and H helices. The earliest detected kinetic events occur on a ns to micros time scale. We show that the early folding kinetics of apomyoglobin may be understood as the association of nascent helices through a network of diffusion-collision-coalescence steps G + H <--> GH + A <--> AGH + B <--> ABGH obtained by solving the diffusion-collision model in a chemical kinetics approximation. Our reproduction of the experimental results indicates that the model is a useful way to analyze folding data. One prediction from our fit is that the nascent A and H helices should be relatively more helix-like before coalescence than the other apomyoglobin helices.

Structural heterogeneity of the various forms of apomyoglobin: implications for protein folding

Temperature-induced denaturation transitions of different structural forms of apomyoglobin were studied monitoring intrinsic tryptophan fluorescence. It was found that the tryptophans are effectively screened from solvent both in native and acid forms throughout most of the temperature range tested. Thus, the tryptophans' surrounding do not show a considerable change in structure where major protein conformational transitions have been found in apomyoglobin using other techniques. At high temperatures and under strong destabilizing conditions, the tryptophans' fluorescence parameters show sigmoidal thermal denaturation. These results, combined with previous studies, show that the structure of this protein is heterogeneous, including native-like (tightly packed) and molten globule-like substructures that exhibit conformation (denaturation) transitions under different conditions of pH and temperature (and denaturants). The results suggest that the folding of this protein proceeds via two "nucleation" events whereby native-like contacts are formed. One of these events, which involves AGH "core" formation, appears to occur very early in the folding process, even before significant hydrophobic collapse in the rest of the protein molecule. From the current studies and other results, a rather detailed picture of the folding of myoglobin is presented, on the level of specific structures and their thermodynamical properties as well as formation kinetics.

The A-kinase anchoring domain of type IIalpha cAMP-dependent protein kinase is highly helical

Subcellular localization of the type II cAMP-dependent protein kinase is controlled by interaction of the regulatory subunit with A-Kinase Anchoring Proteins (AKAPs). This contribution examines the solution structure of a 44-residue region that is sufficient for high affinity binding to AKAPs. The N-terminal dimerization domain of the type IIalpha regulatory subunit of cAMP-dependent protein kinase was expressed to high levels on minimal media and uniformly isotopically enriched with 15N and 13C nuclei. Sequence-specific backbone and side chain resonance assignments have been made for greater than 95% of the amino acids in the free dimerization domain using high resolution multidimensional heteronuclear NMR techniques. Contrary to the results from secondary structure prediction algorithms, our analysis indicates that the domain is highly helical with a single 3-5-residue sequence involved in a beta-strand. The assignments and secondary structure analysis provide the basis for analyzing the structure and dynamics of the dimerization domain both free and complexed with specific anchoring proteins.

Thermodynamic model of secondary structure for alpha-helical peptides and proteins

A thermodynamic model describing formation of alpha-helices by peptides and proteins in the absence of specific tertiary interactions has been developed. The model combines free energy terms defining alpha-helix stability in aqueous solution and terms describing immersion of every helix or fragment of coil into a micelle or a nonpolar droplet created by the rest of protein to calculate averaged or lowest energy partitioning of the peptide chain into helical and coil fragments. The alpha-helix energy in water was calculated with parameters derived from peptide substitution and protein engineering data and using estimates of nonpolar contact areas between side chains. The energy of nonspecific hydrophobic interactions was estimated considering each alpha-helix or fragment of coil as freely floating in the spherical micelle or droplet, and using water/cyclohexane (for micelles) or adjustable (for proteins) side-chain transfer energies. The model was verified for 96 and 36 peptides studied by 1H-nmr spectroscopy in aqueous solution and in the presence of micelles, respectively ([set 1] and [set 2]) and for 30 mostly alpha-helical globular proteins ([set 3]). For peptides, the experimental helix locations were identified from the published medium-range nuclear Overhauser effects detected by 1H-nmr spectroscopy. For sets 1, 2, and 3, respectively, 93, 100, and 97% of helices were identified with average errors in calculation of helix boundaries of 1.3, 2.0, and 4.1 residues per helix and an average percentage of correctly calculated helix-coil states of 93, 89, and 81%, respectively. Analysis of adjustable parameters of the model (the entropy and enthalpy of the helix-coil transition, the transfer energy of the helix backbone, and parameters of the bound coil), determined by minimization of the average helix boundary deviation for each set of peptides or proteins, demonstrates that, unlike micelles, the interior of the effective protein droplet has solubility characteristics different from that for cyclohexane, does not bind fragments of coil, and lacks interfacial area.

Folding propensities of peptide fragments of myoglobin

Myoglobin has been studied extensively as a paradigm for protein folding. As part of an ongoing study of potential folding initiation sites in myoglobin, we have synthetized a series of peptides covering the entire sequence of sperm whale myoglobin. We report here on the conformation preferences of a series of peptides that cover the region from the A helix to the FG turn. Structural propensities were determined using circular dichroism and nuclear magnetic resonance spectroscopy in aqueous solution, trifluoroethanol, and methanol. Peptides corresponding to helical regions in the native protein, namely the B, C, D, and E helices, populate the alpha region of (phi, psi) space in water solution but show no measurable helix formation except in the presence of trifluoroethanol. The F-helix sequence has a much lower propensity to populate helical conformations even in TFE. Despite several attempts, we were not successful in synthesizing a peptide corresponding to the A-helix region that was soluble in water. A peptide termed the AB domain was constructed spanning the A- and B-helix sequences. The AB domain is not soluble in water, but shows extensive helix formation throughout the peptide when dissolved in methanol, with a break in the helix at a site close to the A-B helix junction in the intact folded myoglobin protein. With the exception of one local preference for a turn conformation stabilized by hydrophobic interactions, the peptides corresponding to turns in the folded protein do not measurably populate beta-turn conformations in water, and the addition of trifluoroethanol does not enhance the formation of either helical or turn structure. In contrast to the series of peptides described here, either studies of peptides from the GH region of myoglobin show a marked tendency to populate helical structures (H), nascent helical structures (G), or turn conformations (GH peptide) in water solution. This region, together with the A-helix and part of the B-helix, has been shown to participate in an early folding intermediate. The complete analysis of conformational properties of isolated myoglobin peptides supports the hypothesis that spontaneous secondary structure formation in local regions of the polypeptide may play an important role in the initiation of protein folding.

Acceleration of the folding of hen lysozyme by trifluoroethanol

The refolding of a partially structured state of hen lysozyme formed in 60% (v/v) 2,2,2-trifluoroethanol (TFE) has been studied using hydrogen exchange pulse labelling monitored by 2D 1H NMR, and by stopped flow fluorescence and CD measurements. The results are compared with similar studies of the refolding of the protein denatured in 6 M guanidine hydrochloride (GuHCl). Two conclusions have emerged from these studies. First, provided that the refolding conditions are identical, the two denatured states fold with very similar kinetics, despite the fact the extensive secondary structure is present in the TFE-denatured state but not in the protein denatured in 6 M GuHCl. This arises because of the rapid equilibration of structure in the species formed in the initial stage of folding. Second, whilst addition of GuHCl to the refolding buffer decreases the rate of folding, low concentrations of TFE increase the rate of folding. The result is consistent with slow steps in the refolding of lysozyme being associated primarily with the reorganisation of hydrophobic interactions rather than of hydrogen bonded structure.

Insights into protein folding from NMR

NMR has emerged as an important tool for studies of protein folding because of the unique structural insights it can provide into many aspects of the folding process. Applications include measurements of kinetic folding events and structural characterization of folding intermediates, partly folded states, and unfolded states. Kinetic information on a time scale of milliseconds or longer can be obtained by real-time NMR experiments and by quench-flow hydrogen-exchange pulse labeling. Although NMR cannot provide direct information on the very rapid processes occurring during the earliest stages of protein folding, studies of isolated peptide fragments provide insights into likely protein folding initiation events. Multidimensional NMR techniques are providing new information on the structure and dynamics of protein folding intermediates and both partly folded and unfolded states.

Packing interactions in the apomyglobin folding intermediate

The contribution of specific packing to the stability of the sperm whale apomyoglobin intermediate has been studied by urea denaturation monitored by circular dichroism and fluorescence. Mutations disrupting native packing sites within the subdomain formed by the A, G and H helices destabilize the intermediate, in contrast to the conclusion drawn from earlier studies of pH-induced unfolding. Based on these results, the intermediate is proposed to be stabilized by both partially formed native-like tertiary, and non-specific hydrophobic interactions forming a subdomain folding intermediate. The results help to explain how the intermediate acquires its structure and stability.

Molecular dynamics simulations of N-terminal peptides from a nucleotide binding protein

Molecular dynamics (MD) simulations of N-terminal peptides from lactate dehydrogenase (LDH) with increasing length and individual secondary structure elements were used to study their stability in relation to folding. Ten simulations of 1-2 ns of different peptides in water starting from the coordinates of the crystal structure were performed. The stability of the peptides was compared qualitatively by analyzing the root mean square deviation (RMSD) from the crystal structure, radius of gyration, secondary and tertiary structure, and solvent accessible surface area. In agreement with earlier MD studies, relatively short (< 15 amino acids) peptides containing individual secondary structure elements were generally found to be unstable; the hydrophobic alpha 1-helix of the nucleotide binding fold displayed a significantly higher stability, however. Our simulations further showed that the first beta alpha beta supersecondary unit of the characteristic dinucleotide binding fold (Rossmann fold) of LDH is somewhat more stable than other units of similar length and that the alpha 2-helix, which unfolds by itself, is stabilized by binding to this unit. This finding suggests that the first beta alpha beta unit could function as an N-terminal folding nucleus, upon which the remainder of the polypeptide chain can be assembled. Indeed, simulations with longer units (beta-alpha-beta-alpha and beta-alpha-beta-alpha beta-beta) showed that all structural elements of these units are rather stable. The outcome of our studies is in line with suggestions that folding of the N-terminal portion of LDH in vivo can be a cotranslational process that takes place during the ribosomal peptide synthesis.

Native-like secondary structure in a peptide from the alpha-domain of hen lysozyme

BACKGROUND

To gain insight into the local and nonlocal interactions that contribute to the stability of hen lysozyme, we have synthesized two peptides that together comprise the entire alpha-domain of the protein. One peptide (peptide 1-40) corresponds to the sequence that forms two alpha-helices, a loop region, and a small beta-sheet in the N-terminal region of the native protein. The other (peptide 84-129) makes up the C-terminal part of the alpha-domain and encompasses two alpha-helices and a 3(10) helix in the native protein.

RESULTS

As judged by CD and a range of NMR parameters, peptide 1-40 has little secondary structure in aqueous solution and only a small number of local hydrophobic interactions, largely in the loop region. Peptide 84-129, by contrast, contains significant helical structure and is partially hydrophobically collapsed. More specifically, the region corresponding to helix C in native lysozyme is disordered, whereas regions corresponding to the D and 3(10) helices in the native protein are helical in this peptide. The structure in peptide 84-129 is at least partly stabilized by interactions between residues in the two helical regions, as suggested by further NMR analysis of three short peptides corresponding to the individual helices in this region of the native protein.

CONCLUSIONS

Stabilization of structure in the sequence 1-40 appears to be facilitated predominantly by long-range interactions between this region and the sequence 84-129. In native lysozyme, the existence of two disulphide bonds between the N- and C-terminal halves of the alpha-domain is likely to be a major factor in their stabilization. The data show, however, that native-like secondary structure can be generated in the C-terminal portion of the alpha-domain by nonspecific and nonnative interactions within a partially collapsed state.

Folding of the squash trypsin inhibitor EETI II. Evidence of native and non-native local structural preferences in a linear analogue

A peptide, corresponding to the entire sequence of the squash trypsin inhibitor EETI II (Ecballium elaterium trypsin inhibitor) in which the six cysteines, engaged in three disulphide bridges in native EETI II, have been replaced by six serines, has been synthesised. CD, Fourier-transform infrared spectroscopy (FTIR) and 1H-NMR studies of this peptide revealed that some secondary structures present in native EETI II are still populated in the absence of disulphide bonds. Native-like secondary structures were observed for segments 10-15 (helix), 16-19 and 22-25 (reverse turns) but no native tertiary interaction was detected. However, a non-native local interaction between the aromatic ring of Phe26 and the amide group of Gly28 was observed. It is hypothesised that the 10-15, 16-19 and 22-25 native-like local conformations could play a major role in the early folding of EETI II.

Effect of pH on the conformation and backbone dynamics of a 27-residue peptide in trifluoroethanol. An NMR and CD Study

The C-terminal fragment, residues 385-411, from human fibrinogen gamma-chain, i.e. KIIPFNRLTIGEGQQHHLG-GAKQAGDV, shows multiple turn conformations in aqueous solution (Mayo, K. H., Burke, C., Lindon, J. N., and Kloczewiak, M. (1990) Biochemistry 29, 3277-3286). The present study investigates the effect of pH and trifluoroethanol on the conformation and backbone dynamics of this 27-residue peptide. Both circular dichroism (CD) and 1H-NMR data indicate the normally observed increased helical conformations as a function of increasing trifluoroethanol. 1H-NMR structural studies done in the presence of 40% trifluoroethanol, pH 5.3, yield a network of nuclear Overhauser effects consistent with significant populations of helix-like conformation. Distance geometry calculations based on nuclear Overhauser effect-derived distance constraints yield a family of structures with relatively well defined N- and C-terminal conformations and an ill defined mid-peptide region from Gly397 to Gly403. Similar conformational populations are observed at pH 2.5. CD studies, however, indicate an increase in average alpha-helix content on decreasing the pH from 6 to 2. This apparent conflict between CD and NMR results may be explained by a transition from multiple beta-turn character at pH 5.3 to increased alpha-helix structure at pH 2.5. 13C alpha NMR relaxation data analyzed with the Lipari-Szabo model-free approach provide order parameters that demonstrate little if any influence of pH on backbone motional restrictions within the more flexible mid-peptide domain. At low pH, however, motions become less restricted within N-terminal residues Lys385-Phe389 and more restricted within C-terminal residues Ala405-Val411.