Identification and Evaluation of Clinical Substrates of Organic Anion Transporting Polypeptides 1B1 and 1B3

Rifampin- and Silymarin-Mediated Pharmacokinetic Interactions of Exogenous and Endogenous Substrates in a Transgenic OATP1B Mouse Model

Organic anion-transporting polypeptides (OATP) 1B1 and OATP1B3, encoded by the SLCO gene family of the solute carrier superfamily, are involved in the disposition of many exogenous and endogenous compounds. Preclinical rodent models help assess risks of pharmacokinetic interactions, but interspecies differences in transporter orthologs and expression limit direct clinical translation. An OATP1B transgenic mouse model comprising a rodent Slco1a/1b gene cluster knockout and human SLCO1B1 and SLCO1B3 gene insertions provides a potential physiologically relevant preclinical tool to predict pharmacokinetic interactions. Pharmacokinetics of exogenous probe substrates, pitavastatin and pravastatin, and endogenous OATP1B biomarkers, coproporphyrin-I and coproporphyrin-III, were determined in the presence and absence of known OATP/Oatp inhibitors, rifampin or silymarin (an extract of milk thistle [Silybum marianum]), in wild-type FVB mice and humanized OATP1B mice. Rifampin increased exposure of pitavastatin (4.6- and 2.8-fold), pravastatin (3.6- and 2.2-fold), and coproporphyrin-III (1.6- and 2.1-fold) in FVB and OATP1B mice, respectively, but increased coproporphyrin-I AUC0-24h only (1.8-fold) in the OATP1B mice. Silymarin did not significantly affect substrate AUC, likely because the silymarin flavonolignan concentrations were at or below their reported IC50 values for the relevant OATPs/Oatps. Silymarin increased the Cmax of pitavastatin 2.7-fold and pravastatin 1.9-fold in the OATP1B mice. The data of the OATP1B mice were similar to those of the pitavastatin and pravastatin clinical data; however, the FVB mice data more closely recapitulated pitavastatin clinical data than the data of the OATP1B mice, suggesting that the OATP1B mice are a reasonable, though costly, preclinical strain for predicting pharmacokinetic interactions when doses are optimized to achieve clinically relevant plasma concentrations.

Pharmacokinetic Effects of Different Models of Nonalcoholic Fatty Liver Disease in Transgenic Humanized OATP1B Mice

Organic anion transporting polypeptide (OATP) 1B1 and OATP1B3 (collectively, OATP1B) transporters encoded by the solute carrier organic anion transporter (SLCO) genes mediate uptake of multiple pharmaceutical compounds. Nonalcoholic steatohepatitis (NASH), a severe form of nonalcoholic fatty liver disease (NAFLD), decreases OATP1B abundance. This research characterized the pathologic and pharmacokinetics effects of three diet- and one chemical-induced NAFLD model in male and female humanized OATP1B mice, which comprises knock-out of rodent Oatp orthologs and insertion of human SLCO1B1 and SLCO1B3. Histopathology scoring demonstrated elevated steatosis and inflammation scores for all NAFLD-treatment groups. Female mice had minor changes in SLCO1B1 expression in two of the four NAFLD treatment groups, and pitavastatin (PIT) area under the concentration-time curve (AUC) increased in female mice in only one of the diet-induced models. OATP1B3 expression decreased in male and female mice in the chemical-induced NAFLD model, with a coinciding increase in PIT AUC, indicating the chemical-induced model may better replicate changes in OATP1B3 expression and OATP substrate disposition observed in NASH patients. This research also tested a reported multifactorial pharmacokinetic interaction between NAFLD and silymarin, an extract from milk thistle seeds with notable OATP-inhibitory effects. Males showed no change in PIT AUC, whereas female PIT AUC increased 1.55-fold from the diet alone and the 1.88-fold from the combination of diet with silymarin, suggesting that female mice are more sensitive to pharmacokinetic changes than male mice. Overall, the humanized OATP1B model should be used with caution for modeling NAFLD and multifactorial pharmacokinetic interactions. SIGNIFICANCE STATEMENT: Advanced stages of NAFLD cause decreased hepatic OATP1B abundance and increase systemic exposure to OATP substrates in human patients. The humanized OATP1B mouse strain may provide a clinically relevant model to recapitulate these observations and predict pharmacokinetic interactions in NAFLD. This research characterized three diet-induced and one drug-induced NAFLD model in a humanized OATP1B mouse model. Additionally, a multifactorial pharmacokinetic interaction was observed between silymarin and NAFLD.

