CD and (13)C(alpha)-NMR studies of folding equilibria in a two-stranded coiled coil formed by residues 190-254 of alpha-tropomyosin.

Synthesis and CD and (13)C(alpha)-NMR studies in a near-neutral saline buffer are reported for a 65-residue peptide ((190)Tm(254)) comprising residues 190-254 of the alpha-tropomyosin chain. CD on a version disulfide cross-linked via the N-terminal cysteine side chains indicates that this dimer is highly helical and melts near 48 degrees C. The CD is independent of peptide concentration, showing that association of (190)Tm(254) stops at the two-strand level. Similar studies on the reduced version show much lower helix content at low temperature, melting points below room temperature, and the expected concentration dependence. The observed melting temperature of the reduced peptide is far below (by 27 degrees C) that expected from an extant analysis of calorimetry data on parent tropomyosin that designates (190)Tm(254) as an independently melting "cooperative block." This disagreement and the pronounced nonadditivity seen when data for (190)Tm(254) are combined with extant data for other subsequences argue decisively against the concept of specific independently melting blocks within the tropomyosin chain. The data for (190)Tm(254) also serve to test recent ideas on the sequence determinants of structure and stability in coiled coils. Analysis shows that some ideas, such as the stabilizing effect of leucine in the d heptad position, find support, but others--such as the destabilizing effect of alanine in d, the dimer-disfavoring effect of beta-branching in d and its dimer-favoring effect in a, and the dimer-directing effect of asparagine in a--are more questionable in tropomyosin than in the leucine zipper coiled coils. (13)C(alpha)-NMR data at two labeled sites, L228(d) and V246(a), of (190)Tm(254) display well-separated resonances for folded and unfolded forms at each site, indicating that the transition is slow on the NMR time scale and thus demonstrating the possibility of obtaining thermodynamic and kinetic information on the transition at the residue level.

Temperature dependence of the folding and unfolding kinetics of the GCN4 leucine zipper via 13C(alpha)-NMR.

Studies by one-dimensional NMR are reported on the interconversion of folded and unfolded forms of the GCN4 leucine zipper in neutral saline buffer. The peptide bears 99% 13C(alpha) labels at three sites: V9, L12, and G31. Time-domain 13C(alpha)-NMR spectra are interpreted by global Bayesian lineshape analysis to extract the rate constants for both unfolding and folding as functions of temperature in the range 47-71 degrees C. The data are well fit by the assumption that the same rate constants apply at each labeled site, confirming that only two conformational states need be considered. Results show that 1) both processes require a free energy of activation; 2) unfolding is kinetically enthalpy-opposed and entropy-driven, while folding is the opposite; and 3) the transition state dimer ensemble averages approximately 40% helical. The activation parameters for unfolding, derived from NMR data at the elevated temperatures where both conformations are populated, lead to estimates of the rate constant at low temperatures (5-15 degrees C) that agree with extant values determined by stopped-flow CD via dilution from denaturing media. However, the corresponding estimated values for the folding rate constant are larger by two to three orders of magnitude than those obtained by stopped flow. We propose that this apparent disagreement is caused by the necessity, in the stopped-flow experiment, for initiation of new helices as the highly denaturant-unfolded molecule adjusts to the newly created benign solvent conditions. This must reduce the success rate of collisions in producing the folded molecule. In the NMR determinations, however, the unfolded chains always have a small, but essential, helix content that makes such initiation unnecessary. Support for this hypothesis is adduced from recent extant experiments on the helix-coil transition in single-chain helical peptides and from demonstration that the folding rate constants for coiled coils, as obtained by stopped flow, are influenced by the nature of the denaturant used.

Ultracentrifuge and circular dichroism studies of folding equilibria in a retro GCN4-like leucine zipper.

