| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
* European Media Laboratory, Schloss-Wolfsbrunnenweg 33, D-69118 Heidelberg, Germany; Department of Biochemistry and Molecular Biology, Syddansk Universitet, DK-5230 Odense M, Denmark; and Department of Biological Sciences, Wayne State University, Detroit, Michigan 48202 USA
Correspondence: Address reprint requests to Lars F. Olsen, Dept. of Biochemistry and Molecular Biology, Syddansk Universitet, Campusvej 55, DK-5230 Odense M, Denmark. Tel.: +45-6550-2482; Fax: +45-6550-2467; E-mail: lfo{at}dou.dk.
 |
ABSTRACT |
|---|
 TOP ABSTRACT INTRODUCTIONMATERIALS AND METHODSResultsDISCUSSIONACKNOWLEDGEMENTSREFERENCES |
|---|
|
INTRODUCTION |
|---|
|
TOPABSTRACTINTRODUCTIONMATERIALS AND METHODSResultsDISCUSSIONACKNOWLEDGEMENTSREFERENCES |
|---|
). Electrons are transferred from cytoplasmic NADPH to oxygen on the phagosomal side of the membrane, generating the so-called reactive oxygen species, e.g., superoxide, hydrogen peroxide, hydroxyl radical, and singlet oxygen. The NADPH used is itself produced by the hexose monophosphate shunt. Although there is good evidence that the respiratory burst is associated with NADPH oxidase activity, other proteins also participate in the oxidative destruction of targets. One protein is myeloperoxidase, which is, in fact, the most abundant protein of the neutrophil (Klebanoff, 1991). Myeloperoxidase, which is stored in the azurophil granules from which it is released into the phagosome upon activation, is believed to catalyze the reaction of hydrogen peroxide with Cl- ions to form HOCl, which also has antibacterial activity (Dunford, 1999; Hampton et al., 1998; Klebanoff, 1991).
Recently, it was shown that in migrating neutrophils the concentrations of NAD(P)H and reactive oxygen species oscillate (Petty, 2001; Amit et al., 1999). The mechanism responsible for such oscillations is not known. However, various known biochemical oscillators such as glycolysis and cytosolic calcium (Petty, 2001) have been proposed as sources of the oscillations. It is quite likely that several biochemical oscillators operate together to produce the experimentally observed oscillations. Here we wish to present evidence for another such oscillator that involves myeloperoxidase, melatonin, NADPH, and NADPH oxidase.
In a previous paper (Olsen et al., 2001) it was shown that melatonin activates the reaction:
 |
). In addition to activating the PO reaction. melatonin also induces oscillations in this reaction when the pH is between 4 and 6. Because lactoperoxidase and myeloperoxidase are structurally very similar (Dunford, 1999) and because melatonin has been shown to be a substrate for the classical peroxidase reaction:
 |
