Ten Years' Experience with the CYP2D6 Activity Score: A Perspective on Future Investigations to Improve Clinical Predictions for Precision Therapeutics

Delayed gametocyte clearance in Plasmodium vivax malaria is associated with polymorphisms in the cytochrome P450 reductase (CPR)

Primaquine (PQ) is the main drug used to eliminate dormant liver stages and prevent relapses in Plasmodium vivax malaria. It also has an effect on the gametocytes of Plasmodium falciparum; however, it is unclear to what extent PQ affects P. vivax gametocytes. PQ metabolism involves multiple enzymes, including the highly polymorphic CYP2D6 and the cytochrome P450 reductase (CPR). Since genetic variability can impact drug metabolism, we conducted an evaluation of the effect of CYP2D6 and CPR variants on PQ gametocytocidal activity in 100 subjects with P. vivax malaria. To determine gametocyte density, we measured the levels of pvs25 transcripts in samples taken before treatment (D0) and 72 hours after treatment (D3). Generalized estimating equations (GEEs) were used to examine the effects of enzyme variants on gametocyte densities, adjusting for potential confounding factors. Linear regression models were adjusted to explore the predictors of PQ blood levels measured on D3. Individuals with the CPR mutation showed a smaller decrease in gametocyte transcript levels on D3 compared to those without the mutation (P = 0.02, by GEE). Consistent with this, higher PQ blood levels on D3 were associated with a lower reduction in pvs25 transcripts. Based on our findings, the CPR variant plays a role in the persistence of gametocyte density in P. vivax malaria. Conceptually, our work points to pharmacogenetics as a non-negligible factor to define potential host reservoirs with the propensity to contribute to transmission in the first days of CQ-PQ treatment, particularly in settings and seasons of high Anopheles human-biting rates.

Physiologically Based Pharmacokinetic Modeling to Unravel the Drug-gene Interactions of Venlafaxine: Based on Activity Score-dependent Metabolism by CYP2D6 and CYP2C19 Polymorphisms

BACKGROUND

Venlafaxine (VEN) is a commonly utilized medication for alleviating depression and anxiety disorders. The presence of genetic polymorphisms gives rise to considerable variations in plasma concentrations across different phenotypes. This divergence in phenotypic responses leads to notable differences in both the efficacy and tolerance of the drug.

PURPOSE

A physiologically based pharmacokinetic (PBPK) model for VEN and its metabolite O-desmethylvenlafaxine (ODV) to predict the impact of CYP2D6 and CYP2C19 gene polymorphisms on VEN pharmacokinetics (PK).

METHODS

The parent-metabolite PBPK models for VEN and ODV were developed using PK-Sim® and MoBi®. Leveraging prior research, derived and implemented CYP2D6 and CYP2C19 activity score (AS)-dependent metabolism to simulate exposure in the drug-gene interactions (DGIs) scenarios. The model's performance was evaluated by comparing predicted and observed values of plasma concentration-time (PCT) curves and PK parameters values.

RESULTS

In the base models, 91.1%, 94.8%, and 94.6% of the predicted plasma concentrations for VEN, ODV, and VEN + ODV, respectively, fell within a twofold error range of the corresponding observed concentrations. For DGI scenarios, these values were 81.4% and 85% for VEN and ODV, respectively. Comparing CYP2D6 AS = 2 (normal metabolizers, NM) populations to AS = 0 (poor metabolizers, PM), 0.25, 0.5, 0.75, 1.0 (intermediate metabolizers, IM), 1.25, 1.5 (NM), and 3.0 (ultrarapid metabolizers, UM) populations in CYP2C19 AS = 2.0 group, the predicted DGI AUC0-96 h ratios for VEN were 3.65, 3.09, 2.60, 2.18, 1.84, 1.56, 1.34, 0.61, and for ODV, they were 0.17, 0.35, 0.51, 0.64, 0.75, 0.83, 0.90, 1.11, and the results were similar in other CYP2C19 groups. It should be noted that PK differences in CYP2C19 phenotypes were not similar across different CYP2D6 groups.

CONCLUSIONS

In clinical practice, the impact of genotyping on the in vivo disposition process of VEN should be considered to ensure the safety and efficacy of treatment.

Effect of isotretinoin on CYP2D6 and CYP3A activity in patients with severe acne

AIMS

Previously, retinoids have decreased CYP2D6 mRNA expression in vitro and induced CYP3A4 in vitro and in vivo. This study aimed to determine whether isotretinoin administration changes CYP2D6 and CYP3A activities in patients with severe acne.

METHODS

Thirty-three patients (22 females and 11 males, 23.5 ± 6.0 years old) expected to receive isotretinoin treatment completed the study. All participants were genotyped for CYP2D6 and CYP3A5. Participants received dextromethorphan (DM) 30 mg orally as a dual-probe substrate of CYP2D6 and CYP3A activity at two study timepoints: pre-isotretinoin treatment and with isotretinoin for at least 1 week. The concentrations of isotretinoin, DM and their metabolites were measured in 2-h postdose plasma samples and in cumulative 0-4-h urine collections using liquid chromatography-mass spectrometry.

RESULTS

In CYP2D6 extensive metabolizers, the urinary dextrorphan (DX)/DM metabolic ratio (MR) (CYP2D6 activity marker) was numerically, but not significantly, lower with isotretinoin administration compared to pre-isotretinoin (geometric mean ratio [GMR] [90% confidence interval (CI)] 0.78 [0.55, 1.11]). The urinary 3-hydroxymorphinan (3HM)/DX MR (CYP3A activity marker) was increased (GMR 1.18 [1.03, 1.35]) and the urinary DX-O-glucuronide/DX MR (proposed UGT2B marker) was increased (GMR 1.22 [1.06, 1.39]) with isotretinoin administration compared to pre-isotretinoin.

CONCLUSIONS

Administration of isotretinoin did not significantly reduce CYP2D6 activity in extensive metabolizers, suggesting that the predicted downregulation of CYP2D6 based on in vitro data does not translate into humans. We observed a modest increase in CYP3A activity (predominantly CYP3A4) with isotretinoin treatment. The data also suggest that DX glucuronidation is increased following isotretinoin administration.

Beyond Population-Level Targets for Drug Concentrations: Precision Dosing Needs Individual-Level Targets that Include Superior Biomarkers of Drug Responses

The purpose of precision dosing is to increase the chances of therapeutic success in individual patients. This is achieved in practice by adjusting doses to reach precision dosing targets determined previously in relevant populations, ideally with robust supportive evidence showing improved clinical outcomes compared with standard dosing. But is this implicit assumption of translatable population-level precision dosing targets correct and the best for all patients? In this review, the types of precision dosing targets and how they are determined are outlined, problems with the translatability of these targets to individual patients are identified, and ways forward to address these challengers are proposed. Achieving improved clinical outcomes to support precision dosing over standard dosing is currently hampered by applying population-level targets to all patients. Just as "one-dose-fits-all" may be an inappropriate philosophy for drug treatment overall, a "one-target-fits-all" philosophy may limit the broad clinical benefits of precision dosing. Defining individual-level precision dosing targets may be needed for greatest therapeutic success. Superior future precision dosing targets will integrate several biomarkers that together account for the multiple sources of drug response variability.

Solanidine is a sensitive and specific dietary biomarker for CYP2D6 activity

BACKGROUND

Individual assessment of CYP enzyme activities can be challenging. Recently, the potato alkaloid solanidine was suggested as a biomarker for CYP2D6 activity. Here, we aimed to characterize the sensitivity and specificity of solanidine as a CYP2D6 biomarker among Finnish volunteers with known CYP2D6 genotypes.

RESULTS

Using non-targeted metabolomics analysis, we identified 9152 metabolite features in the fasting plasma samples of 356 healthy volunteers. Machine learning models suggested strong association between CYP2D6 genotype-based phenotype classes with a metabolite feature identified as solanidine. Plasma solanidine concentration was 1887% higher in genetically poor CYP2D6 metabolizers (gPM) (n = 9; 95% confidence interval 755%, 4515%; P = 1.88 × 10-11), 74% higher in intermediate CYP2D6 metabolizers (gIM) (n = 89; 27%, 138%; P = 6.40 × 10-4), and 35% lower in ultrarapid CYP2D6 metabolizers (gUM) (n = 20; 64%, - 17%; P = 0.151) than in genetically normal CYP2D6 metabolizers (gNM; n = 196). The solanidine metabolites m/z 444 and 430 to solanidine concentration ratios showed even stronger associations with CYP2D6 phenotypes. Furthermore, the areas under the receiver operating characteristic and precision-recall curves for these metabolic ratios showed equal or better performances for identifying the gPM, gIM, and gUM phenotype groups than the other metabolites, their ratios to solanidine, or solanidine alone. In vitro studies with human recombinant CYP enzymes showed that solanidine was metabolized mainly by CYP2D6, with a minor contribution from CYP3A4/5. In human liver microsomes, the CYP2D6 inhibitor paroxetine nearly completely (95%) inhibited the metabolism of solanidine. In a genome-wide association study, several variants near the CYP2D6 gene associated with plasma solanidine metabolite ratios.

CONCLUSIONS

These results are in line with earlier studies and further indicate that solanidine and its metabolites are sensitive and specific biomarkers for measuring CYP2D6 activity. Since potato consumption is common worldwide, this biomarker could be useful for evaluating CYP2D6-mediated drug-drug interactions and to improve prediction of CYP2D6 activity in addition to genotyping.

Circulating Biomarkers Instead of Genotyping to Establish Metabolizer Phenotypes

Pharmacogenomics (PGx) enables personalized treatment for the prediction of drug response and to avoid adverse drug reactions. Currently, PGx mainly relies on the genetic information of absorption, distribution, metabolism, and excretion (ADME) targets such as drug-metabolizing enzymes or transporters to predict differences in the patient's phenotype. However, there is evidence that the phenotype-genotype concordance is limited. Thus, we discuss different phenotyping strategies using exogenous xenobiotics (e.g., drug cocktails) or endogenous compounds for phenotype prediction. In particular, minimally invasive approaches focusing on liquid biopsies offer great potential to preemptively determine metabolic and transport capacities. Early studies indicate that ADME phenotyping using exosomes released from the liver is reliable. In addition, pharmacometric modeling and artificial intelligence improve phenotype prediction. However, further prospective studies are needed to demonstrate the clinical utility of individualized treatment based on phenotyping strategies, not only relying on genetics. The present review summarizes current knowledge and limitations.

Risperidone plasma level, and its correlation with CYP2D6 gene polymorphism, clinical response and side effects in chronic schizophrenia patients

BACKGROUND

