Different effects of SLCO1B1 polymorphism on the pharmacokinetics of atorvastatin and rosuvastatin

Association of SLCO1B1 gene variants with angiotensin-converting enzyme inhibitor-induced cough in a Pakistani hypertensive cohort

Background

Angiotensin-converting enzyme inhibitors (ACEIs) are prescribed for individuals with high cardiovascular (CV) risk; however, persistent cough limits the use of ACEIs in a large number of patients. The current study aimed to identify the genetic variants in the SLCO1B1 gene that might be associated with ACEI-related cough in a Pakistani hypertensive population.

Methods

A prospective cohort study was conducted at a tertiary care hospital in Pakistan. A total of 74 patients who had been treated with ACEIs were recruited through a convenient sampling method. The study was approved by the Institutional Review Board & Ethics Committee of the Shifa International Hospital, Islamabad. Patients provided 2 ml of blood for sequencing after signing informed consent. Partial gene sequencing of SLCO1B1 was carried out to find single nucleotide polymorphisms (SNPs) and haplotypes.

Results

It was found, through a structured questionnaire, that thirty-eight (38) patients experienced cough within 2 weeks of ACEI administration and were considered as a case group (cough), and thirty-six (36) patients were considered as a control group (no cough). The incidence of cough was 51%. We found six different SNPs and 9 haplotypes in the partial gene sequences of SLCO1B1. Haplotype H4 was associated significantly with cough after adjusting for sex and smoking status. Other SNPs and haplotypes were not significantly associated with ACE-Is-induced cough.

Conclusion

These findings emphasize the significance of SLCO1B1 genetic variants, specifically H4, as a potential predictor of ACEI-induced cough. It could be included in clinical practice as a possible risk factor for ACEI-induced cough once confirmed in larger clinical trials with bigger sample sizes. The replication of these findings in larger and more diverse populations is likely to contribute to the therapeutic use of ACEIs by predicting ACEI-induced cough.

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.

Validation of a novel direct method to determine reduced adherence to atorvastatin therapy

AIMS

Objective methods to determine statin adherence are requested to improve lipid management. We have recently established a method to detect reduced adherence to atorvastatin therapy with cut-off values based on the sum of atorvastatin and its major metabolites in the blood. We aimed to validate this method in patients with and without cardiovascular disease, and optimize previous cut-off values.

METHODS AND RESULTS

The pharmacokinetic study included 60 participants treated with atorvastatin 20 mg (N = 20), 40 mg (N = 20), and 80 mg (N = 20). Atorvastatin was then stopped and blood samples collected from day zero to day four. Quantification of the parent drug and its metabolites in blood plasma was performed with a liquid chromatography-tandem mass spectrometry assay. The cut-off values for reduced adherence were validated and optimized by calculating diagnostic sensitivity and specificity. Our candidate cut-off value of dose-normalized six-component sum of atorvastatin plus metabolites <0.10 nM/mg provided a sensitivity of 97% and a specificity of 93% for detecting ≥2 omitted doses. An optimized cut-off <0.062 nM/mg provided a sensitivity of 90% and a specificity of 100%. An alternative simplified two-component metabolite sum with a cut-off value <0.05 nM/mg provided a sensitivity of 98% and a specificity of 76%. An optimized cut-off <0.02 nM/mg provided a sensitivity of 97% and a specificity of 98%.

CONCLUSION

This validation study confirms that our direct method discriminates reduced adherence from adherence to atorvastatin therapy with high diagnostic accuracy. The method may improve lipid management in clinical practice and serve as a useful tool in future studies.

Genome-Wide Association Study of Atorvastatin Pharmacokinetics: Associations With SLCO1B1, UGT1A3, and LPP

In a genome-wide association study of atorvastatin pharmacokinetics in 158 healthy volunteers, the SLCO1B1 c.521T>C (rs4149056) variant associated with increased area under the plasma concentration-time curve from time zero to infinity (AUC0-∞) of atorvastatin (P = 1.2 × 10-10), 2-hydroxy atorvastatin (P = 4.0 × 10-8), and 4-hydroxy atorvastatin (P = 2.9 × 10-8). An intronic LPP variant, rs1975991, associated with reduced atorvastatin lactone AUC0-∞ (P = 3.8 × 10-8). Three UGT1A variants linked with UGT1A3*2 associated with increased 2-hydroxy atorvastatin lactone AUC0-∞ (P = 3.9 × 10-8). Furthermore, a candidate gene analysis including 243 participants suggested that increased function SLCO1B1 variants and decreased activity CYP3A4 variants affect atorvastatin pharmacokinetics. Compared with individuals with normal function SLCO1B1 genotype, atorvastatin AUC0-∞ was 145% (90% confidence interval: 98-203%; P = 5.6 × 10-11) larger in individuals with poor function, 24% (9-41%; P = 0.0053) larger in those with decreased function, and 41% (16-59%; P = 0.016) smaller in those with highly increased function SLCO1B1 genotype. Individuals with intermediate metabolizer CYP3A4 genotype (CYP3A4*2 or CYP3A4*22 heterozygotes) had 33% (14-55%; P = 0.022) larger atorvastatin AUC0-∞ than those with normal metabolizer genotype. UGT1A3*2 heterozygotes had 16% (5-25%; P = 0.017) smaller and LPP rs1975991 homozygotes had 34% (22-44%; P = 4.8 × 10-5) smaller atorvastatin AUC0-∞ than noncarriers. These data demonstrate that genetic variation in SLCO1B1, UGT1A3, LPP, and CYP3A4 affects atorvastatin pharmacokinetics. This is the first study to suggest that LPP rs1975991 may reduce atorvastatin exposure. [Correction added on 6 April, after first online publication: An incomplete sentence ("= 0.017) smaller in heterozygotes for UGT1A3*2 and 34% (22%, 44%; P × 10-5) smaller in homozygotes for LPP noncarriers.") has been corrected in this version.].

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.

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.

Quantitative prediction of transporter-mediated drug-drug interactions using the mechanistic static pharmacokinetic (MSPK) model

Guidance/guidelines on drug-drug interactions (DDIs) have been issued in Japan, the United States, and Europe. These guidance/guidelines provide decision trees for conducting metabolizing enzyme-mediated clinical DDI studies; however, the decision trees for transporter-mediated DDIs lack quantitative prediction methods. In this study, the accuracy of a net-effect mechanistic static pharmacokinetics (MSPK) model containing the fraction transported (ft) of transporters was examined to predict transporter-mediated DDIs. This study collected information on 25 oral drugs with new active reagents that were used in clinical DDI studies as perpetrators (42 cases) from drugs approved in Japan between April 2016 and June 2020. The AUCRs (AUC ratios with and without perpetrators) of victim drugs were predicted using the net-effect MSPK model. As a result, 83 and 95% of the predicted AUCRs were within 1.5- and 2-fold error in the observed AUCRs, respectively. In cases where the victims were statins in which pharmacokinetics several transporters are involved, 70 and 91% of the predicted AUCRs were within 1.5- and 2-fold errors, respectively. Therefore, the net-effect MSPK model was applicable for predicting the AUCRs of victims, which are substrates for multiple transporters.

Atorvastatin-associated myotoxicity: A toxicokinetic review of pharmacogenetic associations to evaluate the feasibility of precision pharmacotherapy

Atorvastatin (ATV) and other statins are highly effective in reducing cholesterol levels. However, in some patients, the development of drug-associated muscle side effects remains an issue as it compromises the adherence to treatment. Since the toxicity is dose-dependent, exploring factors modulating pharmacokinetics (PK) appears fundamental. The purpose of this review aims at reporting the current state of knowledge about the singular genetic susceptibilities influencing the risk of developing ATV muscle adverse events through PK modulations. Multiple single nucleotide polymorphisms (SNP) in efflux (ABCB1, ABCC1, ABCC2, ABCC4 and ABCG2) and influx (SLCO1B1, SLCO1B3 and SLCO2B1) transporters have been explored for their association with ATV PK modulation or with statin-related myotoxicities (SRM) development. The most convincing pharmacogenetic association with ATV remains the influence of the rs4149056 (c.521 T > C) in SLCO1B1 on ATV PK and pharmacodynamics. This SNP has been robustly associated with increased ATV systemic exposure and consequently, an increased risk of SRM. Additionally, the SNP rs2231142 (c.421C > A) in ABCG2 has also been associated with increased drug exposure and higher risk of SRM occurrence. SLCO1B1 and ABCG2 pharmacogenetic associations highlight that modulation of ATV systemic exposure is important to explain the risk of developing SRM. However, some novel observations credit the hypothesis that additional genes (e.g. SLCO2B1 or ABCC1) might be important for explaining local PK modulations within the muscle tissue, indicating that studying the local PK directly at the skeletal muscle level might pave the way for additional understanding.

