Hydrogen–deuterium exchange reveals long-range dynamical allostery in soybean lipoxygenase

In lipoxygenases, the topologically conserved C-terminal domain catalyzes the oxidation of polyunsaturated fatty acids, generating an assortment of biologically relevant signaling mediators. Plant and animal lipoxygenases also contain a 100-150-amino acid N-terminal C2-like domain that has been implicated in interactions with isolated fatty acids and at the phospholipid bilayer. These interactions may lead to increased substrate availability and contribute to the regulation of active-site catalysis. Because of a lack of structural information, a molecular understanding of this lipid-protein interaction remains unresolved. Herein, we employed hydrogen–deuterium exchange MS (HDXMS) to spatially resolve changes in protein conformation upon interaction of soybean lipoxygenase with a fatty acid surrogate, oleyl sulfate (OS), previously shown to act at a site separate from the substrate-binding site. Specific, OS-induced conformational changes are detected both at the N-terminal domain and within the substrate portal nearly 30 Å away. Combining previously measured kinetic properties in the presence of OS with its impact on the Kd for linoleic acid substrate binding, we conclude that OS binding brings about an increase in rate constants for both the ingress and egress of substrate. We discuss the role of OS-induced changes in protein flexibility in the context of changes in the mechanism of substrate acquisition.

Lipoxygenases, widely represented among plants, animals, fungi, and bacteria, catalyze the C-H abstraction from longchain polyunsaturated fatty acids in the process of generating a diverse repertoire of (per)oxidized products that play important roles in various signaling pathways (1,2). Select eicosanoid products from human enzymes have been implicated in a number of inflammation-based diseases (examples include arthritis, asthma, atherosclerosis, and cancer) (3); lipoxygenases also represent a "neglected frontier for regulation of chronic pain" (4). Because of the role of lipoxygenases in a wide range of signaling pathways, an understanding of allosteric regulation of active-site chemistry is fundamental to a mechanistic framework of enzyme reactivity as well as the success of future rational drug design efforts.
Lipoxygenases are structurally characterized by a 500 -700amino acid (predominantly ␣-helical) C-terminal domain (Fig.  1). The C-terminal domain contains the catalytic metallocofactor, typically a ferric hydroxide that is responsible for catalysis (Fig. 2). Although there is considerable variability in sequence and loop features, the structural topology of this domain is fairly conserved; however, a 100 -150-amino acid N-terminal polycystin-1/lipoxygenase/␣-toxin (PLAT) 3 domain is only found in animals and plants (5)(6)(7). Of the plant lipoxygenases, soybean lipoxygenase-1 (SLO) has emerged as the paradigm of the family due to its native abundance, ease of purification, and structural stability. Because SLO and human lipoxygenases catalyze similar reactions and share similar structural elements, SLO is often related to the molecular enzymology of human enzymes as well.
Whereas the function of the C-terminal domain is well established, the role for the N-terminal domain has not been fully elucidated. The N-terminal PLAT domain is structurally reminiscent of C2 domains in proteins that have been established to bind calcium and associate with the lipid membrane bilayer for cell signaling (8,9). Studies on the human lipoxygenases have implicated the role of this PLAT domain in Ca 2ϩ -mediated phospholipid membrane binding (10 -14).
In the case of plant enzyme SLO, calcium does not regulate catalysis or mediate membrane binding (15); instead, isolated fatty acids, such as the substrate linoleic acid (LA) at high concentrations, are postulated to associate at a remote site, presumably within the protein's N terminus, and alter the interaction of substrate within the catalytic center (16 -18). Effectors such as oleyl sulfate (OS) and mixed inhibitors, such as oleic acid (OA), have been shown to stimulate the same allosteric behavior (18,19), suggesting that LA, OA, and OS (and perhaps the phospholipid membrane (16)) may all interact at (or near) a remote allosteric site. Biochemical analysis of the product feedback inhibition found in humans, although lacking in SLO, has further supported the presence of an allosteric site for isolated lipid effector interactions within or near the PLAT domain (20,21). Thus, although an ongoing recognized hypothesis is that the N-terminal domain of lipoxygenase interacts with lipid effectors to promote changes in substrate binding within the C-terminal catalytic domain, a network for communication between the allosteric and active sites has yet to be resolved.
To extend the molecular framework for lipid-induced allostery in lipoxygenases, we present HDXMS experiments with SLO in the absence and presence of the effector, OS. HDXMS has emerged as an incisive tool to describe networks of allosteric communication in proteins and to offer molecular insight into the mechanism of enzyme inhibition and/or activation (22)(23)(24)(25)(26)(27). In the context of remote binding effects that alter lipoxygenase behavior, we find that SLO may serve as a tractable system, due to its structural stability and previously well-characterized kinetic impact from effectors such as OS. Further, we recently reported conditions in which HDXMS is capable of spatially resolving remote changes in SLO flexibility that correlate with the introduction of catalysis-impairing active-site mutations (28). Consistent with evidence from lipid interaction studies on SLO and human lipoxygenases (for examples, see Refs. 16 and 21), this current HDXMS study supports a localization of OS, as a lipid surrogate, at the N terminus (within the PLAT domain). Further, six non-overlapping peptides within the catalytic domain display a significant increase in the observed rate of deuterium exchange upon OS addition. This change in apparent rates of HDX supports a model in which lipid interactions at the N terminus are linked to the destabilization of a group of ␣ helices stemming from the N terminus to the substrate entrance portal. Integrating prior kinetic analyses of OS inhibition (18) with measurement of the impact of OS on the K d for linoleic acid substrate has enabled us to estimate the degree to which rates for LA binding and release are increased, implicating a shift in an internal equilibrium constant describing open and closed ES complexes and attributed to the aforementioned helical destabilization. Taken together, the data highlight the ability of HDXMS to interrogate the molecular mechanism of long-range allosteric effects in lipoxygenases.