To explore the influence of CYP2D6 genetic polymorphism on risperidone metabolism, thereby affecting risperidone's effects and safeties in patients with chronic schizophrenia.

METHODS

Sixty-nine subjects with chronic schizophrenia treated with risperidone were recruited. CYP2D6 genotypes was determined using targeted sequencing and translated into phenotype using activity system. Risperidone plasma concentrations were measured using HPLC. Positive and Negative Symptom Scale (PANSS) and Brief Psychiatric Rating Scale (BPRS) were used to evaluate the existence and severity of psychiatric symptoms, Barnes Akathisia Scale (BAS) and Extrapyramidal Symptom Rating Scale (ESRS) for neurological side effects. Metabolic and endocrine status assess were also included.

RESULTS

The plasma drug concentrations varied hugely among individuals. Intermediate metabolizer (IM) group had higher plasma levels of RIP and dose corrected RIP concentration, RIP/9-OH-RIP ratio and C/D ratio than normal metabolizer (NM) group (p < 0.01). There was no statistic difference between responders and non-responders in dose-adjusted plasma concentrations and ratios of RIP/9-OH-RIP and C/D. The occurrence of EPS was related to active moiety levels in 4th week (p < 0.05). The prolactin (PRL) levels in two follow-ups were both significantly higher than baseline (p < 0.01). PRL change from baseline to week 4 and week 8 were both positively associated with active moiety concentration detected in week 4 (p < 0.05).

CONCLUSIONS

The risperidone plasma levels have great inter- and intraindividual variations, and are associated with the CYP2D6 phenotypes, as well as the changes in serum prolactin in patients diagnosed with chronic schizophrenia.

Novel Independent Trans- and Cis-Genetic Variants Associated with CYP2D6 Expression and Activity in Human Livers

Cytochrome P450 2D6 (CYP2D6) is a critical hepatic drug-metabolizing enzyme in humans, responsible for metabolizing approximately 20%-25% of commonly used medications such as codeine, desipramine, fluvoxamine, paroxetine, and tamoxifen. The CYP2D6 gene is highly polymorphic, resulting in substantial interindividual variability in its catalytic function and the pharmacokinetics and therapeutic outcomes of its substrate drugs. Although many functional CYP2D6 variants have been discovered and validated, a significant portion of the variability in the expression and activity of CYP2D6 remains unexplained. In this study, we performed a genome-wide association study (GWAS) to identify novel variants associated with CYP2D6 protein expression in individual human livers, followed by a conditional analysis to control for the effect of functional CYP2D6 star alleles. We also examined their impact on hepatic CYP2D6 activity. Genotyping on a genome-wide scale was achieved using the Illumina Multi-Ethnic Genotyping Array (MEGA). A data-independent acquisition (DIA)-based proteomics method was used to quantify CYP2D6 protein concentrations. CYP2D6 activity was determined by measuring the dextromethorphan O-demethylation in individual human liver s9 fractions. The GWAS identified 44 single nuclear polymorphisms (SNPs) that are significantly associated with CYP2D6 protein expressions with a P value threshold of 5.0 × 10-7 After the conditional analysis, five SNPs, including the cis-variants rs1807493 and rs1062753 and the trans-variants rs4073010, rs729559, and rs80274432, emerged as independent variants significantly correlated with hepatic CYP2D6 protein expressions. Notably, four of these SNPs, except for rs80274432, also exhibited a significant association with CYP2D6 activities in human livers, suggesting their potential as novel and independent cis- and trans-variants regulating CYP2D6. SIGNIFICANT STATEMENT: Using individual human livers, we identified four novel cis- and trans-pQTLs/aQTLs (protein quantitative trait loci/activity quantitative trait loci) of Cytochrome P450 2D6 (CYP2D6) that are independent from known functional CYP2D6 star alleles. This study connects the CYP2D6 gene expression and activity, enhancing our understanding of the genetic variants associated with CYP2D6 protein expression and activity, potentially advancing our insight into the interindividual variability in CYP2D6 substrate medication response.

CYP2D6 Activity Is Correlated with Changes in Plasma Concentrations of Taurocholic Acid during Pregnancy and Postpartum in CYP2D6 Extensive Metabolizers

Cytochrome P450 2D6 (CYP2D6) is involved in the metabolism of >20% of marketed drugs. CYP2D6 expression and activity exhibit high interindividual variability and is induced during pregnancy. The farnesoid X receptor (FXR) is a transcriptional regulator of CYP2D6 that is activated by bile acids. In pregnancy, elevated plasma bile acid concentrations are associated with maternal and fetal risks. However, modest changes in bile acid concentrations may occur during healthy pregnancy, thereby altering FXR signaling. A previous study demonstrated that hepatic tissue concentrations of bile acids positively correlated with the hepatic mRNA expression of CYP2D6. This study sought to characterize the plasma bile acid metabolome in healthy women (n = 47) during midpregnancy (25-28 weeks gestation) and ≥3 months postpartum and to determine if plasma bile acids correlate with CYP2D6 activity. It is hypothesized that during pregnancy, plasma bile acids would favor less hydrophobic bile acids (cholic acid vs. chenodeoxycholic acid) and that plasma concentrations of cholic acid and its conjugates would positively correlate with the urinary ratio of dextrorphan/dextromethorphan. At 25-28 weeks gestation, taurine-conjugated bile acids comprised 23% of the quantified serum bile acids compared with 7% ≥3 months postpartum. Taurocholic acid positively associated with the urinary ratio of dextrorphan/dextromethorphan, a biomarker of CYP2D6 activity. Collectively, these results confirm that the bile acid plasma metabolome differs between pregnancy and postpartum and provide evidence that taurocholic acid may impact CYP2D6 activity during pregnancy. SIGNIFICANCE STATEMENT: Bile acid homeostasis is altered in pregnancy, and plasma concentrations of taurocholic acid positively correlate with CYP2D6 activity. Differences between plasma and/or tissue concentrations of farnesoid X receptor ligands such as bile acids may contribute to the high interindividual variability in CYP2D6 expression and activity.

Pharmacogenomic variation in the Malagasy population: implications for the antimalarial drug primaquine metabolism

Aim: Antimalarial primaquine (PQ) eliminates liver hypnozoites of Plasmodium vivax. CYP2D6 gene variation contributes to PQ therapeutic failure. Additional gene variation may contribute to PQ efficacy. Information on pharmacogenomic variation in Madagascar, with vivax malaria and a unique population admixture, is scanty. Methods: The authors performed genome-wide genotyping of 55 Malagasy samples and analyzed data with a focus on a set of 28 pharmacogenes most relevant to PQ. Results: Mainly, the study identified 110 coding or splicing variants, including those that, based on previous studies in other populations, may be implicated in PQ response and copy number variation, specifically in chromosomal regions that contain pharmacogenes. Conclusion: With this pilot information, larger genome-wide association analyses with PQ metabolism and response are substantially more feasible.

Implementation of CYP2D6-guided opioid therapy at Cincinnati Children's Hospital Medical Center

PURPOSE

We describe the implementation of CYP2D6-focused pharmacogenetic testing to guide opioid prescribing in a quaternary care, nonprofit pediatric academic medical center.

SUMMARY

Children are often prescribed oral opioids after surgeries, for cancer pain, and occasionally for chronic pain. In 2004, Cincinnati Children's Hospital Medical Center implemented pharmacogenetic testing for CYP2D6 metabolism phenotype to inform codeine prescribing. The test and reports were updated to align with changes over time in the testing platform, the interpretation of genotype to phenotype, the electronic health record, and Food and Drug Administration (FDA) guidance. The use of the test increased when a research project required testing and decreased as prescribing of oxycodone increased due to FDA warnings about codeine. Education about the opioid-focused pharmacogenetic test was provided to prescribers (eg, the pain and sickle cell teams) as well as patients and families. Education and electronic health record capability increased provider compliance with genotype-guided postsurgical prescribing of oxycodone, although there was a perceived lack of utility for oxycodone prescribing.

CONCLUSION

The implementation of pharmacogenetic testing to inform opioid prescribing for children has evolved with accumulating evidence and guidelines, requiring changes in reporting of results and recommendations.

Impact of Pregnancy and Vitamin A Supplementation on CYP2D6 Activity

The mechanism of cytochrome P450 2D6 (CYP2D6) induction during pregnancy has not been evaluated in humans. This study assessed the changes in CYP2D6 and CYP3A activities during pregnancy and postpartum, and the effect of vitamin A administration on CYP2D6 activity. Forty-seven pregnant CYP2D6 extensive metabolizers (with CYP2D6 activity scores of 1 to 2) received dextromethorphan (DM) 30 mg orally as a single dose during 3 study windows (at 25 to 28 weeks of gestation, study day 1; at 28 to 32 weeks of gestation, study day 2; and at ≥3 months postpartum, study day 3). Participants were randomly assigned to groups with no supplemental vitamin A (control) or with supplemental vitamin A (10 000 IU/day orally for 3 to 4 weeks) after study day 1. Concentrations of DM and its metabolites, dextrorphan (DX) and 3-hydroxymorphinan (3HM), were determined from a 2-hour post-dose plasma sample and cumulative 4-hour urine sample using liquid chromatography-mass spectrometry. Change in CYP2D6 activity was assessed using DX/DM plasma and urine metabolic ratios. The activity change in CYP3A was also assessed using the 3HM/DM urine metabolic ratio. The DX/DM urine ratio was significantly higher (43%) in pregnancy compared with postpartum (P = .03), indicating increased CYP2D6 activity. The DX/DM plasma ratio was substantially higher in the participants, with an activity score of 1.0 during pregnancy (P = .04) compared with postpartum. The 3HM/DM urinary ratio was significantly higher (92%) during pregnancy, reflecting increased CYP3A activity (P = .02). Vitamin A supplementation did not change CYP2D6 activity during pregnancy; however, plasma all-trans retinoic acid (atRA) concentrations were positively correlated with increased CYP2D6 activity during pregnancy and postpartum. Further research is needed to elucidate the mechanisms of increased CYP2D6 activity during pregnancy.

Physiologically-Based Pharmacokinetic Model Development, Validation, and Application for Prediction of Eliglustat Drug-Drug Interactions

