Kinetics of CO Ligation with Nitric-oxide Synthase by Flash Photolysis and Stopped-flow Spectrophotometry*

(Received for publication, January 30, 1997, and in revised form, March 3, 1997)

Jürgen S. Scheele Dagger , Vladimir G. Kharitonov §, Pavel Martásek , Linda J. Roman , Vijay S. Sharma §, Bettie Sue Siler Masters and Douglas Magde Dagger par

From the Dagger  Department of Chemistry and Biochemistry and the § Department of Medicine, University of California at San Diego, La Jolla, California 92093 and the  Department of Biochemistry, University of Texas Health Sciences Center, San Antonio, Texas 78284-7760

ABSTRACT
INTRODUCTION
MATERIALS AND METHODS
RESULTS
DISCUSSION
FOOTNOTES
ACKNOWLEDGEMENTS
REFERENCES


ABSTRACT

Interaction of CO with hemeproteins has physiological importance. This is especially true for nitric-oxide synthases (NOS), heme/flavoenzymes that produce ·NO and citrulline from L-arginine (Arg) and are inhibited by CO in vitro. The kinetics of CO ligation with both neuronal NOS and its heme domain module were determined in the presence and absence of tetrahydrobiopterin and Arg to allow comparison with other hemeproteins. Geminate recombination in the nanosecond time domain is followed by bimolecular association in the millisecond time domain. Complex association kinetics imply considerable heterogeneity but can be approximated with two forms, one fast (2-3 × 106 M-1 s-1) and another slow (2-4 × 104 M-1 s-1). The relative proportions of the two forms vary with conditions. For the heme domain, fast forms dominate except in the presence of both tetrahydrobiopterin and Arg. In the holoenzyme, slow forms dominate except when both reagents are absent. Geminate recombination is substantial, ~50%, only when fast forms predominate. Stopped-flow mixing found dissociation constants near 0.3 s-1. These data imply an equilibrium constant such that very little CO should bind at physiological conditions unless large CO concentrations are present locally.


INTRODUCTION

Nitric-oxide synthase is a hemoprotein that catalyzes conversion of the amino acid L-arginine to citrulline and nitric oxide (1, 2). The NO produced activates another heme enzyme, guanylate cyclase (see Reaction R1), that catalyzes formation of cyclic guanosine monophosphate, a second messenger that mediates numerous biochemical events including vascular smooth muscle relaxation, platelet disaggregation, photoreceptor cell signaling, ion transport in gastrointestinal cells, and myeloid cell differentiation.


[View Larger Version of this Image (7K GIF file)]

Fig. R1.

At least three isoforms of NOS1 have been purified and characterized. In each, a heme prosthetic group, iron protoporphyrin IX, is attached to the protein by a proximal cysteine. Visible absorption spectroscopy and other physicochemical measurements reveal similarities to cytochromes P450 (3-8). Like that family of proteins, reduced (ferro) NOS forms carboxy derivatives with an absorption maximum at ~446 nm (4). Marletta (9) and Stuehr et al. (10) have proposed mechanisms for the oxidation of Arg catalyzed by NOS, based largely upon the similarity of NOS and P450, and Masters (11) has proposed a scheme involving redox interactions with the flavoprotein module of NOS. A key observation was that enzyme activity is inhibited by CO (4, 6, 9) as well as by NO (12), the physiological significance of which remains uncertain. Enzyme inhibition in vitro required large concentrations of either CO or NO. If those reagents play a role in vivo, the mechanism must involve locally enhanced concentrations. For NO, fast turnover and slow diffusion away from the site of production in unstirred cytoplasm might possibly lead to large concentrations locally and, consequently, to autoregulation of NO production. Recent evidence suggests a role for CO produced by heme oxygenase as a messenger molecule in the same cell types as nNOS (13-15). Could the activity of heme oxygenase produce high concentrations? The CO concentration needed depends upon the equilibrium constant for CO binding. Unfortunately, equilibrium constants for binding CO and NO are not known for NOS. Matsuoka et al. (16) reported CO combination rate constants for neuronal NOS but not CO dissociation data, which could have been combined with the combination rate constant to calculate an equilibrium constant. In hemeproteins, CO dissociation rate constants vary from 6.5 s-1 in P450 to 7.2 × 10-5 s-1 in horseradish peroxidase. Therefore, knowledge of CO association and dissociation is important for determining NOS reactivity relative to myoglobin, guanylate cyclase, and other heme proteins, which potentially could compete with NOS for reactive ligands.

In addition to the interest generated by its physiology, this sensitivity of CO dissociation rate constants to protein structure makes it an excellent parameter for providing qualitative information about the local heme environment. A review of CO dissociation constants reveals two trends. 1) CO off rates are markedly affected by the nature of any base proximal to the heme-CO bond, whereas for a given base, any weakening of the base-heme bond, or its total rupture, dramatically increases the CO dissociation rate. 2) Steric hindrance on the distal side reduces both association and dissociation rates for CO in all cases known (17).

These same issues arise for carboxy guanylate cyclase activation. In vitro activity is significantly affected by 1 atm of CO, but equilibrium constants indicate that in vivo, little, if any, guanylate cyclase activation by CO should be expected. Nevertheless, several biological phenomena seem to be affected by small amounts of CO (15). We have no explanation, but the example points out the difficulty in extrapolating from in vitro to in vivo conditions. Nevertheless, a knowledge of equilibrium and rate parameters is the essential starting point of all discussion.

