Large-scale next generation sequencing based analysis of SLCO1B1 pharmacogenetics variants in the Saudi population.

BACKGROUND: SLCO1B1 plays an important role in mediating hepatic clearance of many different drugs including statins, angiotensin-converting enzyme inhibitors, chemotherapeutic agents and antibiotics. Several variants in SLCO1B1 have been shown to have a clinically significant impact, in relation to efficacy of these medications. This study provides a comprehensive overview of SLCO1B1 variation in Saudi individuals, one of the largest Arab populations in the Middle East. METHODS: The dataset of 11,889 (9,961 exomes and 1,928 pharmacogenetic gene panel) Saudi nationals, was used to determine the presence and frequencies of SLCO1B1 variants, as described by the Clinical Pharmacogenetic Implementation Consortium (CPIC). RESULTS: We identified 141 previously described SNPs, of which rs2306283 (50%) and rs4149056 (28%), were the most common. In addition, we observed six alleles [*15 (24.7%) followed by *20 (8.04%), *14 (5.86%), *5 (3.84%), *31 (0.21%) and *9 (0.03%)] predicted to be clinically actionable. Allele diplotype to phenotype conversion revealed 41 OATP1B1 diplotypes. We estimated the burden of rare, and novel predicted deleterious variants, resulting from 17 such alterations. CONCLUSIONS: The data we present, from one of the largest Arab cohorts studied to date, provides the most comprehensive overview of SLCO1B1 variants, and the subsequent OATP1B1 activity of this ethnic group, which thus far remains relatively underrepresented in available international genomic databases. We believe that the presented data provides a basis for further clinical investigations and the application of personalized statin drug therapy guidance in Arabs.

The reduced function allele SLCO1B1 c.521T>C is of no practical relevance for the renal graft function over the first post-transplant year in patients treated with mycophenolic acid.

OBJECTIVE: It is unclear whether renal transplant recipients treated with mycophenolic acid (MPA) who carry the reduced-function allele at polymorphism SLCO1B1 c.521T>C differ from their wild-type peers regarding renal outcomes and tolerability. We aimed to estimate the effect of this polymorphism on the graft function (estimated glomerular filtration rate, eGFR) over the first 12 post-transplant months in patients on MPA-based maintenance immunosuppression. METHODS: In a 12-month observational cohort study, consecutive adult patients were repeatedly assessed for eGFR. The SLCO1B1 c.521C>T variant allele carriers (exposed) and wild-type subjects (controls) were balanced on a range of demographic, medical, and genetic variables at baseline, and eGFR trajectory was estimated with further adjustment for time-varying covariates. A subset of patients were assessed for exposure to MPA 5-7 days after the transplantation. RESULTS: The adjusted eGFR slopes from day 1 to day 28 (daily), and from day 28 to day 365 (monthly) were practically identical in exposed (n = 86) and control (n = 168) patients [geometric means ratios (GMR) = 0.99, 95% confidence interval (CI) = 0.92-1.06 and GMR = 0.98, 0.94-1.01, respectively]. The rates of adverse renal outcomes and possible MPA-related adverse effects were low, and similar in exposed and controls [rate ratios (RR) = 0.94, 0.49-1.84 and RR = 1.08, 0.74-1.58, respectively]. The pharmacokinetic analysis did not signal meaningful differences regarding exposure to MPA, overall (exposed n = 23, control n = 45), if cotreated with cyclosporine (n = 17 vs. n = 26) or with tacrolimus (n = 8 vs. n = 17). CONCLUSIONS: In patients treated with MPA, variant allele SLCO1B1 c.521T>C appears of no practical relevance regarding the 12-month renal graft function, MPA safety and exposure to MPA at early steady-state.

Precision medicine in cardiovascular therapeutics: Evaluating the role of pharmacogenetic analysis prior to drug treatment.

