Folding-unfolding thermodynamics of a beta-heptapeptide from equilibrium simulations

The thermodynamics of folding and unfolding of a beta-heptapeptide in methanol solution has been studied at four different temperatures, 298 K, 340 K, 350 K, and 360 K, by molecular dynamics simulation. At each of these temperatures, the 50-ns simulations were sufficient to generate an equilibrium distribution between a relatively small number of conformations (approximately 10(2)), showing that, even above the melting temperature (approximately 340 K), the peptide does not randomly sample conformational space. The free energy of folding and the free energy difference between pairs of conformations have been calculated from their relative populations. The experimentally determined folded conformation at 298 K, a left-handed 3(1)-helix, is at each of the four temperatures the predominant conformation, with its probability and average lifetime decreasing with increasing temperature. The most common intermediates of folding and unfolding are also the same at the four temperatures. Paths and rates of interconversion between different conformations have been determined. It has been found that folding can occur through multiple pathways, not necessarily downhill in free energy, although the final step involves a reduced number of intermediates.     

Folding mechanism of FIS, the intertwined, dimeric factor for inversion stimulation

FIS, the factor for inversion stimulation, from Escherichia coli and other enteric bacteria, is an interwined alpha-helical homodimer. Size exclusion chromatography and static light scattering measurements demonstrated that FIS is predominately a stable dimer at the concentrations (1-10 microM monomer) and buffer conditions employed in this study. The folding and unfolding of FIS were studied with both equilibrium and kinetic methods by circular dichroism using urea and guanidinium chloride (GdmCl) as the perturbants. The equilibrium folding is reversible and well-described by a two-state folding model, with stabilities at 10 degrees C of 15.2 kcal mol(-1) in urea and 13.5 kcal mol(-1) in GdmCl. The kinetic data are consistent with a two-step folding reaction where the two unfolded monomers associate to a dimeric intermediate within the mixing time for the stopped-flow instrument (<5 ms), and a slower, subsequent folding of the dimeric intermediate to the native dimer. Fits of the burst phase amplitudes as a function of denaturant showed that the free energy for the formation of the dimeric intermediate constitutes the majority of the stability of the folding (9.6 kcal mol(-1) in urea and 10.5 kcal mol(-1) in GdmCl). Folding-to-unfolding double jump kinetic experiments were also performed to monitor the formation of native dimer as a function of folding delay times. The data here demonstrate that the dimeric intermediate is obligatory and on-pathway. The folding mechanism of FIS, when compared to other intertwined, alpha-helical, homodimers, suggests that a transient kinetic dimeric intermediate may be a common feature of the folding of intertwined, segment-swapped, alpha-helical dimers.     

High- and low-temperature unfolding of human high-density apolipoprotein A-2.

Human plasma apolipoprotein A-2 (apoA-2) is the second major protein of the high-density lipoproteins that mediate the transport and metabolism of cholesterol. Using CD spectroscopy and differential scanning calorimetry, we demonstrate that the structure of lipid-free apoA-2 in neutral low-salt solutions is most stable at approximately 25 degrees C and unfolds reversibly both upon heating and cooling from 25 degrees C. High-temperature unfolding of apoA-2, monitored by far-UV CD, extends from 25-85 degrees C with midpoint Th = 56 +/- 2 degrees C and vant Hoff's enthalpy delta H(Th) = 17 +/- 2 kcal/mol that is substantially lower than the expected enthalpy of melting of the alpha-helical structure. This suggests low-cooperativity apoA-2 unfolding. The apparent free energy of apoA-2 stabilization inferred from the CD analysis of the thermal unfolding, delta G(app)(25 degrees) = 0.82 +/- 0.15 kcal/mol, agrees with the value determined from chemical denaturation. Enhanced low-temperature stability of apoA-2 observed upon increase in Na2HPO4 concentration from 0.3 mM to 50 mM or addition of 10% glycerol may be linked to reduced water activity. The close proximity of the heat and cold unfolding transitions, that is consistent with low delta G(app)(25 degrees), indicates that lipid-free apoA-2 has a substantial hydrophobic core but is only marginally stable under near-physiological solvent conditions. This suggests that in vivo apoA-2 transfer is unlikely to proceed via the lipid-free state. Low delta H(Th) and low apparent delta Cp approximately 0.52 kcal/mol.K inferred from the far-UV CD analysis of apoA-2 unfolding, and absence of tertiary packing interactions involving Tyr groups suggested by near-UV CD, are consistent with a molten globular-like state of lipid-free apoA-2.     

