PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of nihpaAbout Author manuscriptsSubmit a manuscriptHHS Public Access; Author Manuscript; Accepted for publication in peer reviewed journal;
 
Biochemistry. Author manuscript; available in PMC 2010 September 1.
Published in final edited form as:
PMCID: PMC2746415
NIHMSID: NIHMS137242

Local and Global Effects of a Cavity Filling Mutation in a Metastable Serpin

Abstract

The serpins are an unusual class of protease inhibitors which fold to a metastable form and subsequently undergo a massive conformational change to a stable form when they inhibit their target proteases. The driving force for this conformational change has been extensively investigated by site directed mutagenesis, and it has been found that mutations which stabilize the metastable form frequently result in activity deficiency. Here we employ hydrogen/deuterium exchange to probe the effects of a cavity filling mutant of α1AT. The Gly117→Phe substitution fills a cavity between the F-helix and the face of β-sheet A, stabilizes the metastable form of α1AT by ~4 kcal/mole and results in a 60% reduction in inhibitory activity against elastase. Globally, the G117F substitution alters the unfolding mechanism by eliminating the molten globule intermediate that is seen in Wild Type unfolding. Remarkably, this is accomplished primarily by destabilizing the molten globule rather than stabilizing the metastable native state. Locally, conformational flexibility in the native state is reduced in specific regions: the top of the F-helix, β-strands 5A, 1C and 4C, and helix D. Excepting strand 4C, all of these regions mediate or propagate conformational changes. The F-helix and strand 5A must be displaced during protease inhibition, displacement of strand 1C is required for polymer formation, and helix D is a site (in antithrombin) of allosteric regulation. Our results indicate that these functionally important regions form a delocalized network of residues that are dynamically coupled, and that both local and global stability mediate inhibitory activity.

Keywords: Serpin, α1-antitrypsin, hydrogen/deuterium exchange, mass spectrometry, molten globule

The native state of most globular proteins is the conformation with the lowest Gibbs free energy. However, for some proteins an energy barrier prevents them from adopting the lowest free energy state, thus trapping them in a metastable state. Some typical examples are the membrane fusion proteins of some viruses, α-lytic proteases, arrestin and serpins. Serpins are a large and widely distributed family of serine protease inhibitors which include plasma protease inhibitors such as α1-antitrypsin, α1-antichymotrypsin, antithrombin III, plasminogen activator inhibitor-I, C1 inhibitor, as well as non-inhibitory members such as ovalbumin and angiotensinogen (1). Inhibitory serpins share a common tertiary structure consisted of three β sheets (A, B, C) and 8 to 9 helices, with the extended reactive center loop exposed at one end of the molecule for protease binding (Fig. 1a) (2). Upon cleavage by a protease, the amino terminal portion of the cleaved reactive center loop (RCL) of the serpin inserts into the central beta sheet A with the protease covalently linked to it (3). This conformational transition from stressed (S) (also referred to as “native”) form to relaxed (R) form results in a dramatic increase in thermal stability of the molecule (4). Unfavorable interactions in the native form such as side-chain overpacking, buried polar groups and cavities present in serpin molecule have been identified as the structural basis for native metastability (5-8) and such metastability thus turns out to be absolutely crucial for efficient protease inhibitor activity. The metastability of serpins and their ability to undergo controlled conformational changes also rendered these molecule convert to a latent inactive state (3). In the latent form, the RCL of the serpin molecule inserts into β-sheet A, therefore cannot react with the target protease. This also makes them vulnerable to misfolding and polymerization (9).

Figure 1
(A) Structure of the metastable form of α1AT. β-sheet A is shown in black and other functionally significant regions are indicated with arrows. (B) Close up view of the structure of G117F (21). Phe117 is shown in red and surrounding residues ...

Extensive mutagenesis studies have probed the role of local and global stability in serpins (10-14). Single stabilizing mutations near the RCL and in β-sheet A result in decreased activity, indicating that local instability in these regions is important for function (10-12). On the other hand, it has been shown that combinations of multiple stabilizing mutations can compromise activity regardless of location, provided that the degree of stabilization exceeds a threshold of ~13 kcal/mole (14). From these studies it appears that both local and global instability contribute to serpin function. One drawback to these studies is that, while they determined the effect of mutations on the stability of α1AT, they did not investigate their effect on conformational flexibility, which is also expected to play a role in inhibitory function.