Pharmacogenomics is the examination of how genetic variation influences drug metabolism and response, in terms of both efficacy and safety. In cardiovascular disease, patient-specific diplotypes determine phenotypes, thereby influencing the efficacy and safety of drug treatments, including statins, antiarrhythmics, anticoagulants and antiplatelets. Notably, polymorphisms in key genes, such as CYP2C9, CYP2C19, VKORC1 and SLCO1B1, significantly impact the outcomes of treatment with clopidogrel, warfarin and simvastatin. Furthermore, the CYP2C19 polymorphism influences the pharmacokinetics and safety of the novel hypertrophic cardiomyopathy inhibitor, mavacamten. In this review, we critically assess the clinical application of pharmacogenomics in cardiovascular disease and delineate present and future utilization of pharmacogenomics. This includes insights into identifying missing heritability, the integration of whole genome sequencing and the application of polygenic risk scores to enhance the precision of personalized drug therapy. Our discussion encompasses health economic analyses that underscore the cost benefits associated with pre-emptive genotyping for warfarin and clopidogrel treatments, albeit acknowledging the need for further research in this area. In summary, we contend that cardiovascular pharmacogenomic analyses are underpinned by a wealth of evidence, and implementation is already occurring for some of these gene-drug pairs, but as with any area of medicine, we need to continually gather more information to optimize the use of pharmacogenomics in clinical practice.

Genetic and clinical predictors of rifapentine and isoniazid pharmacokinetics in paediatrics with tuberculosis infection.

OBJECTIVES: Twelve weekly doses of rifapentine and isoniazid (3HP regimen) are recommended for TB preventive therapy in children with TB infection. However, they present with variability in the pharmacokinetic profiles. The current study aimed to develop a pharmacokinetic model of rifapentine and isoniazid in 12 children with TB infection using NONMEM. METHODS: Ninety plasma and 41 urine samples were collected at Week 4 of treatment. Drug concentrations were measured using a validated HPLC-UV method. MassARRAY® SNP genotyping was used to investigate genetic factors, including P-glycoprotein (ABCB1), solute carrier organic anion transporter B1 (SLCO1B1), arylacetamide deacetylase (AADAC) and N-acetyl transferase (NAT2). Clinically relevant covariates were also analysed. RESULTS: A two-compartment model for isoniazid and a one-compartment model for rifapentine with transit compartment absorption and first-order elimination were the best models for describing plasma and urine data. The estimated (relative standard error, RSE) of isoniazid non-renal clearance was 3.52 L·h-1 (23.1%), 2.91 L·h-1 (19.6%), and 2.58 L·h-1 (20.0%) in NAT2 rapid, intermediate and slow acetylators. A significant proportion of the unchanged isoniazid was cleared renally (2.7 L·h-1; 8.0%), while the unchanged rifapentine was cleared primarily through non-renal routes (0.681 L·h-1; 3.6%). Participants with the ABCB1 mutant allele had lower bioavailability of rifapentine, while food prolonged the mean transit time of isoniazid. CONCLUSIONS: ABCB1 mutant allele carriers may require higher rifapentine doses; however, this must be confirmed in larger trials. Food did not affect overall exposure to isoniazid and only delayed absorption time.

Mechanistic Account of Distinct Change in Organic Anion Transporting Polypeptide 1B (OATP1B) Substrate Pharmacokinetics during OATP1B-Mediated Drug-Drug Interactions Using Physiologically Based Pharmacokinetic Modeling.