Relationship of plasma 3-carboxy-4-methyl-5-propyl-2-furanpropanoic acid concentration with OATP1B activity in patients with chronic kidney disease

Organic anion-transporting polypeptides (OATP)1B are drug transporters mainly expressed in the sinusoidal membrane. Many studies have suggested that OATP1B activity is affected by genetic factor, the uremic toxin 3-carboxy-4-methyl-5-propyl-2-furanpropanoic acid (CMPF), and inflammatory cytokines, such as tumor necrosis factor-α (TNF-α) and interleukin-6 (IL-6). Coproporphyrin-I (CP-I) is spotlighted as a highly accurate endogenous substrate of OATP1B. We previously reported a positive correlation between plasma CMPF and CP-I concentrations in patients with chronic kidney disease (CKD). The present study evaluated the impact of genetic polymorphisms, CMPF, IL-6, TNF-α, and estimated glomerular filtration rate (eGFR) on individual differences in OATP1B activity in patients with CKD. Seventy-three patients with CKD who received kidney transplant at least 3 months earlier were analyzed. Plasma CP-I concentration was higher in OATP1B1*15 carriers than in non-carriers. In all patients, CP-I did not correlate significantly with CMPF, IL-6, TNF-α, or eGFR. However, when the dataset was cut off at CMPF concentration of 8 and 7 μg/mL, 4 μg/mL, 3 μg/mL or 2 μg/mL, CMPF correlated positively with CP-I, and correlation coefficient tended to be higher as plasma CMPF concentration was lower. In conclusion, OATP1B1*15 impacted OATP1B activity in patients with CKD, but IL-6 and TNF-α did not. However, the impact of CMPF on OATP1B activity was limited to low CMPF concentrations, and the effect could be saturated at high concentrations. When prescribing an OATP1B substrate drug for patients with CKD, the OATP1B1*15 carrier status and plasma CMPF concentration may need to be considered to decide the dose regimen.

SLCO1B1 Genetic Variation Influence on Atorvastatin Systemic Exposure in Pediatric Hypercholesterolemia

This clinical study examined the influence of SLCO1B1 c.521T>C (rs4149056) on plasma atorvastatin concentrations in pediatric hypercholesterolemia. The participants (8-21 years), including heterozygous (c.521T/C, n = 13), homozygous (c.521C/C, n = 2) and controls (c.521T/T, n = 13), completed a single-oral-dose pharmacokinetic study. Similar to in adults, the atorvastatin (AVA) area-under-concentration-time curve from 0 to 24 h (AUC0-24) was 1.7-fold and 2.8-fold higher in participants with c.521T/C and c.521C/C compared to the c.521T/T participants, respectively. The inter-individual variability in AVA exposure within these genotype groups ranged from 2.3 to 4.8-fold, indicating that additional factors contribute to the inter-individual variability in the AVA dose-exposure relationship. A multivariate model reinforced the SLCO1B1 c.521T>C variant as the central factor contributing to AVA systemic exposure in this pediatric cohort, accounting for ~65% of the variability in AVA AUC0-24. Furthermore, lower AVA lactone concentrations in participants with increased body mass index contributed to higher exposure within the c.521T/T and c.521T/C genotype groups. Collectively, these factors contributing to higher systemic exposure could increase the risk of toxicity and should be accounted for when individualizing the dosing of atorvastatin in eligible pediatric patients.

Ticagrelor Increases Exposure to the Breast Cancer Resistance Protein Substrate Rosuvastatin

Ticagrelor and rosuvastatin are often used concomitantly after atherothrombotic events. Several cases of rhabdomyolysis during concomitant ticagrelor and rosuvastatin have been reported, suggesting a drug-drug interaction. We showed recently that ticagrelor inhibits breast cancer resistance protein (BCRP) and organic anion transporting polypeptide (OATP) 1B1, 1B3, and 2B1-mediated rosuvastatin transport in vitro. The aim of this study was to investigate the effects of ticagrelor on rosuvastatin pharmacokinetics in humans. In a randomized, crossover study, 9 healthy volunteers ingested a single dose of 90 mg ticagrelor or placebo, followed by a single 10 mg dose of rosuvastatin 1 hour later. Ticagrelor 90 mg or placebo were additionally administered 12, 24, and 36 hours after their first dose. Ticagrelor increased rosuvastatin area under the plasma concentration-time curve (AUC) and peak plasma concentration 2.6-fold (90% confidence intervals: 1.8-3.8 and 1.7-4.0, P = 0.001 and P = 0.003), and prolonged its half-life from 3.1 to 6.6 hours (P = 0.009). Ticagrelor also decreased the renal clearance of rosuvastatin by 11% (3%-19%, P = 0.032). The N-desmethylrosuvastatin:rosuvastatin AUC0-10h ratio remained unaffected by ticagrelor. Ticagrelor had no effect on the plasma concentrations of the endogenous OATP1B substrates glycodeoxycholate 3-O-glucuronide, glycochenodeoxycholate 3-O-glucuronide, glycodeoxycholate 3-O-sulfate, and glycochenodeoxycholate 3-O-sulfate, or the sodium-taurocholate cotransporting polypeptide substrate taurocholic acid. These data indicate that ticagrelor increases rosuvastatin concentrations more than twofold in humans, probably mainly by inhibiting intestinal BCRP. Because the risk for rosuvastatin-induced myotoxicity increases along with rosuvastatin plasma concentrations, using ticagrelor concomitantly with high doses of rosuvastatin should be avoided.

Influence of a Short Course of Ritonavir Used as Booster in Antiviral Therapies Against SARS-CoV-2 on the Exposure of Atorvastatin and Rosuvastatin

PURPOSE

Early antiviral treatment with nirmatrelvir/ritonavir is recommended for SARS-CoV-2-infected patients at high risk for severe courses. Such patients are usually chronically ill and susceptible to adverse drug interactions caused by ritonavir. We investigated the interactions of short-term low-dose ritonavir therapy with atorvastatin and rosuvastatin, two statins commonly used in this population.

METHOD

We assessed exposure changes (area under the concentration-time curve (AUC∞) and maximum concentration (Cmax)) of a single dose of 10 mg atorvastatin and 10 mg rosuvastatin before and on the fifth day of ritonavir treatment (2 × 100 mg/day) in healthy volunteers and developed a semi-mechanistic pharmacokinetic model to estimate dose adjustment of atorvastatin during ritonavir treatment.

RESULTS

By the fifth day of ritonavir treatment, the AUC∞ of atorvastatin increased 4.76-fold and Cmax 3.78-fold, and concurrently, the concentration of atorvastatin metabolites decreased to values below the lower limit of quantification. Pharmacokinetic modelling indicated that a stepwise reduction in atorvastatin dose during ritonavir treatment with a stepwise increase up to 4 days after ritonavir discontinuation can keep atorvastatin exposure within safe and effective margins. Rosuvastatin pharmacokinetics were only mildly modified; ritonavir significantly increased the Cmax 1.94-fold, while AUC∞ was unchanged.

CONCLUSION

Atorvastatin doses should likely be adjusted during nirmatrelvir/ritonavir treatment. For patients on a 20-mg dose, we recommend half of the original dose. In patients taking 40 mg or more, a quarter of the dose should be taken until 2 days after discontinuation of nirmatrelvir/ritonavir. Patients receiving rosuvastatin do not need to change their treatment regimen.

TRIAL REGISTRATION

EudraCT number: 2021-006634-39. DRKS00027838.

A gene risk score using missense variants in SLCO1B1 is associated with earlier onset statin intolerance

BACKGROUND AND AIMS

The efficacy of statin therapy is hindered by intolerance to the therapy, leading to discontinuation. Variants in SLCO1B1, which encodes the hepatic transporter OATB1B1, influence statin pharmacokinetics, resulting in altered plasma concentrations of the drug and its metabolites. Current pharmacogenetic guidelines require sequencing of the SLCO1B1 gene, which is more expensive and less accessible than genotyping. In this study, we aimed to develop an easy, clinically implementable functional gene risk score (GRS) of common variants in SLCO1B1 to identify patients at risk of statin intolerance.

METHODS AND RESULTS

