Next Article in Journal
Fabrication of Three-Dimensional Composite Scaffold for Simultaneous Alveolar Bone Regeneration in Dental Implant Installation
Next Article in Special Issue
Anti-Inflammatory and Anti-Oxidative Effects of luteolin-7-O-glucuronide in LPS-Stimulated Murine Macrophages through TAK1 Inhibition and Nrf2 Activation
Previous Article in Journal
Characterization and Evaluation of Resistance to Powdery Mildew of Wheat–Aegilops geniculata Roth 7Mg (7A) Alien Disomic Substitution Line W16998
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Application of an Inter-Species Extrapolation Method for the Prediction of Drug Interactions between Propolis and Duloxetine in Humans

1
College of Pharmacy, Chungnam National University, Daejeon 34134, Korea
2
College of Pharmacy, Seoul National University, Gwanak-ro 1, Gwanakgu, Seoul 08826, Korea
3
Department of Pharmaceutics, Ernest Mario School of Pharmacy, Rutgers, The State University of New Jersey, Piscataway, NJ 08854, USA
4
Department of Clinical Pharmacology and Therapeutics, Seoul National University College of Medicine and Hospital, Seoul 03080, Korea
5
Department of Pharmaceutical Sciences, Leslie Dan Faculty of Pharmacy, University of Toronto, Toronto, ON M5S 3M2, Canada
6
College of Pharmacy, Chung-Ang University, Seoul 06974, Korea
*
Authors to whom correspondence should be addressed.
Int. J. Mol. Sci. 2020, 21(5), 1862; https://doi.org/10.3390/ijms21051862
Submission received: 4 December 2019 / Revised: 28 February 2020 / Accepted: 5 March 2020 / Published: 9 March 2020

Abstract

:
Duloxetine (DLX) is a potent drug investigated for the treatment of depression and urinary incontinence. DLX is extensively metabolized in the liver by two P450 isozymes, CYP2D6 and CYP1A2. Propolis (PPL) is one of the popular functional foods known to have effects on activities of CYPs, including CYP1A2. Due to the high probability of using DLX and PPL simultaneously, the present study was designed to investigate the potent effect of PPL on pharmacokinetics (PKs) of DLX after co-administration in humans. A PK study was first conducted in 18 rats (n = 6/group), in which the plasma concentration of DLX and its major metabolite 4-hydroxy duloxetine (4-HD) with or without administration of PPL was recorded. Population PKs and potential effects of PPL were then analyzed using NONMEM software. Lastly, these results were extrapolated from rats to humans using the allometric scaling and the liver blood flow method. PPL (15,000 mg/day) exerts a statistically significant increase in DLX exposures at steady state, with a 20.2% and 24.6% increase in DLX C m a x , s s and the same 28.0% increase in DLX A U C s s when DLX (40 or 60 mg) was administered once or twice daily, respectively. In conclusion, safety issues are required to be attended to when individuals simultaneously use DLX and PPL at high doses, and the possibility of interactions between DLX and PPL might be noted.

Graphical Abstract

1. Introduction

Propolis (PPL) is a resin-like material made by honeybees from the buds of poplar and cone-bearing trees and is used in beehive construction and maintenance. PPL has been used in traditional medicine for centuries worldwide [1,2]. A significant body of evidence has shown that extracts of PPL have anti-inflammatory [3], anticancer [4,5], antioxidant [6], antifungal [7], antibacterial [8], and antivirus [9] activities. These biological activities depend mainly upon the chemical composition of the constituents, which vary considerably depending on the type of plants accessible to the bees as well as the sites of collection. Currently, PPL is widely used as popular functional foods and beverages that found in many health food stores for treating various diseases. One of the most common applications of PPL is to support the immune system for treatment of cold syndromes, including upper respiratory tract infections, the common cold, and flu-like infections. Since PPL is considered to be a health supplement, not a drug, PPL dosage is generally recommended by manufacturers on the product label without any guideline from heath organizations.
Previous literature has reported that PPL has protective effects against chemical-induced liver injury [10,11,12]. For example, PPL pretreatment protected against acetaminophen-induced hepatotoxicity in mice. One of the protective mechanisms was that PPL inhibited the formation of N-acetyl-p-benzoquinoneimine (an active metabolite of acetaminophen), which was produced by cytochrome P450 (CYP) enzymes [10]. Changes in the activity of CYP450 enzymes is one of the reasons for interactions between supplements and drugs that may cause adverse events. Therefore, effects of PPL and its major ingredients on CYP450 activities has been characterized recently [13,14,15,16]. For instance, Naramoto et al. has reported that ethanol extraction of Brazilian Green PPL inhibited enzymatic activities of all the tested CYP enzymes (CYP1A2, 2C9, 2C19, 2D6, and 3A4) in vitro in a concentration-dependent manner [13]. The inhibition effects of PPL on CYP1A2, 2E1, and 2C19 activities in vitro and in vivo have also been reported in our recently study [14]. Interestingly, none of published studies had evaluated the interactions of PPL with clinical drugs in humans, although PPL has been used widely as a health supplement.
Duloxetine (DLX), (+)-(S)-N-methyl-γ-(1-naphthyloxy)-2-thiophenrpropylamine, administered as duloxetine hydrochloride, is a potent and selective serotonin and norepinephrine reuptake inhibitor. This drug is used clinically in many countries for the treatment of disorders related to norepinephrine and serotonin, including major depressive disorder, diabetic peripheral neuropathic pain, and generalized anxiety disorder [17,18]. The drug is rapidly and extensively metabolized by CYP1A2, and to a lesser degree, by CYP2D6, to form multiple oxidative metabolites, which are then conjugated before being excreted in the urine [19]. Correspondingly, the glucuronide conjugate of 4-hydroxy duloxetine (4-HD) and the sulfate conjugate of 5-hydroxy-6-methoxy duloxetine are the two major metabolites of DLX in plasma [19]. A clinical study by Lobo et al. [20] reported that following a single oral DLX dose, the co-administration of fluvoxamine (a strong CYP1A2 inhibitor) resulted in clinically important increases in the DLX A U C 0 (460%) and C m a x   (141%). To explain these significant increases, the authors suggested that fluvoxamine increased the amount of DLX absorbed into systemic circulation by decreasing the extent of first-pass metabolism and, additionally, decreased systemic clearance of DLX by inhibiting CYP1A2-mediated metabolism.
Since only DLX has pharmacological effects, while all its circulating metabolites are inactive, any interaction resulting in a variation in DLX metabolism may alter its therapeutic effects, safety, and/or adverse events [17]. Nowadays, many people, especially elderly people, usually take drugs and supplementary simultaneously. As PPL is widely used for addressing various diseases, there is a high probability that DLX and PPL are administered at the same time. This study was conducted to investigate potential effects of PPL on the pharmacokinetic (PK) profile of DLX when they are co-administered in humans, using population PK model-based approach [21,22]. First, a PK study for drug interactions between DLX and PPL was conducted in rats. A population PK model was subsequently developed using nonlinear mixed-effects modelling (NONMEM) software to describe the PK profiles of DLX and its major metabolite 4-HD in rats. The results were then extrapolated from rats to humans by applying the simple allometric scaling method [23,24,25,26] and the liver blood flow (LBF) method [27,28].

2. Results

2.1. PK Data of DLX and 4-HD in Rats

Plasma concentration of DLX and 4-HD after a single oral dose of 40 mg/kg DLX or a single oral dose of 40 mg/kg DLX with 500 or 1500 mg/kg PPL in rats were determined using a high-performance liquid chromatography-tandem mass spectrometer (HPLC-MS/MS) method. The calibration curves established by plotting the peak area ratio vs. concentration were linear in the range of 5–1000 ng/mL for both of these two analytes (R2 > 0.99). During the method validation, all the test results met the required criteria according to the U.S Food and Drug Administration (FDA) guidelines for bioanalytical method validation [29]. Individual and mean of population of plasma concentration-time curves of DLX and 4-HD are illustrated in Supplementary Materials Figure S1.