The role of transporters in drug clearance is widely acknowledged, directly and indirectly by facilitating tissue/enzyme exposure. Through the latter, transporters also affect volume of distribution. Drug-drug interactions (DDIs) involving organic anion transporting polypeptides (OATPs) 1B1/1B3 and SLCO1B1 pharmacogenetics lead to altered pharmacokinetics of OATP1B substrates; however, several factors may confound direct interpretation of pharmacokinetic parameters from these clinical studies using noncompartmental analysis (NCA). A review of clinical data herein indicates a single dose of OATP1B inhibitor rifampin almost never leads to increased substrate half-life but often a decrease and that most clinical OATP1B substrates are CYP3A4 substrates and/or undergo enterohepatic cycling (EHC). Using hypothetically simple OATP1B substrate physiologically based pharmacokinetic (PBPK) models, simulated effect of rifampin differed from specific OATP1B inhibition due to short rifampin half-life causing dissipation of OATP1B inhibition over time combined with CYP3A4 induction. Calculated using simulated tissue data, volume of distribution indeed decreased with OATP1B inhibition and was expectedly limited to the contribution of liver volume. However, an apparent and counterintuitive effect of rifampin on volume greater than that on clearance resulted for CYP3A4 substrates using NCA. The effect of OATP1B inhibition and rifampin on OATP1B substrate models incorporating EHC plus or minus renal clearance was distinct compared with simpler models. Using PBPK models incorporating reversible lactone metabolism for clinical OATP1B substrates atorvastatin and pitavastatin, DDIs reporting decreased half-life with rifampin were reproduced. These simulations provide an explanation for the distinct change in OATP1B substrate pharmacokinetics observed in clinical studies, including changes in volume of distribution and additional mechanisms. SIGNIFICANCE STATEMENT: Transporters are involved in drug clearance and volume of distribution, and distinct changes in OATP1B substrate pharmacokinetics are observed with OATP1B inhibitor rifampin. Using hypothetical and validated PBPK models and simulations, this study addresses the limitations of single-dose rifampin and complicated clinical OATP1B substrate disposition in evaluating the pharmacokinetic parameters of OATP1B substrates during rifampin drug-drug interactions (DDIs). These models account for change in volume of distribution and identify additional mechanisms underlying apparent pharmacokinetic changes in OATP1B DDIs.

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.

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.

Interplay of Ritonavir-Boosted Oral Cabazitaxel with the Organic Anion-Transporting Polypeptide (OATP) Uptake Transporters and Carboxylesterase 1 in Mice.

Intravenously administered chemotherapeutic cabazitaxel is used for palliative treatment of prostate cancer. An oral formulation would be more patient-friendly and reduce the need for hospitalization. We therefore study determinants of the oral pharmacokinetics of cabazitaxel in a ritonavir-boosted setting, which reduces the CYP3A-mediated first-pass metabolism of cabazitaxel. We here assessed the role of organic anion-transporting polypeptides (OATPs) in the disposition of orally boosted cabazitaxel and its active metabolites, using the Oatp1a/b-knockout and the OATP1B1/1B3-transgenic mice. These transporters may substantially affect plasma clearance and hepatic and intestinal drug disposition. The pharmacokinetics of cabazitaxel and DM2 were not significantly affected by Oatp1a/b and OATP1B1/1B3 activity. In contrast, the plasma AUC0-120 min of DM1 in Oatp1a/b-/- was 1.9-fold (p < 0.05) higher than that in wild-type mice, and that of docetaxel was 2.4-fold (p < 0.05) higher. We further observed impaired hepatic uptake and intestinal disposition for DM1 and docetaxel in the Oatp-ablated strains. None of these parameters showed rescue by the OATP1B1 or -1B3 transporters in the humanized mouse strains, suggesting a minimal role of OATP1B1/1B3. Ritonavir itself was also a potent substrate for mOatp1a/b, showing a 2.9-fold (p < 0.0001) increased plasma AUC0-120 min and 3.5-fold (p < 0.0001) decreased liver-to-plasma ratio in Oatp1a/b-/- compared to those in wild-type mice. Furthermore, we observed the tight binding of cabazitaxel and its active metabolites, including docetaxel, to plasma carboxylesterase (Ces1c) in mice, which may complicate the interpretation of pharmacokinetic and pharmacodynamic mouse studies. Collectively, these results will help to further optimize (pre)clinical research into the safety and efficacy of orally applied cabazitaxel.

