Pyranopterin Coordination Controls Molybdenum Electrochemistry in Escherichia coli Nitrate Reductase*

Background: The role of the pyranopterin component of the mononuclear molybdenum cofactor is largely unknown. Results: Variants of pyranopterin-coordinating amino acid residues were generated, and their effects on electrochemistry/catalysis investigated. Conclusion: The pyranopterin environment modulates molybdenum electrochemistry. Significance: The pyranopterin coordination environment enables redox-tuning of the molybdenum atom, and facilitates molybdoenzyme reactivity toward a broad range of substrates. We test the hypothesis that pyranopterin (PPT) coordination plays a critical role in defining molybdenum active site redox chemistry and reactivity in the mononuclear molybdoenzymes. The molybdenum atom of Escherichia coli nitrate reductase A (NarGHI) is coordinated by two PPT-dithiolene chelates that are defined as proximal and distal based on their proximity to a [4Fe-4S] cluster known as FS0. We examined variants of two sets of residues involved in PPT coordination: (i) those interacting directly or indirectly with the pyran oxygen of the bicyclic distal PPT (NarG-Ser719, NarG-His1163, and NarG-His1184); and (ii) those involved in bridging the two PPTs and stabilizing the oxidation state of the proximal PPT (NarG-His1092 and NarG-His1098). A S719A variant has essentially no effect on the overall Mo(VI/IV) reduction potential, whereas the H1163A and H1184A variants elicit large effects (ΔEm values of −88 and −36 mV, respectively). Ala variants of His1092 and His1098 also elicit large ΔEm values of −143 and −101 mV, respectively. An Arg variant of His1092 elicits a small ΔEm of +18 mV on the Mo(VI/IV) reduction potential. There is a linear correlation between the molybdenum Em value and both enzyme activity and the ability to support anaerobic respiratory growth on nitrate. These data support a non-innocent role for the PPT moieties in controlling active site metal redox chemistry and catalysis.