2.2. Population PK Model for DLX and 4-HD in Rats

The final population PK model used to describe the PKs of DLX and its metabolite 4-HD with characterizing the metabolism of DLX converted to 4-HD was chosen after testing various models based on the objective function value (OFV) changes and goodness of fit (GOF) plots. A summary of building steps for covariate models resulting in statistical significance in the OFV during the model development is listed in Table 1. As demonstrated, PPL showed a significant inhibitory effect (p < 0.000) on the fraction of 4-HD converted from DLX ( F p m ) at the depot, as well as the conversion of DLX to 4-HD ( C L p m ) at the systemic circulation.
( F p m ) and systemic ( C L p m ) metabolisms.
PKs of DLX and 4-HD were well described by a one-compartment disposition model, with first-order elimination and first-order absorption kinetics. The metabolism of DLX converted to 4-HD was described by two metabolism pathways, including the first-pass effect and the systemic metabolism (Figure 1). The covariate effects of PPL on PK of DLX were modeled using a Michaelis Menten equation, as described in the following equations (Equations (1) and (2)) and illustrated in Figure 2.
C L p m = 1.26 × 1 1 × D o s e P P L D o s e P P L   +   276
F p m = 0.589 × 1 0.147 × D o s e P P L D o s e P P L   +   538
By using the final model, only clearance of DLX via all elimination routes except the conversion to 4-HD via the systemic metabolism from the body ( C L p ) need to be modeled with inter-individual variation effect. Estimates of fixed effects, intra- and inter-individual random effects, and their relative standard errors (% RSE) for the estimations are listed in Table 1. As shown in Table 2, the PK parameters of DLX and 4-HD in rats were generally well-estimated, with their % RSE in the range of 8.10%–85.7%, except IC50_ F p m (358%). The IIV estimate for C L p was 7.9% with an RSE of 70%. Accordingly, a decrease by 149% (from 149% to 0%) in IIV of C L p m and a decrease by 10.5% (from 18.4% to 7.9%) in IIV of C L p were observed in comparison to the base model. The proportional residual error was 19.9% (for DLX) and 24.0% (for 4-HD).
The basic GOF plots of the final model of the plasma concentration-time profiles of DLX and 4-HD are shown in Figure 3, suggesting that the individual prediction concentrations were reasonably well fitted by the final covariate model. The population and individual post-hoc predictions were distributed around the line of identity with no systematic bias. Moreover, conditional weighted residuals for the population predicted plasma concentrations were generally distributed around zero and were symmetric.

2.3. Population PK Model Evaluation

Evaluation and selection of the models were based on graphical analysis (visual predictive check, VPC plots) as well as statistical methods (bootstrap analysis). The predicted concentrations from 1000 simulated datasets using the final model for DLX and 4-HD after a single oral dose of 40 mg DLX and a variable dose of PPL in rats are shown in Figure 4. As illustrated in Figure 4, all of the observed concentrations were within the 5th–95th percentiles of the simulated concentrations. In addition, the majority of the observed median values were within the 95% confidence interval of the predicted median values obtained by the simulation. Results of the bootstrap analysis are listed in Table 2, indicating the parameter estimates of the final model were in the 2.5th–97.5th percentiles and were close to the medians estimated from bootstrap replicates. This result demonstrated that the parameter estimates of the final model were not statistically different from the parameters estimated by the bootstrap process.

2.4. Extrapolation of PKs of DLX and 4-HD in Humans

PK parameters of DLX and 4-HD in humans extrapolated from rats and from previous studies, including absorption rate constant K a , elimination rate constant K e , apparent oral clearance ( C L p m / F p and C L p / F p ), volume of distribution (   V p / F p ) of DLX, where F p is bioavailability of DLX after the first-pass effect, and apparent clearance C L m / F m and apparent volume of distribution   V m / F m of 4-HD, where F m is bioavailability of 4-HD, are listed in Table 3. Other PK parameters, including the effects of PPL dose on the PKs of DLX, as well as the IIV and RSE of all PK parameters in humans, were assumed to be the same as in the rats. The predicted total clearance ( C L D L X / F p ) and volume of distribution ( V p / F p ) of DLX was 55.9 L/h and 1022 L, respectively. Following conversion, the predicted apparent total clearance ( C L D L X / F D L X ) and volume of distribution ( V D L X / F D L X ) was 136 L/h and 2487 L, respectively, where F D L X is the total bioavailability of DLX. In other words, the fraction of drug that reaches the systemic circulation from the DLX dose is calculated following the Equation (3) below:
F D L X = 1 F p m × F p
The predicted PK profile of DLX and 4-HD following a single oral administration of DLX at the dose of 40 or 60 mg in humans is illustrated in Figure 5. The maximum concentration ( C m a x ) and area under the concentration-time curves (AUC) from zero to last time ( A U C 0 t ) of DLX are listed in Table 4.

2.5. Extrapolation of Effect of PPL on PKs of DLX in Humans

PK profiles of DLX and 4-HD after administration of multiple DLX 40 or 60 mg doses given once or twice daily co-administered with difference PPL doses (0, 5000, or 15,000 mg/day) in the first day and at the steady state condition are illustrated in Figure 6. The respective PK parameters of DLX are listed in Table 4. As shown in Table 4, DLX t1/2 increased by 1.1 h (from 12.7 to 13.8 h) and 2.6 h (from 12.7 to 15.3 h) when the drug was co-administered with PPL at doses of 5000 or 15,000 mg/day, respectively, comparing to when the drug was administered alone. Correspondingly, when DLX was administered once daily, the presence of PPL resulted in an increase by 11.6% and 28.0% in A U C during one dosing interval ( A U C S S ) and an increase by 8.35% and 20.2% in C m a x ( C m a x , S S ) at the steady state. Additionally, when DLX was administered twice daily, the presence of PPL resulted in an increase by 11.6% and 28.0% in A U C S S and an increase by 10.2% and 24.6% in C m a x , S S at the steady state, respectively.

3. Discussions