G45 and V386 in Transmembrane Domains 1 and 8 Are Critical for the Activation of OATP1B3-Mediated E17ßG Uptake by Clotrimazole.

Organic anion-transporting polypeptides (OATPs) 1B1 and 1B3 are two highly homologous transport proteins. However, OATP1B1- and 1B3-mediated estradiol-17ß-glucuronide (E17ßG) uptake can be differentially affected by clotrimazole. In this study, by functional characterization on chimeric transporters and single mutants, we find that G45 in transmembrane domain 1 (TM1) and V386 in TM8 are critical for the activation of OATP1B3-mediated E17ßG uptake by clotrimazole. However, the effect of clotrimazole on the function of OATP1B3 is substrate-dependent as clotrimazole does not stimulate OATP1B3-mediated uptake of 4',5'-dibromofluorescein (DBF) and rosuvastatin. In addition, clotrimazole is not transported by OATP1B3, but it can efficiently permeate the plasma membrane due to its lipophilic properties. Homology modeling and molecular docking indicate that E17ßG binds in a substrate binding pocket of OATP1B3 through hydrogen bonding and hydrophobic interactions, among which its sterol scaffold forms hydrophobic contacts with V386. In addition, a flexible glycine residue at position 45 is essential for the activation of OATP1B3. Finally, clotrimazole is predicted to bind at an allosteric site, which mainly consists of hydrophobic residues located at the cytoplasmic halves of TMs 4, 5, 10, and 11.

Associations between UGT1A1, SLCO1B1, SLCO1B3, BLVRA and HMOX1 polymorphisms and susceptibility to neonatal severe hyperbilirubinemia in Chinese Han population.

BACKGROUND: Severe neonatal hyperbilirubinemia could lead to kernicterus and neonatal death. This study aimed to analyze the association between single nucleotide polymorphisms in genes involved in bilirubin metabolism and the incidence of severe hyperbilirubinemia. METHODS: A total of 144 neonates with severe hyperbilirubinemia and 50 neonates without or mild hyperbilirubinemia were enrolled in 3 institutions between 2019 and 2020. Twelve polymorphisms of 5 genes (UGT1A1, SLCO1B1, SLCO1B3, BLVRA, and HMOX1) were analyzed by PCR amplification of genomic DNA. Genotyping was performed using an improved multiplex ligation detection reaction technique based on ligase detection reaction. RESULTS: The frequencies of the A allele in UGT1A1-rs4148323 and the C allele in SLCO1B3-rs2417940 in the severe hyperbilirubinemia group (30.2% and 90.6%, respectively) were significantly higher than those in the controls (30.2% vs.13.0%, 90.6% vs. 78.0%, respectively, both p < 0.05). Haplotype analysis showed the ACG haplotype of UGT1A1 were associated with an increased hyperbilirubinemia risk (OR 3.122, p = 0.001), whereas the GCG haplotype was related to a reduced risk (OR 0.523, p = 0.018). CONCLUSION: The frequencies of the A allele in rs4148323 and the C allele in rs2417940 are highly associated with the incidence of severe hyperbilirubinemia in Chinese Han neonates. TRIAL REGISTRATION: Trial registration number:ChiCTR1800020424; Date of registration:2018-12-29.

SLCO1B1 Exome Sequencing and Statin Treatment Response in 64,000 UK Biobank Patients.