A GRS was developed from four common variants in SLCO1B1. In statin users from Tayside, Scotland, UK, those with a high-risk GRS had increased odds across three phenotypes of statin intolerance [general statin intolerance (GSI): ORGSI 2.42; 95% confidence interval (CI): 1.29-4.31, P = 0.003; statin-related myopathy: ORSRM 2.51; 95% CI: 1.28-4.53, P = 0.004; statin-related suspected rhabdomyolysis: ORSRSR 2.85; 95% CI: 1.03-6.65, P = 0.02]. In contrast, using the Val174Ala genotype alone or the recommended OATP1B1 functional phenotypes produced weaker and less reliable results. A meta-analysis with results from adjudicated cases of statin-induced myopathy in the PREDICTION-ADR Consortium confirmed these findings (ORVal174Ala 1.99; 95% CI: 1.01-3.95, P = 0.048; ORGRS 1.76; 95% CI: 1.16-2.69, P = 0.008). For those requiring high-dose statin therapy, the high-risk GRS was more consistently associated with the time to onset of statin intolerance amongst the three phenotypes compared with Val174Ala (GSI: HRVal174Ala 2.49; 95% CI: 1.09-5.68, P = 0.03; HRGRS 2.44; 95% CI: 1.46-4.08, P < 0.001). Finally, sequence kernel association testing confirmed that rare variants in SLCO1B1 are associated with the risk of intolerance (P = 0.02).

CONCLUSION

We provide evidence that a GRS based on four common SLCO1B1 variants provides an easily implemented genetic tool that is more reliable than the current recommended practice in estimating the risk and predicting early-onset statin intolerance.

Pharmacogenetic analysis of inter-ethnic variability in the uptake transporter SLCO1B1 gene in Colombian, Mozambican, and Portuguese populations

BACKGROUND

Statin-induced myopathy is reported to be associated with the solute carrier organic anion transporter family member 1B1 gene single nucleotide polymorphism, c.521 T > C. There is no epidemiologic data on this gene polymorphism in several countries. Therefore, this study aimed at assessing the genotype and allele frequencies of the gene variant in three countries.

METHODS

This study involved healthy individuals from Colombia, Mozambique, and Portugal. Genomic DNA was isolated from blood samples using the Qiamp DNA Extraction Kit (Qiagen). The isolated DNA was genotyped using novel Polymerase Chain Reaction-Restriction Fragment Length Polymorphism. Microstat and GraphPad QuickCal software were used for the Chi-square test and the evaluation of Hardy-Weinberg equilibrium respectively.

RESULTS

A total of 181 individuals' blood samples were analyzed. Overall, the TT (74.0%) genotype was the highest and the CC (7.8%) was the lowest. Country wise genotypic frequencies were Colombia 47(70.2%) TT, 12(17.9%) TC and 8(11.9%) CC; Mozambique 47(88.7%) TT, 5(9.4%) TC, and 1(1.9%) CC; and Portugal 40(65.6%) TT, 16(26.2%) TC, and 5(8.2%) CC. The reference (T) allele was highest among Mozambicans (93.4%) compared to Colombians (79.1%) and Portuguese (78.7%). Mozambicans showed statistically significant genotypic and allelic frequency differences compared to Colombians (p < 0.01) and Portuguese (p < 0.01).

CONCLUSIONS

Overall and country-wise, CC genotype was less frequent and it is relatively high for Colombians and Portuguese populations. This finding may imply statins risk-benefit variability associated with CC genotype among these populations that needs further understanding.

SLCO1B1*5 is protective against non-senile cataracts in cohort prescribed statins: analysis in a British-South Asian cohort

BACKGROUND

Reported association between statin use and cataract risk is controversial. The SLCO1B1 gene encodes a transport protein responsible for statin clearance. The aim of this study was to investigate a possible association between the SLCO1B1*5 reduced function variant and cataract risk in statin users of South Asian ethnicity.

METHODS

The Genes & Health cohort consists of British-Bangladeshi and British-Pakistani participants from East London, Manchester and Bradford, UK. SLCO1B1*5 genotype was assessed with the Illumina GSAMD-24v3-0-EA chip. Medication data from primary care health record linkage was used to compare those who had regularly used statins compared to those who had not. Multivariable logistic regression was used to test for association between statin use and cataracts, adjusting for population characteristics and potential confounders in 36,513 participants. Multivariable logistic regression was used to test association between SLCO1B1*5 heterozygotes or homozygotes and cataracts, in subgroups having been regularly prescribed statins versus not.

RESULTS

Statins were prescribed to 35% (12,704) of participants (average age 41 years old, 45% male). Non-senile cataract was diagnosed in 5% (1686) of participants. An apparent association between statins and non-senile cataract (12% in statin users and 0.8% in non-statin users) was negated by inclusion of confounders. In those prescribed a statin, presence of the SLCO1B1*5 genotype was independently associated with a decreased risk of non-senile cataract (OR 0.7 (CI 0.5-0.9, p 0.007)).

CONCLUSIONS

Our findings suggest that there is no independent association between statin use and non-senile cataract risk after adjusting for confounders. Among statin users, the SLCO1B1*5 genotype is associated with a 30% risk reduction of non-senile cataracts. Stratification of on-drug cohorts by validated pharmacogenomic variants is a useful tool to support or repudiate adverse drug events in observational cohorts.

Real-world pharmacogenetics of statin intolerance: effects of SLCO1B1, ABCG2 , and CYP2C9 variants

OBJECTIVE

The association of SLCO1B1 c.521T>C with simvastatin-induced muscle toxicity is well characterized. However, different statins are subject to metabolism and transport also by other proteins exhibiting clinically meaningful genetic variation. Our aim was to investigate associations of SLCO1B1 c.521T>C with intolerance to atorvastatin, fluvastatin, pravastatin, rosuvastatin, or simvastatin, those of ABCG2 c.421C>A with intolerance to atorvastatin, fluvastatin, or rosuvastatin, and that of CYP2C9*2 and *3 alleles with intolerance to fluvastatin.

METHODS

We studied the associations of these variants with statin intolerance in 2042 patients initiating statin therapy by combining genetic data from samples from the Helsinki Biobank to clinical chemistry and statin purchase data.

RESULTS

We confirmed the association of SLCO1B1 c.521C/C genotype with simvastatin intolerance both by using phenotype of switching initial statin to another as a marker of statin intolerance [hazard ratio (HR) 1.88, 95% confidence interval (CI) 1.08-3.25, P  = 0.025] and statin switching along with creatine kinase measurement (HR 5.44, 95% CI 1.49-19.9, P  = 0.011). No significant association was observed with atorvastatin and rosuvastatin. The sample sizes for fluvastatin and pravastatin were relatively small, but SLCO1B1 c.521T>C carriers had an increased risk of pravastatin intolerance defined by statin switching when compared to homozygous reference T/T genotype (HR 2.11, 95% CI 1.01-4.39, P  = 0.047).

CONCLUSION

The current results can inform pharmacogenetic statin prescribing guidelines and show feasibility for the methodology to be used in larger future studies.

Effect of SLCO1B1 c.521T>C polymorphism on the lipid response to statins in people living with HIV on a boosted protease inhibitor-containing regimen

AIMS

We previously observed that some individuals on HIV boosted protease inhibitor-containing regimen do not achieve their lipid targets despite elevated statin concentrations. This study evaluated whether the common single polymorphism c.521T>C in SLCO1B1, associated with reduced statin uptake in the liver, could explain this observation.

METHODS

People living with HIV in the Swiss HIV Cohort Study were eligible if they were on a boosted protease inhibitor concomitantly with a statin for at least 6 months and if their SLCO1B1 genotype was available. Furthermore, their lipids had to be documented before and after the introduction of the statin. The statin efficacy was defined as % change in total cholesterol, low-density lipoprotein-cholesterol, high-density lipoprotein-cholesterol and triglycerides levels after statin initiation compared to pretreatment levels. Lipid response was adjusted for differences in potency and dose between statins.

RESULTS

In total, 88 people living with HIV were included, of whom 58, 28 and 2 carried the SLCO1B1 TT, TC and CC genotypes, respectively. The change in lipid levels after statin initiation tended to be lower in carriers of the polymorphism although the difference was not statistically significant (TT vs. TC/CC: total cholesterol: -11.7 vs. -4.8%; low-density lipoprotein- cholesterol: -20.6 vs. -7.4%; high-density lipoprotein-cholesterol: 1.6 vs. 0%; triglycerides: -11.5 vs. -7.9%). In the multiple linear regression, change in total cholesterol was inversely correlated with the total cholesterol level prestatin treatment (coefficient -6.60, 95% confidence interval: -9.63 to -3.56, P < .001).

CONCLUSION

The lipid-lowering effect of statins tended to be attenuated by SLCO1B1 polymorphism and progressively declined as total cholesterol under the boosted protease inhibitor treatment decreased.

Rhabdomyolysis during concomitant ticagrelor and rosuvastatin: A breast cancer resistance protein-mediated drug interaction?