Eliglustat is a glucosylceramide synthase inhibitor indicated as a long-term substrate reduction therapy for adults with type 1 Gaucher disease, a lysosomal rare disease. It is primarily metabolized by cytochrome P450 2D6 (CYP2D6), and variants in the gene encoding this enzyme are important determinants of eliglustat pharmacokinetics (PK) and drug-drug interactions (DDIs). The existing drug label addresses the DDIs to some extent but has omitted scenarios where both metabolizing CYPs (2D6 and 3A4) are mildly or moderately inhibited. The objectives of this study were (i) to develop and validate an eliglustat physiologically-based pharmacokinetic (PBPK) model with and without drug interactions, (ii) to simulate untested DDI scenarios, and (iii) to explore potential dosing flexibility using lower dose strength of eliglustat (commercially not available). PK data from healthy adults receiving eliglustat with or without interacting drugs were obtained from literature and used for the PBPK model development and validation. The model-predicted single-dose and steady-state maximum concentration (C_max ) and area under the concentration-time curve (AUC) of eliglustat were within 50-150% of the observed values when eliglustat was administered alone or coadministered with ketoconazole or paroxetine. Then as model-based simulations, we illustrated eliglustat exposure as a victim of interaction when coadministered with fluvoxamine following the US Food and Drug Administration (FDA) dosing recommendations. Second, we showed that with lower eliglustat doses (21 mg, 42 mg once daily) the exposure in participants of intermediate and poor metabolizer phenotypes was within the outlined safety margin (C_max <250 ng/mL) when eliglustat was administered with ketoconazole, where the current recommendation is a contraindication of coadministration (84 mg). The present study demonstrated that patients with CYP2D6 deficiency may benefit from lower doses of eliglustat.

Prediction of Drug-Drug-Gene Interaction Scenarios of (E)-Clomiphene and Its Metabolites Using Physiologically Based Pharmacokinetic Modeling

Clomiphene, a selective estrogen receptor modulator (SERM), has been used for the treatment of anovulation for more than 50 years. However, since (E)-clomiphene ((E)-Clom) and its metabolites are eliminated primarily via Cytochrome P450 (CYP) 2D6 and CYP3A4, exposure can be affected by CYP2D6 polymorphisms and concomitant use with CYP inhibitors. Thus, clomiphene therapy may be susceptible to drug-gene interactions (DGIs), drug-drug interactions (DDIs) and drug-drug-gene interactions (DDGIs). Physiologically based pharmacokinetic (PBPK) modeling is a tool to quantify such DGI and DD(G)I scenarios. This study aimed to develop a whole-body PBPK model of (E)-Clom including three important metabolites to describe and predict DGI and DD(G)I effects. Model performance was evaluated both graphically and by calculating quantitative measures. Here, 90% of predicted C_max and 80% of AUC_last values were within two-fold of the corresponding observed value for DGIs and DD(G)Is with clarithromycin and paroxetine. The model also revealed quantitative contributions of different CYP enzymes to the involved metabolic pathways of (E)-Clom and its metabolites. The developed PBPK model can be employed to assess the exposure of (E)-Clom and its active metabolites in as-yet unexplored DD(G)I scenarios in future studies.

Global perspectives on CYP2D6 associations with primaquine metabolism and Plasmodium vivax radical cure

Clinical trial and individual patient treatment outcomes have produced accumulating evidence that effective primaquine (PQ) treatment of Plasmodium vivax and P. ovale liver stage hypnozoites is associated with genetic variation in the human cytochrome P450 gene, CYP2D6. Successful PQ treatment of individual and population-wide infections by the Plasmodium species that generate these dormant liver stage forms is likely to be necessary to reach elimination of malaria caused by these parasites globally. Optimizing safe and effective PQ treatment will require coordination of efforts between the malaria and pharmacogenomics research communities.

Physiologically based pharmacokinetic (PBPK) modeling of the role of CYP2D6 polymorphism for metabolic phenotyping with dextromethorphan

The cytochrome P450 2D6 (CYP2D6) is a key xenobiotic-metabolizing enzyme involved in the clearance of many drugs. Genetic polymorphisms in CYP2D6 contribute to the large inter-individual variability in drug metabolism and could affect metabolic phenotyping of CYP2D6 probe substances such as dextromethorphan (DXM). To study this question, we (i) established an extensive pharmacokinetics dataset for DXM; and (ii) developed and validated a physiologically based pharmacokinetic (PBPK) model of DXM and its metabolites dextrorphan (DXO) and dextrorphan O-glucuronide (DXO-Glu) based on the data. Drug-gene interactions (DGI) were introduced by accounting for changes in CYP2D6 enzyme kinetics depending on activity score (AS), which in combination with AS for individual polymorphisms allowed us to model CYP2D6 gene variants. Variability in CYP3A4 and CYP2D6 activity was modeled based on in vitro data from human liver microsomes. Model predictions are in very good agreement with pharmacokinetics data for CYP2D6 polymorphisms, CYP2D6 activity as described by the AS system, and CYP2D6 metabolic phenotypes (UM, EM, IM, PM). The model was applied to investigate the genotype-phenotype association and the role of CYP2D6 polymorphisms for metabolic phenotyping using the urinary cumulative metabolic ratio (UCMR), DXM/(DXO + DXO-Glu). The effect of parameters on UCMR was studied via sensitivity analysis. Model predictions indicate very good robustness against the intervention protocol (i.e. application form, dosing amount, dissolution rate, and sampling time) and good robustness against physiological variation. The model is capable of estimating the UCMR dispersion within and across populations depending on activity scores. Moreover, the distribution of UCMR and the risk of genotype-phenotype mismatch could be estimated for populations with known CYP2D6 genotype frequencies. The model can be applied for individual prediction of UCMR and metabolic phenotype based on CYP2D6 genotype. Both, model and database are freely available for reuse.

Effects of CYP2D6*4 polymorphism on the steady-state concentration of paroxetine in patients diagnosed with depressive episode and comorbid alcohol use disorder

INTRODUCTION

Selective serotonin reuptake inhibitors have a common and increasing use for the treatment of patients diagnosed with depressive disorders. Some of them do not respond adequately to therapy, and numerous previous studies have indicated an increased risk of type A adverse drug reactions.

OBJECTIVE

The objective of our study was to evaluate the effect of 1846G>A polymorphism of CYP2D6 on the concentration/dose ratio of paroxetine.

MATERIAL AND METHODS

The study enrolled 267 patients with depressive episode (average age, 40.3 ± 14.3 years). Therapy included paroxetine in an average daily dose of 25.1 ± 9.5 mg per day. The efficacy and safety rates of treatment were evaluated using the international psychometric scales. For genotyping, we performed the real-time polymerase chain reaction. Therapeutic drug monitoring has been performed using high-performance liquid chromatography mass spectrometry (HPLC-MS/MS).

RESULTS

Our study revealed the statistically significant results in terms of treatment efficacy (Hamilton Depression Rating Scale scores): (GG) 2.0 [1.0; 3.0] and (GA) 4.0 [2.0; 5.0], p < 0.001; meanwhile, no statistically significant results were obtained for the safety profile (Udvalg for Kliniske Undersogelser (UKU) Scale scores): (GG) 3.0 [2.0; 3.0] and (GA) 3.0 [3.0; 4.0], p = 0.056. We revealed the statistically significant results for the concentration/dose ratio of paroxetine in patients with different genotypes: (GG) 2.803 [2.154; 4.098] and (GA) 5.098 [3.560; 7.241], p < 0.001.

CONCLUSION

The effect of CYP2D6*4 genetic polymorphism on the efficacy profile of paroxetine was demonstrated in a group of 267 patients with depressive disorder.

A longitudinal study of cytochrome P450 2D6 (CYP2D6) activity during adolescence

CYP2D6 substrates are among the most highly prescribed medications in teenagers and also commonly associated with serious adverse events. To investigate the relative contributions of genetic variation, growth, and development on CYP2D6 activity during puberty, healthy children and adolescents 7-15 years of age at enrollment participated in a longitudinal phenotyping study involving administration of 0.3 mg/kg dextromethorphan (DM) and 4-h urine collection every 6 months for 3 years (7 total visits). At each visit, height, weight, and sexual maturity were recorded, and CYP2D6 activity was determined as the urinary molar ratio of DM to its metabolite dextrorphan (DX). A total of 188 participants completed at least one visit, and 102 completed all seven study visits. Following univariate analysis, only CYP2D6 activity score (p < 0.001), urinary pH (p < 0.001), weight (p = 0.018), and attention-deficit/hyperactivity disorder (ADHD) diagnosis (p < 0.001) were significantly correlated with log(DM/DX). Results of linear mixed model analysis with random intercept, random slope covariance structure revealed that CYP2D6 activity score had the strongest effect on log(DM/DX), with model-estimated average log(DM/DX) being 3.8 SDs higher for poor metabolizers than for patients with activity score 3. A moderate effect on log(DM/DX) was observed for sex, and smaller effects were observed for ADHD diagnosis and urinary pH. The log(DM/DX) did not change meaningfully with age or pubertal development. CYP2D6 genotype remains the single, largest determinant of variability in CYP2D6 activity during puberty. Incorporation of genotype-based dosing guidelines should be considered for CYP2D6 substrates given the prevalent use of these agents in this pediatric age group.

From Croatian Roma to 1000 Genomes: The Story of the CYP2D6 Gene Promoter and Enhancer SNPs

The CYP2D6 gene encodes an enzyme responsible for the metabolism of ~20% of clinically prescribed drugs. In this study, 18 SNPs from the enhancer and promoter regions of CYP2D6 in 323 Roma from Croatia were genotyped, to find out whether the demographic history of Roma affected the distribution of the studied SNPs and their linkage disequilibrium (LD) values, with the major SNPs defining the CYP2D6 star alleles. No differences were found between the three Roma groups in allele and genotype frequencies. The distribution of LD values of Roma was compared with LD values of European and Asian populations. Regulatory CYP2D6 SNPs (rs5758550, rs28624811, rs1080985 and rs1080983) showed similar distribution and the highest LDs with rs16947 from the gene-coding region in all populations. In the promoter region, a complete LD between rs1080989 and rs28588594, and between rs1080983 and rs28624811, was found in Croatian Roma and investigated populations from 1000 genomes. A high LD was also found between rs1080985 from the promoter and rs5758550 from the enhancer region. SNP rs28735595 from the gene promoter region had the highest LD, with two gene region SNPs, rs1058164 and rs1135840. To conclude, the Croatian Roma population shows an LD pattern of the CYP2D6 gene region similar to the 1000 Genomes European and Asian populations.

Physiologically Based Pharmacokinetic Modeling to Describe the CYP2D6 Activity Score-Dependent Metabolism of Paroxetine, Atomoxetine and Risperidone