The studies reported here determined the kinetics of ligation of CO with neuronal nitric-oxide synthase stably expressed in human kidney 293 cells (4) and a holoenzyme (18) and heme module (residues 1-714) (19) expressed in Escherichia coli. Flash photolysis and stopped-flow spectrophotometry were used to determine association and dissociation rate constants. We investigated the effects of tetrahydrobiopterin (BH4) and L-arginine (Arg). For the latter, advantage was taken of the ability in E. coli to express nitric-oxide synthases in the absence of BH4 and with diminished Arg concentrations (18, 20) in the case of the heme domain. In addition to bimolecular association, we also investigated geminate recombination of CO-iron pairs confined within the protein. These studies allowed us to estimate the equilibrium binding constant of CO with NOS and provided insights into the nature of the heme pocket in the enzyme.


MATERIALS AND METHODS

Carbon monoxide (99.8%), nitric oxide (99.9%), high purity argon, and premixed CO in argon were all from Matheson. NO was further purified by passage through a column of fresh KOH. Buffers were prepared (20 mM Tris-HCl, usually pH 7.8, 100 µM EDTA, 100 mM NaCl) in a gas-tight syringe and deoxygenated by bubbling with argon for >40 min. When indicated, buffers also contained 250 µM BH4 and/or 0.5 or 10 mM Arg. Microperoxidase solutions were prepared as described elsewhere (21).

Rat nNOS holoenzyme was expressed in two different systems. Unless otherwise indicated, it was in E. coli, purified and reconstituted as described previously (18). Nitric oxide formation, measured in E. coli-expressed nNOS after Arg and BH4 repletion using a hemoglobin capture assay described previously (22), was found to be 250 nmol min-1 mg-1 (25 °C) or 540 nmol min-1 mg-1 (25 °C). Expression of rat nNOS in human embryonic kidney cell was performed as reported previously (4); nitric oxide formation was measured using an L-[3H]citrulline assay (23) and found to be 420 nmol min-1 mg-1 (25 °C). The amino-terminal heme-binding domain (residues 1-714) of rat nNOS was expressed in E. coli and purified as described previously (19). All proteins were >90% homogeneous as determined by sodium dodecyl sulfate-polyacrylamide gel electrophoresis. Protein solutions were degassed and then reduced using degassed 1% sodium dithionite solution. Finally, CO was added in a premixed ratio with argon by gentle gas flow over the liquid using a large surface-to-volume ratio. Spectra for reduced holoenzyme and heme domain were consistent with previously published data as were the spectra of liganded preparations (7). Spectra were taken before and after flash photolysis measurements to detect any possible degradation.

Flash photolysis used an improved version of an instrument described previously (24). Photolysis laser pulses were 4 ns in duration at 545 nm with up to 4 mJ of energy over an area of 0.1 cm2. The probe had an 8-nm bandwidth selected from a stable tungsten lamp, except that measurements of geminate recombination required a pulsed xenon flash. The photomultiplier was wired with both a large standing current and large interdynode capacitors to optimize linearity. The digitizer had been upgraded to a Lecroy Model 9361 digitizing oscilloscope.

A quantitative treatment of signal-to-noise ratio (S/N) is essential whenever there is a question of complicated kinetics. Heme concentrations were near 1 µM with 1-cm optical path length using a volume of 0.5 cm3 in a 1-cm square cell. (This gave a large surface area for gas exchange after we learned that even mild stirring caused degradation.) Initial absorbance changes Delta A (t = 0) were typically about 0.03-0.04. By accumulating 80-200 recordings, the root mean square uncertainty at each time point was reduced below 0.00008. As a result, we achieved S/N approx  500 or better. Nanosecond geminate recombination had to be measured with a wider electronic bandwidth, resulting in about a 5-fold loss in S/N despite the brighter, pulsed lamp. Good S/N for random scatter exposes a measurement to systematic distortions in the data. These have been reduced through several generations of the instrument. The most convincing short argument that systematic distortion is not responsible for any major feature described below is that, by coincidence, while these measurements were being made we were also measuring kinetics of recombination in HbA2. The latter system gave beautiful single exponential decay curves.

Stopped-flow measurements were carried out at 20 °C using a Durrum instrument following protocols described earlier (25). The CO dissociation rates were determined both by the microperoxidase method (21) and by replacing CO with NO (26). For the microperoxidase method, after mixing, the sample was 0.5 µM protein, 200 µM sodium dithionite, and 5 µM CO. The microperoxidase concentrations varied between 3 and 10 µM. For the NO substitution method, the sample was 0.5 µM protein, 25 µM sodium dithionite, and 150 µM CO. The NO concentrations varied between 175 and 275 µM. Spectra before and after reactions were as expected (27). For both methods, absorption changes were monitored at 443 nm. Measurements were done in duplicate and were proved to be independent of concentration of either microperoxidase or NO.


RESULTS

The kinetics of the reaction of CO with both the complete nNOS holoenzyme and the heme domain fraction of nNOS were characterized, and the rates in the presence and absence of Arg and BH4 were measured. Thus, there are eight cases to distinguish. The effect of Ca2+/calmodulin on rates of CO binding has been previously shown to be negligible (16).