Drug interactions are one of the most common causes of medical adverse events. Therefore, evaluating potential drug interaction is an integral part of drug development. It is also important for making dosing recommendations for co-administered drugs. PK drug interactions may occur if a drug, a supplement, or food affects any one of the following four processes of another drug: absorption, distribution, metabolism, or elimination (ADME) [21,30,31].
PPL has been found to inhibit more than 50% activities of many enzymatic CYPs in a concentration-dependent manner, including CYP1A2. Changes in activity of CYP enzymes is thought to be a factor in interactions between supplements and drugs that may cause adverse events. While DLX has been reported to be a substrate of CYP450, the drug is rapidly and extensively metabolized by CYP1A2 and, to a lesser degree, metabolized by CYP2D6 to form multiple oxidative metabolites [19]. Consequently, there are probability of interactions between DLX and PPL through CYP1A2-mediated inhibition of PPL when they are co-administered in humans.
Approaches to predict human PK profiles from in vitro and/or PK profiles in preclinical species, utilizing both physiological and empirical approaches, are highly desirable [32,33,34,35,36,37,38,39], because this can drastically reduce the time and expense of drug development and dose regiment. Allometric scaling, an empirical method, has been intensively studied and widely applied to predict many important PK parameters, including CL, Vd, and t1/2, due to its simplicity [23,24,25,26]. However, allometric scaling has limitations, especially since it assumes that there are anatomical, physiological, and biochemical similarities among animals. The method may not be suitable for drugs that exhibit species-specific differences; for instance, drugs with high protein binding properties, significant biliary excretion, extensive active renal secretion, active metabolism, or show species-specific binding [40,41,42]. DLX is one of these cases. The drug is primarily cleared by hepatic metabolism via the activity of CYP1A2 and CYP2D6. Therefore, another method instead of a simple allometric scaling method, which can account for the difference in hepatic elimination functions between humans and rats, is required. As a result, an LBF method was selected, because it assumed that a drug clearance is proportional to the LBF of the subject; therefore, it can account for the difference in the liver functions of inter-species. By which methods, the mean apparent total clearance ( C L D L X / F D L X ), volume of distribution ( V D L X / F D L X ), and terminal elimination half-life t1/2 of DLX in humans was estimated to be 136 L/h, 2487 L, and 12.7 h, respectively. They are all in agreement with published literature [17,18,19,20,43,44,45,46,47]. For instance, according to the product information of Cymbalta (duloxetine hydrochloride) by the FDA, the volume of distribution of DLX in humans was 1640 L [17]. In the EMA product information, the apparent plasma clearance of DLX ranges from 33 to 261 L/h (mean 101 L/h) after oral administration of DLX [18].
Up until now, there have been no official criteria for extrapolated data from animals to humans. However, the criteria that the prediction of PK parameters within the range of one-half and twice (also known as two-fold dimensions) of the observed data are determined to be acceptable predictions has been suggested and applied in many published studies [28,48,49,50]. Therefore, we also applied the two-fold dimension criteria to assess the success/failure of our predicted data. Accordingly, the present study succeeded in predicting C L D L X / F D L X and V D L X / F D L X parameters of DLX in humans.
The absorption rate constant of DLX in humans was predicted based on the relationship between K a ,   K e , and time to reach C m a x   ( T m a x ), as described in Equation (4):
T m a x = ln K a ln K e K a K e
Given that DLX is acid labile, it therefore must be formulated with enteric-coated pellets to protect DLX from degradation in the acidic stomach environment by preventing DLX dissolution in the acidic stomach but allowing immediate release and rapid absorption at the small intestine. Consequently, after oral administration, there is a median 2-h lag time ( T l a g ) until absorption begins [17]. Therefore, the T m a x in Equation (4) was calculated as the time from when the absorption begins until the time when the C m a x reaches, Equation (5) as follows:
T m a x = T m a x , p o s t   d o s e T l a g
As a result, K a was predicted to be 0.697 1/h, which also is in agreement with the experimental estimate in humans by Sharma et al., where Ka was around 0.6 1/h [51].
By applying the above-extrapolated PK parameters from rats to humans, the mean C m a x and A U C 0 t of DLX after a single oral administration of 60 mg DLX were predicted to be 19.4 ng/mL and 441 h*ng/mL, respectively, which are reasonably similar to values reported in previous literature [20].
In conclusion, the PKs of DLX were successfully extrapolated from rats to humans by applying a method that accounted for both the pre-systemic and systemic metabolism of DLX to 4-HD. The extrapolation method developed here was further applied to investigate the potential effects of PPL, a CYP1A2 inhibitor, to the PKs of DLX when they are co-administered in humans.
Since PPL is considered to be only a health supplement, not a drug, the PPL dosage is generally recommended by manufacturers without any guidance from health organizations. PPL is designed in a capsule formulation, and one capsule contains up to 500 mg dry extract that is equivalent to up to 5000 mg fresh natural PPL. The recommended dose for this product is one to three capsules daily [52,53,54,55]. The recommended dosing regimen of DLX is 60 mg daily (given as 60 mg once daily or 30 mg twice daily) for major depressive disorder, diabetic peripheral neuropathic pain, and generalized anxiety disorder [17]. For stress urinary incontinence, the recommended dose is 40 mg twice daily. Based on the response of patients to treatment, the dose may be adjusted up to 120 mg daily (given as 60 mg twice daily). To mimic the therapeutic dosage of DLX and PPL, we extrapolated our data at different scenarios, including DLX 40 or 60 mg given once or twice daily and PPL 5000 mg/day (the medium dose) or 15,000 mg/day (the maximum dose). Due to a terminal elimination half-life of 12 h, which allows steady-state conditions to be achieved after three days, PKs of DLX and 4-HD were extracted after seven days post-dose to confirm the steady-state condition.
As a result, PPL (15,000 mg/day) exerts a statistically significant increase in DLX exposures at steady state, with a 20.2% and 24.6% increase in DLX C m a x , s s and the same 28.0% increase in DLX A U C s s when DLX was administered once or twice daily, respectively. However, a previous study has reported that despite higher DLX plasma concentrations during co-administration with fluvoxamine (138% increase in C m a x and 460% increase in A U C 0 ), adverse events were not observed more often than when DLX was used alone [20]. In addition, according to the FDA and European Medicines Agency (EMA), the safety of DLX up to 120 mg has been adequately established for the approved indications [17]. These results suggested that when DLX at the dose less than 120 mg was administered, the increased DLX exposure due to co-administration with PPL (even at the highest PPL dose) are still within the safety range, and it may not cause clinically important safety concerns. However, at the highest DLX dose (60 mg, twice daily; equivalent to 120 mg/day), safety issues need to be concerned, because there has been no safety test performed on doses greater than 120 mg/day. In addition, due to these changes in PKs of DLX, in all cases, we suggest that the probability of interaction between DLX and PPL should to be noticed by manufacturers and health organizations.
The present study had some limitations. The inhibitory CYP1A2 effect of PPL on PKs of DLX was conducted in rats, followed by extrapolation into healthy human subjects directly. The extrapolation was conducted with the assumption that there are anatomical, physiological, and biochemical similarities among rats and humans. However, animal populations may not exactly represent human populations in these conditions. In addition, healthy subjects also may not well represent the patient populations. Therefore, to make future decisions about using a treatment combination of DLX and PPL, for instance, whether the combination should be avoided or patient safety should be monitored during co-administration in humans, future clinical studies should be conducted in diverse human populations, especially patients with liver or renal failures and patients treated with DLX. Importantly, as we mentioned above, the biological activities in PPL depend mainly upon the chemical composition of the constituents, which vary considerably depending on the type of plants accessible to the bees as well as the sites of collection. As the results of the drug interaction between PPL and duloxetine could be different depending on the major ingredients and its respective contents of PPL, its interpretation should be paid attention to.

4. Material and Methods

4.1. Quantitative Determination of Major Ingredients in Propolis Extract

PPL capsules were purchased from Natural Immix Health Ltd. (Port Coquitlam, BC, Canada) and Country Life Ltd. (Hauppauge, NY, USA) and pooled at the ratio 1:1 (w/w). The PPL capsules contained only PPL extract from beehives of the honey bee without additives. PPL capsules were dissolved to 200 μg/mL in 0.1 M phosphate buffer (pH 7.4) containing 0.5% dimethyl sulfoxide (DMSO). Major ingredients in the PPL extract were then determined by a HPLC–MS/MS system as described in our previous report [14]. The quantities of major ingredients are shown in Table 5.

4.2. PK Study Design in Rats and Data Analysis

