Next Article in Journal
Spontaneous Right Ventricular Pseudoaneurysms and Increased Arrhythmogenicity in a Mouse Model of Marfan Syndrome
Previous Article in Journal
Exosome Composition and Seminal Plasma Proteome: A Promising Source of Biomarkers of Male Infertility
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 21
Issue 19
10.3390/ijms21197023
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Prediction of Cyclosporin-Mediated Drug Interaction Using Physiologically Based Pharmacokinetic Model Characterizing Interplay of Drug Transporters and Enzymes

by
Yiting Yang
1,†,
Ping Li
1,†,
Zexin Zhang
1,
Zhongjian Wang
2,
Li Liu
1,* and
Xiaodong Liu
1,*
1
Center of Drug Metabolism and Pharmacokinetics, China Pharmaceutical University, Nanjing 210009, China
2
Jiangsu Key Laboratory of New Drug Research and Clinical Pharmacy, Xuzhou Medical University, Xuzhou 221004, China
*
Authors to whom correspondence should be addressed.
†
These authors contributed equally to this work.
Int. J. Mol. Sci. 2020, 21(19), 7023; https://0-doi-org.brum.beds.ac.uk/10.3390/ijms21197023
Submission received: 25 July 2020 / Revised: 13 September 2020 / Accepted: 18 September 2020 / Published: 24 September 2020
(This article belongs to the Section Biochemistry)
Browse Figures
Review Reports Versions Notes

Abstract

:
Uptake transporter organic anion transporting polypeptides (OATPs), efflux transporters (P-gp, BCRP and MRP2) and cytochrome P450 enzymes (CYP450s) are widely expressed in the liver, intestine or kidney. They coordinately work to control drug disposition, termed as “interplay of transporters and enzymes”. Cyclosporine A (CsA) is an inhibitor of OATPs, P-gp, MRP2, BCRP and CYP3As. Drug–drug interaction (DDI) of CsA with victim drugs occurs via disordering interplay of transporters and enzymes. We aimed to establish a whole-body physiologically-based pharmacokinetic (PBPK) model which predicts disposition of CsA and nine victim drugs including atorvastatin, cerivastatin, pravastatin, rosuvastatin, fluvastatin, simvastatin, lovastatin, repaglinide and bosentan, as well as drug–drug interactions (DDIs) of CsA with nine victim drugs to investigate the integrated effect of enzymes and transporters in liver, intestinal and kidney on drug disposition. Predictions were compared with observations. Most of the predictions were within 0.5–2.0 folds of observations. Atorvastatin was represented to investigate individual contributions of transporters and CYP3As to atorvastatin disposition and their integrated effect. The contributions to atorvastatin disposition were hepatic OATPs >> hepatic CYP3A > intestinal CYP3As ≈ efflux transporters (P-gp/BCRP/MRP2). The results got the conclusion that the developed PBPK model characterizing the interplay of enzymes and transporters was successfully applied to predict the pharmacokinetics of 10 OATP substrates and DDIs of CsA with 9 victim drugs.

1. Introduction

Drug influx transporters, efflux transporters and metabolic enzymes such as cytochrome P450 enzymes (CYP450s) and UDP-glucuronosyltransferases (UGTs) are widely expressed in the liver, intestine or kidney. They coordinately work to control the disposition of the drug in vivo, called as “interplay of transporters and enzymes” [1,2,3,4]. CYP3A4, the main CYP450s, accounts for 30–40% of hepatic CYP450s and mediates 60–70% of drug metabolism [5]. CYP3A4 is also abundant in the intestine, thus being implicated in the first-pass effect of orally administrated drugs. Uptake transporters organic anion transporting polypeptides (OATPs) such as OATP1B1, OATP1B3 and OATP2B1, are mainly expressed at the basolateral surface of hepatocytes and mediate hepatic uptake of their substrates from portal blood, becoming the limited step of hepatic clearance of some drugs such as statins [6]. Efflux transporters mainly include P-glycoprotein (P-gp), multidrug resistance-associated proteins (MRPs) and breast cancer resistance protein (BCRP). They are expressed in the brush border membrane of enterocytes, bile canalicular membrane of hepatocytes and apical membranes of proximal tubules, exporting their substrates out of cells to the intestinal lumen, bile and urine, respectively [7]. Some drugs are often substrates of these transporters and CYP450s. For example, atorvastatin is a substrate of P-gp, BCRP, MRP2, OATPs and CYP3A4, indicating that the disposition of atorvastatin should be attributed to the integrated effects of P-gp, BCRP, MRP2, OATPs and CYP3A4 [7,8,9,10] in the intestine and liver. Another example is cyclosporin A (CsA). CsA is a substrate of P-gp, OATPs and CYP3A4. Moreover, CsA is also a strong inhibitor of P-gp, BCRP, MRP2, OATPs and CYP3A4 [11,12,13]. Accumulating evidences have demonstrated that coadministration of CsA may lead potential drug–drug interactions (DDIs) with victim drug such as atorvastatin [14,15,16], cerivastatin [17], pravastatin [18,19], rosuvastatin [20], fluvastatin [21], simvastatin [22], lovastatin [18], repaglinide [23] and bosentan [24] via inhibiting these transporters or CYP3A4. Hepatic uptake of these drugs is mainly mediated by hepatic OATPs, which leads to high Kp values of drugs, then metabolized by CYP450s (main CYP3A4) and biliary or renal excretion via a P-gp, BCRP or MRP2-dependent mechanism. All these processes may be remarkably inhibited by CsA, indicating that CsA-mediated DDIs should be attributed to integrated effects of inhibiting transporters and CYP3A4 in the intestine, liver and kidney. Physiologically based pharmacokinetic (PBPK) model has been widely used to quantitatively predict the pharmacokinetics of drugs and the magnitude of DDIs [25]. Moreover, PBPK also illustrates individual contributions of transporters and enzymes to drug disposition as well as their integrated effects [26].
There have been researches about the DDI between statin and CsA with the PBPK model [27,28], but few have integrated simultaneously of the liver, intestinal and kidney on drug disposal. The highlight of our work is to investigate the integrated effect of enzymes and transporters in the liver, intestinal and kidney on drug disposition using a whole-body PBPK model characterizing the interplay of the enzyme-transporter in the liver, intestine and kidney. The developed PBPK first was used to predict the disposition of CsA and nine victim drugs (OATP substrates) including atorvastatin, cerivastatin, pravastatin, rosuvastatin, fluvastatin, simvastatin, lovastatin, repaglinide and bosentan, as well as DDIs of CsA with nine victim drugs. The predictions were compared with clinical reports. Atorvastatin, the substrate of uptake transporter (OATP), efflux transporter (BCRP, P-gp and MRP) and enzyme CYP3A, was represented to further investigate individual contributions of transporters and CYP3A4 to the disposition of drugs and to document sensitivity analysis.

2. Results

2.1. Collection of Drug Parameters Used for PBPK Prediction

Physiological parameters of humans such as organ volume, blood flow rate, metabolic and transport parameters of CsA and victim drugs were cited from publications PubMed. Ratios of drug concentration in tissue-to-plasma were cited from previous reports or estimated using the previously reported concentration in tissues and plasma or calculated using the previous method [29] based on tissue composition and physicochemical parameters of drugs. These parameters were list in Table 1, Table 2, Table 3 and Table 4. Based on the DDIs for criterions, nine victim drugs were selected including atorvastatin (substrate of P-gp, BCRP, MRP2, OATPs and CYP3A4) [8], cerivastatin (OATPs, P-gp, CYP3A4 and CYP2C8) [30,31], pravastatin (OATP and MRP2) [9,32], rosuvastatin (P-gp, BCRP and MRP2, OATPs) [8], fluvastatin (substrate of OATPs, P-gp, BCRP, MRP2, CYP2C9 and CYP3A4) [8,33], simvastatin (substrate of P-gp, OATPs and CYP3A4) [31,32], lovastatin (substrate of MRP2, OATPs and CYP3A4) [34], repaglinide (P-gp, BCRP, OATPs, CYP2C8 and CYP3A4) [35] and bosentan (MRP2, OATPs, CYP3A4 and CYP2C9) [36,37]. The CsA-mediated DDI data came from fifteen clinic reports as showed in follow-up sections.

2.2. Quantitative Prediction of CsA and Victim Drugs Using the Developed PBPK Model

Pharmacokinetic profiles (Figure 1) and main pharmacokinetic parameters (Table 5) of perpetrator CsA and nine victim drugs (atorvastatin, cerivastatin, pravastatin, rosuvastatin, fluvastatin, simvastatin, lovastatin, repaglinide and bosentan) were predicted using the developed PBPK model and parameters listed in Table 1, Table 2, Table 3 and Table 4. Model performance was assessed using absolute average fold error (AAFE) around each data point, as well as the ratio between observed and predicted peak concentration (Cmax) and area under the curve (AUC). The results showed that 62% predicted concentrations (24/39) and 96% predicted exposure parameters (75/78) fell within 0.5–2.0-fold of observations, inferring that the developed PBPK model was suitable for illustrating pharmacokinetic behaviors of CsA and the tested victim drugs.

2.3. Quantitative DDIs of CsA with Victim Drugs

CsA is a strong inhibitor of P-gp, BCRP, MRP2, OATP1B1 and CYP3A4, whose Ki values were reported to be 0.895 μM [91] 0.28 μM [55], 4.1 μM [27], 0.014 [27] μM and 2 μM [27], respectively. The DDIs of CsA co-administrated with nine victim drugs were predicted using the developed PBPK and compared with clinic observations. The DDIs were indexed as the ratio of AUC (AUCR) or Cmax (Cmax R) for victim drugs with CsA to without CsA. The dosage of CsA was set to be the highest dosage reported in clinic reports. The Pharmacokinetic profiles for coadministration of perpetrator CsA and victim drugs (atorvastatin, cerivastatin, pravastatin, rosuvastatin, fluvastatin, simvastatin, lovastatin, repaglinide and bosentan) were showed in Figure 2, Figure S1 and Table 6. The detailed interpretation was given one by one for the pharmacokinetic of each drug administrated with and without CsA.

2.3.1. DDI of CsA with Atorvastatin

Atorvastatin is a substrate of OATPs, CYP3A4, BCRP, P-gp and MRP2. Renal efflux transporters (BCRP, P-gp and MRP2) also mediated renal excretion of atorvastatin, although its contribution to total clearance was minor. Coadministration of CsA induced DDI with atorvastatin via inhibiting CYP3A4 and above transporters. Two clinical studies for DDIs of CsA and atorvastatin were included in the study.
Report 1 [14,16]. Eighteen renal transplant patients received CsA-based immunosuppressive therapy. The CsA dosage was 5.20 (2.03–12.66) mg/kg/d. The patients orally received daily 10 mg atorvastatin for four weeks. Pharmacokinetics of atorvastatin was performed following a four-week treatment. Eight age-matched healthy subjects were selected for controls and administrated daily 10 mg atorvastatin for one week [14]. The results showed that the simulated pharmacokinetic profile of atorvastatin following oral multi-dose of 40 mg atorvastatin was comparable to observed data (Figure 2A). Predicted and observed steady-state AUC0–24 and Cmax of atorvastatin in healthy subjects and in renal transplant patients are show in Table 6. Dosages of CsA were set to be 886 mg/70 kg/d. The predictions were within the fold error range compared with clinical observations. The predicted AUCR and CmaxR values were in line with observations.
Report 2 [15]. Thirteen healthy subjects orally received 40 mg atorvastatin for six days. On the evening of day 6, the first dose (2.5 mg/kg) of CsA was administrated. On day 7, CsA (2.5 mg/kg) was co-administrated with atorvastatin in the morning. Pharmacokinetics of atorvastatin after morning dose was implemented on day 6 and day 7. CsA dosage was set to be 175 mg/70 kg. The pharmacokinetic profiles and corresponding pharmacokinetic parameters of atorvastatin before and after the second dose of CsA were predicted using the developed PBPK model (Figure S1A). The results showed that the simulated pharmacokinetic profile of atorvastatin following oral multi-dose of 40 mg atorvastatin was comparable to observed data. Predicted steady-state AUC0–6 and Cmax following multi-dose of atorvastatin to healthy subjects were within 0.5–2.0 folds of clinic observations (data are shown in Table 6). Coadministration of CsA increased plasma exposure to atorvastatin. The predicted AUC0–6 and Cmax of atorvastatin following co-administration of CsA were lower than the clinic observations. The predicted AUCR and CmaxR were also lower than observed data.

2.3.2. DDIs of CsA with Cerivastatin

Cerivastatin is a substrate of P-gp, BCRP, OATP1B1 and CYP2C8. CsA-mediated DDI with cerivastatin is attributed to inhibiting OATP1B1-mediated hepatic uptake and P-gp/BCRP-mediated intestinal efflux. A clinic report [17] showed that 12 renal transplant recipients received CsA-based immunosuppressive therapy. The dosage regimen of CsA ranged from 75 mg to 225 mg twice daily. The patients were administrated 0.2 mg cerivastatin daily for seven days. Pharmacokinetics of cerivastatin were implemented on day 1 and day 7. Twelve healthy subjects were administrated a single dose of cerivastatin, serving as the control group. The predicted Cmax and AUC of cerivastatin following oral 0.2 mg of cerivastatin were good agreement with the measured data in healthy subjects. The dosage of CsA was assumed 225 mg/d. The predicted Cmax and AUC of cerivastatin following oral 0.2 mg of cerivastatin to renal transplant recipients on day 1 and day 7 were fell within 0.5–2.0 folds of observations, shown in Figure 2B and Table 6. Compared with healthy subjects, the predicted AUCR on day 7 was 5.2, falling within 0.5–2.0 folds of observation (4.8), but predicted CmaxR was 2.3, less than 0.5-fold observations (5.0).

