Characterization of Bile Acid Sulfate Conjugates as Substrates of Human Organic Anion Transporting Polypeptides

Drug interactions involving the inhibition of hepatic organic anion transporting polypeptides (OATPs) 1B1 and OATP1B3 are considered important. Therefore, we sought to study various sulfated bile acids (BA-S) as potential clinical OATP1B1/3 biomarkers. It was determined that BA-S [e.g., glycochenodeoxycholic acid 3-O-sulfate (GCDCA-S) and glycodeoxycholic acid 3-O-sulfate (GDCA-S)] are substrates of OATP1B1, OATP1B3, and sodium-dependent taurocholic acid cotransporting polypeptide (NTCP) transfected into human embryonic kidney 293 cells, with minimal uptake evident for other solute carriers (SLCs) like OATP2B1, organic anion transporter 2, and organic cation transporter 1. It was also shown that BA-S uptake by plated human hepatocytes (PHH) was inhibited (≥96%) by a pan-SLC inhibitor (rifamycin SV), and there was greater inhibition (≥77% versus ≤12%) with rifampicin (OATP1B1/3-selective inhibitor) than a hepatitis B virus myristoylated-preS1 peptide (NTCP-selective inhibitor). Estrone 3-sulfate was also used as an OATP1B1-selective inhibitor. In this instance, greater inhibition was observed with GDCA-S (76%) than GCDCA-S (52%). The study was expanded to encompass the measurement of GCDCA-S and GDCA-S in plasma of SLCO1B1 genotyped subjects. The geometric mean GDCA-S concentration was 2.6-fold (90% confidence interval 1.6, 4.3; P = 2.1 × 10^-4) and 1.3-fold (1.1, 1.7; P = 0.001) higher in individuals homozygous and heterozygous for the SLCO1B1 c.521T > C loss-of-function allele, respectively. For GCDCA-S, no significant difference was noted [1.2-fold (0.8, 1.7; P = 0.384) and 0.9-fold (0.8, 1.1; P = 0.190), respectively]. This supported the in vitro data indicating that GDCA-S is a more OATP1B1-selective substrate (versus GCDCA-S). It is concluded that GCDCA-S and GDCA-S are viable plasma-based OATP1B1/3 biomarkers, but they are both less OATP1B1-selective when compared to their corresponding 3-O-glucuronides (GCDCA-3G and GDCA-3G). Additional studies are needed to determine their utility versus more established biomarkers, such as coproporphyrin I, for assessing inhibitors with different OATP1B1 (versus OATP1B3) inhibition signatures.

Effect of Hepatic Impairment on OATP1B Activity: Quantitative Pharmacokinetic Analysis of Endogenous Biomarker and Substrate Drugs

Hepatic impairment (HI) is known to modulate drug disposition and may lead to elevated plasma exposure. The aim of this study was to quantitate the in vivo OATP1B-mediated hepatic uptake activity in populations with varying degrees of HI. First, we measured baseline levels of plasma coproporphyrin-I, an endogenous OATP1B biomarker, in an open-label, parallel cohort study in adult subjects with normal liver function and mild, moderate, and severe HI (n = 24, 6/cohort). The geometric mean plasma concentrations of coproporphyrin-I were 1.66-fold, 2.81-fold (P < 0.05), and 7.78-fold (P < 0.0001) higher in mild, moderate, and severe impairment than those healthy controls. Second, we developed a dataset of 21 OATP1B substrate drugs with HI data extracted from literature. Median disease-to-healthy plasma area under the curve (AUC) ratios for substrate drugs were ~ 1.4, 3.0, and 6.4 for mild, moderate, and severe HI, respectively. Additionally, significant linear relationship was noted between AUC ratios of substrate drugs without and with co-administration of rifampin, a prototypic OATP1B inhibitor, and AUC ratios in moderate (P < 0.01) and severe (P < 0.001) HI. Third, a physiologically-based pharmacokinetic model analysis was conducted with 10 substrate drugs following estimation of relative OATP1B functional activity in virtual disease population models using coproporphyrin-I data and verification of substrate models with rifampin drug-drug interaction data. This approach adequately predicted plasma AUC change particularly in moderate (9 of 10 within 2-fold) and severe (5 of 5 within 2-fold) HI. Collective findings indicate progressive reduction, by as much as 90-92%, in OATP1B activity in the HI population.

Generation of a VeroE6 Pgp gene knock out cell line and its use in SARS-CoV-2 antiviral study