DLX and 4-HD were purchased from Waterstone Technology (Carmel, IN, USA) and Santa Cruz Biotechnology (Santa Cruz, CA, USA), respectively. Carbamazepine, used as an internal standard (IS) for analyzing DLX and 4-HD, was purchased from Sigma Chemical (St Louis, MO, USA).
The present study was approved by the Animal Ethics Committee of Chungnam National University (no. CNU-00721). It was conducted in 18 SD rats divided into 3 groups (G1, G2, and G3; n = 6 rats/group). Each was orally administered the same DLX dose (40 mg/kg) and a different PPL dose (G1, 0 mg/kg; G2, 500 mg/kg; and G3, 1500 mg/kg). PPL doses in the present study were the equivalent amounts of natural PPL. Blood samples (0.3 mL) were collected in Vacutainer® tubes (Becton Dickinson and Company, Franklin Lakes, NJ, USA) at pre-dose (0 h) and at various predetermined time-points (0.25, 0.5, 1, 2, 4, 8, 12, and 24 h) post-dose. Samples were centrifuged; the plasma was then transferred to polyethylene tubes and stored at −70 °C until analysis.
DLX and 4-HD were extracted from the plasma by a liquid-liquid extraction method, followed by analysis with a validated HPLC–MS/MS method as described in our previous report [56]. In detail, DLX and 4-HD were ionized in the positive electrospray ionization mode and detected through multiple reaction monitoring transitions at 298.00→44.00 and 314.13→154.10 m/z, respectively. Calibration curves were linear over a concentration range of 5–1000 ng/mL for both these two analytes (R2 > 0.99). Separation was performed using an Atlantis C18 (50 mm × 4.6 mm, 3 mm; Waters, Dublin, Ireland) column operated at 40 ± 5 °C. The mobile phase was a mixture of methanol with 5 nM ammonium acetate at a volume ratio of 60:40 and a flow rate of 0.3 mL/min.

4.3. Rat Population PK Model Development

PK analysis was conducted using a software package for nonlinear mixed-effects modeling (NONMEM version 7.3.0, ICON Development Solutions, Hanover, MD, USA) [57]. The PK model was estimated using the first-order conditional estimation (FOCE) method with interaction. First, a base model was developed without characterization of the effect of PPL on the PKs of DLX. The absorption and elimination were assumed to follow first-order kinetics. The lognormal distribution was used to model the IIV in the PK parameters, as described by Equation (6).
Pi = PTV × exp(ƞi)
where Pi is the parameter value for the ith individual, PTV is the typical population parameter, and ƞi (IIV) is the difference between Pi and PTV that is assumed to be normally distributed with a mean of zero and a variance of ωƞ2. A proportional error model was applied to describe the residual random error (ε) that contains contributions from IIV, assay error, and model misspecification, as expressed by Equation (7):
Cij = Cpred,ij × (1 + εpro,ij)
where Cij and Cpred,ij present the jth observed and predicted concentrations for the ith individual, respectively; εpro,ij is the proportional component that has normally distributed residual random effects with a mean of zero and a variance of σ2pro.
The proposed disposition compartment model for DLX and its metabolite, 4-HD, is schematically presented in Figure 1. This model suggests two pathways of 4-HD formation from DLX: (1) pre-systemic metabolism and (2) systemic metabolism. Correspondingly, the estimated population PK parameters for 4-HD were a fraction of the metabolite converted from DLX at the depot ( F p m ), which was exerted by the first-pass liver metabolism, and the rate of the metabolite converted from DLX circulating in the plasma ( K p m ), which exerted by systemic metabolism.
The estimated population PK parameters for DLX were apparent oral clearance ( C L p m / F p and C L p / F p ), volume of distribution ( V p / F p ), and the absorption rate constant K a . Where C L p m is the clearance of DLX via conversion into 4-HD; C L p is the clearance of DLX via all routes except for conversion into 4-HD; F p is bioavailability of DLX after the first-pass effect. In other words, F p is the fraction of the drug that reaches the systemic circulation in relation to the amount of DLX ( A m o u n t D L X ) that “survives” after passing the liver from the depot. The A m o u n t D L X is calculated following the Equation (8) below:
A m o u n t D L X = D o s e D L X × 1 F p m
Second, a covariate model to describe the possible effect of PPL on the PKs of DLX through its interaction with CYP1A2 was developed based on the base model. Covariate model building was accomplished by mixed stepwise forward addition (p < 0.05) and stepwise backward elimination (p < 0.01) based on change in the OFV using a likelihood ratio test within NONMEM, as well as reductions in IIV of PK parameters, GOF plots, and precision of estimates. Effects of PPL on the metabolism of DLX were modeled using a Michaelis Menten equation. For example, the effects of PPL on C L p m is described by the Equation (9) as follows:
C L p m = T V C L p m × 1 E m a x C L p m × D o s e P P L D o s e P P L   +   I C 50 C L p m
where T V C L p m is the typical value of C L p m in the population without co-administration with PPL, E m a x C L p m is the maximum inhibition effect of PPL on C L p m , and I C 50 C L p m is the co-administered dose (mg/kg) of PPL at which PPL has 50% of the maximum inhibition effect.

4.4. Rat Population PK Model Evaluation

Evaluation and selection of the models were based on graphical analysis as well as statistical methods. A nonparametric bootstrap analysis and a VPC were performed by simulating 1000 replicates each group using NONMEM. The final model assumed the same study design with same parameter estimates and their variances. Concentrations predicted from the VPC simulations were plotted against time for a 95% CI of the median, 5th, and 95th percentiles and overlaid with observed data. One thousand bootstrap replications were performed in PsN on the final model to obtain the median and 95% nonparametric CI (the 2.5th and 97.5th percentiles) for all parameters to assess parameter uncertainty.

4.5. Extrapolating Population PK of DLX in Humans

Assuming that an adult human has an average body weight of 70 kg, Vd of DLX and 4-HD in humans was extrapolated from the PK parameters in rats using a simple inter-species allometric scaling method [23,24,25,26], which assumes the similarities among animals can be generalized and expressed mathematically by an allometric relationship, as presented in Equation (10):
Y   = a × W b
where Y is the parameter of interest, W is the average body weight of a species, and a and b are the allometric coefficient and exponent of the equation, respectively. In the case of Vd, a fixed value of 1.0 was used for the exponent in the equation [23,24,25,26]. Accordingly, the Vd in humans was predicted from the rat PK parameter, as presented in Equation (11):
V h u m a n = V r a t × B W h u m a n B W r a t 1.00
Clearance (CL) in humans was extrapolated from rats following a method described by Ward and Smith [27], wherein CL is considered to be proportional to LBF. Accordingly, human CL was predicted following the Equation (12) below:
C L h u m a n = C L r a t × L B F h u m a n L B F r a t
where LBF values for rats and humans are 85 and 21 mL/min per kg of body weight, respectively [27].
The absorption rate constant was predicted based on an equation that describes the relationship between K a , elimination rate constant K e , and time to reach C m a x ( T max ) , as expressed by Equation (13):
T max = ln K a ln K e K a K e
where   T max was extracted from literature reports for PKs of DLX in humans and K e was calculated following the Equation (14) below:
K e h u m a n   = C L h u m a n V h u m a n = C L p m + C L p V p
To determine whether the extrapolation from rats to humans was successful or not, the extrapolated PK parameters were compared directly with the corresponding values estimated experimentally in humans. The extrapolation was assessed to be successful when the extrapolated human PK parameters lay within two-fold dimensions of the experimental values. Otherwise, it was assessed to be unsuccessful [28,48,49,50].

4.6. Simulation of PPL Effects on PKs of DLX in Humans

A potential effect of co-administered PPL on PKs of DLX in humans was performed by simulating 1000 replicates using the developed population PK model in rats with DLX and 4-HD. PK parameters in humans were obtained from extrapolation, as described in Section 2.4. The IIV between subjects in the human population was assumed to be the same as those in the rat population. In addition, the effects of PPL on the PKs of DLX were assumed to be constant across species.
Six scenarios for DLX and PPL co-administrations were simulated. For DLX, the drug was administered as follows: (1) 40 or 60 mg once daily and (2) 40 or 60 mg twice daily. For PPL, the supplement was co-administered with DLX at a total dose of (1) 0 mg, (2) 5000 mg, or (3) 15,000 mg of PPL per day. The potential effects of PPL on the PKs of DLX were assessed by comparing AUC and C m a x of DLX with co-administration (PPL = 5000 or 15,000 mg) and without co-administration (PPL = 0) of PPL.