Association

Flash photolysis measurements were made at 23 ± 0.5 °C. They were carried out at both pH 7.0 and pH 7.8, with identical results. There was also no discernible difference between nNOS from kidney cells and that from E. coli. There was no irreversible photochemistry. Transient recovery was >99.5% after each flash, and spectra did not change during a series of flashes. Association rates were determined by monitoring recombination after photodissociation in the presence of excess CO. Kinetics were quite complex. In general, we observed three phases on different time scales. There was a very fast process with a half-life on the order of 100 ns that was unaffected by changing CO concentration and was assigned as geminate recombination of protein-caged heme-ligand pairs. All slower processes varied systematically with [CO]. Solutions were equilibrated with premixed gas mixtures at 1 atm (105 pascals) total pressure. We used CO in argon at molar ratios of 5, 10, 20, and 50% as well as 100% CO. Molar concentrations were calculated assuming that the solubility of CO in aqueous solution at 25 °C is 930 µM atm-1.

We must distinguish a kinetic "phase," as the term was used above, from an exponential component of a fit. The two words are frequently used as synonyms. Here, however, a phase may be composed of two or even three exponentials that are quite close in characteristic time constant. The group as a whole differs from other rates by at least a factor of 10, and the components behave similarly as CO concentration is changed.

Best fits to the data were calculated using our own implementation of the classic Marquardt algorithm for a nonlinear least squares fit to a sum of exponentials,
&Dgr;A(t)=s <LIM><OP>∑</OP></LIM> a<SUB>i</SUB>e<SUP><UP>−</UP>k<SUB>i</SUB>t</SUP> (Eq. 1)
where s is an overall scale factor, ai are amplitudes, which sum to unity, and ki is the rate constant. The nanosecond geminate phase involved only a single exponential, characterized by agem and kgem. The bimolecular association required as many as four additional components, but they could be grouped into fast and slow phases characterized by kf or ks, each defined as a sum over one, two, or, occasionally, three exponentials as follows.
k<SUP><UP>−</UP>1</SUP><SUB>f</SUB>=<LIM><OP>∑</OP></LIM> a<SUB>j</SUB>k<SUP><UP>−</UP>1</SUP><SUB>j</SUB>/<LIM><OP>∑</OP></LIM> a<SUB>j</SUB> (Eq. 2)
Fig. 1 displays a typical transient absorption measurement with a multiexponential fit. The definition of the mean rate as the mean of reciprocal rates (see Ref. 2) mimics what a single exponential fit to the same data yields.


Fig. 1. Top panel, semilog display of transient absorbance change at 445 nm after flash photolysis of holoenzyme in the presence of tetrahydrobiopterin and arginine and 20% atm CO, 80% atm argon. a, the best single exponential fit (Delta A = 0.0539 e-5.11t) might be thought to fit several half-lives, but it clearly deviates at long times and, although difficult to see in the semilog plot, also at short times. b, a triple exponential fit of the major slow phase plus a small single exponential for a fast phase (Delta A(t) = 0.0566 (0.0177 e-763t + 0.0864 e-14.8t + 0.5883 e-5.75t + 0.3074 e-3.74t)) appears coincident with the data. The weighted mean rate of the dominant, slower phase is 5.17, only 1% different from the best single exponential. Bottom panel, a linear display of the residuals between fits and data shows clearly that the single exponential (a) fails miserably (reduced chi 2 = 18 for root mean square noise delta Delta A = 9.1 × 10-5), whereas the triple exponential (b) fits well (reduced chi 2 = 1.1).
[View Larger Version of this Image (18K GIF file)]

To determine the bimolecular association constants ka,f or ka,s, measurements were carried out for several [CO] to prove the linear dependence on [CO]. The phases were distinguished according to Equations 1 and 2, and the rates ka,f and ka,s were plotted versus [CO]. The association rate constants are the slopes of those graphs.

There are three reasons for simplifying to two bimolecular phases. First, it permits easier comparison with other experiments that do not have the S/N to resolve all components. Second, fitting exponentials with similar rates is a difficult statistical problem. Parameters are highly correlated, and even though heterogeneity is clear, little significance should be attached to a particular expression of that heterogeneity. It is probable that there is a continuous distribution of rates. Third, since so little is known at present about NOS, it is preferable to simplify the discussion for the present while documenting whatever detail is measured for its potential future significance.

To confirm that the complex kinetics were not an instrumental artifact, we carried out a number of control experiments. 1) We investigated pulse energies from 0.4 to 4 mJ to confirm that there was no effect on the measured decay curves. 2) In preliminary measurements, we recorded kinetic traces at 5-nm intervals from 380 to 450 nm and obtained a transient spectrum very similar to that reported previously (16). In particular, we confirmed the extremes of the difference spectrum, the location of the isosbestic point, and the near equality of the magnitude of peak transient absorption and bleach. 3) We confirmed that we had a well defined isosbestic point at 430-432 nm (slightly shifted between the heme domain and the holoenzyme), which persisted throughout the entire time range. Subsequently, in each measurement traces were collected both at the maximum transient bleach, 445 nm, and the maximum transient absorption, 405 nm. Kinetics reported are averages of the two. The heme domain decay curves were often indistinguishable at the two wavelengths. There may be a small difference in some cases for the holoprotein, presumably due to very small, time-dependent shifts in the transient spectrum related to conformational relaxation of the holoprotein.

