Equilibrium structure and folding of a helix-forming peptide: circular dichroism measurements and replica-exchange molecular dynamics simulations

We have performed experimental measurements and computer simulations of the equilibrium structure and folding of a 21-residue alpha-helical heteropeptide. Far ultraviolet circular dichroism spectroscopy is used to identify the presence of helical structure and to measure the thermal unfolding curve. The observed melting temperature is 296 K, with a folding enthalpy of -11.6 kcal/mol and entropy of -39.6 cal/(mol K). Our simulations involve 45 ns of replica-exchange molecular dynamics of the peptide, using eight replicas at temperatures between 280 and 450 K, and the program CHARMM with a continuum solvent model. In a 30-ns simulation started from a helical structure, conformational equilibrium at all temperatures was reached after 15 ns. This simulation was used to calculate the peptide melting curve, predicting a folding transition with a melting temperature in the 330-350 K range, enthalpy change of -10 kcal/mol, and entropy change of -30 cal/(mol K). The simulation results were also used to analyze the peptide structural fluctuations and the free-energy surface of helix unfolding. In a separate 15-ns replica-exchange molecular dynamics simulation started from the extended structure, the helical conformation was first attained after approximately 2.8 ns, and equilibrium was reached after 10 ns of simulation. These results showed a sequential folding process with a systematic increase in the number of hydrogen bonds until the helical state is reached, and confirmed that the alpha-helical state is the global free-energy minimum for the peptide at low temperatures.     

Unfolding of a leucine zipper is not a simple two-state transition

Temperature-induced unfolding of the leucine zipper, an alpha-helical, double-stranded, coiled-coil, was studied by circular dichroism spectroscopy, spectrofluorimetry and heat capacity scanning calorimetry. It is shown that this process does not represent a simple two-state transition, as previously believed, but consists of several stages. The first transition starts at the very beginning of heating from 0 degrees C and proceeds with significant heat absorption and decrease of ellipticity. This transition does not depend on the concentration of protein and is sensitive to modification of the N terminus; it is therefore associated with unfolding or fraying of this part of the leucine zipper. The second transition takes place at a considerably higher temperature; it is more pronounced than the first one and does not depend on the concentration of protein, i.e. it is unimolecular. This transition is sensitive to modification of both termini of the leucine zipper and affects the optical properties of a tryptophan residue placed in the central part of the zipper. It therefore involves the whole dimer but does not result in its dissociation, presumably being associated with some repacking of the coiled-coil. This second transition is followed at higher temperatures by the concentration-dependent cooperative unfolding/dissociation of the two strands. The enthalpy and entropy of the temperature-induced structural changes of the leucine zipper that take place before its cooperative unfolding/dissociation comprises almost 40% of the total enthalpy and entropy of unfolding of the completely folded coiled-coil, the state in which it appears to be below 0 degrees C. Comparison of the total enthalpy of leucine zipper unfolding with that of a single-stranded alpha-helix shows that their temperature-dependence correlates with the extent of intramolecular non-polar contacts and allows an assessment of the enthalpy of hydrogen bonding in alpha-helices, which appears to be about 3.3kJmol(-1) at 20 degrees C.     

Role of backbone hydration and salt-bridge formation in stability of alpha-helix in solution

We test molecular level hypotheses for the high thermal stability of alpha-helical conformations of alanine-based peptides by performing detailed atomistic simulations of a 20-amino-acid peptide with explicit treatment of water. To assess the contribution of large side chains to alpha-helix stability through backbone desolvation and salt-bridge formation, we simulate the alanine-rich peptide, Ac-YAEAAKAAEAAKAAEAAKAF-Nme, referred to as the EK peptide, that has three pairs of "i, i + 3" glutamic acid(-) and lysine(+) substitutions. Efficient configurational sampling of the EK peptide over a wide temperature range enabled by the replica exchange molecular dynamics technique allows characterization of the stability of alpha-helix with respect to heat-induced unfolding. We find that near ambient temperatures, the EK peptide predominately samples alpha-helical configurations with 80% fractional helicity at 300 K. The helix melts over a broad range of temperatures with melting temperature, T(m), equal to 350 K, that is significantly higher than the T(m) of a 21-residue polyalanine peptide, A(21). Salt-bridges between oppositely charged Glu(-) and Lys(+) side chains can, in principle, provide thermal stability to alpha-helical conformers. For the specific EK peptide sequence, we observe infrequent formation of Glu-Lys salt-bridges (with approximately 10-20% probability) and therefore we conclude that salt-bridge formation does not contribute significantly to the EK peptide's helical stability. However, lysine side chains are found to shield specific "i, i + 4" backbone hydrogen bonds from water, indicating that large side-chain substituents can play an important role in stabilizing alpha-helical configurations of short peptides in aqueous solution through mediation of water access to backbone hydrogen bonds. These observations have implications on molecular engineering of peptides and biomolecules in the design of their thermostable variants where the shielding mechanism can act in concert with other factors such as salt-bridge formation, thereby increasing thermal stability considerably.     