5. Conclusions

After a single oral dose of DLX 40 mg/kg in rats, PKs of DLX and its metabolite 4-HD were modeled well by a one-compartment disposition model with first-order absorption and first-order elimination. While the metabolism from DLX to 4-HD was described by both pre-systemic (first-pass effect) and systemic metabolic pathways, a Michaelis-Menten equation was used to model the CYP1A2 inhibition kinetics of PPL on the metabolism of DLX converted to 4-HD. The PK profiles of DLX, as well as the interaction of PPL with the PKs of DLX in humans, were successfully extrapolated from the corresponding PK parameters in rats by applying an allometric scaling method (for predicting volume of distribution), an LBF method (for predicting clearance), and other PK parameters in humans extracted from experimental literature, with the final developed population model in rats. In conclusion, at the DLX doses less than 120 mg, the presence of PPL exerts a statistically significant increase in DLX exposures but does not seem to exert clinically significant changes in the PKs of DLX, even at the maximum dose of PPL (15,000 mg/day). At the maximum therapeutic dose of DLX (60 mg twice daily), safety issues are required to be attended to, and the possibility of interaction between DLX and PPL might be noted. Future studies further investigating these findings should be conducted in diverse human populations to provide a better representation.

Supplementary Materials

Supplementary Materials are available online at https://www.mdpi.com/1422-0067/21/5/1862/s1.

Author Contributions