Principal data for all measurements are collected in Table I, but the following sections add information about certain details, mainly the multiple exponentials included in certain phases.

Table I. Kinetic rate constants for binding CO to nNOS heme domain and nNOS holoenzyme

Rates listed here do not always represent single exponential decays. Rates listed here do not always represent single exponential decays.
Systema agem/Delta A(0) kgem af/(af + as) ka,f ka,s kd

s-1 M-1s-1 M-1s-1 s-1
Dom 0.35  ± 0.10 2.0  ± 0.5 × 107 0.95  ± 0.03 2.1  ± 0.2 × 106 3.5  ± 0.5 × 104 0.15
Dom + Arg 0.40  ± 0.10 1.6  ± 0.5 × 107 0.90  ± 0.04 3.5  ± 0.2 × 106 2.2  ± 1.5 × 104 0.32
Dom + BH4 0.55  ± 0.10 1.8  ± 0.3 × 107 0.95  ± 0.02 3.0  ± 0.4 × 106 3.0  ± 0.3 × 104 0.16
Dom + Arg + BH4 0.08  ± 0.03 1.8  ± 0.5 × 107 0.30  ± 0.05 1.5  ± 0.1 × 106 2.1  ± 0.1 × 104 0.34
Holo 0.45  ± 0.05 1.3  ± 0.5 × 107 0.94  ± 0.02 2.8  ± 0.2 × 106 4.2  ± 0.5 × 104
Holo + Arg 0.09  ± 0.03 1.2  ± 0.4 × 107 0.10  ± 0.03 2.0  ± 0.2 × 106 1.8  ± 0.2 × 104
Holo + BH4 0.10  ± 0.03 1.8  ± 0.5 × 107 0.14  ± 0.03 1.9  ± 0.2 × 106 2.2  ± 0.2 × 104 0.16 and 15
Holo + Arg + BH4 <0.03 0.03  ± 0.01 ~2-3 × 106 2.0  ± 0.1 × 104 0.33

a Dom, domain; Holo, holoenzyme.

Heme Domain with neither L-Arginine nor BH4

The fast phase could be fit fairly well using two exponentials of equal amplitude with rates differing by a factor of 3, but at very high S/N it became clear that two exponentials were not completely adequate. The steady-state absorption spectrum showed a small amount of extra absorption near 420 nm, consistent with that reported and explained by Wang et al. (28). We also noticed additional complexity in the transient spectra, particularly around the usual isosbestic point. This could be pursued, if one were interested in the so-called P420 form of NOS, but for our present purpose we believe that by choosing wavelengths characteristic of the P450-like spectrum, we can include this case in the comparison study, with minimal contamination by the P420 form.

Heme Domain with L-Arginine but Not BH4

The fast phase required at least two exponentials with rates differing by about a factor of 3.

Heme Domain with BH4 but Not L-Arginine

The fast phase required two exponentials differing by a factor of 4-5, slightly more than in the cases above.

Heme Domain with both L-Arginine and BH4

Adding both reagents had a striking effect, not observed with either alone. Geminate recombination was greatly reduced so that the slow phase dominated bimolecular association. The slow phase could be well fit using two exponentials differing in rate by a factor of 3, but instead of the amplitudes being approximately equal, the slower contributed 90% to the phase.

Holoenzyme with neither BH4 nor L-Arginine

The bimolecular combination could be described very well using three exponentials, of which two could be grouped into a fast phase. The fastest is 5-6 times faster than the middle rate, which is about 15 times the slowest. The fastest also provides most of the amplitude, being about 10 times larger than the middle rate, which is somewhat larger than the slowest. The primary data were obtained using a preparation that had been prepared entirely in the absence of Arg or BH4. This was possible by expressing nNOS holoenzyme in E. coli in which BH4 is not biosynthesized (18). Similar results were obtained for a sample originally prepared with BH4 and Arg but then dialyzed 48 h against buffer lacking both. That, however, is rather harsh treatment, required because BH4 is strongly bound.

Holoenzyme with L-Arginine but Not BH4

The slow phase was well fit using two exponentials differing by a factor of only 2, with the slower carrying the majority of the amplitude.

Holoenzyme with BH4 but Not L-Arginine

The slow phase required two exponentials differing in rate by a factor of about 3, with the slower carrying most (about 70%) of the weight.

Holoenzyme with L-Arginine and BH4

Protein with 0.25 mM BH4 and either 0.5 or 10 mM Arg showed no detectable geminate recombination (<3%), and the bimolecular association was almost exclusively (>95%) a slow phase. That phase, however, was not at all a single exponential (see Fig. 1). It required at least two and often three exponentials, distributed in rate over a range of 5-10. The small, fast phase was at least as fast as the fast phase observed in the absence of Arg. Although a small fast phase at the beginning of a decay can be quite apparent, it is difficult to characterize since it is affected by the random noise in the much larger slow phase.

Dissociation