Results
The conditions of the described HDXMS experiments are at temperatures ranging from 10 to 40°C and effective pD 8 (see "Experimental procedures" for details), where the rate-limiting step for k cat has been shown to involve the chemical C-H bond cleavage (29) and HDX has been shown (28) to proceed via the EX-2 exchange mechanism (30,31). EX-2 conditions permit a temperature-dependent analysis of regional conformational flexibility (28,32,33). To streamline data analysis, 45 non-overlapping peptides (Fig. 3) were selected from the nearly 300 available for this 839-amino acid protein, providing 63 and 89% sequence coverage of the N-and C-terminal domains, respectively. For HDXMS experiments with OS, the peptide sequence map is identical to that for SLO in the absence of effector. HDXMS was performed for both sets of peptides as a function of time (10 s to 4 h) and temperature (10, 20, 30, and 40°C); previous controls showed that wildtype SLO is structurally stable, contains nearly complete iron occupancy, and maintains high activity (Ͼ90%) during the course of these incubations (28). For each peptide, the percentage of deuterium exchange (or incorporation) was determined by calculating the difference between the mass of the peptide before and after a given incubation time; the data in the presented HDXMS traces (Figs. 4 and 6 and Fig. S1) have been corrected for peptide-specific back-exchange during work-up.
Two sets of experiments were conducted, both on wildtype SLO, in which one set of samples was incubated in the presence of 10 M OS and another set was incubated without OS. Although OA can act as an allosteric effector (17), we chose to focus on OS to preclude competition between OA and LA at the substrate-binding site. As shown previously, OS does not compete with substrate binding and exclusively targets an alternate site in SLO (19). The concentration of OS was chosen to be  Left, reaction of model enzyme, soybean lipoxygenase, with substrate linoleic acid as an example. The metallocofactor, often a ferric hydroxide, although a manganese center is employed by fungal enzymes (63), abstracts a hydrogen atom from polyunsaturated fatty acids, in a regio-and stereospecific manner, through a proton-coupled electron transfer (PCET) mechanism (49). Molecular oxygen inserts into the delocalized radical intermediates, forming the product upon reverse proton-coupled electron transfer. The carbon numbering of the LA substrate is shown for reference. Relevant microscopic rate constant, k 2 , associated with the rate-determining C-H abstraction, is labeled. The enzyme-substrate and enzyme-LA radical complexes are labeled as ESЈ and ES ⅐ , respectively. Right, the structures of substrates, as well as the allosteric effectors OA and OS, are shown for reference. The reactive carbon for each substrate is designated by an asterisk.