The solute carrier organic anion transporter family member 1B1 (SLCO1B1) encodes the organic anion-transporting polypeptide 1B1 (OATP1B1 protein) that transports statins to liver cells. Common genetic variants in SLCO1B1, such as *5, cause altered systemic exposure to statins and therefore affect statin outcomes, with potential pharmacogenetic applications; yet, evidence is inconclusive. We studied common and rare SLCO1B1 variants in up to 64,000 patients from UK Biobank prescribed simvastatin or atorvastatin, combining whole-exome sequencing data with up to 25-year routine clinical records. We studied 51 predicted gain/loss-of-function variants affecting OATP1B1. Both SLCO1B1*5 alone and the SLCO1B1*15 haplotype increased LDL during treatment (beta*5 = 0.08 mmol/L, p = 6 × 10-8; beta*15 = 0.03 mmol/L, p = 3 × 10-4), as did the likelihood of discontinuing statin prescriptions (hazard ratio*5 = 1.12, p = 0.04; HR*15 = 1.05, p = 0.04). SLCO1B1*15 and SLCO1B1*20 increased the risk of General Practice (GP)-diagnosed muscle symptoms (HR*15 = 1.22, p = 0.003; HR*20 = 1.25, p = 0.01). We estimated that genotype-guided prescribing could potentially prevent 18% and 10% of GP-diagnosed muscle symptoms experienced by statin patients, with *15 and *20, respectively. The remaining common variants were not individually significant. Rare variants in SLCO1B1 increased LDL in statin users by up to 1.05 mmol/L, but replication is needed. We conclude that genotype-guided treatment could reduce GP-diagnosed muscle symptoms in statin patients; incorporating further SLCO1B1 variants into clinical prediction scores could improve LDL control and decrease adverse events, including discontinuation.

Sex-Hormone-Binding Globulin Gene Polymorphisms and Breast Cancer Risk in Caucasian Women of Russia.

In our work, the associations of GWAS (genome-wide associative studies) impact for sex-hormone-binding globulin (SHBG)-level SNPs with the risk of breast cancer (BC) in the cohort of Caucasian women of Russia were assessed. The work was performed on a sample of 1498 women (358 BC patients and 1140 control (non BC) subjects). SHBG correlated in previously GWAS nine polymorphisms such as rs780093 GCKR, rs17496332 PRMT6, rs3779195 BAIAP2L1, rs10454142 PPP1R21, rs7910927 JMJD1C, rs4149056 SLCO1B1, rs440837 ZBTB10, rs12150660 SHBG, and rs8023580 NR2F2 have been genotyped. BC risk effects of allelic and non-allelic SHBG-linked gene SNPs interactions were detected by regression analysis. The risk genetic factor for BC developing is an SHBG-lowering allele variant C rs10454142 PPP1R21 ([additive genetic model] OR = 1.31; 95%CI = 1.08-1.65; pperm = 0.024; power = 85.26%), which determines 0.32% of the cancer variance. Eight of the nine studied SHBG-related SNPs have been involved in cancer susceptibility as part of nine different non-allelic gene interaction models, the greatest contribution to which is made by rs10454142 PPP1R21 (included in all nine models, 100%) and four more SNPs-rs7910927 JMJD1C (five models, 55.56%), rs17496332 PRMT6 (four models, 44.44%), rs780093 GCKR (four models, 44.44%), and rs440837 ZBTB10 (four models, 44.44%). For SHBG-related loci, pronounced functionality in the organism (including breast, liver, fibroblasts, etc.) was predicted in silico, having a direct relationship through many pathways with cancer pathophysiology. In conclusion, our results demonstrated the involvement of SHBG-correlated genes polymorphisms in BC risk in Caucasian women in Russia.

[Clinical characteristics and genetic analysis of four children with Rotor syndrome].

OBJECTIVE: To explore the characteristics of SLCO1B1/SLCO1B3 gene variants among children with Rotor syndrome (RS). METHODS: Four children who were admitted to the Department of Hepatology of Hunan Children's Hospital between January 2019 and January 2022 were selected as the study subjects. Trio-whole exome sequencing was carried out for the four families, and gel electrophoresis was used to verify an insertional variant of long-interspersed element-1 (LINE-1). RESULTS: Genetic testing has identified three variants of the SLCO1B1 gene, including c.1738C>T (p.R580*), c.757C>T (p.R253*) and c.1622A>C (p.Q541P), and two variants of the SLCO1B3 gene, including c.481+22insLINE-1 and c.1747+1G>A among the children. Three of them were found to harbor homozygous variants of the SLCO1B1/SLCO1B3 genes, and one has harbored compound heterozygous variants. Sanger sequencing confirmed the existence of all variants, and gel electrophoresis has confirmed the existence of the LINE-1 insertional variant of about 6 kb within intron 6 of the SLCO1B3 gene in all children. CONCLUSION: The pathogenesis of the RS among the four children may be attributed to the variants of the SLCO1B1/SLCO1B3 genes. The LINE-1 insertion variant of the SLCO1B3 gene may be common among Chinese RS patients.