Plots of association rates versus [CO] were linear and could be extrapolated to the y intercept, but the plots were not very useful for estimating dissociation constants. Only upper limits could be assigned, typically <20 s-1 for the fast phase and <2 s-1 for the slow phase. Consequently, we resorted to stopped-flow methods carried out at 20 °C. With one exception, all stopped-flow measurements were well fit by single exponentials over the time scale that could be monitored. A typical trace is shown in Fig. 2. The absorbance change of 0.07 is what was expected for the change from reactants to products, so that any additional component cannot be very large. Results are included in Table I. The holoenzyme lacking Arg showed about 15-20% faster phase with rate constant 15 ± 2. Perhaps surprisingly, the heme domain did not show any fast process under any conditions.


Fig. 2. Reaction time course for CO dissociation from the carboxy form of the heme domain protein of NOS, measured by NO replacement. lambda , 445 nm; T, 20 °C; pH, 7.8; Tris, 20 mM; EDTA, 100 µM; sodium dithionite, 50 µM; CO, 150 µM; NO, 100-250 µM; Arg, 10 mM. Symbols are data points, and the continuous line is the best fit to a single exponential, with parameters as in Table I. Note that this is total absorbance, not absorbance change, so it plateaus at a nonzero value at long times.
[View Larger Version of this Image (15K GIF file)]


DISCUSSION

Heterogeneous Kinetics

Bimolecular association required three or four exponentials to fit, but the exponentials tend to group into two distinct time regimes. Fig. 1 displays a large amplitude slow phase requiring three exponentials varying in rate by only a factor of 4 between the extreme values together, with a minority fast phase that is 150 times faster. We attribute the different rates to different protein conformations, each reacting with its own rate constant. This is based on two considerations. 1) Protein expression should produce proteins with a common sequence identity. 2) Lack of spectral evidence for intermediates appears to rule out sequential mechanisms. (Additional spectral features can be detected under certain circumstances but not the conditions discussed here.) We postulate a quasi-continuous distribution of conformations, which is bimodal, clustering about the two most common forms. The effect of cofactors, substrates, and other reagents is to bias the probability distribution in ways that can be described approximately (but not completely) as shifting the protein between fast reacting forms and slow reacting forms. The one previous report of combination kinetics in NOS distinguished the two major phases, but did not recognize additional complexity (16).

Tentatively assuming that kinetic heterogeneity is due to multiple conformational states, the question remains as to how much is intrinsic and how much may be an artifact of in vitro preparations. The much-studied myoglobins exhibit no such heterogeneity. The hemoglobins do have distinct R and T quaternary conformations with different kinetics (29, 30). There are also alpha  and beta  chain differences in Hb (30-32), but those reflect different primary sequences. On the other hand, NOS is compared not with Mb or Hb but with P450, and a similar multiplicity of kinetic phases is observed in human cytochrome P450 3A4 (33). Three phases were recognized over the range 105-107 M-1 s-1, and it was suggested there may be further sub-phases. Heterogeneity may be common for enzymes, even in vivo, although rare for the simpler oxygen transport proteins. On the other hand, the globins purified by simpler methods are less likely to introduce artifacts.

Bimolecular Association-Dissociation

Despite the complexity of association kinetics, Table I shows that the behavior can be modeled approximately using a fast phase with ka,f ~ 2-3 × 106 M-1 s-1 and a slow phase with ka,f ~ 2-4 × 104 M-1 s-1 in various ratios. Table I also shows that such a two-state model is an approximation. When the fast dominates, it becomes faster; when the slow state dominates, it is at its slowest. The full story apparently involves a complex redistribution among multiple states.

For the holoenzyme in the presence of both BH4 and Arg, the slow phase dominates but never quite reaches 100%. The opposite extreme occurs in the heme domain in the absence of BH4 and Arg. All other cases lie in between. The holoenzyme lacking both BH4 and Arg is similar to the heme domain, but adding either reagent is sufficient to convert it largely to the slow reacting form even though both are needed for the full effect. For the heme domain, either reagent alone has a small effect, but both together have a dramatic effect. The limiting cases are reproducible. For the intermediate cases, different preparations vary slightly in the percentage of each phase, again suggesting heterogeneity. Thus, it is not surprising that there is some slight difference between our numbers and those reported previously, which we believe pertain to two intermediate permutations (16).

It is known that both the BH4 and the Arg binding sites are located in the heme domain (19, 34). Our results confirm that model. It has also been reported that addition of Arg to the native enzyme is accompanied by changes in the UV-visible spectrum that are consistent with conversion of a low-spin, six-coordinate iron to a high-spin, five-coordinate iron (35). This implies that the Arg does not bind directly to the iron but perturbs or displaces some as yet unidentified ligand (or protein side chain) coordinated to the iron (35, 36). Steric interference by Arg of the ligand binding site was suggested (37). It is not surprising, therefore, that in the presence of both BH4 and Arg the CO association rate is reduced by a factor of 100 below the fast phase. Lack of any profound effect by Arg alone in the heme domain suggests that Arg can be incorporated into the heme domain, but its spatial orientation is different and directed away from the ligand binding site. In the holoenzyme, Arg alone is only slightly different from Arg + BH4, showing that the options for spatial orientation are severely limited. These conformational constraints may also correlate with quaternary changes, that is, the formation of homodimer, which is facilitated by BH4 and may be more favored in the holoenzyme.