Conceptualization, T.L.N., C.-H.L., H.-Y.Y., W.K. and J.-W.C; Data curation, T.L.N. and C.-H.L.; Formal analysis, T.L.N., C.-H.L., H.-M.B. and S.-J.R.; Investigation, T.L.N. and K.N.; Methodology, T.L.N., C.-H.L., N.H., H.-M.B., H.-Y.Y. and J.-W.C; Supervision, W.K. and J.-W.C; Validation, T.L.N., N.H., S.-J.R. and K.N.; Writing – original draft, T.L.N.; Writing – review & editing, T.L.N., H.-Y.Y., W.K. and J.-W.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the National Research Foundation of Korea (NRF- 2018R1C1B5085278).

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Ghisalberti, E.L. Propolis: A Review. Bee World 1979, 60, 59–84. [Google Scholar] [CrossRef]
  2. Castaldo, S.; Capasso, F. Propolis, an old remedy used in modern medicine. Fitoterapia 2002, 73, S1–S6. [Google Scholar] [CrossRef]
  3. Paulino, N.; Teixeira, C.; Martins, R.; Scremin, A.; Dirsch, V.; Vollmar, A.; Abreu, S.; de Castro, S.; Marcucci, M. Evaluation of the Analgesic and Anti-Inflammatory Effects of a Brazilian Green Propolis. Planta Med. 2006, 72, 899–906. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. El-khawaga, O.-A.Y.; Salem, T.A.; Elshal, M.F. Protective role of Egyptian propolis against tumor in mice. Clin. Chim. Acta 2003, 338, 11–16. [Google Scholar] [CrossRef]
  5. Ozkul, Y.; Silici, S.; Eroğlu, E. The anticarcinogenic effect of propolis in human lymphocytes culture. Phytomedicine 2005, 12, 742–747. [Google Scholar] [CrossRef]
  6. Ahn, M.-R.; Kunimasa, K.; Kumazawa, S.; Nakayama, T.; Kaji, K.; Uto, Y.; Hori, H.; Nagasawa, H.; Ohta, T. Correlation between antiangiogenic activity and antioxidant activity of various components from propolis. Mol. Nutr. Food Res. 2009, 53, 643–651. [Google Scholar] [CrossRef] [PubMed]
  7. Shehu, A.; Ismail, S.; Rohin, M.; Harun, A.; Aziz, A.; Haque, M. Antifungal Properties of Malaysian Tualang Honey and Stingless Bee Propolis against Candida Albicans and Cryptococcus Neoformans. J. Appl. Pharm. Sci. 2016, 044–050. [Google Scholar] [CrossRef] [Green Version]
  8. Kartal, M.; Yıldız, S.; Kaya, S.; Kurucu, S.; Topçu, G. Antimicrobial activity of propolis samples from two different regions of Anatolia. J. Ethnopharmacol. 2003, 86, 69–73. [Google Scholar] [CrossRef]
  9. Amoros, M.; Lurton, E.; Boustie, J.; Girre, L.; Sauvager, F.; Cormier, M. Comparison of the Anti-Herpes Simplex Virus Activities of Propolis and 3-Methyl-but-2-enyl Caffeate. J. Nat. Prod. 1994, 57, 644–647. [Google Scholar] [CrossRef]
  10. Seo, K.W.; Park, M.; Song, Y.J.; Kim, S.-J.; Yoon, K.R. The protective effects of Propolis on hepatic injury and its mechanism. Phytother. Res. 2003, 17, 250–253. [Google Scholar]
  11. González, R.; Corcho, I.; Remirez, D.; Rodriguez, S.; Ancheta, O.; Merino, N.; González, A.; Pascual, C. Hepatoprotective effects of propolis extract on carbon tetrachloride-induced liver injury in rats. Phyther. Res. 1995, 9, 114–117. [Google Scholar] [CrossRef]
  12. Lin, S.C.; Lin, Y.H.; Chen, C.F.; Chung, C.Y.; Hsu, S.H. The hepatoprotective and therapeutic effects of propolis ethanol extract on chronic alcohol-induced liver injuries. Am. J. Chin. Med. 1997, 25, 325–332. [Google Scholar] [CrossRef] [PubMed]
  13. Naramoto, K.; Kato, M.; Ichihara, K. Effects of an Ethanol Extract of Brazilian Green Propolis on Human Cytochrome P450 Enzyme Activities in Vitro. J. Agric. Food Chem. 2014, 62, 11296–11302. [Google Scholar] [CrossRef] [PubMed]
  14. Ryu, C.S.; Oh, S.J.; Oh, J.M.; Lee, J.-Y.; Lee, S.Y.; Chae, J.; Kwon, K.; Kim, S.K. Inhibition of Cytochrome P450 by Propolis in Human Liver Microsomes. Toxicol. Res. 2016, 32, 207–213. [Google Scholar] [CrossRef] [Green Version]
  15. Sasaki, T.; Sato, Y.; Kumagai, T.; Yoshinari, K.; Nagata, K. Effect of health foods on cytochrome P450-mediated drug metabolism. J. Pharm. Health Care Sci. 2017, 3, 14. [Google Scholar] [CrossRef]
  16. Šarić Mustapić, D.; Debeljak, Ž.; Maleš, Ž.; Bojić, M. The Inhibitory Effect of Flavonoid Aglycones on the Metabolic Activity of CYP3A4 Enzyme. Molecules 2018, 23, 2553. [Google Scholar] [CrossRef] [Green Version]
  17. European Medicines Agency. Cymbalta-Product Information; European Medicines Agency: Amsterdam, The Netherlands, 2004. [Google Scholar]
  18. U.S. Food and Drug Administration. Cymbalta-Product Information; U.S. Food and Drug Administration: Silver Spring, MD, USA, 2004. [Google Scholar]
  19. Lantz, R.J.; Gillespie, T.A.; Rash, T.J.; Kuo, F.; Skinner, M.; Kuan, H.-Y.; Knadler, M.P. Metabolism, excretion, and pharmacokinetics of duloxetine in healthy human subjects. Drug Metab. Dispos. 2003, 31, 1142–1150. [Google Scholar] [CrossRef] [Green Version]
  20. Lobo, E.D.; Bergstrom, R.F.; Reddy, S.; Quinlan, T.; Chappell, J.; Hong, Q.; Ring, B.; Knadler, M.P. In vitro and in vivo evaluations of cytochrome P450 1A2 interactions with duloxetine. Clin. Pharmacokinet. 2008, 47, 191–202. [Google Scholar] [CrossRef]
  21. U.S. Food and Drug Administration. Clinical Drug Interaction Studies—Study Design, Data Analysis, and Clinical Implications Guidance for Industry; U.S. Food and Drug Administration: Silver Spring, MD, USA, 2017. [Google Scholar]
  22. European Medicines Agency Guideline on the Investigation of Drug Interactions; European Medicines Agency: Amsterdam, The Netherlands, 2012.
  23. Boxenbaum, H. Interspecies scaling, allometry, physiological time, and the ground plan of pharmacokinetics. J. Pharmacokinet. Biopharm. 1982, 10, 201–227. [Google Scholar] [CrossRef]
  24. Tang, H.; Hussain, A.; Leal, M.; Mayersohn, M.; Fluhler, E. Interspecies Prediction of Human Drug Clearance Based on Scaling Data from One or Two Animal Species. Drug Metab. Dispos. 2007, 35, 1886–1893. [Google Scholar] [CrossRef] [Green Version]
  25. Caldwell, G.W.; Masucci, J.A.; Yan, Z.; Hageman, W. Allometric scaling of pharmacokinetic parameters in drug discovery: Can human CL, Vss and t1/2 be predicted fromin-vivo rat data? Eur. J. Drug Metab. Pharmacokinet. 2004, 29, 133–143. [Google Scholar] [CrossRef] [PubMed]
  26. Dedrick, R.L. Animal scale-up. J. Pharmacokinet. Biopharm. 1973, 1, 435–461. [Google Scholar] [CrossRef]
  27. Ward, K.W.; Smith, B.R. A comprehensive quantitative and qualitative evaluation of extrapolation of intravenous pharmacokinetic parameters from rat, dog, and monkey to humans. I. Clearance. Drug Metab. Dispos. 2004, 32, 603–611. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Nagilla, R.; Ward, K.W. A comprehensive analysis of the role of correction factors in the allometric predictivity of clearance from rat, dog, and monkey to humans. J. Pharm. Sci. 2004, 93, 2522–2534. [Google Scholar] [CrossRef] [PubMed]
  29. U.S. Food and Drug Administration Bioanalytical Method Validation: Guidance for Industry. Available online: https://www.fda.gov/files/drugs/published/Bioanalytical-Method-Validation-Guidance-for-Industry.pdf (accessed on 15 June 2018).
  30. Huang, S.-M.; Temple, R.; Throckmorton, D.C.; Lesko, L.J. Drug Interaction Studies: Study Design, Data Analysis, and Implications for Dosing and Labeling. Clin. Pharmacol. Ther. 2007, 81, 298–304. [Google Scholar] [CrossRef]
  31. Palleria, C.; Di Paolo, A.; Giofrè, C.; Caglioti, C.; Leuzzi, G.; Siniscalchi, A.; De Sarro, G.; Gallelli, L. Pharmacokinetic drug-drug interaction and their implication in clinical management. J. Res. Med. Sci. 2013, 18, 601–610. [Google Scholar]
  32. Kim, S.H.; Kim, W.B.; Lee, M.G. Interspecies pharmacokinetic scaling of a new carbapenem, DA-1131, in mice, rats, rabbits and dogs, and prediction of human pharmacokinetics. Biopharm. Drug Dispos. 1998, 19, 231–235. [Google Scholar] [CrossRef]
  33. Shim, H.J.; Kim, Y.C.; Lee, J.H.; Kwon, J.W.; Kim, W.B.; Kim, Y.G.; Kim, S.H.; Lee, M.G. Interspecies pharmacokinetic scaling of DA-8159, a new erectogenic, in mice, rats, rabbits and dogs, and prediction of human pharmacokinetics. Biopharm. Drug Dispos. 2005, 26, 269–277. [Google Scholar] [CrossRef]
  34. Lave, T.; Dupin, S.; Schmitt, M.; Kapps, M.; Meyer, J.; Morgenroth, B.; Chou, R.C.; Jaeck, D.; Coassolo, P. Interspecies scaling of tolcapone, a new inhibitor of catechol- O -methyltransferase (COMT). Use of in vitro data from hepatocytes to predict metabolic clearance in animals and humans. Xenobiotica 1996, 26, 839–851. [Google Scholar] [CrossRef]
  35. Mahmood, I.; Balian, J.D. Interspecies Scaling: Predicting Pharmacokinetic Parameters of Antiepileptic Drugs in Humans from Animals with Special Emphasis on Clearance. J. Pharm. Sci. 1996, 85, 411–414. [Google Scholar] [CrossRef]
  36. Mahmood, I. Application of allometric principles for the prediction of pharmacokinetics in human and veterinary drug development. Adv. Drug Deliv. Rev. 2007, 59, 1177–1192. [Google Scholar] [CrossRef] [PubMed]
  37. Yassen, A.; Olofsen, E.; Kan, J.; Dahan, A.; Danhof, M. Animal-to-Human Extrapolation of the Pharmacokinetic and Pharmacodynamic Properties of Buprenorphine. Clin. Pharmacokinet. 2007, 46, 433–447. [Google Scholar] [CrossRef] [PubMed]
  38. Richter, W.F.; Gallati, H.; Schiller, C.D. Animal pharmacokinetics of the tumor necrosis factor receptor-immunoglobulin fusion protein lenercept and their extrapolation to humans. Drug Metab. Dispos. 1999, 27, 21–25. [Google Scholar] [PubMed]
  39. Ward, K.W.; Coon, D.J.; Magiera, D.; Bhadresa, S.; Struharik, M.; Lawrence, M.S. Exploration of the African green monkey as a preclinical pharmacokinetic model: Oral pharmacokinetic parameters and drug-drug interactions. Xenobiotica 2009, 39, 266–272. [Google Scholar] [CrossRef]
  40. Mahmood, I. Integration of in Vitro Data and Brain Weight in Allometric Scaling to Predict Clearance in Humans: Some Suggestions. J. Pharm. Sci. 1998, 87, 527–529. [Google Scholar] [CrossRef]
  41. Lavé, T.; Coassolo, P.; Reigner, B. Prediction of hepatic metabolic clearance based on interspecies allometric scaling techniques and in vitro-in vivo correlations. Clin. Pharmacokinet. 1999, 36, 211–231. [Google Scholar]
  42. Huang, Q.; Riviere, J.E. The application of allometric scaling principles to predict pharmacokinetic parameters across species. Expert Opin. Drug Metab. Toxicol. 2014, 10, 1241–1253. [Google Scholar] [CrossRef]
  43. Tianmei, S.; Knadler, M.P.; Lim, M.T.; Yeo, K.P.; Teng, L.; Liang, S.; Pan, A.X.; Lobo, E.D. Pharmacokinetics and tolerability of duloxetine following oral administration to healthy Chinese subjects. Clin. Pharmacokinet. 2007, 46, 767–775. [Google Scholar] [CrossRef]
  44. Lobo, E.D.; Loghin, C.; Knadler, M.P.; Quinlan, T.; Zhang, L.; Chappell, J.; Lucas, R.; Bergstrom, R.F. Pharmacokinetics of duloxetine in breast milk and plasma of healthy postpartum women. Clin. Pharmacokinet. 2008, 47, 103–109. [Google Scholar] [CrossRef]
  45. Skinner, M. Duloxetine is both an inhibitor and a substrate of cytochrome P4502D6 in healthy volunteers. Clin. Pharmacol. Ther. 2003, 73, 170–177. [Google Scholar] [CrossRef]
  46. Li, H.; Li, T.; Li, Y.; Shen, Y. Pharmacokinetics and Safety of Duloxetine Enteric-coated Tablets in Chinese Healthy Volunteers: A Randomized, Open-label, Single- and Multiple-dose Study. Clin. Psychopharmacol. Neurosci. 2013, 11, 28–33. [Google Scholar] [CrossRef] [PubMed]
  47. Chan, C.; Yeo, K.P.; Pan, A.X.; Lim, M.; Knadler, M.P.; Small, D.S. Duloxetine pharmacokinetics are similar in Japanese and Caucasian subjects. Br. J. Clin. Pharmacol. 2007, 63, 310–314. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  48. Jolivette, L.J.; Ward, K.W. Extrapolation of human pharmacokinetic parameters from rat, dog, and monkey data: Molecular properties associated with extrapolative success or failure. J. Pharm. Sci. 2005, 94, 1467–1483. [Google Scholar] [CrossRef] [PubMed]
  49. Evans, C.A.; Jolivette, L.J.; Nagilla, R.; Ward, K.W. Extrapolation of preclinical pharmacokinetics and molecular feature analysis of “discovery-like” molecules to predict human pharmacokinetics. Drug Metab. Dispos. 2006, 34, 1255–1265. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  50. Feng, M.R.; Lou, X.; Brown, R.R.; Hutchaleelaha, A. Allometric pharmacokinetic scaling: Towards the prediction of human oral pharmacokinetics. Pharm. Res. 2000, 17, 410–418. [Google Scholar] [CrossRef]
  51. Sharma, A.; Goldberg, M.J.; Cerimele, B.J. Pharmacokinetics and safety of duloxetine, a dual-serotonin and norepinephrine reuptake inhibitor. J. Clin. Pharmacol. 2000, 40, 161–167. [Google Scholar] [CrossRef]
  52. Ferngrove Pharmaceuticals Australia Pty Ltd. Black Propolis 5000 mg. Available online: https://global.rakuten.com/en/store/the-natural/item/plo5000-100c/ (accessed on 14 April 2019).
  53. Ever Health EverHealth Propolis 5000, 120 Capsules. Available online: http://www.everhealth.co.nz/portfolio-item/propolis5000/ (accessed on 14 April 2019).
  54. My Natural Health Limited Health Vibrant Propolis 5000. Available online: https://www.mynaturalhealth.co.nz/products/view/health-vibrant-propolis-5000 (accessed on 25 February 2019).
  55. True Blue Premium Dark Bee Propolis 5000mg 250 Capsules Australian Made. Available online: https://www.amazon.com/True-Blue-Propolis-Capsules-Australian/dp/B0768SLYK8 (accessed on 15 February 2019).
  56. Chae, J.; Baek, H.; Kim, S.K.; Kang, H.; Kwon, K. Quantitative determination of duloxetine and its metabolite in rat plasma by HPLC-MS/MS. Biomed. Chromatogr. 2013, 27, 953–955. [Google Scholar] [CrossRef]
  57. Boeckmann, A.J.; Beal, S.L.; Sheiner, L.B. NONMEM User’s Guide, Part I; NONMEM Project Group, University of California at San Francisco: San Francisco, CA, USA, 1998. [Google Scholar]