Recently, the local distribution of both stability and conformational flexibility in α1AT was investigated using hydrogen/deuterium exchange and mass spectrometry (15-16). Conformational flexibility was found to vary widely between different regions, for example, the top of the F-helix is extremely labile while most of β-sheets A and B are rigid (15). In contrast, stability was much more uniformly distributed. The entire molecule undergoes a concerted transition to a molten globule state at 0.8 M GdnHCl, while the ΔG of unfolding for all regions clustered around a value of 5 kcal/ mole (16).

In the present study we have employed HXMS to study the local dynamics and local/global unfolding in a mutant of α1AT in which glycine 117 has been substituted with phenylalanine. Gly 117 is on β-sheet 2A and the Cα of the residue is surrounded by Phe-119, Tyr-160, and val-185 on the outer face of β-sheet A (Fig. 1b). This site forms a cavity that can accommodate larger side chains without appreciably altering the main structure. Substitution of Gly by the larger hydrophobic residue Phe fills a cavity between helix F and β-sheet A, and this substitution was shown to both increase the thermodynamic stability of α1AT and to reduce its inhibitory activity toward elastase (12). As the rate of association with elastase was unchanged, it was assumed that activity loss resulted from reduced efficiency of RCL translocation and insertion. Using this mutant, we tried to further elucidate the respective roles of conformational flexibility and stability in serpin function.

Materials and Methods

Protein expression, purification, and activity of G117F

The mutant protein was expressed and purified as described earlier (15). Activity assay of α1-AT was performed against trypsin (sigma) and elastase as dscribed (12, 15). Because active site titrants of PPE are not readily available, absolute activity was not determined. However, Wild Type and G117F assays were performed simultaneously against the same preparation of trypsin or PPE, and so the relative activity values are accurate.

Peptide mapping by HPLC MS/MS

Peptide mapping was performed as described. In short, 5 μg of G117F in 10 mM sodium phosphate (pH 7.8), 50 mM NaCl (Buffer A) was digested with 5 μg porcine pepsin dissolved in 0.05% (v/v) TFA and incubated for 5 min on ice. The digested protein was injected into a micropeptide trap (Michrom Bioresources) connected to a C18 HPLC column (5 cm × 1 mm, Altech) coupled to a Finnigan LCQ quadrupole ion-trap mass spectrometer (ThermoElectron) for the sequencing of each peptic peptide.

Equilibrium unfolding by H/D exchange

5 μg of the protein was incubated in buffer A containing different concentrations of GdnHCl for 1 hr at room temperature. The sample was then tenfold diluted with 10 mM sodium phosphate (pD 7.8) and 50 mM NaCl in D2O buffer containing the same concentration of Gdn2HCl as the pre-incubation sample. The exchange reaction was allowed to proceed for 10 s at room temperature followed by addition of 100 mM NaH2PO4 (pH 2.4) as the quench buffer. The final protein concentration was 0.5 μM. The quenched sample was then quickly frozen in a dry ice/ ethanol bath and stored at −80°C. The concentration of GdnHCl and GdnH2Cl were determined by refractrometry.

Conformational dynamics study by H/D exchange

5 μg of the protein in buffer A was 24 times diluted with the same buffer prepared in D2O at room temperature to label the sample. H/D exchange was terminated at different time points by the addition of the quench buffer (100 mM NaH2PO4 pH 2.4) and quickly frozen in a dry ice/ ethanol bath. The final protein concentration was same as that in equilibrium study. The samples were kept at −80°C until use.

Isotope analysis of HPLC electrospray ionization (ESI)

The frozen samples were quickly thawed and digested with 5 μg of porcine pepsin for 5 min. The digested sample was quickly injected into a micropeptide trap connected to a C18 HPLC column coupled to a Finnigan LCQ quadrupole ion-trap mass spectrometer. Peptides were eluted in 12 min using a gradient of 10-45% acetonitrile at a flow rate of 50 μl/ min. The micropeptide trap and column were all immersed in ice during the experiment to minimize back exchange.

Analyses of the unfolding curves were done according to the method described in 16. Percent deuterium exchange at different points was calculated as described previously (15).

Equilibrium unfolding in GdnHCl by CD spectroscopy

CD spectra were obtained using an Aviv 215 series CD spectrometer at 25°C with 1nm/10 s signal averaging from 210 nm to 250 nm using a 1mm path-length cell. The protein was incubated in buffer A containing different concentrations of optical-grade GdnHCl (Pierce) for 1 h at 25°C. The equilibrium unfolding curve was determined from the signal at 222 nm as a function of the concentration of denaturant. The data were fitted to a two state unfolding equation.

Results

The inhibitory activity of G117F against bovine trypsin was found to be 20% less than that of Wild Type. This is a smaller reduction than the 55-60% decrease seen against porcine pancreatic elastase that is reported in this work and elsewhere (12), but larger than the negligible change in activity seen against chymotrypsin (17).