Vero cells are widely used for antiviral tests and virology research for SARS-CoV-2 as well as viruses from various other families. However, Vero cells generally express high levels of multi-drug resistance 1 (MDR1) or Pgp protein, the efflux transporter of foreign substances including many antiviral compounds, affecting the antiviral activity as well as interpretation of data. To address this, a Pgp gene knockout VeroE6 cell line (VeroE6-Pgp-KO) was generated using CRISPR-CAS9 technology. These cells no longer expressed the Pgp protein as indicated by flow cytometry analysis following staining with a Pgp-specific monoclonal antibody. They also showed significantly reduced efflux transporter activity in the calcein acetoxymethyl ester (calcein AM) assay. The VeroE6-Pgp-KO cells and the parental VeroE6 cells were each infected with SARS-CoV-2 to test antiviral activities of remdesivir and nirmatrelvir, two known Pgp substrates, in the presence or absence of a Pgp inhibitor. The compounds showed antiviral activities in VeroE6-Pgp-KO cells similar to that observed in the presence of the Pgp inhibitor. Thus, the newly established VeroE6-Pgp-KO cell line adds a new in vitro virus infection system for SARS-CoV-2 and possibly other viruses to test antiviral therapies without a need to control the Pgp activity. Removal of the Pgp inhibitor for antiviral assays will lead to less data variation and prevent failed assays.

Biomarker-Informed Model-Based Risk Assessment of Organic Anion Transporting Polypeptide 1B Mediated Drug-Drug Interactions

Quantitative prediction of drug-drug interactions (DDIs) involving organic anion transporting polypeptide (OATP)1B1/1B3 inhibition is limited by uncertainty in the translatability of experimentally determined in vitro inhibition potency (half-maximal inhibitory concentration (IC₅₀ )). This study used an OATP1B endogenous biomarker-informed physiologically-based pharmacokinetic (PBPK) modeling approach to predict the effect of inhibitor drugs on the pharmacokinetics (PKs) of OATP1B substrates. Initial static analysis with about 42 inhibitor drugs, using in vitro IC₅₀ values and unbound liver inlet concentrations (I_in,max,u ), suggested in vivo OATP1B inhibition risk for drugs with R-value (1+ I_in,max,u /IC₅₀ ) above 1.5. A full-PBPK model accounting for transporter-mediated hepatic disposition was developed for coproporphyrin I (CP-I), an endogenous OATP1B biomarker. For several inhibitors (cyclosporine, diltiazem, fenebrutinib, GDC-0810, itraconazole, probenecid, and rifampicin at 3 different doses), PBPK models were developed and verified against available CP-I plasma exposure data to obtain in vivo OATP1B inhibition potency-which tend to be lower than the experimentally measured in vitro IC₅₀ by about 2-fold (probenecid and rifampicin) to 37-fold (GDC-0810). Models verified with CP-I data are subsequently used to predict DDIs with OATP1B probe drugs, rosuvastatin and pitavastatin. The predicted and observed area under the plasma concentration-time curve ratios are within 20% error in 55% cases, and within 30% error in 89% cases. Collectively, this comprehensive study illustrates the adequacy and utility of endogenous biomarker-informed PBPK modeling in mechanistic understanding and quantitative predictions of OATP1B-mediated DDIs in drug development.

Quantitative prediction of breast cancer resistant protein mediated drug-drug interactions using physiologically-based pharmacokinetic modeling

Quantitative assessment of drug‐drug interactions (DDIs) involving breast cancer resistance protein (BCRP) inhibition is challenged by overlapping substrate/inhibitor specificity. This study used physiologically‐based pharmacokinetic (PBPK) modeling to delineate the effects of inhibitor drugs on BCRP‐ and organic anion transporting polypeptide (OATP)1B‐mediated disposition of rosuvastatin, which is a recommended BCRP clinical probe. Initial static model analysis using in vitro inhibition data suggested BCRP/OATP1B DDI risk while considering regulatory cutoff criteria for a majority of inhibitors assessed (25 of 27), which increased rosuvastatin plasma exposure to varying degree (~ 0–600%). However, rosuvastatin area under plasma concentration‐time curve (AUC) was minimally impacted by BCRP inhibitors with calculated G‐value (= gut concentration/inhibition potency) below 100. A comprehensive PBPK model accounting for intestinal (OATP2B1 and BCRP), hepatic (OATP1B, BCRP, and MRP4), and renal (OAT3) transport mechanisms was developed for rosuvastatin. Adopting in vitro inhibition data, rosuvastatin plasma AUC changes were predicted within 25% error for 9 of 12 inhibitors evaluated via PBPK modeling. This study illustrates the adequacy and utility of a mechanistic model‐informed approach in quantitatively assessing BCRP‐mediated DDIs.