Figure 1. Schematic representation of the base population pharmacokinetic (PK) model for DLX and its metabolite, 4-HD, and strategies to investigate the effect of PPL on the PKs of DLX after a single oral dose of 40 mg/kg DLX in rats. Notations are described in the main text. Pre-systemic ( F p m ) and systemic ( C L p m ) metabolisms.
Figure 1. Schematic representation of the base population pharmacokinetic (PK) model for DLX and its metabolite, 4-HD, and strategies to investigate the effect of PPL on the PKs of DLX after a single oral dose of 40 mg/kg DLX in rats. Notations are described in the main text. Pre-systemic ( F p m ) and systemic ( C L p m ) metabolisms.
Ijms 21 01862 g001
Figure 2. Predicted conversion fractions of DLX to 4-HD via the pre-systemic ( F p m ) and systemic ( C L p m ) metabolisms according to the different administered doses of PPL.
Figure 2. Predicted conversion fractions of DLX to 4-HD via the pre-systemic ( F p m ) and systemic ( C L p m ) metabolisms according to the different administered doses of PPL.
Ijms 21 01862 g002
Figure 3. Goodness of fit plots for the final model for predictions of the PK profiles of (A) DLX and (B) 4-HD in rats.
Figure 3. Goodness of fit plots for the final model for predictions of the PK profiles of (A) DLX and (B) 4-HD in rats.
Ijms 21 01862 g003
Figure 4. Visual predictive check plots for the final model for predictions of the PK profiles of (upper) DLX and (lower) 4-HD in rats. Open circles: observed concentrations, solid red line: median of the observed concentrations, dark grey area: 95% confidence interval (2.5th–97.5th percentiles) of the median predicted concentrations, and light grey area: 90% prediction intervals (5.0th to 95.0th percentiles) of the predicted concentrations.
Figure 4. Visual predictive check plots for the final model for predictions of the PK profiles of (upper) DLX and (lower) 4-HD in rats. Open circles: observed concentrations, solid red line: median of the observed concentrations, dark grey area: 95% confidence interval (2.5th–97.5th percentiles) of the median predicted concentrations, and light grey area: 90% prediction intervals (5.0th to 95.0th percentiles) of the predicted concentrations.
Ijms 21 01862 g004
Figure 5. Extrapolated concentration-time curve profiles of DLX and 4-HD after a single oral administration of 40 mg or 60 mg in humans. Data are presented as means (solid square, DLX and open circles, 4-DH) and standard deviations (error bars).
Figure 5. Extrapolated concentration-time curve profiles of DLX and 4-HD after a single oral administration of 40 mg or 60 mg in humans. Data are presented as means (solid square, DLX and open circles, 4-DH) and standard deviations (error bars).
Ijms 21 01862 g005
Figure 6. (A) Extrapolated concentration-time curve profiles of DLX and 4-HD at the first dosing interval (dashed line) and at the steady-state condition (solid line) following multiple DLX doses given once daily according to difference doses of PPL in humans: (black) PPL 0 mg, (red) PPL 5000 mg, and (blue) PPL 15,000 mg. (Upper) 40 mg and (lower) 60 mg; (B) Extrapolated concentration-time curve profiles of DLX and 4-HD at the first dosing interval (dashed line) and at the steady-state condition (solid line) following multiple DLX doses given twice daily according to difference doses of PPL in humans: (black) PPL 0 mg, (red) PPL 5000 mg, and (blue) PPL 15,000 mg. (Upper) 40 mg and (lower) 60 mg.
Figure 6. (A) Extrapolated concentration-time curve profiles of DLX and 4-HD at the first dosing interval (dashed line) and at the steady-state condition (solid line) following multiple DLX doses given once daily according to difference doses of PPL in humans: (black) PPL 0 mg, (red) PPL 5000 mg, and (blue) PPL 15,000 mg. (Upper) 40 mg and (lower) 60 mg; (B) Extrapolated concentration-time curve profiles of DLX and 4-HD at the first dosing interval (dashed line) and at the steady-state condition (solid line) following multiple DLX doses given twice daily according to difference doses of PPL in humans: (black) PPL 0 mg, (red) PPL 5000 mg, and (blue) PPL 15,000 mg. (Upper) 40 mg and (lower) 60 mg.
Ijms 21 01862 g006
Table 1. Summary of covariate model building steps.
Table 1. Summary of covariate model building steps.
Model No.Model Descriptiondf.Change in OFVCompared withSignificance (p-Value)
1Base model
2Model 1 with PPL as covariate for Fpm1−38.4Model 1<0.000
3Model 1 with PPL as covariate for CLpm2−28.7Model 1<0.000
4 aModel 1 with PPL as covariate for Fpm and CLp3−47.8Model 1<0.000
4 aModel 1 with PPL as covariate for Fpm and CLp2−10.6Model 2<0.005
OFV, objective function value; df., degrees of freedom; and a final model. Pre-systemic (Fpm) and systemic ( C L p m ) metabolisms.
Table 2. Parameter estimates from the final model and results of bootstrap validation for DLX and 4-HD after a single oral administration of DLX at dose of 40 mg without or with co-administration with PPL 500 or 1500 mg/kg in rats.
Table 2. Parameter estimates from the final model and results of bootstrap validation for DLX and 4-HD after a single oral administration of DLX at dose of 40 mg without or with co-administration with PPL 500 or 1500 mg/kg in rats.
ParametersUnitEstimatesRSE (%)Shrinkage (%)Bootstrap Replicates (n = 1000)
Median95% CI
K a 1/h1.358.10 1.351.241.45
Emax_ F p m 0.14785.7 0.1980.04350.844
IC50_ F p m mg/kg538358 80645.75478
F p m 0.58915.9 0.5740.4510.746
C L p L/h/kg1.9722.6 2.101.082.80
V p L/kg14.626.4 15.19.0119.6
Emax_ C L p m 1.00 1.001.001.00
IC50_ C L p m mg/kg27677.2 27060.1596
C L p m / F p L/h/kg1.2630.7 1.240.7701.77
C L m / F m L/h/kg12.313.6 11.99.4516.0
V m / F m L/kg84.221.6 81.364.3110
Inter-individual variability (IIV, %)
IIV C L p 7.9070.014.07.103.2110.2
Residual variability (%)
Prop_p 19.913.4 19.818.021.9
Prop_m 24.010.3 23.520.826.2
RSE, relative standard error; IIV, inter-individual variation; CI, confidence interval (the 2.5th and 97.5th percentiles); Prop_p and Prop_m are the proportional residual errors for DLX and 4-HD, respectively; F p m , conversion fractions of DLX to 4-HD via the pre-systemic metabolism; C L p m , clearance of DLX via the conversion to 4-HD via the systemic metabolism; C L p , clearance of DLX via all elimination routes except the conversion to 4-HD via the systemic metabolism; K a , V p , and F p : absorption rate constant, volume of distribution, and bioavailability of DLX after the pre-systemic metabolism, respectively; and C L m , V m , and F m : clearance, volume of distribution, and bioavailability of 4-HD, respectively.
Table 3. Predicted PK parameters of DLX and 4-HD in humans extrapolated from rats and experimental human parameters in literatures.
Table 3. Predicted PK parameters of DLX and 4-HD in humans extrapolated from rats and experimental human parameters in literatures.
ParametersUnitPredicted Value
  K a 1/h0.687
  C L p / F p L/h34.1
  C L p m / F p L/h21.8
  V p / F p L1022
  K e p 1/h0.0547
  C L m / F m L/h212
  C L m / F m L/h5894