2.3.3. DDIs of CsA with Pravastatin

Pravastatin is a substrate of MRP2 and OATPs. The drug mainly eliminates via MRP2-mediated biliary and MRP2-mediated renal secretion, indicating that CsA leads to DDI with pravastatin via affecting OATP-mediated hepatic uptake, MRP2-mediated intestinal efflux, MRP2-mediated biliary excretion and MRP2-mediated renal secretion. Two studies for DDIs of CsA with pravastatin were included in the study.
Report 1 [18]. Twenty-three stable kidney allograft recipients undergoing CsA-based immunosuppressive therapy received a single 20 mg oral daily dose of pravastatin for 28 days. Pharmacokinetics of pravastatin were implemented on day 1 and day 28 following pravastatin administration. CsA was administered as an oral capsule formulation twice daily, whose dosage was set to be 420 mg/d. The pharmacokinetic profiles of pravastatin following multidose of 20 mg to patients on day 1 and day 28 were predicted and compared with observations (Figure 2C). The Cmax and AUC0–24 of pravastatin co-administrated with CsA were predicted in Table 6. The predictions fell within 0.5–2.0 folds of clinic observations. The mean values in Table 5 were normalized by 20 mg, as controls, AUCR and CmaxR on day 28 were calculated to 3.1 and 3.6, respectively, within the fold error range compared with predictions (3.7 and 4.1 respectively).
Report 2 [19]. Eleven heart-transplant recipients underwent CsA-based immunosuppressive therapy. The patients received a daily dose of 40 mg/d pravastatin for the first eight days, then was reduced to 10 mg/d, administered until day 29. Pharmacokinetics of pravastatin were carried out on day 1, 8 and 29. Pharmacokinetics of pravastatin in eight healthy subjects following a single dose of 60 mg pravastatin was also investigated. The dosage of CsA in patients was set to 400 mg/d. The predicted Cmax and AUC0–24 values normalized by 10 mg pravastatin in healthy subjects were in line with clinic observations. The predicted 10 mg normalized Cmax values of pravastatin co-administrated CsA on days 1, 8 and 29 were also within the fold error range compared with observations. The predicted 10 mg-normalized AUC0–24 values in heart-transplant on day 1, 8 and 29 were lower than observations. The predicted CmaxR fell within 0.5–2.0 folds observations but estimated AUCR were lower than observations.

2.3.4. DDI of CsA with Rosuvastatin

Rosuvastatin is a substrate of P-gp, BCRP, MRP2 and OATPs. Rosuvastatin mainly eliminates via BCRP-dependent biliary and renal secretion. Thus, DDI of CsA with rosuvastatin may be due to inhibition of these transporter function. A report [20] demonstrated heart transplanted patients underwent a standard scheme of CsA-based therapy and 10 healthy subjects received a single dose of 10 mg rosuvastatin, followed by once-daily oral dose of 10 mg rosuvastatin for 10 days. Another group might receive a dose of 20 mg rosuvastatin, followed by a once-daily oral dose of 20 mg rosuvastatin for 10 days. Pharmacokinetics of rosuvastatin were documented at the first and last dose for patients and healthy subjects. CsA dosage was set to be 200 mg twice a day. The pharmacokinetics of rosuvastatin in healthy subjects and patients were predicted and compared with observation (Figure 2D and Table 6). The predicted Cmax and AUC of rosuvastatin following multi-dose of 10 mg rosuvastatin to healthy subjects were agreed with observations. The predicted Cmax and AUC of single-dose and multi-dose of 10 mg to patients were lower than that for measured geometric mean, the same for single-dose and multi-dose of 20 mg rosuvastatin to patients. The predicted AUCR and CmaxR following multi-doses to patients were 2.4 and 2.9, respectively, which were also lower than observations.

2.3.5. DDIs of CSA with Fluvastatin

Fluvastatin is metabolized by CYP2C9. Fluvastatin is also a substrate of P-gp, BCRP, MRP2 and OATPs. Liver BCRP also mediates biliary excretion of fluvastatin. Thus DDI of CsA with fluvastatin mainly results from inhibiting the above transporters. Two clinic reports were cited.
Report 1 [21]. Ten heart transplant patients under CsA-based therapy and 10 healthy subjects received 40 mg fluvastatin daily for 28 days. Pharmacokinetics was documented on day 1 and day 28. CsA dosage was assumed to be 200 mg/d. The plasma concentrations of fluvastatin in patients and healthy subjects were predicted and compared with observations (Figure 2E and Table 6). The predicted Cmax and AUC0–24 on day 1 and day 28 following oral doses of 40 mg fluvastatin to healthy subjects were within the 2-fold error range compared with the observations. The predicted Cmax and AUC0–24 on day 1 and day 28 following oral 40 mg of fluvastatin to heart transplant patients were also fell within 0.5-2.0 folds of clinical observations. On day 28, predicted CmaxR and AUCR were 3.4 and 3.2, respectively, which were in line with observations (6.0 and 3.1, respectively).
Report 2 [92]. Twenty hypercholesterolemic renal transplant recipients receiving CsA in combination with azathioprine and methylprednisolone administrated 20 mg fluvastatin daily for 4–6 weeks. CsA dosage was set to 200 mg/d. The predicted Cmax value was falling within 0.5–2.0 folds of observation, but the predicted AUC was twice higher than observation.

2.3.6. DDI of CsA with Simvastatin

Simvastatin is a substrate of OATPs and CYP3A4, which mainly eliminates via the CYP3A4-mediated metabolism mechanism. CsA leads to DDI with simvastatin via inhibiting OATP-mediated hepatic uptake, CYP3A4-mediated metabolism in the intestine and liver. Two clinic reports were included in the simulation.
Report 1 [93]. Pharmacokinetics of simvastatin following an oral dose of 20 mg was implemented in five transplant patients receiving CsA (1.1~3.8 mg/kg) and 5 patients without CsA. The dosage of CsA was set to be 266 mg/70 kg. The pharmacokinetic profiles of simvastatin following oral 20 mg simvastatin to patients treated with CsA or without CsA were predicted and compared with observations (Figure 2F and Table 6). The predicted AUC0–24 and Cmax values in patients treated with CsA were fell within ranges of observed AUC0–24 and Cmax. The predicted AUC0–24 and Cmax values in patients treated without CsA were lower than clinic observations. Predicted AUCR and CmaxR were 3.6 and 7.3, respectively, which were higher than observations (2.7 and 2.08). However, using observed mean values normalized by 20 mg in Table 5, AUCR values of Cmax and AUC were calculated to 6.5 and 6.9, respectively, within the 2-fold error range compared with predictions (7.3 and 3.6. respectively).
Report 2 [94]. Seven hypercholesterolemic heart transplant patients under CsA-based immunosuppressive therapy and seven hypercholesterolemic non-heart transplant patients (served as the control) were included. Both groups undergoing 10 mg simvastatin treatment for 6 weeks, attended the clinic in the morning after an overnight fast and received a double dose (20 mg). The dosage of CsA was set to be 200 mg twice daily. The predicted concentrations of simvastatin acid at 1, 2 and 3 h following oral 20 mg of simvastatin to heart transplant patients were higher than those in non-heart transplants. The observed mean plasma concentrations of simvastatin acid in heart transplant patients at 1, 2, and 3 h following oral 20 mg simvastatin were higher than those in non-heart transplants. The plasma concentrations of simvastatin acid in non-heart transplants were lower than 2 ng/mL. The predictions were higher than the observations.

2.3.7. DDI of CsA with Lovastatin

Lovastatin is a substrate of OATPs, MRP2 and CYP3A4, which mainly eliminates via CYP3A4-mediated metabolism mechanism. DDI of CsA with lovastatin occurs via inhibiting OATP-mediated hepatic uptake, MRP2-mediated intestinal efflux and CYP3A4-mediated metabolism. Three clinic reports were cited in the study.
Report 1 [18]. Twenty-one stable kidney allograft recipients undergoing CsA-based immunosuppressive therapy received a single 20 mg oral daily dose of lovastatin for 28 days. CsA dosage was between 2 and 6 mg/kg/day. Pharmacokinetic profiles of lovastatin were implemented on day 1 and day 28. CsA was administered as an oral capsule formulation twice daily. Pharmacokinetic profiles of lovastatin co-administrated 420 mg CsA were simulated and compared with observations (Figure 2G and Table 6). The predicted Cmax and AUC of lovastatin on day 1 and day 28 were lower than the observed medium values. Using the data list in Table 5 as controls, the calculated CmaxR and AUCR were within the 2-fold error range compared with predicted.
Report 2 [95]. Five types of patients were investigated as follows: (1) six heart transplant recipients receiving CsA (273 mg)-based triple immunosuppressive regimen; (2) five kidney transplant recipients receiving CsA (230 mg)-based triple immunosuppressive regimen; (3) five patients with psoriasis treated with CsA (290 mg) monotherapy; (4) five kidney transplant recipients receiving prednisolone and azathioprine; and (5) eight hypercholesterolemic patients serving as “control”. These patients orally received 10 mg lovastatin once daily for 10 days. Pharmacokinetics of lovastatin was implemented on day 10 (Figure S1B). The simulated AUC0–8 values in heart transplant recipients receiving CsA, kidney transplant recipients receiving CsA, patients with psoriasis treated with CsA, and hypercholesterolemic patients were 31.33, 27.87, 32.64, and 3.89 ng·h/mL. The measured mean AUC0–8 values in heart transplant recipients receiving CsA, kidney transplant recipients receiving CsA, patients with psoriasis treated with CsA and hypercholesterolemic patients were 175, 110, 110 and 26 ng·h/mL. Compared with hypercholesterolemic patients, AUCRs of lovastatin following co-administration of CsA were calculated to be 6.7, 4.2 and 4.2, respectively. The predicated AUCRs of lovastatin in heart and kidney transplant recipients and patients with psoriasis were 8.1, 7.2 and 8.4, respectively. All predictions were fell within 0.5–2.0 folds of observations.
Report 3 [96]. Twenty-four (14 cardiac and 10 renal) transplanted patients under CsA based-triple immunosuppressive regimen, received lovastatin treatment. Dosage of CsA was 300 (150–400) mg/d for heart transplanted patients and 255 (200–400) mg/d for renal transplanted patients. Five milligrams lovastatin was then given once daily in the morning for three weeks. The dose was thereafter increased by 5 mg every third week until 20 mg once daily, which was continued until week 18. Blood levels of lovastatin were measured prior to and 2 h after oral CsA and the first dose of lovastatin at each dose level. Blood levels of lovastatin in hypercholesterolaemic non-transplanted patients prior to and 2 h after 10 mg lovastatin were served as control. The dosages of CsA were set to be 300 mg. The predicted concentrations of lovastatin at 2 h following lovastatin 5 mg at weeks 0 and 3, 10 mg at week 6, 15 mg at week 9, and 20 mg at weeks 12 and 18 were 5.09, 5.93, 12.27, 18.91, 24.77 and 25.70 ng/mL, respectively. They were within 0.5–2.0-fold of observations. The predicted concentrations at 2 h following lovastatin 10 mg to hypercholesterolaemic non-transplanted patients was 1 ng/mL, lower than the observation concentration 4.6 ng/mL. Predicted AUCR using concentration at 2 h following 10 mg to transplanted patients at 6 weeks was 12, which was higher than observation (4.6).

2.3.8. DDI of CsA with Repaglinide

Repaglinide eliminates mainly via CYP2C8. To some extent, CYP3A4 also mediates repaglinide metabolism. Repaglinide is also a substrate of OATP1B1 and P-gp, indicating that DDI of CsA with repaglinide is involved in the inhibition of CYP3A4, P-gp, and OATPs.
A report [23] showed that 12 healthy subjects at PM 8 on day 1 and at 8 AM on day 2, received 100 mg CsA or placebo. At 9 a.m. on day 2, they ingested a 0.25 mg dose of repaglinide. Pharmacokinetics of repaglinide was documented. The plasma concentrations of repaglinide following 0.25 mg dose of repaglinide alone or coadministration of CsA were simulated and compared with observations (Figure 2H and Table 6). The predicted Cmax and AUC of repaglinide when co-administrated CsA were 2.4- and 2.5-fold higher than those following repaglinide alone. The measured Cmax and AUC of repaglinide when co-administrated CsA was 1.7 folds and 2.4 folds higher than those following repaglinide alone. All predictions fell within 0.5–2.0 folds of observations.

2.3.9. DDI of CsA with Bosentan

Bosentan is mainly eliminated via the metabolic mechanism, which is partly attributed to CYP3A4. Bosentan is also a substrate of MRP2 and OATPs and renal MRP2 may be involved in renal excretion of bosentan. These results inferred that DDI of CsA with bosentan is attributed to inhibition of CYP3A4-mediated metabolism, OATPs-mediated hepatic uptake and MRP2-mediated renal secretion.
A report [24] demonstrated that bosentan was given to 10 healthy subjects twice daily at a dose of 500 mg for 8 days. The subjects received 300 mg CsA twice daily starting from day 1 evening for 8 days. On days 1 and 8, pharmacokinetics of bosentan was documented following the morning dose. The plasma concentrations of bosentan were stimulated (Figure 2I and Table 6). The predicted Cmax and AUC of bosentan following coadministration of CsA were 2.2-fold and 3.76-fold higher than those before co-administration of CsA. The measured Cmax and AUC of bosentan following coadministration of CsA were higher than those before coadministration of CsA. Predicted CmaxR and AUCR were within the fold error range compared with the observed values.

2.4. Sensitivity Analysis of Model Parameters