We present 3 patients diagnosed with rhabdomyolysis 1-6 months after the initiation of concomitant rosuvastatin and ticagrelor medication. A literature review and Food and Drug Administration adverse event reporting system revealed >40 reports of rhabdomyolysis during concomitant ticagrelor and rosuvastatin, including 3 with a fatal outcome. We show that ticagrelor inhibits breast cancer resistance protein-, organic anion transporting polypeptide (OATP) 1B1-, 1B3- and 2B1-mediated transport of rosuvastatin in vitro with half-maximal unbound inhibitory concentrations of 0.36, 4.13, 7.5 and 3.26 μM, respectively. A static drug interaction model predicted that ticagrelor may inhibit intestinal breast cancer resistance protein and thus increase rosuvastatin plasma exposure 2.1-fold, whereas the OATP-mediated hepatic uptake of rosuvastatin should not be inhibited due to relatively low portal ticagrelor concentrations. Taken together, concomitant use of ticagrelor with rosuvastatin may increase the systemic exposure to rosuvastatin and the risk of rosuvastatin-induced rhabdomyolysis. Further studies are warranted to investigate the potential pharmacokinetic interaction between ticagrelor and rosuvastatin in humans.

The frequency of major ABCG2, SLCO1B1 and CYP2C9 variants in Asian, Native Hawaiian and Pacific Islander women subgroups: implications for personalized statins dosing

Aim: The frequencies of SLCO1B1*5 and CYP2C9*2 and *3 in specific Asian, Native Hawaiian and Pacific Islander (NHPI) subgroups are unknown. Patients & methods: Repository DNA samples from 1064 women self-identifying as Filipino, Korean, Japanese, Native Hawaiian, Marshallese or Samoan and aged 18 years or older were used for targeted sequencing of three genetic variants (rs4149056, rs1799853 and rs1057910). Results: SLCO1B1*5 was significantly less frequent in NHPI women (0.5-6%) than in Europeans (16%). Except for Koreans, CYP2C9*2 (0-1.4%) and *3 (0.5-3%) were significantly less frequent in all subgroups than in Europeans (8 and 12.7%, respectively). Prior reports showed that Asian and NHPI individuals have significantly higher ABCG2 Q141K allele frequency (13-46%) than Europeans (9.4%). Combined phenotype rates for rosuvastatin and fluvastatin revealed that Filipinos and Koreans had the highest frequencies of statin-associated myopathy symptoms risk alleles. Conclusion: Differences in ABCG2, SLCO1B1 and CYP2C9 allele frequencies among different racial and ethnic subgroups highlight the need for increased diversity in pharmacogenetic research. Risk alleles for statin-associated myopathy symptoms are more prevalent in Filipinos, underscoring the importance of genotype-based statin dosing.

Atorvastatin Metabolite Pattern in Skeletal Muscle and Blood from Patients with Coronary Heart Disease and Statin-Associated Muscle Symptoms

Self-perceived statin-associated muscle symptoms (SAMS) are prevalent, but only a minority is drug-dependent. Diagnostic biomarkers are not yet identified. The local statin exposure in skeletal muscle tissue may correlate to the adverse effects. We aimed to determine whether atorvastatin metabolites in blood reflect the corresponding metabolite levels in skeletal muscle, and whether genetic variants of statin transporters modulate this relationship. We also addressed atorvastatin metabolites as potential objective biomarkers of SAMS. Muscle symptoms were examined in patients with coronary disease and self-perceived SAMS during 7 weeks of double-blinded treatment with atorvastatin 40 mg/day and placebo in randomized order. A subset of 12 patients individually identified with more muscle symptoms on atorvastatin than placebo (confirmed SAMS) and 15 patients with no difference in muscle symptom intensity (non-SAMS) attended the present follow-up study. All received 7 weeks of treatment with atorvastatin 40 mg/day followed by 8 weeks without statins. Biopsies from the quadriceps muscle and blood plasma were collected after each treatment period. Strong correlations (rho > 0.7) between muscle and blood plasma concentrations were found for most atorvastatin metabolites. The impact of the SLCO1B1 c.521T>C (rs4149056) gene variant on atorvastatin's systemic pharmacokinetics was translated into muscle tissue. The SLCO2B1 c.395G>A (rs12422149) variant did not modulate the accumulation of atorvastatin metabolites in muscle tissue. Atorvastatin pharmacokinetics in patients with confirmed SAMS were not different from patients with non-SAMS. In conclusion, atorvastatin metabolite levels in skeletal muscle and plasma are strongly correlated, implying that plasma measurements are suitable proxies of atorvastatin exposure in muscle tissue. The relationship between atorvastatin metabolites in plasma and SAMS deserves further investigation.

PharmVar GeneFocus: SLCO1B1

The Pharmacogene Variation Consortium (PharmVar) is now providing star (*) allele nomenclature for the highly polymorphic human SLCO1B1 gene encoding the organic anion transporting polypeptide 1B1 (OATP1B1) drug transporter. Genetic variation within the SLCO1B1 gene locus impacts drug transport, which can lead to altered pharmacokinetic profiles of several commonly prescribed drugs. Variable OATP1B1 function is of particular importance regarding hepatic uptake of statins and the risk of statin-associated musculoskeletal symptoms. To introduce this important drug transporter gene into the PharmVar database and serve as a unified reference of haplotype variation moving forward, an international group of gene experts has performed an extensive review of all published SLCO1B1 star alleles. Previously published star alleles were self-assigned by authors and only loosely followed the star nomenclature system that was first developed for cytochrome P450 genes. This nomenclature system has been standardized by PharmVar and is now applied to other important pharmacogenes such as SLCO1B1. In addition, data from the 1000 Genomes Project and investigator-submitted data were utilized to confirm existing haplotypes, fill knowledge gaps, and/or define novel star alleles. The PharmVar-developed SLCO1B1 nomenclature has been incorporated by the Clinical Pharmacogenetics Implementation Consortium (CPIC) 2022 guideline on statin-associated musculoskeletal symptoms.

The frequency of rs2231142 in ABCG2 among Asian subgroups: implications for personalized rosuvastatin dosing

Statins are widely used medications for the primary and secondary prevention of cardiovascular diseases. Statin-induced musculoskeletal symptoms are the primary adverse drug events contributing to poor adherence to lipid-lowering therapy. Rosuvastatin is characterized by interindividual differences in systemic exposure among different patient population subgroups. The missense variant Q141K within ABCG2, highly prevalent in some Asian subgroups, results in decreased transporter efflux function and increased exposure to rosuvastatin. We aim to highlight the implications of ABCG2 genotype in prescribing rosuvastatin and the ramifications of interpopulation differences in Q141K frequencies in the starting dose of rosuvastatin in major Asian subgroups, using the most recent genetic-based guidelines. The high frequency of Q141K in Filipinos could warrant a lower starting rosuvastatin dose versus non-Filipinos. The Q141K genotype frequencies in Asian subgroups suggest significant interpopulation differences, reinforcing the need to move beyond race-based to genotype-based rosuvastatin dosing.

A comprehensive pharmacogenomic study indicates roles for SLCO1B1, ABCG2 and SLCO2B1 in rosuvastatin pharmacokinetics

AIMS

The aim was to comprehensively investigate the effects of genetic variability on the pharmacokinetics of rosuvastatin.

METHODS

We conducted a genome-wide association study and candidate gene analyses of single dose rosuvastatin pharmacokinetics in a prospective study (n = 159) and a cohort of previously published studies (n = 88).

RESULTS

In a genome-wide association meta-analysis of the prospective study and the cohort of previously published studies, the SLCO1B1 c.521 T > C (rs4149056) single nucleotide variation (SNV) associated with increased area under the plasma concentration-time curve (AUC) and peak plasma concentration of rosuvastatin (P = 1.8 × 10^-12 and P = 3.2 × 10^-15 ). The candidate gene analysis suggested that the ABCG2 c.421C > A (rs2231142) SNV associates with increased rosuvastatin AUC (P = .0079), while the SLCO1B1 c.388A > G (rs2306283) and SLCO2B1 c.1457C > T (rs2306168) SNVs associate with decreased rosuvastatin AUC (P = .0041 and P = .0076). Based on SLCO1B1 genotypes, we stratified the participants into poor, decreased, normal, increased and highly increased organic anion transporting polypeptide (OATP) 1B1 function groups. The OATP1B1 poor function phenotype associated with 2.1-fold (90% confidence interval 1.6-2.8, P = 4.69 × 10^-5 ) increased AUC of rosuvastatin, whereas the OATP1B1 highly increased function phenotype associated with a 44% (16-62%; P = .019) decreased rosuvastatin AUC. The ABCG2 c.421A/A genotype associated with 2.2-fold (1.5-3.0; P = 2.6 × 10^-4 ) increased AUC of rosuvastatin. The SLCO2B1 c.1457C/T genotype associated with 28% decreased rosuvastatin AUC (11-42%; P = .01).