An average of human body weight of 70 kg was assumed. K a was predicted based on the relationship between   K a ,   K e , and T m a x as: T m a x   =   ln K a ln K e K a K e , where   T m a x = 4 h and   K e   =   C L p m + C L p V p .
Table 4. Pharmacokinetic (PK) parameters of DLX in humans according to different co-administered doses of PPL.
Table 4. Pharmacokinetic (PK) parameters of DLX in humans according to different co-administered doses of PPL.
ScenarioParameterPPL 0 mgPPL 5000 mgDifference a (%)PPL 15,000 mgDifference b (%)
t 1 / 2 12.7 ± 0.60213.8 ± 0.699 15.3 ± 0.867
Single dose40 mg C m a x 12.9 ± 0.0928
A U C 0 t 294 ± 13.8
60 mg C m a x , S S   19.4 ± 0.139
A U C 0 t 441 ± 20.7
Multiple dose, once daily40 mg C m a x , S S   18.2 ± 0.55719.7 ± 0.6718.3521.9 ± 0.86820.2
A U C S S 294 ± 14.0328 ± 16.811.6376 ± 21.528.0
60 mg C m a x , S S   27.3 ± 0.83629.6 ± 1.018.3532.8 ± 1.3020.2
A U C S S 441 ± 21.0492 ± 25.211.6561 ± 32.328.0
Multiple dose, twice daily40 mg C m a x , S S   28.7 ± 1.1631.6 ± 1.3910.235.7 ± 1.7824.6
A U C S S 294 ± 14.0328 ± 16.811.6376 ± 21.528.0
60 mg C m a x , S S   43.0 ± 1.7447.4 ± 2.0810.253.6 ± 2.6724.6
A U C S S 441 ± 21.0492 ± 25.211.6564 ± 32.328.0
C m a x and C m a x , S S , maximum concentration of DLX after the first dose and at steady state condition, respectively; A U C 0 t ,area under the concentration-time curve from zero to last time after single DLX dose; A U C S S area under the concentration-time curve during one dosing interval at the steady state condition; the percentage difference in DLX PK parameters when drug was administered alone (PPL 0 mg) and (a) with PPL 5000 mg or (b) with PPL 15,000 mg, respectively; and unit: t 1 / 2 h; C m a x and C m a x , S S ng/mL; and A U C 0 t and A U C S S h*ng/mL.
Table 5. Major ingredients and their respective contents in PPL extract.
Table 5. Major ingredients and their respective contents in PPL extract.
IngredientContent (µg/mg) *
Chrysin23.57±3.02
Galangin7.45±0.51
Kaempferide4.30±0.57
Kaempferol4.19±0.27
Caffeic acid phenethyl ester0.35±0.01
Apigenin0.19±0.10
Artepillin C0.05±0.04
p-Coumaric acid0.03±0.01
Caffeic acid0.01±0.00
* Each value represents the mean ± standard deviation (n = 3 samples).

Share and Cite

MDPI and ACS Style

Ngo, T.L.; Lee, C.-H.; Han, N.; Back, H.-M.; Rhee, S.-J.; Noh, K.; Yun, H.-Y.; Kang, W.; Chae, J.-W. Application of an Inter-Species Extrapolation Method for the Prediction of Drug Interactions between Propolis and Duloxetine in Humans. Int. J. Mol. Sci. 2020, 21, 1862. https://doi.org/10.3390/ijms21051862

AMA Style

Ngo TL, Lee C-H, Han N, Back H-M, Rhee S-J, Noh K, Yun H-Y, Kang W, Chae J-W. Application of an Inter-Species Extrapolation Method for the Prediction of Drug Interactions between Propolis and Duloxetine in Humans. International Journal of Molecular Sciences. 2020; 21(5):1862. https://doi.org/10.3390/ijms21051862

Chicago/Turabian Style

Ngo, Thi Lien, Chung-Hee Lee, Nayoung Han, Hyun-Moon Back, Su-Jin Rhee, Keumhan Noh, Hwi-Yeol Yun, Wonku Kang, and Jung-Woo Chae. 2020. "Application of an Inter-Species Extrapolation Method for the Prediction of Drug Interactions between Propolis and Duloxetine in Humans" International Journal of Molecular Sciences 21, no. 5: 1862. https://doi.org/10.3390/ijms21051862

APA Style

Ngo, T. L., Lee, C. -H., Han, N., Back, H. -M., Rhee, S. -J., Noh, K., Yun, H. -Y., Kang, W., & Chae, J. -W. (2020). Application of an Inter-Species Extrapolation Method for the Prediction of Drug Interactions between Propolis and Duloxetine in Humans. International Journal of Molecular Sciences, 21(5), 1862. https://doi.org/10.3390/ijms21051862

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop