|Home | About | Journals | Submit | Contact Us | Français|
Harmaline is a β-carboline alkaloid showing neuroprotective and neurotoxic properties. Our recent studies have revealed an important role for cytochrome P450 2D6 (CYP2D6) in harmaline O-demethylation. This study, therefore, aimed to delineate the effects of CYP2D6 phenotype/genotype on harmaline metabolism, pharmacokinetics (PK) and pharmacodynamics (PD), and to develop a pharmacogenetics mechanism-based compartmental PK model. In vitro kinetic studies on metabolite formation in human CYP2D6 extensive metabolizer (EM) and poor metabolizer (PM) hepatocytes indicated that harmaline O-demethylase activity (Vmax/Km) was about 9-fold higher in EM hepatocytes. Substrate depletion showed mono-exponential decay trait, and estimated in vitro harmaline clearance (CLint, μL/min/106 cells) was significantly lower in PM hepatocytes (28.5) than EM hepatocytes (71.1). In vivo studies in CYP2D6-humanized and wild-type mouse models showed that wild-type mice were subjected to higher and longer exposure to harmaline (5 and 15 mg/kg; i.v. and i.p.), and more severe hypothermic responses. The PK/PD data were nicely described by our pharmacogenetics-based PK model involving the clearance of drug by CYP2D6 (CLCYP2D6) and other mechanisms (CLother), and an indirect response PD model, respectively. Wild-type mice were also more sensitive to harmaline in marble-burying tests, as manifested by significantly lower ED50 and steeper Hill slope. These findings suggest that distinct CYP2D6 status may cause considerable variations in harmaline metabolism, PK and PD. In addition, the pharmacogenetics-based PK model may be extended to define PK difference caused by other polymorphic drug-metabolizing enzyme in different populations.
Cytochrome P450 2D6 (CYP2D6) is one of the most important polymorphic phase I drug-metabolizing enzymes that is involved in the biotransformation of 20–30% of marketed drugs and some endogenous substrates [1, 2]. Individuals with impaired or deficient CYP2D6 activity could be more sensitive to CYP2D6-metabolized drugs, which is likely dependent on the fraction of CYP2D6-mediated metabolism in total drug clearance and the range between therapeutic and toxic doses of the drug. CYP2D6 also metabolizes a variety of drugs of abuse such as amphetamines [3–6], indoleamines [7–9] and many other designer drugs [10, 11]. Delineation of the impact of CYP2D6 on the metabolism, pharmacokinetics (PK) and pharmacodynamics (PD) of these substances would provide increased understanding of individual vulnerability to and/or protection from illicit drugs of abuse .
Harmaline, 3,4-dihydro-7-methoxy-1-methyl-β-carboline, is a major form of alkaloid in the seeds of Peganum harmala (Syria rue) that has been used in traditional medicine and for recreational purpose . Harmaline is also a typical psychotropic ingredient in the recreational beverage Ayahuasca, which shows potential in treating many psychiatric disorders , whereas the therapeutic potential and medical uses are hampered by many safety issues including abuse liability. Harmaline readily modulates the levels of amine neurotransmitters via inhibition of monoamine oxidase. It may directly bind to specific receptors including 5-HT, benzodiazepine and imidazoline receptors. With diverse biochemical properties, harmaline has shown neuroprotective, anxiolytic, antinociceptive, antioxidant and vasorelaxant properties in animal models [15–19]. On the other hand, harmaline alone or along with other tryptamine substances may result in severe or even fatal toxicities in animals [20–22] and in humans [23–25]. Of particular note, harmaline acts on the olivo-cerebellar system, causes rhythmic activity, and leads to postural tremors in animal models. Harmaline-induced tremor has many similar characteristics as essential tremor, a common movement disorders in humans. Therefore, harmaline tremor has been widely used as a model to understand the pathogenesis of essential tremor and to develop therapeutic agents for generalized tremor [21, 22, 26].
Following an acute administration to rats, harmaline is rapidly distributed into different organs and mainly metabolized via O-demethylation . Recently we have shown that harmaline O-demethylation is primarily catalyzed by cytochrome P450 2D6 (CYP2D6) and CYP1A2 in human liver microsomes . This study, therefore, aimed to 1) define the effects of CYP2D6 status on harmaline metabolism, pharmacokinetics (PK) and pharmacodynamics (PD), and 2) use harmaline as a model drug to develop pharmacogenetics-based PK compartmental model to provide increased understanding of the PK difference caused by CYP2D6 metabolism.
Harmaline hydrochloride dihydrate, harmalol hydrochloride dihydrate and β-glucuronidase were purchased from Sigma-Aldrich (St. Louis, MO). Cryopreserved human hepatocytes, InVitroGRO HT Medium, KHB buffer and Torpedo Antibiotic Mix were purchased from Celsis (Chicago, IL) or BD Biosciences (San Jose, CA). The lot numbers of CYP2D6 PM hepatocyte samples were VTA (dextromethorphan O-demethylase activity, 1 pmol/min/106 cells; CYP2D6 genotype, *3/*4), ETR (1 pmol/min/106 cells; *4/*4), PFM (1 pmol/min/106 cells; *1/*7), and HH183 (bufuralol 1′-hydroxylase activity, not detectable), and EM hepatocyte samples were GNG (17 pmol/min/106 cells; *1/*4), KRM (36 pmol/min/106 cells; *1/*4), KRJ (40 pmol/min/106 cells; *1/*4), VEN (69 pmol/min/106 cells; *1/*1), AGR (96 pmol/min/106 cells; genotyped not determined), and HH169 (bufuralol 1′-hydroxylase activity, 60 pmol/min/106 cells). Dulbecco’s Modified Eagle’s Medium (DMEM) was bought from Mediatech, Inc (Manassas, VA). All other reagents or organic solvents used were either analytical or high-performance liquid chromatography (HPLC) grade.
Age-matched (7–10 weeks old) male wild-type FVB/N and CYP2D6-humanized (Tg-CYP2D6) mice  weighing 25–32 g were used in the study. Animals were housed in an animal care facility maintained at 20 ± 2°C on a 12-h light/dark cycle. Food and water were provided ad libitum. All animal procedures were approved by the Institutional Animal Care and Use Committee at University at Buffalo, The State University of New York.
Cryopreserved human hepatocytes were thawed in a 37°C water bath, and suspended in pre-warmed suspension media (DMEM or InVitroGRO HT Medium supplemented with Torpedo Antibiotic Mix, 49/1, v/v). Cell suspension was mixed gently, centrifuged at 500 rpm for 5 min at 4°C. After the supernatant was removed, cells were resuspended in incubation media, and cell viability was determined by the Trypan Blue exclusion test.
Reaction was initiated by mixing the same volume of drug solution (in KHB buffer) and resuspended hepatocytes, and the mixture was incubated at 37°C, 5% carbon dioxide and saturating humidity. In metabolite production kinetic studies, incubations were conducted with 96-well plate, and all reactions were carried out in triplicate. Each well consisted of 0.25×106 cells/mL of hepatocytes and desired concentrations of harmaline (0.5–150 μM) in a total volume of 200 μL. Reaction mixtures were snap-frozen in liquid nitrogen after 30 min incubation. In substrate depletion kinetic studies, incubations were carried out using 6-well plate and all reactions were performed in duplicate. Each well consisted of 2 μM of harmaline and 0.5×106 cells/mL of hepatocytes in a total volume of 1.6 mL. Forty μL of media were aliquoted at different time point (0, 15, 30, 45, 60 and 90 min) and snap-frozen in liquid nitrogen.
Adult male Tg-CYP2D6 and wild-type mice were randomized to receive either 5 or 15 mg/kg of harmaline intravenously (i.v.) or intraperitoneally (i.p.). Blood samples were collected from retro-orbital plexus at different time points (0–120 min, N = 3 per time point) after drug administration. Serum was prepared with a serum separator (BD Biosciences, San Jose, CA), and stored at −80°C until analysis.
Upon analyses, hepatocyte samples were thawed on ice, and cells were removed by centrifugation. Forty μL of sample were digested with β-glucuronidase (80 U in 40 μL of 0.2% sodium chloride solution) in 0.2 M sodium acetate buffer (pH 4.75, 80 μL) at 37°C for 3 h. The mixture was then deproteinized with 5 μL of ice-cold perchloric acid (70%) and centrifuged at 14,000 g for 10 min. The supernatant was transferred into new vials and injected for HPLC analysis. Likewise, serum sample (40 μL) was digested with β-glucuronidase, and the mixture was processed by liquid-liquid extraction with 2 mL of methyl t-butyl ether/ethyl acetate (v/v, 1:2). The organic layer was transferred into new tubes, evaporated to dryness, reconstituted with 40 μL of 50% methanol containing 0.1% formic acid, and then analyzed by tandem mass spectrometry (MS/MS).
HPLC quantification of harmaline and harmalol in hepatocyte samples was performed as reported  on an Agilent 1100 series HPLC system consisting of the online vacuum degasser, quaternary pump, autosampler, thermostat controlled column compartment, fluorescence detector and diode-array detector. The calibration curve was linear from 0.2 to 500 pmol for both harmaline and harmalol. Assay accuracy was 95–110%, and intra-day and inter-day variations were both less than 10%.
The LC-MS/MS system used for the analyses of serum samples consisted of a Shimadzu prominence HPLC (Kyoto, Japan) coupled to an API 3000 turbo ionspray ionization triple-quadrupole mass spectrometer (Applied Biosystems, Foster City, CA). A Phenomenex phenyl-hexyl column (50 × 4.6 mm, 3μm; Torrance, CA) was used for the analyses. A gradient elution was employed to separate the analytes, which included an initial 1-min isocratic elution with 75% of Buffer A (0.1% formic acid in water) and 25% of Buffer B (0.1% formic acid in methanol), an increase of Buffer B to 90% at 1.1 min and then elution with 90% of Buffer B for 3 min, followed by the initial condition from 4.1 min. The mass spectrometer was operated in positive mode with transitions m/z 215.2→174.2, 201.2→160.2, and 203.2→158.3 for harmaline, harmalol, and 5-methyl-N,N-dimethyltryptamine (internal standard), respectively. An online motorized six-port divert valve was used to introduce the LC eluent to the mass spectrometer over the period of 2.5–6 min for data acquisition, whereas eluent of 0–2.5 min was diverted to the waste. Calibration curve was linear from 12.3 to 27,000 nmol/L for harmaline, and linear from 37 to 27,000 nmol/L for harmalol, respectively. Assay accuracy was 91–106%, and intra-day and inter-day variations were both less than 10%.
After i.p. administration of harmaline (5 or 15 mg/kg) or vehicle control, mouse body temperature was monitored at different time points (0, 5, 10, 15, 20, 30, 45, 60, 90, 120, 150 and 180 min) using a High Precision Digital Thermometer with a microprobe (Barnstead International, Dubuque, IA) that was inserted about 1.5 cm into the rectum .
Thirty minutes after i.p. administration of harmaline (0–10 mg/kg) or drug vehicle, mouse was placed in a 28 × 17 × 16 cm transparent plastic cage with 5 cm layer of sawdust (Harlan Sani-Chips, Harlan-Teklad, Indianapolis, IN) on the floor and 20 glass marbles (1.5 cm diameter) distributed evenly against the wall of cage. Mouse was removed from the cage after 30 min, and the number of marbles buried (at least 2/3 in the sawdust) was recorded and normalized with vehicle control.
All values were expressed as mean ± SD when experiments were carried out using different samples, or mean ± SEM when using the same sample. Enzyme kinetic parameters were estimated by nonlinear regression using GraphPad Prism 5 (San Diego, CA). Specifically, metabolite formation data were fit to one-enzyme or two-enzyme Michaelis-Menten equation, and intrinsic clearance (CLint) was calculated as Vmax/Km. Substrate depletion data were fit to mono-exponential decay model, C(t) = C0·e−kt, where C0 is the initial drug concentration in the incubation, and k is the terminal depletion rate constant. In vitro half-life (T1/2) was calculated by dividing 0.693 by the rate constant. In vitro clearance (CLint) was calculated by dividing initial dose (D0) by the area under the concentration-time curve (AUC), and normalized with cell variability.
PK/PD modeling was conducted using WinNonLin (version 5.0, Pharsight, Mountain View, CA). Serum drug concentration data were first modeled by non-compartmental analysis, essentially the same as described . A pharmacogenetics-based PK model (Fig. 1) was then constructed to illustrate harmaline PK difference in Tg-CYP2D6 and wild-type mice. The initial estimates were obtained from non-compartmental analysis. Parameters describing the pharmacogenetics-based PK model include the clearance contributed by CYP2D6-mediated metabolism (CLCYP2D6) and clearance contributed by other elimination mechanisms (CLother), volume of distribution of the central compartment (VC) and peripheral compartment (VT), absorption rate constant (Ka), distribution clearance between central and peripheral compartments (CLD), and bioavailability (F). Note that harmaline is eliminated by CLCYP2D6 and CLother in Tg-CYP2D6 mice (or EMs), i.e. CL = CLCYP2D6 + CLother, whereas the drug is eliminated by CLother only in wild-type control mice (or PMs), i.e. CL = CLother. In contrast, drug absorption (Ka) and distribution (CLD) properties are the same in the two genotyped subjects. Once the pharmacogenetics-based PK model was finalized, the associated PK parameters were set as constants and linked to the PD model (Fig. 1), in which the time course of hypothermic response was sequentially analyzed by the basic indirect response (IDR) model IV . Change of mouse body temperature (T) is due to heat production (kin) and loss (kout) nature of temperature control.
The differential equations derived from the final PK/PD model (Fig. 1) were as follows:
Equations (1), (2) and (6) were used to fit the i.p. data, and (4) and (6) fit the i.v. data of Tg-CYP2D6 mice, respectively. Likewise, equations (1), (3) and (6) were used to fit the i.p. data, and (5) and (6) fit the i.v. data of wild-type mice, respectively. Equation (7) was used to fit mouse body temperature data, where Cc is drug concentration, Smax reflects the maximum fractional drug effect, SC50 is the drug concentration producing 50% of the maximum effect, and kin can be defined as the initial condition of body temperature (T0) times kout. The appropriate PK/PD model was selected based on the goodness-of-fit criteria including coefficient variation of the estimates (CV%), R-square, Akaike’s Information Criterion and Schwarz Criterion.
Mouse marble-burying data were fit to the empirical dose-response model (Normalized response = 100/(1+10^((LogIC50-LogDose)*Hill Slope)) (GraphPad Prism 5).
Depending on the number of groups and variances, data were analyzed by unpaired Student’s t-test, one-way ANOVA analysis followed by Dunnett’s test, or two-way ANOVA analysis (GraphPad Prism 5). Significant difference was considered if the probability was less than 0.05 (P < 0.05).
To delineate the effects of CYP2D6 status on harmaline metabolism, we first compared harmalol production between CYP2D6 EM and PM hepatocytes in vitro. Interestingly, all EM hepatocytes (N = 5) showed monophasic Michaelis-Menten enzyme kinetics, whereas PM hepatocytes (N = 4) exhibited biphasic kinetics (Fig. 2). Estimated enzyme efficiency (Vmax/Km, μL/min/106 cells) was about 9-fold lower in CYP2D6 PM hepatocytes (0.71 ± 0.44) than that in CYP2D6 EM hepatocytes (6.32 ± 2.35) (Table 1). The results suggest that PMs with CYP2D6 deficiency may have a lower capacity to catalyze harmaline O-demethylation metabolism.
We then investigated the impact of CYP2D6 status on overall harmaline depletion in human hepatocytes. The data showed that harmaline depletion was significantly slower in CYP2D6 PM hepatocytes (N = 4) than EM hepatocytes (N = 6) (Fig. 3A). Consistently, harmalol levels were much higher in CYP2D6 EM hepatocytes (Fig. 3B). Kinetic parameters were estimated for harmaline depletion using mono-exponential decay model. Harmaline in vitro half-life (T1/2, min) was much longer in CYP2D6 PM hepatocytes (111 ± 28) than EM hepatocytes (46.1 ± 15.8), which was consistent with a significantly lower in vitro intrinsic clearance (CLint, μL/min/106 cells) in PM hepatocytes (28.5 ± 6.2) than that in EM hepatocytes (71.1 ± 23.1) (Table 2). The results suggest that CYP2D6 status could be an important factor controlling harmaline metabolic elimination.
We thus employed Tg-CYP2D6 and wild-type mouse models to examine the effects of CYP2D6 status on harmaline PK. Serum harmaline concentration versus time courses were monitored (Fig. 4) in the two genotyped mice after administration of harmaline (5 and 15 mg/kg; i.v. and i.p.). Non-compartmental analyses of harmaline PK profiles showed that wild-type mice lacking of CYP2D6 had higher exposure (AUC) to harmaline (Table 3). For example, AUC values were about 2-fold higher in wild-type mice than Tg-CYP2D6 mice dosed i.p. with harmaline. Consistently, the elimination of harmaline was slower in wild-type mice, as indicated by a lower systemic clearance (CL) and longer half-life (T1/2) (Table 3). Additionally, wild-type mice had significantly lower blood metabolite (harmalol) concentrations than Tg-CYP2D6 mice (data not shown), indicating the difference in harmaline O-demethylation capacity between the two genotyped mice in vivo. The results suggest that transgenic mice expressing functional CYP2D6 have enhanced capacity to eliminate harmaline, leading to lower and shorter exposure to the drug.
Because Tg-CYP2D6 and wild-type mice differ only in CYP2D6-mediated drug metabolism [1, 29, 33], we further built a new pharmacogenetics mechanism-based compartmental PK model (Fig. 1) to elucidate harmaline PK difference between the two groups of mice. The mechanistic difference in drug clearance between Tg-CYP2D6 and wild-type mice (or human CYP2D6 EM and PM subjects) is that the drug is eliminated via CLCYP2D6 and CLother in Tg-CYP2D6 mice, whereas only via CLother in wild-type mice. This model nicely described harmaline PK data in the two genotyped mice, and the results indicated that harmaline clearance by CLCYP2D6 and CLother was 60.8 and 96.2 mL/min/kg, respectively (Table 4). In turn, this pharmacogenetics-based PK model can be applied to understand potential difference in PK of CYP2D6-metabolized drugs between human CYP2D6 PM and EM subjects.
To investigate if the difference in harmaline PK could be translated into significant difference in drug response, we compared harmaline-induced hypothermic effects between wild-type and Tg-CYP2D6 mice. Administration of harmaline resulted in a rapid decrease in body temperature in both wild-type and Tg-CYP2D6 mice (Fig. 5), whereas injection of drug vehicle had no significant effects within the period (data not shown). Importantly, the hypothermia induced by 15 mg/kg of harmaline was more severe in wild-type mice than Tg-CYP2D6 mice. In particular, wild-type mice recovered from hypothermia at a lower rate, which was associated with lower drug clearance and higher exposure to harmaline (Fig. 4 and Table 3). The body temperature data were fit to the indirect response PD model with stimulation of kout (IDRIV), which possibly involves the serotonergic action of harmaline . As expected, there was no difference in PD parameters (Table 5) estimated for the two genotyped mice, suggesting that the observed difference in hypothermia was caused by the differences in harmaline PK between wild-type and Tg-CYP2D6 mice.
In addition, we employed marble-burying test [34–36] to evaluate the difference in harmaline drug effects between wild-type and Tg-CYP2D6 mice. Of note, the behavior of burying harmless marbles by mice is considered as an animal model of obsessive-compulsive disorder, and may indicate the anxiolytic activity of a drug. The results showed that Tg-CYP2D6 mice dosed i.p. with 5.0 or 7.5 mg/kg of harmaline buried a significantly greater number of marbles than wild-type mice treated with the same dose of harmaline (Fig. 6). Overall, harmaline dose responses in the two genotyped mice were significantly different. Particularly, wild-type mice had a significantly lower ED50 value (5.00 ± 1.04 mg/kg) and a steeper Hill slope (−5.79 ± 1.58) than Tg-CYP2D6 mice (ED50 6.75 ± 1.05 mg/kg and Hill slope −4.28 ± 0.84, respectively). These findings suggest that wild-type mice lacking of CYP2D6 are more sensitive to harmaline.
Humans are exposed to endogenously-synthesized, dietary or recreational β-carbolines, an important group of indoleamines that have a wide spectrum of biochemical, pharmacological and toxicological properties. Harmaline represents a β-carboline alkaloid that shows both neurotoxic and neuroprotective effects, and cases of harmaline-related intoxication have been documented in recent years [23–25]. Our previous study has revealed an important contribution of polymorphic CYP2D6 to harmaline O-demethylation metabolism , which is the predominant pathway in harmaline clearance in rat models . The data from current study, using human hepatocytes and transgenic mouse models, clearly demonstrated that CYP2D6 status had significant impact on harmaline total clearance and drug responses, suggesting that CYP2D6 polymorphism could be an additional risk factor in harmaline intoxication.
CYP2D6 genetic polymorphism is an important factor causing inter-patient variations in PK, and is well recognized for its clinical significance . Defect of CYP2D6 gene accurately predicts CYP2D6 PM phenotype, in comparison to EM carrying functional CYP2D6 alleles. Of note, 7–10% of Caucasians are CYP2D6 PMs who, taking standard or higher dose of CYP2D6-inactivated drug, may have higher blood drug concentrations. While harmaline is mainly eliminated via O-demethylation metabolism  and CYP2D6 significantly contributes to this biotransformation , CYP2D6 deficiency led to the impaired harmaline O-demethylase capacity (Fig. 2 and Table 1) and slower substrate depletion (Fig. 3 and Table 2). The less dramatic difference in substrate depletion (2.5-fold) than metabolite formation (9-fold) between EM and PM hepatocytes may reflect the actual contribution of CYP2D6-mediated metabolism to total harmaline depletion. In comparison to Tg-CYP2D6 mice, wild-type mice lacking of CYP2D6 activity indeed had higher and longer exposure to harmaline (Fig. 4 and Table 3). Note that the Tg-CYP2D6 and wild-type mice may be valuable animal models for human CYP2D6 EMs and PMs, respectively [29, 33], in assessing the significance of CYP2D6-mediated drug metabolism at the systemic level. Our findings suggest that CYP2D6 PMs would have lower harmaline hepatic clearance, leading to potential risk of higher exposure to the drug. It is also documented that people taking harmine, another β-carboline substance metabolized by CYP2D6, CYP1A and CYP2C , show distinct slow and fast metabolizer phenotypes . Together, the presence of metabolic pharmacogenetics may raise additional concerns for use of these psychotropic β-carboline xenobiotics.
Assessing the risk of polymorphic metabolism for a drug is important not only in clinical pharmacotherapy but also in drug discovery and development. A quantitative relationship between in vitro CYP2D6 metabolism and in vivo PK has been constructed for evaluation of the impact of CYP2D6 polymorphic metabolism . The pharmacogenetics-based compartmental PK model (Fig. 1) described in current study was developed from the two genotyped animal models that differ in CYP2D6 drug-metabolizing activity [29, 33]. This model may be readily applied to human subjects (e.g., PMs and EMs) with different metabolic capacity caused by polymorphic drug-metabolizing enzyme. When linear drug clearance is assumed, the function of AUC ratio between PMs and EMs (RAUC = AUCPM/AUCEM = 1 + CLCYP2D6/CLother) derived from our model (Fig. 1) is essentially the same as that (RAUC = 1/(1 − fmCYP2D6)) in Gibbs model , in which fmCYP2D6, the fraction of drug clearance by CYP2D6 enzyme, equals to CLCYP2D6 divided by the sum of CLCYP2D6 and CLother. One may notice that fmCYP2D6 determined by in vitro assays with human liver microsomes or hepatocytes  is rather a fraction of hepatic clearance than total clearance. Therefore, use of in vitro fmCYP2D6 value, which overlooks renal or other drug clearance mechanisms, might provide an overestimation of the actual in vivo fmCYP2D6 and an exaggerated prediction of the impact of polymorphic metabolism on PK. Meanwhile, quantitative prediction of human PK from transgenic mouse data remains a challenge  unless reliable mouse-to-human scaling factors exist for drug clearance . Nevertheless, clinical studies is indispensable to gain a clear understanding of the impact of polymorphic metabolism on drug exposure, and our PK model (Fig. 1) may provide a helpful insight.
It is noteworthy that significant difference in PK may not be translated into considerable difference in PD, due to distinct PD relationship (e.g., the absence of linear dose-effect relationship)  and sensitivity of the assay for drug response. Harmaline is able to induce a variety of biochemical, physiological and pharmacological effects. Current study investigated the change of mouse body temperature and marble-burying behavior following the administration of harmaline. Although exposure to 5 mg/kg i.p. harmaline was significantly higher in wild-type mice than Tg-CYP2D6 mice (Table 3), the change of body temperature was not significantly different between the two genotyped mice (Fig. 5). In contrast, the anxiolytic activity of 5 mg/kg i.p. harmaline shown in mouse marble-burying test was significantly higher in wild-type mice (Fig. 6). Furthermore, wild-type mice exposed to15 mg/kg of harmaline showed more severe hypothermia than transgenic mice (Fig. 5), which is consistent with drug exposure data (Fig. 3 and Table 3). These findings suggest that PK/PD relationship need be understood for precise assessment of polymorphic metabolism on drug responses.
In summary, CYP2D6 drug-metabolizing enzyme plays an important role in harmaline clearance. Compared to human CYP2D6 EM hepatocytes, PM hepatocytes showed defective harmaline O-demethylase capacity, leading to lower rate of substrate depletion. Compared to Tg-CYP2D6 mice, wild-type mice had significantly higher and longer exposure to harmaline xenobiotics, which was translated into more severe hypothermia and largely-altered marble-burying behavior. Our findings suggest that distinct CYP2D6 status can cause considerable variations in harmaline metabolism, PK and PD, and CYP2D6 deficiency may lead to increased risk when high dose of harmaline is used. In addition, the pharmacogenetics-based PK model can be generalized to understand and assess the degree of PK difference caused by other polymorphic drug-metabolizing enzyme in different populations.
The project was supported by Award Number R01DA021172 from the National Institute On Drug Abuse (NIDA), National Institutes of Health (NIH). The authors also thank the Pharmaceutical Sciences Instrumentation Facility at University at Buffalo for use of LC-MS system that was obtained with Shared Instrumentation Grants S10RR014592 from the National Center for Research Resources, NIH.
Publisher's Disclaimer: This is a PDF file of an unedited manuscript that has been accepted for publication. As a service to our customers we are providing this early version of the manuscript. The manuscript will undergo copyediting, typesetting, and review of the resulting proof before it is published in its final citable form. Please note that during the production process errors may be discovered which could affect the content, and all legal disclaimers that apply to the journal pertain.