Equilibrium unfolding of Wild Type α1AT by acid or chemical denaturants follows a three state transition model as revealed by CD spectroscopy (18-20). The biphasic nature of the Wild Type unfolding curve in Figure 2 indicates the presence of an intermediate. The unfolding profile for the mutant G117F was found to be different from Wild Type. The far UV CD spectral data showed a two state transition with the transition midpoint at 1.4 M GdnHCl (Fig. 2). The nature of the unfolding curve and transition midpoint of G117F were in good agreement with that found previously using fluorescence spectroscopy (12). The transition midpoint of 1.4 M is slightly higher than the midpoint of ~1.2 M found in (12), but this small difference may be due to the difference in pH (pD 7.8 in this study, pH 6.5 in (12)) as well as to the fact that unfolding was monitored by different techniques. For Wild Type α1AT intrinsic tryptophan fluorescence spectroscopy indicates two state unfolding (12) while unfolding monitored by far UV CD spectroscopy clearly reveals three state unfolding. In the case of G117F, both CD spectroscopy and fluorescence provide no evidence for an intermediate in the unfolding pathway.

Figure 2
GdnHCl induced equilibrium unfolding curves of Wild Type (●) and G117F (■) α1-AT as measured by CD spectroscopy. The far-UV CD signals were measured for samples equilibrated in varying concentration of GdnHCl. Data for Wild Type ...

To study site specific equilibrium unfolding, G117F was incubated in various concentrations of GdnHCl at pH 7.8 for 1 hour to attain equilibrium. After that, the protein was pulsed for 10 s with D2O buffer (25°C, pD 7.8) containing the same concentration of GdnCl as the pre-incubated samples. Under these conditions, amide hydrogens in unfolded regions of the molecule will undergo nearly complete exchange while hydrogens in folded regions will remain largely protected. This type of pulse labeling experiment has been shown to be an effective method for monitoring site specific folding and unfolding in proteins.

Fig. 3 shows the MS spectra of peptide 38-60 in the presence of 0.5 M, 1.4 M and 3.4 M guanidine. At 0.5 M GdnCl only one mass envelope was noted. At 1.4 M GdnCl two peaks of equal intensity were observed, one at lower and other at higher m/z value. The appearance of such bimodal peaks indicates that both folded and unfolded states are populated. As the concentration of the denaturant was increased, intensity of the peak at lower m/z started to decrease and disappeared above 3 M. Therefore, G117F was completely unfolded above 3 M guanidine. The unfolded population at each denaturant concentration was calculated as described in Methods and plotted against denaturant concentration. The transition midpoint obtained from the unfolding curve, fitted to a two state transition model (Fig. 4), matched well with that obtained from CD spectroscopic measurements. We performed similar analysis on all peptides obtained from peptic digests of G117F (Table 1). Due to noise caused by the presence of denaturant, some peptide signals were suppressed, and we present results for 13 peptides covering 42% of G117F. These peptides are, however, well distributed throughout the structure. The unfolding curves for all peptides followed a two state transition with midpoints clustered very tightly around 1.4 M GdnCl, in good agreement with the midpoint of 1.4 M obtained from fitting the CD unfolding curve of G117F, but significantly higher then the midpoint of 0.8 M for the native→molten globule transition in Wild Type. The coincidence of the unfolding curves determined from HXMS and CD indicates that, unlike Wild Type, G117F unfolds via a single cooperative transition in which both stable hydrogen bonding and secondary structure are lost simultaneously. Stabilization against chemical denaturation is increased quite uniformly throughout the structure. Denaturation midpoints are clustered around a median GuHCl concentration of 1.4 M with a standard deviation of 0.3 M. This is 0.55 M higher than the average transition midpoint of 0.85 M observed for Wild Type (16). A note about these transition midpoints is in order. We are comparing denaturation midpoints determined by HXMS. However, it was shown previously (16) that, for Wild Type α1AT, HXMS detects only the first of the two transitions seen by CD. Thus, we are here comparing the Wild Type native→molten globule transition with the G117F native→fully unfolded transition. Despite this complication, the fact that all peptides studied show a remarkably similar increase in their denaturation midpoints indicates that the effects of the G117F substitution are distributed fairly uniformly throughout the structure. Using the denaturation m-value of 7.9 kcal/mol*M from (12) we can calculate ΔΔG values for all peptides: these cluster around a median value of 4.2 kcal/mol. As with the shifts in transition midpoints, these ΔΔG values should be interpreted with caution, as they represent differences in ΔG between two qualitatively different transitions (native→molten globule and native→fully unfolded).