CONCLUSION

These data suggest roles for SLCO1B1, ABCG2 and SLCO2B1 in rosuvastatin pharmacokinetics. Poor SLCO1B1 or ABCG2 function genotypes may increase the risk of rosuvastatin-induced myotoxicity. Reduced doses of rosuvastatin are advisable for patients with these genotypes.

Pharmacogenomics of lipid-lowering agents: the impact on efficacy and safety

Hyperlipidemia is a significant risk factor for cardiovascular disease morbidity and mortality. The lipid-lowering drugs are considered the cornerstone of primary and secondary prevention of atherosclerotic cardiovascular disease. Unfortunately, the lack of efficacy and associated adverse effects, ranging from mild-to-moderate to potentially life-threatening, lead to therapy discontinuation. Numerous reports support the role of gene polymorphisms in drugs' pharmacokinetic parameters and their associated adverse reactions. Therefore, this study aims to understand the pharmacogenomics of lipid-lowering drugs and the impact of genetic variants of key genes on the drugs' efficacy and toxicity. Indeed, genetically guided lipid-lowering therapy enhances overall safety, improves drug adherence and achieves long-term therapy.

The Farnesoid X Receptor as a Master Regulator of Hepatotoxicity

The nuclear receptor farnesoid X receptor (FXR, NR1H4) is a bile acid (BA) sensor that links the enterohepatic circuit that regulates BA metabolism and elimination to systemic lipid homeostasis. Furthermore, FXR represents a real guardian of the hepatic function, preserving, in a multifactorial fashion, the integrity and function of hepatocytes from chronic and acute insults. This review summarizes how FXR modulates the expression of pathway-specific as well as polyspecific transporters and enzymes, thereby acting at the interface of BA, lipid and drug metabolism, and influencing the onset and progression of hepatotoxicity of varying etiopathogeneses. Furthermore, this review article provides an overview of the advances and the clinical development of FXR agonists in the treatment of liver diseases.

Sinusoidal Organic Anion-Transporting Polypeptide 1B1/1B3 and Bile Canalicular Multidrug Resistance-Associated Protein 2 Play an Essential Role in the Hepatobiliary Disposition of a Synthetic Cyclic Dinucleotide (STING Agonist)

The liver is central to the elimination of many drugs from the body involving multiple processes and understanding of these processes is important to quantitively assess hepatic clearance of drugs. The synthetic STING (STimulator of INterferon Genes protein) agonist is a new class of drugs currently being evaluated in clinical trials as a potential anticancer therapy. In this study, we used ML00960317 (synthetic STING agonist) to investigate the hepatobiliary disposition of this novel molecular entity. A bile-duct cannulated (BDC) rat study indicated that biliary excretion is the major route of elimination for ML00960317 (84% of parent dose in bile). The human biliary clearance using in vitro sandwich cultured human hepatocyte model predicted significant biliary excretion of ML00960317 (biliary excretion index (BEI) of 47%). Moreover, the transport studies using transporter expressing cell lines, hepatocytes, and membrane vesicles indicated that ML00960317 is a robust substrate of OATP1B1, OATP1B3, and MRP2. Using relative expression factor approach, the combined contribution of OATP1B1 (fraction transported (ft) = 0.62) and OATP1B3 (ft = 0.31) was found to be 93% of the active uptake clearance of ML00960317 into the liver. Furthermore, OATP1B1 and OATP1B3-mediated uptake of ML00960317 was inhibited by rifampicin with IC50 of 6.5 and 2.3 μM, respectively indicating an in vivo DDI risk (R value of 1.5 and 2.5 for OATP1B1 and OATP1B3, respectively). These results highlighted an important role of OATP1B1, OATP1B3, and MRP2 in the hepatobiliary disposition of ML00960317. These pathways may act as rate-determining steps in the hepatic clearance of ML00960317 thus presenting clinical DDI risk.

Drug dosing in children with obesity: a narrative updated review

Childhood obesity and its associated comorbidities are highly prevalent diseases that may add to any other possible health problem commonly affecting the pediatric age. Uncertainties may arise concerning drug dosing when children with obesity need pharmacologic therapies. In general, in pediatric practice, there is a tendency to adapt drug doses to a child's total body weight. However, this method does not consider the pharmacological impact that a specific drug can have under a two-fold point of view, that is, across various age and size groups as well. Moreover, there is a need for a therapeutic approach, as much as possible tailored considering relevant interacting aspects, such as modification in metabolomic profile, drug pharmacokinetics and pharmacodynamics. Taking into account the peculiar differences between children with overweight/obesity and those who are normal weight, the drug dosage in the case of obesity, cannot be empirically determined solely by the per kg criterion. In this narrative review, we examine the pros and cons of several drug dosing methods used when dealing with children who are affected also by obesity, focusing on specific aspects of some of the drugs most frequently prescribed in real-world practice by general pediatricians and pediatric subspecialists.

Frequency of functional exonic single-nucleotide polymorphisms and haplotype distribution in the SLCO1B1 gene across genetic ancestry groups in the Qatari population

Organic anion transporting polypeptides (OATP), which are encoded by SLCO genes, participate in the hepatic elimination of drugs and xenobiotics. SLCO1B1 is an important pharmacogenomic gene (encoding OATP1B1) associated with response to the uptake of endogenous compounds, such as statin and bilirubin. Ethnicity of the patient modulates the response to these drugs; the frequency and haplotype data for SLCO1B1 genetic variants in the Arab population is lacking. Therefore, we determined the frequencies of two well-characterized SLCO1B1 single nucleotide polymorphisms (SNP) and haplotypes that affect the OATP1B1 drugs transportation activity in Qatari population. Genotyping data for two SLCO1B1 SNPs (c.388A > G, c.521 T > C) were extracted from whole exome data of 1050 Qatari individuals, who were divided into three ancestry groups, namely Bedouins, Persians/South Asians, and Africans. By way of using Fisher's exact and Chi-square tests, we evaluated the differences in minor allele frequency (MAF) of the two functional SNPs and haplotype frequencies (HF) among the three ancestry groups. The OATP1B1 phenotypes were assigned according to their function by following the guidelines from the Clinical Pharmacogenetics Implementation Consortium for SLCO1B1 and Simvastatin-Induced Myopathy.The MAF of SLCO1B1:c.388A > G was higher compared to that of SLCO1B1:c.521 T > C in the study cohort. It was significantly high in the African ancestry group compared with the other two groups, whereas SLCO1B1:c.521 T > C was significantly low in the African ancestry group compared with the other two groups. The SLCO1B1 *15 haplotype had the highest HF, followed by *1b, *1a, and *5. Only the SLCO1B1 *5 haplotype showed no significant difference in frequency across the three ancestry groups. Furthermore, we observed that the OATP1B1 normal function phenotype accounted for 58% of the Qatari individuals, the intermediate function phenotype accounted for 35% with significant differences across the ancestry groups, and the low function phenotype accounted for 6% of the total Qatari individuals with a higher trend observed in the Bedouin group.The results indicate that the phenotype frequencies of the OATP1B1 intermediate and low function in the Qatari population appear at the higher end of the frequency range seen worldwide. Thus, a pharmacogenetic screening program for SLCO1B1 variants may be necessary for the Qatari population.

Novel Insights in Drug Transporter Sciences: the Year 2021 in Review

On behalf of the team I am pleased to present the second annual 'novel insights into drug transporter sciences review' focused on peer-reviewed articles that were published in the year 2021. In compiling the articles for inclusion, preprints available in 2021 but officially published in 2022 were considered to be in scope. To support this review the contributing authors independently selected one or two articles that were thought to be impactful and of interest to the broader research community. A similar approach as published last year was adopted whereby key observations, methods and analysis of each paper is concisely summarized in the synopsis followed by a commentary highlighting the impact of the paper in understanding of drug transporters' role in drug disposition.As the goal of this review is not to provide a comprehensive overview of each paper but rather highlight important findings that are well supported by the data, the reader is encouraged to consult the original articles for additional information. Further, and keeping in line with the goals of this review, it should be noted that all authors actively contributed by writing synopsis and commentary for individual papers and no attempt was made to standardize language or writing styles. In this way, the review article is reflective of not only the diversity of the articles but also that of the contributors. I extend my thanks to the authors for their continued support and also welcome Diane Ramsden and Pallabi Mitra as contributing authors for this issue.

Functional in vitro characterization of SLCO1B1 variants and simulation of the clinical pharmacokinetic impact of impaired OATP1B1 function