rate, trimethylamine-N-oxide, and dimethyl sulfoxide (1, 6 -8). A plausible hypothesis explaining the presence of the second pyranopterin in the DMSOR family is that it adds additional facets of fine tuning to facilitate reactivity toward the observed broad range of substrates. The pyranopterin coordination environments may play a critical role in modulating dithiolene chelate electrostatics via differences in hydrogen bonding environment, conformation, and oxidation state (8, 12-15, 20, 21). Although much attention has been paid to the role of the immediate metal environment in defining catalysis (10,(22)(23)(24), little is known about the role of the pyranopterin coordination environment.
The molybdenum cofactor found in the DMSOR family is a molybdo-bis(pyranopterin guanine dinucleotide) (Mo-bisPGD) (1, 6 -8). The pyranopterins are referred to as proximal and distal based on their positions relative to a conserved [4Fe-4S] cluster, known as FS0, which is present in almost every member of the DMSOR enzyme family. In enzymes for which high-resolution structures are available, the proximal pyranopterin has a more distorted conformation than that of the distal pyranopterin, consistent with the former being in a fully reduced tetrahydro-form; and with the latter having accessibility to the partially oxidized 10,10a-dihydro form (15). In two members of the family, Escherichia coli respiratory nitrate reductase (NarGHI) (25) and Aromatoleum aromaticum ethylbenzene dehydrogenase (EbdABC) (26), high-resolution protein structures reveal that the distal pyranopterin has a bicyclic structure with an open pyran ring. 3 These observations raise the question of the role of the alternate open pyran ring structure in defining active site redox chemistry and catalysis. In the case of NarGHI, protein crystallography reveals that the open pyran ring oxygen participates in hydrogen bonding interactions with two conserved residues in NarG: NarG-His 1163 and NarG-Ser 719 (3) (Fig. 1A). The distances between the Ser 719 OG proton and the pyran oxygen and between the His 1163 NE2 proton and the pyran oxygen are both ϳ2 Å. Closer inspection reveals the presence of a second conserved His residue (NarG-His 1184 ) that completes a charge-transfer relay connecting the pyran oxygen with three structurally conserved water molecules at its distal end (Fig. 1). It would therefore be of interest to generate variants of NarG-Ser 719 , -His 1163 , and -His 1184 and determine their effects on molybdenum electrochemistry and catalysis.
We recently examined the coordination environment of the pyranopterin piperazine nitrogens in members of the DMSOR family for which high-resolution structures are available (8,15), revealing that all members of the family contain a His or Arg residue that bridges the two piperazine N-5 atoms. With the exception of the Thermus thermophilus polysulfide reductase (27), all DMSOR-fold enzyme structures obtained to date contain an additional residue (a His/Ser/Gln) that stabilizes the N-5 of the proximal pyranopterin in an sp 3 hybridized state. Fig. 2A illustrates piperazine nitrogen coordination of the NarGHI pyranopterins. The N-10 atoms of both piperazines are hydrogenbond donors to backbone amide oxygens (of NarG-Thr 259 and NarG-Ser 720 ). The N-5 atoms are bridged by the imidazole of NarG-His 1092 , with the N-5 atom of the proximal pyranopterin participating in an additional H-bond with the NE2 nitrogen of NarG-His 1098 . These two residues are referred to as the "bridging" and "stabilizing" histidines, respectively, and their H-bonding contacts are summarized in Fig. 2B. That almost all of the proximal N-5 atoms in the DMSOR family have two hydrogen bonding contacts favors a tetrahydro oxidation state for the proximal pyranopterin. The single H-bonding contact observed for the N-5 atoms of all the distal pyranopterins provides further support for the suggestion that these have accessibility to the 10,10a-dihydro oxidation state (8,15). In some members of the DMSOR family (e.g. E. coli formate dehydrogenase N (FdnGHI (28)) and the periplasmic nitrate reductase from Cupriavidus necator (NapAB (29))), the bridging ligand is an Arg residue. It would therefore be of interest to generate variants of His 1092 and His 1098 and determine their effects on molybdenum electrochemistry and catalysis.
In this paper, we tested the hypothesis that the pyranopterins of the DMSOR family of molybdoenzymes have a non-innocent role in defining molybdenum electrochemistry and catalysis. If the role of the pyranopterin dithiolene chelates is merely to contribute to a molybdenum-binding scaffold at the active site, then variants of the target residues studied herein should have little or no effect on enzyme function. If, however, the noninnocence reported for model compounds extends beyond the dithiolene chelate and into the heterocyclic ring system of the pyanopterin (13)(14)(15)20), variants of pyranopterin-coordinating residues should have a dramatic effect on molybdenum electrochemistry and catalysis.
Site-directed Mutagenesis-Mutant plasmids were generated using the QuikChange site-directed mutagenesis kit from Stratagene. DpnI was purchased from Invitrogen, and DNA purification kits were purchased from Qiagen. Mutants were verified by DNA sequencing (DNA Core Facility, Department of Biochemistry, University of Alberta). Preparation of competent cells and their transformations with plasmids were carried out as previously described (30).
Generation of a NarG-S719A Variant-A 4.5-kbp SacI-AvaI fragment of pVA700 was subcloned into pTZ18R and the resultant plasmid was used for site-directed mutagenesis. A 3.3kbp SacI-NcoI fragment was cloned back into the pVA700 expression vector, creating pVA700/S719A.
Growth of Cells-Cells were grown microaerobically in 2-liter batch cultures of Terrific Broth (30) at 30°C in the presence of 100 g ml Ϫ1 streptomycin and ampicillin. A 10% innoculum of a stationary phase culture was used and NarGHI overexpression was induced by addition of 0.2 mM isopropyl 1-thio-␤-Dgalactopyranoside. Following addition of the innoculum and isopropyl 1-thio-␤-D-galactopyranoside, cells were grown overnight with gentle shaking, and were harvested by centrifugation, and subsequently washed in a buffer containing 100 mM MOPS and 5 mM EDTA (pH 7.0). Cells were resuspended in buffer, and phenylmethylsulfonyl fluoride was added to a final concentration of 0.2 mM. Cell lysis was achieved by three passages through an Emulsiflex C3 microfluidizer (Avestin) at a pressure of 17,000 p.s.i. Crude membranes were prepared by differential centrifugation as previously described (31). These were resuspended in buffer and layered on top of a 55% (w/v) sucrose step (made up in buffer). Following ultracentrifugation at 40,000 rpm for 1.5 h, the floating band enriched in cytoplasmic membrane vesicles was collected. Excess sucrose was removed by two dilution and ultracentrifugation steps. The final membrane pellet was resuspended in buffer to a concentration of ϳ30 mg ml Ϫ1 , and flash frozen in liquid nitrogen prior to being stored at Ϫ70°C until use. Where appropriate, buffer exchange prior to EPR analysis was achieved by dilution and re-centrifugation. Overexpression of NarGHI was evaluated by polyacrylamide gel electrophoresis (32).

Bacterial Growth on Glycerol-Nitrate Minimal Medium-
Anaerobic growth of E. coli harboring LCB79/pVA700 and mutant derivatives was assessed in a glycerol-nitrate minimal medium essentially as previously described (33,34). Nitrate was added (as KNO 3 ) to a final concentration of 40 mM. Growth was evaluated at 37°C, and culture turbidity was measured using a Klett-Summerson spectrophotometer equipped with a number 66 filter. Maximal growth rates were calculated as described by Zwietering et al. (35), and were expressed as the parameter m in units of (Klett units) h Ϫ1 .
Form A Fluorescence of the Pyranopterin Cofactor-Pyranopterin was assayed by generation of its Form A fluorescent derivative following protein denaturation (36,37), essentially as previously described (33). As starting material, 10 mg of membrane protein was used.
Redox Potentiometry and EPR Spectroscopy-Redox titrations were carried out under argon at 25°C as previously described (38,39) in 100 mM Tricine and 5 mM EDTA (pH 8.0). The protein concentration used was ϳ30 mg ml Ϫ1 . The following redox mediators were used at a concentration of 25 M: quinhydrone, 2,6-dichloroindophenol, 1,2-naphthoquinone, toluene blue, phenazine methosulfate, thionine, methylene blue, resorufin, indigotrisulfonate, indigocarmine, anthraquinone-2-sulfonic acid, phenosafranine, and neutral red. All samples were prepared in quartz EPR tubes with a 3-mm internal diameter, and were rapidly frozen in liquid nitrogen-chilled ethanol and stored at 77 K until use. EPR spectra were recorded at 150 K using a Bruker Elexsys E500 spectrometer equipped FIGURE 1. A, charge-transfer relay in NarG. Three water molecules are conserved across a range of NarGHI structures, including those described by Protein Data Bank codes 1Q16, 1R27, 1YZ4 and 1Y5I (25,60,61). These waters are linked to the pyran oxygen of the bicyclic distal pyranopterin by NarG-His 1184 and NarG-His 1163 . The OG oxygen of NarG-Ser 719 is also within H-bonding distance of the pyran oxygen. The image was generated using the PyMol molecular visualization package (version 1.4.1, Schrödinger LLC). B, mechanism of charge-transfer between the pyran oxygen and the conserved water molecules. The indicated charge-transfer results in protonation of the alkoxide anion form of the pyran oxygen. The image was generated using the MarvinSketch software package (www.chemaxon.com). In both panels, some predicted hydrogens have been added for clarity.
with a Bruker SHQE microwave cavity and a Bruker ER4131 Variable Temperature Unit using liquid nitrogen as a cryogen. EPR spectra were plotted as signal intensity versus g value, with the latter being calculated from the microwave frequency and the field intensity for each sample studied. This allowed for direct comparison of the g values of the individual spectral features and corrected for minor variations in microwave frequency (Fig. 3). Potentiometric titration data were analyzed by plotting the intensity of the molybdenum peak-trough versus E h and fitting the data to two E m values (E 1 and E 2 ) (40). Mo(V) stability constants (K stab ) were calculated as described by Ohnishi et al. (41). Fitting was performed using the Solver for Nonlinear Programming extension of Libre Office Calc. Where appropriate the statistical significance of differences in E m values was evaluated using a t test.
Protein Assays-Protein concentrations were assayed as described by Lowry et al. (42) modified by the inclusion of 1% (w/v) sodium dodecyl sulfate in the incubation mixture to solubilize membrane proteins (43).
Enzyme Assays-Nitrate reductase activities were measured using the nonspecific electron donor benzyl viologen (44). Assays were carried out at pH 8.0 in a buffer containing 100 mM Tricine, 0.23 mM benzyl viologen, 5 mM EDTA, and 4 mM KNO 3 . Benzyl viologen was reduced by adding an excess of sodium dithionite (ϳ0.4 mM), and the reaction was initiated by addition of a catalytic amount of enzyme.