Hepatic clearance of these tested agents is highly dependent on OATP-mediated uptake [26]. The empirical scaling factor (SF) is often introduced to gaps in translating in vitro transporter kinetic data. The Ki.OATP of CsA and clinic dosage of CsA largely varied. Five parameters, SF, CLint,uptake, Ki,OATP, CsA dosage and transition rate (Kt) on DDIs of CsA with atorvastatin were investigated. Variabilities of Ki,OATP SF, CLint,uptake, CsA dosage and Kt were set to be 0.5, 1 and 2.
We set atorvastatin as a typical drug to investigate the contribution of five parameters mentioned above. With the same structure of the PBPK model, we have reason to believe that other substances may with a similar effect by these factors. CsA dosage was set to be 200 mg/d. The results showed that the tested four factors remarkably affected DDIs of CsA with atorvastatin (Figure 3A–F), whose extents were CLint,uptake ≈ SF > CsA dosage > Ki,OATP > Kt. Individual contributions of hepatic CYP3A, integrated effects of the liver (i.e., hepatic CYP3A4-mediated metabolism/hepatic OATP-mediated hepatic uptake) and integrated effects of the intestine (i.e., intestinal CYP3A-mediated metabolism/efflux transporters) to DDI of CsA with atorvastatin were also investigated (Figure 3F). The results demonstrated that no considering inhibitions of CsA on CYP3A-mediated metabolism (non-CYP3A4) and integrated effects occurring in the liver (non-liver) or intestine (non-gut) led to remarkably lower concentrations of atorvastatin compared with actual concentrations of atorvastatin when co-administrated CsA. These results indicated that DDI of CsA with atorvastatin should be attributed to integrated effects of CYP3A4, OATP, P-gp, MRP2 and BCRP inhibitions. The pharmacokinetics parameters variation made by different factors listed in Figure 3 are shown in Table 7 (DDI Part).
Atorvastatin was used as mode drug to further investigate individual contributions of CYP3A, OATPs, P-gp, MRP and BCRP to atorvastatin disposition (Figure 4). The results showed that individual contributions to atorvastatin disposition following an oral dose were hepatic OATPs >> hepatic CYP3A > intestinal CYP3A ≈ intestinal P-gp ≈ intestinal MRP2 ≈ intestinal BCRP. The integrated contribution of intestinal CYP3A, intestinal P-gp, intestinal MRP2 and intestinal BCRP to atorvastatin disposition following oral dose contributions was about 20%. The pharmacokinetics parameters variation made by different factors listed in Figure 4 are shown in Table 7 (Alone Part).

3. Discussion

Drug disposition is highly dependent on the interplay of drug metabolism enzymes (such as CYP450 and UGTs) and transporters (such as OATPs, P-gp, BCRP and MRP2). These enzymes and transporters are widely expressed in the liver, intestine, or kidney. Some drugs, such as atorvastatin, are substrates of these enzymes and transporters, some of which (such as CsA) are also inhibitors of these enzymes and transporters. When they were co-administered, DDIs occur due to the inhibition of these enzymes and transporters. The extents of DDIs should be attributed to integrated effects on enzymes and transporters in these tissues. The aim of the study was to establish a PBPK model characterizing the interplay of enzymes and transporters in the intestine, liver and kidney, to predict DDIs of CsA with victim drugs.
Nine victim drugs were selected. They are all substrates of OATPs, some of which are also substrates of CYP3A4, BCRP, P-gp or MRP2. The developed PBPK model was first validated using nine victim drugs and perpetrator CsA using in vitro parameters of drug metabolism and transport. The lag time of simvastatin can be interpreted as the time for simvastatin (lactone) converse to simvastatin acid. The results showed that the developed PBPK model was successfully applied to predict pharmacokinetic profiles of nine victim drugs and CsA. Of predications, 96% pharmacokinetics parameters fell within 0.5–2.0 folds of observation. Following validation, the developed PBPK model was successfully scaled to DDIs of CsA with nine victim drugs. Except for rosuvastatin, all DDI predictions had the same trend with observations, demonstrating that the CsA-mediated DDIs may be predicted using the developed PBPK model involving enzymes and transporters in the liver, intestine and kidney.
Atorvastatin was served as a model drug for sensitivity analysis and individual contributions of enzymes and transporters to drug disposition. A sensitivity analysis that CLint,up, Ki,OATP, SF, the dosage of CsA and intestinal transit remarkably affected DDI of CsA with atorvastatin. The effect of SF, CLint,up was the strongest, followed by CsA dosage and Ki,OATP. The Ki values for CsA were reported to vary 6.3-folds when different OATP1B1 substrates were used in the in vitro transporter inhibition assays [97], showing a possible substrate-dependent inhibition of cyclosporine on OATP1B1 functions. In the clinic, the CsA dosage was varied, which may explain large variances in plasma concentrations of victim drugs when co-administered CsA. The empirical scaling factor (SF) was frequently applied during PBPK model development to capture the observed in vivo clearance [98]. The present study showed that each of the substrates had itself SF values and 300-fold varied (from CsA 0.1 to atorvastatin 32.6) among the tested agents. For a special substrate, the reported SF value was often different. For example, during simulating pharmacokinetics of atorvastatin. We also estimated SF using data from 10 reports. The results showed SF values had 5.4 times variations (ranging from 10.6 to 57). The medium value (32.6) was selected for simulating atorvastatin disposition, and good predictions were obtained. All these indicated that unreasonable SF value may lead to mispredictions. The ratio of drug concentration in tissue to plasma (Kp), which is used to describe drug distribution, is one of the most important factors for the PBPK model. We used a different method to get Kp value and selected the most reasonable one to continue our research showed in Table 2.
The simulation also demonstrated that intestinal BCRP, P-gp, MRP2 and CYP3A, hepatic OATPs and CYP3A differently contributed to atorvastatin disposition, hepatic OATPs > hepatic CYP3A > intestinal CYP3A ≈ intestinal BCRP ≈ intestinal P-gp ≈ intestinal MRP2. Systemic exposure of atorvastatin was more affected by hepatic OATP1B1 inhibition than by CYP3A4 inhibition, supporting that OATP1B1 was the rate-limiting step for the hepatic clearance of atorvastatin. All these results also demonstrated that DDI of CsA with atorvastatin should be mainly attributed to hepatic OATP inhibition, which may explain the findings that the exposure of atorvastatin increased 12-fold in the presence of the OATP1B1 inhibitor rifampin, but did not change in the presence of the strong CYP3A4 inhibitor itraconazole [99]. Stimulation also demonstrated that contribution of intestinal integrated effects (CYP3A-mediated metabolism and P-gp/BCRP/MRP2-mediated efflux) to atorvastatin disposition following oral dose was about 20%, leading to low intestinal bioavailability. High exposure to CsA following an oral dose of CsA occurred, inferring that extent of DDI following oral administration was larger than that following intravenous dose
However, the study also has some limitations. DDIs of CsA with statins were mainly documented in organ transplant patients under CsA-based therapy. The predictions of DDIs using the PBPK model were dependent on the assumption that diseases little affected pharmacokinetics of the indicated drugs. In fact, these diseases and the co-medicine themselves may affect pharmacokinetics. For example, rosuvastatin is mainly eliminated via the kidney. Heart failure was reported to also contribute to the decline of kidney function [100], which seemed to explain under- prediction of plasma exposure to rosuvastatin in heart transplanted patients. Moreover, OATP1B1 polymorphisms also affect the pharmacokinetics of statins [101] and repaglinide [101]. All these may lead to mispredictions, which needed further investigation.
In conclusion, the PBPK model characterizing the interplay of enzymes and transporters was successfully applied to predict the pharmacokinetics of OATP substrates and DDIs of CsA with 9 victim drugs. Most of the predictions fell within 0.5–2.0-fold of observations.

4. Method

4.1. Collection of Data

Clinical DDI studies involving CsA and OATP substrates were collected from publications on Pubmed. The CsA-mediated DDIs were based on the following criterions: (1) data might come from healthy subjects or patients following a single-dose or multi-dose of victim drugs when co-administrated with CsA. What is more, DDI data and the data for drugs used alone might come from different reports. (2) CsA was orally administrated to subjects in an immediate-release formulation. (3) Pharmacokinetic parameters such as peak concentration (Cmax) or area under the curve (AUC) or pharmacokinetic profiles were shown.

4.2. Development of a Whole-Body PBPK Model Charactering Interplay of Enzymes and Transporters in Intestine, Liver and Kidney

A whole-body PBPK model characterizing the interplay of enzymes and transporters in the intestine, liver and kidney (Figure 5) was developed to predict pharmacokinetic behaviors of victim drugs and CsA as well as DDIs of CsA with victim drugs. The PBPK model consisted of the lung, heart, brain, muscle, kidney skin, spleen, adipose tissues, other body tissues, gastrointestinal systems and liver [26,50,102].
In gastrointestinal systems
Gastrointestinal systems were divided into the stomach, gut lumen and gut wall. Assuming that absorption or metabolism of drug did not occur in stomach, drug amount (A0) in stomach was controlled by transition rate (Kti), i.e.,
dA 0 dt =   −   K ti × A 0
where the ith stand for different parts of gastrointestinal systems (i = stomach/duodenum/jejunum/ileum/cecum/colon).
The gut lumen was divided into five segments including duodenum, jejunum, ileum, caecum and colon according to their physiological and anatomical characteristics. Drug absorption was assumed to only occur at duodenum, jejunum and ileum.
The disposition of the drug in the ith gut lumen (i = duodenum/jejunum/ileum) was illustrated by
dA i dt = K ti   −   1 × A i   −   1   −   K ti × A i   −   k a , i × A i + k b , i × A gwi × T sf , i
where Kti, ka,i and Agwi were gut transit rate constant, absorption rate constant in the ith gut segment and the drug amount in the ith enterocyte compartment. kb,i is the efflux rate constant of the drug from enterocytes into the lumen. The expressions of intestinal P-gp, BCRP and MRP2 are regional [103], thus, transporter-mediated Peff, B-A in the ith gut segment was corrected by a relative transporter scaling factor (Tsf,i). The expressions of P-gp, BCRP and MRP2 in ileum were assumed to be 1. According to a previous report [103], the Tsf,i values of P-gp, BCRP and MRP2 in duodenum, jejunum and ileum were calculated as follows: 0.16:1.2:1.0 for P-gp, 0.68:1.45:1.0 for BCRP and 1.27:1.77:1.0 for MRP2. ka,i and kb,i may be calculated from human effective permeability in apical to basolateral direction (Peff,man,A-B) and in basolateral to apical direction (Peff,man,B-A) across the intestinal wall, i.e.,
k a , i = 2 × P eff , A - B / r i
k b , i = 2 × P eff ,   B - A / r i
where ri is the intestinal radius and the ri values for duodenum, jejunum and ileum were set to be 2.0, 1.63 and 1.45 cm [102], respectively. Peff, B-A was determined by intestinal P-gp, BCRP and MRP2, which should be attributed to the summed effect of these efflux transporters (P-gp, BCRP and MRP2). In the presence of CsA, Peff, B-A was rewritten as
P eff , B - A = P eff , B - A , P - gp 1 + ( A gwi I / V gwi ) / K i , P - gp + P eff , B - A , BCRP 1 + ( A gwi I / V gwi ) / K i , BCRP + P eff , B - A , MRP 2 1 + ( A gwi I / V gwi ) / K i , MRP 2
where Vgwi and A gwi I are the volume of the ith intestinal segment and amount of inhibitor in the ith intestinal segment. Peff, B-A, P-gp, Peff, B-A, BCRP, and Peff, B-A, MRP2 are efflux parameters mediated by P-gp, BCRP and MRP2, respectively. Ki, P-gp, Ki, BCRP and Ki, MRP2 are inhibition constants of corresponding transporters by CsA.
The drug amount (Agwi) in the ith gut wall is expressed as:
dA gwi dt =   Q gwi × A art V art   +   k a , i × A i   −   k b , i × A gwi × T sf , i   −   ( Q gwi + f ub   ×   CL int , gwi 1   +   ( A gwi I / V gwi ) / K i ) × A gwi / V gwi / K g : b
where Qgwi is the blood flow of the ith gut wall. Aart, Vart and Kg:b are the amounts of drug in artery blood, volume of artery blood and ratio of drug concentrations in intestine to blood, respectively. CLint, gwi is intrinsic clearance in the ith gut wall. Superscript “I” indicated perpetrator, CsA. The expression of intestinal CYP3A is also regional, expressions of CYP3A4 in duodenum, jejunum and ileum were 9.7 nmol, 38.4 nmol and 22.4 nmol [53], respectively. CLint,gut values of atorvastatin, CsA, simvastatin and lovastatin in intestinal microsomes were reported to be 0.03, 0.0277, 1.86 and 2.143 mL/mg microsome protein, respectively. From the report [104], there were 3155.2 mg microsomal protein/intestine for humans. According to the volume of the intestinal segment [105], median regional wet weights of duodenum, jejunum and ileum were 79 (0.098), 411 (0.39) and 319 (0.51) g, respectively. Assume that intestinal metabolic activity mainly comes from the jejunum. CLint, jejunum was calculated as 0.51 × 3155.2 × CLint, gut. While, with the reported expression of regional CYP3A4 [53], convert jejunal metabolism to CYP3A metabolism and we can get CLint, duodenum and CLint, ileum, which are listed in Table 4.
In artery blood (art):
dA art dt =   Q total × ( ( A lu / V lu ) / K lu : b − A art / V art )
where Qtotal, Vlu, Alu and Klu:b were cardiac output, lung volume, amount of drug in lung and ratio of drug concentration in lung to blood, respectively.
In venous blood (ven):
dA ven dt = ∑ Q i × ( A i / V i ) / K tissue : b −   Q total   ×   A ven / V ven
where Ktissue:b is ratio of drug concentration in different tissues to blood.
In liver:
The liver compartment was divided into hepatic blood and hepatocytes. Hepatic uptake of drugs was mainly controlled by OATPs and biliary excretion of drug was mediated by P-gp, BCRP or MRP2.
In hepatic blood (Ah, b):
dA h , b dt = ∑ Q gwi × A gwi V gwi × K g : b + Q sp × A sp V sp × K s : b + Q ha × A h , b V h   −   Q h × A h , b V h −   ( SF act × CL int , up , OATP 1 + f ub I × A h , b I V h , b × K i , OATP + CL int , pd ) ×   f ub × A h , b V h , b + f ub × CL int , pd × A h , c V h , c × K h : b
where Vh, b, Vh, c, Kh:p and Qh (Qh = Qha + Qsp + Qgwi) are the volume of hepatic blood and hepatocytes, ratio of drug concentrations in liver to blood, and blood flow in liver venous. CLint, up, OATP and CLint, pd are OATP-mediated intrinsic clearance of uptake and passive clearance, respectively. SP indicates the spleen. SFact represents the empirical scaling factor which was applied to scale-up the in vitro clearance to the in vivo hepatic clearance to make up for the underpredicted overall hepatic intrinsic clearance [27]. SFact, the value of atorvastatin, cerivastatin, pravastatin, rosuvastatin, repaglinide and bosentan, was cited from reference, while SFact, the value of cyclosporine, fluvastatin, simvastatin and lovastatin, was estimated from the software based on their pharmacokinetic data, which were shown in Table S1. At least three sets of clinical reports were used for SFact estimation, and their means were used for PBPK simulation. The estimated and cited SFact values are listed in Table 3.
In hepatocytes (Ah, c):
dA h , c dt = ( SF act × CL int , up , OATP 1 + f ub I   × A h , b I / V h , b / K i , OATP + CL int , pd ) ×   f ub × A h V h − ( CL int , pd + CL int , bile 1 + f ub I × A h , c I / V h , c / K i , P - gp / BCRP / MRP 2 + CL int , met , CYP 3 A 1 + f ub I × A h , c I / V h , c / K i , CYP 3 A + CL int , met , non - CYP 3 A )   ×   f ub × A h , c V h , c × K h : b
where CLint, met, CYP3A and CLint, met, non-CYP3A are intrinsic clearances of CYP3A-mediated metabolism and non-CYP3A-mediated metabolism, respectively. CLint, bile is the biliary excretion clearance by P-gp, BCRP or MRP2. Ki, OATP, Ki, CYP3A and Ki, P-gp/BCRP/MRP2 are inhibition constants of OATP-mediated uptake, CYP3A-mediated metabolism and P-gp, BCRP-or MRP2-mediated biliary excretion by CsA, respectively.
In kidney:
Renal excretion of drug mainly occurs via glomerular filtration and renal secretion, the renal clearance (CLren) is illustrated by
CL ren   =   f ub × GFR   +   CL sec
where GFR (set to be 120 mL/min/70 kg) and CLsec are glomerular filtration rate and renal secretory clearance by P-gp, BCRP or MRP2, respectively. CLsec was calculated using the equation:
CL sec   =   CL ren   −   f ub × GFR
Assuming that DDI occurred at secretion process, amount of drug in kidney is
dA k dt   =   Q k × A art V art −   A k   ×   f ub V k × K k : b × ( CL int , GFR + CL int , sec 1 + f ub I × A k I V k × K i , P - gp , BCRP , MRP 2 )
The CLint,GFR and CLint, sec were calculated using equation [106]
CL GFR ( or   CL sec ) =   f ub × CL int ,   GFR ( or   CL int ,   sec )   ×   Q k f ub × CL int ,   GFR ( or   CL int ,   sec )   +   Q k
where Qk is blood flow in kidney and CLGFR = fub × GFR. The CLren of cyclosporine A and bosentan were 3.44 and 2.4 mL/min. CLint, GFR of pravastatin and rosuvastatin were 132.9 mL/min. CLint,sec of pravastatin (MRP2) and rosuvastatin (BCRP) were 1538.96 mL/min and 420.3 mL/min, respectively.
The amount (Ai) of drug in other tissue is illustrated by the following formulation.
dA i dt = Q i × ( A art V art − A i V i × K tissue : b )
where Qi is blood flow in other tissues.
Phoenix WinNonlin 8.1 (Pharsight, St. Louis, MO, USA) was used for coding and solving of the PBPK model as well as estimation of kinetic parameters (Cmax and AUC).