; Allegra et al., 2001, Ximenes et al., 2001), it is tempting to suggest that a myeloperoxidase-catalyzed PO reaction could also be involved in the oscillations in neutrophils. Myeloperoxidase has been shown to catalyze the PO reaction (Fujimoto et al., 1993). Furthermore, peroxidases do not discriminate between NADH and NADPH as the reducing substrate (Nakamura et al., 1969). Finally, it has been shown that leukocytes may synthesize melatonin (Finocchiaro et al., 1991). Here we first study experimentally the PO reaction catalyzed by horseradish peroxidase enclosed in liposomes while the substrates NADH and O2 are continuously supplied to the external medium. Melatonin is added to the reaction mixture in order to mediate electron trafficking between NADH and peroxidase. We then use the experimental results to develop a model of the PO reaction catalyzed by myeloperoxidase in neutrophils. The computational results are tested by performing key experiments on neutrophils confirming the main results of the simulations.
|
MATERIALS AND METHODS |
|---|
|
TOPABSTRACTINTRODUCTIONMATERIALS AND METHODSResultsDISCUSSIONACKNOWLEDGEMENTSREFERENCES |
|---|
-phosphatidylcholine (based on palmitoyl-oleoyl phosphatidylcholine) (>99%) was purchased from Avanti Polar Lipids (Alabaster, AL).
Liposomes
Liposomes were made of palmitoyl-oleoyl phosphatidylcholine. The lipid was first dissolved in chloroform and placed in a vacuum chamber overnight. The resulting lipid film was then hydrated by adding an aqueous solution containing 20 µM horseradish peroxidase and 0.1 M sodium acetate buffer, pH 5.1. The final concentration of lipid in the suspension was 10 mM. The lipid suspension was placed in a heating bath at 35°C for at least 2 h. The suspension was vigorously shaken from time to time. Unilamellar vesicles containing horseradish peroxidase were made by extrusion (MacDonald et al., 1991). The lipid suspension was passed through a Millipore HTTP (Millipore, Glostrup, Denmark) filter with a pore size of 0.4 µM under an N2 pressure of 10 atmospheres 11 times. The liposomes were separated from the peroxidase-containing solution by elution through a Sephacryl S-200 HR (Amersham Biosciences, Hørsholm, Denmark) column using 0.1 M sodium acetate, pH 5.1, as the eluent.
PO reaction
Experiments were conducted at 28 (±0.1)°C in a 2.0 x 2.0 x 4.3 cm quartz cuvette fitted with a thermostating jacket. The cuvette was connected to a Zeiss Specord S 10 diode array spectrophotometer (Carl Zeiss, Jena, Germany) through optical fibers. Oxygen in the solution was measured with a Clark-type oxygen electrode (Microelectrodes, Bedford, NH). The reaction mixture consisted of an 8 ml well-stirred homogeneous aqueous solution containing 0.1 M sodium acetate, pH 5.1, 200 µM melatonin and peroxidase-containing liposomes (see above). O2 is entering the reaction mixture from a 1.05% (v/v) O2/N2 gas mixture above the liquid. The rate of diffusion into the liquid follows Fick's law and the oxygen transfer constant was measured as 5.6 x 10-3 (±0.1 x 10-3) s-1 corresponding to a stirring rate of 900 (±10) rpm. NADH was supplied to the reaction mixture by infusion of a 0.1 M stock solution through a capillary whose tip was below the surface of the liquid. The infusion was mediated by a Harvard Apparatus (Holliston, MA) model 22 infusion pump and the infusion rate was 35 µl/h.
Cells
Human peripheral blood neutrophils were prepared as previously described (Amit et al., 1999). Cells were suspended in Hanks' balanced salt solution containing glucose but not phenol red. Morphologically polarized cells with well-defined uropods, which are characteristics of cell motility, were examined.
Detection of NAD(P)H oscillations
Metabolic oscillations were detected using NAD(P)H microscopy (Amit et al., 1999; Rosenspire et al. 2000). NAD(P)H autofluorescence was detected using narrow bandpass interference filters (AF series, Omega Optical, Brattleboro, VT) and a dichroic mirror. An iris diaphragm was adjusted to exclude light from neighboring cells. A cooled photomultiplier tube held in a Products for Research (Danvers, MA) housing attached to a Zeiss microscope was used (Rosenspire et al., 2000).
Imaging spectrophotometry of chemiluminescence
Cells were observed microscopically at 37°C. Imaging spectrophotometry was performed using an Axiovert 135 fluorescence microscope with a quartz condenser and quartz objective (Carl Zeiss, New York). To increase light collection efficiency, the microscope's bottom port was employed. This port was fiberoptically coupled to the input of an Acton 150 (Acton Instruments, Acton, MA) imaging spectrophotometer (Petty et al., 2000). The exit side was connected to a liquid N2-cooled intensifier attached to a Peltier-cooled IMAX-512 camera (-20°C) (Princeton Instruments, Trenton, NJ). The camera was controlled by a high-speed Princeton ST-133 interface and a Stanford Research Systems (Sunnyvale, CA) DG-535 delay gate generator (Petty et al., 2000). A Dell Precision 410 workstation with an 800 MHz processor, 1.0 Gb RAM, 10 Mb onboard cache, and a high-speed Lava Dual PCI enhanced port (Rexdale, CA) was used. Winspec (Princeton) software was employed.
|
Results |
|---|
|
TOPABSTRACTINTRODUCTIONMATERIALS AND METHODSResultsDISCUSSIONACKNOWLEDGEMENTSREFERENCES |
|---|
; Kummer et al., 1997; Hauser and Olsen, 1998). However, unlike melatonin, which is uncharged at pH 5.1, 4-hydroxybenzoic acid is a polar molecule with a net negative charge under the present experimental conditions. The latter could prevent 4-hydroxybenzoic acid from diffusing across the liposomal membrane and thus inducing an acceleration in the PO reaction. Consequently, our results suggest that peroxidase is shielded from NADH in the solution by a hydrophobic barrier.
Next we investigated whether the peroxidase-containing liposomes could catalyze an oscillating PO reaction in the presence of melatonin. Fig. 1 shows a time series of the recorded absorbance changes at 360 nm and 418 nm plus the oxygen concentration when NADH is infused into a reaction mixture containing liposomes and melatonin. The absorbance changes at the two wavelengths are characteristic for NADH and peroxidase compound III, respectively. We note that the reaction oscillates, although the oscillations are strongly damped. Using an extinction coefficient for compound III at 418 nm of 1.16 x 105 M-1 cm-1 (Olsen et al., 2001), we can estimate a concentration of peroxidase in the suspension of 84 nM. This corresponds to a concentration of peroxidase in the liposome of 20 µM as expected if peroxidase is evenly distributed between the external solution and the liposome aqueous fraction after the preparation of the liposomes.
|
Construction of the model
In the construction of a model for a myeloperoxidase-catalyzed PO reaction in neutrophils, one has to consider that the enzyme and one of its substrates (NADPH) are not located in the same compartment. Although NADPH is formed in the cytoplasm, the enzyme is situated either in the azurophil granules or in the phagosome. In the previous section we showed that it is possible for peroxidase to catalyze the PO reaction even when the substrate NADH and the enzyme are separated by a lipid membrane, provided that melatonin (or another nonpolar activator of peroxidase) is also present. Thus, we break our model of the PO reaction in neutrophils up into reactions that can occur in the cytosol and those that occur in the phagosome. Most of the reactions are well-documented and the in vitro rate constants are known (see Scheeline et al., 1997, for a review). The reactions and diffusion terms are all listed in Table 1 and the equation numbers listed below all refer to the numbers in the table.
|
 | (7) |
The reaction is slow at pH around 7, but it may be catalyzed by a number of organic substances, e.g., methylene blue (Sevcik and Dunford, 1991). Here the reaction is modeled as a second order reaction with a very low rate constant (of the order of 1 M-1 s-1).
Another reaction likely to take place in the cytosol is the reduction of molecular oxygen by the NADP· radical:
 | (8) |