Competing interactions contributing to alpha-helical stability in aqueous solution.

The stability of a 15-residue peptide has been investigated using CD spectroscopy and molecular simulation techniques. The sequence of the peptide was designed to include key features that are known to stabilize alpha-helices, including ion pairs, helix dipole capping, peptide bond capping, and aromatic interactions. The degree of helicity has been determined experimentally by CD in three solvents (aqueous buffer, methanol, and trifluoroethanol) and at two temperatures. Simulations of the peptide in the aqueous system have been performed over 500 ps at the same two temperatures using a fully explicit solvent model. Consistent with the CD data, the degree of helicity is decreased at the higher temperature. Our analysis of the simulation results has focused on competition between different side-chain/side-chain and side-chain/main-chain interactions, which can, in principle, stabilize the helix. The unfolding in aqueous solution occurs at the amino terminus because the side-chain interactions are insufficient to stabilize both the helix dipole and the peptide hydrogen bonds. Loss of capping of the peptide backbone leads to water insertion within the first peptide hydrogen bond and hence unfolding. In contrast, the carboxy terminus of the alpha-helix is stable in both simulations because the C-terminal lysine residue stabilizes the helix dipole, but at the expense of an ion pair.     

Molecular dynamics simulations of helix denaturation.

An understanding of the structural transitions that an alpha-helix undergoes will help to elucidate such motions in proteins and their role in protein folding. We present the results of molecular dynamics simulations to investigate these transitions in a short polyalanine peptide (13 residues) both in vacuo and in the presence of solvent. The denaturation of this peptide was monitored as a function of temperature (ranging from 5 to 200 degrees C). In vacuo, the helical state predominated at all temperatures, whereas in solution the helix melted with increasing temperature. The peptide was predominantly helical at low temperature in solution, while at intermediate temperatures the peptide spent the bulk of the time fluctuating between different conformations with intermediate amounts of helix, e.g. not completely helical nor entirely non-helical. Many of these conformations consisted of short helical segments with intervening non-helical residues. At high temperature the peptide unfolded and adopted various collapsed unstructured states. The intrahelical hydrogen bonds that break at high temperature were not fully compensated by hydrogen bonds with water molecules in the partially unfolded forms of the peptide. Increases in temperature disrupted both the helical structure and the peptide-water interactions. Water played a major but indirect role in facilitating unfolding, as opposed to specifically competing for the intrapeptide hydrogen bonds. The implications of our results to protein folding are discussed.     

Role of water on unfolding kinetics of helical peptides studied by molecular dynamics simulations

Molecular dynamics simulations have been carried out with four polypeptides, Ala13, Val(13), Ser13, and Ala4Gly5Ala4, in vacuo and with explicit hydration. The unfolding of the polypeptides, which are initially fully alpha-helix in conformation, has been monitored during trajectories of 0.3 ns at 350 K. A rank of Ala < Val < Ser < Gly is found in the order of increasing rate of unwinding. The unfolding of Ala13 and Val(13) is completed in hundreds of picoseconds, while that of Ser13 is about one order of magnitude faster. The helix content of the peptide containing glycine residues falls to zero within a few picoseconds. Ramachandran plots indicate quite distinct equilibrium distributions and time evolution of dihedral angles in water and in vacuum for each residue type. The unfolding of polyalanine and polyvaline helices is accelerated due to solvation. In contrast, polyserine is more stable in water compared to vacuum, because its side chains can form intramolecular hydrogen bonds with the backbone more readily in vacuum, which disrupts the helix. Distribution functions of the spatial and angular position of water molecules in the proximity of the polypeptide backbone polar groups reveal the stabilization of the coiled structures by hydration. The transition from helix to coil is characterized by the appearance of a new peak in the probability distribution at a specific location characteristic of hydrogen bond formation between water and backbone polar groups. No significant insertion of water molecules is observed at the precise onset of unwinding, while (i, i+3) hydrogen bond formation is frequently detected at the initiation of alpha-helix unwinding.     