Results and Discussion
Residues Targeted for Site-directed Mutagenesis-During nitrate reduction, the molybdenum atom of the Mo-bisPGD cycles through the IV, V, and VI oxidation states to catalyze the two-electron reduction of nitrate to nitrite (45). To investigate the effect of pyranopterin coordination on molybdenum electrochemistry, we generated variants of two sets of pyranopterin-coordinating residues in NarG: (i) those interacting directly or indirectly with the oxygen of the open pyran ring of the distal pyranopterin (NarG-Ser 719 , NarG-His 1163 , and NarG-His 1184 ; Fig. 1); and (ii) those involved in bridging the proximal and distal pyranopterins and stabilizing the oxidation state of the proximal pyranopterin (NarG-His 1092 and NarG-His 1098 ; Fig. 2). As judged by SDS-PAGE and Form A pyranopterin assays, each variant was expressed in a Mo-bisPGD-containing form at levels comparable with that of the wild-type. We used these variants to test the hypothesis that residues coordinating the pyranopterin moiety of the Mo-bisPGD cofactor play a critical role in defining molybdenum electrochemistry and substrate reactivity.
Impact of the Variants on the NarG Mo(V) EPR Spectrum-Of the three accessible molybdenum oxidation states, only the intermediate Mo(V) state is paramagnetic and EPR visible. It can be observed in spectra recorded of samples poised at appropriate reduction potentials. Because the NarGHI Mo(V) signal is sensitive to pH in an anion-dependent way (46,47), we recorded spectra of samples poised at pH 8.0, wherein the socalled "high pH" form is dominant. This form exhibits a simple rhombic EPR spectrum with an absence of resolved hyperfine ( 1 H) couplings (46), simplifying subsequent analyses of potentiometric data (see below). Fig. 3A shows spectra of redox-poised samples containing the wild-type enzyme and variants of residues implicated in coordinating the open pyran oxygen of the distal pyranopterin (S719A, H1163A, and H1184A). These spectra are similar to those previously observed for the Mo(V) species of NarGHI (46,47), wherein a rhombic Mo(V) species is observed with g 1,2,3 values of ϳ1.988, 1.982, and 1.962. As was the case in previous studies, features consistent with the presence of a small amount of the "low pH" form of the Mo(V) are observed immediately downfield from the high-pH g 1 (46,47), but these are inconsistently observed in the variants studied herein. Also in agreement with previous work (46 -48), no hyperfine coupling was resolved with the isotopes of molybdenum having a nuclear spin of 5/2 ( 95 Mo and 97 Mo, with natural abundances of ϳ16 and 10%, respectively). These data indicate that the S719A, H1163A, and H1184A variants do not significantly impact the molybdenum coordination sphere. A, residues defining pyranopterin piperazine ring coordination. The image was generated using the NarGHI structure described by Protein Data Bank code 1Q16 (25). His 1092 functions to bridge the two pyranopterins, with its ND1 nitrogen functioning as an H-bond donor to the proximal piperazine N-5 atom, and its NE2 nitrogen functioning as an H-bond donor to the distal piperazine N-5 atom. Note that His 1092 is shown in its protonated form to facilitate H-bond donation to both piperazine N-5 atoms. For further details, see text. B, proposed H-bonding network around the piperazine rings of the two pyranopterins. The proximal pyranopterin (labeled P) is shown in its tetrahydro form, with this redox state stabilized by both the bridging His 1092 and the stabilizing His 1098 , whereas the piperazine of the distal pyranopterin (labeled D) is shown in a form equivalent to the 10,10a-dihydro pyranopterin. In both panels, some predicted hydrogens have been added for clarity. Images were generated as described in the legend to Fig. 1. Fig. 3B shows spectra of redox-poised samples of variants of residues implicated in bridging the two pyranopterins (H1092A and H1092R), or stabilizing the proximal pyranopterin in its tetrahydro form (H1098A). Of these, the H1092A variant has the most dramatic effect on the Mo(V) EPR spectrum, shifting g 2 upfield from approximately g ϭ 1.982 (wild-type) to approximately g ϭ 1.962 (H1092A). The H1092R variant exhibits an EPR spectrum wherein heterogeneity is observed, especially in the g 3 feature. Finally, the H1098A variant exhibits a spectrum with clearly defined g 1,2,3 values, with greater rhombicity than that of the wild-type. Overall, Mo(V) species are clearly visible in all of the NarGHI variants studied herein, with the H1092A variant having the greatest effect on the Mo(V) EPR spectrum (see below).
Influence of the NarGHI Variants on Molybdenum Electrochemistry-Molybdenum redox cycling is defined by two reduction potentials, E 1 and E 2 , corresponding to the Mo(VI/V) and Mo(V/IV) reduction potentials, with the average of these being defined as the overall reduction potential (E m ). The difference between E 1 and E 2 is a reflection of the stability of the intermediate Mo(V) species, and can be used to calculate its stability constant (K stab ) (41). We performed potentiometric titrations on all the variants of NarGHI studied herein, and in each case titrations were carried out on samples generated from 3 to 4 independent biological replicates. Fig. 4A shows representative titrations for the variants of residues interacting with the O-1 pyran oxygen (Fig. 1). Fig. 4B shows plots of the residuals between the experimental data and the fits. For the wild-type enzyme, we estimate values of E m and K stab of 142 mV and 28, respectively (summarized in Fig. 4C). These values are in reasonable agreement with those previously obtained (38, 49 -51). The S719A variant has a statistically insignificant effect on the E m (⌬E m ϭ 4 mV, p ϭ 0.57) and decreases the K stab from 28 to 8, whereas the H1163A variant elicits a large decrease in the E m (⌬E m ϭ Ϫ88 mV, Table 1) and a concomitant decrease in the K stab from 28 to 1. The H1184A variant has an intermediate effect, eliciting a ⌬E m of Ϫ36 mV, and decreases the K stab , again from 28 to 1. Thus, the charge-transfer relay depicted in Fig. 1B appears to play two critical roles: (i) modulating the molybdenum reduction potentials (see below for further discussion of this); and (ii) stabilizing the Mo(V) intermediate. Fig. 5A shows representative potentiometric titrations of the H1092A, H1092R, and H1098A variants. The H1092R variant replaces the bridging His residue with an Arg, and mimics the pyranopterin-coordination environment observed in the formate dehydrogenases and periplasmic nitrate reductases (e.g. FdnGHI and NapAB, respectively (8)). The H1092R variant elicits a modest increase of marginal statistical significance in the overall molybdenum E m (⌬E m ϭ 18 mV, p ϭ 0.07; Table 1), consistent with there being no change in overall charge of the bridging residue when it is changed from a His to an Arg (the pK a of the Arg guanidinium group is ϳ12.5). This suggests that the imidazole pK a is significantly increased from ϳ6.0 in aqueous solution to Ͼ8.0 in the bridging environment between the two pyranopterins. The H1092A variant elicits a ⌬E m on the molybdenum center of approximately Ϫ143 mV along with a  modest decrease in its K stab (from 28 to 19, Table 1). A possible explanation for the large negative ⌬E m is that the void between the two pyranopterins contains waters in the H1092A variant, and that these waters lack an overall positive charge, thus decreasing the overall molybdenum reduction potential. Alternatively, given that the H1092A variant has the most severe effect on the Mo(V) EPR spectrum, loss of "scaffolding" by the bridging imidazole moiety may result in a significant alteration of the geometry of dithiolene coordination that may cause the large observed ⌬E m .
The H1098A variant elicits a ⌬E m of Ϫ101 mV, along with a significant increase in Mo(V) stability, with the K stab increasing from 28 to 1822 (Fig. 5, A and C, Table 1). Fig. 2B shows a model for how His 1092 and His 1098 define the H-bonding environment of the proximal pyranopterin piperazine N-5 atom. As suggested above, formal protonation of the imidazole of His 1092 is supported by the relative lack of effect of the H1092R variant on the overall molybdenum E m . One of the H-bonds is eliminated in the H1098A variant, which may result in transfer of the proton from the ND1 atom to the piperazine N-5, resulting in loss of a positive charge and the observed decrease in the overall molybdenum E m .
We proposed that in the DMSOR family of enzymes, the role of the stabilizing residue (His 1098 in NarG) is to maintain the N-5 atom of the proximal pyranopterin in an sp 3 -hybridized form, thus stabilizing the tetrahydro oxidation state (8,15). This suggests a possible explanation for the observed increase in Mo(V) stability in the H1098A variant. Removal of the hydrogen-bonding contact provided by the His 1098 imidazole may result in the 10,10a-dihydro oxidation state becoming available to the proximal pyranopterin. If the 10,10a-dihydro to tetrahydro reduction potential occurs within the range in which Mo(V) is extant, then E 1 (the Mo(VI/V) reduction potential) could be shifted to the observed value of 139 mV (Table 1; see Ref. 15).
Correlation between Enzyme Activity, Cell Growth, and Molybdenum Electrochemistry-How do changes in molybdenum electrochemistry impact enzyme activity and anaerobic respiratory growth on nitrate? Specific activities in enriched membranes for the variants studied herein range from 1 to 54 mol min Ϫ1 mg Ϫ1 for the H1092A and S719A variants, respectively (Table 1). With the exception of the H1092R variant, a strong correlation exists between overall molybdenum reduction potential and enzyme activity (r ϭ 0.99, Fig. 6). A correlation also exists between respiratory growth rate and molybdenum reduction potential across all variants (r ϭ 0.95, inset to Fig. 6). Thus, catalysis and growth are convincingly correlated with increased molybdenum E m . A possible explanation for this is that catalytic turnover depends on the energetics of electron transfer from FS0 (E m ϭ Ϫ55 mV (49)) to the molybdenum.
It is notable that the NarGHI catalytic activity appears to be dependent on the overall molybdenum E m , but not on the Mo(V) K stab (Table 1). Protein film voltammetry studies of NarGHI (45) and two other members of the DMSOR family (the periplasmic nitrate reductase from Rhodobacter sphaeroides, NapAB (52), and the Me 2 SO reductase from E. coli (53)) demonstrated that enzyme turnover decreases below a reduction potential referred to as E switch . It was therefore proposed that observation of an E switch arises because substrate binds preferentially to the Mo(V) form of the cofactor (45,52,53). In the case of NarGHI, this is corroborated by the observation of a stable nitrate adduct to the Mo(V) (47). If the mechanism of nitrate reduction proceeds via substrate binding to a Mo(V) species, we would expect a variant that increases the Mo(V) K stab to retain significant catalytic activity. The H1098A variant exhibits a very large Mo(V) K stab (ϳ1822; Table 1), but has low specific activity and growth rate compared with the wild-type (Fig. 6).
As previously noted (15), an alternative explanation for the observed E switch is that it controls a reduction of the distal pyranopterin from an oxidation state equivalent to the 10,10a-dihydro form to one equivalent to the tetrahydro form. This provides a possible explanation for why the H1092R variant supports respiratory growth on nitrate, but has decreased enzyme activity in our in vitro assay (Fig. 6). It is possible that a significant positive ⌬E m is elicited on the E switch , resulting in  Table 1. Colored squares indicate the mean E m values and the horizontal lines indicate the E 1 and E 2 values in the cases where E 1 Ϫ E 2 Ͼ 0 (K stab Ն 1).