Organic Anion Transporting Polypeptide 1B1 is important to the hepatic elimination and distribution of many drugs. If OATP1B1 function is decreased, it can increase plasma exposure of e.g. several statins leading to increased risk of muscle toxicity. First, we examined the impact of three naturally occurring rare variants and the frequent SLCO1B1 c.388A>G variant on in vitro transport activity with cellular uptake assay using two substrates: 2', 7'-dichlorofluorescein (DCF) and rosuvastatin. Secondly, LC-MS/MS based quantitative targeted absolute proteomics measured the OATP1B1 protein abundance in crude membrane fractions of HEK293 cells over-expressing these SNVs. Additionally, we simulated the effect of impaired OATP1B1 function on rosuvastatin pharmacokinetics to estimate the need for genotype-guided dosing. R57Q impaired DCF and rosuvastatin transport significantly yet did not change protein expression considerably, while N130D and N151S did not alter activity but increased protein expression. R253Q did not change protein expression but reduced DCF uptake and increased rosuvastatin K_m. Based on pharmacokinetics simulations, doses of 30 mg (with 50% OATP1B1 function) and 20 mg (with 0% OATP1B1 function) result in plasma exposure similar to 40 mg dose (with 100% OATP1B1 function). Therefore dose reductions might be considered to avoid increased plasma exposure caused by function-impairing OATP1B1 genetic variants, such as R57Q.

Pharmacogenomic Study of Statin-Associated Muscle Symptoms in the ODYSSEY OUTCOMES Trial

Statin-associated muscle symptoms (SAMS) are the most frequently reported adverse events for statin therapies. Previous studies have reported an association between the p.Val174Ala missense variant in SLCO1B1 and SAMS in simvastatin-treated subjects; however, evidence for genetic predictors of SAMS in atorvastatin- or rosuvastatin-treated subjects is currently lacking.

Atorvastatin population pharmacokinetics in a real-life setting: Influence of genetic polymorphisms and association with clinical response

The purpose of this study was to investigate the potential clinical relevance of estimating the apparent clearance (CL/F) of atorvastatin through population pharmacokinetic (PopPK) modeling with samples collected in a real-life setting in a cohort of ambulatory patients at risk of cardiovascular disease by using an opportunistic sampling strategy easily accessible in clinical routine. A total of 132 pharmacokinetic (PK) samples at a maximum of three visits were collected in the 70 included patients. The effects of demographic, genetic, and clinical covariates were also considered. With the collected data, we developed a two-compartment PopPK model that allowed estimating atorvastatin CL/F relatively precisely and considering the genotype of the patient for SLCO1B1 c.521T>C single-nucleotide polymorphism (SNP). Our results indicate that the estimation of the CL/F of atorvastatin through our PopPK model might help in identifying patients at risk of myalgia. Indeed, we showed that a patient presenting a CL/F lower than 414.67 L h-1 is at risk of suffering from muscle discomfort. We also observed that the CL/F was correlated with the efficacy outcomes, suggesting that a higher CL/F is associated with a better drug efficacy (i.e., a greater decrease in total and LDL-cholesterol levels). In conclusion, our study demonstrates that PopPK modeling can be useful in daily clinics to estimate a patient' atorvastatin clearance. Notifying the clinician with this information can help in identifying patients at risk of myalgia and gives indication about the potential responsiveness to atorvastatin therapy.

Ten Years of Experience Support Pharmacogenetic Testing to Guide Individualized Drug Therapy

Precision medicine utilizing the genetic information of genes involved in the metabolism and disposition of drugs can not only improve drug efficacy but also prevent or minimize adverse events. Polypharmacy is common among multimorbid patients and is associated with increased adverse events. One of the main objectives in health care is safe and efficacious drug therapy, which is directly correlated to the individual response to treatment. Precision medicine can increase drug safety in many scenarios, including polypharmacy. In this report, we share our experience utilizing precision medicine over the past ten years. Based on our experience using pharmacogenetic (PGx)-informed prescribing, we implemented a five-step precision medicine protocol (5SPM) that includes the assessment of the biological-clinical characteristics of the patient, current and past prescription history, and the patient's PGx test results. To illustrate our approach, we present cases highlighting the clinical relevance of precision medicine with a focus on patients with a complex history and polypharmacy.

A phase I/IIa trial of atorvastatin in Japanese patients with acute Kawasaki disease with coronary artery aneurysm: Study protocol of a multicenter, single-arm, open-label trial

Background

Kawasaki disease (KD) is a systemic vasculitis complicated with coronary artery abnormalities (CAAs). Intravenous immunoglobulin reduces the occurrence of CAAs, but significant number of KD patients with CAAs still exists. Thus, new approaches to prevent and attenuate CAAs are warranted. Atorvastatin has been shown to promote endothelial cell homeostasis and suppress vascular inflammation and has received enthusiasm as a potentially new candidate treatment for KD. In the United States, a phase I/IIa dose-escalation study of atorvastatin in KD patients with CAAs demonstrated the safety and pharmacokinetic data of atorvastatin. However, due to the uncertainty in the application of these results to other populations, we aim to examine the tolerability and generate pharmacokinetics data in Japanese KD patients.

Methods

This is a multicenter, single-arm, open-label, phase I/IIa study of atorvastatin in acute KD patients with CAAs in Japan. A minimum of 9 and a maximum of 18 KD patients (2 years–17 years old) will be recruited for a 3 + 3 dose-escalation study of a 6-week course of atorvastatin (0.125–0.5 mg/kg/day). The primary outcome will be safety of atorvastatin. The secondary outcomes will be pharmacokinetics of atorvastatin, activity of atorvastatin and echocardiographic assessment of CAAs. The activity of atorvastatin will include assessment of C-reactive protein or high sensitivity C-reactive protein and white blood cell levels.

Discussion

This study will provide evidence of the safety, tolerability, and pharmacokinetics of atorvastatin in Japanese KD patients and may lead new standard therapy for acute-phase KD associated with CAA complications.

Trial registration

Japan Registry of Clinical Trials (JRCTs031180057). Registered December 19, 2018, https://jrct.niph.go.jp/en-latest-detail/jRCTs031180057.

Serum Predose Metabolic Profiling for Prediction of Rosuvastatin Pharmacokinetic Parameters in Healthy Volunteers

Rosuvastatin is a well-known lipid-lowering agent generally used for hypercholesterolemia treatment and coronary artery disease prevention. There is a substantial inter-individual variability in the absorption of statins usually caused by genetic polymorphisms leading to a variation in the corresponding pharmacokinetic parameters, which may affect drug therapy safety and efficacy. Therefore, the investigation of metabolic markers associated with rosuvastatin inter-individual variability is exceedingly relevant for drug therapy optimization and minimizing side effects. This work describes the application of pharmacometabolomic strategies using liquid chromatography coupled to mass spectrometry to investigate endogenous plasma metabolites capable of predicting pharmacokinetic parameters in predose samples. First, a targeted method for the determination of plasma concentration levels of rosuvastatin was validated and applied to obtain the pharmacokinetic parameters from 40 enrolled individuals; then, predose samples were analyzed using a metabolomic approach to search for associations between endogenous metabolites and the corresponding pharmacokinetic parameters. Data processing using machine learning revealed some candidates including sterols and bile acids, carboxylated metabolites, and lipids, suggesting the approach herein described as promising for personalized drug therapy.

Determining the potential clinical value of panel-based pharmacogenetic testing in patients with chronic pain or gastroesophageal reflux disease

We aimed to determine the potential value of panel-based pharmacogenetic (PGx) testing in patients with chronic pain or gastroesophageal reflux disease (GERD) who underwent single-gene PGx testing to guide opioid or proton pump inhibitor (PPI) therapy, respectively. Of 448 patients included (chronic pain, n = 337; GERD, n = 111), mean age was 57 years, 68% were female, and 73% were white. Excluding opiates for the pain cohort and PPIs for the GERD cohort, 76.6% of patients with pain and 71.2% with GERD were prescribed at least one additional medication with a high level of PGx evidence, most commonly ondansetron or selective serotonin reuptake inhibitors. The most common genes that could inform PGx drug prescribing were CYP2C19, CYP2D6, CYP2C9, and SLCO1B1. Our findings suggest that patients with chronic pain or GERD are commonly prescribed drugs with a high level of evidence for a PGx-guided approach, supporting panel-based testing in these populations.

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.

Loss of function polymorphisms in SLCO1B1 (c.521T>C, rs4149056) and ABCG2 (c.421C>A, rs2231142) genes are associated with adverse events of rosuvastatin: a case-control study

PURPOSE

The study aims to evaluate relationship between polymorphisms associated with a reduced function of two transporter proteins resulting in increased exposure to rosuvastatin - organic anion transporter 1B1 (OATP1B1) (SLCO1B1 c.521T>C) and ATP binding cassette subfamily G member 2 (ABCG2) (ABCG2 c.421C>A) and occurrence of rosuvastatin related myotoxicity/hepatotoxicity.