Helix-forming tendencies of amino acids depend on the restrictions of side-chain rotamer conformations: crystal structure of the tripeptide GAI in two crystalline forms.

In our attempts to design crystalline alpha-helical peptides, we synthesized and crystallized GAI (C11H21N3O4) in two crystal forms, GAI1 and GAI2. Form 1 (GAI1) Gly-L-Ala-L-Ile (C11H21N3O4.3H2O) crystals are monoclinic, space group P2(1) with a = 8.171(2), b = 6.072(4), c = 16.443(4) A, beta = 101.24(2) degrees, V = 800 A3, Dc = 1.300 g cm-3 and Z = 2, R = 0.081 for 482 reflections. Form 2 (GAI2) Gly-L-Ala-L-Ile (C11H21N3O4.1/2H2O) is triclinic, space group P1 with a = 5.830(1), b = 8.832(2), c = 15.008(2) A, alpha = 102.88(1), beta = 101.16(2), gamma = 70.72(2) degrees, V = 705 A3, Z = 2, Dc = 1.264 g cm-3, R = 0.04 for 2582 reflections. GAI1 is isomorphous with GAV and forms a helix, whereas GAI2 does not. In GAI1, the tripeptide molecule is held in a near helical conformation by a water molecule that bridges the NH3+ and COO- groups, and acts as the fourth residue needed to complete the turn by forming two hydrogen bonds. Two other water molecules form intermolecular hydrogen bonds in stabilizing the helical structure so that the end result is a column of molecules that looks like an incipient alpha-helix. GAI2 imitates a cyclic peptide and traps a water molecule. The conformation angles chi 11 and chi 12 for the side chain are (-63.7 degrees, 171.1 degrees) for the helical GAI1, and (-65.1 degrees, 58.6 degrees) and (-65.0 degrees, 58.9 degrees) for the two independent nonhelical molecules in GAI2; in GAI1, both the C gamma atoms point away from the helix, whereas in GAI2 the C gamma atom with the g+ conformation points inward to the helix and causes sterical interaction with atoms in the adjacent peptide plane. From these results, it is clear that the helix-forming tendencies of amino acids correlate with the restrictions of side-chain rotamer conformations. Both the peptide units in GAI1 are trans and show significant deviation from planarity [omega 1 = -168(1) degrees; omega 2 = -171(1) degrees] whereas both the peptide units in both the molecules A and B in GAI2 do not show significant deviation from planarity [omega 1 = 179.3(3) degrees; omega 2 = -179.3(3) degrees for molecule A and omega 1 = 179.5(3) degrees; omega 2 = -179.4(3) degrees for molecule B], indicating that the peptide planes in these incipient alpha-helical peptides are considerably bent.     

Thermodynamics of the alpha-helix-coil transition of amphipathic peptides in a membrane environment: implications for the peptide-membrane binding equilibrium

Amphipathic alpha-helices are the membrane binding motif in many proteins. The corresponding peptides are often random coil in solution but are folded into an alpha-helix upon interaction with the membrane. The energetics of this ubiquitous folding process are still a matter of conjecture. Here, we present a new method to quantitatively analyze the thermodynamics of peptide folding at the membrane interface. We have systematically varied the helix content of a given amphipathic peptide when bound to the membrane and have correlated the thermodynamic binding parameters determined by isothermal titration calorimetry with the alpha-helix content obtained by circular dichroism spectroscopy. The peptides investigated were the antibiotic magainin 2 amide and three analogs in which two adjacent amino acid residues were substituted by their d-enantiomers. The thermodynamic parameters controlling the alpha-helix formation were found to be linearly related to the helicity of the membrane-bound peptides. Helix formation at the membrane surface is characterized by an enthalpy change of DeltaH(helix) approximately -0.7 kcal/mol per residue, an entropy change of DeltaS(helix) approximately -1.9 cal/molK residue and a free energy change of DeltaG(helix)=-0.14 kcal/mol residue. Helix formation is a strong driving force of peptide insertion into the membrane and accounts for about 50 % of the free energy of binding. An increase in temperature entails an unfolding of the membrane-bound helix. The temperature dependence can be described with the Zimm-Bragg theory and the enthalpy of unfolding agrees with that deduced from isothermal titration calorimetry.     