Thermal unfolding and conformational stability of the recombinant domain II of glutamate dehydrogenase from the hyperthermophile Thermotoga maritima

Domain II (residues 189-338, M(r) = 16 222) of glutamate dehydrogenase from the hyperthermophilic bacterium Thermotoga maritima was used as a model system to study reversible unfolding thermodynamics of this hyperthermostable enzyme. The protein was produced in large quantities in E.COLI: using a T7 expression system. It was shown that the recombinant domain is monomeric in solution and that it comprises secondary structural elements similar to those observed in the crystal structure of the hexameric enzyme.The recombinant domain is thermostable and undergoes reversible and cooperative thermal unfolding in the pH range 5.90-8.00 with melting temperatures between 75.1 and 68.0 degrees C. Thermal unfolding of the protein was studied using differential scanning calorimetry and circular dichroism spectroscopy. Both methods yielded comparable values. The analysis revealed an unfolding enthalpy at 70 degrees C of 70.2 +/- 4.0 kcal/mol and a DeltaC(p) value of 1.4 +/- 0.3 kcal/mol K. Chemical unfolding of the recombinant domain resulted in m values of 3.36 +/- 0.10 kcal/mol M for unfolding in guanidinium chloride and 1.46 +/- 0.04 kcal/mol M in urea. The thermodynamic parameters for thermal and chemical unfolding equilibria indicate that domain II from T.MARITIMA: glutamate dehydrogenase is a thermostable protein with a DeltaG(max) of 3.70 kcal/mol. However, the thermal and chemical stabilities of the domain are lower than those of the hexameric protein, indicating that interdomain interactions must play a significant role in the stabilization of T. MARITIMA: domain II glutamate dehydrogenase.     

Multifrequency calorimetry of the folding/unfolding transition of cytochrome c.

The folding-unfolding transition of Fe(III) cytochrome c has been studied with the new technique of multifrequency calorimetry. Multifrequency calorimetry is aimed at measuring directly the dynamics of the energetic events that take place during a thermally induced transition by measuring the frequency dispersion of the heat capacity. This is done by modulating the folding/unfolding equilibrium using a variable frequency, small oscillatory temperature perturbation (approximately 0.05-0.1 degrees C) centered at the equilibrium temperature of the system. Fe(III) cytochrome c at pH 4 undergoes a fully reversible folding/unfolding transition centered at 67.7 degrees C and characterized by an enthalpy change of 81 kcal/mol and heat capacity difference between unfolded and folded states of 0.9 kcal/K*mol. By measuring the temperature dependence of the frequency dispersion of the heat capacity in the frequency range of 0.1-1 Hz it has been possible to examine the time regime of the enthalpic events associated with the transition. The multifrequency calorimetry results indicate that approximately 85% of the excess heat capacity associated with the folding/unfolding transition relaxes with a single relaxation time of 326 +/- 68 ms at the midpoint of the transition region. This is the first time that the time regime in which heat is absorbed and released during protein folding/unfolding has been measured.     

Comparing the thermodynamic stabilities of a related thermophilic and mesophilic enzyme