Figure 3
MS spectra of a representative peptic fragment 38-60 deuterated in 0 M (top), 1.4 M (middle), and 3.4 M GdnCl.
Figure 4
Equilibrium unfolding curve for peptic fragment 38-60 determined by H/D exchange mass spectrometry. The curve was fit to a two-state transition model.
Table 1
Denaturation midpoint of each peptic fragment determined by hydrogen-deuterium exchange mass spectrometry.

In addition examining the effects of the G117F substitution on the unfolding mechanism, we also probed the effects of the Gly→Phe substitution on the conformational dynamics of the native state. We investigated the distribution of flexibility and rigidity in the mutant by H/D exchange mass spectrometry in order to assess the structural differences from the Wild Type. H/D exchange was performed over period of 3000s at pD 7.8 at 25°C followed by HPLC-MS to quantify the mass of peptic fragments. Clear differences in rates of deuterium uptake are evident in several regions of α1AT. Figure 5 shows normalized deuterium uptake in the peptic fragment consisting of residues 190-208 for both Wild Type (open circles) and G117F (filled squares). By 500 seconds of incubation, the extent of deuterium uptake in Wild Type is ~35% greater than in G117F.

Figure 5
Raw data showing the percent deuterium exchange at different time points for the peptide fragment 190-208 derived from the Wild Type (○) and G117F (■) α1-AT.

Differences in normalized deuterium uptake between G117F and Wild Type at 10, 100, and 1000 seconds are shown for all peptides in Figure 6. Noticeable differences in flexibility were observed in peptides 160-171 and 160-172 containing six and seven residues in the C-terminal region of F helix. In the Wild Type this region showed exchange behavior comparable to that observed for unstructured surface loops such as the RCL. Exchange reached more than 80% in 100s time for the Wild Type, whereas it took 1000 s to show 85% deuterium exchange in G117F. Changes in the flexibility of the F-helix were anticipated as the mutation fills a cavity between the F helix and β strand A. Significant changes were also seen in other regions of the molecule. β-strands 5A (325-338), 1C (352-372), and 4C (190-208) all showed reduced exchange compared to Wild Type. The top of strand 5A, together with the top of strand 3A, comprises the “breach” region, which serves as the initial site of RCL insertion. Unfortunately, we lack coverage in strand 3A and so are unable to determine whether the entire breach region is stabilized in G117F.

Figure 6
Difference in normalized deuterium exchange in G117F relative to Wild Type. For each peptic fragment, differences are shown for 10 s (light blue bars), 100 s (purple bars) and 1000 s (dark blue bars) incubation in D2O.

In addition to helix F, parts of helix I (304-317) and helix D (85-99) also exhibited lower deuterium uptake than Wild Type. There are also some secondary structures in the mutant showing greater exchange behavior compared to the Wild Type. Residues 228-237 (comprised of β1B and some unstructured surface loops) and 64-77 (comprised of a loop and a portion of helix C) showed higher exchange behavior compared to the Wild Type at longer time points, although the showed less exchange at shorter time points. A plausible explanation for this behavior is that the least protected amide hydrogens in these regions (those that chiefly contribute to exchange at short labeling times) have been stabilized by the G117F substitution, while those amide hydrogens that are the most strongly protected against exchange are somewhat destabilized. Residues 385-395 also show slightly increased exchange at both short and long labeling times. Fig. 7 represents the different structural regions of α1AT which show significant changes in native state dynamics due to mutation. Conformational dynamics in other parts of the molecule remained more or less similar to that of the Wild Type.

Figure 7
Relative deuterium levels at 100 seconds in G117F compared to Wild Type mapped onto the three dimensional structure of α1AT. Both front (A) and back (B) views of the structure are shown. The figures were prepared using PyMol (37).

Discussion