The radical NADP· is in turn regenerated in the reaction:
 | (10) |
Superoxide (O2-) formed in reaction (8) either disproportionates:
 | (9) |
 | (18) |
As for the disproportionation of superoxide, we have assumed second order kinetics, but we have also set the second order rate constant more than 50 times as high as it should be. This is done to account for the superoxide dismutase activity in the cytoplasm. As for the diffusion of superoxide from cytoplasm to phagosome the rate is probably extremely slow. We consider the rate of this process to be essentially zero in our simulations.
We also assume that the NADP· radical dismutates to form a (NADP)2 dimer (reaction (11)). Even though the rate constant for this reaction is fairly high, the actual rate of (NADP)2 formation is low because the concentration of the radical is always very low. We have not included any oxidation of the dimer to NADP+. Some enzymes such as peroxidase are able to use (NADP)2 as a substrate (Kirkor and Scheeline, 2000).
NADPH production is assumed to come exclusively from the hexose monophosphate shunt (Henderson and Chappell, 1996). In the absence of any further knowledge about the dynamics of this reaction pathway, we assume that NADPH production is constant (reaction (12)). The entry of O2 from the neighboring tissue into the cytoplasm of the neutrophil is considered to be due to simple diffusion:
 | (13) |
Here we have taken a very conservative view on the maximum rate of oxygen diffusion into the cell (
0.75 mM per minute), which is small compared to the rate of oxygen consumption in cardiac tissue (McCord, 2000).
Reactions occurring in the phagosome
These include the classical peroxidase cycle:
 | (1) |
 | (2) |
 | (3) |
Note that reaction (1), in contrast to the corresponding reaction involving horseradish peroxidase, is reversible (Furtmüller et al., 2000). The rate constants for reactions (1), k1 and k-1 listed in Table 1 are taken from Furtmüller et al. (2000), whereas the rate constants k2 and k3 are taken from Allegra et al. (2001).
The free ferric enzyme may form compound III by reaction with superoxide:
 | (4) |
In the absence of any reported rate constant for myeloperoxidase we have used the corresponding rate constant k4 measured for horseradish peroxidase (Scheeline et al., 1997). In homogeneous solution there is an alternative way of forming compound III through reduction of ferric peroxidase to ferrous peroxidase by the NAD(P)· radical and subsequent oxygenation of ferrous peroxidase (Bronnikova et al., 1995). However, as the NAD· radical is probably not able to enter the phagosome, these reactions are not included in the present model. As in the cytosol superoxide may disproportionate (reaction (5)), but in this case the reaction is not considered to be catalyzed by superoxide dismutase. In most detailed models compound III is reduced to compound I through reaction with the NAD(P)· radical. For the same reason as given above for the formation of ferrous peroxidase this reaction is probably not occurring in the phagosome. Instead we have included the reaction:
 | (6) |
. In the absence of detailed knowledge about the magnitude of k6, we have chosen a reasonably low value for this reaction.
Diffusion terms
The present model assumes exchange of electrons between components in two different cellular compartments. The following species must be able to permeate the lipid membrane, either by simple diffusion or by mediated transport: oxygen and hydrogen peroxide, melatonin and its corresponding free radical, and superoxide. The latter molecule will probably cross the lipid membrane as the uncharged protonated form (HO2·). However, we cannot exclude that the phagosome membrane, like the plasma membrane, contains an anion channel through which the unprotonated form can enter the phagosome. In any case the rate of diffusion of O2- across the phagosome membrane will be considered to be very slow. The diffusional terms are computed from Fick's law of diffusion:
 | (6) |
 | (6) |
) suggesting a first order constant of 20–50 s-1. We can easily assume that the permeability coefficient of O2 is no less than that. As for the permeability coefficients of melatonin and its free radical form, we have assumed values around 10-4 cm/s resulting in pseudo first order constants of the order of 10 s-1. It was previously shown that melatonin can easily cross lipid bilayers, suggesting a high permeability coefficient (Costa et al., 1995).
NADPH oxidase
Upon activation of neutrophils, this membrane-associated multisubunit enzyme is recruited and assembled in the phagosome membrane (Hampton et al., 1998; Henderson and Chappell, 1996; Jones et al., 2000). The function of the enzyme is to produce superoxide radicals in the phagosome. The electrons needed for the reduction of oxygen come from NADPH, which again is formed by the hexose monophosphate shunt. The stoichiometry of the reaction is:
 | (6) |
NADPH oxidase is a very complex enzyme, the activity of which requires the involvement of at least two G-proteins (Jones et al., 2000). These proteins act as molecular switches that turn the enzyme's activity on and off when they alternate between binding GTP in the active form and GDP in the inactive form. Thus the activity of NADPH oxidase undoubtedly has a very complex kinetics with respect to its two substrates, NADPH and O2. Here we have taken a simple nonlinear dependence of the enzyme kinetics of NADPH oxidase, namely that the enzyme shows cooperativeness with respect to the binding of NADPH and simple Michaelis-Menten behavior with respect to O2 (see Table 1). The dissociation constant for the binding of NADPH is set to 60 µM (Henderson and Chappell, 1996), and the Michaelis constant for oxygen is assumed to be around 1.5 µM. The maximal rate of the oxidase was estimated at 288 µM/s in accordance with values of around 10 nmol/106 cells per 10 min reported in Levy et al. (1994). However, it should be emphasized that the exact kinetic expression for NADPH oxidase and the constants assumed for the binding of its substrates as well as its maximal turnover number are not critical for the behavior of the model.
The differential equations corresponding to the model presented in Table 1 are listed in Table 2. The concentration changes in the phagosome are indicated with a subscript p whereas the concentration changes in the cytosol are indicated with a subscript c. The factor f refers to the fractional volume of the phagosome over the cytosol volume and is set to 1:10 in all simulations. However, both higher and lower values of f yield oscillating solutions of the equations in Table 2.
|
). The melatonin free radical also oscillates with the same frequency as the oscillations shown in Fig. 2.
|
 |