Several models have been proposed to explain the high temperatures required to denature enzymes from thermophilic organisms; some involve greater maximum thermodynamic stability for the thermophile, and others do not. To test these models, we reversibly melted two analogous protein domains in a two-state manner. E2cd is the isolated catalytic domain of cellulase E2 from the thermophile Thermomonospora fusca. CenAP30 is the analogous domain of the cellulase CenA from the mesophile Cellulomonas fimi. When reversibly denatured in a common buffer, the thermophilic enzyme E2cd had a temperature of melting (Tm) of 72.2 degrees C, a van't Hoff enthalpy of unfolding (DeltaHVH) of 190 kcal/mol, and an entropy of unfolding (DeltaSu) of 0.55 kcal/(mol*K); the mesophilic enzyme CenAP30 had a Tm of 56.4 degrees C, a DeltaHVH of 107 kcal/mol, and a DeltaSu of 0. 32 kcal/(mol*K). The higher DeltaHVH and DeltaSu values for E2cd suggest that its free energy of unfolding (DeltaGu) has a steeper dependence on temperature at the Tm than CenAP30. This result supports models that predict a greater maximum thermodynamic stability for thermophilic enzymes than for their mesophilic counterparts. This was further explored by urea denaturation. Under reducing conditions at 30 degrees C, E2cd had a concentration of melting (Cm) of 5.2 M and a DeltaGu of 11.2 kcal/mol; CenAP30 had a Cm of 2.6 M and a DeltaGu of 4.3 kcal/mol. Under nonreducing conditions, the Cm and DeltaGu of CenAP30 were increased to 4.5 M and 10.8 kcal/mol at 30 degrees C; the Cm for E2cd was increased to at least 7.4 M at 32 degrees C. We were unable to determine a DeltaGu value for E2cd under nonreducing conditions due to problems with reversibility. These data suggest that E2cd attains its greater thermal stability (DeltaTm = 15.8 degrees C) through a greater thermodynamic stability (DeltaDeltaGu = 6.9 kcal/mol) compared to its mesophilic analogue CenAP30.     

Equilibrium unfolding mechanism of chicken muscle triose phosphate isomerase

Triose phosphate isomerase (TIM) was prepared and purified from chicken breast muscle. The equilibrium unfolding of TIM by urea was investigated by following the changes of intrinsic fluorescence and circular dichroism spectroscopy, and the equilibrium thermal unfolding by differential scanning calorimetry (DSC). Results show that the unfolding of TIM in urea is highly cooperative and no folding intermediate was detected in the experimental conditions used. The thermodynamic parameters of TIM during its urea induced unfolding were calculated as DeltaG degrees =3.54 kcal.mol(-1), and m(G) = 0.67 kcal.mol(-1)M(-1), which just reflect the unfolding of dissociated folded monomer to fully unfolded monomer transition, while the dissociation energy of folded dimer to folded monomer is probe silence. DSC results indicate that TIM unfolding follows an irreversible two-state step with a slow aggregation process. The cooperative unfolding ratio, DeltaH(cal)/DeltaH(vH), was measured close to 2, indicating that the two subunits of chicken muscle TIM unfold independently. The van't Hoff enthalpy, DeltaH(vH), was estimated as about 200 kcal.mol(-1). These results support the unfolding mechanism with a folded monomer formation before its tertiary structure and secondary structure unfolding.     

Molecular mechanism for stabilizing a short helical peptide studied by generalized-ensemble simulations with explicit solvent

We study the folding mechanism of an analog of the C-peptide of ribonuclease A in explicit water by a replica-exchange multicanonical molecular dynamics simulation based on all-atom models. The multicanonical weight factor was determined by the combined use of the multicanonical replica-exchange method and the replica-exchange multicanonical algorithm. Using statistically reliable data thus obtained, we have examined the free-energy landscape of the peptide system. The global-minimum free-energy state in the landscape at room temperature has an alpha-helix structure with a distortion near the N-terminus. The state also has a salt bridge between Glu(-)-2 and Arg(+)-10 and an aromatic-aromatic interaction between Phe-8 and His(+)-12, both of which have been observed in x-ray and other experimental measurements. Principal component analysis clearly shows the different roles of these side-chain interactions in the peptide folding. The side-chain interaction between Phe-8 and His(+)-12 greatly enhances the stability of helical structure toward the C-terminal end, whereas the salt bridge between Glu(-)-2 and Arg(+)-10 mainly works as a restraint to prevent the alpha-helix structure from extending to the N-terminus. The free-energy landscape of C-peptide reveals a funnel-like shape where all of these interactions consistently exist only in the global-minimum state. This is the major reason why the native structure of the short helical peptide shows significant stability at low temperatures.     