It was previously demonstrated that the G117F substitution increases the stability of the metastable form of α1AT. Employing HXMS, we have found that this substitution also completely alters the unfolding mechanism. It is remarkable that a single amino acid substitution is sufficient to abolish the molten globule intermediate during equilibrium unfolding from the metastable state to the unfolded state. While position 117 has relatively low solvent accessibility, it cannot be considered as occupying the core of α1AT. Position 117 lies on the outside of β-sheet A, and the phenylalanine side chain will face outward, interacting primarily with the F-helix as revealed by the crystal structure (21). Further, the G117F substitution alters the unfolding mechanism not primarily by stabilizing the native state (although there is some stabilization relative to the native→molten globule transition of Wild Type) but by destabilizing the unfolding intermediate (Fig. 2). At 2.5 M GuHCl, for example, G117F is already fully unfolded as judged by CD, but Wild Type still retains considerable secondary structure. This is especially surprising since the introduction of a large nonpolar Phe side chain would, if anything, be expected to stabilize a compact denatured state. How this destabilization of the molten globule is achieved is unclear. Presumably the Phe side chain must disrupt stabilizing interactions in the intermediate state ensemble. From H/D exchange, we know that the intermediate lacks both stable secondary structure and a well packed hydrophobic core, but this does not rule out local or even long ranged side chain interactions. Lacking detailed information on the structure of the intermediate, we cannot say which specific interactions are disrupted by the G117F substitution. We are aware of one other instance in which 3 state unfolding is eliminated through destabilization of the intermediate rather than stabilization of the native state, and it is notable that this instance involves another member of the serpin family. The thermophilic serpin thermopin contains a C terminal extension not found in eukaryotic serpins. Chemical denaturation of thermopin shows 3 state unfolding. Deletion of the 5 C terminal residues eliminates the unfolding intermediate without substantial stabilization of the native state (22).

Seo et al. showed that the accumulation of multiple stabilizing mutations in α1AT results in activity deficient variants (14). They further found that substantial activity loss is accompanied by the disappearance of the unfolding intermediate (seen in Wild Type and modestly stabilized mutants) and the onset of two-state rather than three state unfolding. Having demonstrated using H/D exchange and other techniques that the intermediate seen in the three state unfolding of Wild Type α1AT is a molten globule, we proposed that the ability to populate a molten globule is correlated with, and perhaps necessary for, efficient protease inhibition by serpins (16). We suggested that the reason for this correlation is that α1AT must undergo transient disruptions throughout the structure, and not just at the loop insertion site, during inhibition. This is also supported by the work of Baek et al where they used H/D exchange mass spectrometry to show that α1AT undergoes transient unfolding in several structural regions – including regions distant from the loop insertion site – during protease inhibition (23). Now with G117F we have also demonstrated that decreased activity due to a stabilizing mutation is accompanied by loss of the ability to populate the molten globule state at low concentrations of denaturant. If large scale transient disruption of the folded structure is indeed required for RCL insertion during inhibition, then the altered unfolding observed here could influence activity through two mechanisms. On the one hand, the G117F mutation introduces favorable interactions that stabilize the folded structure. On the other hand, this mutation appears to destabilize the molten globule by an unknown mechanism. Both effects will increase the energetic distance between the folded state and the compact denatured state, thus increasing the energy barrier to be surmounted during the inhibitory conformational change (if in fact α1AT must pass through a compact denatured state during inhibition).

We note that the effects of the G117F mutation are not identical for all denaturants. Gooptu et al. found that primary effects of G117F on urea denaturation were to shift the onset of unfolding to higher concentrations (from 0.7 M for Wild Type to 2.0 M for G117F) and to shift the midpoint of transition (from 3.4 M to 4.3 M) (21). There was no clear indication that an unfolding intermediate had been eliminated. This is not surprising given that the unfolding behavior of Wild Type α1AT is quite different in urea and GuHCl. Monitored by CD spectroscopy, GuHCl denaturation of α1AT shows a very clear unfolding intermediate in which nearly 70% of the native ellipticity at 222 nm is retained (19). In contrast, the unfolding intermediate in urea is much less distinctly resolved by CD (it is more apparent to tryptophan fluorescence), and retains at most ~25% of its native ellipticity at 222 nm (24). The unfolding mechanism of α1AT (and perhaps other serpins) thus shows unusual plasticity: the intermediates populated (or not) and the degree of structure they retain can be significantly altered either by changing the denaturant employed or by amino acid substitutions.

In addition to its effects on the unfolding mechanism, the G117F substitution also affects local conformational dynamics. Based on the location of position 117 and the expected contacts between the F helix and the Phe side chain, the significant reduction in flexibility seen in the C terminal half of the F-helix is not surprising. Based on crystal structures, the F-helix blocks the path traversed by the RCL during translocation and must be physically displaced in order for complete loop insertion to occur. Therefore, it has an important role in serpin function (25-31). Previous H/D exchange data on Wild Type α1AT showed that the top of the F-helix is highly dynamic in solution, indicating weak interactions with beta-sheet A (15). It was suggested that these weak interactions would facilitate easy displacement of the F-helix and thus allow efficient loop translocation. The recently published 3.2 Å crystal structure of G117F showed that the F-helix is shifted towards the bottom of the β sheet 3A by half a turn. Also, introduction of the extra aromatic ring (Phe 117) into the helix F- β sheet A interface increases favorable packing interactions ultimately leading to the stabilization of helix F (21). This might also retard loop translocation (as F helix is more difficult to move aside) in G117F, thus leading to decreased activity.