Energetic contribution of solvent-exposed ion pairs to alpha-helix structure.

Understanding the role of amino acid side-chain interactions in forming secondary structure in proteins is useful for deciphering how proteins fold and for predicting folded structures of proteins from their sequence. Analysis of the secondary structure as a function of pH in two designed synthetic peptides with identical composition but different sequences, affords a quantitative estimate of the free energy contribution of a single ion pair to the stability of an isolated alpha-helix. One peptide contains repeated blocks of Glu4Lys4. The second has repeated blocks of Glu2Lys2. The former contains significant helical structure at neutral pH while the latter has none, based on ultraviolet light circular dichroism measurements and 1H nuclear magnetic resonance spectroscopy. The difference is attributed to formation of helix-stabilizing salt-bridges between Glu- and Lys+ spaced at i, i + 4 intervals in the former peptide. The free energy of formation of a single Glu(-)-Lys+ salt-bridge can be evaluated by using a statistical model of the helix-coil transition that explicitly includes salt-bridges: the result is -0.50(+/- 0.05) kcal/mol at 4 degrees C and neutral pH in 10 mM salt, in agreement with a value derived for a single salt-bridge in a helix on the surface of a globular protein.     

Increase in the stability and helical content of estrogen receptor alpha in the presence of the estrogen response element: analysis by circular dichroism spectroscopy

Ligand-dependent stabilization of the estrogen receptor (ER) is often postulated, with limited support from experimental data. We studied the thermal unfolding of recombinant ERalpha by circular dichroism (CD) spectroscopy. The T(M) of unfolding of ERalpha was 38 +/- 2.4 degrees C, and the van't Hoff enthalpy of unfolding was 31.7 +/- 3.4 kcal/mol in the absence of ligands. Addition of estradiol (E(2)) increased the T(M) to 43.6 +/- 2.3 degrees C, while addition of E(2) and an oligonucleotide harboring the estrogen response element (ERE) increased the T(M) to 47.9 +/- 1.6 degrees C. Addition of the antiestrogen 4-hydroxytamoxifen (HT) alone did not increase the T(M); however, a combination of HT and the ERE increased the T(M) to 48.9 +/- 1.0 degrees C. The ERE alone increased the T(M) to 46.1 +/- 0.9 degrees C. Addition of E(2) alone had no effect on the apparent enthalpy of unfolding; however, the ERE alone increased the apparent enthalpy from 31.7 to 36.1 kcal/mol. ERalpha samples containing the ERE also exhibited an increase in the negative ellipticity at 208 and 222 nm, relative to that of ligand-free ERalpha, suggesting a stabilization of the alpha-helix. CD data analysis further showed that the presence of the ERE caused a large increase in alpha-helical content of ERalpha in both the presence and absence of the ligands. This increase in alpha-helical content of ERalpha was not observed in the presence of a nonspecific oligonucleotide. These results show that the ERE can increase the thermal stability of ERalpha, enhance its alpha-helical content, and facilitate the cooperativity of the folding transition.     

Thermodynamic consequences of burial of polar and non-polar amino acid residues in the protein interior