Identification of Glycochenodeoxycholate 3-O-Glucuronide and Glycodeoxycholate 3-O-Glucuronide as Highly Sensitive and Specific OATP1B1 Biomarkers

The aim of this study was to investigate the sensitivity and specificity of endogenous glycochenodeoxycholate and glycodeoxycholate 3-O-glucuronides (GCDCA-3G and GDCA-3G) as substrates for organic anion transporting polypeptide 1B1 (OATP1B1) in humans. We measured fasting levels of plasma GCDCA-3G and GDCA-3G using liquid chromatography-tandem mass spectrometry in 356 healthy volunteers. The mean plasma levels of both compounds were ~ 50% lower in women than in men (P = 2.25 × 10^-18 and P = 4.73 × 10^-9 ). In a microarray-based genome-wide association study, the SLCO1B1 rs4149056 (c.521T>C, p.Val174Ala) variation showed the strongest association with the plasma GCDCA-3G (P = 3.09 × 10^-30 ) and GDCA-3G (P = 1.60 × 10^-17 ) concentrations. The mean plasma concentration of GCDCA-3G was 9.2-fold (P = 8.77 × 10^-31 ) and that of GDCA-3G was 6.4-fold (P = 2.45x10^-13 ) higher in individuals with the SLCO1B1 c.521C/C genotype than in those with the c.521T/T genotype. No other variants showed independent genome-wide significant associations with GCDCA-3G or GDCA-3G. GCDCA-3G was highly efficacious in detecting the SLCO1B1 c.521C/C genotype with an area under the receiver operating characteristic curve of 0.996 (P < 0.0001). The sensitivity (98-99%) and specificity (100%) peaked at a cutoff value of 180 ng/mL for men and 90 ng/mL for women. In a haplotype-based analysis, SLCO1B1*5 and *15 were associated with reduced, and SLCO1B1*1B, *14, and *35 with increased OATP1B1 function. In vitro, both GCDCA-3G and GDCA-3G showed at least 6 times higher uptake by OATP1B1 than OATP1B3 or OATP2B1. These data indicate that the hepatic uptake of GCDCA-3G and GDCA-3G is predominantly mediated by OATP1B1. GCDCA-3G, in particular, is a highly sensitive and specific OATP1B1 biomarker in humans.

Physiologically-Based Pharmacokinetic Modeling Approach to Predict Rifampin-Mediated Intestinal P-Glycoprotein Induction

Physiologically‐based pharmacokinetic (PBPK) modeling is a powerful tool to quantitatively describe drug disposition profiles in vivo, thereby providing an alternative to predict drug–drug interactions (DDIs) that have not been tested clinically. This study aimed to predict effects of rifampin‐mediated intestinal P‐glycoprotein (Pgp) induction on pharmacokinetics of Pgp substrates via PBPK modeling. First, we selected four Pgp substrates (digoxin, talinolol, quinidine, and dabigatran etexilate) to derive in vitro to in vivo scaling factors for intestinal Pgp kinetics. Assuming unbound Michaelis‐Menten constant (K_m) to be intrinsic, we focused on the scaling factors for maximal efflux rate (J_max) to adequately recover clinically observed results. Next, we predicted rifampin‐mediated fold increases in intestinal Pgp abundances to reasonably recover clinically observed DDI results. The modeling results suggested that threefold to fourfold increases in intestinal Pgp abundances could sufficiently reproduce the DDI results of these Pgp substrates with rifampin. Hence, the obtained fold increases can potentially be applicable to DDI prediction with other Pgp substrates.

Quantitative Contribution of Six Major Transporters to the Hepatic Uptake of Drugs: "SLC-Phenotyping" Using Primary Human Hepatocytes