The Effect of Original Russian Neurotropic Drugs on Organic Anion Transporting Polypeptides OATP1B1 and OATP1B3.

In experiments on HepG2 cells, we studied the effect of the original domestic neurotropic drugs omberacetam, fabomotizole, and ethylmethylhydroxypyridine succinate (EMHPS) (1-500 µM) on the activity and content of organic anion transporting polypeptides OATP1B1 and OATP1B3. It was shown that omberacetam (500 µM) increased the content of OATP1B1 and OATP1B3, fabomotizole did not affect the level of both transporters, and EMHPS (500 µM) increased the content of OATP1B1 compared to the control and did not affect the level of OATP1B3. The tested substances also reduced the OATP1B1/OATP1B3 ratio, as evidenced by a decrease in the penetration of atorvastatin, a substrate of the transporters, into HepG2 cells in the presence of omberacetam (100-500 µM), fabomotizole (500 µM), and EMHPS (10-500 µM). Evaluation of clinical significance of the obtained results, according to the FDA approach based on the calculation of the Cmax/IC50 ratio, showed that the effect of the tested substances on OATP1B1/OATP1B3 is clinically insignificant.

Physiologically Based Pharmacokinetic Model of OATP1B Substrates with a Nonlinear Mixed Effect Approach: Estimating Empirical In Vitro-to-In Vivo Scaling Factors.

BACKGROUND AND OBJECTIVE: Physiologically based pharmacokinetic (PBPK) models are valuable for translating in vitro absorption, distribution, metabolism, and excretion (ADME) data to predict clinical pharmacokinetics, and can enable discovery and early clinical stages of pharmaceutical research. However, in predicting pharmacokinetics of organic anion transporting polypeptide (OATP) 1B substrates based on in vitro transport and metabolism data, PBPK models typically require additional empirical in vitro-to-in vivo scaling factors (ESFs) in order to accurately recapitulate observed clinical profiles. As model simulation is very sensitive to ESFs, a critical evaluation of ESF estimation is prudent. Previously studies have applied classic 'two-stage' and 'naïve pooled data' approaches in identifying a set of compound independent ESFs. However, the 'two-stage' approach has the parameter identification issue in separately fitting data for individual compounds, while the 'naïve pooled data' approach ignores interstudy variability, leading to potentially biased ESF estimates. METHODS: In this study, we have applied a nonlinear mixed-effect approach in estimating ESF of the PBPK model and incorporated additional data from 86 runs of in vitro uptake assay and 49 clinical studies of 12 training compounds in model development to further enhance the translation of in vitro data to predict the pharmacokinetics of OATP1B substrate drugs. To test predication accuracy of the model, a 'leave-one-out' analysis has been performed. RESULTS: The established model can reasonably describe the clinical observations, with both mean values and interstudy variabilities quantified for ESF and volume of distribution parameters. The mean estimates are largely consistent with values in the previous reports. The interstudy variabilities of these parameters are estimated to be at least 50% (as coefficient of variation). Most compounds can be reasonably predicted in the 'leave-one-out' analysis. CONCLUSION: This study improves the confidence in predicting the pharmacokinetics of OATP1B substrates in individual studies of small sample sizes, and quantifies the variability associated with the prediction.

Discovery of TP0628103: A Highly Potent and Selective MMP-7 Inhibitor with Reduced OATP-Mediated Clearance Designed by Shifting Isoelectric Points.