Effects of amino acid substitutions at four fully buried sites of the ubiquitin molecule on the thermodynamic parameters (enthalpy, Gibbs energy) of unfolding were evaluated experimentally using differential scanning calorimetry. The same set of substitutions has been incorporated at each of four sites. These substitutions have been designed to perturb packing (van der Waals) interactions, hydration, and/or hydrogen bonding. From the analysis of the thermodynamic parameters for these ubiquitin variants we conclude that: (i) packing of non-polar groups in the protein interior is favorable and is largely defined by a favorable enthalpy of van der Waals interactions. The removal of one methylene group from the protein interior will destabilize a protein by approximately 5 kJ/mol, and will decrease the enthalpy of a protein by 12 kJ/mol. (ii) Burial of polar groups in the non-polar interior of a protein is highly destabilizing, and the degree of destabilization depends on the relative polarity of this group. For example, burial of Thr side-chain in the non-polar interior will be less destabilizing than burial of Asn side-chain. This decrease in stability is defined by a large enthalpy of dehydration of polar groups upon burial. (iii) The destabilizing effect of dehydration of polar groups upon burial can be compensated if these buried polar groups form hydrogen bonding. The enthalpy of this hydrogen bonding will compensate for the unfavorable dehydration energy and as a result the effect will be energetically neutral or even slightly stabilizing.     

Molecular dynamics simulations of the unfolding of an alpha-helical analogue of ribonuclease A S-peptide in water

Molecular dynamics simulations of the S-peptide analogue AETAAAKFLREHMDS have been conducted in aqueous solution for 300 ps at 278 K and for 500 ps in two different runs at 358 K. The results show agreement with experimental observations in that at low temperature, 5 degrees C, the helix is stable, while unfolding is observed at 85 degrees C. In the low-temperature simulation a solvent-separated ion pair was formed between Glu-2 and Arg-10, and the side chain of His-12 reoriented toward the C-terminal end of the alpha-helix. Detailed analyses of the unfolding pathways at high temperature have also revealed that the formation or disappearance of main-chain helical hydrogen bonds occurs frequently through an alpha in equilibrium with 3(10) in equilibrium with no hydrogen bond sequence.     

Protein recognition motifs: design of peptidomimetics of helix surfaces

Helices represent one of the most common recognition motifs in proteins. The design of nonpeptidic scaffolds, such as the 3,2',2'-tris-substituted terphenyl, that can imitate the side-chain orientation along one face of an alpha-helix potentially provides an effective means to modulate helix-recognition functions. Here, based on theoretical arguments, we described novel alpha-helix mimetics which are more effective than the terphenyl at constraining the aryl-aryl torsion angles to those associated with structures suitable for mimicking the alpha-helical twist for side-chain orientation and for superimposing the side chains of residues i, i + 3 or i + 4, i + 7 when compared with the alpha-beta side-chain vectors of the regular alpha-helix with an improved root mean square deviation (RMSD) of approximately 0.5 A. In addition, this study suggests that rotamer distributions around the C(alpha)--C(beta) bonds of these helix mimetics are similar to those of alpha-helices, except that these rotamer distributions show an approximately 60 degrees shift compared to those of alpha-helices when the mimetic axis is superimposed upon the helix axis. This change in rotamer orientation complicates mimicry of the helix surface.     

Differential stability of the triple helix of (Pro-Pro-Gly)10 in H2O and D2O: thermodynamic and structural explanations

(Pro-Pro-Gly)10 [(PPG10)], a collagen-like polypeptide, forms a triple-helical, polyproline-II structure in aqueous solution at temperatures somewhat lower than physiological, with a melting temperature of 24.5 degrees C. In this article, we present circular dichroism spectra that demonstrate an increase of the melting temperature with the addition of increasing amounts of D2O to an H2O solution of (PPG)10, with the melting temperature reaching 40 degrees C in pure D2O. A thermodynamic analysis of the data demonstrates that this result is due to an increasing enthalpy of unfolding in D2O vs. H2O. To provide a theoretical explanation for this result, we have used a model for hydration of (PPG)10 that we developed previously, in which inter-chain water bridges are formed between sterically crowded waters and peptide bond carbonyls. Energy minimizations were performed upon this model using hydrogen bond parameters for water, and altered hydrogen bond parameters that reproduced the differences in carbonyl oxygen-water oxygen distances found in small-molecule crystal structures containing oxygen-oxygen hydrogen bonds between organic molecules and H2O or D2O. It was found that using hydrogen bond parameters that reproduced the distance typical of hydrogen bonds to D2O resulted in a significant lowering of the potential energy of hydrated (PPG)10. This lowering of the energy involved energetic terms that were only indirectly related to the altered hydrogen bond parameters, and were therefore not artifactual; the intra-(PPG10) energy, plus the water-(PPG10) van der Waals energy (not including hydrogen bond interactions), were lowered enough to qualitatively account for the lower enthalpy of the triple-helical conformation, relative to the unfolded state, in D2O vs. H2O. This result indicates that the geometry of the carbonyl-D2O hydrogen bonds allows formation of good hydrogen bonds without making as much of an energetic sacrifice from other factors as in the case of hydration by H2O.     