METHODS

In a case-control study, cases (rosuvastatin treated patients developing myotoxicity or hepatotoxicity) and controls (concurrent rosuvastatin treated patients free of adverse events) were prospectively recruited over a 2 year period in a single tertiary center specialized in treatment of metabolic disorders. Subjects were evaluated for clinical, comorbidity, and comedication characteristics and for genotype predicted metabolizing phenotypes regarding cytochrome P450 enzymes CYP2C9 and CYP2C19. Standard regression analysis and analysis in matched sets of cases and controls (optimal full matching) were undertaken by fitting frequentist and Bayesian models (covariates/matching variables: age, sex, diabetes, liver/renal disease, hypertension, CYP2C9 and C19 phenotype, use of CYP or transporter inhibitors, non evaluated transporter genotype).

RESULTS

A total of 88 cases (81 with myotoxicity, 6 with hepatotoxicity, 1 with both) and 129 controls were recruited. Odds of variant SLCO1B1 c.521T>C allele were 2.2-2.5 times higher in cases than in controls (OR = 2.45, 95% CI 1.34-4.48; Bayesian OR = 2.59, 95% CrI 1.42-4.90 in regression analysis; OR = 2.20, 1.10-4.42; Bayesian OR = 2.26, 1.28-4.41 in matched analysis). Odds of variant ABCG2 c.421C>A allele were 2.1-2.3 times higher in cases than in controls (OR = 2.24, 1.04-4.83; Bayesian OR = 2.35, 1.09-4.31 in regression analysis; OR = 2.10, 0.83-5.31; Bayesian OR = 2.17, 1.07-4.35 in matched analysis).

CONCLUSION

Loss of function polymorphisms in SLCO1B1 c.521T>C and ABCG2 c.421C>A genes are associated with the presence of rosuvastatin related myotoxicity and/or hepatotoxicity.

The Effect of Single Nucleotide Variations in the Transmembrane Domain of OATP1B1 on in vitro Functionality

Purpose

Organic Anion Transporting Polypeptide 1B1 (OATP1B1) mediates hepatic influx and clearance of many drugs, including statins. The SLCO1B1 gene is highly polymorphic and its function-impairing variants can predispose patients to adverse effects. The effects of rare genetic variants of SLCO1B1 are mainly unexplored. We examined the impact of eight naturally occurring rare variants and the well-known SLCO1B1 c.521C > T (V174A) variant on in vitro transport activity, cellular localization and abundance.

Methods

Transport of rosuvastatin and 2,7-dichlorofluorescein (DCF) in OATP1B1 expressing HEK293 cells was measured to assess changes in activity of the variants. Immunofluorescence and confocal microscopy determined the cellular localization of OATP1B1 and LC–MS/MS based quantitative targeted absolute proteomics analysis quantified the amount of OATP1B1 in crude membrane fractions.

Results

All studied variants, with the exception of P336R, reduced protein abundance to varying degree. V174A reduced protein abundance the most, over 90% compared to wild type. Transport function was lost in G76E, V174A, L193R and R580Q variants. R181C decreased activity significantly, while T345M and L543W retained most of wild type OATP1B1 activity. P336R showed increased activity and H575L decreased the transport of DCF significantly, but not of rosuvastatin. Decreased activity was interrelated with lower absolute protein abundance in the studied variants.

Conclusions

Transmembrane helices 2, 4 and 11 appear to be crucial for proper membrane localization and function of OATP1B1. Four of the studied variants were identified as loss-of-function variants and as such could make the individual harboring these variants susceptible to altered pharmacokinetics and adverse effects of substrate drugs.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11095-021-03107-8.

Protective Effects of Statin Therapy in Cirrhosis Are Limited by a Common SLCO1B1 Transporter Variant

Complications of cirrhosis and portal hypertension (PH) can be reduced by statin therapy. The common loss-of-function variant p.V174A in the solute carrier organic anion transporter gene 1B1 (SLCO1B1) gene encoding the organic anion transporting polypeptide 1B1 results in decreased hepatic uptake of statins. Our specific aim was to assess the impact of this variant in patients with cirrhosis and statin treatment while controlling for the stage of cirrhosis and other potential confounders with propensity score matching (PSM), availing of a large cohort of genotyped study patients. In total, from 1,088 patients with cirrhosis in two German academic medical centers, PSM yielded 154 patients taking statins and 154 matched controls. The effect on PH was assessed by the liver stiffness-spleen size-to-platelet score (LSPS), and complications of cirrhosis were retrospectively recorded applying consensus criteria. As hypothesized, patients on statin treatment presented less frequently with signs of PH: Esophageal varices (41% vs. 62%; P < 0.001) were less common, and LSPS (4.8 ± 11.5 vs. 5.6 ± 6.4; P = 0.01) was reduced. Correspondingly, decompensation events were also reduced in patients on statins (odds ratio [OR] = 0.54, 95% confidence interval [CI] 0.32-0.90; P = 0.02). When the variant in SLCO1B1 was present in patients on statins, esophageal varices (OR = 2.68, 95% CI 1.24-5.81; P = 0.01) and bacterial infections (OR = 2.50, 95% CI 1.14-5.47; P = 0.02) were more common as compared with wild type carriers on statins. Conclusion: In this cohort, signs and complications of PH were reduced in patients with cirrhosis treated with statins. Notably, this effect was diminished by the common loss-of-function variant in SLCO1B1. Further prospective studies in independent cohorts are warranted to confirm these genotype-specific observations.

Association of SLCO1B1 c.521T>C (rs4149056) with discontinuation of atorvastatin due to statin-associated muscle symptoms

The most common adverse drug reaction from statins are statin-associated muscle symptoms (SAMS), characterized by myopathy (weakness), myalgia (muscle pain), and commonly elevation in serum creatine kinase. All statins are substrates of the organic anion transporter 1B1 (OATP1B1; gene: SLCO1B1), albeit to different degrees. A genetic polymorphism in SLCO1B1, c.521T>C (rs4149056), markedly decreases OATP1B1 function. The literature is currently unclear as to whether SLCO1B1 c.521T>C is significantly associated with discontinuation of atorvastatin specifically due to SAMS. Our hypothesis was that individuals carrying the SLCO1B1 decreased function 521C allele are more likely to discontinue atorvastatin due to SAMS. This was a retrospective analysis of survey data from 379 Caucasians genotyped for rs4149056 and treated with atorvastatin for at least 12 months. Crude and multivariable logistic regression, adjusted for established risk factors for SAMS, determined the association of SLCO1B1 c.521T>C with discontinuation of atorvastatin due to SAMS (SLCO1B1 521T-homozygotes vs. 521C-carriers). The sample was 51% male, with a mean age of 57 years (SD = 11). Sixty-one percent of participants reported discontinuing atorvastatin due to SAMS, and 32% overall carried the 521C allele. SLCO1B1 521C-carrier status was not a significant predictor of atorvastatin discontinuation in any model: crude OR = 1.07; 95% CI, 0.68-1.66; P = 0.78 and adjusted OR = 1.07; 95% CI, 0.68-1.69; P = 0.76. The results were similar in a sub-group of participants treated with higher doses of atorvastatin (>20 mg). In summary, SLCO1B1 c.521T>C was not significantly associated with discontinuation of atorvastatin therapy due to SAMS.

Current Evidence, Challenges, and Opportunities of Physiologically Based Pharmacokinetic Models of Atorvastatin for Decision Making

Atorvastatin (ATS) is the gold-standard treatment worldwide for the management of hypercholesterolemia and prevention of cardiovascular diseases associated with dyslipidemia. Physiologically based pharmacokinetic (PBPK) models have been positioned as a valuable tool for the characterization of complex pharmacokinetic (PK) processes and its extrapolation in special sub-groups of the population, leading to regulatory recognition. Several PBPK models of ATS have been published in the recent years, addressing different aspects of the PK properties of ATS. Therefore, the aims of this review are (i) to summarize the physicochemical and pharmacokinetic characteristics involved in the time-course of ATS, and (ii) to evaluate the major highlights and limitations of the PBPK models of ATS published so far. The PBPK models incorporate common elements related to the physicochemical aspects of ATS. However, there are important differences in relation to the analyte evaluated, the type and effect of transporters and metabolic enzymes, and the permeability value used. Additionally, this review identifies major processes (lactonization, P-gp contribution, ATS-Ca solubility, simultaneous management of multiple analytes, and experimental evidence in the target population), which would enhance the PBPK model prediction to serve as a valid tool for ATS dose optimization.

