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.     

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.     

An engineered leucine zipper a position mutant with an unusual three-state unfolding pathway

The leucine zipper is a dimeric coiled-coil structural motif consisting of four to six heptad repeats, designated (abcdefg)(n). In the GCN4 leucine zipper, a position 16 in the third heptad is occupied by an Asn residue whereas the other a positions are Val residues. Recently, we have constructed variants of the GCN4 leucine zipper in which the a position Val residues were replaced by Ile. The folding and unfolding of the wild-type GCN4 leucine zipper and the Val to Ile variant both adhere to a simple two-state mechanism. In this study, another variant of the GCN4 leucine zipper was constructed by moving the single Asn residue from a position 16 to a position 9. This switch causes the thermal unfolding of the GCN4 leucine zipper to become three state. The unfolding pathway of this variant was determined by thermal denaturation, limited proteinase K digestion, and sedimentation equilibrium analysis. Our data are consistent with a model in which the variant first unfolds from its N terminus and changes the oligomerization specificity from a native dimer to a partially unfolded intermediate containing a mixture of dimers and trimers and then completely unfolds to unstructured monomers.     

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.     

T-jump infrared study of the folding mechanism of coiled-coil GCN4-p1

Partially folded intermediates have been frequently observed in equilibrium and kinetic protein folding studies. However, folding intermediates that exist at the native side of the rate-limiting step are rather difficult to study because they often evade detection by conventional folding kinetic methods. Here, we demonstrated that a laser-induced temperature-jump method can potentially be used to identify the existence of such post-transition or hidden intermediates. Specifically, we studied two cross-linked variants of GCN4-p1 coiled-coil. The GCN4 leucine zipper has been studied extensively and most of these studies have regarded it as a two-state folder. Our static circular dichroism and infrared data also indicate that the thermal unfolding of these two monomeric coiled-coils can be adequately described by an apparent two-state model. However, their temperature-jump-induced relaxation kinetics exhibit non-monoexponential behavior, dependent upon sequence and temperature. Taken together, our results support a folding mechanism wherein at least one folding intermediate populates behind the main rate-limiting step.     

Folding of a designed simple ankyrin repeat protein

Ankyrin repeats (AR) are 33-residue motifs containing a beta-turn, followed by two alpha-helices connected by a loop. AR occur in tandem arrangements and stack side-by-side to form elongated domains involved in very different cellular tasks. Recently, consensus libraries of AR repeats were constructed. Protein E1_5 represents a member of the shortest library, and consists of only a single consensus repeat flanked by designed N- and C-terminal capping repeats. Here we present a biophysical characterization of this AR domain. The protein is compactly folded, as judged from the heat capacity of the native state and from the specific unfolding enthalpy and entropy. From spectroscopic data, thermal and urea-induced unfolding can be modeled by a two-state transition. However, scanning calorimetry experiments reveal a deviation from the two-state behavior at elevated temperatures. Folding and unfolding at 5 degrees C both follow monoexponential kinetics with k(folding) = 28 sec(-1) and k(unfolding) = 0.9 sec(-1). Kinetic and equilibrium unfolding parameters at 5 degrees C agree very well. We conclude that E1_5 folds in a simple two-state manner at low temperatures while equilibrium intermediates become populated at higher temperatures. A chevron-plot analysis indicates that the protein traverses a very compact transition state along the folding/unfolding pathway. This work demonstrates that a designed minimal ankyrin repeat protein has the thermodynamic and kinetic properties of a compactly folded protein, and explains the favorable properties of the consensus framework.     

Non-linear effects of temperature and urea on the thermodynamics and kinetics of folding and unfolding of hisactophilin

Extensive measurements and analysis of thermodynamic stability and kinetics of urea-induced unfolding and folding of hisactophilin are reported for 5-50 degrees C, at pH 6.7. Under these conditions hisactophilin has moderate thermodynamic stability, and equilibrium and kinetic data are well fit by a two-state transition between the native and the denatured states. Equilibrium and kinetic m values decrease with increasing temperature, and decrease with increasing denaturant concentration. The betaF values at different temperatures and urea concentrations are quite constant, however, at about 0.7. This suggests that the transition state for hisactophilin unfolding is native-like and changes little with changing solution conditions, consistent with a narrow free energy profile for the transition state. The activation enthalpy and entropy of unfolding are unusually low for hisactophilin, as is also the case for the corresponding equilibrium parameters. Conventional Arrhenius and Eyring plots for both folding and unfolding are markedly non-linear, but these plots become linear for constant DeltaG/T contours. The Gibbs free energy changes for structural changes in hisactophilin have a non-linear denaturant dependence that is comparable to non-linearities observed for many other proteins. These non-linearities can be fit for many proteins using a variation of the Tanford model, incorporating empirical quadratic denaturant dependencies for Gibbs free energies of transfer of amino acid constituents from water to urea, and changes in fractional solvent accessible surface area of protein constituents based on the known protein structures. Noteworthy exceptions that are not well fit include amyloidogenic proteins and large proteins, which may form intermediates. The model is easily implemented and should be widely applicable to analysis of urea-induced structural transitions in proteins.     

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.     

Stabilization of coiled-coil peptide domains by introduction of trifluoroleucine

Substitution of leucine residues by 5,5,5-trifluoroleucine at the d-positions of the leucine zipper peptide GCN4-p1d increases the thermal stability of the coiled-coil structure. The midpoint thermal unfolding temperature of the fluorinated peptide is elevated by 13 degrees C at 30 microM peptide concentration. The modified peptide is more resistant to chaotropic denaturants, and the free energy of folding of the fluorinated peptide is 0.5-1.2 kcal/mol larger than that of the hydrogenated form. A similarly fluorinated form of the DNA-binding peptide GCN4-bZip binds to target DNA sequences with affinity and specificity identical to those of the hydrogenated form, while demonstrating enhanced thermal stability. Molecular dynamics simulation on the fluorinated GCN4-p1d peptide using the Surface Generalized Born implicit solvation model revealed that the coiled-coil binding energy is 55% more favorable upon fluorination. These results suggest that fluorination of hydrophobic substructures in peptides and proteins may provide new means of increasing protein stability, enhancing protein assembly, and strengthening receptor-ligand interactions.     

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.     

Stabilization of native protein fold by intein-mediated covalent cyclization

A mutant version of the N-terminal domain of Escherichia coli DnaB helicase was used as a model system to assess the stabilization against unfolding gained by covalent cyclization. Cyclization was achieved in vivo by formation of an amide bond between the N and C termini with the help of a split mini-intein. Linear and circular proteins were constructed to be identical in amino acid sequence. Mutagenesis of Phe102 to Glu rendered the protein monomeric even at high concentration. A difference in free energy of unfolding, DeltaDeltaG, between circular and linear protein of 2.3(+/-0.5) kcal mol(-1) was measured at 10 degrees C by circular dichroism. A theoretical estimate of the difference in conformational entropy of linear and circular random chains in a three-dimensional cubic lattice model predicted DeltaDeltaG=2.3 kcal mol(-1), suggesting that stabilization by protein cyclization is driven by the reduced conformational entropy of the unfolded state. Amide-proton exchange rates measured by NMR spectroscopy and mass spectrometry showed a uniform, approximately tenfold decrease of the exchange rates of the most slowly exchanging amide protons, demonstrating that cyclization globally decreases the unfolding rate of the protein. The amide proton exchange was found to follow EX1 kinetics at near-neutral pH, in agreement with an unusually slow refolding rate of less than 4 min(-1) measured by stopped-flow circular dichroism. The linear and circular proteins differed more in their unfolding than in their folding rates. Global unfolding of the N-terminal domain of E.coli DnaB is thus promoted strongly by spatial separation of the N and C termini, whereas their proximity is much less important for folding.     

The magnitude of the backbone conformational entropy change in protein folding

The magnitude of the conformational entropy change experienced by the peptide backbone upon protein folding was investigated experimentally and by computational analysis. Experimentally, two different pairs of mutants of a 33 amino acid peptide corresponding to the leucine zipper region of GCN4 were used for high-sensitivity microcalorimetric analysis. Each pair of mutants differed only by having alanine or glycine at a specific solvent-exposed position under conditions in which the differences in stability could be attributed to differences in the conformational entropy of the unfolded state. The mutants studied were characterized by different stabilities but had identical heat capacity changes of unfolding (ΔC_p), identical solvent-related entropies of unfolding (ΔS_solv), and identical enthalpies of unfolding (ΔH) at equivalent temperatures. Accordingly, the differences in stability between the different mutants could be attributed to differences in conformational entropy. The computational studies were aimed at generating the energy profile of backbone conformations as a function of the main chain dihedral angles ϕ and ϱ. The energy profiles permit a direct calculation of the probability distribution of different conformers and therefore of the conformational entropy of the backbone. The experimental results presented in this paper indicate that the presence of the methyl group in alanine reduces the conformational entropy of the peptide backbone by 2.46 ± 0.2 cal/K · mol with respect to that of glycine, consistent with a 3.4-fold reduction in the number of allowed conformations in the alanine-containing peptides. Similar results were obtained from the energy profiles. The computational analysis also indicates that the addition of further carbon atoms to the side chain had only a small effect as long as the side chains were unbranched at position β. A further reduction with respect to Ala of only 0.61 and 0.81 cal/K · mol in the backbone entropy was obtained for leucine and lysine, respectively. β-branching (Val) produces the largest decrease in conformational entropy (1.92 cal/K · mol less than Ala). Finally, the backbone entropy change associated with the unfolding of an α-helix is 6.51 cal/K · mol for glycine. These and previous results have allowed a complete estimation of the conformational entropy changes associated with protein folding. © 1996 Wiley-Liss, Inc.     

Effect of ethanol on folding of hen egg-white lysozyme under acidic condition

The equilibrium and kinetics of folding of hen egg-white lysozyme were studied by means of CD spectroscopy in the presence of varying concentrations of ethanol under acidic condition. The equilibrium transition curves of guanidine hydrochloride-induced unfolding in 13 and 26% (v/v) ethanol have shown that the unfolding significantly deviates from a two-state mechanism. The kinetics of denaturant-induced refolding and unfolding of hen egg-white lysozyme were investigated by stopped-flow CD at three ethanol concentrations: 0, 13, and 26% (v/v). Immediately after dilution of the denaturant, the refolding curves showed a biphasic time course in the far-UV region, with a burst phase with a significant secondary structure and a slower observable phase. However, when monitored by the near-UV CD, the burst phase was not observed and all refolding kinetics were monophasic. To clarify the effect of nonnative secondary structure induced by the addition of ethanol on the folding/unfolding kinetics, the kinetic m values were estimated from the chevron plots obtained for the three ethanol concentrations. The data indicated that the folding/unfolding kinetics of hen lysozyme in the presence of varying concentrations of ethanol under acidic condition is explained by a model with both on-pathway and off-pathway intermediates of protein folding.