Hepatic uptake transporters [solute carriers (SLCs)], including organic anion transporting polypeptide (OATP) 1B1, OATP1B3, OATP2B1, sodium-dependent taurocholate cotransporting polypeptide (NTCP), and organic anion (OAT2) and organic cation (OCT1) transporters, play a key role in determining the systemic and liver exposure of chemically diverse drugs. Here, we established a phenotyping approach to quantify the contribution of the six SLCs, and passive diffusion, to the overall uptake using plated human hepatocytes (PHHs). First, selective inhibitor conditions were identified by screening about 20 inhibitors across the six SLCs using single-transfected human embryonic kidney 293 cells. Data implied rifamycin SV (20 µM) inhibits three OATPs, while rifampicin (5 µM) inhibits OATP1B1/1B3 only. Further, hepatitis B virus myristoylated-preS1 peptide (0.1 µM), quinidine (100 µM), and ketoprofen (100-300 µM) are relatively selective against NTCP, OCT1, and OAT2, respectively. Second, using these inhibitory conditions, the fraction transported (f_t ) by the individual SLCs was characterized for 20 substrates with PHH. Generally, extended clearance classification system class 1A/3A (e.g., warfarin) and 1B/3B compounds (e.g., statins) showed predominant OAT2 and OATP1B1/1B3 contribution, respectively. OCT1-mediated uptake was prominent for class 2/4 compounds (e.g., metformin). Third, in vitro f_t values were corrected using quantitative proteomics data to obtain "scaled f_t " Fourth, in vitro-in vivo extrapolation of the scaled OATP1B1/1B3 f_t was assessed, leveraging statin clinical drug-drug interaction data with rifampicin as the perpetrator. Finally, we outlined a novel stepwise strategy to implement phenotypic characterization of SLC-mediated hepatic uptake for new molecular entities and drugs in a drug discovery and development setting.

Simultaneous Assessment of Transporter-Mediated Drug-Drug Interactions Using a Probe Drug Cocktail in Cynomolgus Monkey

We aim to establish an in vivo preclinical model to enable simultaneous assessment of inhibition potential of an investigational drug on clinically relevant drug transporters, organic anion-transporting polypeptide (OATP)1B, breast cancer resistance protein (BCRP), P-glycoprotein (P-gp), and organic anion transporter (OAT)3. Pharmacokinetics of substrate cocktail consisting of pitavastatin (OATP1B substrate), rosuvastatin (OATP1B/BCRP/OAT3), sulfasalazine (BCRP), and talinolol (P-gp) were obtained in cynomolgus monkey-alone or in combination with transporter inhibitors. Single-dose rifampicin (30 mg/kg) significantly (P < 0.01) increased the plasma exposure of all four drugs, with a marked effect on pitavastatin and rosuvastatin [area under the plasma concentration-time curve (AUC) ratio ∼21-39]. Elacridar, BCRP/P-gp inhibitor, increased the AUC of sulfasalazine, talinolol, as well as rosuvastatin and pitavastatin. An OAT1/3 inhibitor (probenecid) significantly (P < 0.05) impacted the renal clearance of rosuvastatin (∼8-fold). In vitro, rifampicin (10 µM) inhibited uptake of pitavastatin, rosuvastatin, and sulfasalazine by monkey and human primary hepatocytes. Transport studies using membrane vesicles suggested that all probe substrates, except talinolol, are transported by cynoBCRP, whereas talinolol is a cynoP-gp substrate. Elacridar and rifampicin inhibited both cynoBCRP and cynoP-gp in vitro, indicating potential for in vivo intestinal efflux inhibition. In conclusion, a probe substrate cocktail was validated to simultaneously evaluate perpetrator impact on multiple clinically relevant transporters using the cynomolgus monkey. The results support the use of the cynomolgus monkey as a model that could enable drug-drug interaction risk assessment, before advancing a new molecular entity into clinical development, as well as providing mechanistic insights on transporter-mediated interactions.

Application of Physiologically Based Pharmacokinetic Modeling in Understanding Bosutinib Drug-Drug Interactions: Importance of Intestinal P-Glycoprotein

Bosutinib is an orally available Src/Abl tyrosine kinase inhibitor indicated for the treatment of patients with Ph+ chronic myelogenous leukemia at a clinically recommended dose of 500 mg once daily. Clinical results indicated that increases in bosutinib oral exposures were supraproportional at the lower doses (50-200 mg) and approximately dose-proportional at the higher doses (200-600 mg). Bosutinib is a substrate of CYP3A4 and P-glycoprotein and exhibits pH-dependent solubility with moderate intestinal permeability. These findings led us to investigate the factors influencing the underlying pharmacokinetic mechanisms of bosutinib with physiologically based pharmacokinetic (PBPK) models. Our primary objectives were to: 1) refine the previously developed bosutinib PBPK model on the basis of the latest oral bioavailability data and 2) verify the refined PBPK model with P-glycoprotein kinetics on the basis of the bosutinib drug-drug interaction (DDI) results with ketoconazole and rifampin. Additionally, the verified PBPK model was applied to predict bosutinib DDIs with dual CYP3A/P-glycoprotein inhibitors. The results indicated that 1) the refined PBPK model adequately described the observed plasma concentration-time profiles of bosutinib and 2) the verified PBPK model reasonably predicted the effects of ketoconazole and rifampin on bosutinib exposures by accounting for intestinal P-glycoprotein inhibition/induction. These results suggested that bosutinib DDI mechanism could involve not only CYP3A4-mediated metabolism but also P-glycoprotein-mediated efflux on absorption. In summary, P-glycoprotein kinetics could constitute an element in the PBPK models critical to understanding the pharmacokinetic mechanism of dual CYP3A/P-glycoprotein substrates, such as bosutinib, that exhibit nonlinear pharmacokinetics owing largely to a saturation of intestinal P-glycoprotein-mediated efflux.