Pyranopterin Coordination in E. coli Nitrate Reductase A (NarGHI)
OCTOBER 9, 2015 • VOLUME 290 • NUMBER 41 lower activity measured using benzyl viologen as electron donor (E m ϭ Ϫ374 mV (54)). NarGHI is able to function in vivo with either ubiquinol or menaquinol as electron donor (55), and these have reduction potentials of ϩ110 and Ϫ80 mV, respectively (2), rendering them both less likely to reduce the distal pyranopterin than the artificial electron donor benzyl viologen.

Role of the His 1163 /His 1184 Charge-transfer Relay in NarGHI
Maturation and in Modulation of Molybdenum Electrochemistry-NarGHI maturation is a carefully orchestrated process involving the NarJ system-specific chaperone (56) and the twin arginine translocase (57). The process ensures assembly of correctly folded enzyme with a complete complement of iron-sulfur clusters and hemes (56,58). Insertion of the Mo-bisPGD cofactor depends on the presence of the FS0 [4Fe-4S] cluster (56,58). With the exception of EbdABC (26), NarGHI is the only enzyme examined to date that has a bicyclic distal pyranopterin (25). Closer examination of the charge-transfer relay comprising His 1163 and His 1184 (Fig. 1) suggests a possible mechanism for pyran ring opening (Fig. 7), which proceeds as follows. The His 1163 /His 1184 charge-transfer relay catalyzes elimination of the C-4a proton (Fig. 7A), generating an open ring-form with a piperazine ring in a form equivalent to the 5,10-dihydro form (Fig. 7B) in a mechanism essentially identi-  . Correlation between enzyme activity and overall molybdenum reduction potential. Excluding the NarG-H1092R variant, a correlation exists between enzyme activity and overall molybdenum reduction potential (r ϭ 0.99). Inset, correlation between growth rate and overall molybdenum reduction potential. Including the NarG-H1092R variant, a correlation exists between the maximal rate of respiratory growth on nitrate and overall molybdenum reduction potential (r ϭ 0.95). Growth rates were calculated using the method of Zwietering et al. (35). cal to that previously described (59). This form rearranges to the lowest energy tautomer (equivalent to the 10,10a-dihydro form (15,59)) shown in Fig. 7C. A mechanism essentially identical to that presented in Fig. 7, A-C, can be proposed for ringopening of the distal pyranopterin in EbdABC, with the distinction that ring opening would be catalyzed by a guanidinium side chain rather than by an imidazole. The equilibrium between the structures shown in Fig. 7, C and D, illustrates how the formal change on the pyran oxygen (and the overall molybdenum reduction potential) can be modulated by the His 1163 /His 1184 charge-transfer relay. Stabilization of a the pyran oxygen in its alkoxide form would decrease the overall molybdenum reduction potential, whereas formal protonation to the hydroxyl form would increase the overall molybdenum reduction potential.
In the context of the model presented in Fig. 7, C and D, for charge-transfer-induced molybdenum reduction potential modulation, it is notable that the H1163A variant elicits a large ⌬E m that is approximately double that elicited by the H1184A variant. This is consistent with the entire charge-transfer relay being involved in molybdenum reduction potential modulation. It is possible to speculate that the H1163A variant may impact an equilibrium between the open and closed pyran ring forms of the distal pyranopterin, and we are currently addressing this issue using protein crystallography.
Conclusions-We have demonstrated the importance of the pyranopterin coordination environment in defining molybde-num electrochemistry and substrate reactivity, with enzyme activity correlating with increasing overall molybdenum reduction potential. Our studies demonstrate the importance of factors beyond the immediate molybdenum coordination sphere in defining enzyme activity.
Author Contributions-S. W. designed, performed, and analyzed the experiments presented in the manuscript. R. A. R. wrote the paper and contributed to the data analyses. J. H. W. conceived and coordinated the study.