Allostery in soybean lipoxygenase
saturating (K d ϳ0.6 M) while remaining below the critical micelle concentration (Ն18 M) (18). A thorough inspection of the impact of OS on HDXMS, using pattern recognition (28), revealed that for the majority (38 of the 45 peptides), protein flexibility was unaffected by OS (see Fig. S1 for the complete set of HDXMS traces). Among the remaining peptides, two regionally defined behaviors emerged from this data analysis, with OS observed either to increase the extent of exchange at the N terminus ( Fig. 4) or to increase the rate of exchange in the C-terminal domain (see Fig. 6 and Table 1).

Conformational change in the N-terminal PLAT domain
We first searched for changes in the extent of the deuterium exchange (at 4 h) in the presence of OS with the expectation of observing a canonical "protection" effect (i.e. decreased HDX) due to formation of a complex of OS with protein. In contrast, at all temperatures examined, OS displayed a significantly increased extent of exchange at a single, N-terminal peptide, 1-14 (Fig. 4). The difference in percent exchange, between SLO/OS and SLO alone, shows its largest impact at elevated temperatures (30 and 40°C). Based on the percent exchange and number of exchangeable amide linkages (13 for this peptide), we calculate that OS induces a perturbation to 1-3 amide hydrogen bonds (depending on temperature) within the N-terminal peptide, 1-14. Comparison with the next sequential peptide, 15-44, showed no noticeable change to the extent or rate of exchange. We further examined peptides that overlap with 1-14, namely 1-20 and 1-32 (Fig. 4). As anticipated, the extents of exchange are increased in the presence of OS. The differences in percent exchange are smaller for these longer peptides, although the number of altered amide hydrogen bonds remains the same (i.e. 1-3, depending on temperature). These results suggest that OS binds to or in the vicinity of residues 1-14 and introduces a distinct change in conformation within the SLO N-terminal domain (localized around peptide [1][2][3][4][5][6][7][8][9][10][11][12][13][14]. Our inability to detect any protection from HDX by OS at its actual binding site is probably due to relatively poor peptide coverage within the N-terminal region of protein, only 63% compared with 89% coverage within the C-terminal catalytic domain. From analysis of the X-ray-derived structure (Fig. 5A), the first six amino acids (MFSAGH) could not be modeled, suggesting a high degree of structural mobility in crystallo. The rest of this peptide sequence is predominantly hydrophobic, aside from two lysine residues (KIKGTVVL), and forms ␤-sheet secondary structure. As noted above, further localization of the OS interaction was not feasible because peptides within the ␤-barrel that sandwich the peptide 1-14 through interstrand hydrogen bonding are not covered (Fig. 5A, gray schematic) in our HDXMS measurements. Attempts by our group to soak crystals or co-crystalize with OS were unsuccessful; no electron density could be assigned to OS (or even OA) at the N terminus or anywhere else in the protein structure.
It is possible that the negative charge on the sulfate group of OS alone is sufficient to interact with one or both of the above

Allostery in soybean lipoxygenase
highlighted lysine residues to elicit a conformational change. However, these residues are not conserved among plant and animal lipoxygenases, and OS has been shown to bind with comparable affinity to the human 15-lipoxygenase (18). Fur-ther, previous stopped-flow kinetic experiments with palmitoyl sulfate, a compound structurally similar to OS with 2 fewer carbons, did not exhibit a high specificity for the allosteric site (K i ϳ140 M, compared with 0.6 M for OS) (19). Both oleic (18:1) and linoleic (18:2) acid also bind at the allosteric site (K d Ն 20 M), although not nearly as tightly as OS (18). The length of the molecule and the structure of the headgroup probably play a role in the tight association of OS with the allosteric site.

Destabilization of C-terminal ␣ helices
The second revealing pattern in the data is the observation of six non-overlapping SLO peptides whose apparent rate of exchange is increased upon incubation with OS ( Fig. 6). HDX is typically regarded as a multiexponential process with slow, intermediate, and fast kinetic exchange rate constants. Whereas many reports generally use three exponentials to fit data within the dynamic range, we have found, using various multiexponential fitting functions, that the SLO HDXMS traces can be accurately modeled with two or even one exponential kinetic process without over-or underfitting data. In the case of peptides with two exponentials, the reported apparent rate constants represent their weighted average. A complete comparison of apparent rate constants can be found in Table S1.  Table 1 Apparent average rates (min ؊1 ) of exchange, k HDX , from fits to HDXMS data (at 10°C) for peptides affected by OS Apparent average rate constants were determined from one-or two-exponential fits to the data within the dynamic range of the measurement. The rates are represented within one S.D. from the exponential fits to the data. For complete rate constant analysis across all temperatures, see Table S1; Arrhenius plots are shown in Fig. 7

Allostery in soybean lipoxygenase
Comparing rate constants across all temperatures revealed three peptides with rates of exchange increased by at least 3-fold: 297-305, 306 -316, and 317-334 (cf. Table S1). The apparent rates of exchange for these peptides were previously shown to be sensitive to active-site mutations (28). A more thorough inspection of rate constants showed that three additional peptides exhibit altered rates, but only at lower incubation temperatures: 239 -256, 257-273, and 751-761 (Table 1). These peptides have previously been characterized as having rapid exchange approaching nearly full extent. Indeed, at elevated temperatures (30 and 40°C), the exchange rates for the latter three peptides (Table S1) are rapid and comparable for SLO in the presence and absence of effector (Ն3 min Ϫ1 ). However, at 10°C, the exchange process for the SLO protein in the absence of effector is slowed, sufficiently to be monitored versus time; thus, we focus on the HDXMS results at 10°C, where the largest exchange differences are observed. In the presence of OS, a substantial increase in the apparent rate of exchange (by Ն8-fold) is observed for all six peptides at 10°C (Fig. 6 and Table 1). Peptide 257-273 (and perhaps 751-761) is the most impacted by the presence of OS. Overall, the increase in the rates of exchange for all peptides shown in Fig. 6 is consistent with a destabilization of their secondary structure (34).
Besides the peptide 317-334, which is modeled primarily as a loop structural element in the X-ray structure, the remaining five C-terminal peptides affected by OS are ␣ helices (Fig. 6). Our collective HDXMS data support a model of OS binding at the N-terminal PLAT domain that triggers the destabilization of these ␣ helices in the catalytic domain. Previously, Tatulian et al. (16) provided evidence from the vibrational shifts in the peptidyl amide I vibrational modes that linoleic acid induces destabilization of ␣ helices (i.e. to a more open structure). It has been shown that at elevated concentrations, substrate LA and mixed inhibitor, OA, associate at both an allosteric site and the active site in SLO. As stated above, OS exclusively interacts at the allosteric site of SLO and alters the binding of LA at the active site. Our HDXMS results, when combined with studies of OS binding (18,19), support the location of an allosteric site for SLO in the PLAT domain that leads to destabilization of ␣ helices within the catalytic domain.
Assessment of the rates of HDX for the peptides affected by OS as a function of the four temperatures studied here further revealed striking changes in the apparent activation energies for exchange; the Arrhenius plots are displayed in Fig. 7. A change in the Arrhenius slope is related to an altered activation energy of exchange. For EX-2 conditions, the observed rate of exchange (k HDX ) will represent the product between the equilibrium constant for opening (k open ) and closing (k close ) processes, K op ϭ k open /k close , and the intrinsic chemical exchange rate, k int : k HDX ϭ K op ϫ k int (30,31,35). The latter term represents the exchange rate on a fully exposed amide bond and depends on inductive effects from side-chain composition. In the absence of OS (Fig. 7, black dots/lines), the apparent activation energies of HDX for peptides shown (15 Ϯ 2 Ͻ E aHDX Ͻ 20 Ϯ 2 kcal/mol) are within error of the enthalpy for the intrinsic exchange (⌬H ‡ int ϭ 17 kcal/mol) (64). Upon OS addition, the first clear trend is that the exchange rates for N-terminal peptide, 1-14, are unaffected across all temperatures (Fig. 7A). The same is true for all other N-terminal domain peptides (Table S1). The second feature to emerge

Allostery in soybean lipoxygenase
is a drastic change in the Arrhenius slopes for 239 -256 and 257-273 (Fig. 7, B and C). For these two peptides in the presence of OS, E aHDX(avg) values are considerably reduced; apparent activation energies of exchange decreased from the intrinsic value of 17 kcal/mol, requiring negative values for K op that can be interpreted in the context of increased disorder within this region upon remote binding of OS (36,37). For the remaining three peptides affected by OS, 297-305, 306 -316, and 317-334 (Fig. 7, D-F), Arrhenius slopes are seen to be only modestly altered by OS when compared with the behavior of 239 -256 and 257-273. The temperature dependence of the exchange rates for peptide 317-334 is perhaps the least affected. This behavior is in contrast to a previous observation from activesite mutational impact, in which this loop displayed a linear relationship between the enthalpic barrier for HDX and the corresponding barrier for catalytic PCET (28). An analysis of the temperature dependence of 751-761 has been omitted, as reliable rate constants were not obtainable at temperatures above 20°C.
We point out that analyses of HDXMS traces at a single temperature, for example 30°C in which most of the k cat , K m , and D k cat kinetic parameters are often reported for SLO, would have revealed only changes to the loop region including peptides within the 297-334 segment. The impact on conformational flexibility at peptides 239 -256, 257-273, and 751-761 is masked at this temperature because the exchange rates are practically identical for both conditions interrogated here (absence versus presence of OS). Whereas peptides including loop 317-334 and those spatially adjacent (297-305 and 306 -316) are seen to have slightly impacted temperature dependences for the apparent rate constants, the change in the temperature dependence is most pronounced for 239 -256 and 257-273. These results illustrate the advantage of performing HDXMS experiments at multiple temperatures.
HDXMS experiments conducted with OA (collected at 30°C and with 45 M OA) display trends similar to those for OS at the single temperature analyzed (Fig. S2). OA has also been shown to bind at the active site (38); however, OA is poorly reactive with ferric SLO (39) and non-reactive with the ferrous state of the enzyme (40) studied here. Whereas OA is expected to interact at both the allosteric and active sites, HDXMS traces at 30°C for OA and OS are nearly indistinguishable (Fig. S2). Most notably, OA causes a similar conformational change at the N-terminal peptide, 1-14; these data are consistent with association of lipids, both OS and OA, at the proposed allosteric site located within the N-terminal PLAT domain. The observation of similar HDXMS behavior between free and substrate-bound protein has been documented previously for other systems (41,42). HDXMS experiments with SLO were not performed in the presence of substrate, LA, because of the complicated results that were expected from structural alterations occurring during the catalytic cycle.

Discussion
HDXMS is proving to be a vital tool for the interrogation of protein flexibility associated with long-range allosteric communication, capable of detecting subtle changes in backbone mobility that are often invisible with X-ray crystallography. This work represents the HDXMS interrogation of lipid-induced allostery in lipoxygenase. Here, two distinct patterns of HDX have emerged from interactions of the lipid-like allosteric effector, OS, with the model plant enzyme, SLO. First, in the N-terminal PLAT domain, found only in humans and plants, OS induces a significant increase in the overall extent of HDX at long incubation times, consistent with the additional exchange of 1-3 amide hydrogens. We interpret this effect in the context of a conformational change, associated with OS binding nearby. Further, six peptides in the catalytic domain, where substrate binds and is oxidized, are seen to undergo dramatic increases in flexibility, supporting destabilization of their secondary structure.
Typically, protein ligand-inhibitor interactions show patterns of rates and/or extents of HDX that are reduced from the native protein (35,43,44). This canonical ligand-induced behavior is often attributed to enhanced intramolecular hydrogen bonding (i.e. structural stabilization). Protection of amide hydrogen bonds from exchange, as read out in HDXMS experiments as decreased exchange (extent and/or rates), is widely accepted as a hallmark for protein-protein and protein-ligand interactions. However, protein-ligand binding can exhibit enhanced rates and/or extents, as documented previously (for examples, see Refs. [45][46][47], and these patterns are characteristic of a different type of HDX behavior (48). This latter type of HDX behavior in which ligand binding increases HDX rates has been suggested to arise when the ligand binds following (or

. Arrhenius-like plots of the weighted average exchange rates, ln(k HDX(avg) ), for SLO in the absence (black) or presence (gray) of OS.
Apparent rates of HDXMS were determined from fits to the traces in Fig. S1. All rate constants can be found in Table S1.

Allostery in soybean lipoxygenase
followed by) 4 a protein structural change (48). Such may be the case for OS association at the SLO PLAT domain, where the presence of this effector shifts a protein conformational equilibrium toward a more solvent-exposed state within both the N terminus and catalytic domains.
The previously characterized kinetic ramifications from OS binding provide a mechanistic framework for the HDXMS results presented herein (18,50). Select kinetic parameters for the SLO reactions on substrate LA in the absence and presence of OS are summarized below in Table 2. Importantly, an increase in the D k cat /K m was shown upon the addition of OS; this has been linked to an increased microscopic rate constant for the substrate off rate, k off , by 60-fold compared with the absence of OS (see Table 2 for details). From the estimated off rates, the relative substrate on rates, k on , can be determined from the dissociation constants, K d ϭ k off /k on . Because C-H bond cleavage is fully rate-limiting for deuterated substrate (29), K d is equivalent to the Michaelis-Menten constant for the deuterated substrate (i.e. K m (D) ϭ K d ) measured herein. Estimated k on values, reported in Table 2, are seen to be accelerated by ϳ12-fold in the presence of OS. These microscopic rate constants suggest that substrate diffusion is no longer partially rate-limiting under these conditions. This effect is corroborated by the demonstrated elimination of the solvent viscosity dependence ((k cat /K m°) /(k cat /K m )) in the presence of OS (see Table 2) (18). We note that the summarized kinetic properties of SLO (Table 2) represent a pH value slightly different from that of the HDXMS analyses (pH of 9 versus pD 8, respectively). From previous studies, the k cat and D k cat parameters are seen to be identical at pH 8 and 9, whereas k cat /K m undergoes only a modest change; further, the value of D k cat /K m at pH 9 is within experimental error of the value at pH 8 (29). For these reasons, it has been possible to quantify the kinetic impact of OS at pH 9, where the majority of kinetic data have been collected for SLO, and relate its behavior to HDXMS at pD 8, conditions where EX2 behavior has been demonstrated (28).
To explain the altered kinetics for substrate binding, in which k off is increased to a larger extent than k on , we propose a two-step mechanism for productive enzyme-substrate complex formation with an internal reorganization ES^ESЈ associated with the equilibrium constant, K eq (Fig. 8). In the native system without OS, K eq favors ESЈ so that the rate-determining steps for k cat /K m include k 2 (chemistry) and substrate diffusional control. In the presence of OS, substrate diffusion no longer (significantly) contributes to the rate-determining steps (i.e. there is no viscogen dependence, (k cat /K m°) /(k cat /K m ) ϭ 1). Because k on and k off are both dependent on K eq (i.e. k on ϭ k 1 [S] ϫ K eq ; k off ϭ k Ϫ1 ϫ 1/K eq ; see Fig. 8) and OS addition impacts the off rate 5 times more than the on rate, the data imply that K eq is reduced in the presence of OS, shifting this internal equilibrium toward a "more open" form (denoted as ES).
We propose that the enhanced substrate on and off rates may be facilitated by the destabilization of peptide 257-273, which is anticipated to line the substrate portal in SLO (51,52). This putative entrance for substrate is located over 30 Å from the N-terminal peptide, 1-14, whose OS-induced conformational change is linked to enhanced substrate on and off rates (Fig.  9A). This raises the question: What is the link between the OS-induced conformational change at the N-terminal PLAT domain and the destabilization of ␣ helices in the catalytic domain? A conserved cationinteraction between Trp-130 and Arg-242 (Fig. 9, A and B; SLO numbering) has been previously proposed to serve as a mediator of Ca 2ϩ activation in mammalian enzymes (13). Further, the homologous tryptophan in the human 5-lipoxygenase lies within a segment (FPCYRW 102 ) that has been speculated to interact with the 4 It is important to note here that the exact order of these events (binding of molecules induces conformational change, or conformational change induces ligand binding) is difficult to clarify because the HDXMS experiment as assessed here is a thermally averaged measurement under EX-2 conditions. The data are consistent with a conformational change that occurs with ligand (i.e. lipid) binding.  The addition of effector, OS, causes a relative reduced barrier for the first binding step and a change in the internal reorganization. K eq refers to this proposed internal reorganization from the initial enzyme-substrate complex (ES) to ESЈ structure, with the latter representing the productive binding for chemistry. The substrate on rate (k on ϭ k 1 [S] ϫ K eq ) and off rate (k off ϭ k Ϫ1 ϫ 1/K eq ) are both a product of this internal equilibrium constant.

Allostery in soybean lipoxygenase
phospholipid membrane and an activating protein, coactosinlike protein, and elicit conformational allostery (53,54). The HDXMS-identified network of allosteric communication in SLO, stemming from the PLAT domain to the substrate binding portal, may also be facilitated by this structural feature.
The current HDXMS and kinetic study was carried out at pH values elevated in relation to physiological conditions. The rationale for the conditions used here was to develop a molecular framework for the regulation of substrate diffusion by allosteric effectors. The key finding from the current work is the implied localization of lipid (surrogate) binding interactions at the N-terminal PLAT domain together with long-range allostery that arises from the lipoxygenase-lipid interactions. Beyond the present results, the most substantiated evidence to implicate the association of N-terminal PLAT domains with the free lipids or the lipid membranes can be attributed to N-terminal deletions. Examples from select plant and mammalian lipoxygenases suggest that removal of the PLAT domain impairs stability and catalytic proficiency and affects allosteric regulation (see Refs. 21, 55, and 56 and references within).
At the time of submission of this work, an HDXMS study of human 15-lipoxygenase, 15-LOX-2, appeared in the literature (14). This study showed little impact of either Ca 2ϩ or lipid nanodiscs on HDXMS, whereas substrate was shown to yield a pattern of protection within the protein overall. Important differences between that study and the current work are the reduced coverage of the protein sequence for 15-LOX-2 in relation to SLO and the focus in the 15-LOX-2 study on a single temperature and different time regimes. As we show herein, changes in temperature can have a dramatic impact on observed patterns of HDXMS (Figs. 4 and 7), expanding our capacity to infer mechanistic features.
It should be noted that whereas Ca 2ϩ triggers translocation of human LOX to the phospholipid membrane, calcium does not regulate membrane binding for SLO (15). Human enzymes also exhibit product feedback inhibition, whereas SLO does not. Together with the fact that SLO is 20% larger than human LOXs, the data suggest that 15-LOX-2 and SLO behave as different model enzymes.
To conclude, an important feature of the presented work is the focus on the interactions of SLO with OS, a molecule that neither turns over with SLO nor binds at the active site. These properties have allowed us to isolate the allosteric properties of a lipid effector as monitored by HDX and steady-state kinetics. The enhanced HDX observed for the SLO PLAT domain provides strong support for OS binding to the N terminus of SLO, whereas detailed analyses of HDXMS at multiple temperatures have uncovered a long-range interaction between the remote PLAT domain and the substrate binding portal of SLO. We believe that these HDXMS results provide the first robust structural evidence for the interaction of lipids at the N terminus of an intact SLO, with a resulting communication of dynamical effects to the catalytic domain.

SLO protein expression and purification
SLO was expressed using the T7-7 plasmid in BL21(DE3) Codon Plus RIL cells and purified as described previously (58). Purity was assessed by SDS-PAGE to be 90%. LC-MS (electrospray ionization) observed: 94,411 Da, calculated: 94,411. Iron analysis using the ferrozine assay resulted in 0.9 Ϯ 0.05 Fe/SLO. The quality of the enzyme preparation was also confirmed by kinetic assays. Before freezing purified aliquots, the protein was dialyzed in 0.1 M borate, pH 9, and concentrated to 100 M. Steady-state kinetic measurements with D 2 -LA were performed (as described previously (58)) in triplicate with 0.1 M borate, pH 9.0, buffer in the presence and absence of effector, OS.

HDX sample preparation
Purified aliquots of wildtype SLO were thawed and diluted to 50 M. Protein aliquots (5 l) from this stock sample were  is colored pale blue. Peptides, whose backbone HDX properties are impacted by interactions with OS, are colored as follows. Yellow, increased extent of exchange (not rate); salmon, increased rate of exchange. A, proposed communication network for OS-mediated allostery includes the conserved cation-interaction (green sticks). The iron cofactor is represented as a dark gray sphere. The modeled substrate is depicted in spheres (orange, carbon; white, hydrogen). Side chains Thr-259 and Leu-541 at the substrate entrance are displayed in spheres for reference. The green highlighted area suggests the binding site for OS. B, a focused view of the cation-interaction between Trp-130 and Arg-242. The conserved carboxylate, Asp-243, also forms hydrogen bonds to Trp-130 backbone.

Liquid chromatography-mass spectrometry for hydrogen-deuterium exchange measurements
Peptide assignments (from the digestion in water) were previously determined for WT SLO (28). Deuterated, pepsin-digested samples of SLO from hydrogen-deuterium exchange experiments were analyzed using a 1200 series liquid chromatograph (LC; Agilent) that was connected in-line with an LTQ Orbitrap XL mass spectrometer (Thermo). The LC was equipped with a reversed-phase analytical column (Viva C8, 30-mm length ϫ 1.0-mm inner diameter, 5-m particles; Restek) and guard precolumn (C8, Restek). Solvent A was 99.9% water, 0.1% formic acid, and solvent B was 99.9% acetonitrile, 0.1% formic acid (v/v). Each sample was thawed immediately before injection onto the column. The elution program consisted of a linear gradient from 5 to 10% B over 1 min, a linear gradient to 40% B over 5 min, a linear gradient to 100% B over 4 min, isocratic conditions at 100% B for 3 min, a linear gradient to 5% B over 0.5 min, and isocratic conditions at 5% B for 5.5 min, at a flow rate of 300 l/min. The column compartment was maintained at 4°C and lined with towels to absorb atmospheric moisture condensation. The column exit was connected to the electrospray ionization source of the mass spectrometer using PEEK tubing (0.005-inch inner diameter ϫ 1 ⁄16-inch outer diameter; Agilent). Mass spectra were acquired in the positive ion mode over the range m/z ϭ 350 -1800 using the Orbitrap mass analyzer, in profile format, with a mass resolution setting of 100,000 (at m/z ϭ 400). Data acquisition was controlled using Xcalibur software (version 2.0.7, Thermo).

HDX data analysis
Mass spectral data acquired for HDX measurements were analyzed using the software HDX WorkBench (62). The percentage of deuterium incorporation was calculated for each of these peptides, taking into account the number of amide linkages (excluding proline residues) and the calculated number of deuterons incorporated. The values were normalized for 100% D 2 O and corrected for peptide-specific back-exchange, as determined previously (28). The data were plotted as deute-rium exchange versus time using Igor Pro software. The rates and extents of exchange (at 4 h) were determined from oneor two-exponential fits to the analyzed time-resolved HDX data.