The cytochrome P450 2D6 (CYP2D6) genotype is the single most important determinant of CYP2D6 activity as well as interindividual and interpopulation variability in CYP2D6 activity. Here, the CYP2D6 activity score provides an established tool to categorize the large number of CYP2D6 alleles by activity and facilitates the process of genotype-to-phenotype translation. Compared to the broad traditional phenotype categories, the CYP2D6 activity score additionally serves as a superior scale of CYP2D6 activity due to its finer graduation. Physiologically based pharmacokinetic (PBPK) models have been successfully used to describe and predict the activity score-dependent metabolism of CYP2D6 substrates. This study aimed to describe CYP2D6 drug–gene interactions (DGIs) of important CYP2D6 substrates paroxetine, atomoxetine and risperidone by developing a substrate-independent approach to model their activity score-dependent metabolism. The models were developed in PK-Sim^®, using a total of 57 plasma concentration–time profiles, and showed good performance, especially in DGI scenarios where 10/12, 5/5 and 7/7 of DGI AUC_last ratios and 9/12, 5/5 and 7/7 of DGI C_max ratios were within the prediction success limits. Finally, the models were used to predict their compound’s exposure for different CYP2D6 activity scores during steady state. Here, predicted DGI AUC_ss ratios were 3.4, 13.6 and 2.0 (poor metabolizers; activity score = 0) and 0.2, 0.5 and 0.95 (ultrarapid metabolizers; activity score = 3) for paroxetine, atomoxetine and risperidone active moiety (risperidone + 9-hydroxyrisperidone), respectively.

Quantitative Proteomics in Translational Absorption, Distribution, Metabolism, and Excretion and Precision Medicine

A reliable translation of in vitro and preclinical data on drug absorption, distribution, metabolism, and excretion (ADME) to humans is important for safe and effective drug development. Precision medicine that is expected to provide the right clinical dose for the right patient at the right time requires a comprehensive understanding of population factors affecting drug disposition and response. Characterization of drug-metabolizing enzymes and transporters for the protein abundance and their interindividual as well as differential tissue and cross-species variabilities is important for translational ADME and precision medicine. This review first provides a brief overview of quantitative proteomics principles including liquid chromatography-tandem mass spectrometry tools, data acquisition approaches, proteomics sample preparation techniques, and quality controls for ensuring rigor and reproducibility in protein quantification data. Then, potential applications of quantitative proteomics in the translation of in vitro and preclinical data as well as prediction of interindividual variability are discussed in detail with tabulated examples. The applications of quantitative proteomics data in physiologically based pharmacokinetic modeling for ADME prediction are discussed with representative case examples. Finally, various considerations for reliable quantitative proteomics analysis for translational ADME and precision medicine and the future directions are discussed. SIGNIFICANCE STATEMENT: Quantitative proteomics analysis of drug-metabolizing enzymes and transporters in humans and preclinical species provides key physiological information that assists in the translation of in vitro and preclinical data to humans. This review provides the principles and applications of quantitative proteomics in characterizing in vitro, ex vivo, and preclinical models for translational research and interindividual variability prediction. Integration of these data into physiologically based pharmacokinetic modeling is proving to be critical for safe, effective, timely, and cost-effective drug development.

Physiologically-based pharmacokinetic modeling of dextromethorphan to investigate interindividual variability within CYP2D6 activity score groups

This study provides a whole‐body physiologically‐based pharmacokinetic (PBPK) model of dextromethorphan and its metabolites dextrorphan and dextrorphan O‐glucuronide for predicting the effects of cytochrome P450 2D6 (CYP2D6) drug‐gene interactions (DGIs) on dextromethorphan pharmacokinetics (PK). Moreover, the effect of interindividual variability (IIV) within CYP2D6 activity score groups on the PK of dextromethorphan and its metabolites was investigated. A parent‐metabolite‐metabolite PBPK model of dextromethorphan, dextrorphan, and dextrorphan O‐glucuronide was developed in PK‐Sim and MoBi. Drug‐dependent parameters were obtained from the literature or optimized. Plasma concentration‐time profiles of all three analytes were gathered from published studies and used for model development and model evaluation. The model was evaluated comparing simulated plasma concentration‐time profiles, area under the concentration‐time curve from the time of the first measurement to the time of the last measurement (AUC_last) and maximum concentration (C_max) values to observed study data. The final PBPK model accurately describes 28 population plasma concentration‐time profiles and plasma concentration‐time profiles of 72 individuals from four cocktail studies. Moreover, the model predicts CYP2D6 DGI scenarios with six of seven DGI AUC_last and seven of seven DGI C_max ratios within the acceptance criteria. The high IIV in plasma concentrations was analyzed by characterizing the distribution of individually optimized CYP2D6 k_cat values stratified by activity score group. Population simulations with sampling from the resulting distributions with calculated log‐normal dispersion and mean parameters could explain a large extent of the observed IIV. The model is publicly available alongside comprehensive documentation of model building and model evaluation.

Influence of CYP2D6, CYP3A4 and CYP2C19 Genotypes on Recurrence of Plasmodium vivax

Background

The influence of the CYPs (cytochrome P-450) in the success of antimalarial therapy remains uncertain. In this study, the association of CYP2D6, CYP2C19 and CYP3A4 polymorphisms and predicted phenotypes with malaria recurrence was investigated.

Methods

After diagnosis of vivax malaria, individuals treated at a reference center in Manaus were followed up for 180 days. Patients were separated into two groups: a recurrence group and a non-recurrence group. Genotyping of CYP2D6, CYP2C19 and CYP3A4 was performed using a TaqMan™ assay and real-time PCR.

Findings

The frequencies of decreased-function and normal-function alleles and phenotypes for all CYPs were similar between the groups, except for the CYP2D6*2xN allele (p=0.047) and the CYP2D6 gUM phenotype (p=0.057), which were more frequent in individuals without recurrence. Despite this, the CYP2D6, CYP2C19 and CYP3A4 genotypes had no association with an increased risk of recurrence. CYPs polymorphisms also had no influence in parasite clearance, neither in the time nor the number of recurrence episodes. MAIN

Conclusion

This prospective cohort study demonstrated that CYP2D6, CYP2C19 and CYP3A4 polymorphisms have no influence on malaria recurrence. Nonetheless, our findings suggest that the CYP2D6 predicted ultrarapid phenotype was less susceptible to recurrence, and that patients with the CYP2D6 gUM phenotype are less susceptible to primaquine failure. Additional investigation of pharmacogenetics and pharmacokinetics are needed before implementing CYP analysis to better orientate individualized radical treatment of vivax malaria in reference centers that treat patients with multiple recurrences.

Antipsychotic medication for women with schizophrenia spectrum disorders

There are significant differences between men and women in the efficacy and tolerability of antipsychotic drugs. Here, we provide a comprehensive overview of what is currently known about the pharmacokinetics and pharmacodynamics of antipsychotics in women with schizophrenia spectrum disorders (SSDs) and translate these insights into considerations for clinical practice. Slower drug absorption, metabolism and excretion in women all lead to higher plasma levels, which increase the risk for side-effects. Moreover, women reach higher dopamine receptor occupancy compared to men at similar serum levels, since oestrogens increase dopamine sensitivity. As current treatment guidelines are based on studies predominantly conducted in men, women are likely to be overmedicated by default. The risk of overmedicating generally increases when sex hormone levels are high (e.g. during ovulation and gestation), whereas higher doses may be required during low-hormonal phases (e.g. during menstruation and menopause). For premenopausal women, with the exceptions of quetiapine and lurasidone, doses of antipsychotics should be lower with largest adjustments required for olanzapine. Clinicians should be wary of side-effects that are particularly harmful in women, such as hyperprolactinaemia which can cause oestrogen deficiency and metabolic symptoms that may cause cardiovascular diseases. Given the protective effects of oestrogens on the course of SSD, oestrogen replacement therapy should be considered for postmenopausal patients, who are more vulnerable to side-effects and yet require higher dosages of most antipsychotics to reach similar efficacy. In conclusion, there is a need for tailored, female-specific prescription guidelines, which take into account adjustments required across different phases of life.

The Impact of the CYP2D6 "Enhancer" Single Nucleotide Polymorphism on CYP2D6 Activity

rs5758550 has been associated with enhanced transcription and suggested to be a useful marker of CYP2D6 activity. As there are limited and inconsistent data regarding the utility of this distant "enhancer" single nucleotide polymorphism (SNP), our goal was to further assess the impact of rs5758550 on CYP2D6 activity toward two probe substrates, atomoxetine (ATX) and dextromethorphan (DM), using in vivo urinary metabolite (DM; n = 188) and pharmacokinetic (ATX; n = 70) and in vitro metabolite formation (ATX and DM; n = 166) data. All subjects and tissues were extensively genotyped, the "enhancer" SNP phased with established CYP2D6 haplotypes either computationally or experimentally, and the impact on CYP2D6 activity investigated using several linear models of varying complexity to determine the proportion of variability in CYP2D6 activity captured by each model. For all datasets and models, the "enhancer" SNP had no or only a modest impact on CYP2D6 activity prediction. An increased effect, when present, was more pronounced for ATX than DM suggesting potential substate-dependency. In addition, CYP2D6*2 alleles with the "enhancer" SNP were associated with modestly higher metabolite formation rates in vitro, but not in vivo; no effect was detected for CYP2D6*1 alleles with "enhancer" SNP. In summary, it remains inconclusive whether the small effects detected in this investigation are indeed caused by the "enhancer" SNP or are rather due to its incomplete linkage with other variants within the gene. Taken together, there does not appear to be sufficient evidence to warrant the "enhancer" SNP be included in clinical CYP2D6 pharmacogenetic testing.

Measurements of 5,6 orthoquinone, surrogate for presumed active primaquine metabolite 5-hydroxyprimaquine, in the urine of Cambodian adults