The dissociation rate for CO from either NOS-CO itself or the heme domain subunit is close to 0.3 s-1 in the presence of Arg and 0.16 s-1 in its absence. Either number is high compared with the oxygen-binding hemeproteins but substantially less than in P450. Consider: five-coordinate carboxyheme (a protein-free model compound) in cetyltrimethylammonium bromide detergent has a CO dissociation rate constant of 400 s-1. Coordination of proximal imidazole trans to CO in either a protein or a model system reduces the CO dissociation constant to between 0.1 and 0.01 s-1, a trans effect of about 104. In six-coordinate bacterial cytochrome P450cam, in the absence of substrate, the CO dissociation rate constant is 6.3 s-1 (38) upon extrapolation to room temperature, which represents a much smaller trans effect of only 64. In the presence of substrate, CO dissociation from carboxy-P450 decreases to 1.4 s-1, but either with or without substrate the proximal cysteine thiolate in P450 exerts a much smaller positive trans effect than does imidazole. The dissociation rate for NOS-CO, which also has a trans-cysteine, is intermediate: lower than in bacterial cytochrome P450cam but much greater than when histidine (imidazole) is trans. The smaller trans effect of thiolate is probably due to its ionic nature; ions generally have lower affinity for ferroheme. In fact, for Mb mutants in which the proximal histidine was replaced by cysteine, carboxy derivative formation was accompanied by rupture of the proximal bond, a related consequence of a small (or negative) trans effect (39, 40).

We noted above the evidence for some steric crowding at the ligand binding site. As mentioned in the introduction, such crowding seems always to decrease rate constants for both association and dissociation of CO. Consequently, we can infer from kinetic measurements that the pocket in NOS is probably more constrained than in P450cam. The is no significant sequence homology between NOS and P450 and, according to computer modeling, no structural similarity. Consequently, differences in their heme pockets should be expected.

From the rate constants, one may calculate association equilibrium constants for the fast and slow phases of ~107 and ~105. (Conditions that generate some fraction of faster dissociation will reduce those numbers even further.) This affinity is 300 times less than in human Mb. It is unlikely that CO binding will be significant in vivo, where low CO concentrations are expected, unless there is much local enhancement and an absence of nearby traps, such as Hb or Mb. Ingi et al. (15) predicted, on the basis of different evidence, that a 1 mM concentration of CO would be needed to inhibit NOS activity. The equilibrium constant reported here suggests that their value is too high, but we concur with the conclusion that CO inhibition is unlikely. In fact, a recent study suggests just the opposite; far from being an inhibitor, small amounts of CO, 1000 ppm, may increase NO production (41). If that observation is confirmed, it will have important implications for the mechanism of NOS action. One possibility is that small amounts of CO may stabilize the reduced, ferro intermediate of iron in NOS in the reaction scheme proposed by Marletta (9).

Geminate Recombination

After bond dissociation by flash photolysis, CO may remain in the heme pocket for some time. The simplest model for such a geminate pair temporarily trapped in a protein pocket is as follows.
<UP>Fe—CO</UP> <AR><R><C>k<SUB>1</SUB></C></R><R><C>⇄</C></R><R><C>k<SUB>2</SUB></C></R></AR> [<UP>Fe</UP>⋯<UP>CO</UP>] <AR><R><C>k<SUB>3</SUB></C></R><R><C>⇄</C></R><R><C>k<SUB>4</SUB></C></R></AR> <UP>Fe</UP>+<UP>CO</UP> (Eq. 3)
In this model, admittedly oversimplified (24), the rate of disappearance of the caged pair intermediate is kgem = k2 + k3, whereas the geminate yield (the fraction recombining directly from the geminate pair) is Y = k2/(k2 + k3). For NOS, Table I shows that kgem is about 1.5 × 107 s-1 and varies by less than a factor of 2 among our eight conditions, while Y changes by a factor of 10, from <5% to almost 50%. This implies that k3 remains approximately constant, while k2, the rate of bond formation proper, changes significantly.

Changes in k2 may be due to changes either in the distal environment or in strain on the proximal ligand, both of which are influenced by protein conformation. Hence, geminate recombination provides direct evidence that in addition to other consequences, when the two domains are combined or BH4 or Arg is added, one effect is a substantial perturbation at the heme. Resonance Raman spectroscopy also probes the heme and its environment. Binding Arg causes a frequency shift of the Fe-NO stretch (37). Furthermore, that mode is broad, suggesting multiple conformations (7). All this echoes what is inferred from bimolecular combination rates.

The rates kgem are similar in all hemeproteins. In P450, two components were detected, one 2-3 times faster than we observe in NOS, the other 2-3 times smaller (42). For three different hemoglobins, kgem varied from approximately equal to about half the rate measured in NOS (43). Part of that variation is due to the fact that kgem = k2 + k3 is faster when k2 is greater. Escape from the protein k3 is surprisingly constant for different hemeproteins.

The geminate yield Y varies greatly among hemeproteins. For CO, Y is no more than about 5% in Mb, 20-50% in various HbR, depending on pH, and somewhat less in Hb under conditions that favor T-like structures (42). In P450cam, Y is reported to be as large as 90% in the absence of substrate, reduced to only 2% with camphor bound (41). NOS is like P450 in having very small Y with substrate bound. In the absence of substrate, P450 is like NOS without BH4 or like the heme domain and unlike the competent holoenzyme lacking only substrate. This echoes the conclusion from measurements by magnetic circular dichroism that "the active sites of NOS and P450 may share some common structural features, but significant distinctions exist" (44).