Equilibrium ultracentrifuge and circular dichroism (CD) studies of a retropeptide of a GCN4-like leucine zipper in neutral saline buffer are reported as functions of temperature. Ultracentrifuge results indicate the presence of three oligomeric species: monomer, dimer, and tetramer, in quantifiable amounts, and the data provide values for the standard DeltaG, DeltaH, and DeltaS for interconversion. CD at 222 nm displays the strong concentration dependence characteristic of dissociative unfolding, but also shows a helicity far below that of the parent propeptide. Remarkably enough, the CD at 222 nm shows an extremum in the region between 0 and 20 degrees C. At higher T, the usual cooperative unfolding is observed. Comparable data are presented for a mutant retropeptide, in which a single asparagine residue is restored to the characteristic heptad position it occupies in the propeptide. The mutant shows marked differences from its unmutated relative in both thermodynamic properties and CD, although the oligomeric ensemble also comprises monomers, dimers, and tetramers. The mutant is closer in helicity to the parent propeptide but is less stable. These findings do not support either of the extant views on retropeptides. The behavior seen is consistent neither with the view that retropeptides should have the same structure as propeptides nor with the view that they should have the same structure but opposite chirality. The simultaneous availability of oligomeric population data and CD allows the latter to be dissected into individual contributions from monomers, dimers, and tetramers. This dissection yields explanations for the observed extrema in curves of CD (222 nm) versus T and reveals that the dimer population in both retropeptides undergoes "cold denaturation."

Thermodynamics and kinetics of a folded-folded' transition at valine-9 of a GCN4-like leucine zipper.

Spin inversion transfer (SIT) NMR experiments are reported probing the thermodynamics and kinetics of interconversion of two folded forms of a GCN4-like leucine zipper near room temperature. The peptide is 13Calpha-labeled at position V9(a) and results are compared with prior findings for position L13(e). The SIT data are interpreted via a Bayesian analysis, yielding local values of T1a, T1b, kab, kba, and Keq as functions of temperature for the transition FaV9 right arrow over left arrow FbV9 between locally folded dimeric forms. Equilibrium constants, determined from relative spin counts at spin equilibrium, agree well with the ratios kab/kba from the dynamic SIT experiments. Thermodynamic and kinetic parameters are similar for V9(a) and L13(e), but not the same, confirming that the molecular conformational population is not two-state. The energetic parameters determined for both sites are examined, yielding conclusions that apply to both and are robust to uncertainties in the preexponential factor (kT/h) of the Eyring equation. These conclusions are 1) the activation free energy is substantial, requiring a sparsely populated transition state; 2) the transition state's enthalpy far exceeds that of either Fa or Fb; 3) the transition state's entropy far exceeds that of Fa, but is comparable to that of Fb; 4) "Arrhenius kinetics" characterize the temperature dependence of both kab and kba, indicating that the temperatures of slow interconversion are not below that of the glass transition. Any postulated free energy surface for these coiled coils must satisfy these constraints.

Site-specific thermodynamics and kinetics of a coiled-coil transition by spin inversion transfer NMR.

A 33-residue pseudo-wild-type GCN4 leucine zipper peptide is used to probe the equilibrium conformational population in proteins. 13Calpha-NMR shows that chain sites differ in structural content at a given temperature, and that two dimeric folded forms are evident at many sites. Spin inversion transfer experiments are reported bearing on the thermodynamics and kinetics of interconversion of the two dimeric folded forms (Fa <--> Fb) at the 13Calpha-labeled position L13. At each temperature, at conditions wherein the population of unfolded chains is quite small, inversion of the Fa spins via a tuned Gaussian pi-pulse is followed by a time interval (tau), interrogation, and recording of the free induction decay. Fifteen such inversions, with varying tau, provide the time course for recovery of equilibrium magnetization after inversion. Similar experiments follow inversion of the Fb spins. Re-equilibration is known to be modulated by four first-order rate constants: two (T1a(-1) and T1b(-1)) for spin-lattice relaxation intrinsic to the respective sites, and two (kab and kba) for the conformational change. All four follow from joint, Bayesian analysis of all the data at each temperature. The equilibrium constant at each temperature for this local transition, determined simply from the equilibrium relative magnetizations at Fa and Fb sites, agrees well with the kinetic ratio kab/kba. The standard Gibbs energies, enthalpy, and entropy follow. Activation parameters, both ways, are accessible from the rate constants and suggest a transition state with high Gibbs energy and enthalpy, but with entropy between those of Fa and Fb.

Thermal unfolding in a GCN4-like leucine zipper: 13C alpha NMR chemical shifts and local unfolding curves.