Structural dynamics of neuropeptide hPYY

We have employed a combination of experiment and simulation to characterize the ensemble of structures sampled by human Peptide YY (hPYY), an important member of the neuropeptide Y family. Experimental structural characterization carried out with far UV circular dichroism spectroscopy and Fourier Transform-Infrared measurements confirmed that the major feature of the secondary structure of hPYY is the α-helix, encompassing about half the peptide residues, with smaller contributions from turn and β-sheet like structures. The peptide undergoes thermal denaturation characterized by a melting temperature of 48°C with an enthalpy change of -24.5 kcal/mol and entropy change of -76.2 cal/(mol K). In our computational studies, based on a 4-μsec MD trajectory generated with the AMBER03 potential, we found excellent agreement of the predicted features with experimental data, including a stable C-terminal helix, a central turn and conservation of about 80% of measured long-range NOE contacts. The main structural fluctuations involved partial helix unwinding and large-scale motions of the N-terminal. Our joint experimental/computational approach leads to several insights into the biological function of PYY. We conclude that the C-terminal helix is crucial for the structural integrity of PYY. The structures and motions found in the simulations suggest microscopic explanations for observed changes in biological activity of the peptide upon mutation and truncation. We also performed microsecond-length MD and replica-exchange simulations of hPYY with the OPLS-AA force field, for which computed structures did not agree well with experimental data, predicting significant loss of helicity and NOE contacts.     

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.     

Thermodynamics and folding pathway of tetraloop receptor-mediated RNA helical packing

Little is known about the thermodynamic forces that drive the folding pathways of higher-order RNA structure. In this study, we employ calorimetric [isothermal titration calorimetry (ITC) and differential scanning calorimetry (DSC)] and spectroscopic (NMR and UV) methods to characterize the thermodynamics of the GAAA tetraloop-receptor interaction, utilizing a previously described bivalent construct. ITC studies indicate that the bivalent interaction is enthalpy driven and highly stable, with a binding constant (K(obs)) of 5.5x10(6) M(-1) and enthalpy (DeltaH(obs)(o)) of -33.8 kcal/mol at 45 degrees C in 20 mM KCl and 2 mM MgCl(2). Thus, we derive the DeltaH(obs)(o) for a single tetraloop-receptor interaction to be -16.9 kcal/mol at these conditions. UV absorbance data indicate that an increase in base stacking quality contributes to the enthalpy of complex formation. These highly favorable thermodynamics are consistent with the known critical role for the tetraloop-receptor motif in the folding of large RNAs. Additionally, a significant heat capacity change (DeltaC(p,obs)(o)) of -0.24 kcal mol(-1) K(-1) was determined by ITC. DSC and UV-monitored thermal denaturation experiments indicate that the bivalent tetraloop-receptor construct follows a minimally five-state unfolding pathway and suggest the observed DeltaC(p,obs)(o) for the interaction results from a temperature-dependent unbound receptor RNA structure.     

Thermal unfolding of the archaeal DNA and RNA binding protein Ssh10

The reversible thermal unfolding of the archaeal histone-like protein Ssh10b from the extremophile Sulfolobus shibatae was studied using differential scanning calorimetry and circular dichroism spectroscopy. Analytical ultracentrifugation and gel filtration showed that Ssh10b is a stable dimer in the pH range 2.5–7.0. Thermal denaturation data fit into a two-state unfolding model, suggesting that the Ssh10 dimer unfolds as a single cooperative unit with a maximal melting temperature of 99.9 °C and an enthalpy change of 134 kcal/mol at pH 7.0. The heat capacity change upon unfolding determined from linear fits of the temperature dependence of ΔH^cal is 2.55 kcal/(mol K). The low specific heat capacity change of 13 cal/(mol K residue) leads to a considerable flattening of the protein stability curve (ΔG (T)) and results in a maximal ΔG of only 9.5 kcal/mol at 320 K and a ΔG of only 6.0 kcal/mol at the optimal growth temperature of Sulfolobus.     