Geminate yield is related to bimolecular association. According to equation 3, 
k<SUB>a</SUB>=k<SUB>4</SUB>Y (Eq. 4)
where k4 is a bimolecular rate constant that describes entry into the protein and involves both diffusion to the periphery of the protein and a steric factor describing the fraction of encounters in which the ligand penetrates the protein. After entry, only a fraction Y of the ligands bind to the iron; the rest escape and must try again. From Table I, a correlation is evident between Y and the fraction of the bimolecular fast phase. The fast binding forms have substantial, ~50%, geminate yield, whereas slow binding conformations show very little geminate recombination. The 100-fold ratio in association rates between fast and slow phases can be accounted for largely by their different Y.

Geminate recombination strongly suggests that conformational effects on ligand binding are mediated not by changes in ligand entry into the protein but rather by interactions related to binding at the iron, involving strain from the proximal cysteine or polar or steric effects from the distal side.


FOOTNOTES

*   This work was supported in part by National Institutes of Health Grants HL13581 (to V. S. S.), HL40818 (to D. M.), GM52419, and HL30050 (to B. S. S. M.) and the Robert Welch Foundation (to B. S. S. M.).The costs of publication of this article were defrayed in part by the payment of page charges. The article must therefore be hereby marked "advertisement" in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.
par    To whom correspondence should be addressed.
1   The abbreviations used are: NOS, nitric-oxide synthase; nNOS, neuronal NOS; BH4, tetrahydrobiopterin; P450, cytochrome(s) P450; S/N, signal-to-noise ratio.

ACKNOWLEDGEMENTS

We thank T. Shea for excellent technical assistance and C. S. Roman for critical comments.