Matrix metalloproteinase-7 (MMP-7) has been shown to play an important role in pathophysiological processes such as cancer and fibrosis. We previously discovered selective MMP-7 inhibitors by molecular hybridization and structure-based drug design. However, the systemic clearance (CLtot) of the biologically active lead compound was very high. Because our studies revealed that hepatic uptake by organic anion transporting polypeptide (OATP) was responsible for the high CLtot, we found a novel approach to reducing their uptake based on isoelectric point (IP) values as an indicator for substrate recognition by OATP1B1/1B3. Our "IP shift strategy" to adjust the IP values culminated in the discovery of TP0628103 (18), which is characterized by reduced in vitro OATP-mediated hepatic uptake and in vivo CLtot. Our in vitro-in vivo extrapolation of OATP-mediated clearance and the "IP shift strategy" provide crucial insights for a new medicinal chemistry approach to reducing the systemic clearance of OATP1B1/1B3 substrates.

Non-targeted metabolomics for the identification of plasma metabolites associated with organic anion transporting polypeptide 1B1 function.

Our aim was to evaluate biomarkers for organic anion transporting polypeptide 1B1 (OATP1B1) function using a hypothesis-free metabolomics approach. We analyzed fasting plasma samples from 356 healthy volunteers using non-targeted metabolite profiling by liquid chromatography high-resolution mass spectrometry. Based on SLCO1B1 genotypes, we stratified the volunteers to poor, decreased, normal, increased, and highly increased OATP1B1 function groups. Linear regression analysis, and random forest (RF) and gradient boosted decision tree (GBDT) regressors were used to investigate associations of plasma metabolite features with OATP1B1 function. Of the 9152 molecular features found, 39 associated with OATP1B1 function either in the linear regression analysis (p < 10-5) or the RF or GBDT regressors (Gini impurity decrease > 0.01). Linear regression analysis showed the strongest associations with two features identified as glycodeoxycholate 3-O-glucuronide (GDCA-3G; p = 1.2 × 10-20 for negative and p = 1.7 × 10-19 for positive electrospray ionization) and one identified as glycochenodeoxycholate 3-O-glucuronide (GCDCA-3G; p = 2.7 × 10-16). In both the RF and GBDT models, the GCDCA-3G feature showed the strongest association with OATP1B1 function, with Gini impurity decreases of 0.40 and 0.17. In RF, this was followed by one GDCA-3G feature, an unidentified feature with a molecular weight of 809.3521, and the second GDCA-3G feature. In GBDT, the second and third strongest associations were observed with the GDCA-3G features. Of the other associated features, we identified with confidence two representing lysophosphatidylethanolamine 22:5. In addition, one feature was putatively identified as pregnanolone sulfate and one as pregnenolone sulfate. These results confirm GCDCA-3G and GDCA-3G as robust OATP1B1 biomarkers in human plasma.

Interactions of oleanane pentacyclic triterpenoids with human organic anion transporting polypeptide 1B1 and 1B3.

Oleanane pentacyclic triterpenoids have been widely used in clinical practice. However, studies on their interactions with hepatic transporters remain limited. In this study, we systematically investigated the inhibitory effects of 14 oleanane pentacyclic triterpenoids on organic anion transporting polypeptide 1B1 and 1B3 (OATP1B1 and OATP1B3), two liver-specific uptake transporters. Through fluorescence-based cellular uptake assays, we identified three potent OATP1B1 inhibitors (saikosaponin B1, saikosaponin A and 18ß-glycyrrhetinic acid) and five potent OATP1B3 inhibitors (echinocystic acid, 3-oxo-16α-hydroxy-olean-12-en-28ß-oic acid, chikusetsu saponin IVa, saikosaponin B1 and 18ß-glycyrrhetinic acid). Structural analysis revealed that free oleanane triterpenoids inhibited OATP1B1/1B3 more potently than triterpene glycosides. Despite their similar structures, 18ß-glycyrrhetinic acid exhibited much stronger inhibition on OATP1B1/1B3 than 18α-glycyrrhetinic acid, while both were substrates of OATP1B3. Interestingly, OATP1B3 overexpression significantly increased reactive oxygen species (ROS) levels in HepG2 cells after treatment with 18ß-glycyrrhetinic acid. To conclude, this study highlights the potential interactions of oleanane pentacyclic triterpenoids with OATP1B1/1B3, and provides novel insights into the anti-cancer activity of 18ß-glycyrrhetinic acid.