A calorimetric study of the folding-unfolding of an alpha-helix with covalently closed N and C-terminal loops.

The thermal melting of a dicyclic 29-residue peptide, having helix-stabilizing side-chain to side-chain covalent links at each terminal, has been studied by circular dichroism spectropolarimetry (CD) and differential scanning calorimetry (DSC). The CD spectra for this dicyclic peptide indicate that it is monomeric, almost fully alpha-helical at -10 degrees C, and undergoes a reversible transition from the folded to the disordered state with increasing temperature. The temperature dependencies of the ellipticity at 222 nm and the excess heat capacity measured calorimetrically are well fit by a two-state model, which indicates a cooperative melting transition that is complete within the temperature ranges of these experiments (from -10 degrees C to 100 degrees C). This allows a complete analysis of the thermodynamics of helix formation. The helix unfolding is found to proceed with a small positive heat-capacity increment, consistent with the solvation of some non-polar groups upon helix unfolding. It follows that the hydrogen bonds are not the only factors responsible for the formation of the alpha-helix, and that hydrophobic interactions are also playing a role in its stabilization. At 30 degrees C, the calorimetric enthalpy and entropy values are estimated to be 650(+/-50) cal mol(-1)and 2.0(+/-0.2) cal K(-1)mole(-1), respectively, per residue of this peptide. Comparison with the thermodynamic characteristics obtained for the unfolding of double-stranded alpha-helical coiled-coils shows that at that temperature the enthalpic contribution of non-polar groups to the stabilization of the alpha-helix is insignificant and the estimated transition enthalpy can be assigned to the hydrogen bonds. With increasing temperature, the increasing magnitude of the negative enthalpy of hydration of the exposed polar groups should decrease the helix-stabilizing enthalpy of the backbone hydrogen bonds. However, the helix-stabilizing negative entropy of hydration of these groups should also increase in magnitude with increasing temperature, offsetting this effect.     

Helix-coil transition of alanine peptides in water: force field dependence on the folded and unfolded structures

The force fields used in classical modeling studies are semiempirical in nature and rely on their validation by comparison of simulations with experimental data. The all-atom replica-exchange molecular dynamics (REMD) methodology allows us to calculate the thermodynamics of folding/unfolding of peptides and small proteins, and provides a way of evaluating the reliability of force fields. We apply the REMD to obtain equilibrium folding/unfolding thermodynamics of a 21-residue peptide containing only alanine residues in explicit aqueous solution. The thermodynamics of this peptide is modeled with both the OPLS/AA/L and the A94/MOD force fields. We find that the helical content and the values for the helix propagation and nucleation parameters for this alanine peptide are consistent with measurements on similar peptides and with calculations using the modified AMBER force field (A94/MOD). The nature of conformations, both folded and unfolded, that contributes to the helix-coil transition profile, however, is quite different between these two force fields.     

Folding and stability of the three-stranded beta-sheet peptide Betanova: insights from molecular dynamics simulations