REFERENCES

  1. Palmer, R. M. J., and Moncada, S. (1989) Biochem. Biophys. Res. Commun. 158, 348-353 [Medline] [Order article via Infotrieve]
  2. Knowles, R. G., Palacios, M., Palmer, R. M., and Moncada, S. (1989) Proc. Natl. Acad. Sci. U. S. A. 86, 5159-5162 [Abstract]
  3. Mayer, B., Klatt, P., Böhme, E., and Schmidt, K. (1992) J. Neurochem. 59, 2024-2029 [Medline] [Order article via Infotrieve]
  4. McMillan, K., Bredt, D. S., Hirsch, D. J., Snyder, S. H., Clark, J. E., and Masters, B. S. S. (1992) Proc. Natl. Acad. Sci. U. S. A. 89, 11141-11145 [Abstract]
  5. Stuehr, D. J., and Ikeda-Saito, M. (1992) J. Biol. Chem. 267, 20547-20550 [Abstract/Free Full Text]
  6. White, K. A., and Marletta, M. A. (1992) Biochemistry 31, 6627-6631 [Medline] [Order article via Infotrieve]
  7. Wang, J., Stuehr, D. J., Ikeda-Saito, M., and Rousseau, D. L. (1993) J. Biol. Chem. 268, 22255-22258 [Abstract/Free Full Text]
  8. Masters, B. S. S., McMillan, K., Sheta, E. A., Nishimura, J. S., Roman, L. J., and Martásek, P. (1996) FASEB 10, 552-558 [Abstract/Free Full Text]
  9. Marletta, M. A. (1993) J. Biol. Chem. 268, 12231-12234 [Free Full Text]
  10. Stuehr, D. J., Kwon, N. S., Nathan, C. F., Griffith, O. W., Feldman, P. L., and Wiseman, J. (1991) J. Biol. Chem. 266, 6259-6263 [Abstract/Free Full Text]
  11. Masters, B. S. S. (1994) Annu. Rev. Nutr. 14, 131-145 [CrossRef][Medline] [Order article via Infotrieve]
  12. Abu-Soud, H. M., Wang, J., Rousseau, D. L., Fukuto, J. M., Ignarro, L. J., and Stuehr, D. J. (1995) J. Biol. Chem. 270, 22997-23006 [Abstract/Free Full Text]
  13. Dawson, T. M., and Snyder, S. H. (1994) J. Neurosci. 14, 5147-5150 [Abstract]
  14. Christodoulides, N., Durante, W., Kroll, M. H., and Shafer, A. I. (1995) Circulation 91, 2306-2309 [Abstract/Free Full Text]
  15. Ingi, T., Cheng, J., and Ronnett, G. V. (1996) Neuron 16, 835-842 [Medline] [Order article via Infotrieve]
  16. Matsuoka, A., Stuehr, D. J., Olson, J. S., Clark, P., and Ikeda-Saito, M. (1994) J. Biol. Chem. 269, 20335-20339 [Abstract/Free Full Text]
  17. Springer, B. A., Sligar, S. G., Olson, J. S., and Phillips, G. N., Jr. (1994) Chem. Rev. 94, 699-714
  18. Roman, L. J., Sheta, E. A., Martásek, P., Gross, S. S., Liu, Q., and Masters, B. S. S. (1995) Proc. Natl. Acad. Sci. U. S. A. 92, 8428-8432 [Abstract]
  19. McMillan, K., and Masters, B. S. S. (1995) Biochemistry 34, 3686-3693 [Medline] [Order article via Infotrieve]
  20. Martásek, P., Liu, Q., Roman, L. J., Gross, S. S., Sessa, W. C., and Masters, B. S. S. (1996) Biochem. Biophys. Res. Commun. 219, 359-365 [CrossRef][Medline] [Order article via Infotrieve]
  21. Sharma, V. S., Schmidt, M. R., and Ranney, H. M. (1976) J. Biol. Chem. 251, 4267-4272 [Abstract]
  22. Sheta, E. A., McMillan, K., and Masters, B. S. S. (1994) J. Biol. Chem. 269, 15147-15153 [Abstract/Free Full Text]
  23. Bredt, D. S., and Snyder, S. H. (1990) Proc. Natl. Acad. Sci. U. S. A. 87, 682-685 [Abstract]
  24. Walda, K. N., Liu, X. Y., Sharma, V. S., and Magde, D. (1994) Biochemistry 33, 2198-2209 [Medline] [Order article via Infotrieve]
  25. Sharma, V. S., Isaacson, R. A., John, M. E., Waterman, M. R., and Chevion, M. (1983) Biochemistry 22, 3897-3902 [Medline] [Order article via Infotrieve]
  26. Olson, J. S. (1981) in Hemoglobins: Methods in Enzymology (Antonini, E., Rossi-Bernardi, L., and Chiancone, E., eds), Vol. 76, pp. 631-651, Academic Press, New York
  27. Hurshman, A. R., and Marletta, M. A. (1995) Biochemistry 34, 5627-5634 [Medline] [Order article via Infotrieve]
  28. Wang, J., Stuehr, D. J., and Rousseau, D. L. (1995) Biochemistry 34, 7080-7087 [Medline] [Order article via Infotrieve]
  29. Sawicki, C. A., and Gibson, Q. H. (1976) J. Biol. Chem. 251, 1533-1542 [Abstract]
  30. Moore, E. G., and Gibson, Q. H. (1976) J. Biol. Chem. 251, 2788-2794 [Abstract]
  31. Sharma, V. S., Vedvick, T. S., Magde, D., Luth, R., Friedman, D., Schmidt, M. R., and Ranney, H. M. (1980) J. Biol. Chem. 255, 5879-5884 [Abstract/Free Full Text]
  32. Vandegriff, K. D., Le Tellier, Y. C., Winslow, R. M., Rohlfs, R. J., and Olson, J. S. (1991) J. Biol. Chem. 266, 17049-17059 [Abstract/Free Full Text]
  33. Koley, A. P., Buters, J. T. M., Robinson, R. C., Markowitz, A., and Friedman, F. K. (1995) J. Biol. Chem. 270, 5014-5018 [Abstract/Free Full Text]
  34. Nishimura, J. S., Mártasek, P., McMillan, K., Salerno, J. C., Liu, Q., Gross, S. S., and Masters, B. S. S. (1995) Biochem. Biophys. Res. Commun. 210, 288-294 [CrossRef][Medline] [Order article via Infotrieve]
  35. McMillan, K., and Masters, B. S. S. (1993) Biochemistry 32, 9875-9880 [Medline] [Order article via Infotrieve]
  36. Salerno, J. C., Frey, C., McMillan, K., Williams, R. F., Masters, B. S. S., and Griffith, O. W. (1995) J. Biol. Chem. 270, 27423-27428 [Abstract/Free Full Text]
  37. Wang, J., Rousseau, D. L., Abu-Soud, H. M., and Stuehr, D. J. (1994) Proc. Natl. Acad. Sci. U. S. A. 91, 10512-10516 [Abstract/Free Full Text]
  38. Peterson, J. A., and Griffin, B. W. (1972) Arch. Biochem. Biophys. 151, 427-433 [Medline] [Order article via Infotrieve]
  39. Hildebrand, D. P., Ferrer, J. C., Tang, H.-L., Smith, M., and Mauk, A. G. (1995) Biochemistry 34, 11598-11605 [Medline] [Order article via Infotrieve]
  40. Matsui, T., Nagano, S., Ishimori, K., Watanabe, Y., and Morishimi, I. (1996) Biochemistry 35, 13118-13124 [CrossRef][Medline] [Order article via Infotrieve]
  41. Meilin, S., Rogatsky, G. G., Thom, S. R., Zarchin, T. N., Guggenheimer-Furman, E., and Mayevsky, A. (1996) J. Appl. Physiol. 81, 1078-1083 [Abstract/Free Full Text]
  42. Tian, W. D., Wells, A. V., Champion, P. M., Di Primo, C., Gerber, G., and Sligar, S. G. (1995) J. Biol. Chem. 270, 8673-8679 [Abstract/Free Full Text]
  43. Campbell, B. F., Magde, D., and Sharma, V. S. (1985) J. Biol. Chem. 260, 2752-2756 [Abstract]
  44. Sono, M., Stuehr, D. J., Ikeda-Saito, M., and Dawson, J. H. (1995) J. Biol. Chem. 270, 19943-19948 [Abstract/Free Full Text]

©1997 by The American Society for Biochemistry and Molecular Biology, Inc.