13C alpha chemical shifts and site-specific unfolding curves are reported for 12 sites on a 33-residue, GCN4-like leucine zipper peptide (GCN4-lzK), ranging over most of the chain and sampling most heptad positions. Data were derived from NMR spectra of nine synthetic, isosequential peptides bearing 99% 13C alpha at sites selected to avoid spectral overlap in each peptide. At each site, separate resonances appear for unfolded and folded forms, and most sites show resonances for two folded forms near room temperature. The observed chemical shifts suggest that 1) urea-unfolded GCN4-lzK chains are randomly coiled; 2) thermally unfolded chains include significant transient structure, except at the ends; 3) the coiled-coli structure in the folded chains is atypical near the C-terminus; 4) only those interior sites surrounded by canonical interchain salt bridges fail to show two folded forms. Local unfolding curves, obtained from integrated resonance intensities, show that 1) sites differ in structure content and in melting temperature, so the equilibrium population must comprise more than two molecular conformations; 2) there is significant end-fraying, even at the lowest temperatures, but thermal unfolding is not a progressive unwinding from the ends; 3) residues 9-16 are in the lowest melting region; 4) heptad position does not dictate stability; 5) significant unfolding occurs below room temperature, so the shallow, linear decline in backbone CD seen there has conformational significance. It seems that only a relatively complex array of conformational states could underlie these findings.

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.

Observation via one-dimensional 13Calpha NMR of local conformational substates in thermal unfolding equilibria of a synthetic analog of the GCN4 leucine zipper.

Synthesis of a 33-residue, capped leucine zipper analogous to that in GCN4 is reported. Histidine and arginine residues are mutated to lysine to reduce the unfolding temperature. CD and ultracentrifugation studies indicate that the molecule is a two-stranded coiled coil under benign conditions. Versions of the same peptide are made with 99% 13Calpha at selected sites. One-dimensional 13C NMR spectra are assigned by inspection and used to study thermal unfolding equilibria over the entire transition from 8 to 73 degrees C. Spectra at the temperature extremes establish the approximate chemical shifts for folded and unfolded forms at each labeled site. Resonances for each amino acid appear at both locations at intermediate T, indicating that folded and unfolded forms interconvert slowly (> >2 ms) on the NMR time scale. Moreover, near room temperature, the structured form's resonance is double at several, but not all, sites, indicating at least two slowly interconverting, structured, local conformational substates. Analysis of the dynamics is possible. For example, near room temperature at the Val-9, Ala-24, and Gly-31 positions, the equilibrium constant for interconversion of the two structured forms is near unity and the time scale is > or= 10-20 ms.

The use of spectral decomposition via the convex constraint algorithm in interpreting the CD-observed unfolding transitions of coiled coils.

Decomposition of CD spectra for the unfolding of both coiled-coil and single-helical molecules is carried out via the convex constraint algorithm (CCA) [A. Perczel, M. Hollósi, G. Tusnády, and G. D. Fasman (1991) Protein Engineering, Vol. 4, pp. 669-679]. Examined are (1) our thermal unfolding data for rabbit alpha alpha-tropomyosin and chicken gizzard gamma gamma-tropomyosin coiled coils, and for a 35-residue, tropomyosin-model peptide that forms single helices, not coiled coils; (2) extent pH-induced unfolding data for 50- and 400-residue poly-L-glutamic acid. Each set of spectra shows a sharp isodichroic point near 203 nm. We find here that the CCA is of sharply limited use for analyzing such data. The component spectra obtained for a given substance not only depend on the particular experimental spectra included and on the chosen number of component spectra, but all pass through the experimental isodichroic point. The latter is physically unlikely for more than three component spectra, and physically impossible for conformers, such as beta structures, having known isodichroic points elsewhere. Our conclusions are in contrast to those of an extant decomposition via CCA of thermal spectra for rabbit alpha alpha-tropomyosin [N. J. Greenfield and S. E. Hitchcock-DeGregori (1993) Protein Science, Vol. 2, pp. 1263-1273] that postulates the existence of five conformers, including beta structures, in the unfolding. Moreover, an extant diagnostic based on the theta 222/theta 208 ratio and allegedly distinguishing between spectra for coiled coil and for single alpha-helix is shown here to be unreliable.