Conformational intermediates in the folding of a coiled-coil model peptide of the N-terminus of tropomyosin and alpha alpha-tropomyosin.

Circular dichroism was used to study the folding of alpha alpha-tropomyosin and AcTM43, a 43-residue peptide designed to serve as a model for the N-terminal domain of tropomyosin. The sequence of the peptide is AcMDAIKKKMQMLKLDVENLLDRLEQLEADLKALEDRYKQLEGGC. The peptide appeared to form a coiled coil at low temperatures (< 25 degrees C) in buffers with physiological ionic strength and pH. The folding and unfolding of the peptide, however, were noncooperative. When CD spectra were examined as a function of temperature, the apparent degree of folding differed when the ellipticity was followed at 222, 208, and 280 nm. Deconvolution of the spectra suggested that at least three component curves contributed to the CD in the far UV. One component curve was similar to the CD spectrum of the coiled-coil alpha-helix of native alpha alpha-tropomyosin. The second curve resembled the spectrum of single-stranded short alpha-helical segments found in globular proteins. The third was similar to that of polypeptides in the random coil conformation. These results suggested that as the peptide folded, the alpha-helical content increased before most of the coiled coil was formed. When the CD spectrum of striated muscle alpha alpha-tropomyosin was examined as a function of temperature, the unfolding was also not totally cooperative. As the temperature was raised from 0 to 25 degrees C, there was a decrease in the coiled coil and an increase in the conventional alpha-helix type spectrum without formation of random coil. The major transition, occurring at 40 degrees C, was a cooperative transition characterized by the loss of all of the remaining coiled coil and a concomitant increase in random coil.     

Helix nucleation kinetics from molecular simulations in explicit solvent

We study the reversible folding/unfolding of short Ala and Gly-based peptides by molecular dynamics simulations of all-atom models in explicit water solvent. A kinetic analysis shows that the formation of a first alpha-helical turn occurs within 0.1-1 ns, in agreement with the analyses of laser temperature jump experiments. The unfolding times exhibit Arrhenius temperature dependence. For a rapidly nucleating all-Ala peptide, the helix nucleation time depends only weakly on temperature. For a peptide with enthalpically competing turn-like structures, helix nucleation exhibits an Arrhenius temperature dependence, corresponding to the unfolding of enthalpic traps in the coil ensemble. An analysis of structures in a "transition-state ensemble" shows that helix-to-coil transitions occur predominantly through breaking of hydrogen bonds at the helix ends, particularly at the C-terminus. The temperature dependence of the transition-state ensemble and the corresponding folding/unfolding pathways illustrate that folding mechanisms can change with temperature, possibly complicating the interpretation of high-temperature unfolding simulations. The timescale of helix formation is an essential factor in molecular models of protein folding. The rapid helix nucleation observed here suggests that transient helices form early in the folding event.     

Thermodynamics of barnase unfolding.

The thermodynamics of barnase denaturation has been studied calorimetrically over a broad range of temperature and pH. It is shown that in acidic solutions the heat denaturation of barnase is well approximated by a 2-state transition. The heat denaturation of barnase proceeds with a significant increase of heat capacity, which determines the temperature dependencies of the enthalpy and entropy of its denaturation. The partial specific heat capacity of denatured barnase is very close to that expected for the completely unfolded protein. The specific denaturation enthalpy value extrapolated to 130 degrees C is also close to the value expected for the full unfolding. Therefore, the calorimetrically determined thermodynamic characteristics of barnase denaturation can be considered as characteristics of its complete unfolding and can be correlated with structural features--the number of hydrogen bonds, extent of van der Waals contacts, and the surface areas of polar and nonpolar groups. Using this information and thermodynamic information on transfer of protein groups into water, the contribution of various factors to the stabilization of the native structure of barnase has been estimated. The main contributors to the stabilization of the native state of barnase appear to be intramolecular hydrogen bonds. The contributions of van der Waals interactions between nonpolar groups and those of hydration effects of these groups are not as large if considered separately, but the combination of these 2 factors, known as hydrophobic interactions, is of the same order of magnitude as the contribution of hydrogen bonding.     