Despite the fact that the crystal structure of G117F indicates no significant conformational changes other than those in and around the F-helix, we find that distant regions show altered dynamics in G117F. This is consistent with the proposition that allosteric communication between distant protein regions can occur without a discrete conformational change (32). Interestingly, the effects of the substitution are not uniformly distributed in the structure, despite the global effect on the unfolding mechanism, and in much of the structure, native state dynamics are not significantly altered by the mutation. Although in the crystal structure no significant changes were seen in β sheet A or β sheet C (21), regions showing significant changes in H/D exchange behavior include strand 5A, strands 1 and 4C and helix D. They all show reduced flexibility/increased stability in G117F. Strand 5A makes contacts with residues on the loop C terminal to the F helix which, in turn, makes contacts with the F helix itself. Stabilization could propagate through these interactions. More surprising is the significant stabilization of beta-strands 1C and 4C. These are ≥40 angstroms from position 117 and make no contacts with helix F or strand 5 A.

Excepting strand 4 C, all of the significantly effected regions are involved in regulating or propagating conformational changes, either in α1AT or in other serpins. Strand 5 A is part of the breach region and must separate from strand 3 A during RCL insertion. The release of strand 1 C from the rest of sheet C is critical for both the native→latent transition and for polymerization (33). G117F shows retarded polymerization in addition to decreased activity against some proteases (21). The most obvious explanation for this is strengthened interaction between the F-helix and the face of β-sheet A, as seen in the crystal structure. Our results suggest that stabilization of strand 1C may also contribute to retarded polymerization. Additionally, the observed stabilization of strand 5A is also consistent with retarded polymerization. Increased stability in strand 5A might discourage polymerization by making it more difficult to separate from strand 3A (in the loop sheet model of polymerization) or by making it more difficult to expel strand 5 from sheet A entirely (in the recently proposed domain swap model of polymerization (34). In serpins such as antithrombin, helix D plays a key role in allosteric activation; the affects of heparin binding to helix D are transmitted to the breach region and result in expulsion of the partially inserted RCL (4). Finally, large displacements are required of helix F during RCL insertion (30). The fact that the dynamical effects of the G117F substitution are concentrated in these functionally important regions suggests that they may form a network of dynamically and energetically coupled residues. Such coupling would facilitate the propagation of conformational changes between different regions of the molecule.

Proteins in solution exist as statistical ensembles, sampling many different conformations more or less frequently based on their relative energies. Hilser et al. have proposed that both local stability/flexibility and long range cooperative interactions in proteins are mediated by the global protein ensemble (35). In the context of ligand binding and allostery, this approach suggests that instability is critically important for the propagation of binding induced structural changes to distant regions in a protein's structure. Ensemble based models indicates that binding effects are propagated to distant regions most efficiently if a portion of the binding site is highly unstable in the absence of ligand and is stabilized upon ligand binding (36). This ensemble based view may aid in understanding the dramatic global affects of the G117F substitution inα1AT. Previous H/D exchange studies of Wild Type α1AT found that the top of the F-helix is one of the most unstable regions of the molecule, exchanging at rates comparable to an unstructured loop. The affect of introducing a Phe side chain at position 117 can be thought of as analogous to ligand binding to an unstable binding site in that it forms contacts with the F-helix and dramatically stabilizes it. Such stabilization would lead to a redistribution of the native state ensemble in a manner similar to that proposed for allosteric ligand binding. Applying an ensemble based computational algorithm to a large set of enzymes, they also found that functional residues showed stronger energetic coupling to the rest of the molecule than non-functional residues (37). α1AT is not an enzyme, and therefore the “functional residues” are not easily defined. α1AT's function requires a massive conformational change involving much of the molecule, and the intensive mutagenesis efforts of Yu et al. have established that “functional” residues are in fact distributed throughout the structure. This unusual property of α1AT (and perhaps of inhibitory serpins in general) may explain the unusually strong energetic coupling whereby local stabilizing interactions between sheet A and the F-helix propagate to distant regions and eliminate cooperative unfolding to the molten globule intermediate at low concentrations of denaturant.

This perspective in which energetic coupling between local regions is mediated by the overall protein ensemble suggests that both local and global changes in stability contribute to decreased activity against trypsin and elastase. Locally, strengthened interactions between the F-helix and β-sheet A will make displacing the F-helix more difficult, thus retarding loop insertion. However, this strengthened F-helix:sheet A interaction modulates the entire ensemble of α1AT, increasing the thermodynamic stability of the active state and apparently decreasing the stability of the compact intermediate. If distributed transient unfolding is required during the inhibitory transition, as the work of Baek et al. suggests, then this reordering of the relative stability of states will also retard loop insertion by increasing the energy barrier that must be overcome.

ACKNOWLEDGMENT

We thank Dr. Anne Gershenson (Brandeis University) for providing the plasmid encoding the G117F mutant of α1AT.

This work was supported by NIH grant R01HL085469.

ABBREVIATIONS

α1AT
α1-antitrypsin
HXMS
hydrogen/deuterium exchange/mass spectrometry
CD
circular dichroism
G117F
Gly117→Phe

REFERENCES

1. Potempa J, Korzus E, Travis J. The serpin superfamily of proteinase inhibitors: structure, function, and regulation. J. Biol. Chem. 1994;269:15957–15960. [PubMed]
2. Huber R, Carrell RW. Implications of the three-dimensional structure of alpha 1- antitrypsin for structure and function of serpins. Biochemistry. 1989;28:8951–8966. [PubMed]
3. Stratikos E, Gettins PG. Formation of the covalent serpin-proteinase complex involves translocation of the proteinase by more than 70 Å and full insertion of the reactive center loop into β-sheet A. Proc. Natl. Acad. Sci. 1999;96:4808–4813. [PubMed]
4. Gettins PG. Serpin structure, mechanism, and function. Chem. Rev. 2002;102:4751–4804. [PubMed]
5. Lee KN, Park SD, Yu MH. Probing the native strain in alpha 1-antitrypsin. Nat. Struct. Biol. 1996;3:497–500. [PubMed]
6. Lee KN, Im H, Kang SW, Yu MH. Characterization of a human alpha 1- antitrypsin variant that is as stable as ovalbumin. J. Biol. Chem. 1998;273:2509–2516. [PubMed]
7. Im H, Seo EJ, Yu MH. Metastability in the inhibitory mechanism of human alpha 1-antitrypsin. J. Biol. Chem. 1999;274:11072–11077. [PubMed]
8. Ryu SE, Choi HJ, Kwon KS, Lee KN, Yu MH. The native strain in the hydrophobic core and flexible reactive loop of a serine protease inhibitor: crystal structure of a uncleaved alpha 1-antitrypsin at 2.7 Å Structure. 1996;4:1181–1192. [PubMed]
9. Lomas DA, Evans DL, Finch JT, Carrell RW. The accumulation of Z alpha 1-antitrypsin in liver. Nature. 1992;357:605–607. [PubMed]
10. Seo EJ, Im H, Maeng JS, Kim KE, Yu MH. Distribution of the native strain in human alpha 1-antitrypsin and its association with protease inhibitor function. J. Biol. Chem. 2000;275:16904–16909. [PubMed]
11. Im H, Yu MH. Role of Lys 335 in the metastability and function of inhibitory serpins. Protein Sci. 2000;9:934–941. [PubMed]
12. Lee C, Park SH, Lee MY, Yu MH. Regulation of protein function by native metastability. Proc. Natl. Acad. Sci. 2000;97:7727–7731. [PubMed]
13. Lee C, Maeng JS, Kocher JP, Lee B, Yu MH. Cavities of alpha (1)-antitrypsin that play structural and functional roles. Protein Sci. 2001;10:1446–1453. [PubMed]
14. Seo EJ, Lee C, Yu MH. Concerted regulation of inhibitory activity of alpha 1-antitrypsin by the native strain distributed throughout the molecule. J. Biol. Chem. 2002;277:14216–14220. [PubMed]
15. Tsutsui Y, Liu L, Gershenson A, Wintrode PL. The conformational dynamics of a metastable serpin studied by hydrogen exchange and mass spectrometry. Biochemistry. 2006;45:6561–6569. [PubMed]
16. Tsutsui Y, Wintrode PL. Cooperative unfolding of a metastable serpin to a molten globule suggests a link between functional and folding energy landscapes. J. Mol. Biol. 2007;371:245–255. [PubMed]
17. Parfrey H, Mahadeva R, Ravenhill NA, Zhou A, Dafforn TR, Foreman RC, Lomas DA. Targeting a surface cavity of α1-antitrypsin to prevent conformational disease. J. Biol. Chem. 2003;278:33060–33066. [PubMed]
18. Devlin GL, Chow MK, Howlett GJ, Bottomley SP. Acid denaturation of alpha 1-antitrypsin: characterization of a novel mechanism of serpin polymerization. J. Mol. Biol. 2002;324:859–870. [PubMed]
19. James EL, Whisstock JC, Gore MG, Bottomley SP. Probing the unfolding pathway of alpha 1-antitrypsin. J. Biol. Chem. 1999;274:9482–9488. [PubMed]
20. Tew DJ, Bottomley SP. Probing the equilibrium denaturation of the serpin alpha 1-antitrypsin with single tryptophan mutants; evidence for structure in the urea unfolded state. J. Mol. Biol. 2001;313:1161–1169. [PubMed]
21. Gooptu B, Miranda E, Nobeli I, Mallya M, Purkiss A, Leigh Brown SC, Summers C, Phillips RL, Lomas DA, Barrett TE. Crystallographic and cellular characterization of two mechanisms stabilizing the native fold of alpha (1)-antitrypsin: implications for disease and drug design. J. Mol. Biol. 2009 Epub ahead of print. [PMC free article] [PubMed]
22. Irving JA, Cabrita LD, Rossjohn J, Pike RN, Bottomley SP, Whisstock JC. The 1.5 Å crystal structure of a prokaryote serpin: controlling conformational change in a heated environment. Structure. 2003;11:387–397. [PubMed]
23. Baek J, Yang WS, Lee C, Yu MH. Functional unfolding of α1-antitrypsin probed by hydrogen-deuterium exchange coupled with mass spectrometry. Mol. Cell. Proteomics. 2009;8:1072–1081. [PMC free article] [PubMed]
24. Tew DJ, Bottomley SP. Probing the equilibrium denaturation of the serpin α1-antitrypsin with single tryptophan mutants: evidence for structure in the urea unfolded state. J. Mol. Biol. 2001;313:1161–1169. [PubMed]
25. Vleugels N, Gils A, Bijnens AP, Knockaert I, Declerck PJ. The importance of helix F in plasminogen activator inhibitor-1. Biochim. Biophys. Acta. 2000;1476:20–26. [PubMed]
26. Bijnens AP, Gils A, Knockaert I, Stassen JM, Declerck PJ. Importance of the hinge region between alpha-helix F and the main part of serpins, based upon identification of the epitope of plasminogen activator inhibitor type1 neutralizing antibodies. J. Biol. Chem. 2000;275:6375–6380. [PubMed]
27. Cabrita LD, Whisstock JC, Bottomley SP. Probing the role of the F-helix in serpin stability through a single tryptophan substitution. Biochemistry. 2002;41:4575–4581. [PubMed]
28. Cabrita LD, Dai W, Bottomley SP. Different conformational changes within the F-helix occur during serpin folding, polymerization, and proteinase inhibition. Biochemistry. 2004;43:9834–9839. [PubMed]
29. Baek JH, Im H, Kang UB, Seong KM, Lee C, Kim J, Yu MH. Probing the local conformational change of α1-antitrypsin. Protein Sci. 2007;16:1842–1850. [PubMed]
30. Gettins PG. The F-helix plays an essential, active role in the proteinase inhibition mechanism. FEBS Letters. 2002;523:2–6. [PubMed]
31. Zheng X, Wintrode PL, Chance MR. Complementary structural mass spectrometry techniques reveal local dynamics in functionally important regions of a metastable serpin. Structure. 2008;16:38–51. [PubMed]
32. Tsai CJ, Del Sol A, Nussinov R. Allostery: absence of a change in shape does not imply that allostery is not at play. J. Mol. Biol. 2008;378:1–11. [PMC free article] [PubMed]
33. Chang W.S. Whisstock J., Hopkins P.C. Lesk A.M., Carrell RW, Wardell MR. Importance of the release of strand 1C to the polymerization mechanism of inhibitory serpins. Protein Sci. 1997;6:89–98. [PubMed]
34. Yamasaki M, Wei L, Johnson D, Huntington JA. Crystal structure of a stable dimer reveals the molecular basis of serpin polymerization. Nature. 2008;455:1255–1258. [PubMed]
35. Hilser VJ, Garcίa-Moreno E,B, Oas TG, Kapp G, Whitten ST. A statistical thermodynamic model of the protein ensemble. Chem. Rev. 2006;106:1545–1558. [PubMed]
36. Freire E. The propagation of binding interactions to remote sites in proteins: analysis of the binding of the monoclonal antibody D1.3 to lysozyme. Proc. Natl. Acad. Sci. 1999;96:10118–10122. [PubMed]
37. Liu T, Whitten ST, Hilser VJ. Functional residues serve a dominant role in mediating the cooperativity of the protein ensemble. Proc. Natl. Acad. Sci. 2007;104:4347–4352. [PubMed]
38. Delano Scientific . The PUMOL molecular graphics system. San Carlos, CA: 2002. http://www.pymol.org.