The dynamics of the three-stranded beta-sheet peptide Betanova has been studied at four different temperatures (280, 300, 350, and 450 K by molecular dynamics simulation techniques, in explicit water. Two 20-ns simulations at 280 K indicate that the peptide remains very flexible under "folding" conditions sampling a range of conformations that together satisfy the nuclear magnetic resonance (NMR)-derived experimental constraints. Two simulations at 300 K (above the experimental folding temperature) of 20 ns each show partial formation of "native"-like structure, which also satisfies most of the NOE constraints at 280 K. At higher temperature, the presence of compact states, in which a series of hydrophobic contacts remain present, are observed. This is consistent with experimental observations regarding the role of hydrophobic contacts in determining the peptide's stability and in initiating the formation of turns and loops. A set of different structures is shown to satisfy NMR-derived distance restraints and a possible mechanism for the folding of the peptide into the NMR-determined structure is proposed.     

Folding study of an Aib-rich peptide in DMSO by molecular dynamics simulations.

To evaluate the ability of molecular dynamics (MD) simulations using atomic force-fields to correctly predict stable folded conformations of a peptide in solution, we show results from MD simulations of the reversible folding of an octapeptide rich in alpha-aminoisobutyric acid (2-amino-2-methyl-propanoic acid, Aib) solvated in di-methyl-sulfoxide (DMSO). This solvent generally prevents the formation of secondary structure, whereas Aib-rich peptides show a high propensity to form secondary structural elements, in particular 3(10)- and alpha-helical structures. Aib is, moreover, achiral, so that Aib-rich peptides can form left- or right-handed helices depending on the overall composition of the peptide, the temperature, and the solvation conditions. This makes the system an interesting case to study the ensembles of peptide conformations as a function of temperature by MD simulation. Simulations involving the folding and unfolding of the peptide were performed starting from two initial structures, a right-handed alpha-helical structure and an extended structure, at three temperatures, 298 K, 340 K, and 380 K, and the results are compared with experimental nuclear magnetic resonance (NMR) data measured at 298 K and 340 K. The simulations generally reproduce the available experimental nuclear Overhauser effect (NOE) data, even when a wide range of conformations is sampled at each temperature. The importance of adequate statistical sampling in order to reliably interpret the experimental data is discussed.     

Probing the conformation of a human apolipoprotein C-1 by amino acid substitutions and trimethylamine-N-oxide.

To test, at the level of individual amino acids, the conformation of an exchangeable apolipoprotein in aqueous solution and in the presence of an osmolyte trimethylamine-N-oxide (TMAO), six synthetic peptide analogues of human apolipoprotein C-1 (apoC-1, 57 residues) containing point mutations in the predicted alpha-helical regions were analyzed by circular dichroism (CD). The CD spectra and the melting curves of the monomeric wild-type and plasma apoC-1 in neutral low-salt solutions superimpose, indicating 31 +/- 4% alpha-helical structure at 22 degrees C that melts reversibly with T(m,WT) = 50 +/- 2 degrees C and van't Hoff enthalpy deltaH(v,WT)(Tm) = 18 +/- 2 kcal/mol. G15A substitution leads to an increased alpha-helical content of 42 +/- 4% and an increased T(m,G15A) = 57 +/- 2 degrees C, which corresponds to stabilization by delta deltaG(app) = +0.4 +/- 1.5 kcal/mol. G15P mutant has approximately 20% alpha-helical content at 22 degrees C and unfolds with low cooperativity upon heating to 90 degrees C. R23P and T45P mutants are fully unfolded at 0-90 degrees C. In contrast, Q31P mutation leads to no destabilization or unfolding. Consequently, the R23 and T45 locations are essential for the stability of the cooperative alpha-helical unit in apoC-1 monomer, G15 is peripheral to it, and Q31 is located in a nonhelical linker region. Our results suggest that Pro mutagenesis coupled with CD provides a tool for assigning the secondary structure to protein groups, which should be useful for other self-associating proteins that are not amenable to NMR structural analysis in aqueous solution. TMAO induces a reversible cooperative coil-to-helix transition in apoC-1, with the maximal alpha-helical content reaching 74%. Comparison with the maximal alpha-helical content of 73% observed in lipid-bound apoC-1 suggests that the TMAO-stabilized secondary structure resembles the functional lipid-bound apolipoprotein conformation.     

Stabilisation of alpha-helices by site-directed mutagenesis reveals the importance of secondary structure in the transition state for acylphosphatase folding

The effects of stabilising mutations on the folding process of common-type acylphosphatase have been investigated. The mutations were designed to increase the helical propensity of the regions of the polypeptide chain corresponding to the two alpha-helices of the native protein. Various synthetic peptides incorporating the designed mutations were produced and their helical content estimated by circular dichroism. The most substantial increase in helical content is found for the peptide carrying five mutations in the second alpha-helix. Acylphosphatase variants containing the corresponding mutations display, to different extents, enhanced conformational stabilities as indicated by equilibrium urea denaturation experiments monitored by changes of intrinsic fluorescence. All the protein variants studied here refold with apparent two-state kinetics. Mutations in the first alpha-helix are responsible for a small increase in the refolding rate, accompanied by a marked decrease in the unfolding rate. On the other hand, multiple mutations in the second helix result in a considerable increase in the refolding rate without any significant effect on the unfolding rate. Addition of trifluoroethanol was found to accelerate the folding of the acylphosphatase variants, the extent of the acceleration being inversely proportional to the intrinsic rate of folding of the corresponding mutant. The trifluoroethanol-induced acceleration is far less marked for those variants whose alpha-helical structure is efficiently stabilised by amino acid replacements. This observation suggests that trifluoroethanol acts in a similar manner to the stabilising mutations in promoting native-like secondary structure. Analysis of the kinetic data indicates that the second helix is fully consolidated in the transition state for folding of acylphosphatase, whereas the first helix is only partially formed. These data suggest that the second helix is an important element in the folding process of the protein.     

A scanning calorimetric study of unfolding equilibria in homodimeric chicken gizzard tropomyosins.

Using both circular dichroism (CD) and differential scanning calorimetry (DSC), several laboratories find that the thermal unfolding transitions of alpha alpha and beta beta homodimeric coiled coils of rabbit tropomyosin are multistate and display an overall unfolding enthalpy of near 300 kcal (mol dimer)(-1). In contrast, an extant CD study of beta beta and gamma gamma species of chicken gizzard tropomyosin concludes that their unfolding transitions are simple two-state transitions, with much smaller overall enthalpies (98 kcal mol(-1) for beta beta and 162 kcal mol(-1) for gamma gamma). However, these smaller enthalpies have been questioned, because they imply a concentration dependence of the melting temperatures that is far larger than observed by CD. We report here DSC studies of the unfolding of both beta beta and gamma gamma chicken gizzard homodimers. The results show that these transitions are very similar to those in rabbit tropomyosins in that 1) the overall unfolding enthalpy is near 300 kcal mol(-1); 2) the overall delta C(rho) values are significantly positive; 3) the various transitions are multistate, requiring at least two and as many as four domains to fit the DSC data. DSC studies are also reported on these homodimeric species of chicken gizzard tropomyosin with a single interchain disulfide cross-link. These results are also generally similar to those for the correspondingly cross-linked rabbit tropomyosins.     