Energetics of protein folding

The energetics of protein folding determine the 3D structure of a folded protein. Knowledge of the energetics is needed to predict the 3D structure from the amino acid sequence or to modify the structure by protein engineering. Recent developments are discussed: major factors are reviewed and auxiliary factors are discussed briefly. Major factors include the hydrophobic factor (burial of non-polar surface area) and van der Waals interactions together with peptide hydrogen bonds and peptide solvation. The long-standing model for the hydrophobic factor (free energy change proportional to buried non-polar surface area) is contrasted with the packing-desolvation model and the approximate nature of the proportionality between free energy and apolar surface area is discussed. Recent energetic studies of forming peptide hydrogen bonds (gas phase) are reviewed together with studies of peptide solvation in solution. Closer agreement is achieved between the 1995 values for protein unfolding enthalpies in vacuum given by Lazaridis-Archontis-Karplus and Makhatadze-Privalov when the solvation enthalpy of the peptide group is taken from electrostatic calculations. Auxiliary factors in folding energetics include salt bridges and side-chain hydrogen bonds, disulfide bridges, and propensities to form alpha-helices and beta-structure. Backbone conformational entropy is a major energetic factor which is discussed only briefly for lack of knowledge.     

Thermal unfolding of myosin rod and light meromyosin: circular dichroism and tryptophan fluorescence studies

Rabbit skeletal myosin rod, which is the coiled-coil alpha-helical portion of myosin, contains two tryptophan residues located in the light meromyosin (LMM) portion whose fluorescence contributes 27% to the fluorescence of the entire myosin molecule. The temperature dependence of several fluorescence parameters (quantum yield, spectral position, polarization) of the rod and its LMM portion was compared to the thermal unfolding of the helix measured with circular dichroism. Rod unfolds with three major helix unfolding transitions: at 43, 47, and 53 degrees C, with the 43 and 53 degrees C transitions mainly located in the LMM region and the 47 degrees C transition mainly located in the subfragment 2 region. The fluorescence study showed that the 43 degrees C transition does not involve the tryptophan-containing region and that the 47 degrees C transition produces an intermediate with different fluorescence properties from both the completely helical and fully unfolded states. That is, although the fluorescence of the 47 degrees C intermediate is markedly quenched, the tryptophyl residues do not become appreciably exposed to solvent until the 53 degrees C transition. It is suggested that although the intermediate that is formed in the 47 degrees C transition contains an extensive region which is devoid of alpha-helix, the unfolded region is not appreciably solvated or flexible. It appears to have the properties of a collapsed nonhelical state rather than a classical random coil.     

Reversible unfolding of beta-sheets in membranes: a calorimetric study

The hexapeptide acetyl-Trp-Leu(5) (AcWL(5)) has the remarkable ability to assemble reversibly and spontaneously into beta-sheets on lipid membranes as a result of monomer partitioning followed by cooperative assembly. This system provides a unique opportunity to study the thermodynamics of protein folding in membranes, which we have done using isothermal titration calorimetry (ITC) and differential scanning calorimetry (DSC). The results, which may represent the first example of reversible thermal unfolding of peptides in membranes, help to define the contribution of hydrogen bonding to the extreme thermal stability of membrane proteins. ITC revealed that the enthalpy change for partitioning of monomeric, unstructured AcWL(5) from water into membranes was zero within experimental error over the temperature range of 5 degrees C to 75 degrees C. DSC showed that the beta-sheet aggregates underwent a reversible, endothermic, and very asymmetric thermal transition with a concentration-dependent transition temperature (T(m)) in the range of 60 degrees C to 80 degrees C. A numerical model of nucleation and growth-dependent assembly of oligomeric beta-sheets, proposed earlier to describe beta-sheet formation in membranes, recreated remarkably well the unusual shape and concentration-dependence of the transition peaks. The enthalpy for thermal unfolding of AcWL(5) beta-sheets in the membrane was found to be about 8(+/-1)kcal mol(-1), or about 1.3(+/-0.2)kcal mol(-1) per residue.