4.3. Model Validation

The criteria for selecting victim drugs: 1. the victim drugs selected should be used clinically at present; 2. they all must be substances of OATP, some of which are substances of CYP3A; 3. they all must have clinical data of interactions with CsA, some of which may have strong DDI effect with CsA. Only the ones that meet all conditions were selected. The selected drugs may be combined with cyclosporine in certain patients clinically.
There were two steps in this process. Pharmacokinetic profiles and their parameters (AUC and Cmax) of nine victim drugs and CsA following oral administration to humans were first predicted using the developed PBPK models and compared with clinic observations. Following validation, the developed PBPK model was separately used to predict DDIs of CsA with victim drugs and confirmed with clinic data. The dosage of CsA was set to be the highest dosage reported in clinic reports. The extent of DDI was assessed as the ratio of AUC (AUCR) or Cmax (CmaxR) with CsA to without CsA and absolute average fold error (AAFE) around each data point. The predictions were considered successful if the ratio of predication to observation and AAFE fell between 0.5 and 2.0.
AAFE =   10 1 n ∑ i = 1 n | log 10 C Pred C Obs | .

4.4. Sensitivity Analysis of Model Parameters

Hepatic uptakes of the tested victims are mainly mediated by OATPs. Several pieces of evidence showed that the inhibitory effects of cyclosporin A were long-lasting [107], whose Ki showed about a 100-fold difference between pre-incubation non-pre-incubation. Moreover, the clinical dosage of CsA often varied. During in vitro–in vivo extrapolation of transporter-mediated clearance, the empirical scaling factor (SF) is frequently introduced to capture the observed in vivo clearance [98]. Sensitivity analysis was operated on CLint,up, Ki,OATP, SF, CsA dosage and rate of intestinal transit.
Atorvastatin is a substrate of OATPs, P-gp, BCRP and CYP3A. CsA is a strong inhibitor of OATPs, P-gp, BCRP, MRP2 and CYP3A. The individual contributions of BCRP, P-gp, MRP2, CYP3A, OATPs to DDI of atorvastatin with CsA were investigated.

Supplementary Materials

Supplementary Materials can be found at https://0-www-mdpi-com.brum.beds.ac.uk/1422-0067/21/19/7023/s1.

Author Contributions

Y.Y. and P.L. with the same contributed to Investigation, Methodology and Writing original draft. Z.Z. contributed to the Data curation and Z.W. contributed to Conceptualization. L.L. contributed to Funding acquisition and Software. X.L. contributed to Funding acquisition, Formal analysis and Writing original draft. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China, No. 81872930 and 81673505; the “Double First-Class” university project, No. CPU2018GY22.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

CYP450sCytochrome P450 enzymes
UGTsUDP-glucuronosyltransferases
OATPsOrganic Anion transporting polypeptides
P-gpP-glycoprotein
MRPMultidrug resistance-associated protein
BCRPBreast cancer resistance protein
CsACyclosporin A
DDIsDrug–Drug Interactions
PBPKPhysiologically based pharmacokinetic
fubunbound fraction of drug in blood
Rbblood/plasma ratio
KtTransition rate
SFEmpirical scaling factor
AUCArea under the curve
CmaxPeak concentration
AUCRRatio of AUC for victim drugs with CsA to without CsA
CmaxRRatio of Cmax for victim drugs with CsA to without CsA