The "cratic correction" and related fallacies.

The "cratic correction" for removing the contribution of translational motion from standard free energies and entropies of chemical reactions is examined. Although removal of translational effects is useful in interpreting free energy and entropy changes in molecular terms, this particular correction does not do so. No basis for it is found in thermodynamics or statistical mechanics. The appropriate correction for center-of-mass translational motion is derived from classical statistical mechanics. However, the importance of recognizing the proper molecular meaning of the remaining free energy is stressed. In any case, this latter correction, however legitimate, is shown to be futile, since it has long been recognized from classical statistical mechanics that all terms involving momenta (thermal wavelength terms)-be they classified as translational, rotational, vibrational, or whatever-necessarily cancel out in forming the standard free energy change of any chemical reaction from the constituent standard chemical potentials.

Some properties of acid-reassembled tropomyosin.

The preference for parallelism of the two chains in tropomyosin coiled coils is thought to result from interchain salt bridges. To examine this idea, studies are presented of tropomyosin molecules reassembled from chaotropic solvents in acid solution, where cross-links cannot exist. The acid-reassembled molecules are appreciably less disulfide cross-linkable in acid than native molecules, a result explainable if some antiparallel dimers indeed form at low pH. Physical studies (backbone- and tyrosine-region CD and intrinsic viscosity) indicate that refolding in acid yields a molecular population demonstrably different in tyrosine-region CD from native, but having comparable (but not identical) helix content, thermal stability, and dimensions. Moreover, the refolding in acid after either thermal or chaotropic-solvent denaturation yields the same final state, arguing that it is an equilibrium state. All these results are consistent with, but do not prove, that the acid-reassembled population includes an appreciable fraction (2/3) of antiparallel coiled-coil dimers.

Structural stability of short subsequences of the tropomyosin chain.

The native tropomyosin molecule is a parallel, registered, alpha-helical coiled coil made from two 284-residue chains. Long excised subsequences (> or = 95 residues) form the same structure with comparable thermal stability. Here, we investigate local stability using shorter subsequences (20-50 residues) that are chemically synthesized or excised from various regions along the protein chain. Thermal unfolding studies of such shorter peptides by CD in the same solvent medium used in extant studies of the parent protein indicate very low helix content, almost no coiled-coil formation, and high thermal lability of such secondary structure as does form. This behavior is in stark contrast to extant data on leucine-zipper peptides and short "designed" synthetic peptides, many of which have high alpha-helix content and form highly stable coiled coils. The existence of short coiled coils calls into question the older idea that short subsequences of a protein have little structure. The present study supports the older view, at least in its application to tropomyosin. The intrinsic local alpha-helical propensity and helix-helix interaction in this prototypical alpha-helical protein is sufficiently weak as to require not only dimerization, but macro-molecular amplification in order to attain its native conformation in common benign media near neutral pH.

Thermal unfolding equilibria in homodimeric chicken gizzard tropomyosin coiled coils.

CD studies are presented on thermal unfolding of coiled-coil homodimers of two genetic variant chains of chicken gizzard tropomyosin (CG-Tm). The experiments include the effects of cross-linking both isoforms and the dependence on protein concentration of unfolding in both reduced isoforms, variables not examined in extant work. The general shapes of the unfolding curves for singly cross-linked species depend on whether the cross-link is at C190 (its site on one isoform) or at C36 (its site on the other). These curves are compared with extant ones for various cross-linked species of rabbit tropomyosin. The comparison supports the view that the unfolding behavior of cross-linked species results from a complex interaction of strain at the cross-link, local variations in structural stability, and loop entropy. The observed concentration dependence of the transition temperature for the uncross-linked (reduced) species of CG-Tm is very small (2.9 degrees C) for one variant homodimer and unobservably small (< 2 degrees C) for the other in the 100-fold concentration range (approximately 0.01-1.0 mg/mL) accessible here. These experimental values of delta Tm are much smaller than are predicted from extant values of the van't Hoff transition enthalpies, calling the latter into question.

Does Flory-Huggins theory help in interpreting solute partitioning experiments?