Thermodynamics and mechanism of alpha helix initiation in alanine and valine peptides

We used molecular dynamics simulations to study the folding/unfolding of one of turn of an alpha helix in Ac-(Ala)3-NHMe and Ac-(Val)3-NHMe. Using specialized sampling techniques, we computed free energy surfaces as functions of a conformational coordinate that corresponds to alpha helices at small values and to extended conformations at large values. Analysis of the peptide conformations populated during the simulations showed that alpha helices, reverse turns, and extended conformations correspond to minima on the free energy surfaces of both peptides. The free energy difference between alpha helix and extended conformations, determined from the equilibrium constants for helix unfolding, is approximately -1 kcal/mol for Ac-(Ala)3-NHMe and -5 kcal/mol for Ac-(Val)3-NHMe. The mechanism observed in our simulations, which includes reverse turns as important intermediates along the helix folding/unfolding pathway, is consistent with a mechanism proposed previously. Our results predict that both peptides (but especially the Ala peptide) have a much larger equilibrium constant for helix initiation than is predicted by the helix-coil transition theory with the host-guest parameters. We also predict a much greater difference in the equilibrium constants than the theory predicts. Insofar as helix initiation is concerned, our results suggest that the large difference between the helical propensities of Ala and Val cannot be explained by simple concepts such as side-chain rotamer restriction or unfavorable steric interactions. Rather, the origin of the difference appears to be quite complicated because it involves subtle differences in the solvation of the two peptides. The two peptides have similar turn-extended equilibria but very different helix-turn equilibria, and the difference in helical propensities reflects the fact that the helix-turn equilibrium strongly favors the turns in Ac-(Val)3-NHMe, while it favors the helices in Ac-(Ala)3-NHMe. We also computed thermodynamic decompositions of the free energy surfaces, and these revealed that the helix-turn equilibria are vastly different primarily because the changes in peptide-water interactions that accompany helix-to-turn conformational changes are qualitatively different for the two peptides.