The active metabolites of primaquine, in particular 5-hydroxyprimaquine, likely responsible for clearance of dormant hypnozoites, are produced through the hepatic CYP450 2D6 (CYP2D6) enzymatic pathway. With the inherent instability of 5-hydroxyprimaquine, a stable surrogate, 5,6 orthoquinone, can now be detected and measured in the urine as part of primaquine pharmacokinetic studies. This study performed CYP450 2D6 genotyping and primaquine pharmacokinetic testing, to include urine 5,6 orthoquinone, in 27 healthy adult Cambodians, as a preliminary step to prepare for future clinical studies assessing primaquine efficacy for Plasmodium vivax infections. The CYP2D6 *10 reduced activity allele was found in 57% of volunteers, and the CYP2D6 genotypes were dominated by *1/*10 (33%) and *10/*10 (30%). Predicted phenotypes were evenly split between Normal Metabolizer (NM) and Intermediate Metabolizer (IM) except one volunteer with a gene duplication and unclear phenotype, classifying as either IM or NM. Median plasma PQ area under the curve (AUC) was lower in the NM group (460 hr*ng/mL) compared to the IM group (561 hr*ng/mL), although not statistically significant. Similar to what has been found in the US study, no 5,6 orthoquinone was detected in the plasma. The urine creatinine-corrected 5,6 orthoquinone AUC in the NM group was almost three times higher than in the IM group, with peak measurements (T max ) at 4 hours. Although there is variation among individuals, future studies examining the relationship between the levels of urine 5,6 orthoquinone and primaquine radical cure efficacy could result in a metabolism biomarker predictive of radical cure.

An update on pharmacogenetic factors influencing the metabolism and toxicity of artemisinin-based combination therapy in the treatment of malaria

INTRODUCTION

Artemisinin-based combination therapies (ACTs) are recommended first-line antimalarials for uncomplicated Plasmodium falciparum malaria. Pharmacokinetic/pharmacodynamic variation associated with ACT drugs and their effect is documented. It is accepted to an extent that inter-individual variation is genetically driven, and should be explored for optimized antimalarial use.

AREAS COVERED

We provide an update on the pharmacogenetics of ACT antimalarial disposition. Beyond presently used antimalarials, we also refer to information available for the most notable next-generation drugs under development. The bibliographic approach was based on multiple Boolean searches on PubMed covering all recent publications since our previous review.

EXPERT OPINION

The last 10 years have witnessed an increase in our knowledge of ACT pharmacogenetics, including the first clear examples of its contribution as an exacerbating factor for drug-drug interactions. This knowledge gap is still large and is likely to widen as a new wave of antimalarial drug is looming, with few studies addressing their pharmacogenetics. Clinically useful pharmacogenetic markers are still not available, in particular, from an individual precision medicine perspective. A better understanding of the genetic makeup of target populations can be valuable for aiding decisions on mass drug administration implementation concerning region-specific antimalarial drug and dosage options.

Genetic Variation of G6PD and CYP2D6: Clinical Implications on the Use of Primaquine for Elimination of Plasmodium vivax

Primaquine, an 8-aminoquinoline, is the only medication approved by the World Health Organization to treat the hypnozoite stage of Plasmodium vivax and P. ovale malaria. Relapse, triggered by activation of dormant hypnozoites in the liver, can occur weeks to years after primary infection, and provides the predominant source of transmission in endemic settings. Hence, primaquine is essential for individual treatment and P. vivax elimination efforts. However, primaquine use is limited by the risk of life-threatening acute hemolytic anemia in glucose-6-phosphate dehydrogenase (G6PD) deficient individuals. More recently, studies have demonstrated decreased efficacy of primaquine due to cytochrome P450 2D6 (CYP2D6) polymorphisms conferring an impaired metabolizer phenotype. Failure of standard primaquine therapy has occurred in individuals with decreased or absent CYP2D6 activity. Both G6PD and CYP2D6 are highly polymorphic genes, with considerable geographic and interethnic variability, adding complexity to primaquine use. Innovative strategies are required to overcome the dual challenge of G6PD deficiency and impaired primaquine metabolism. Further understanding of the pharmacogenetics of primaquine is key to utilizing its full potential. Accurate CYP2D6 genotype-phenotype translation may optimize primaquine dosing strategies for impaired metabolizers and expand its use in a safe, efficacious manner. At an individual level the current challenges with G6PD diagnostics and CYP2D6 testing limit clinical implementation of pharmacogenetics. However, further characterisation of the overlap and spectrum of G6PD and CYP2D6 activity may optimize primaquine use at a population level and facilitate region-specific dosing strategies for mass drug administration. This precision public health approach merits further investigation for P. vivax elimination.

Prediction of the CYP2D6 enzymatic activity based on investigating of the CYP2D6 genotypes around the vivax malaria patients in Yunnan Province, China

Background

In recent years, the incidence rate of vivax malaria recurrence still had 3.1% in Yunnan Province population after eradication therapy using primaquine (PQ). In order to understand the specific failure reasons for preventing vivax malaria relapses, a preliminary exploration on the CYP2D6 enzyme activity was carried out in the vivax malaria patients in Yunnan Province population by analysing mutational polymorphism in the coding region of CYP2D6 gene.

Methods

Blood samples were collected from vivax malaria patients with suspected relapse (SR) and non-relapsed (NR) malaria in Yunnan Province. The DNA fragments containing 9 exons regions of human CYP2D6 gene were amplified by performing PCR and sequenced. The sequencing results were aligned by using DNAStar 11.0 to obtain the coding DNA sequence (CDS) of CYP2D6 gene. DnaSP 6.11.01 software was used to identify mutant polymorphisms and haplotypes of the CDS chain. The waterfall function of GenVisR package in R was utilized to visualize the mutational landscape. The alleles of CYP2D6 gene were identified according to the criteria prescribed by Human Cytochrome P450 (CYP) Allele Nomenclature Committee Database and the CYP2D6 enzyme activity was predicted based on diploid genotype.

Results

A total of 320 maternal CDS chains, including 63 from SR group and 257 from NR group, were obtained. Twelve mutant loci, including c.31 (rs769259), c.100 (rs1065852), c.271 (rs28371703), c.281 (rs28371704), c.294 (rs28371705), c.297 (rs200269944), c.336 (rs1081003), c.408 (rs1058164), c.505 (rs5030865), c.801 (rs28371718), c.886 (rs16947), and c.1,457 (rs1135840) were observed on the 640 CDS chains (including 320 maternal and 320 paternal chains). The high-frequency mutation at rs1135840 (0.703) and low-frequency mutation, such as rs28371703, were detected only in the SR group. The frequency of mutant rs1058164 and rs1135840 were significantly increased in the SR group (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${x}^{2}$$\end{document}x2= 4.468, 5.889, P < 0.05), as opposed to the NR group. Of the 23 haplotypes (from Hap_1 to Hap_23), the nomenclatures of 11 allelic forms could be found: Hap_3 was non-mutant, Hap_2 accounted for the highest frequency (36.9%, 236/640), and Hap_9 had the most complex sequence structure, containing 7 loci mutations. Allele *10 was the most frequent among these genotypes (0.423). Among the allele *10 standard named genotypes, *1/*10, *1/*1 and *2/*10 were significantly more frequent in the NR group (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${x}^{2}$$\end{document}x2= 3.911, P < 0.05) and all showed uncompromised enzyme activity; the impaired genotype *10/*39 was more frequent in the SR group (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${x}^{2}$$\end{document}x2= 10.050, P < 0.05), and genotype *4/*4was detected only in the SR group.

Conclusion

In the patients receiving PQ dosage in Yunnan Province population, both rs1135840 single nucleotide polymorphism and *10 allele form was common in the CYP2D6 gene. Low-frequency mutation sites, such as rs28371703, were only presented in patients with vivax malaria relapse.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12936-021-03988-5.

Novel models for the prediction of drug-gene interactions

INTRODUCTION

Adverse drug reactions (ADRs) are among the leading causes of death, and frequently associated with drug-gene interactions (DGIs). In addition to pharmacogenomic programs for implementation of genetic preemptive testing into clinical practice, mathematical modeling can help to understand, quantify and predict the effects of DGIs in vivo. Moreover, modeling can contribute to optimize prospective clinical drug trial activities and to reduce DGI-related ADRs.

AREAS COVERED

Approaches and challenges of mechanistical DGI implementation and model parameterization are discussed for population pharmacokinetic and physiologically based pharmacokinetic models. The broad spectrum of published DGI models and their applications is presented, focusing on the investigation of DGI effects on pharmacology and model-based dose adaptations.

EXPERT OPINION

Mathematical modeling provides an opportunity to investigate complex DGI scenarios and can facilitate the development process of safe and efficient personalized dosing regimens. However, reliable DGI model input data from in vivo and in vitro measurements are crucial. For this, collaboration among pharmacometricians, laboratory scientists and clinicians is important to provide homogeneous datasets and unambiguous model parameters. For a broad adaptation of validated DGI models in clinical practice, interdisciplinary cooperation should be promoted and qualification toolchains must be established.

The Influence of CYP2D6 and CYP2C19 Genetic Variation on Diabetes Mellitus Risk in People Taking Antidepressants and Antipsychotics

CYP2D6 and CYP2C19 enzymes are essential in the metabolism of antidepressants and antipsychotics. Genetic variation in these genes may increase risk of adverse drug reactions. Antidepressants and antipsychotics have previously been associated with risk of diabetes. We examined whether individual genetic differences in CYP2D6 and CYP2C19 contribute to these effects. We identified 31,579 individuals taking antidepressants and 2699 taking antipsychotics within UK Biobank. Participants were classified as poor, intermediate, or normal metabolizers of CYP2D6, and as poor, intermediate, normal, rapid, or ultra-rapid metabolizers of CYP2C19. Risk of diabetes mellitus represented by HbA1c level was examined in relation to the metabolic phenotypes. CYP2D6 poor metabolizers taking paroxetine had higher Hb1Ac than normal metabolizers (mean difference: 2.29 mmol/mol; p < 0.001). Among participants with diabetes who were taking venlafaxine, CYP2D6 poor metabolizers had higher HbA1c levels compared to normal metabolizers (mean differences: 10.15 mmol/mol; p < 0.001. Among participants with diabetes who were taking fluoxetine, CYP2D6 intermediate metabolizers and decreased HbA1c, compared to normal metabolizers (mean difference -7.74 mmol/mol; p = 0.017). We did not observe any relationship between CYP2D6 or CYP2C19 metabolic status and HbA1c levels in participants taking antipsychotic medication. Our results indicate that the impact of genetic variation in CYP2D6 differs depending on diabetes status. Although our findings support existing clinical guidelines, further research is essential to inform pharmacogenetic testing for people taking antidepressants and antipsychotics.

Model-based assessment of the safety of community interventions with primaquine in sub-Saharan Africa

BACKGROUND

Single low-dose primaquine (SLD-PQ) is recommended in combination with artemisinin-based combination therapy to reduce Plasmodium falciparum transmission in areas threatened by artemisinin resistance or aiming for malaria elimination. SLD-PQ may be beneficial in mass drug administration (MDA) campaigns to prevent malaria transmission but uptake is limited by concerns of hemolysis in glucose-6-phosphate dehydrogenase (G6PD)-deficient individuals. The aim of this study was to improve the evidence on the safety of MDA with SLD-PQ in a sub-Saharan African setting.

METHODS

A nonlinear mixed-effects model describing the pharmacokinetics and treatment-induced hemolysis of primaquine was developed using data from an adult (n = 16, G6PD deficient) and pediatric study (n = 38, G6PD normal). The relationship between primaquine pharmacokinetics and hemolysis was modeled using an established erythrocyte lifespan model. The safety of MDA with SLD-PQ was explored through Monte Carlo simulations for SLD-PQ at 0.25 or 0.4 mg/kg using baseline data from a Tanzanian setting with detailed information on hemoglobin concentrations and G6PD status.

RESULTS

The predicted reduction in hemoglobin levels following SLD-PQ was small and returned to pre-treatment levels after 25 days. G6PD deficiency (African A- variant) was associated with a 2.5-fold (95% CI 1.2-8.2) larger reduction in hemoglobin levels. In the Tanzanian setting where 43% of the population had at least mild anemia (hemoglobin < 11-13 g/dl depending on age and sex) and 2.73% had severe anemia (hemoglobin < 7-8 g/dl depending on age and sex), an additional 3.7% and 6.0% of the population were predicted to develop at least mild anemia and 0.25% and 0.41% to develop severe anemia after 0.25 and 0.4 mg/kg SLD-PQ, respectively. Children < 5 years of age and women ≥ 15 years of age were found to have a higher chance to have low pre-treatment hemoglobin.

CONCLUSIONS

This study supports the feasibility of MDA with SLD-PQ in a sub-Saharan African setting by predicting small and transient reductions in hemoglobin levels. In a setting where a substantial proportion of the population had low hemoglobin concentrations, our simulations suggest treatment with SLD-PQ would result in small increases in the prevalence of anemia which would most likely be transient.

Oral Yohimbine as a New Probe Drug to Predict CYP2D6 Activity: Results of a Fixed-Sequence Phase I Trial

Objective

Yohimbine pharmacokinetics were determined after oral administration of a single oral dose of yohimbine 5 mg and a microdose of yohimbine 50 µg in relation to different cytochrome P450 (CYP) 2D6 genotypes. The CYP2D6 inhibitor paroxetine was used to investigate the influence on yohimbine pharmacokinetics. Microdosed midazolam was applied to evaluate a possible impact of yohimbine on CYP3A activity and the possibility of combining microdosed yohimbine and midazolam to simultaneously determine CYP2D6 and CYP3A activity.

Methods

In a fixed-sequence clinical trial, 16 healthy volunteers with a known CYP2D6 genotype [extensive (10), intermediate (2) and poor (4) metaboliser] received an oral dose of yohimbine 50 µg, yohimbine 5 mg at baseline and during paroxetine as a CYP2D6 inhibitor. Midazolam (30 µg) was co-administered to determine CYP3A activity at each occasion. Plasma concentrations of yohimbine, its main metabolite 11-OH-yohimbine, midazolam and paroxetine were quantified using validated liquid chromatography-tandem mass spectrometry assays.

Results

Pharmacokinetics of yohimbine were highly variable and a CYP2D6 genotype dependent clearance was observed. After yohimbine 5 mg, the clearance ranged from 25.3 to 15,864 mL/min and after yohimbine 50 µg, the clearance ranged from 39.6 to 38,822 mL/min. A more than fivefold reduction in clearance was caused by paroxetine in CYP2D6 extensive metabolisers, while the clearance in poor metabolisers was not affected. Yohimbine did not alter CYP3A activity as measured by microdosed midazolam.

Conclusions

The pharmacokinetics of yohimbine were highly correlated with CYP2D6, which was further supported by the clearance inhibition caused by the CYP2D6 inhibitor paroxetine. With these data, yohimbine is proposed to be a suitable probe drug to predict CYP2D6 activity. In addition, the microdose can be used in combination with microdosed midazolam to simultaneously evaluate CYP2D6 and CYP3A activity without any interaction between the probe drugs and because the microdoses exert no pharmacological effects.

Clinical Trial Registration

EudraCT2017-001801-34.

Electronic supplementary material

The online version of this article (10.1007/s40262-020-00862-6) contains supplementary material, which is available to authorized users.

Comparison of Multidrug Use in the General Population and among Persons with Diabetes in Denmark for Drugs Having Pharmacogenomics (PGx) Based Dosing Guidelines

Background: This study measures the use of drugs within the therapeutic areas of antithrombotic agents (B01), the cardiovascular system (C), analgesics (N02), psycholeptics (N05), and psychoanaleptics (N06) among the general population (GP) in comparison to persons with diabetes in Denmark. The study focuses on drugs having pharmacogenomics (PGx) based dosing guidelines for CYP2D6, CYP2C19, and SLCO1B1 to explore the potential of applying PGx-based decision-making into clinical practice taking drug–drug interactions (DDI) and drug–gene interactions (DGI) into account. Methods: This study is cross-sectional, using The Danish Register of Medicinal Product Statistics as the source to retrieve drug consumption data. Results: The prevalence of use in particular for antithrombotic agents (B01) and cardiovascular drugs (C) increases significantly by 4 to 6 times for diabetic users compared to the GP, whereas the increase for analgesics (N02), psycoleptics, and psychoanaleptics (N06) was somewhat less (2–3 times). The five most used PGx drugs, both in the GP and among persons with diabetes, were pantoprazole, simvastatin, atorvastatin, metoprolol, and tramadol. The prevalence of use for persons with diabetes compared to the GP (prevalence ratio) increased by an average factor of 2.9 for all PGx drugs measured. In addition, the prevalence of use of combinations of PGx drugs was 4.6 times higher for persons with diabetes compared to GP. In conclusion, the findings of this study clearly show that a large fraction of persons with diabetes are exposed to drugs or drug combinations for which there exist PGx-based dosing guidelines related to CYP2D6, CYP2C19, and SLCO1B1. This further supports the notion of accessing and accounting for not only DDI but also DGI and phenoconversion in clinical decision-making, with a particular focus on persons with diabetes.

Analysis Approaches to Identify Pharmacogenetic Associations With Pharmacodynamics

Pharmacogenetics (PGx) seeks to enable selection of the right dose of the right drug for each patient to optimize therapeutic outcomes. Most PGx focuses on pharmacokinetics (PKs), due to our relatively advanced understanding of the genes involved in PKs and the causative effects of variants in those genes. Genetic variants can also affect pharmacodynamics (PDs), but relatively few PGx-PD associations have been identified. This is partially due to a more limited understanding of the relevant genes and the consequences of genetic variation, but is also due in part to the potential confounding of PK variability in assessments of clinical outcomes that have a contribution from both PKs and PDs. For example, it is challenging to confirm the effect of mu opioid receptor (OPRM1) genetic variation on opioid response due to the contribution of CYP2D6 genotype to bioactivation of some opioid drugs (i.e., codeine and tramadol). The objectives of this mini-review are to describe several recent efforts to discover and validate PGx-PD that disentangle the influence of PK variability and propose potential approaches that could be used in future PGx-PD analyses. We use the effect of OPRM1 genetics on opioid response to illustrate how these analyses could be conducted and conclude by discussing how PGx-PD could be translated into clinical practice to improve therapeutic outcomes.

Impact of CYP2D6 Pharmacogenomic Status on Pain Control Among Opioid-Treated Oncology Patients

BACKGROUND

Several opioids have pharmacogenomic associations impacting analgesic efficacy. However, germline pharmacogenomic testing is not routinely incorporated into supportive oncology. We hypothesized that CYP2D6 profiling would correlate with opioid prescribing and hospitalizations.

MATERIALS AND METHODS

We analyzed 61,572 adult oncology patients from 2012 to 2018 for opioid exposures. CYP2D6 metabolizer phenotype (ultra-rapid [UM], normal metabolizer [NM], intermediate [IM], or poor [PM]), the latter two of which may cause inefficacy of codeine, tramadol, and standard-dose hydrocodone, was determined for patients genotyped for reasons unrelated to pain. The primary endpoint was number of opioid medications received during longitudinal care (IM/PMs vs. NMs). Secondary endpoint was likelihood of pain-related hospital encounters.

RESULTS

Most patients with cancer (n = 34,675, 56%) received multiple opioids (average 2.8 ± 1.6/patient). Hydrocodone was most commonly prescribed (62%), followed by tramadol, oxycodone, and codeine. In the CYP2D6 genotyped cohort (n = 105), IM/PMs received a similar number of opioids (3.4 ± 1.4) as NMs (3.3 ± 1.9). However, IM/PMs were significantly more likely to experience pain-related hospital encounters compared with NMs, independent of other variables (odds ratio [OR] = 5.4; 95% confidence interval [CI], 1.2-23.6; p = .03). IM/PMs were also more likely to be treated with later-line opioids that do not require CYP2D6 metabolism, such as morphine and hydromorphone (OR = 3.3; 95% CI, 1.1-9.8; p = .03).

CONCLUSION

CYP2D6 genotype may identify patients with cancer at increased risk for inadequate analgesia when treated with typical first-line opioids like codeine, tramadol, or standard-dose hydrocodone. Palliative care considerations are an integral part of optimal oncology care, and these findings justify prospective evaluation of preemptive genotyping as a strategy to improve oncology pain management.

IMPLICATIONS FOR PRACTICE

Genomic variation in metabolic enzymes can predispose individuals to inefficacy when receiving opioid pain medications. Patients with intermediate and/or poor CYP2D6 metabolizer status do not adequately convert codeine, tramadol, and hydrocodone into active compounds, with resulting increased risk of inadequate analgesia. This study showed that patients with cancer frequently receive CYP2D6-dependent opioids. However, patients with CYP2D6 intermediate and poor metabolizer status had increased numbers of pain-related hospitalizations and more frequently required the potent non-CYP2D6 opioids morphine and hydromorphone. This may reflect inadequate initial analgesia with the common "first-line" CYP2D6-metabolized opioids. Preemptive genotyping to guide opioid prescribing during cancer care may improve pain-related patient outcomes.

Ontogeny of Drug-Metabolizing Enzymes

Almost 50% of prescription drugs lack age-appropriate dosing guidelines and therefore are used "off-label." Only ~10% drugs prescribed to neonates and infants have been studied for safety or efficacy. Immaturity of drug metabolism in children is often associated with drug toxicity. This chapter summarizes data on the ontogeny of major human metabolizing enzymes involved in oxidation, reduction, hydrolysis, and conjugation of drugs. The ontogeny data of individual drug-metabolizing enzymes are important for accurate prediction of drug pharmacokinetics and toxicity in children. This information is critical for designing clinical studies to appropriately test pharmacological hypotheses and develop safer pediatric drugs, and to replace the long-standing practice of body weight- or surface area-normalized drug dosing. The application of ontogeny data in physiologically based pharmacokinetic model and regulatory submission are discussed.

Increased primaquine total dose prevents Plasmodium vivax relapses in patients with impaired CYP2D6 activity: report of three cases

Background

The relapsing nature of Plasmodium vivax infection is a major barrier to its control and elimination. Factors such as adequate dosing, adherence, drug quality, and pharmacogenetics can impact the effectiveness of radical cure of P. vivax and need to be adequately evaluated. CYP2D6 pathway mediates the activation of primaquine (primaquine) into an active metabolite(s) in hepatocytes, and impaired activity has been linked to a higher risk of relapse.

Cases presentation

Three patients diagnosed with P. vivax malaria presented repeated relapses after being initially treated with chloroquine (25 mg/kg) and primaquine (3.5 mg/kg in 14 days) at a non-endemic travel clinic. Recurring episodes were subsequently treated with a higher dose of primaquine (7 mg/kg in 14 days), which prevented further relapses in two patients. However, one patient still presented two episodes after a higher primaquine dose and was prescribed 300 mg of chloroquine weekly to prevent further episodes. Impaired CYP2D6 function was observed in all of them.

Conclusion

Lack of response to primaquine was associated with impaired CYP2D6 activity in three patients presenting multiple relapses followed in a non-endemic setting. Higher primaquine dosage was safe and effectively prevented relapses in two patients and should be further investigated as an option in Latin America. It is crucial to investigate the factors associated with unsuccessful radical cures and alternative therapeutic options.

METABOLOMICS REVEALS BIOMARKERS IN HUMAN URINE AND PLASMA TO PREDICT CYP2D6 ACTIVITY

Background and Purpose

Individualized assessment of cytochrome P450 2D6 (CYP2D6) activity is usually performed through phenotyping following administration of a probe drug to measure the enzyme's activity. To avoid any iatrogenic harm (allergic drug reaction, dosing error) related to the probe drug, the development of non‐burdensome tools for real‐time phenotyping of CYP2D6 could significantly contribute to precision medicine. This study focuses on the identification of markers of the CYP2D6 enzyme in human biofluids using an LC‐high‐resolution mass spectrometry‐based metabolomic approach.

Experimental Approach

Plasma and urine samples from healthy volunteers were analysed before and after intake of a daily dose of paroxetine 20 mg over 7 days. CYP2D6 genotyping and phenotyping, using single oral dose of dextromethorphan 5 mg, were also performed in all participants.

Key Results

We report four metabolites of solanidine and two unknown compounds as possible novel CYP2D6 markers. Mean relative intensities of these features were significantly reduced during the inhibition session compared with the control session (n = 37). Semi‐quantitative analysis showed that the largest decrease (−85%) was observed for the ion m/z 432.3108 normalized to solanidine (m/z 398.3417). Mean relative intensities of these ions were significantly higher in the CYP2D6 normal–ultrarapid metabolizer group (n = 37) compared with the poor metabolizer group (n = 6). Solanidine intensity was more than 15 times higher in CYP2D6‐deficient individuals compared with other volunteers.

Conclusion and Implications

The applied untargeted metabolomic strategy identified potential novel markers capable of semi‐quantitatively predicting CYP2D6 activity, a promising discovery for personalized medicine.

Toward predicting CYP2D6-mediated variable drug response from CYP2D6 gene sequencing data

Pharmacogenomics is a key component of personalized medicine that promises safer and more effective drug treatment by individualizing drug choice and dose based on genetic profiles. In clinical practice, genetic biomarkers are used to categorize patients into *-alleles to predict CYP450 enzyme activity and adjust drug dosages accordingly. However, this approach leaves a large part of variability in drug response unexplained. Here, we present a proof-of-concept approach that uses continuous-scale (instead of categorical) assignments to predict enzyme activity. We used full CYP2D6 gene sequences obtained with long-read amplicon-based sequencing and cytochrome P450 (CYP) 2D6-mediated tamoxifen metabolism data from a prospective study of 561 patients with breast cancer to train a neural network. The model explained 79% of interindividual variability in CYP2D6 activity compared to 54% with the conventional *-allele approach, assigned enzyme activities to known alleles with previously reported effects, and predicted the activity of previously uncharacterized combinations of variants. The results were replicated in an independent cohort of tamoxifen-treated patients (model R ² adjusted = 0.66 versus *-allele R ² adjusted = 0.35) and a cohort of patients treated with the CYP2D6 substrate venlafaxine (model R ² adjusted = 0.64 versus *-allele R ² adjusted = 0.55). Human embryonic kidney cells were used to confirm the effect of five genetic variants on metabolism of the CYP2D6 substrate bufuralol in vitro. These results demonstrate the advantage of a continuous scale and a completely phased genotype for prediction of CYP2D6 enzyme activity and could potentially enable more accurate prediction of individual drug response.

Validation of a Population Pharmacokinetic Model of Vortioxetine Using Therapeutic Drug Monitoring Data

Introduction

Vortioxetine is an antidepressant primarily metabolized by the polymorphic enzyme cytochrome P450 (CYP) 2D6. A population pharmacokinetic (popPK) model of vortioxetine and its CYP2D6-dependent metabolite was recently published.

Objective

The aim of the current study was to assess the predictive performance of the popPK model using vortioxetine concentration measurements from a clinical setting. Furthermore, the study aimed to evaluate the ability of different CYP2D6 phenotype classification systems to provide accurate concentration predictions.

Methods

Overall, 1388 patients receiving vortioxetine treatment were identified from a therapeutic drug monitoring (TDM) database in Oslo, Norway; 334 CYP2D6-genotyped patients with 502 serum concentrations of vortioxetine, analysed by a validated ultra-high-performance liquid chromatography–high-resolution mass spectrometry (UHPLC-HRMS) method, were retrospectively included. The performance of the vortioxetine popPK model was tested on the clinical data from the TDM database.

Results

Overall, the model had a good ability to predict vortioxetine concentrations measured in clinical practice, with a slight tendency to overpredict concentrations. Using simulation-based diagnostics, 76% of the prediction-corrected TDM concentrations were within the 90% prediction interval based on 1000 simulated data sets. Prediction-based diagnostics showed the best performance for CYP2D6 poor and ultrarapid metabolizers, with a median prediction error (MDPE) of 12% and 23%, respectively, while the poorest performance was observed for normal metabolizers, with an MDPE of 66%. In the comparison of different CYP2D6 phenotype classification systems, the use of differentiated activity scores for decreased function alleles did not improve the concentration predictions. Grouping the CYP2D6 genotypes into the four conventional phenotype groups provided predictions closest to the TDM measured concentrations.

Conclusion

TDM data provide a unique insight into real-world clinical practice with vortioxetine. The tendency of the popPK model to overpredict vortioxetine concentrations measured in TDM may be attributed to several factors, including poor treatment compliance for some patients and, to a lesser extent, lack of information on patient characteristics and misspecified CYP2D6 alleles. To optimize personalized therapy with vortioxetine, real-world clinical data sets originating from different ethnicities need to be studied in the future.

Supplementary Information

The online version contains supplementary material available at 10.1007/s40262-021-01029-7.

Evaluation of the CYP2D6 Haplotype Activity Scores Based on Metabolic Ratios of 4,700 Patients Treated With Three Different CYP2D6 Substrates

The metabolic activity of the polymorphic CYP2D6 enzyme is dependent on the CYP2D6 genotype; however, the guidelines for translating the genotype into phenotype, which are of relevance for adequate drug dose personalization, are ambiguous. In the present study, retrospective therapeutic drug monitoring data from 4,700 CYP2D6 genotyped patients treated with risperidone, venlafaxine, and/or aripiprazole were analyzed to quantify the effect of CYP2D6 genotype on the CYP2D6 metabolic activities, as measured by metabolic ratios of these substrates. The patients were categorized into diplotypes based on the presence of normal function (CYP2D6Norm), nonfunctional (CYP2D6Nonf), and decreased function (CYP2D6Decr; i.e., CYP2D6*9, CYP2D6*10, and CYP2D6*41) CYP2D6 haplotypes. Significant correlations between the metabolic ratios were observed in patients (n = 77-103) cotreated with risperidone and venlafaxine, risperidone and aripiprazole, or venlafaxine and aripiprazole (ρ = 0.874, 0.785, and 0.644, respectively; P < 0.001 for all). Relative metabolic CYP2D6 diplotype activity was calculated based on that the metabolic ratios, where median values for CYP2D6Nonf/Nonf and CYP2D6Norm/Norm subgroups were set to 0% and 100%, respectively. The relative CYP2D6 activities were: 7.0% for CYP2D6Nonf/*41, 16.7% for CYP2D6Nonf/*9-10, 13.2% for CYP2D6*41/*41, 24.9% for CYP2D6*41/*9-10, 33.1% for CYP2D6*9-10/*9-10, 41.3% for CYP2D6Nonf/Norm, 55.0% for CYP2D6*41/Norm, 58.9% for CYP2D6*9-10/Norm, and 149.2% for CYP2D6Norm/Normx2. Compared with the CYP2D6Norm alleles, the activity scores of CYP2D6*41 and CYP2D6*9-10 alleles were estimated to be one sixth and one third, respectively. The results of this highly powered study provide a solid basis for the translation of the CYP2D6 genotype into a drug metabolic phenotype.

CYP2D6 Reduced Function Variants and Genotype/Phenotype Translations of CYP2D6 Intermediate Metabolizers: Implications for Personalized Drug Dosing in Psychiatry

Genetic differences in cytochrome P450 (CYP)-mediated metabolism have been known for several decades. The clinically most important polymorphic CYP enzyme is CYP2D6, which plays a key role in the metabolism of many antidepressants and antipsychotics, along with a range of non-psychiatric medications. Dose individualization based on CYP2D6 genotype to improve the effect and safety of drug treatment has been an ambition for a long time. Clinical use of CYP2D6 genotyping is steadily increasing; however, for pre-emptive genotyping to be successful in predicting individual dose requirements, high precision of genotype-to-phenotype translations are required. Recently, guidelines for assigning CYP2D6 enzyme activity scores of CYP2D6 variant alleles, and subsequent diplotype-to-phenotype translations, were published by the Clinical Pharmacogenetics Implementation Consortium (CPIC) and the Dutch Pharmacogenetics Working Group. Consensus on assigning activity scores of CYP2D6 variant alleles and translating diplotype scores into CYP2D6 poor, intermediate, normal, or ultrarapid metabolizer groups were obtained by consulting 37 international experts. While assigning enzyme activities of non-functional (score 0) and fully functional (score 1) alleles are straightforward, reduced function variant alleles are more complex. In this article, we present data showing that the assigned activity scores of reduced function variant alleles in current guidelines are not of sufficient precision; especially not for CYP2D6*41, where the guideline activity score is 0.5 compared to 0.05–0.15 in pharmacogenetic studies. Due to these discrepancies, CYP2D6 genotypes with similar guidelinediplotype scores exhibit substantial differences in CYP2D6 metabolizer phenotypes. Thus, it is important that the guidelines are updated to be valid in predicting individual dose requirements of psychiatric drugs and others metabolized by CYP2D6.

Co-expression of drug metabolizing cytochrome P450 enzymes and estrogen receptor alpha (ESR1) in human liver: racial differences and the regulatory role of ESR1

OBJECTIVES

The function and expression of cytochrome P450 (CYP) drug metabolizing enzymes is highly variable, greatly affecting drug exposure, and therapeutic outcomes. The expression of these enzymes is known to be controlled by many transcription factors (TFs), including ligand-free estrogen receptor alpha (ESR1, in the absence of estrogen). However, the relationship between the expression of ESR1, other TFs, and CYP enzymes in human liver is still unclear.

METHODS

Using real-time PCR, we quantified the mRNA levels of 12 CYP enzymes and nine TFs in 246 human liver samples from European American (EA, n = 133) and African American (AA, n = 113) donors.

RESULTS

Our results showed higher expression levels of ESR1 and six CYP enzymes in EA than in AA. Partial least square regression analysis showed that ESR1 is the top-ranking TF associating with the expression of eight CYP enzymes, six of which showed racial difference in expression. Conversely, four CYP enzymes without racial difference in expression did not have ESR1 as a top-ranking TF. These results indicate that ESR1 may contribute to variation in CYP enzyme expression between these two ancestral backgrounds.

CONCLUSIONS

These results are consistent with our previous study showing ESR1 as a master regulator for the expression of several CYP enzymes. Therefore, factors affecting ESR1 expression may have broad influence on drug metabolism through altered expression of CYP enzymes.

Estimated impact of tafenoquine for Plasmodium vivax control and elimination in Brazil: A modelling study

BACKGROUND

Despite recent intensification of control measures, Plasmodium vivax poses a major challenge for malaria elimination efforts. Liver-stage hypnozoite parasites that cause relapsing infections can be cleared with primaquine; however, poor treatment adherence undermines drug effectiveness. Tafenoquine, a new single-dose treatment, offers an alternative option for preventing relapses and reducing transmission. In 2018, over 237,000 cases of malaria were reported to the Brazilian health system, of which 91.5% were due to P. vivax.

METHODS AND FINDINGS

We evaluated the impact of introducing tafenoquine into case management practices on population-level transmission dynamics using a mathematical model of P. vivax transmission. The model was calibrated to reflect the transmission dynamics of P. vivax endemic settings in Brazil in 2018, informed by nationwide malaria case reporting data. Parameters for treatment pathways with chloroquine, primaquine, and tafenoquine with glucose-6-phosphate dehydrogenase deficiency (G6PDd) testing were informed by clinical trial data and the literature. We assumed 71.3% efficacy for primaquine and tafenoquine, a 66.7% adherence rate to the 7-day primaquine regimen, a mean 5.5% G6PDd prevalence, and 8.1% low metaboliser prevalence. The introduction of tafenoquine is predicted to improve effective hypnozoite clearance among P. vivax cases and reduce population-level transmission over time, with heterogeneous levels of impact across different transmission settings. According to the model, while achieving elimination in only few settings in Brazil, tafenoquine rollout in 2021 is estimated to improve the mean effective radical cure rate from 42% (95% uncertainty interval [UI] 41%-44%) to 62% (95% UI 54%-68%) among clinical cases, leading to a predicted 38% (95% UI 7%-99%) reduction in transmission and over 214,000 cumulative averted cases between 2021 and 2025. Higher impact is predicted in settings with low transmission, low pre-existing primaquine adherence, and a high proportion of cases in working-aged males. High-transmission settings with a high proportion of cases in children would benefit from a safe high-efficacy tafenoquine dose for children. Our methodological limitations include not accounting for the role of imported cases from outside the transmission setting, relying on reported clinical cases as a measurement of community-level transmission, and implementing treatment efficacy as a binary condition.

CONCLUSIONS

In our modelling study, we predicted that, provided there is concurrent rollout of G6PDd diagnostics, tafenoquine has the potential to reduce P. vivax transmission by improving effective radical cure through increased adherence and increased protection from new infections. While tafenoquine alone may not be sufficient for P. vivax elimination, its introduction will improve case management, prevent a substantial number of cases, and bring countries closer to achieving malaria elimination goals.

The association of CYP2D6 gene polymorphisms in the full-length coding region with higher recurrence rate of vivax malaria in Yunnan Province, China

Background

Accumulating evidence suggest that compromised CYP2D6 enzyme activity caused by gene mutation could contribute to primaquine failure for the radical cure of vivax malaria. The current study aims to preliminarily reveal the association between the recurrence of vivax malaria in Yunnan Province and CYP2D6 gene mutation by analysing polymorphisms in the entire coding region of human CYP2D6 gene.

Methods

Blood samples were collected from patients with vivax malaria, who received "chloroquine and 8-day course of primaquine therapy" in Yunnan Province. The suspected relapsed cases were determined by epidemiological approaches and gene sequence alignment. PCR was conducted to amplify the CYP2D6 gene in the human genome, and the amplified products were then sequenced to compare with the non-mutation “reference” sequence, so as to ensure correct sequencing results and to determine 9 exon regions. Subsequently, the DNA sequences of 9 exons were spliced into the coding DNA sequence (CDS), which, by default, is known as maternal CDS. The paternal CDS was obtained by adjusting the bases according to the sequencing peaks. The mutation loci, haplotypes (star alleles), genotypes and odds ratios (OR) of all the CDSs were analysed.

Results

Of the119 maternal CDS chains in total with 1491 bp in length, 12 mutation sites in the 238 maternal and paternal CDS chains were detected. The c.408G > C mutation was most frequently detected in both suspected relapsed group (SR) and non-relapsed group (NR), reaching 85.2% (75/88) and 76.0% (114/150), respectively. The c.886C > T mutation was most closely related to the recurrence of vivax malaria (OR = 2.167, 95% CI 1.104–4.252, P < 0.05). Among the 23 haplotypes (Hap_1 ~ Hap_23), Hap_3 was non-mutant, and the sequence structure of Hap_9 was the most complicated one. Five star alleles, including *1, *2, *4, *10 and *39, were confirmed by comparison, and CYP2D6*10 allele accounted for the largest percentage (45.4%, 108/238). The frequency of CYP2D6*2 allele in the SR group was significantly higher than that in the NR group (Χ² = 16.177, P < 0.05). Of the defined 24 genotypes, 8 genotypes, including *4/*4, *4/*o, *2/*39, *39/*m, *39/*x, *1/*r, *1/*n, and *v/*10, were detected only in the SR group.

Conclusion

Mutation of CYP2D6*10 allele accounts for the highest proportion of vivax malaria cases in Yunnan Province. The mutations of c. 886C > T and CYP2D6*2 allele, which correspond to impaired PQ metabolizer phenotype, are most closely related to the relapse of vivax malaria. In addition, the genotype *4/*4 with null CYP2D6 enzyme function was only detected in the SR group. These results reveal the risk of defected CYP2D6 enzyme activity that diminishes the therapeutic effect of primaquine on vivax malaria.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12936-021-03685-3.

Whole Genome Interpretation for a Family of Five

Although best practices have emerged on how to analyse and interpret personal genomes, the utility of whole genome screening remains underdeveloped. A large amount of information can be gathered from various types of analyses via whole genome sequencing including pathogenicity screening, genetic risk scoring, fitness, nutrition, and pharmacogenomic analysis. We recognize different levels of confidence when assessing the validity of genetic markers and apply rigorous standards for evaluation of phenotype associations. We illustrate the application of this approach on a family of five. By applying analyses of whole genomes from different methodological perspectives, we are able to build a more comprehensive picture to assist decision making in preventative healthcare and well-being management. Our interpretation and reporting outputs provide input for a clinician to develop a healthcare plan for the individual, based on genetic and other healthcare data.

Potential implications of DMET ontogeny on the disposition of commonly prescribed drugs in neonatal and pediatric intensive care units

Introduction: Pediatric patients, especially neonates and infants, are more susceptible to adverse drug events as compared to adults. In particular, immature small molecule drug metabolism and excretion can result in higher incidences of pediatric toxicity than adults if the pediatric dose is not adjusted.Area covered: We reviewed the top 29 small molecule drugs prescribed in neonatal and pediatric intensive care units and compiled the mechanisms of their metabolism and excretion. The ontogeny of Phase I and II drug metabolizing enzymes and transporters (DMETs), particularly relevant to these drugs, are summarized. The potential effects of DMET ontogeny on the metabolism and excretion of the top pediatric drugs were predicted. The current regulatory requirements and recommendations regarding safe and effective use of drugs in children are discussed. A few representative examples of the use of ontogeny-informed physiologically based pharmacokinetic (PBPK) models are highlighted.Expert opinion: Empirical prediction of pediatric drug dosing based on body weight or body-surface area from the adult parameters can be inaccurate because DMETs are not mature in children and the age-dependent maturation of these proteins is different. Ontogeny-informed-PBPK modeling provides a better alternative to predict the pharmacokinetics of drugs in children.

StellarPGx: A Nextflow Pipeline for Calling Star Alleles in Cytochrome P450 Genes

Bioinformatics pipelines for calling star alleles (haplotypes) in cytochrome P450 (CYP) genes are important for the implementation of precision medicine. Genotyping CYP genes using high throughput sequencing data is complicated, e.g., by being highly polymorphic, not to mention the structural variations especially in CYP2D6, CYP2A6, and CYP2B6. Genome graph-based variant detection approaches have been shown to be reliable for genotyping HLA alleles. However, their application to enhancing star allele calling in CYP genes has not been extensively explored. We present StellarPGx, a Nextflow pipeline for accurately genotyping CYP genes by combining genome graph-based variant detection, read coverage information from the original reference-based alignments, and combinatorial diplotype assignments. The implementation of StellarPGx using Nextflow facilitates its portability, reproducibility, and scalability on various user platforms. StellarPGx is currently able to genotype 12 important pharmacogenes belonging to the CYP1, 2, and 3 families. For purposes of validation, we use CYP2D6 as a model gene owing to its high degree of polymorphisms (over 130 star alleles defined to date, including complex structural variants) and clinical importance. We applied StellarPGx and three existing callers to 109 whole genome sequenced samples for which the Genetic Testing Reference Material Coordination Program (GeT-RM) has recently provided consensus truth CYP2D6 diplotypes. StellarPGx had the highest CYP2D6 diplotype concordance (99%) with GeT-RM compared with Cyrius (98%), Aldy (82%), and Stargazer (84%). This exemplifies the high accuracy of StellarPGx and highlights its importance for both research and clinical pharmacogenomics applications. The StellarPGx pipeline is open-source and available from https://github.com/SBIMB/StellarPGx.