References

  1. Lu, C.; Di, L. In vitro and in vivo methods to assess pharmacokinetic drug- drug interactions in drug discovery and development. Biopharm. Drug Dispos. 2020, 41, 3–31. [Google Scholar] [CrossRef] [Green Version]
  2. Bteich, M.; Poulin, P.; Haddad, S. The potential protein-mediated hepatic uptake: Discussion on the molecular interactions between albumin and the hepatocyte cell surface and their implications for the in vitro-to-in vivo extrapolations of hepatic clearance of drugs. Expert Opin. Drug Metab. Toxicol. 2019, 15, 633–658. [Google Scholar] [CrossRef]
  3. Liu, X. SLC Family Transporters. Adv. Exp. Med. Biol. 2019, 1141, 101–202. [Google Scholar] [CrossRef] [PubMed]
  4. Yang, Y.; Liu, X. Imbalance of Drug Transporter-CYP450s Interplay by Diabetes and Its Clinical Significance. Pharmaceutics 2020, 12, 348. [Google Scholar] [CrossRef] [PubMed]
  5. Guengerich, F.P. Cytochrome P-450 3A4: Regulation and role in drug metabolism. Ann. Rev. Pharm. Toxicol. 1999, 39, 1–17. [Google Scholar] [CrossRef] [PubMed]
  6. Liu, X. ABC Family Transporters. Adv. Exp. Med. Biol. 2019, 1141, 13–100. [Google Scholar] [CrossRef] [PubMed]
  7. Lennernas, H. Clinical pharmacokinetics of atorvastatin. Clin. Pharm. 2003, 42, 1141–1160. [Google Scholar] [CrossRef] [PubMed]
  8. Li, J.; Volpe, D.A.; Wang, Y.; Zhang, W.; Bode, C.; Owen, A.; Hidalgo, I.J. Use of transporter knockdown Caco-2 cells to investigate the in vitro efflux of statin drugs. Drug Metab. Dispos. 2011, 39, 1196–1202. [Google Scholar] [CrossRef] [Green Version]
  9. Ellis, L.C.; Hawksworth, G.M.; Weaver, R.J. ATP-dependent transport of statins by human and rat MRP2/Mrp2. Toxicol. Appl. Pharm. 2013, 269, 187–194. [Google Scholar] [CrossRef]
  10. Kalliokoski, A.; Niemi, M. Impact of OATP transporters on pharmacokinetics. Br. J. Pharm. 2009, 158, 693–705. [Google Scholar] [CrossRef] [Green Version]
  11. Lechner, C.; Reichel, V.; Moenning, U.; Reichel, A.; Fricker, G. Development of a fluorescence-based assay for drug interactions with human Multidrug Resistance Related Protein (MRP2; ABCC2) in MDCKII-MRP2 membrane vesicles. Eur. J. Pharm. Biopharm. 2010, 75, 284–290. [Google Scholar] [CrossRef] [PubMed]
  12. Shitara, Y. Clinical importance of OATP1B1 and OATP1B3 in drug-drug interactions. Drug Metab. Pharm. 2011, 26, 220–227. [Google Scholar] [CrossRef] [PubMed]
  13. Qadir, M.; O’Loughlin, K.L.; Fricke, S.M.; Williamson, N.A.; Greco, W.R.; Minderman, H.; Baer, M.R. Cyclosporin A is a broad-spectrum multidrug resistance modulator. Clin. Cancer Res. 2005, 11, 2320–2326. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Hermann, M.; Asberg, A.; Christensen, H.; Holdaas, H.; Hartmann, A.; Reubsaet, J.L. Substantially elevated levels of atorvastatin and metabolites in cyclosporine-treated renal transplant recipients. Clin. Pharm. 2004, 76, 388–391. [Google Scholar] [CrossRef] [PubMed]
  15. Lemahieu, W.P.; Hermann, M.; Asberg, A.; Verbeke, K.; Holdaas, H.; Vanrenterghem, Y.; Maes, B.D. Combined therapy with atorvastatin and calcineurin inhibitors: No interactions with tacrolimus. Am. J. Transpl. 2005, 5, 2236–2243. [Google Scholar] [CrossRef] [PubMed]
  16. Asberg, A.; Hartmann, A.; Fjeldsa, E.; Bergan, S.; Holdaas, H. Bilateral pharmacokinetic interaction between cyclosporine A and atorvastatin in renal transplant recipients. Am. J. Transpl. 2001, 1, 382–386. [Google Scholar] [CrossRef]
  17. Muck, W.; Mai, I.; Fritsche, L.; Ochmann, K.; Rohde, G.; Unger, S.; Johne, A.; Bauer, S.; Budde, K.; Roots, I.; et al. Increase in cerivastatin systemic exposure after single and multiple dosing in cyclosporine-treated kidney transplant recipients. Clin. Pharm. 1999, 65, 251–261. [Google Scholar] [CrossRef]
  18. Olbricht, C.; Wanner, C.; Eisenhauer, T.; Kliem, V.; Doll, R.; Boddaert, M.; O’Grady, P.; Krekler, M.; Mangold, B.; Christians, U. Accumulation of lovastatin, but not pravastatin, in the blood of cyclosporine-treated kidney graft patients after multiple doses. Clin. Pharm. 1997, 62, 311–321. [Google Scholar] [CrossRef]
  19. Park, J.W.; Siekmeier, R.; Merz, M.; Krell, B.; Harder, S.; Marz, W.; Seidel, D.; Schuler, S.; Gross, W. Pharmacokinetics of pravastatin in heart-transplant patients taking cyclosporin A. Int. J. Clin. Pharm. 2002, 40, 439–450. [Google Scholar] [CrossRef]
  20. Simonson, S.G.; Raza, A.; Martin, P.D.; Mitchell, P.D.; Jarcho, J.A.; Brown, C.D.; Windass, A.S.; Schneck, D.W. Rosuvastatin pharmacokinetics in heart transplant recipients administered an antirejection regimen including cyclosporine. Clin. Pharm. 2004, 76, 167–177. [Google Scholar] [CrossRef]
  21. Park, J.W.; Siekmeier, R.; Lattke, P.; Merz, M.; Mix, C.; Schuler, S.; Jaross, W. Pharmacokinetics and pharmacodynamics of fluvastatin in heart transplant recipients taking cyclosporine A. J. Cardiovasc. Pharm. 2001, 6, 351–361. [Google Scholar] [CrossRef] [PubMed]
  22. Ichimaru, N.; Takahara, S.; Kokado, Y.; Wang, J.D.; Hatori, M.; Kameoka, H.; Inoue, T.; Okuyama, A. Changes in lipid metabolism and effect of simvastatin in renal transplant recipients induced by cyclosporine or tacrolimus. Atherosclerosis 2001, 158, 417–423. [Google Scholar] [CrossRef]
  23. Kajosaari, L.I.; Niemi, M.; Neuvonen, M.; Laitila, J.; Neuvonen, P.J.; Backman, J.T. Cyclosporine markedly raises the plasma concentrations of repaglinide. Clin. Pharm. 2005, 78, 388–399. [Google Scholar] [CrossRef] [PubMed]
  24. Binet, I.; Wallnofer, A.; Weber, C.; Jones, R.; Thiel, G. Renal hemodynamics and pharmacokinetics of bosentan with and without cyclosporine A. Kidney Int. 2000, 57, 224–231. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Saeheng, T.; Na-Bangchang, K.; Karbwang, J. Utility of physiologically based pharmacokinetic (PBPK) modeling in oncology drug development and its accuracy: A systematic review. Eur. J. Clin. Pharm. 2018, 74, 1365–1376. [Google Scholar] [CrossRef]
  26. Wang, Z.; Yang, H.; Xu, J.; Zhao, K.; Chen, Y.; Liang, L.; Li, P.; Chen, N.; Geng, D.; Zhang, X.; et al. Prediction of Atorvastatin Pharmacokinetics in High-Fat Diet and Low-Dose Streptozotocin-Induced Diabetic Rats Using a Semiphysiologically Based Pharmacokinetic Model Involving Both Enzymes and Transporters. Drug Metab. Dispos. 2019, 47, 1066–1079. [Google Scholar] [CrossRef]
  27. Varma, M.V.; Bi, Y.A.; Kimoto, E.; Lin, J. Quantitative prediction of transporter- and enzyme-mediated clinical drug-drug interactions of organic anion-transporting polypeptide 1B1 substrates using a mechanistic net-effect model. J. Pharm. Exp. 2014, 351, 214–223. [Google Scholar] [CrossRef]
  28. Varma, M.V.; Lai, Y.; Feng, B.; Litchfield, J.; Goosen, T.C.; Bergman, A. Physiologically based modeling of pravastatin transporter-mediated hepatobiliary disposition and drug-drug interactions. Pharm. Res. 2012, 29, 2860–2873. [Google Scholar] [CrossRef]
  29. Rodgers, T.; Rowland, M. Physiologically based pharmacokinetic modelling 2: Predicting the tissue distribution of acids, very weak bases, neutrals and zwitterions. J. Pharm. Sci. 2006, 95, 1238–1257. [Google Scholar] [CrossRef]
  30. Fujino, H.; Saito, T.; Tsunenari, Y.; Kojima, J.; Sakaeda, T. Metabolic properties of the acid and lactone forms of HMG-CoA reductase inhibitors. Xenobiotica 2004, 34, 961–971. [Google Scholar] [CrossRef]
  31. Kunze, A.; Poller, B.; Huwyler, J.; Camenisch, G. Application of the extended clearance concept classification system (ECCCS) to predict the victim drug-drug interaction potential of statins. Drug Metab. Pers. 2015, 30, 175–188. [Google Scholar] [CrossRef] [PubMed]
  32. Neuvonen, P.J.; Backman, J.T.; Niemi, M. Pharmacokinetic comparison of the potential over-the-counter statins simvastatin, lovastatin, fluvastatin and pravastatin. Clin. Pharm. 2008, 47, 463–474. [Google Scholar] [CrossRef] [PubMed]
  33. Toda, T.; Eliasson, E.; Ask, B.; Inotsume, N.; Rane, A. Roles of different CYP enzymes in the formation of specific fluvastatin metabolites by human liver microsomes. Basic Clin. Pharm. Toxicol. 2009, 105, 327–332. [Google Scholar] [CrossRef] [PubMed]
  34. Wang, R.W.; Kari, P.H.; Lu, A.Y.; Thomas, P.E.; Guengerich, F.P.; Vyas, K.P. Biotransformation of lovastatin. IV. Identification of cytochrome P450 3A proteins as the major enzymes responsible for the oxidative metabolism of lovastatin in rat and human liver microsomes. Arch. Biochem. Biophys. 1991, 290, 355–361. [Google Scholar] [CrossRef]
  35. Kudo, T.; Goda, H.; Yokosuka, Y.; Tanaka, R.; Komatsu, S.; Ito, K. Estimation of the Contribution of CYP2C8 and CYP3A4 in Repaglinide Metabolism by Human Liver Microsomes Under Various Buffer Conditions. J. Pharm. Sci. 2017, 106, 2847–2852. [Google Scholar] [CrossRef]
  36. Srinivas, N.R. Clinical drug-drug interactions of bosentan, a potent endothelial receptor antagonist, with various drugs: Physiological role of enzymes and transporters. Gen. Physiol. Biophys. 2016, 35, 243–258. [Google Scholar] [CrossRef]
  37. Fahrmayr, C.; Konig, J.; Auge, D.; Mieth, M.; Munch, K.; Segrestaa, J.; Pfeifer, T.; Treiber, A.; Fromm, M. Phase I and II metabolism and MRP2-mediated export of bosentan in a MDCKII-OATP1B1-CYP3A4-UGT1A1-MRP2 quadruple-transfected cell line. Br. J. Pharm. 2013, 169, 21–33. [Google Scholar] [CrossRef] [Green Version]
  38. Davies, B.; Morris, T. Physiological parameters in laboratory animals and humans. Pharm. Res. 1993, 10, 1093–1095. [Google Scholar] [CrossRef]
  39. DeSesso, J.M.; Jacobson, C.F. Anatomical and physiological parameters affecting gastrointestinal absorption in humans and rats. Food Chem. Toxicol. 2001, 39, 209–228. [Google Scholar] [CrossRef]
  40. Perdaems, N.; Blasco, H.; Vinson, C.; Chenel, M.; Whalley, S.; Cazade, F.; Bouzom, F. Predictions of metabolic drug-drug interactions using physiologically based modelling: Two cytochrome P450 3A4 substrates coadministered with ketoconazole or verapamil. Clin. Pharm. 2010, 49, 239–258. [Google Scholar] [CrossRef]
  41. Mikkaichi, T.; Nakai, D.; Yoshigae, Y.; Imaoka, T.; Okudaira, N.; Izumi, T. Liver-selective distribution in rats supports the importance of active uptake into the liver via organic anion transporting polypeptides (OATPs) in humans. Drug Metab. Pharm. 2015, 30, 334–340. [Google Scholar] [CrossRef]
  42. Masuda, N.; Akasaka, I.; Ohtawa, M. Metabolic Fate of Fluvastatin, an Inhibitor of HMG-CoA Reductase (1): Absorption, Distribution and Excretion of [14C] Fluvastatin after Single Administration in Rats. Drug Metab. Pharmacokinet. 1995, 513–528. [Google Scholar]
  43. Liu, Y.; Zeng, B.H.; Shang, H.T.; Cen, Y.Y.; Wei, H. Bama miniature pigs (Sus scrofa domestica) as a model for drug evaluation for humans: Comparison of in vitro metabolism and in vivo pharmacokinetics of lovastatin. Comp. Med. 2008, 58, 580–587. [Google Scholar] [PubMed]
  44. Modry, D.L.; Stinson, E.B.; Oyer, P.E.; Jamieson, S.W.; Baldwin, J.C.; Shumway, N.E. Acute rejection and massive cyclosporine requirements in heart transplant recipients treated with rifampin. Transplantation 1985, 39, 313–314. [Google Scholar] [CrossRef] [PubMed]
  45. Germershausen, J.I.; Hunt, V.M.; Bostedor, R.G.; Bailey, P.J.; Karkas, J.D.; Alberts, A.W. Tissue selectivity of the cholesterol-lowering agents lovastatin, simvastatin and pravastatin in rats in vivo. Biochem. Biophys. Res. Commun. 1989, 158, 667–675. [Google Scholar] [CrossRef]
  46. Camenisch, G.; Umehara, K. Predicting human hepatic clearance from in vitro drug metabolism and transport data: A scientific and pharmaceutical perspective for assessing drug-drug interactions. Biopharm. Drug Dispos. 2012, 33, 179–194. [Google Scholar] [CrossRef]
  47. Varma, M.V.; Lin, J.; Bi, Y.A.; Kimoto, E.; Rodrigues, A.D. Quantitative Rationalization of Gemfibrozil Drug Interactions: Consideration of Transporters-Enzyme Interplay and the Role of Circulating Metabolite Gemfibrozil 1-O-beta-Glucuronide. Drug Metab. Dispos. 2015, 43, 1108–1118. [Google Scholar] [CrossRef] [Green Version]
  48. Kato, M.; Shitara, Y.; Sato, H.; Yoshisue, K.; Hirano, M.; Ikeda, T.; Sugiyama, Y. The quantitative prediction of CYP-mediated drug interaction by physiologically based pharmacokinetic modeling. Pharm. Res. 2008, 25, 1891–1901. [Google Scholar] [CrossRef]
  49. Gertz, M.; Harrison, A.; Houston, J.B.; Galetin, A. Prediction of human intestinal first-pass metabolism of 25 CYP3A substrates from in vitro clearance and permeability data. Drug Metab. Dispos. 2010, 38, 1147–1158. [Google Scholar] [CrossRef]
  50. Kim, S.J.; Toshimoto, K.; Yao, Y.; Yoshikado, T.; Sugiyama, Y. Quantitative Analysis of Complex Drug-Drug Interactions Between Repaglinide and Cyclosporin A/Gemfibrozil Using Physiologically Based Pharmacokinetic Models With In Vitro Transporter/Enzyme Inhibition Data. J. Pharm. Sci. 2017, 106, 2715–2726. [Google Scholar] [CrossRef] [Green Version]
  51. Sato, M.; Toshimoto, K.; Tomaru, A.; Yoshikado, T.; Tanaka, Y.; Hisaka, A.; Lee, W.; Sugiyama, Y. Physiologically Based Pharmacokinetic Modeling of Bosentan Identifies the Saturable Hepatic Uptake As a Major Contributor to Its Nonlinear Pharmacokinetics. Drug Metab. Dispos. 2018, 46, 740–748. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  52. Cetkovic, Z.; Cvijic, S.; Vasiljevic, D. In vitro/in silico approach in the development of simvastatin-loaded self-microemulsifying drug delivery systems. Drug Dev. Ind. Pharm. 2018, 44, 849–860. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Gertz, M.; Houston, J.B.; Galetin, A. Physiologically based pharmacokinetic modeling of intestinal first-pass metabolism of CYP3A substrates with high intestinal extraction. Drug Metab. Dispos. 2011, 39, 1633–1642. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Morse, B.L.; Alberts, J.J.; Posada, M.M.; Rehmel, J.; Kolur, A.; Tham, L.S.; Loghin, C.; Hillgren, K.M.; Hall, S.D.; Dickinson, G.L. Physiologically-Based Pharmacokinetic Modeling of Atorvastatin Incorporating Delayed Gastric Emptying and Acid-to-Lactone Conversion. Cpt. Pharmacomet. Syst. Pharm. 2019, 8, 664–675. [Google Scholar] [CrossRef] [PubMed]
  55. Duan, P.; Zhao, P.; Zhang, L. Physiologically Based Pharmacokinetic (PBPK) Modeling of Pitavastatin and Atorvastatin to Predict Drug-Drug Interactions (DDIs). Eur. J. Drug Metab. Pharm. 2017, 42, 689–705. [Google Scholar] [CrossRef]
  56. Jones, H.M.; Barton, H.A.; Lai, Y.; Bi, Y.A.; Kimoto, E.; Kempshall, S.; Tate, S.C.; El-Kattan, A.; Houston, J.B.; Galetin, A.; et al. Mechanistic pharmacokinetic modeling for the prediction of transporter-mediated disposition in humans from sandwich culture human hepatocyte data. Drug Metab. Dispos. 2012, 40, 1007–1017. [Google Scholar] [CrossRef] [Green Version]
  57. Nishimuta, H.; Sato, K.; Yabuki, M.; Komuro, S. Prediction of the intestinal first-pass metabolism of CYP3A and UGT substrates in humans from in vitro data. Drug Metab. Pharm. 2011, 26, 592–601. [Google Scholar] [CrossRef] [Green Version]
  58. Varma, M.V.; Lai, Y.; Kimoto, E.; Goosen, T.C.; El-Kattan, A.F.; Kumar, V. Mechanistic modeling to predict the transporter- and enzyme-mediated drug-drug interactions of repaglinide. Pharm. Res. 2013, 30, 1188–1199. [Google Scholar] [CrossRef]
  59. Wu, X.; Whitfield, L.R.; Stewart, B.H. Atorvastatin transport in the Caco-2 cell model: Contributions of P-glycoprotein and the proton-monocarboxylic acid co-transporter. Pharm. Res. 2000, 17, 209–215. [Google Scholar] [CrossRef]
  60. Kivisto, K.T.; Zukunft, J.; Hofmann, U.; Niemi, M.; Rekersbrink, S.; Schneider, S.; Luippold, G.; Schwab, M.; Eichelbaum, M.; Fromm, M.F. Characterisation of cerivastatin as a P-glycoprotein substrate: Studies in P-glycoprotein-expressing cell monolayers and mdr1a/b knock-out mice. Naunyn Schmiedebergs Arch. Pharm. 2004, 370, 124–130. [Google Scholar] [CrossRef]
  61. Li, J.; Wang, Y.; Zhang, W.; Huang, Y.; Hein, K.; Hidalgo, I.J. The role of a basolateral transporter in rosuvastatin transport and its interplay with apical breast cancer resistance protein in polarized cell monolayer systems. Drug Metab. Dispos. 2012, 40, 2102–2108. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  62. Winiwarter, S.; Bonham, N.M.; Ax, F.; Hallberg, A.; Lennernas, H.; Karlen, A. Correlation of human jejunal permeability (in vivo) of drugs with experimentally and theoretically derived parameters. A multivariate data analysis approach. J. Med. Chem. 1998, 41, 4939–4949. [Google Scholar] [CrossRef] [PubMed]
  63. Yaghoobian, M.; Haeri, A.; Bolourchian, N.; Shahhosseini, S.; Dadashzadeh, S. An Investigation into the Role of P-Glycoprotein in the Intestinal Absorption of Repaglinide: Assessed by Everted Gut Sac and Caco-2 Cell Line. Iran J. Pharm. Res. 2019, 18, 102–110. [Google Scholar] [PubMed]
  64. Guo, T.; Qin, H.-x.; Xia, D.-y. Pharmacokinetics and Bioequivalence of Three Preparations of Cyclosporine A in Chinese Healthy Volunteers. Chin. Pharmaceut. J. Beijing 2008, 43, 535. [Google Scholar]
  65. He, J.-C.; Feng, E.-F.; Zhang, Q.; Yin, Y.-Q.; He, H.-J.; Xu, G.-L. Study on bioequivalence of ciclosporin microemulsion oral solution in healthy volunteers. Pharm. Care Res. 2009, 9, 137–139. [Google Scholar]
  66. de la Pena, A.; Cui, X.; Geiser, J.; Loghin, C. No Dose Adjustment is Recommended for Digoxin, Warfarin, Atorvastatin or a Combination Oral Contraceptive When Coadministered with Dulaglutide. Clin. Pharm. 2017, 56, 1415–1427. [Google Scholar] [CrossRef]
  67. Guo, C.X.; Pei, Q.; Yin, J.Y.; Peng, X.D.; Zhou, B.T.; Zhao, Y.C.; Wu, L.X.; Meng, X.G.; Wang, G.; Li, Q.; et al. Effects of Ginkgo biloba extracts on pharmacokinetics and efficacy of atorvastatin based on plasma indices. Xenobiotica 2012, 42, 784–790. [Google Scholar] [CrossRef]
  68. Hausner, H.; Derving Karsbol, J.; Holst, A.G.; Jacobsen, J.B.; Wagner, F.D.; Golor, G.; Anderson, T.W. Effect of Semaglutide on the Pharmacokinetics of Metformin, Warfarin, Atorvastatin and Digoxin in Healthy Subjects. Clin. Pharm. 2017, 56, 1391–1401. [Google Scholar] [CrossRef] [Green Version]
  69. Hoch, M.; Hoever, P.; Theodor, R.; Dingemanse, J. Almorexant effects on CYP3A4 activity studied by its simultaneous and time-separated administration with simvastatin and atorvastatin. Eur J. Clin. Pharm. 2013, 69, 1235–1245. [Google Scholar] [CrossRef]
  70. Kubitza, D.; Becka, M.; Roth, A.; Mueck, W. Absence of clinically relevant interactions between rivaroxaban—An oral, direct Factor Xa inhibitor—And digoxin or atorvastatin in healthy subjects. J. Int Med Res. 2012, 40, 1688–1707. [Google Scholar] [CrossRef] [Green Version]
  71. DeGorter, M.K.; Urquhart, B.L.; Gradhand, U.; Tirona, R.G.; Kim, R.B. Disposition of atorvastatin, rosuvastatin, and simvastatin in oatp1b2-/- mice and intraindividual variability in human subjects. J. Clin. Pharm. 2012, 52, 1689–1697. [Google Scholar] [CrossRef]
  72. Blonk, M.; van Beek, M.; Colbers, A.; Schouwenberg, B.; Burger, D. Pharmacokinetic Drug-Drug Interaction Study Between Raltegravir and Atorvastatin 20 mg in Healthy Volunteers. J. Acquir. Immune Defic. Syndr. 2015, 69, 44–51. [Google Scholar] [CrossRef] [PubMed]
  73. Khalilieh, S.; Yee, K.L.; Sanchez, R.I.; Triantafyllou, I.; Fan, L.; Maklad, N.; Jordan, H.; Martell, M.; Iwamoto, M. Results of a Doravirine-Atorvastatin Drug-Drug Interaction Study. Antimicrob. Agents Chemother. 2017, 61. [Google Scholar] [CrossRef] [Green Version]
  74. Ayalasomayajula, S.; Pan, W.; Han, Y.; Yang, F.; Langenickel, T.; Pal, P.; Zhou, W.; Yuan, Y.; Rajman, I.; Sunkara, G. Assessment of Drug-Drug Interaction Potential Between Atorvastatin and LCZ696, A Novel Angiotensin Receptor Neprilysin Inhibitor, in Healthy Chinese Male Subjects. Eur. J. Drug Metab. Pharm. 2017, 42, 309–318. [Google Scholar] [CrossRef] [PubMed]
  75. Shin, K.H.; Kim, T.E.; Kim, S.E.; Lee, M.G.; Song, I.S.; Yoon, S.H.; Cho, J.Y.; Jang, I.J.; Shin, S.G.; Yu, K.S. The effect of the newly developed angiotensin receptor II antagonist fimasartan on the pharmacokinetics of atorvastatin in relation to OATP1B1 in healthy male volunteers. J. Cardiovasc. Pharm. 2011, 58, 492–499. [Google Scholar] [CrossRef] [PubMed]
  76. Muck, W. Clinical pharmacokinetics of cerivastatin. Clin. Pharm. 2000, 39, 99–116. [Google Scholar] [CrossRef]
  77. Backman, J.T.; Kyrklund, C.; Neuvonen, M.; Neuvonen, P.J. Gemfibrozil greatly increases plasma concentrations of cerivastatin. Clin. Pharm. 2002, 72, 685–691. [Google Scholar] [CrossRef]
  78. Bauer, S.; Mwinyi, J.; Stoeckle, A.; Gerloff, T.; Roots, I. Quantification of pravastatin in human plasma and urine after solid phase extraction using high performance liquid chromatography with ultraviolet detection. J. Chromatogr. B Anal. Technol. Biomed. Life Sci. 2005, 818, 257–262. [Google Scholar] [CrossRef]
  79. Neuvonen, P.J.; Kantola, T.; Kivisto, K.T. Simvastatin but not pravastatin is very susceptible to interaction with the CYP3A4 inhibitor itraconazole. Clin. Pharm. 1998, 63, 332–341. [Google Scholar] [CrossRef]
  80. Kantola, T.; Backman, J.T.; Niemi, M.; Kivisto, K.T.; Neuvonen, P.J. Effect of fluconazole on plasma fluvastatin and pravastatin concentrations. Eur. J. Clin. Pharm. 2000, 56, 225–229. [Google Scholar] [CrossRef]
  81. Cooper, K.J.; Martin, P.D.; Dane, A.L.; Warwick, M.J.; Schneck, D.W.; Cantarini, M.V. Effect of itraconazole on the pharmacokinetics of rosuvastatin. Clin. Pharm. 2003, 73, 322–329. [Google Scholar] [CrossRef]
  82. Keskitalo, J.E.; Kurkinen, K.J.; Neuvonen, M.; Backman, J.T.; Neuvonen, P.J.; Niemi, M. No significant effect of ABCB1 haplotypes on the pharmacokinetics of fluvastatin, pravastatin, lovastatin, and rosuvastatin. Br. J. Clin. Pharm. 2009, 68, 207–213. [Google Scholar] [CrossRef] [Green Version]
  83. Kalafsky, G.; Smith, H.T.; Choc, M.G. High-performance liquid chromatographic method for the determination of fluvastatin in human plasma. J. Chromatogr. 1993, 614, 307–313. [Google Scholar] [CrossRef] [PubMed]
  84. Ayalasomayajula, S.; Han, Y.; Langenickel, T.; Malcolm, K.; Zhou, W.; Hanna, I.; Alexander, N.; Natrillo, A.; Goswami, B.; Hinder, M.; et al. In vitro and clinical evaluation of OATP-mediated drug interaction potential of sacubitril/valsartan (LCZ696). J. Clin. Pharm. 2016, 41, 424–431. [Google Scholar] [CrossRef]
  85. Harvey, R.D.; Aransay, N.R.; Isambert, N.; Lee, J.S.; Arkenau, T.; Vansteenkiste, J.; Dickinson, P.A.; Bui, K.; Weilert, D.; So, K.; et al. Effect of multiple-dose osimertinib on the pharmacokinetics of simvastatin and rosuvastatin. Br. J. Clin. Pharm. 2018, 84, 2877–2888. [Google Scholar] [CrossRef] [Green Version]
  86. Chu, N.N.; Chen, W.L.; Xu, H.R.; Li, X.N. Pharmacokinetics and safety of ezetimibe/simvastatin combination tablet: An open-label, single-dose study in healthy Chinese subjects. Clin. Drug Investig. 2012, 32, 791–798. [Google Scholar] [CrossRef]
  87. Yin, O.Q.; Mak, V.W.; Hu, M.; Fok, B.S.; Chow, M.S.; Tomlinson, B. Impact of CYP2D6 polymorphisms on the pharmacokinetics of lovastatin in Chinese subjects. Eur. J. Clin. Pharm. 2012, 68, 943–949. [Google Scholar] [CrossRef]
  88. Tornio, A.; Vakkilainen, J.; Neuvonen, M.; Backman, J.T.; Neuvonen, P.J.; Niemi, M. SLCO1B1 polymorphism markedly affects the pharmacokinetics of lovastatin acid. Pharm. Genom. 2015, 25, 382–387. [Google Scholar] [CrossRef]
  89. Jia, M.M.; Zhou, Y.; He, X.M.; Wu, Y.L.; Li, H.Q.; Chen, H.; Li, W.Y. Development of a liquid chromatography-tandem mass spectrometry method for determination of butoconazole nitrate in human plasma and its application to a pharmacokinetic study. J. Huazhong Univ. Sci. Technol. Med. Sci. 2014, 34, 431–436. [Google Scholar] [CrossRef]
  90. Weber, C.; Schmitt, R.; Birnboeck, H.; Hopfgartner, G.; van Marle, S.P.; Peeters, P.A.; Jonkman, J.H.; Jones, C.R. Pharmacokinetics and pharmacodynamics of the endothelin-receptor antagonist bosentan in healthy human subjects. Clin. Pharm. 1996, 60, 124–137. [Google Scholar] [CrossRef]
  91. Li, M.; de Graaf, I.A.; de Jager, M.H.; Groothuis, G.M. P-gp activity and inhibition in the different regions of human intestine ex vivo. Biopharm. Drug Dispos. 2017, 38, 127–138. [Google Scholar] [CrossRef] [PubMed]
  92. Goldberg, R.; Roth, D. Evaluation of fluvastatin in the treatment of hypercholesterolemia in renal transplant recipients taking cyclosporine. Transplantation 1996, 62, 1559–1564. [Google Scholar] [CrossRef] [PubMed]
  93. Arnadottir, M.; Eriksson, L.O.; Thysell, H.; Karkas, J.D. Plasma concentration profiles of simvastatin 3-hydroxy-3-methyl-glutaryl-coenzyme A reductase inhibitory activity in kidney transplant recipients with and without ciclosporin. Nephron 1993, 65, 410–413. [Google Scholar] [CrossRef] [PubMed]
  94. Campana, C.; Iacona, I.; Regazzi, M.B.; Gavazzi, A.; Perani, G.; Raddato, V.; Montemartini, C.; Vigano, M. Efficacy and pharmacokinetics of simvastatin in heart transplant recipients. Ann. Pharm. 1995, 29, 235–239. [Google Scholar] [CrossRef]
  95. Gullestad, L.; Nordal, K.P.; Berg, K.J.; Cheng, H.; Schwartz, M.S.; Simonsen, S. Interaction between lovastatin and cyclosporine A after heart and kidney transplantation. Transpl. Proc. 1999, 31, 2163–2165. [Google Scholar] [CrossRef]
  96. Gullestad, L.; Nordal, K.P.; Forfang, K.; Ihlen, H.; Hostmark, A.; Berg, K.J.; Cheng, H.; Schwartz, M.S.; Geiran, O.; Simonsen, S. Post-transplant hyperlipidaemia: Low-dose lovastatin lowers atherogenic lipids without plasma accumulation of lovastatin. J. Intern Med. 1997, 242, 483–490. [Google Scholar] [CrossRef] [PubMed]
  97. Izumi, S.; Nozaki, Y.; Komori, T.; Maeda, K.; Takenaka, O.; Kusano, K.; Yoshimura, T.; Kusuhara, H.; Sugiyama, Y. Substrate-dependent inhibition of organic anion transporting polypeptide 1B1: Comparative analysis with prototypical probe substrates estradiol-17beta-glucuronide, estrone-3-sulfate, and sulfobromophthalein. Drug Metab. Dispos. 2013, 41, 1859–1866. [Google Scholar] [CrossRef] [Green Version]
  98. Harwood, M.D.; Neuhoff, S.; Carlson, G.L.; Warhurst, G.; Rostami-Hodjegan, A. Absolute abundance and function of intestinal drug transporters: A prerequisite for fully mechanistic in vitro-in vivo extrapolation of oral drug absorption. Biopharm. Drug Dispos. 2013, 34, 2–28. [Google Scholar] [CrossRef]
  99. Maeda, K.; Ikeda, Y.; Fujita, T.; Yoshida, K.; Azuma, Y.; Haruyama, Y.; Yamane, N.; Kumagai, Y.; Sugiyama, Y. Identification of the rate-determining process in the hepatic clearance of atorvastatin in a clinical cassette microdosing study. Clin. Pharm. 2011, 90, 575–581. [Google Scholar] [CrossRef]
  100. Tuegel, C.; Bansal, N. Heart failure in patients with kidney disease. Heart 2017, 103, 1848–1853. [Google Scholar] [CrossRef]
  101. Lee, H.H.; Ho, R.H. Interindividual and interethnic variability in drug disposition: Polymorphisms in organic anion transporting polypeptide 1B1 (OATP1B1; SLCO1B1). Br. J. Clin. Pharm. 2017, 83, 1176–1184. [Google Scholar] [CrossRef] [Green Version]
  102. Guo, H.; Liu, C.; Li, J.; Zhang, M.; Hu, M.; Xu, P.; Liu, L.; Liu, X. A mechanistic physiologically based pharmacokinetic-enzyme turnover model involving both intestine and liver to predict CYP3A induction-mediated drug-drug interactions. J. Pharm. Sci. 2013, 102, 2819–2836. [Google Scholar] [CrossRef] [PubMed]
  103. Englund, G.; Rorsman, F.; Ronnblom, A.; Karlbom, U.; Lazorova, L.; Grasjo, J.; Kindmark, A.; Artursson, P. Regional levels of drug transporters along the human intestinal tract: Co-expression of ABC and SLC transporters and comparison with Caco-2 cells. Eur. J. Pharm. Sci. 2006, 29, 269–277. [Google Scholar] [CrossRef]
  104. Peters, S.A.; Jones, C.R.; Ungell, A.L.; Hatley, O.J. Predicting Drug Extraction in the Human Gut Wall: Assessing Contributions from Drug Metabolizing Enzymes and Transporter Proteins using Preclinical Models. Clin. Pharm. 2016, 55, 673–696. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  105. Paine, M.F.; Khalighi, M.; Fisher, J.M.; Shen, D.D.; Kunze, K.L.; Marsh, C.L.; Perkins, J.D.; Thummel, K.E. Characterization of interintestinal and intraintestinal variations in human CYP3A-dependent metabolism. J. Pharm. Exp. 1997, 283, 1552–1562. [Google Scholar]
  106. Mathialagan, S.; Piotrowski, M.A.; Tess, D.A.; Feng, B.; Litchfield, J.; Varma, M.V. Quantitative Prediction of Human Renal Clearance and Drug-Drug Interactions of Organic Anion Transporter Substrates Using In Vitro Transport Data: A Relative Activity Factor Approach. Drug Metab. Dispos. 2017, 45, 409–417. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  107. Gertz, M.; Cartwright, C.M.; Hobbs, M.J.; Kenworthy, K.E.; Rowland, M.; Houston, J.B.; Galetin, A. Cyclosporine inhibition of hepatic and intestinal CYP3A4, uptake and efflux transporters: Application of PBPK modeling in the assessment of drug-drug interaction potential. Pharm. Res. 2013, 30, 761–780. [Google Scholar] [CrossRef]