Physiologically based pharmacokinetic (PBPK) modeling of pitavastatin in relation to SLCO1B1 genetic polymorphism.

Pitavastatin, a potent 3-hydroxymethylglutaryl coenzyme A reductase inhibitor, is indicated for the treatment of hypercholesterolemia and mixed dyslipidemia. Hepatic uptake of pitavastatin is predominantly occupied by the organic anion transporting polypeptide 1B1 (OATP1B1) and solute carrier organic anion transporter family member 1B1 (SLCO1B1) gene, which is a polymorphic gene that encodes OATP1B1. SLCO1B1 genetic polymorphism significantly alters the pharmacokinetics of pitavastatin. This study aimed to establish the physiologically based pharmacokinetic (PBPK) model to predict pitavastatin pharmacokinetics according to SLCO1B1 genetic polymorphism. PK-Sim® version 10.0 was used to establish the whole-body PBPK model of pitavastatin. Our pharmacogenomic data and a total of 27 clinical pharmacokinetic data with different dose administration and demographic properties were used to develop and validate the model, respectively. Physicochemical properties and disposition characteristics of pitavastatin were acquired from previously reported data or optimized to capture the plasma concentration-time profiles in different SLCO1B1 diplotypes. Model evaluation was performed by comparing the predicted pharmacokinetic parameters and profiles to the observed data. Predicted plasma concentration-time profiles were visually similar to the observed profiles in the non-genotyped populations and different SLCO1B1 diplotypes. All fold error values for AUC and Cmax were included in the two fold range of observed values. Thus, the PBPK model of pitavastatin in different SLCO1B1 diplotypes was properly established. The present study can be useful to individualize the dose administration strategy of pitavastatin in individuals with various ages, races, and SLCO1B1 diplotypes.

The impact of SLCO1B1 rs4149056 on LDL-C target achievement after lipid lowering therapy optimization in men and women with familial hypercholesterolemia.

Background and aims: FH women are less likely to receive intensive statin treatment and to obtain a 50% reduction of LDL-C from baseline compared to men with FH. SLCO1B1 rs4149056 might influence statin therapy compliance and thus LDL-C target achievement. Our aim was to evaluate the impact of SLCO1B1 rs4149056 on LDL-C target achievement after lipid lowering therapy (LLT) optimization in men and women with FH. Methods: This was a retrospective observational study involving 412 FH subjects with a probable or defined clinical diagnosis of FH who had had genetic analysis from June 2016 to September 2022. Biochemical analysis was obtained from all subjects at baseline and at the last follow-up after LLT optimization. Results: After LLT optimization the percentage of FH subjects on high-intensity statins decreased from the M/SLCO1B1- group to the W/SLCO1B1+ group and the same was found in LDL-C target distribution (for both p for trend < 0.01). The prevalence of SASE fear increased from the M/SLCO1B1- group to the W/SLCO1B1+ group and the same was observed in reported myalgia distribution (for both p for trend < 0.01). Logistic regression analysis showed that the W/SCLO1B1-, M/SCLO1B1+ and W/SCLO1B1+ groups were inversely associated with LDL-C target achievement (p for trend < 0.001) and the W/SCLO1B1+ group exhibited the strongest association. Conclusion: A low prevalence of FH women with SLCO1B1 rs4149056 were on high intensity statins and they rarely achieved LDL-C target. The genotype effect of SLCO1B1 rs4149056 could be more pronounced in FH women than men.