Evaluation of an automated in vitro micronucleus assay in CHO-K1 cells

In this paper, we describe the evaluation of an automated in vitro micronucleus assay using CHO-K1 cells in 96-well plates. CHO-K1 cells were pre-loaded with a cell dye that stains the cytoplasm, after which the cells were treated with the test compounds for either 3h (for the +S9 condition) or 24h (for the -S9 condition). A total of 10 concentrations were tested, of which the top five concentrations were scored (limited by either cytotoxicity or solubility). At the end of the incubation period the cells were fixed and their DNA was stained with Hoechst. The visualization and scoring of the cells was done using an automated fluorescent microscope coupled with proprietary automated image analysis software provided by Cellomics (Pittsburg, PA). A total of 46 compounds were used in this evaluation, including 8 aneugens and 25 clastogens with varied mechanisms of action. Thirteen non-genotoxic compounds were also included. The automated scoring had a sensitivity of 88% and a specificity of 100%, with a predictive value positive of 100% and a predictive value negative of 76%, compared to data from the literature that was obtained with manual scoring. We also describe the incorporation of a metabolic activation system using rat liver S9 homogenates, and the use of cell number counts as a cytotoxicity index which is complementary to the CBPI- (cytokinesis-block proliferation index) based index. Finally, we also discuss the potential for artefactual findings due to fluorescent precipitate, which should be carefully monitored to prevent false positive results. In conclusion, the automated in vitro micronucleus scoring is a valid alternative to the manual scoring of slides, and it has the advantage of generating data in a rapid and consistent manner, and with low compound requirements, which makes it well suited as a screening assay in the early stages of compound development.

QSAR modeling of in vitro inhibition of cytochrome P450 3A4

We report the QSAR modeling of cytochrome P450 3A4 (CYP3A4) enzyme inhibition using four large data sets of in vitro data. These data sets consist of marketed drugs and drug-like compounds all tested in four assays measuring the inhibition of the metabolism of four different substrates by the CYP3A4 enzyme. The four probe substrates are benzyloxycoumarin, testosterone, benzyloxyresorufin, and midazolam. We first show that using state-of-the-art QSAR modeling approaches applied to only one of these four data sets does not lead to predictive models that would be useful for in silico filtering of chemical libraries. We then present the development and the testing of a multiple pharmacophore hypothesis (MPH) that is formulated as a conceptual extension of the traditional QSAR approach to modeling the promiscuous binding of a large variety of drugs to CYP3A4. In the simplest form, the MPH approach takes advantage of the multiple substrate data sets and identifies the binding of test compounds as either proximal or distal relative to that of a given substrate. Application of the approach to the in silico filtering of test compounds for potential inhibitors of CYP3A4 is also presented. In addition to an improvement in the QSAR modeling for the inhibition of CYP3A4, the results from this modeling approach provide structural insights into the drug-enzyme interactions. The existence of multiple inhibition data sets in the BioPrint database motivates the original development of the concept of a multiple pharmacophore hypothesis and provides a unique opportunity for formulating alternative strategies of QSAR modeling of the inhibition of the in vitro metabolism of CYP3A4.

Lower serum osteocalcin levels in pregnant drug users and their newborns at the time of delivery

Osteocalcin was evaluated by radioimmunoassay at the time of delivery in mothers and in the umbilical arteries of newborns in a group of pregnant drug users (eight heroin users and seven cocaine users) and compared with findings from a group of normal mothers and their newborns (N = 18). Drug users had lower osteocalcin values than did the normal women (1.3 +/- 0.7 versus 2.7 +/- 0.8 ng/mL, P less than .001); and infants of drug users had lower values than normal infants (14.1 +/- 3.8 versus 19.0 +/- 4.0 ng/mL, P less than .005). The birth weights of drug users' infants were smaller (3160 +/- 402 versus 3591 +/- 374 g, P less than .05) and there was a significant negative correlation (P less than .001) between osteocalcin and drug intake during pregnancy, but no changes in osteocalcin dependent on the type of drug used. These results suggest a toxic effect of these drugs on the osteoblast, which could account for the lower birth weights and skeletal alterations reported in the infants of drug users.