Relationship of hemoglobin level and plasma coproporphyrin-I concentrations as an endogenous probe for phenotyping OATP1B

Plasma coproporphyrin‐I (CP‐I) concentration is used as a sensitive and selective endogenous probe for phenotyping organic anion transporting polypeptides 1B (OATP1B) activity in many studies. CP‐I is produced in the process of heme synthesis, but the relationship between plasma CP‐I concentrations and heme synthesis activity is unknown. In this study, we evaluated the relationship between plasma CP‐I concentration and hemoglobin level as a biomarker of heme synthesis activity. The data of 391 subjects selected from the Japanese general population were analyzed. One hundred twenty‐six participants had OATP1B1*15 allele, 11 of whom were homozygous (OATP1B1*15/*15). Multiple regression analysis identified hemoglobin level as an independent variable associated with plasma CP‐I concentration (p < 0.0001). A significant positive correlation was observed between hemoglobin level and plasma CP‐I concentration in participants without OATP1B1*15 allele (n = 265; r_s = 0.35, p < 0.0001) and with OATP1B1*15 allele (n = 126; r_s =0.27, p = 0.0022). However, Kruskal–Wallis test showed no large difference in Kruskal–Wallis statistics between the distribution of plasma CP‐I concentrations and that of ratio of plasma CP‐I to hemoglobin among six OATP1B1 polymorphism groups. These findings suggest that the hemoglobin level seems to reflect biosynthesis of CP‐I. However, correction by hemoglobin level is not required when using basal plasma CP‐I concentration for phenotyping OATP1B activity.

Pharmacokinetic/Pharmacodynamic Modelling to Describe the Cholesterol Lowering Effect of Rosuvastatin in People Living with HIV

BACKGROUND

Rosuvastatin is a lipid-lowering agent widely prescribed in people living with HIV, which is actively transported into the liver, making it a potential victim of drug-drug interactions with antiretroviral agents.

OBJECTIVES

The aims of this study were to characterise the pharmacokinetic profile of rosuvastatin and to describe the relationship between rosuvastatin concentrations and non-high-density lipoprotein (HDL)-cholesterol levels in people living with HIV.

METHODS

A population pharmacokinetic model (NONMEM) was developed to quantify the influence of demographics, clinical characteristics and comedications on rosuvastatin pharmacokinetics. This model was combined with an indirect effect model to describe non-HDL-cholesterol measurements.

RESULTS

A two-compartment model with sequential zero- and first-order absorption best fitted the 154 rosuvastatin concentrations provided by 65 people living with HIV. None of the tested covariates significantly influenced rosuvastatin pharmacokinetics. A total of 403 non-HDL cholesterol values were available for pharmacokinetic-pharmacodynamic modelling. Baseline non-HDL cholesterol decreased by 14% and increased by 12% with etravirine and antiretroviral drugs with a known impact on the lipid profile (i.e. protease inhibitors, efavirenz, cobicistat), respectively. The baseline value was surprisingly 43% lower in people living with HIV aged 80 years compared with those aged 40 years. Simulations based on the covariate-free model predicted that, under standard rosuvastatin dosages of 5 mg and 20 mg once daily, 31% and 64% of people living with HIV would achieve non-HDL-cholesterol targets, respectively.

CONCLUSIONS

The high between-subject variability that characterises both rosuvastatin pharmacokinetic and pharmacodynamic profiles remained unexplained after the inclusion of usual covariates. Considering its limited potential for drug-drug interactions with antiretroviral agents and its potent lipid-lowering effect, rosuvastatin prescription appears safe and effective in people living with HIV with hypercholesterolaemia.

CLINICAL TRIAL REGISTRATION NO

NCT03515772.

The Role of Pharmacogenomics in Contemporary Cardiovascular Therapy: A position statement from the European Society of Cardiology Working Group on Cardiovascular Pharmacotherapy

There is a strong and ever-growing body of evidence regarding the use of pharmacogenomics to inform cardiovascular pharmacology. However, there is no common position taken by international cardiovascular societies to unite diverse availability, interpretation and application of such data, nor is there recognition of the challenges of variation in clinical practice between countries within Europe. Aside from the considerable barriers to implementing pharmacogenomic testing and the complexities of clinically actioning results, there are differences in the availability of resources and expertise internationally within Europe. Diverse legal and ethical approaches to genomic testing and clinical therapeutic application also require serious thought. As direct-to-consumer genomic testing becomes more common, it can be anticipated that data may be brought in by patients themselves, which will require critical assessment by the clinical cardiovascular prescriber. In a modern, pluralistic and multi-ethnic Europe, self-identified race/ethnicity may not be concordant with genetically detected ancestry and thus may not accurately convey polymorphism prevalence. Given the broad relevance of pharmacogenomics to areas such as thrombosis and coagulation, interventional cardiology, heart failure, arrhythmias, clinical trials, and policy/regulatory activity within cardiovascular medicine, as well as to genomic and pharmacology subspecialists, this position statement attempts to address these issues at a wide-ranging level.

Statin-Induced Myopathy: Translational Studies from Preclinical to Clinical Evidence

Statins are the most prescribed and effective drugs to treat cardiovascular diseases (CVD). Nevertheless, these drugs can be responsible for skeletal muscle toxicity which leads to reduced compliance. The discontinuation of therapy increases the incidence of CVD. Thus, it is essential to assess the risk. In fact, many studies have been performed at preclinical and clinical level to investigate pathophysiological mechanisms and clinical implications of statin myotoxicity. Consequently, new toxicological aspects and new biomarkers have arisen. Indeed, these drugs may affect gene transcription and ion transport and contribute to muscle function impairment. Identifying a marker of toxicity is important to prevent or to cure statin induced myopathy while assuring the right therapy for hypercholesterolemia and counteracting CVD. In this review we focused on the mechanisms of muscle damage discovered in preclinical and clinical studies and highlighted the pathological situations in which statin therapy should be avoided. In this context, preventive or substitutive therapies should also be evaluated.

[Mitochondrial G12630A variation is associated with statin-induced myalgia in Chinese patients with coronary artery disease]

OBJECTIVE

To identify mitochondrial gene variants associated with statin-induced myalgia in Chinese patients with coronary artery disease (CHD).

METHODS

This study was conducted in a cohort of 403 patients with CHD receiving rosuvastatin therapy, among whom 341 patients had complete follow-up data concerning myalgia and 389 patients had documented measurements of plasma creatine kinase (CK) level. All these patients underwent genetic analysis using GSA chip for detecting mitochondria gene variants associated with myalgia. A logistic regression model was used to assess the association between 69 mitochondrial single-nucleotide polymorphisms (SNPs) and myopathy in 341 patients. The impact of these mutation sites on CK levels in 389 patients was evaluated by linear regression analysis.

RESULTS

G12630A variant was identified to correlate with an increased risk of myalgia in CHD patients (OR: 8.689, 95% CI: 1.586-47.6; P=0.01273), but CK levels did not differ significantly between patients with different genotypes of G12630A (P > 0.05). SNPs at T12285C and A13105G were found to significantly correlate with CK levels in these patients (P < 0.05).

CONCLUSIONS

Mitochondrial G12630A variation is associated with statin-induced myalgia in patients with CHD, indicating the necessity of different treatment strategies for patients who carry this risk allele.

Multi-factorial pharmacokinetic interactions: unraveling complexities in precision drug therapy

Introduction: Precision drug therapy requires accounting for pertinent factors in pharmacokinetic (PK) inter-individual variability (i.e., pharmacogenetics, diseases, polypharmacy, and natural product use) that can cause sub-therapeutic or adverse effects. Although each of these individual factors can alter victim drug PK, multi-factorial interactions can cause additive, synergistic, or opposing effects. Determining the magnitude and direction of these complex multi-factorial effects requires understanding the rate-limiting redundant and/or sequential PK processes for each drug.Areas covered: Perturbations in drug-metabolizing enzymes and/or transporters are integral to single- and multi-factorial PK interactions. Examples of single factor PK interactions presented include gene-drug (pharmacogenetic), disease-drug, drug-drug, and natural product-drug interactions. Examples of multi-factorial PK interactions presented include drug-gene-drug, natural product-gene-drug, gene-gene-drug, disease-natural product-drug, and disease-gene-drug interactions. Clear interpretation of multi-factorial interactions can be complicated by study design, complexity in victim drug PK, and incomplete mechanistic understanding of victim drug PK.Expert opinion: Incorporation of complex multi-factorial PK interactions into precision drug therapy requires advances in clinical decision tools, intentional PK study designs, drug-metabolizing enzyme and transporter fractional contribution determinations, systems and computational approaches (e.g., physiologically-based pharmacokinetic modeling), and PK phenotyping of progressive diseases.