Toward improved predictions of pharmacokinetics of transported drugs in hepatic impairment: Insights from the extended clearance model

Hepatic impairment (HI) moderately (<5-fold) affects the systemic exposure (i.e., area under the plasma concentration-time curve [AUC]) of drugs that are substrates of the hepatic sinusoidal organic anion transporting polypeptide (OATP) transporters and are excreted unchanged in the bile and/or urine. However, the effect of HI on their AUC is much greater (>10-fold) for drugs that are also substrates of cytochrome P450 (CYP) 3A enzymes. Using the extended clearance model, through simulations, we identified the ratio of sinusoidal efflux clearance (CL) over the sum of metabolic and biliary CLs as important in predicting the impact of HI on the AUC of dual OATP/CYP3A substrates. Because HI may reduce hepatic CYP3A-mediated CL to a greater extent than biliary efflux CL, the greater the contribution of the former versus the latter, the greater the impact of HI on drug AUC ratio (AUCRHI ). Using physiologically-based pharmacokinetic modeling and simulation, we predicted relatively well the AUCRHI of OATP substrates that are not significantly metabolized (pitavastatin, rosuvastatin, valsartan, and gadoxetic acid). However, there was a trend toward underprediction of the AUCRHI of the dual OATP/CYP3A4 substrates fimasartan and atorvastatin. These predictions improved when the sinusoidal efflux CL of these two drugs was increased in healthy volunteers (i.e., before incorporating the effect of HI), and by modifying the directionality of its modulation by HI (i.e., increase or decrease). To accurately predict the effect of HI on AUC of hepatobiliary cleared drugs it is important to accurately predict all hepatobiliary pathways, including sinusoidal efflux CL.

CNS Drug Delivery in Stroke: Improving Therapeutic Translation From the Bench to the Bedside

Drug development for ischemic stroke is challenging as evidenced by the paucity of therapeutics that have advanced beyond a phase III trial. There are many reasons for this lack of clinical translation including factors related to the experimental design of preclinical studies. Often overlooked in therapeutic development for ischemic stroke is the requirement of effective drug delivery to the brain, which is critical for neuroprotective efficacy of several small and large molecule drugs. Advancing central nervous system drug delivery technologies implies a need for detailed comprehension of the blood-brain barrier (BBB) and neurovascular unit. Such knowledge will permit the innate biology of the BBB/neurovascular unit to be leveraged for improved bench-to-bedside translation of novel stroke therapeutics. In this review, we will highlight key aspects of BBB/neurovascular unit pathophysiology and describe state-of-the-art approaches for optimization of central nervous system drug delivery (ie, passive diffusion, mechanical opening of the BBB, liposomes/nanoparticles, transcytosis, intranasal drug administration). Additionally, we will discuss how endogenous BBB transporters represent the next frontier of drug delivery strategies for stroke. Overall, this review will provide cutting edge perspective on how central nervous system drug delivery must be considered for the advancement of new stroke drugs toward human trials.

Physiologically based pharmacokinetic (PBPK) modeling to predict the pharmacokinetics of irbesartan in different CYP2C9 genotypes