The NADPH produced is mainly used for the reduction of oxygen to reactive oxygen species, i.e., hydrogen peroxide, superoxide, and hydroxyl radical. We have simulated the flux of production of NADPH as a single constant, k12, the value of which is increasing as the neutrophils are activated. Fig. 3 shows the effect on the oscillation period when k12 is increased from 23.5 to 34 µM/s. Outside this interval the system displays only steady-state behavior. We note that as k12 is increased, the oscillation period decreases from around 200 s to between 15 and 20 s. These periods are in the same range as those observed experimentally in migrating neutrophils (Amit et al., 1999).
|
|
|
|
):
 | (1) |
 | (2a) |
 | (2b) |
Obviously, the reactions 2 a and 2 b compete with reactions 2 and 3 of our model. Furthermore, it was shown (Hampton et al., 1998) that a substantial fraction of the H2O2 formed in the phagosome is converted into HOCl. In order to check if these reactions have a profound impact on the modeling results, we included them in the model and studied the resulting simulations. Using the experimentally determined kinetic rate constants for reactions 2 a and 2 b (Furtmüller et al., 2000), we observed that including these reactions in the model does indeed inhibit the oscillatory behavior. However, with the kinetic parameters used, almost 100% of the oxygen taken up by the neutrophils would be incorporated into HOCl. This does not match corresponding experimental data showing that only 20–70 percent of H2O2 is incorporated. An explanation for this could be that the parameters that were determined experimentally in vitro might be different in the in vivo system. We therefore reduced the rate constant for the forward reaction in Eq. 2 a to such a degree that up to 70% of the oxygen consumed in the phagosome is used for the production of HOCl. Under these circumstances, the oscillatory behavior is restored (data not shown) and the above described modeling results remain valid.
Experiments with cells
We next sought to test certain assumptions of the model described above and its predictions based upon the simulation results. The model replicates key features of an activated neutrophil, such as the reliance on the hexose monophosphate shunt (HMS), activation of the NADPH oxidase, and mobilization of myeloperoxidase. In the first series of cell experiments, we sought to establish the contribution of the HMS in NAD(P)H oscillations. To accomplish this goal, we employed the reagent 6-aminonicotinamide (6-AN), which is known to block the activity of glucose-6-phosphate dehydrogenase and 6-phosphogluconate dehydrogenase; both of these enzymes are specific for the HMS (White et al., 1978).
Adherent neutrophils were incubated in the presence or absence of 6-AN at 10 mM for 30 min at 37°C; previous studies have shown that these conditions effectively block neutrophil activation and the respiratory burst (e.g., Bender and Van Epps, 1985). In the absence of added reagents, NAD(P)H oscillations with a period of 20 s are observed. However, the activating chemotactic peptide N-formyl-methionyl-leucyl-phenylalanine (FMLP, 100 nM), which induces a vigorous respiratory burst, triggers the appearance of a 10 s period in NAD(P)H oscillations (Fig. 7), which has been previously observed (Petty, 2001; Amit et al., 1999). The reagent 6-AN has no effect on the 20 s oscillations, but eliminates the acquisition of the 10 s NAD(P)H oscillation in the presence of FMLP (Fig. 7). Parallel results were obtained using the somewhat less specific reagent dexamethasone (Fig. 7), which also inhibits glucose-6-phosphate dehydrogenase (Whitcomb and Schwartz, 1985; Rom and Harkin, 1991; Mohan and Jacobson, 1993). As previously suggested using various neutrophil activating substances (Petty, 2001; Amit et al., 1999), these findings indicate that the 20 s oscillation is associated with cell adherence/motility, but not activation or utilization of the HMS, whereas the 10 s oscillation is only found in activated cells and draws upon the HMS.
|
). Cells were incubated with DPI at various concentrations for 45 min at 37°C. In the range of 0–20 µg/ml DPI had no effect on NAD(P)H oscillations (Fig. 8 A). At a higher dose of 25–35 µg/ml DPI, FMLP had no effect on NAD(P)H oscillations of neutrophils (Fig. 8 B). No NAD(P)H oscillations were observed at doses above 40 µg/ml DPI (Fig. 8 C). Under the conditions we have employed (adherent cells, HBSS, etc.), concentrations in the range of 100 µg/ml appeared to be toxic to the cells. Our findings confirm the role of the NADPH oxidase in establishing the 10 s NAD(P)H oscillation period.
|
). Therefore, neutrophils were incubated with 50 µM melatonin for 30 min at 37°C followed by addition of 50 nM FMLP. The oxidation of melatonin by peroxidase is accompanied by chemiluminescence via a dioxetane intermediate (Silva et al., 2000). The chemiluminescence properties of melatonin in neutrophils were examined using an imaging spectrophotometer system. The chemiluminescence spectrum was centered at 515 nm (Fig. 9 C). Using an intensifier and camera cooled to -20°C, the chemiluminescence was imaged by accumulating data for 10 s. These images showed that FMLP-treated cells exhibited chemiluminescence within the anterior region of the cell body and lamellipodium (Fig. 9 B), which are areas rich in intracellular granules and peroxidase. When the chemiluminescence was quantitated as a function of time, it was found to oscillate (Fig. 9 D), as does hydrogen peroxide, a cofactor in its oxidation. Therefore, melatonin undergoes chemical reactions within living human neutrophils.
|
; Amit et al., 1999). The level of ROM production is enhanced by activating agents such as FMLP and further enhanced by the inclusion of melatonin (Fig. 10, B, C, E, and F). This accounts for the "priming" effect of melatonin reported in the literature (Recchioni et al., 1998, Morrey et al., 1994); melatonin is incapable of initiating a respiratory burst, but is capable of enhancing ROM production in the presence of a second activating signal.
|
|
DISCUSSION |
|---|
|
TOPABSTRACTINTRODUCTIONMATERIALS AND METHODSResultsDISCUSSIONACKNOWLEDGEMENTSREFERENCES |
|---|
Previous experimental studies have shown that polarized neutrophils exhibit NAD(P)H oscillations with a period of 20 s for untreated cells and 10 s for activated neutrophils (Petty, 2001; Amit et al., 1999). These two distinct periods apparently correspond to two different biochemical pathways. The 20 s oscillatory period presumably reflects the glycolytic apparatus (the nonHMS, 6-AN-insensitive component (Fig. 7)) whereas the 10 s oscillation corresponds to the HMS/NADPH oxidase/myeloperoxidase-dependent (6-AN-sensitive) component of metabolism for activated neutrophils, as modeled above.
The model predicts metabolic oscillations whose period is consistent with experimental observations. The central feature of the model is the cooperation of the NADPH oxidase and myeloperoxidase. The requirement of the NADPH oxidase to stabilize the oscillatory pattern has been confirmed by inactivating the oxidase using DPI, an inhibitor of the oxidase (Cross and Jones, 1986). As predicted by the modeling results, higher concentrations of DPI blocked the 10 s oscillations. Interestingly, at a DPI dose of 25–35 mg/ml, the 20 s oscillation was observed in the presence of FMLP. We suggest that this concentration range is associated with inhibition of the stable oscillations linked with NADPH oxidase-myeloperoxidase cooperation, but that the glycolytic pathway is not significantly affected. At higher DPI concentrations, additional metabolic pathways might be affected.
The model predicted that melatonin enhances the amplitude of NAD(P)H oscillations. The melatonin effect was subsequently tested in neutrophils. Our results showed that melatonin, indeed, enhanced the amplitude of the 10 s NAD(P)H oscillations of activated neutrophils, but had no effect on the 20 s NAD(P)H oscillations. Not only does this confirm the model described above, but it is also consistent with the association of the 10 and 20 s NAD(P)H oscillations with two distinct metabolic pathways. Thus, the model has successfully accounted for prior experimental observations and predicted the effects of DPI and melatonin on living cells.
Oscillations in NAD(P)H are accompanied by oscillations in superoxide and hydrogen peroxide production (Fig. 2). Changes in NAD(P)H oscillation frequency or amplitude are paralleled by changes in oxidant production. The low levels of oxidant production shown above for nonactivated neutrophils displaying the 20 s oscillation have been previously observed by several laboratories (Demaurex et al., 1996; Ginis and Tauber, 1990; Kindzelskii et al., 1998). These oscillations are apparently independent of the HMS. Activation of the HMS is accompanied by higher frequency oscillations and enhanced oxidant production.
The modeling studies described above characterize the properties of activated neutrophils; hence, they do not address the 20 s NAD(P)H oscillations. Also, the modeling results do not account for the dominance of the 10 s oscillations in the presence of neutrophil activating stimuli. We believe that the disappearance of the 20 s component may be due to a greater flux of carbon through the HMS relative to glycolysis. Further analysis of the glycolytic pathway in neutrophils should provide further insights into the coupling of these two oscillators. Importantly, combined modeling/experimental studies should help explain cell-mediated damage due to high glucose levels in diabetes and the unusual metabolic changes in neutrophils during polytrauma (Oehler et al., 2000).
The effect of melatonin on activated leukocytes is interesting for several reasons. Melatonin is a pineal hormone and a biosynthetic product of leukocytes (Finocchiaro et al., 1991). It has an unexplained ability to prime neutrophils and monocytes for oxidant production, cytotoxicity, and other functions (Recchioni et al., 1998; Morrey et al., 1994). Melatonin, a priming agent, enhances oxidant production in response to activating stimuli, such as lipopolysaccharide and phorbol myristate acetate (PMA), but does not induce a response in the absence of an activating stimulus. Previous studies have demonstrated the priming of oxidant production using flow cytometric and superoxide dismutase-inhibitable cytochrome c reduction methodologies (Recchioni et al., 1998; Morrey et al., 1994). In the present study we confirm these findings using single cell fluorescence techniques and correlate this enhancement with an increase in the amplitude of NAD(P)H oscillations. Other studies have shown that activated neutrophils oxidize melatonin (Silva et al., 2000; Recchioni et al., 1998), which suggests that chemical reactions of melatonin participate in leukocyte priming.
The modeling studies also suggest a likely chemical mechanism of melatonin-induced neutrophil priming. The enhanced NADPH and oxidant oscillatory amplitudes predicted by the model and observed experimentally explain why the levels of oxidant production are higher in the presence of melatonin. Because the model is only applicable to activated neutrophils, nonactivated neutrophils are not susceptible to melatonin-mediated enhancement. This is mainly due to the fact that the model predicts that melatonin enhances the myeloperoxidase catalyzed oxidation of NAD(P)H. Myeloperoxidase is only released from granules upon activation of the neutrophils. Therefore, as a priming agent, it cannot enhance oxidant production in the absence of an activating stimulus.
The ability of melatonin to enhance oxidant production by leukocytes may seem paradoxical due to its ability to scavenge reactive oxygen and nitrogen species (Matuszak et al., 1997; Zhang et al., 1998). In a physiological context, however, these contrasting abilities may be advantageous. In the local environment of a neutrophil, melatonin would promote higher oxidant levels, thereby enhancing host defense. However, at some distance from a neutrophil, melatonin's ability to act as a scavenger may minimize collateral damage to healthy host tissues.
Enhanced amplitudes of NAD(P)H and superoxide oscillations have also been observed in neutrophils in response to the biological factors interferon-
(Adachi et al., 1999; Rosenspire et al., 2000) and interleukin-12 (unpublished) and other stimulatory conditions provided by PMA and electric fields (Kindzelskii and Petty, 2000). Not surprisingly, at appropriate doses interferon- and PMA prime cells for oxidant production (Tyagi et al., 1988; Berton et al., 1986). Interferon- and interleukin-12 transmit intracellular signals via calcium (Aas et al., 1998; Aas et al., 1999; Collison et al., 1998). Furthermore, PMA acts on the calcium signaling circuitry (Tauber, 1987). Thus, similar metabolic responses can be achieved by several calcium-related stimuli. One possibility to account for this is that calcium signals may promote mixing of the NADPH oxidase and myeloperoxidase into the same compartment via membrane fusion whereas melatonin may promote electron trafficking in the absence of vesicle fusion. Alternatively, because interferon- induces the synthesis of melatonin (Finocchiaro et al., 1988), it is possible that melatonin is the mediator responsible for the interferon--mediated enhancement of NAD(P)H and oxidant oscillations. Because melatonin is relatively mobile, one might discriminate between these two possibilities by mixing interferon--primed cells with untreated cells then assaying for amplitude modulation in the untreated cells.
One problem with the conventional reductionist approach to rational drug development is that a specific compound designed to affect a specific enzyme often has just the opposite of the desired effect at the cellular level or in patient care. A goal of biocomplexity research is the identification of important interactions, pathways, feedback loops, etc. such that biochemical pathways and living tissues can be intelligently perturbed. Due to their relative simplicity and profound clinical importance, neutrophils are an excellent system for computational analysis. Our model of neutrophil activation was found to be consistent with earlier experimental findings and successfully predicted the novel experimental results described above (Figs. 7–10). The apparent resolution of the mechanism of melatonin-mediated neutrophil priming is a useful example of the importance of computer simulations in biology. The interplay of in silico and in vitro research will have a major impact on modern biology and medicine.
|
ACKNOWLEDGEMENTS |
|---|
|
TOPABSTRACTINTRODUCTIONMATERIALS AND METHODSResultsDISCUSSIONACKNOWLEDGEMENTSREFERENCES |
|---|
Simulations were performed using the Rosenbrock routine for the numerical solution of stiff differential equations. The software used was MADONNA (University of California, Berkeley, CA).
Submitted on March 19, 2002; accepted for publication September 4, 2002.
|
REFERENCES |
|---|
|
TOPABSTRACTINTRODUCTIONMATERIALS AND METHODSResultsDISCUSSIONACKNOWLEDGEMENTSREFERENCES |
|---|
Aas, V., K. Larsen, and L. G. Iversen. 1999. Interferon-gamma elicits a G-protein-dependent Ca2+ signal in human neutrophils after depletion of intracellular Ca2+ stores. Cell. Sign. 11:101–110.[Medline]
Adachi, Y., A. L. Kindzelskii, N. Ohno, T. Yadomae, and H. R. Petty. 1999. Amplitude and frequency modulation of metabolic signals in leukocytes: synergistic role in interferon- in interleukin-6 and interleukin-2-mediated cell activation. J. Immunol. 163:4367–4374.
Allegra, M., P. G. Furtmüller, G. Regelsberger, M. L. Turco-Liveri, L. Tesoriere, M Peretti, M. A. Livrea, and C. Obinger. 2001. Mechanism of the reaction of melatonin with human myeloperoxidase. Biochem. Biophys. Res. Commun. 282:380–386.[Medline]
Amit, A., A. L. Kindzelskii, J. Zanoni, J. N. Jarvis, and H. R. Petty. 1999. Complement deposition on immune complexes reduces the frequencies of metabolic, proteolytic, and superoxide oscillations of migrating neutrophils. Cell. Immunol. 194:47–53.[Medline]
Baker, C. J., K. Deahl, J. Domek, and E. W. Orlandi. 2000. Scavenging of H2O2 and production of oxygen by horseradish peroxidase. Arch. Biochem. Biophys. 382:232–237.[Medline]
Bender, J. G., and D. E. Van Epps. 1985. Inhibition of human neutrophil function by 6-aminonicotinamide: the role of the hexose monophosphate shunt in cell activation. Immunopharm. 10:191–199.
Berton, G., L. Zeni, M. A. Cassatella, and F. Rossi. 1986. Gamma interferon is able to enhance the oxidative metabolism of human neutrophils. Biochem. Biophys. Res. Commun. 138:1276–1282.[Medline]
Bronnikova, T. V., V. R. Fed'kina, W. M. Schaffer, and L. F. Olsen. 1995. Period-doubling bifurcations and chaos in a detailed model of the peroxidase-oxidase reaction. J. Phys. Chem. 99:9309–9312.
Collison, K., S. Saleh, R. Parhar, B. Meyer, A. Kwaasi, K. Al-Hussein, S. Al-Sedairy, and F. Al-Mohanna. 1998. Evidence that IL-12-activated Ca2+ and tyrosine signaling pathways in human neutrophils. J. Immunol. 161:3737–3745.
Costa, E. J. X., R. H. Lopes, and M. T. Lamy-Freund. 1995. Permeability of pure lipid bilayers to melatonin. J. Pineal Res. 19:123–126.[Medline]
Cross, A. R., and O. T. G. Jones. 1986. The effect of the inhibitor diphenylene iodonium on the superoxide-generating system of neutrophils. Biochem. J. 237:111–116.[Medline]
Demaurex, N., G. P. Downey, T. K. Waddell, and S. Grinstein. 1996. Intracellular pH regulation during spreading of human neutrophils. J. Cell Biol. 133:1391–1402.
Dunford, H. B. 1999. Heme Peroxidases. John Wiley & Sons, New York.
Finocchiaro, L. M. E., E. S. Arzt, S. Fernandez-Castelo, M. Criscuolo, S. Finkielman, and V. E. Nahmod. 1988. Serotonin and melatonin synthesis in peripheral blood mononuclear cells: Stimulation by interferon- as part of an immunomodulatory pathway. J. Interferon Res. 8:705–716.[Medline]
Finocchiaro, L. M. E., V. E. Nahmod, and J. M. Launay. 1991. Melatonin biosynthesis and metabolism in peripheral-blood mononuclear leukocytes. Biochem. J. 280:727–731.[Medline]
Fujimoto, S., N. Kawakami, and N. Ohara. 1993. Formation of hydroxyl radical by the myeloperoxidase NADH oxygen system. Biol. Pharm. Bull. 16:525–528.[Medline]
Furtmüller, P. G., C. Obinger, Y. Hsuanyu, and H. B. Dunford. 2000. Mechanism of reaction of myeloperoxidase with hydrogen peroxide and chloride ion. Eur. J. Biochem. 267:5858–5864.[Medline]
Ginis, I., and A. I. Tauber. 1990. Activation mechanisms of adherent human neutrophils. Blood. 76:1233–1239.
Hampton, M. B., A. J. Kettle, and C. C. Winterbourn. 1998. Inside the neutrophil phagosome: oxidants, myeloperoxidase, and bacterial killing. Blood. 92:3007–3017.
Hauser, M. J. B., and L. F. Olsen. 1998. The role of naturally occurring phenols in inducing oscillations in the peroxidase-oxidase reaction. Biochemistry. 37:2458–2469.[Medline]
Henderson, L. M., and J. B. Chappell. 1996. NADPH oxidase of neutrophils. Biochim. Biophys. Acta. 1273:87–107.[Medline]
Henzler, T., and E. Steudle. 2000. Transport and metabolic degradation of hydrogen peroxide in Chara corallina: model calculation and measurements with the pressure probe suggest transport of H2O2 across water channels. J. Exp. Bot. 51:2053–2066.
Jones, R. D., J. T. Hancock, and A. H. Morice. 2000. NADPH oxidase: a universal oxygen sensor? Free Radic. Biol. Med. 29:416–424.[Medline]
Kindzelskii, A. L., M. J. Zhou, R. P. Haugland, and H. R. Petty. 1998. Oscillatory pericellular proteolysis and oxidant deposition during neutrophil migration. Biophys. J. 74:90–97.
Kindzelskii, A. L., and H. R. Petty. 2000. Extremely low frequency electric fields promote neutrophil extension, metabolic resonance and DNA damage during migration. Biochim. Biophys. Acta. 1495:90–111.[Medline]
Kirkor, E. S., and A. Scheeline. 2000. Nicotinamide adenine dinucleotide species in the horseradish peroxidase-oxidase oscillator. Eur. J. Biochem. 267:5014–5022.[Medline]
Klebanoff, S. J. 1991. Myeloperoxidase: occurrence and biological function. In Peroxidases in Chemistry and Biology. J. Everse, K. E. Everse and M. B. Grisham, editors. CRC Press, Boca Raton, FL. 1–35.
Kummer, U., M. J. B. Hauser, K. Wegmann, L. F. Olsen, and G. Baier. 1997. Oscillations and complex dynamics in the peroxidase-oxidase reaction induced by naturally occurring aromatic substrates. J. Am. Chem. Soc. 119:2084–2087.
Levy, R., O. Nagauker, E. Sikuler, T. L. Leto, and F. Schlaeffer. 1994. Elevated NADPH-oxidase activity in neutrophils from bile-duct-ligated rats; changes in the kinetic parameters and in the oxidase cytosolic factor p47. Biochim. Biophys. Acta. 1220:261–265.[Medline]
MacDonald, R. C., R. I. MacDonald, B. P. Menco, K. Takeshita, N. K. Subbarao, and L. R. Hu. 1991. Small-volume extrusion apparatus for preparation of large unilamellar vesicles. Biochim. Biophys. Acta. 1061:297–303.[Medline]
Matuszak, Z., K. Reszka, and C. F. Chignell. 1997. Reaction of melatonin and related indoles with hydroxyl radicals: EPR and spin trapping investigations. Free Rad. Biol. Med. 23:367–372.[Medline]
McCord, J. M. 2000. The evolution of free radicals and oxidative stress. Am. J. Med. 108:652–659.[Medline]
Mohan, R. F., and M. S. Jacobson. 1993. Inhibition of macrophage superoxide generation by dehydroepiandrosterone. Am. J. Med. Sci. 306:10–15.[Medline]
Morrey, K. M., J. A. McLachlan, C. D. Serkin, and O. Bakouche. 1994. Activation of human monocytes by the pineal hormone melatonin. J. Immunol. 153:2671–2680.[Abstract]
Nakamura, S., K. Yokota, and I. Yamazaki. 1969. Sustained oscillations in a lactoperoxidase, NADPH and O2-system. Nature. 222:794.[Medline]
Oehler, R., G. Weingartmann, N. Manhart, U. Salzer, M. Meissner, W. Schlegel, A. Spittler, M. Bergmann, D. Kandioler, S. Oismuller, H. M. Struse, and E. Roth. 2000. Polytrauma induces increased expression of pyruvate kinase in neutrophils. Blood. 95:1086–1092.
Olsen, L. F., A. Lunding, F. R. Lauritsen, and M. Allegra. 2001. Melatonin activates the peroxidase-oxidase reaction and promotes oscillations. Biochem. Biophys. Res. Commun. 284:1071–1076.[Medline]
Petty, H. R., R. G. Worth, and A. L. Kindzelskii. 2000. Imaging sustained dissipative patterns in the metabolism of individual living cells. Phys. Rev. Lett. 84:2754–2757.[Medline]
Petty, H. R. 2001. Neutrophil oscillations: temporal and spatiotemporal aspects of cell behavior. Immunol. Res. 23:85–94.[Medline]
Recchioni, R., F. Marcheselli, F. Moroni, R. Gaspar, S. Damjanovich, and C. Pieri. 1998. Melatonin increases the intensity of respiratory burst and prevents L-selectin shedding in human neutrophils in vitro. Biochem. Biophys. Res. Commun. 252:20–24.[Medline]
Rom, W. N., and T. Harkin. 1991. Dehydroepiandrosterone inhibits the spontaneous release of superoxide radical by alveolar macrophages in vitro in asbestosis. Environ. Res. 55:145–156.[Medline]
Rosenspire, A. J., A. L. Kindzelskii, and H. R. Petty. 2000. Interferon- and sinusoidal electric fields signal by modulating NAD(P)H oscillations in polarized neutrophils. Biophys. J. 79:3001–3008.
Scheeline, A., D. L. Olson, E. P. Williksen, G. A. Horras, M. L. Klein, and R. Larter. 1997. The peroxidase-oxidase oscillator and its constituent chemistries. Chem. Rev. 97:739–756.[Medline]
Sevcik, P., and H. B. Dunford. 1991. Kinetics of the oxidation of NADH by methylene blue in a closed system. J. Phys. Chem. 95:2411–2415.
Silva, S. O., V. F. Ximenes, L. H. Catalani, and A. Campa. 2000. Myeloperoxidase-catalyzed oxidation of melatonin by activated neutrophils. Biochem. Biophys. Res. Commun. 279:657–662.[Medline]
Sørensen, O., and N. Borregaard. 1999. Methods for quantitation of human neutrophil proteins: a survey. J. Immunol. Methods. 232:179–190.[Medline]
Tauber, A. I. 1987. Protein kinase C and the activation of the human neutrophil NADPH oxidase. Blood. 69:711–720.
Tyagi, S. R., M. Tamura, D. N. Burnham, and J. D. Lambeth. 1988. Phorbol myristate acetate (PMA) augments chemoattractant-induced diglyceride generation in human neutrophils but inhibits phosphoinositide hydrolysis. Implications for the mechanism of PMA priming of the respiratory burst. J. Biol. Chem. 263:13191–13198.
Whitcomb, J. M., and A. G. Schwartz. 1985. Dehydroepiandrosterone and 16 alpha-Br-epiandrosterone inhibit 12–0-tetracanolyphorbol-13-acetate stimulation of superoxide radical production by human polymorphonuclear leukocytes. Carcinogen. 6:333–335.
White, A., P. Handler, E. L. Smith, R. L. Hill, and I. R. Lehman. 1978. Principles of Biochemistry. McGraw-Hill, New York.
Ximenes, V. F., L. H. Catalani, and A. Campa. 2001. Oxidation of melatonin and tryptophan by an HRP cycle involving compound III. Biochem. Biophys. Res. Commun. 298:130–134.
Zhang, H., G. Squadrito, and W. A. Pryor. 1998. The reaction of melatonin with peroxynitrite: formation of melatonin radical cation and absence of stable nitrated products. Biochem. Biophys. Res. Commun. 251:83–87.[Medline]
This article has been cited by other articles:
 |
|
 |
|
  U. Kummer, J. Zobeley, J. C. Brasen, R. Fahmy, A. L. Kindzelskii, A. R. Petty, A. J. Clark, and H. R. Petty Elevated Glucose Concentrations Promote Receptor-Independent Activation of Adherent Human Neutrophils: An Experimental and Computational Approach Biophys. J., April 1, 2007; 92(7): 2597 - 2607. [Abstract] [Full Text] [PDF]
|
|