Ijms 21 07023 g001 550
Figure 1. Predicted plasma concentrations (line) of the tested agents using the PBPK model and observed plasma concentrations (points) following a single dose to healthy subjects. Different points represented different dosages for drugs in the corresponding figure. (A) CsA; (B) atorvastatin; (C) cerivastatin; (D) pravastatin; (E) rosuvastatin; (F) fluvastatin; (G) simvastatin; (H) lovastatin; (I) repaglinide; (J) bosentan. Observed data were cited from reports [10,20,23,64,65,66,67,68,69,70,71,72,73,74,75,76,77,78,79,80,81,82,83,84,85,86,87,88,89,90].
Figure 1. Predicted plasma concentrations (line) of the tested agents using the PBPK model and observed plasma concentrations (points) following a single dose to healthy subjects. Different points represented different dosages for drugs in the corresponding figure. (A) CsA; (B) atorvastatin; (C) cerivastatin; (D) pravastatin; (E) rosuvastatin; (F) fluvastatin; (G) simvastatin; (H) lovastatin; (I) repaglinide; (J) bosentan. Observed data were cited from reports [10,20,23,64,65,66,67,68,69,70,71,72,73,74,75,76,77,78,79,80,81,82,83,84,85,86,87,88,89,90].
Ijms 21 07023 g001
Ijms 21 07023 g002 550
Figure 2. Predicted plasma concentrations (line) of the victim drugs agents using the PBPK model and observed plasma concentrations (points) following oral administration to subjects when co-administered with CsA (twice daily). (A) atorvastatin; (B) cerivastatin; (C) pravastatin; (D) rosuvastatin; (E) fluvastatin; (F) simvastatin; (G) lovastatin; (H) repaglinide; (I) bosentan. Observed data were cited from clinic reports [14,15,16,17,18,19,20,21,23,24,92,93,94,95,96].
Figure 2. Predicted plasma concentrations (line) of the victim drugs agents using the PBPK model and observed plasma concentrations (points) following oral administration to subjects when co-administered with CsA (twice daily). (A) atorvastatin; (B) cerivastatin; (C) pravastatin; (D) rosuvastatin; (E) fluvastatin; (F) simvastatin; (G) lovastatin; (H) repaglinide; (I) bosentan. Observed data were cited from clinic reports [14,15,16,17,18,19,20,21,23,24,92,93,94,95,96].
Ijms 21 07023 g002
Ijms 21 07023 g003 550
Figure 3. Contributions of (A) inhibition constant on hepatic organic anion transporting polypeptides (OATP) by CsA (Ki,OATP), (B) OATP-mediated hepatic uptake (CLint,up), (C) CsA dosage (100 mg, 200 mg and 400 mg), (D) transition rate (Kt) and (E) empirical scaling factor (SF) to plasma concentrations of atorvastatin following oral 40 mg of atorvastatin. CsA dosage was set to be 200 mg/d. Contributions of (F) inhibitions on CYP3A4 (non-CYP3A4), integrated effects of hepatic CYP3A-mediated metabolism/hepatic OATP-mediated uptake (non-liver) and integrated effects of intestinal CYP3A4-mediated metabolism/transporter-mediated intestinal efflux (non-gut) to the increased plasma concentrations of atorvastatin by CsA.
Figure 3. Contributions of (A) inhibition constant on hepatic organic anion transporting polypeptides (OATP) by CsA (Ki,OATP), (B) OATP-mediated hepatic uptake (CLint,up), (C) CsA dosage (100 mg, 200 mg and 400 mg), (D) transition rate (Kt) and (E) empirical scaling factor (SF) to plasma concentrations of atorvastatin following oral 40 mg of atorvastatin. CsA dosage was set to be 200 mg/d. Contributions of (F) inhibitions on CYP3A4 (non-CYP3A4), integrated effects of hepatic CYP3A-mediated metabolism/hepatic OATP-mediated uptake (non-liver) and integrated effects of intestinal CYP3A4-mediated metabolism/transporter-mediated intestinal efflux (non-gut) to the increased plasma concentrations of atorvastatin by CsA.
Ijms 21 07023 g003
Ijms 21 07023 g004 550
Figure 4. Effects of (A) hepatic CYP3A4 (non-liver CYP3A4), hepatic BCRP (non-liver BCRP), intestinal CYP3A4 (non-gut CYP3A4), intestinal P-gp (non-gut P-gp), intestinal BCRP (non-gut-BCRP), intestinal MRP2 (non-gut-MRP2) and integrated effects of the intestine (intestinal CYP3A4, intestinal P-gp, intestinal MRP2 and intestinal BCRP, non-gut) to plasma concentrations of atorvastatin following oral 40 mg of atorvastatin, which was compared with actual concentrations (Whole), (B) hepatic OATPs (non-liver-OATP).
Figure 4. Effects of (A) hepatic CYP3A4 (non-liver CYP3A4), hepatic BCRP (non-liver BCRP), intestinal CYP3A4 (non-gut CYP3A4), intestinal P-gp (non-gut P-gp), intestinal BCRP (non-gut-BCRP), intestinal MRP2 (non-gut-MRP2) and integrated effects of the intestine (intestinal CYP3A4, intestinal P-gp, intestinal MRP2 and intestinal BCRP, non-gut) to plasma concentrations of atorvastatin following oral 40 mg of atorvastatin, which was compared with actual concentrations (Whole), (B) hepatic OATPs (non-liver-OATP).
Ijms 21 07023 g004
Ijms 21 07023 g005 550
Figure 5. Schematic diagram of the whole PBPK model involving both enzyme and transporter turnover. Qi indicated blood flow in the corresponding compartment. CL is the intrinsic clearance. Gwi is the the i-th enterocyte compartment. Kti, ka,i and kb,i represented the transit rate constant, drug absorption rate constant and efflux rate constant from enterocyte to gut lumen, respectively. CLint, up, CLint, pd, CLint, bile, CLint, gwi and CLint,met were uptake clearance, passive clearance, efflux clearance to bile, metabolism clearance in the gut and metabolic clearance in hepatocytes.
Figure 5. Schematic diagram of the whole PBPK model involving both enzyme and transporter turnover. Qi indicated blood flow in the corresponding compartment. CL is the intrinsic clearance. Gwi is the the i-th enterocyte compartment. Kti, ka,i and kb,i represented the transit rate constant, drug absorption rate constant and efflux rate constant from enterocyte to gut lumen, respectively. CLint, up, CLint, pd, CLint, bile, CLint, gwi and CLint,met were uptake clearance, passive clearance, efflux clearance to bile, metabolism clearance in the gut and metabolic clearance in hepatocytes.
Ijms 21 07023 g005
Table
Table 1. Physiological parameters for humans used in the physiologically based pharmacokinetic PBPK model.
Table 1. Physiological parameters for humans used in the physiologically based pharmacokinetic PBPK model.
TissueV (mL) [38,39]Q (mL/min) [38]Transit Rate Constant (min−1) [40]
Lung11705600/
Kidney2801240/
Heart310240/
Liver16901518/
Muscle35,000750/
Skin7800300/
Brain1450700/
Adipose10,000260/
Rob5100592/
Spleen19080/
Artery17305600/
Venous34705600/
Stomach160380.0462
Duodenum701180.0462
Jejunum2094130.012
Ileum1392440.0058
Cecum116440.004
Colon11162810.0013
Table
Table 2. Ratio of drug concentration in tissue to plasma for cyclosporin A (CsA), atorvastatin (Ato), cerivastatin (Cer), pravastatin (Pra), rosuvastatin (Ros), fluvastatin (Flu), simvastatin (Sim), bosentan (Bos), repaglinide (Rep) and lovastatin (Lov).
Table 2. Ratio of drug concentration in tissue to plasma for cyclosporin A (CsA), atorvastatin (Ato), cerivastatin (Cer), pravastatin (Pra), rosuvastatin (Ros), fluvastatin (Flu), simvastatin (Sim), bosentan (Bos), repaglinide (Rep) and lovastatin (Lov).
TissueCsA aAto aCer [41]Pra aRos [41]Flu [42]Sim [41]Lov [43]Rep [41]Bos [41]
Adipose253.351810.43 a0.140.820.3390.39298.060.20.52
Liver4.13 [44]4.575156 [45]10.6034.1322.725.6315.851
Muscle3.672.640.20.210.200.1470.2194.690.3020.23
Lung7.755.61.130.360.400.4070.4269.060.6041.01
Kidney5.794.182.9348.5 [45]1.121.592.874.251.481.14
Brain11.778.440.20.220.010.0940.2070.880.0780.35
Heart5.223.7820.320.400.820.61211.630.7020.64
Intestine5.99 [44]8.970.2 a0.210.18 a8.220.2 a36.20.07 a0.47 a
Skin18.2813.160.5330.380.380.3130.3921.460.2090.48
Spleen3.362.430.7330.260.230.230.2920.310.3460.52
Stomach3.672.640.090.210.20 a1.2220.2 a62.690.07 a0.16 a
Rest5.2410.010.110.010.010.010.010.010.01
fub b0.05 [27]0.084 [46]0.017 [47]0.56 [28]0.174 [27]0.014 [27]0.107 [48]0.03 [49] 0.011 [50]0.02 [51]
Rb b1.36 [27]0.61 [27]0.76 [47]0.839 [28]0.69 [27]0.57 [27]0.56 [52]0.57 [53]0.62 [50]0.6 [51]
a calculated from PKa and logP [29]. b fub is unbound fraction of drug in blood. Rb is blood/plasma ratio.
Table
Table 3. Intrinsic metabolic/transport parameters and empirical scaling factor (SF) of cyclosporin A (CsA), atorvastatin (Ato), cerivastatin (Cer), pravastatin (Pra), rosuvastatin (Ros), fluvastatin (Flu), simvastatin (Sim), lovastatin (Lov), repaglinide (Rep) and bosentan (Bos) in the liver.
Table 3. Intrinsic metabolic/transport parameters and empirical scaling factor (SF) of cyclosporin A (CsA), atorvastatin (Ato), cerivastatin (Cer), pravastatin (Pra), rosuvastatin (Ros), fluvastatin (Flu), simvastatin (Sim), lovastatin (Lov), repaglinide (Rep) and bosentan (Bos) in the liver.
DrugCLint metCLint UptakeSFCLint, pdCLint, bile
CYP3A4Other P-gpBCRPMRP2
mL/minmL/min mL/minmL/min
CsA5432 [46]/10,857 [46]0.1 a2933 [46]637 [46]//
Ato3469.5 [27]612.3 [27]6374.7 [54]32.6 [27]3916.79 [46]/302.4 [55]/
Cer979.8 [27]1197.6 [27]1942.8 [56]12.5 [27]3541.5 [27]40.5 [27]//
Pra//283.3 [27]19.4 [27]80.9 [27]//80.9 [27]
Ros//1841.6 [27]9.2 [27]242.8 [27]/566.6 [27]/
Flu2386.3 [56]2801.3 [56]9106.6 [56]21 a4047.4 [56]/3440 [56]/
Sim231,391 [57]/9234.1 [31]25 a23,695.5 [31]101.2 [31]//
Lov6893.96 [43]/13,129.8 [31]4 a11,569.5 [31]///
Rep2503.7 [27]6438 [27]7184 [27]16.9 [58]4452.14 [27]/20.24 [27]/
Bos1365.14 [27]/7184 [27]1.1 [27]2023.7 [27]1472.6 [51]//
a SF is the active uptake scaling factor from the simulation whose source has been shown in Table S1 in Supplementary Material.
Table
Table 4. Intrinsic metabolic/transport parameters of cyclosporin A (CsA), atorvastatin (Ato), cerivastatin (Cer), pravastatin (Pra), rosuvastatin (Ros), fluvastatin (Flu), simvastatin (Sim), lovastatin (Lov), repaglinide (Rep) and bosentan (Bos) in the intestine.
Table 4. Intrinsic metabolic/transport parameters of cyclosporin A (CsA), atorvastatin (Ato), cerivastatin (Cer), pravastatin (Pra), rosuvastatin (Ros), fluvastatin (Flu), simvastatin (Sim), lovastatin (Lov), repaglinide (Rep) and bosentan (Bos) in the intestine.
DrugkaPeff, A-BPeff, B-ACLint, gut [53]
min −1cm/min cm/min mL/min
P-gpBCRPMRP2DuodenumJejunumIleum
CsA0.025 [28]/0.01 [28]//11.2544.625.98
Ato/0.0094 [59]0.0185 [8]0.0118 [8]0.011 [8]12.248.27328.22
Cer/0.0133 [47]0.0099 [60]0.0065 [60]0.0118 [60]///
Pra0.021 [28]///43.7 a///
Ros0.0022 [61]/0.012 [8]0.013 [8]0.011 [8]///
Flu/0.038 [62]0.035 [8]0.026 [8]0.029 [8]///
Sim/0.01 [49]0.0093 [49]//755.22991.31745
Lov/0.016 [49] /0.014 [49]87189.82011
Rep0.04 [58] 0.0148 [63] ///
Bos/0.014 [51] 5.84 b [51]///
a The unit of value is mL/min for the efflux of MRP2. b The value with the unit mL/min was calculated from Vmax and Km.
Table
Table 5. The predicted (Pred) pharmacokinetic parameters and the mean clinic observations (Obs) of cyclosporin A (CsA), atorvastatin (Ato), cerivastatin (Cer), pravastatin (Pra), rosuvastatin (Ros), fluvastatin (Flu), simvastatin (Sim), lovastatin (Lov), repaglinide (Rep) and bosentan (Bos).
Table 5. The predicted (Pred) pharmacokinetic parameters and the mean clinic observations (Obs) of cyclosporin A (CsA), atorvastatin (Ato), cerivastatin (Cer), pravastatin (Pra), rosuvastatin (Ros), fluvastatin (Flu), simvastatin (Sim), lovastatin (Lov), repaglinide (Rep) and bosentan (Bos).
DrugRef.DoseCmaxPred/ObsAUCPred/ObsAAFE
mgng/mL ng·h/mL
PredObs PredObs
CsA[64]200681.17947.160.72 5456.127113.470.77 1.5
[64]200681.17758.860.90 5456.126771.870.81 1.5
[64]200681.17938.710.73 5456.127253.950.75 1.6
[65]5001702.9324800.69 11,984.716,1000.74 1.7
Ato[66]4014.5919.50.75 71.6982.80.87 2.3
[67]4014.5920.220.72 71.7787.480.82 1.8
[68]4013.59.451.43 73.1889.660.82 2.1
[69]4014.5911.61.26 66.7069.40.96 1.6
[70]206.754.71.44 35.7847.40.75 2.5
[71]207.298.90.82 32.6534.50.95 1.7
[72]207.299.280.79 34.2243.70.78 1.7
[73]206.757.870.86 37.2539.90.93 1.7
[74]8029.1838.070.77 136.39182.960.75 1.5
[75]8029.1888.90.33 143.5302.50.47 1.9
Cer[76]0.11.131.280.88 6.066.490.93 1.3
[76]0.22.262.610.87 13.1313.11.00 1.2
[77]0.33.43.21.06 21.0720.91.01 1.7
[76]0.44.533.661.24 26.2524.61.07 2.3
Pra[78]4050.0753.230.95129.86133.630.971.3
[79]4050.0736.91.37127.9893.691.371.4
[80]4050.0749.51.02118.621091.091.2
Ros[20]106.34.581.3834.0440.10.852.1
[81]106.35.81.0934.0445.90.742.6
[81]8050.4253.50.94272.753970.692.0
Flu[82]2091.03121.60.75 264.62199.51.33 2.7
[82]2091.03112.90.81 264.62192.51.37 2.1
[83]2091.03119.060.76 247.18144.621.71 1.9
Sim[84]403.8711.50.34 16.4826.50.62 4.2
[85]403.874.30.90 16.4829.30.56 7.7
[86]403.873.371.15 16.4832.500.51 4.1
Lov[87]404.915.000.98 20.2029.20.69 2.7
[82]202.462.21.12 10.4415.50.67 1.9
[88]404.914.061.21 22.4034.10.66 2.6
Rep[89]215.715.281.03 31.1934.770.90 1.4
[23]0.251.973.990.49 4.154.440.93 1.4
[10]0.53.925.220.75 7.966.361.25 1.5
Bos[90]1001689.9717860.95 11,365.7481801.39 1.9
[90]3005074.9450001.01 33,714.3018,4501.83 2.8
[90]60010,163.5099871.02 67,476.6341,4801.63 2.6
Table
Table 6. The predicted (Pred) pharmacokinetic parameters and the mean clinic observations (Obs) of atorvastatin (Ato), cerivastatin (Cer), pravastatin (Pra), rosuvastatin (Ros), fluvastatin (Flu), simvastatin (Sim), lovastatin (Lov), repaglinide (Rep) and bosentan (Bos) combined with cyclosporine (CsA).
Table 6. The predicted (Pred) pharmacokinetic parameters and the mean clinic observations (Obs) of atorvastatin (Ato), cerivastatin (Cer), pravastatin (Pra), rosuvastatin (Ros), fluvastatin (Flu), simvastatin (Sim), lovastatin (Lov), repaglinide (Rep) and bosentan (Bos) combined with cyclosporine (CsA).
Drug.DoseCsA DoseCmax-withCmax-withoutCmaxRAUC-withAUC-withoutAUCR
mgmgng/mLng/mL ng·h/mLng·h/mL
PredObsPredObsPredObsPredObsPredObsPredObs
Ato10 [14,16]886 BID 28-day28.537.33.313.58.61 10.66 123.2222617.21267.16 8.69
40 [15]17574.336214.626.55.10 13.66 204.78102649.84674.11 15.31
Cer0.2 [17]225 BID 1 day4.027.82.261.561.78 5.00 27.436.29.59.532.88 3.80
225 BID 7-day5.117.822.261.562.26 5.01 49.2745.39.59.535.19 4.75
Pra20 [18]420 BID 1 day91.28425.3523.273.60 3.61 211.3324962.7466.93.37 4.44
420 BID 28-day95.448025.3523.273.76 3.44 243.0224162.7466.93.87 4.30
10 [19]400 BID 1 day53.969612.6813.74.26 7.00 126.39307.0531.3726.254.03 11.70
400 BID 8-day55.829812.6813.74.40 7.15 141.86303.5231.3726.254.52 11.56
400 BID 29-day56.3115.112.6813.74.44 8.40 143.95345.931.3726.254.59 13.18
Ros10 [20]200 BID 1 day15.8639.86.34.582.52 8.69 62.561973440.11.84 4.91
200 BID 10-day18.4748.76.34.582.93 10.63 81.92843440.12.41 7.06
20 [20]200 BID 1 day31.766.5\\\\125.11308\\\\
200 BID 10-day36.9483.0\\\\163.81424\\\\
Flu40 [21]200 BID 1 day622.96869.4182.06211.93.42 4.10 1667.41948.8494.4549.43.37 3.55
200 BID 28-day646.871530190.14254.73.40 6.01 2062.232615.3640841.83.22 3.11
20 [92]200 BID 28-day314155\\\\1056373\\\\
Sim20 [93]266 BID 13.9120.61.99.97.32 2.08 28.841018.0839.63.57 2.55
10 [94]1h 200 BID a18.776.81.652 b11.38 3.40
2 h 200 BID a11.568.41.832 b6.32 4.20
3 h 200 BID a4.5612.11.442 b3.17 6.05
Lov 20 [18]420 BID 1 day30.58462.462.212.43 20.91 113.9824310.4415.510.92 15.68
420 BID 28-day33.12752.462.213.46 34.09 164.8945910.4415.515.79 29.61
10 [95]273 BID 10-day c 31.331753.89268.05 6.73
230 BID 10-day c 27.871103.89267.16 4.23
290 BID 10-day c 32.641103.89268.39 4.23
5 [96]300 d5.094.0514.65.09 0.88
5300 d5.9311.4714.65.93 2.49
10300 d12.2721.0814.612.27 4.58
15300 d18.9125.2514.618.91 5.49
20300 d24.7719.7914.624.77 4.30
20300 d25.719.0514.625.70 4.14
Rep0.25 [23]1004.646.71.973.92.36 1.72 10.3810.824.154.442.50 2.44
Bos e500 [24]300 BID18.557.928.474.742.19 1.67 189.4948.950.4424.783.76 1.97
a The concentration of simvastatin acid showed in Cmax rows are measured at point 1, 2 and 3 h. b The observed plasma concentrations of simvastatin acid without were lower than 2 ng/mL. c Different doses of cyclosporine were given to different patients. See detail in Section 2.3.7 “DDI of CsA with lovastatin”. d The concentrations in Cmax rows were at point 2 h following lovastatin 5 mg at weeks 0 and 3, 10 mg at week 6, 15 mg at week 9, and 20 mg at weeks 12 and 18. Control group was given lovastatin 10 mg. e The unit of Bosentan for Cmax and AUC were μg/mL and μg·h/mL, respectively.
Table
Table 7. Pharmacokinetics parameters of different factors on plasma concentrations of atorvastatin following 40 mg orally.
Table 7. Pharmacokinetics parameters of different factors on plasma concentrations of atorvastatin following 40 mg orally.
FactorsTmaxCmaxAUC
hng/mLng·h/mL
DDI Part
Whole1.17 75.99 264.82
Ki (2-fold)1.17 53.65185.75
Ki (0.5-fold)1.17 103.21378.30
CLint,up (2-fold)1.17 48.13 154.69
CLint,up (0.5-fold)1.33 107.23 419.24
CsA dosage (100 mg)1.17 51.48 178.84
CsA dosage (400 mg)1.33 106.77 392.31
Kt (2-fold)0.83 82.94 228.64
Kt (0.5-fold)1.67 60.31 266.78
SF (2-fold)1.17 48.13 154.69
SF (0.5-fold)1.33 107.23 419.24
non-CYP3A41.17 68.91 240.63
non-liver1.33 19.23 77.52
non-gut1.17 55.96 214.06
Alone Part
Whole1.5014.5966.88
non-liver CYP3A42.3325.23172.17
non-liver BCRP1.5014.9669.11
non-gut CYP3A41.5015.9174.17
non-gut P-gp1.5015.7669.35
non-gut-BCRP1.5015.4868.67
non-gut-MRP21.3315.6368.80
non-gut1.3320.3682.16
non-liver-OATP1.67136.131077.21

Share and Cite

MDPI and ACS Style

Yang, Y.; Li, P.; Zhang, Z.; Wang, Z.; Liu, L.; Liu, X. Prediction of Cyclosporin-Mediated Drug Interaction Using Physiologically Based Pharmacokinetic Model Characterizing Interplay of Drug Transporters and Enzymes. Int. J. Mol. Sci. 2020, 21, 7023. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms21197023

AMA Style

Yang Y, Li P, Zhang Z, Wang Z, Liu L, Liu X. Prediction of Cyclosporin-Mediated Drug Interaction Using Physiologically Based Pharmacokinetic Model Characterizing Interplay of Drug Transporters and Enzymes. International Journal of Molecular Sciences. 2020; 21(19):7023. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms21197023

Chicago/Turabian Style

Yang, Yiting, Ping Li, Zexin Zhang, Zhongjian Wang, Li Liu, and Xiaodong Liu. 2020. "Prediction of Cyclosporin-Mediated Drug Interaction Using Physiologically Based Pharmacokinetic Model Characterizing Interplay of Drug Transporters and Enzymes" International Journal of Molecular Sciences 21, no. 19: 7023. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms21197023

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