Chemokines are small secreted signaling proteins that direct leukocyte migration toward areas of inflammation during an immune response. This task is dependent upon two functional roles associated with chemokines: 1) high affinity binding and activation of specific G-protein-coupled receptors (GPCRs) resulting in cell chemotaxis (
1), and 2) binding to glycosaminoglycans (GAGs) found in the extracellular matrix resulting in establishment of a signaling gradient toward the area of injury or infection (
1). Chemokines are classified into four subfamilies (CX3C, CXC, CC and C) based on the configuration of disulfide bonds within the N terminus (
2). Lymphotactin (Ltn/XCL1) is the defining member of the C-type subfamily, which contains only one of the two disulfide bonds found in all other chemokines. Ltn is unique in the way it bridges these two functional aspects of chemokine physiology, transforming between a monomeric canonical chemokine fold that can stimulate its cognate GPCR XCR1 (
3), and a dimeric β-sheet configuration capable of binding extracellular GAGs with high affinity (
4) (). Depending on solution conditions Ltn adopts the conserved chemokine fold, termed Ltn10 (10 °C, 200 mM NaCl), or the unique dimer configuration, termed Ltn40 (40 °C, 0mM NaCl) (
5). At equilibrium under near-physiological conditions (37 °C, 150 mM NaCl) both conformations are equally populated and each is required for proper biological function (
4). The unique structural interconversion of Ltn has led to further classification as a ‘metamorphic protein’. These newly discovered proteins are capable of adopting different folded conformations for the same primary sequence. However unlike other members of this emerging class, Ltn is the only naturally occurring protein known to populate two distinct folds in a native-state equilibrium that involves a complete restructuring of core residues (
6).
While many of the distinctive structural and functional properties of each Ltn conformation have been characterized (
3–
5,
7,
8), the exact mechanism of structural interconversion has remained elusive. Because native-state interconversion disrupts every hydrogen bond defining the secondary structure and replaces it with a different pairing (
4) we speculated that the transition proceeds through a chiefly disordered state. As an alternative to global unfolding, a structured intermediate could lower the free energy barrier separating Ltn10 and Ltn40 on the conformational energy landscape. Recently, a specific transition-state intermediate resembling a Ltn40 monomer with some Ltn10 tertiary contacts was proposed on the basis of coarse-grained molecular dynamic simulations of the Ltn10-Ltn40 equilibrium (
9). As a first step toward a mechanistic description of this process, we studied the temperature dependence of Ltn interconversion and unfolding by stopped-flow fluorescence to determine if the Ltn structural transformation is thermodynamically similar to protein unfolding, or employs a low energy transition-state structure.
Ltn contains one tryptophan (W55) that can report on the conformational state of the protein (
8). For instance in the Ltn10 conformation W55 is buried within a hydrophobic pocket formed by residues A49, A53, and V59, and displays maximum fluorescence intensity at 330nm. In contrast W55 in Ltn40 is solvent exposed resulting in reduced, red-shifted fluorescence emission relative to the Ltn10 profile. It has been shown that conformational equilibrium can be shifted to favor formation of Ltn10 or Ltn40 by changing solution conditions through addition of NaCl or low molecular weight heparin, respectively (
8). By measuring the change in fluorescence intensity in Ltn samples subjected to rapid mixing with either salt or heparin solutions we were able to monitor conversion to each Ltn state ().
The resulting time-dependent fluorescence intensities were fit to single exponential curves to yield kinetic rates. Analysis revealed a Ltn40 → Ltn10 conversion rate (k
for) of 0.4 s
−1 at 308 K, and a Ltn10 → Ltn40 conversion rate (k
rev) of 1.4 s
−1 at 308 K, in excellent agreement with previously determined values of k
for = 0.4 s
−1 and k
rev = 1.2 s
−1 derived from equilibrium measurements of interconversion by 2D exchange NMR (
8). Based on the large difference in protein concentration employed in fluorescence and NMR experiments (30µM vs. 1mM, respectively) and the absence of heparin in 2D exchange measurements, the results indicated that conformational change is likely a rate limiting process with rate constants independent of total protein or heparin concentrations. This analysis is further supported by estimates of the Ltn dimer K
d in the low micromolar range (
5) that would predict a k
off rate > 10s
−1 (assuming diffusion controlled association ~10
7 M
−1s
−1), which is considerably faster than the measured conversion rate of 0.4s
−1.
The temperature dependence of both conversion rates was then studied to probe the kinetic barrier between the two Ltn conformers. A complete description of experimental conditions is presented in the
supplemental material. The rate constants of interconversion were related to the free energy difference between the corresponding state (Ltn10 or Ltn40) and the transition-state by the Eyring equation:
that can be rewritten as,
where
kB is Boltzmann’s constant,
h is Planck’s constant,
R is the gas constant, and
T is temperature. While we acknowledge the limited validity of Eyring analysis as it applies to protein thermodynamics as previously discussed by d’Avignon et al (
10), this approach has been used extensively in the context of protein interconversion (
11), as well as protein folding and unfolding (
12–
15). Plots of ln(k/T) vs. 1/T () for the forward and reverse reactions were separately fit to
equation 2 allowing estimation of the enthalpic and entropic contributions to the interconversion transition-state.
Between 308 K and 323 K, kfor exhibited a linear Eyring relationship with values ranging from 0.4 s−1 to 1.6 s−1, respectively (). Eyring analysis revealed that a large enthalpic barrier (17.7 kcal/mol) and a modest entropic barrier (−0.8 kcal/mol at 298 K) must be overcome during Ltn40 to Ltn10 interconversion. A list of all thermodynamic parameters determined is displayed in .
| Table 1Thermodynamic parameters for Ltn interconversion and unfolding determined by stopped-flow fluorescence. |
Similar analysis of the krev rates measured between 291 K and 308K also resulted in a linear Eyring dependence with rates between 0.3 s−1 and 1.4 s−1, respectively (). Eyring analysis again indicated large enthalpic (15.6 kcal/mol) and modest entropic (−2.1 kcal/mol at 298 K) barriers associated with the Ltn10 → Ltn40 interconversion transition state at 298K. The temperature dependence of interconversion rates therefore indicates that similar enthalpic and entropic contributions are required for both the Ltn10→Ltn40 and Ltn40→Ltn10 transformations. The large enthalpic barrier (ΔH‡) clearly reflected disruption of non-covalent bonding interactions. However it was unclear if these values represented partial loss of secondary and tertiary structures or a global unfolding event. Entropy changes of the transition-state (ΔS‡) arise from two effects that are typically counterbalanced in globular proteins: conformational freedom of the polypeptide chain, and ordering of solvent molecules around nonpolar residues (the ‘hydrophobic effect’). Because Ltn interconversion is a global structural change with a large ΔH‡, it is unlikely that the individual entropic factors would be unchanged in the transition state. Thus, we reasoned that the small ΔS‡ might correspond to large opposing changes in both solvent exposure and chain entropy, which appears to be most consistent with a global unfolding event. Therefore we investigated the temperature dependence of Ltn10 and Ltn40 unfolding to determine if the resulting transition states were thermodynamically similar to the interconversion process.
Unfolding rates were determined by stopped-flow fluorescence for both Ltn forms using urea as the chemical denaturant. Final kinetic rates for each unfolding transition (k
u10 or k
u40) used in the analysis were extrapolated to 0 M denaturant based on a two-state unfolding model (
Supplemental material, Figure S1). Selection of each Ltn conformer present at the start of denaturant mixing experiments was accomplished by adjustment of temperature and solution conditions. Ltn10 unfolding was measured by buffering solutions in 20 mM NaH
2PO
4, with 200 mM NaCl, at temperatures below 298K, while the Ltn40 unfolding was measured by buffering solutions in 20mM NaH
2PO
4 with no added salt at temperatures above 313 K. The Ltn10 unfolding rates were measured by single exponential fitting of decay curves. Eyring analysis revealed a linear dependence with extrapolated unfolding rates of k
u10 = 3.4 to 4.7 sec
-−1 between 283 K and 291 K, yielding an enthalpic barrier of 23.2 kcal/mol and an entropic contribution of 9.1 kcal/mol calculated at 298 K (
Equation 2 and ). Analysis of Ltn40 unfolding rates were conducted by single exponential fitting of unfolding curves collected between 313 K and 323 K. Extrapolated unfolding rates of k
u40 = 0.5 to 1.3 sec
-−1 were determined between 313 K and 323 K, resulting in an enthalpic barrier of Ltn40
unfold of 17.9 kcal/mol and entropic contribution of −0.7 kcal/mol calculated at 298 K (
Equation2 and ).
Comparison of Ltn interconversion parameters () reveals very similar barrier heights that correspond to a ΔG
Ltn40-Ltn10 of 0.8 kcal/mol separating the two conformational states. This small free energy difference is consistent with NMR measurements of the Ltn conformational equilibrium indicating that K
eq ~ 1 at 298 K (
8). Both conversion events reveal that a large enthalpic barrier must be crossed, consistent with the disruption of hydrogen bonding networks that must accompany Ltn structural rearrangement (
5). Furthermore the slightly negative entropy change likely describes a reduction in water entropy from solvent exposure of hydrophobic core residues that exceeds the increase in chain entropy associated with the loss of tertiary structure, based on the virtually identical thermodynamic signature that characterize both Ltn40→Ltn10 conversion and Ltn40 unfolding. Thermodynamic parameters derived from Ltn10 unfolding also a show a large enthalpic signature, but reveal a favorable entropic contributions leading to a free energy barrier height that is less than structural interconversion. We attribute this disparity in barrier height to the most obvious difference in reaction pathways for the two processes: Ltn10 ↔ Ltn40 interconversion requires formation and disassembly of both tertiary and quaternary contacts, while Ltn10 unfolding reports only on the disruption of tertiary structure. Hence the difference in barrier heights would be expected given that additional free energy is required to bring two monomer subunits together during the Ltn10 → Ltn40 interconversion reaction.
Given the similar thermodynamic profiles of the transition states investigated in this study, we conclude that Ltn interconversion likely requires a global unfolding transition and does not pass through a structured low energy intermediate. This scenario is consistent with a previous H/D exchange study conducted by 2-D
1H-
15N- HSQC analysis (
8), that showed regardless of temperature or salt concentration Ltn protein samples lost all detectable amide signals within the experimental dead time (5 min.), reflecting an absence of stable hydrogen bonding networks within the timescale of interconversion (k
ex ~1.5 sec
−1). Interestingly, a recent report on the structural rearrangement between two β-strand registers of the ARNT PAS-B mutant Y456T (
11) also concluded that the protein interconverts by unfolding to a chiefly disordered state, and refolds into either conformation. An analogous process is likely at work in Ltn interconversion, however a complete mechanistic understanding is still being elucidated. Given that chemokine-mediated recruitment of cells
in vivo depends on distinct functions (GAG binding and GPCR activation) (
16) that are supplied independently by each of the folded Ltn native state structures, our results suggest that the ability to access the unfolded state is an essential feature of its biological activity. Considering the close structural similarity between the two ARNT PAS-B mutant forms (
11) and the disparate structural states of Ltn, these results also suggest that perhaps global unfolding is a common mechanism when structural rearrangement requires disruption of a large number of hydrogen bonds.