The development of Flory-Huggins (FH) theory is reviewed, particularly with regard to the molecular significance of the interaction parameter that scales the contact interaction of solute and solvent. The chemical potential given by FH theory for an "idealute" solute is then compared with that provided by a more general, statistical thermodynamic approach. It is found that the FH contact term does not directly correspond to the solvation free energy. The significance of this result for the interpretation of free energies of transfer of a solute from one solvent to another is examined. It is shown that neither the earlier recommended standard free energy change for the process (using the infinitely dilute reference state, mole fraction units) nor the recently recommended FH-corrected standard free energy change provides the solvation energy desired. Instead, the standard free energy using the infinitely dilute reference state and molarity units, as long advocated by Ben-Naim, provides the desired solvation free energy. Correction of extant values, based on mole fraction units, is easily made. However, application of such results to problems of protein folding is not straightforward.

Further characterization of the kinetic folding intermediate of alpha alpha-tropomyosin and of its 142-281 subsequence.

The backbone CD spectrum from 250 to 212 nm for the kinetic folding intermediate of alpha alpha-tropomyosin (alpha alpha-Tm) and nonpolymerizable alpha alpha-Tm was obtained. The spectrum shows that the intermediate is indeed alpha-helical with about 70% of the equilibrium alpha-helix content. Subsequence 142Tm281 of the alpha-tropomyosin chain has five tyrosine residues (at positions 162, 214, 221, 261, 267). Stopped flow CD at the negative peak in the tyrosine spectral region (280 nm) shows that any tyrosine residues that contribute to the spectrum in the region have already reached their final state in the fast phase of folding (< 0.04 s).

Alpha-helix to random coil transitions: determination of peptide concentration from the CD at the isodichroic point.

A method is presented for determining the concentrations of peptides and proteins having isodichroic points near 203 nm. The existence of an isodichroic point for a given substance indicates a local two-state (alpha-helix, random coil) population. The mean residue ellipticity at the isodichroic point, [theta lambda i], is, of course, independent of helix content. For a wide variety of synthetic and natural peptides, including both single helices and coiled coils, it is shown that [theta lambda i] is also essentially independent of substance and of whether the transition is induced by temperature, ionic strength, pH, chain length changes, amino acid substitution, or solvent perturbation. Averaging [theta lambda i] values culled from various laboratories gives -151 +/- 16 (SD, 7 sources) deg.cm2.mmol-1. In our laboratory, nonpolymerizable rabbit alpha-tropomyosin and two alpha-tropomyosin subsequences yield -135 +/- 10 (SD, 190 values) deg.cm2.mmol-1. Thus, given [theta lambda i] for a peptide of known concentration, it is possible to estimate the concentration of any other peptide provided that it has an isodichroic point at which the ellipticity is accurately measurable. It is then possible to calculate [theta lambda] at any other wavelength for which theta is known. It is advisable to determine [theta lambda i] for the best known peptide in one's own laboratory, since it depends on absolute instrument and cell calibrations and an absolute concentration determination.

Kinetics of folding and unfolding of beta beta-tropomyosin.

The kinetics of folding from random coils to two-chain coiled coils of beta beta-tropomyosin was studied by stopped-flow CD (SFCD) in the backbone region (222 nm). Two species were studied: the reduced form and the doubly disulfide cross-linked form. The proteins were totally unfolded in 6M urea-saline buffer, then refolded by tenfold dilution into benign buffer. In the refolding medium, they spontaneously recover the two-chain coiled-coil structure. Reduced beta beta refolds in at least two stages: one or more fast phases (< 0.04 s), in which an intermediate with 71% of the equilibrium ellipticity forms, followed by a slower time-resolvable phase that completes the folding. The slow phase is first order, signifying that dimerization occurs in the fast phase. The time constant of the slow phase is 2 s at 20 degrees C and requires activation parameters of delta S not equal to = -7 +/- 0.3 cal/mol.K, delta H not equal to = 15 +/- 1 kcal/mol. These results are very similar to those previously found for the reduced genetic variant alpha alpha-tropomyosin. In contrast, refolding of doubly disulfide cross-linked beta beta is complete within the dead time (< 0.04 s), whereas the singly cross-linked alpha alpha species also displays a slow phase. The opposite process, unfolding reduced beta beta from the coiled-coil state, is complete within the dead time, as in the alpha alpha variant.