Irbesartan, a potent and selective angiotensin II type-1 (AT1) receptor blocker (ARB), is one of the representative medications for the treatment of hypertension. Cytochrome P450 (CYP) 2C9 is primarily involved in the oxidation of irbesartan. CYP2C9 is highly polymorphic, and genetic polymorphism of this enzyme is the leading cause of significant alterations in the pharmacokinetics of irbesartan. This study aimed to establish the physiologically based pharmacokinetic (PBPK) model to predict the pharmacokinetics of irbesartan in different CYP2C9 genotypes. The irbesartan PBPK model was established using the PK-Sim® software. Our previously reported pharmacogenomic data for irbesartan was leveraged in the development of the PBPK model and collected clinical pharmacokinetic data for irbesartan was used for the validation of the model. Physicochemical and ADME properties of irbesartan were obtained from previously reported data, predicted by the modeling software, or optimized to fit the observed plasma concentration-time profiles. Model evaluation was performed by comparing the predicted plasma concentration-time profiles and pharmacokinetic parameters to the observed results. Predicted plasma concentration-time profiles were visually similar to observed profiles. Predicted AUCinf in CYP2C9*1/*3 and CYP2C9*1/*13 genotypes were increased by 1.54- and 1.62-fold compared to CYP2C9*1/*1 genotype, respectively. All fold error values for AUC and Cmax in non-genotyped and CYP2C9 genotyped models were within the two-fold error criterion. We properly established the PBPK model of irbesartan in different CYP2C9 genotypes. It can be used to predict the pharmacokinetics of irbesartan for personalized pharmacotherapy in individuals of various races, ages, and CYP2C9 genotypes.

Serum metabolomic analyses reveal the potential metabolic biomarkers for prediction of amatoxin poisoning

Amatoxin poisoning leads to over 90% of deaths in mushroom poisoning. The objective of present study was to identify the potential metabolic biomarkers for early diagnosis of amatoxin poisoning. Serum samples were collected from 61 patients with amatoxin poisoning and 61 healthy controls. An untargeted metabolomics analysis was performed using the ultra-high-performance liquid chromatography-quadrupole time-of-flight tandem mass spectrometry (UPLC-QTOF-MS/MS). Multivariate statistical analysis revealed that the patients with amatoxin poisoning could be clearly separated from healthy controls on the basis of their metabolic fingerprints. There were 33 differential metabolites including 15 metabolites up-regulated metabolites and 18 down-regulated metabolites in patients with amatoxin poisoning compared to healthy controls. These metabolites mainly enriched in the lipid metabolism and amino acid metabolism pathways, such as Glycerophospholipid metabolism, Sphingolipid metabolism, Phenylalanine tyrosine and typtophan biosynthesis, Tyrosine metabolism, Arginine and proline metabolism, which may serve important roles in the amatoxin poisoning. Among the differential metabolites, a total of 8 significant metabolic markers were identified for discriminating patients with amatoxin poisoning from healthy controls, including Glycochenodeoxycholate-3-sulfate (GCDCA-S), 11-Oxo-androsterone glucuronide, Neomenthol-glucuronide, Dehydroisoandrosterone 3-glucuronide, Glucose 6-phosphate (G6P), Lanthionine ketimine, Glycerophosphocholine (GPC) and Nicotinamide ribotide, which achieved satisfactory diagnostic accuracy (AUC>0.8) in both discovery and validation cohorts. Strikingly, the Pearson's correlation analysis indicated that 11-Oxo-androsterone glucuronide, G6P and GCDCA-S were positively correlated with the liver injury induced by amatoxin poisoning. The findings of the current study may provide insight into the pathological mechanism of amatoxin poisoning and screened out the reliable metabolic biomarkers to contribute the clinical early diagnosis of amatoxin poisoning.

Neither Gastric Bypass Surgery Nor Diet-Induced Weight-Loss Affect OATP1B1 Activity as Measured by Rosuvastatin Oral Clearance

INTRODUCTION

Rosuvastatin pharmacokinetics is mainly dependent on the activity of hepatic uptake transporter OATP1B1. In this study, we aimed to investigate and disentangle the effect of Roux-en-Y gastric bypass (RYGB) and weight loss on oral clearance (CL/F) of rosuvastatin as a measure of OATP1B1-activity.

METHODS

Patients with severe obesity preparing for RYGB (n = 40) or diet-induced weight loss (n = 40) were included and followed for 2 years, with four 24-hour pharmacokinetic investigations. Both groups underwent a 3-week low-energy diet (LED; < 1200 kcal/day), followed by RYGB or a 6-week very-low-energy diet (VLED; < 800 kcal/day).

RESULTS

A total of 80 patients were included in the RYGB group (40 patients) and diet-group (40 patients). The weight loss was similar between the groups following LED and RYGB. The LED induced a similar (mean [95% CI]) decrease in CL/F in both intervention groups (RYGB: 16% [0, 31], diet: 23% [8, 38]), but neither induced VLED resulted in any further changes in CL/F. At Year 2, CL/F had increased by 21% from baseline in the RYGB group, while it was unaltered in the diet group. Patients expressing the reduced function SLCO1B1 variants (c.521TC/CC) showed similar changes in CL/F over time compared with patients expressing the wild-type variant.

CONCLUSIONS

Neither body weight, weight loss nor RYGB per se seem to affect OATP1B1 activity to a clinically relevant degree. Overall, the observed changes in rosuvastatin pharmacokinetics were minor, and unlikely to be of clinical relevance.

Comparative Modelling of Organic Anion Transporting Polypeptides: Structural Insights and Comparison of Binding Modes

To better understand the functionality of organic anion transporting polypeptides (OATPs) and to design new ligands, reliable structural data of each OATP is needed. In this work, we used a combination of homology model with molecular dynamics simulations to generate a comprehensive structural dataset, that encompasses a diverse set of OATPs but also their relevant conformations. Our OATP models share a conserved transmembrane helix folding harbouring a druggable binding pocket in the shape of an inner pore. Our simulations suggest that the conserved salt bridges at the extracellular region between residues on TM1 and TM7 might influence the entrance of substrates. Interactions between residues on TM1 and TM4 within OATP1 family shown their importance in transport of substrates. Additionally, in transmembrane (TM) 1/2, a known conserved element, interact with two identified motifs in the TM7 and TM11. Our simulations suggest that TM1/2-TM7 interaction influence the inner pocket accessibility, while TM1/2-TM11 salt bridges control the substrate binding stability.

Toxicological Implications of Biological Heterogeneity

From a micro to macro scale of biological organization, macromolecular diversity and biological heterogeneity are fundamental properties of biological systems. Heterogeneity may result from genetic, epigenetic, and non-genetic characteristics (e.g., tissue microenvironment). Macromolecular diversity and biological heterogeneity are tolerated as long as the sustenance and propagation of life are not disrupted. They also provide the raw materials for microevolutionary changes that may help organisms adapt to new selection pressures arising from the environment. Sequence evolution, functional divergence, and positive selection of gene and promoter dosage play a major role in the evolution of life's diversity including complex metabolic networks, which is ultimately reflected in changes in the allele frequency over time. Robustness in evolvable biological systems is conferred by functional redundancy that is often created by macromolecular diversity and biological heterogeneity. The ability to investigate biological macromolecules at an increasingly finer level has uncovered a wealth of information in this regard. Therefore, the dynamics of biological complexity should be taken into consideration in biomedical research.

Physiologically Based Pharmacokinetic Modeling of Rosuvastatin to Predict Transporter-Mediated Drug-Drug Interactions

Purpose

To build a physiologically based pharmacokinetic (PBPK) model of the clinical OATP1B1/OATP1B3/BCRP victim drug rosuvastatin for the investigation and prediction of its transporter-mediated drug-drug interactions (DDIs).

Methods

The Rosuvastatin model was developed using the open-source PBPK software PK-Sim®, following a middle-out approach. 42 clinical studies (dosing range 0.002–80.0 mg), providing rosuvastatin plasma, urine and feces data, positron emission tomography (PET) measurements of tissue concentrations and 7 different rosuvastatin DDI studies with rifampicin, gemfibrozil and probenecid as the perpetrator drugs, were included to build and qualify the model.

Results

The carefully developed and thoroughly evaluated model adequately describes the analyzed clinical data, including blood, liver, feces and urine measurements. The processes implemented to describe the rosuvastatin pharmacokinetics and DDIs are active uptake by OATP2B1, OATP1B1/OATP1B3 and OAT3, active efflux by BCRP and Pgp, metabolism by CYP2C9 and passive glomerular filtration. The available clinical rifampicin, gemfibrozil and probenecid DDI studies were modeled using in vitro inhibition constants without adjustments. The good prediction of DDIs was demonstrated by simulated rosuvastatin plasma profiles, DDI AUC_last ratios (AUC_last during DDI/AUC_last without co-administration) and DDI C_max ratios (C_max during DDI/C_max without co-administration), with all simulated DDI ratios within 1.6-fold of the observed values.

Conclusions

A whole-body PBPK model of rosuvastatin was built and qualified for the prediction of rosuvastatin pharmacokinetics and transporter-mediated DDIs. The model is freely available in the Open Systems Pharmacology model repository, to support future investigations of rosuvastatin pharmacokinetics, rosuvastatin therapy and DDI studies during model-informed drug discovery and development (MID3).

Supplementary Information

The online version contains supplementary material available at 10.1007/s11095-021-03109-6.

Evaluation of the Effects of Repeat-Dose Dabrafenib on the Single-Dose Pharmacokinetics of Rosuvastatin (OATP1B1/1B3 Substrate) and Midazolam (CYP3A4 Substrate)

Dabrafenib is an oral BRAF kinase inhibitor approved for the treatment of various BRAF V600 mutation–positive solid tumors. In vitro observations suggesting cytochrome P450 (CYP) 3A induction and organic anion transporting polypeptide (OATP) inhibition prompted us to evaluate the effect of dabrafenib 150 mg twice daily on the pharmacokinetics of midazolam 3 mg (CYP3A substrate) and rosuvastatin 10 mg (OATP1B1/1B3 substrate) in a clinical phase 1, open‐label, fixed‐sequence study in patients with BRAF V600 mutation–positive tumors. Repeat dabrafenib dosing resulted in a 2.56‐fold increase in rosuvastatin maximum observed concentration (C_max), an earlier time to C_max, but only a 7% increase in area under the concentration‐time curve from time 0 (predose) extrapolated to infinite time. Midazolam C_max and AUC extrapolated to infinite time decreased by 47% and 65%, respectively, with little effect on time to C_max. No new safety findings were reported. Exposure of drugs that are CYP3A4 substrates is likely to decrease when coadministered with dabrafenib. Concentrations of medicinal products that are sensitive OATP1B1/1B3 substrates may increase during the absorption phase.

The Emerging Role of the SLCO1B3 Protein in Cancer Resistance

Currently, chemotherapy is one of the mainstays of oncologic therapies. But the efficacy of chemotherapy is often limited by drug resistance and severe side effects. Consequently, it is becoming increasingly important to investigate the underlying mechanism and overcome the problem of anticancer chemotherapy resistance. The solute carrier organic anion transporter family member 1B3 (SLCO1B3), a functional transporter normally expressed in the liver, transports a variety of endogenous and exogenous compounds, including hormones and their conjugates as well as some anticancer drugs. The extrahepatic expression of SLCO1B3 has been detected in different cancer cell lines and cancer tissues. Recently, accumulating data indicates that the abnormal expression and function of SLCO1B3 are involved in resistance to anticancer drugs, such as taxanes, camptothecin and its analogs, SN-38, and Androgen Deprivation Therapy (ADT) in breast, prostate, lung, hepatic, and colorectal cancer, respectively. Thus, more investigations have been implemented to identify the potential SLCO1B3-related mechanisms of cancer drug resistance. In this review, we focus on the emerging roles of SLCO1B3 protein in the development of cancer chemotherapy resistance and briefly discuss the mechanisms of resistance. Elucidating the function of SLCO1B3 in chemoresistance may bring out novel therapeutic strategies for cancer treatment.

Alteration in the Plasma Concentrations of Endogenous Organic Anion-Transporting Polypeptide 1B Biomarkers in Patients with Non-Small Cell Lung Cancer Treated with Paclitaxel

Paclitaxel has been considered to cause OATP1B-mediated drug-drug interactions at therapeutic doses; however, its clinical relevance has not been demonstrated. This study aimed to elucidate in vivo inhibition potency of paclitaxel against OATP1B1 and OATP1B3 using endogenous OATP1B biomarkers. Paclitaxel is an inhibitor of OATP1B1 and OATP1B3, with K_i of 0.579 ± 0.107 and 5.29 ± 3.87 μM, respectively. Preincubation potentiated its inhibitory effect on both OATP1B1 and OATP1B3, with K_i of 0.154 ± 0.031 and 0.624 ± 0.183 μM, respectively. Ten patients with non-small cell lung cancer who received 200 mg/m² of paclitaxel by a 3-hour infusion were recruited. Plasma concentrations of 10 endogenous OATP1B biomarkers-namely, coproporphyrin I, coproporphyrin III, glycochenodeoxycholate-3-sulfate, glycochenodeoxycholate-3-glucuronide, glycodeoxycholate-3-sulfate, glycodeoxycholate-3-glucuronide, lithocholate-3-sulfate, glycolithocholate-3-sulfate, taurolithocholate-3-sulfate, and chenodeoxycholate-24-glucuronide-were determined in the patients with non-small cell lung cancer on the day before paclitaxel administration and after the end of paclitaxel infusion for 7 hours. Paclitaxel increased the area under the plasma concentration-time curve (AUC) of the endogenous biomarkers 2- to 4-fold, although a few patients did not show any increment in the AUC ratios of lithocholate-3-sulfate, glycolithocholate-3-sulfate, and taurolithocholate-3-sulfate. Therapeutic doses of paclitaxel for the treatment of non-small cell lung cancer (200 mg/m²) will cause significant OATP1B1 inhibition during and at the end of the infusion. This is the first demonstration that endogenous OATP1B biomarkers could serve as surrogate biomarkers in patients. SIGNIFICANCE STATEMENT: Endogenous biomarkers can address practical and ethical issues in elucidating transporter-mediated drug-drug interaction (DDI) risks of anticancer drugs clinically. We could elucidate a significant increment of the plasma concentrations of endogenous OATP1B biomarkers after a 3-hour infusion (200 mg/m²) of paclitaxel, a time-dependent inhibitor of OATP1B, in patients with non-small cell lung cancer. The endogenous OATP1B biomarkers are useful to assess the possibility of OATP1B-mediated DDIs in patients and help in appropriately designing a dosing schedule to avoid the DDIs.