![]() |
|
|
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
University at Buffalo, State University of New York, Department of Pharmacology and Toxicology, Buffalo, New York
(Received August 9, 2006; accepted October 25, 2006)
| Abstract |
|---|
|
|
|---|
Activation of OPs has been attributed to the cytochrome P450 (P450) family of enzymes. P450s have also been shown to carry out direct detoxification of the parent compounds through a dearylation reaction (Mutch et al., 2003
; Poet et al., 2003
). Humans, with their ability to bioactivate OPs, primarily through CYP1A2, 2B6, 2C19, and 3A4, are particularly sensitive to the actions of OPs (Mutch et al., 1999
; Sams et al., 2000
; Buratti et al., 2003
).
Parathion and chlorpyrifos often have been used as model compounds, the research base for parathion spanning decades (for review see Knaak et al., 2004
). Although parathion use is banned within the United States, it is still used worldwide because of its potency, effectiveness, and low cost. Chlorpyrifos is available for use in the United States and is found in formulations carrying brand names such as Dursban and Lorsban; it is often the insecticide of choice when it comes to cockroach control. Paraoxon and chlorpyrifos-oxon represent the activated forms of parathion and chlorpyrifos, respectively, whereas O,O-diethyl phosphate and p-nitrophenol (PNP), from parathion, or O,O-diethyl phosphate and 3,4,5-trichloropyrindinol (TCP), from chlorpyrifos, represent the more readily cleared detoxification products. The balance between activation and detoxification of OPs determines their risk to humans.
Physiologically based pharmacokinetic/pharmacodynamic (PBPK/PD) models attempt to better quantify exposure, transportation, activation, detoxification, and clearance of OPs to determine safe levels of human exposure (Timchalk et al., 2002
; Knaak et al., 2004
). PBPK/PD models are dependent on the available kinetic parameters for metabolism, which are often not available and/or consistent (for review, see Knaak et al. (2004
).
PBPK/PD models, which use previously published kinetic data for OP metabolism generated from mouse, rat, or human liver microsomes have proven to be unstable, nonreproducible, and under-represent interindividual variability (Knaak et al., 2004
). To circumvent this issue, modifying PBPK/PD models to use individual P450 activities combined with their respective hepatic P450 content may prove to be more reproducible and stable. To have reliable modeling, accurate and complete kinetic data need to be generated. Kinetic data (Km and Vmax) on the metabolism of OPs by specific P450s along with P450-specific content (pmol of P450/mg of microsomal protein) will generate a PBPK/PD model that can better adjust for age, sex, genetic polymorphisms, or other factors, which may alter P450 content and activity, which, in turn, contribute directly to the risk OPs pose to individuals.
We report in this study the identification of specific human P450s that mediate parathion and chlorpyrifos metabolism, as well as report, for the first time, their respective Km and Vmax values for activation and detoxification. These P450-specific kinetic parameters can then be linked to specific hepatic P450 content values as a function of age for future inclusion in PBPK/PD models.
| Materials and Methods |
|---|
|
|
|---|
Experimental Conditions. OP stock solutions were prepared in methanol and stored at 20°C when not in use. Incubations with either human liver microsomes (0.5 mg of protein/ml final concentration) or recombinant P450s (0.030.06 nmol P450 final concentration) were carried out in buffer (100 mM Tris-HCl, and 5 mM MgCl2, 1 mM EDTA, and 50 µM iso-OMPA, pH 7.4 at 37°C) in a final volume of 0.5 ml. Reactions were initiated with the addition of 1 mM NADPH and incubated for 2 min at 37°C. EDTA was included to inhibit A-esterases, whereas iso-OMPA inhibited B-esterases (Reiner et al., 1993
). The reaction was quenched with 1 volume of ice-cold methanol with 0.1% phosphoric acid and centrifuged; then, the supernatant was transferred to HPLC vials and capped for analysis. All reactions were replicated (n = 34) for statistical purposes.
Metabolite Detection. OPs and their respective metabolites were analyzed by reverse-phase HPLC (C18, 5-µm particle size, 25 cm x 4.6 mm i.d; Supelco, Bellefonte, PA) using a Hewlett Packard (Palo Alto, CA) model 1100 high-performance liquid chromatograph with a model 1046A diode array detector. Methanol (buffer A) and 94.99% water/5% acetonitrile/0.01% phosphoric acid (buffer B) were used in a gradient elution consisting of 40% buffer A/60% buffer B to 100% buffer A over 30 min at a 1 ml/min flow rate. Chemical detection was determined at the UV wavelengths of 275 nm for parathion and paraoxon, 320 nm for PNP, 290 nm for chlorpyrifos and chlorpyrifos-oxon, 300 nm for TCP, and 320 nm for diethylthiophosphate. The minimum level of detection of these compounds was 2.5 ng.
Kinetic Calculations. Vmax and Km were determined by nonlinear regression analysis [enzyme kinetics module of SigmaPlot V9.01 (SyStat Software Inc., Point Richmond, CA)] of hyperbolic plots (i.e., velocity versus [S]) obeying Michaelis-Menten kinetics. OP parent compound concentration (micromolar) was set as the independent variable, whereas rate of product formation (pmol/min/nmol P450 or pmol/min/mg protein) was the dependent variable. The kinetic module handles multiple data sets and determines average values along with standard deviations from multiple experiments.
| Results |
|---|
|
|
|---|
Since the use of pooled human liver microsomes did not reflect the interindividual variability in OP metabolism, multiple single-donor hepatic microsomal specimens were assayed. As expected, the rate and amount of OP product formation varied for each donor. Although multiple single donors were assayed, one sample, SD101, demonstrated Km and Vmax values for paraoxon and PNP formation of 8.7 µM, 731 pmol/min/mg protein and 15.2 µM, 630 pmol/min/mg protein, respectively. Km and Vmax values for chlorpyrifos-oxon and TCP formation were 4.5 µM, 611 pmol/min/mg protein and 6.6 µM, 705 pmol/min/mg protein, respectively. In contrast, for SD112, Km and Vmax values for paraoxon and PNP formation were 39.1 µM, 850 pmol/min/mg protein and 55.9 µM, 683 pmol/min/mg protein, respectively. Km and Vmax values for chlorpyrifos-oxon and TCP formation in this specimen were 24.7 µM, 490 pmol/min/mg protein and 27.1 µM, 402 pmol/min/mg protein, respectively. Gentest reported the P450-specific enzyme activities for SD101 as 2.2, 0.34, 0.37, and 4.2 nmol/min/mg protein for CYP1A2, 2B6, 2C19, and 3A4, respectively. In contrast, P450-specific enzyme activities for SD112 were 0.89, 0.11, 0.059, and 8.9 nmol/min/mg protein for CYP1A2, 2B6, 2C19, and 3A4, respectively.
Studies with individual recombinant human P450s were also conducted to identify specific P450s that biotransform these OPs and generate P450-specific kinetic data to predict the P450s that will play a prominent role in vivo. CYP1A2, 2B6, 2C9, 2C19, 3A4, 3A5, and 3A7 were found to metabolize parathion, whereas no detectable metabolites were produced with CYP1A1, 2D6, or 2E1. Figure 1 illustrates Michaelis-Menten plots for three of the most metabolically active P450s (CYP2B6, 2C19, and 3A4) for chlorpyrifos, although CYP1A2, 3A5, and 3A7 were also found to be capable of limited metabolism of chlorpyrifos. The Hill equation was also used for determining the kinetics of CYP3A4 in Fig. 1 (plot not shown) because more than 1 mol of substrate was bound to the enzyme (n > 1); however, both equations resulted in the same Km and Vmax values.
|
|
| Discussion |
|---|
|
|
|---|
Although several studies have assessed the in vitro metabolism of OPs by human liver microsomes, the kinetic values may vary widely because of the marked variability of P450 content and activity in procured human specimens and differences in incubation conditions and analytical methodology (for review, see Knaak et al., 2004
). Differences in experimental methods include detection of only oxon metabolites, the use of excessive substrate concentrations (supersaturated conditions), and nonphysiological pH ranges. Complications and extensive variability in the procurement of human liver specimens for in vitro studies are a major concern. Often it is not possible to limit warm-ischemic time in liver tissue obtained from organ donors. It is necessary to limit warm-ischemic time before proper cryostorage to have accurate measures of enzyme activity. Furthermore, exposure to drugs and other substances that modulate enzyme activity is often not controlled for in specimens obtained from tissue donors. Thus, there remains considerable uncertainty regarding whether variability in metabolism kinetics reflects true age and genetic variability or is an artifact of human liver specimen procurement and storage. The metabolism of selected OPs has also been assessed using specific recombinant human P450s; however, these earlier studies did not report kinetic parameters (Km, Vmax) for all reactions (Sams et al., 2000
, 2004
; Tang et al., 2001
; Mutch et al., 2003
).
In the present study, pooled human liver microsomes demonstrated the ability to metabolize both parathion and chlorpyrifos. Kinetic parameters for the two model OPs varied, as did the ratio of dearylation/desulfation products. Although the results represent an average, the Km and Vmax values using nine substrate concentrations, the results from pooled microsomes did not estimate the interindividual variability that occurs within a population.
Interindividual variability in hepatic metabolism was illustrated in reactions conducted with characterized single-donor human liver microsomes. Marked differences in Km and Vmax values for OP metabolism exist and may be due to variable levels of various P450s within a given microsomal sample from a given donor liver. The higher Km values for specimen SD112 and the high Km values for CYP3A4 (Table 1) suggest that CYP3A4 is playing a primary role in the metabolism of the model OPs by this sample of human liver microsomes. The relative activities of CYP1A2, CYP2B6, and CYP2C19 were lower, whereas CYP3A4 activity was higher in SD112 compared with SD101, supporting the dominant role of CYP3A4 in the metabolism of the two model OPs in SD112.
To better address the limitations of using human liver microsomes for establishing kinetic values, OP metabolism was investigated with recombinant human P450s (summarized in Table 1). As expected, since each P450 isozyme has substrate binding sites that differ, each enzyme demonstrated different Km and Vmax values along with different ratios of product formation. Although CYP3A4 showed the greatest Vmax value for parathion, this enzyme's lack of specificity resulted in the highest Km value, which minimizes the role of CYP3A4 in metabolism at lower OP exposures. CYP2B6 and 2C19 have the lowest Km values for parathion and chlorpyrifos, supporting the major role these forms play in metabolism at low-level real-world exposures.
Due to the low Km values of CYP2B6, its activity is of interest. CYP2B6 activity can be expressed through its Vmax/Km values of 7.87, 2.45, 15.56, and 0.74 for the formation of paraoxon, PNP, chlorpyrifos-oxon, and TCP, respectively (Table 1). Thus, CYP2B6 has greater activity in the formation of the oxon and, thus, may be more important in assessing risk because it preferentially forms the toxic compound over the direct detoxification of the OPs at low level exposures. Conversely, CYP2C19 has the highest Vmax/Km in the formation of PNP and TCP, and thus plays a prominent role in the dearylation (detoxification) reaction.
Current PBPK/PD models for OPs utilize kinetic values from rat liver microsomal metabolism studies that do not reflect human enzymes. The building of kinetic models that use human P450-specific kinetic parameters and specific hepatic P450 content should prove to be more accurate and more easily modified to address factors such as gender, age, or polymorphisms, which may affect P450 protein expression and activity. Human hepatic P450 levels are known to vary across different age groups (Tateishi et al., 1997
). In addition, polymorphism in CYP2B6 and 2C19 is known to alter protein expression and/or activity, and, thus, these P450s could potentially serve as biomarkers of susceptibility to these agents. In current PBPK/PD models, which vary parameters based on body weight, infants are treated as small adults, and these models may not be accurate or protective of the most sensitive segment of the population. As OPs are the most widely used class of pesticides in the world, safe levels of exposure need to be set to protect the most sensitive individuals. Determining human P450-specific kinetic parameters for OP metabolism is the first step in constructing more refined models, which will include age-dependent P450 content to better assess risk associated with OP exposure in infants, children, or other sensitive subgroups.
| Footnotes |
|---|
Part of this work was presented at the 2006 Society of Toxicology Annual Meeting, San Diego, CA [The Toxicologist (supplement to Toxicological Sciences) 90:S-1, March 2006].
Article, publication date, and citation information can be found at http://dmd.aspetjournals.org.
ABBREVIATIONS: OP, organophosphorus pesticide; P450, cytochrome P450; PBPK/PD, physiologically based pharmacokinetic/pharmacodynamic; TCP, 3,4,5-trichloropyrindinol; PNP, p-nitrophenol; iso-OMPA, tetraisopropyl pyrophosphoramide; HPLC, high-performance liquid chromatography.
Address correspondence to: Dr. James R. Olson, 102 Farber Hall, 3435 Main St., Buffalo, NY 14214. E-mail: jolson{at}buffalo.edu
| References |
|---|
|
|
|---|
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||