Polymorphism of Human Cytochrome P450 2D6 and Its Clinical Significance

Optimisation of pharmacotherapy in psychiatry through therapeutic drug monitoring, molecular brain imaging and pharmacogenetic tests: focus on antipsychotics

BACKGROUND

For psychotic disorders (i.e. schizophrenia), pharmacotherapy plays a key role in controlling acute and long-term symptoms. To find the optimal individual dose and dosage strategy, specialized tools are used. Three tools have been proven useful to personalize drug treatments: therapeutic drug monitoring (TDM) of drug levels, pharmacogenetic testing (PG), and molecular neuroimaging.

METHODS

In these Guidelines, we provide an in-depth review of pharmacokinetics, pharmacodynamics, and pharmacogenetics for 50 antipsychotics. Over 30 international experts in psychiatry selected studies that have measured drug concentrations in the blood (TDM), gene polymorphisms of enzymes involved in drug metabolism, or receptor/transporter occupancies in the brain (positron emission tomography (PET)).

RESULTS

Study results strongly support the use of TDM and the cytochrome P450 (CYP) genotyping and/or phenotyping to guide drug therapies. Evidence-based target ranges are available for titrating drug doses that are often supported by PET findings.

CONCLUSION

All three tools discussed in these Guidelines are essential for drug treatment. TDM goes well beyond typical indications such as unclear compliance and polypharmacy. Despite its enormous potential to optimize treatment effects, minimize side effects and ultimately reduce the global burden of diseases, personalized drug treatment has not yet become the standard of care in psychiatry.

Influence of CYP2C19 and CYP2D6 on side effects of aripiprazole and risperidone: A systematic review

Variability in hepatic cytochrome P450 (CYP) enzymes such as 2C19 and 2D6 may influence side-effect and efficacy outcomes for antipsychotics. Aripiprazole and risperidone are two commonly prescribed antipsychotics, metabolized primarily through CYP2D6. Here, we aimed to provide an overview of the effect of CYP2C19 and CYP2D6 on side-effects of aripiprazole and risperidone, and expand on existing literature by critically examining methodological issues associated with pharmacogenetic studies. A PRISMA compliant search of six electronic databases (Pubmed, PsychInfo, Embase, Central, Web of Science, and Google Scholar) identified pharmacogenetic studies on aripiprazole and risperidone. 2007 publications were first identified, of which 34 were included. Quality of literature was estimated using Newcastle-Ottowa Quality Assessment Scale (NOS) and revised Cochrane Risk of Bias tool. The average NOS score was 5.8 (range: 3-8) for risperidone literature and 5 for aripiprazole (range: 4-6). All RCTs on aripiprazole were rated as high risk of bias, and four out of six for risperidone literature. Study populations ranged from healthy volunteers to inpatient individuals in psychiatric units and included adult and pediatric samples. All n = 34 studies examined CYP2D6. Only one study genotyped for CYP2C19 and found a positive association with neurological side-effects of risperidone. Most studies did not report any relationship between CYP2D6 and any side-effect outcome. Heterogeneity between and within studies limited the ability to synthesize data and draw definitive conclusions. Studies lacked statistical power due to small sample size, selective genotyping methods, and study design. Large-scale randomized trials with multiple measurements, providing robust evidence on this topic, are suggested.

Association of CYP2D6 and CYP2C19 metabolizer status with switching and discontinuing antidepressant drugs: an exploratory study

BACKGROUND

Tailoring antidepressant drugs (AD) to patients' genetic drug-metabolism profile is promising. However, literature regarding associations of ADs' treatment effect and/or side effects with drug metabolizing genes CYP2D6 and CYP2C19 has yielded inconsistent results. Therefore, our aim was to longitudinally investigate associations between CYP2D6 (poor, intermediate, and normal) and CYP2C19 (poor, intermediate, normal, and ultrarapid) metabolizer-status, and switching/discontinuing of ADs. Next, we investigated whether the number of perceived side effects differed between metabolizer statuses.

METHODS

Data came from the multi-site naturalistic longitudinal cohort Netherlands Study of Depression and Anxiety (NESDA). We selected depression- and/or anxiety patients, who used AD at some point in the course of the 9 years follow-up period (n = 928). Medication use was followed to assess patterns of AD switching/discontinuation over time. CYP2D6 and CYP2C19 alleles were derived using genome-wide data of the NESDA samples and haplotype data from the PharmGKB database. Logistic regression analyses were conducted to investigate the association of metabolizer status with switching/discontinuing ADs. Mann-Whitney U-tests were conducted to compare the number of patient-perceived side effects between metabolizer statuses.

RESULTS

No significant associations were observed of CYP metabolizer status with switching/discontinuing ADs, nor with the number of perceived side effects.

CONCLUSIONS

We found no evidence for associations between CYP metabolizer statuses and switching/discontinuing AD, nor with side effects of ADs, suggesting that metabolizer status only plays a limited role in switching/discontinuing ADs. Additional studies with larger numbers of PM and UM patients are needed to further determine the potential added value of pharmacogenetics to guide pharmacotherapy.

Polymorphisms in drug-metabolizing genes and urinary bladder cancer susceptibility and prognosis: Possible impacts and future management

Epidemiological studies have shown the association of genetic variants with risks of occupational and environmentally induced cancers, including bladder (BC). The current review summarizes the effects of variants in genes encoding phase I and II enzymes in well-designed studies to highlight their contribution to BC susceptibility and prognosis. Polymorphisms in genes codifying drug-metabolizing proteins are of particular interest because of their involvement in the metabolism of exogenous genotoxic compounds, such as tobacco and agrochemicals. The prognosis between muscle-invasive and non-muscle-invasive diseases is very different, and it is difficult to predict which will progress worse. Web of Science, PubMed, and Medline were searched to identify studies published between January 1, 2010, and February 2023. We included 73 eligible studies, more than 300 polymorphisms, and 46 genes/loci. The most studied candidate genes/loci of phase I metabolism were CYP1B1, CYP1A1, CYP1A2, CYP3A4, CYP2D6, CYP2A6, CYP3E1, and ALDH2, and those in phase II were GSTM1, GSTT1, NAT2, GSTP1, GSTA1, GSTO1, and UGT1A1. We used the 46 genes to construct a network of proteins and to evaluate their biological functions based on the Reactome and KEGG databases. Lastly, we assessed their expression in different tissues, including normal bladder and BC samples. The drug-metabolizing pathway plays a relevant role in BC, and our review discusses a list of genes that could provide clues for further exploration of susceptibility and prognostic biomarkers.

Genotype-guided new approach for dose optimisation of hydroxychloroquine administration in Chinese patients with SLE

OBJECTIVES

The study aims to investigate the impact of gene polymorphisms on blood hydroxychloroquine (HCQ) concentrations in patients with SLE and provide guidelines for individualised care.

METHODS

489 Chinese patients with SLE taking HCQ for more than 3 months were collected in this study. The blood HCQ, desethylhydroxychloroquine (DHCQ) and desethylchloroquine concentrations were measured. The optimal blood concentration of HCQ was determined by receiver operating characteristic curve analysis. Single nucleotide polymorphisms of metabolic enzymes involved in HCQ metabolism were genotyped and the associations with treatment effects were investigated.

RESULTS

The cut-off value of HCQ was 559.67 ng/mL, with sensitivity and specificity values of 0.51 and 0.89, respectively. The TC and CC genotypes of CYP2C8 (rs7910936) were significantly related to the increase in blood HCQ concentrations, and the CYP2C8 (rs10882521) TT genotype was associated with lower blood HCQ concentrations. The DHCQ:HCQ ratio was highest in patients with the GG genotype of the CYP2D6*10 (rs1065852) polymorphism and lowest in those with the AA genotype. Patients with the CYP2C8 (rs7910936) CC genotype were more likely to achieve the optimal blood concentration (p=0.030) in HCQ 200 mg/day group and patients with the CYP2D6*10 (rs1065852) GG genotype were more likely to reach the optimal blood concentration (p=0.049) in 400 mg/day group.

CONCLUSIONS

Our results suggest that the optimal blood concentration of HCQ measured approximately 12-18 hours after the last dosage may be between 500 and 600 ng/mL in Chinese patients with SLE. The observed variations in HCQ concentrations between individuals can potentially be attributed to genetic polymorphisms in CYP2D6*10 (rs1065852) and CYP2C8 (rs7910936 and rs10882521). Genotypical testing of patients and regular monitoring of blood levels are recommended for optimising HCQ dosage management in Chinese patients with SLE.

TRIAL REGISTRATION NUMBER

ChiCTR2300070628.

Effect of apatinib on the pharmacokinetics of tramadol and O-desmethyltramadol in rats

Since the combination of anticancer drugs and opioids is very common, apatinib and tramadol are likely to be used in combination clinically. This study evaluated the effects of apatinib on the pharmacokinetics of tramadol and its main metabolite O-desmethyltramadol in Sprague-Dawley (SD) rats and the inhibitory effects of apatinib on tramadol in rat liver microsomes (RLMs), human liver microsomes (HLMs) and recombinant human CYP2D6.1. The samples were determined by ultra-performance liquid chromatography-tandem mass spectrometry (UPLC-MS/MS). The in vivo results showed that compared with the control group, apatinib increased the AUC(0-t), AUC(0-∞) and Cmax values of tramadol and O-desmethyltramadol, and decreased the values of VZ/F and CLz/F. In addition, the MRT(0-t), MRT(0-∞) values of O-desmethyltramadol were increased. In vitro, apatinib inhibited the metabolism of tramadol by a mixed way with IC50 of 1.927 µM in RLMs, 2.039 µM in HLMs and 15.32 µM in CYP2D6.1. In summary, according to our findings, apatinib has a strong in vitro inhibitory effect on tramadol, and apatinib can increase the analgesic effect of tramadol and O-desmethyltramadol in rats.

Multivariate QSAR, similarity search and ADMET studies based in a set of methylamine derivatives described as dopamine transporter inhibitors

The dopamine transporter (DAT), responsible for the regulation of dopaminergic neurotransmission, is implicated in the etiology of several neuropsychiatric disorders which, in turn, have contributed to high rates of disability and numerous deaths in recent years, significantly impacting the global health system. Although the research for new drugs for the treatment of neuropsychiatric disorders has evolved in recent years, the availability of DAT-selective drugs that do not generate the same psychostimulant effects observed in drugs of abuse remains scarce. Therefore, we performed a QSAR study based on a dataset of 36 methylamine derivatives described as DAT inhibitors. The model was obtained based only in descriptors derived from 2D structures, and it was validated and generated satisfactory results considering the metrics used for internal and external validation. Subsequently, a virtual screening step also based on 2D similarity was performed, where it was possible to identify a total of 1157 compounds. After a series of reductions of the set using toxicity filters, applicability domain evaluation, and pharmacokinetic properties in silico assessment, seven hit compounds were selected as the most promising to be used, in future studies, as new scaffolds for the development of new DAT inhibitors.

Tramadol poisoning and its management and complications: a scoping review

Tramadol (TR) is an opioid agonist (μ-opioid receptors) that also affects the noradrenergic and serotonergic systems. TR is a synthetic analgesic substance with opioid agonist properties that has been approved for pain management by affecting the central nervous system.

OBJECTIVE

The primary objective of this scoping review was to determine the poisoning consequences of TR and its associated symptom management techniques, as well as its effect on opioid and non-opioid receptors, metabolism, and complications on various organ systems.

DISCUSSION

Typically, acute overdose of TR is not considered life threatening, and most fatalities are attributed to polysubstance overdose. TR can cause drowsiness, headaches, seizures, respiratory depression, low heart rate, coma, and even death. Additionally, the prolonged use of TR can lead to addiction, with withdrawal resulting in both common and atypical symptoms.

CONCLUSIONS

Therefore, it is recommended that patients with TR poisoning receive close monitoring of their cardiovascular system, along with a comprehensive management plan for their levels of consciousness and respiratory function.

How to Choose the Optimal Starting Dose of Clomiphene Citrate (50 or 100 mg per Day) for a First Cycle of Ovulation Induction in Anovulatory PCOS Women?

RESEARCH QUESTION

Clomiphene citrate (CC) is one of the first-line treatments for ovulation induction in women with anovulatory polycystic ovary syndrome (PCOS). However, nearly 1 out of 2 women is resistant to 50 mg/day of CC. The objective of this study is to investigate the clinical, biological, and/or ultrasound factors that may predict the resistance to 50 mg/day of CC in the first cycle of treatment in women with anovulatory PCOS. This would make it possible to identify PCOS patients to whom the dose of 100 mg/day would be offered as of the first cycle.

DESIGN

A retrospective and monocentric study was conducted on 283 women with anovulatory PCOS who required the use of ovulation induction with CC (903 cycles).

RESULTS

During the first cycle of treatment, 104 patients (36.8%) were resistant to 50 mg/day of CC. Univariate regression analysis showed that patients who resisted 50 mg/day of CC had significantly higher BMI, waist circumference, serum levels of AMH, total testosterone, Δ4-androstenedione, 17-OHP, and insulin (p < 0.05), compared to patients ovulating with this dose. Serum levels of SHBG were significantly lower in patients resistant to 50 mg/day (p < 0.05). After multivariate analysis, only AMH and SHBG remained statistically significant (p = 0.01 and p = 0.001, respectively). However, areas under the ROC curves were weak (0.59 and 0.68, respectively).

CONCLUSION

AMH and SHBG are the only two parameters significantly associated with the risk of resistance to 50 mg/day of CC. However, no satisfactory thresholds have been established to predict resistance to 50 mg CC.

Case report: Use of therapeutic drug monitoring and pharmacogenetic testing as opportunities to individualize care in a case of flecainide toxicity after fetal supraventricular tachycardia

Flecainide is a class IC antiarrhythmic utilized in prophylaxis of refractory paroxysmal supraventricular tachycardias in pediatric populations. Despite being a highly effective agent, its narrow therapeutic index increases the risk of toxicity and proarrhythmic events, including wide-complex tachycardia. In the absence of direct plasma sampling in the fetus to quantitate flecainide systemic concentrations, clinicians typically make drug dosing decisions from maternal plasma concentrations and QRS duration on maternal ECGs. There remains a paucity of standard guidelines and data to inform the timing and frequency of the aforementioned test in pregnancy and timing of flecainide discontinuation prior to childbirth. Flecainide primarily undergoes metabolism via cytochrome P450 (CYP). Given the variance of CYP-mediated metabolism at the level of the individual patient, pharmacogenomics can be considered in patients who present with flecainide toxicity to determine the maternal vs. fetal factors as an etiology for the event. Finally, pharmacogenetic testing can be utilized as an adjunct to guide flecainide dosing decisions, but must be done with caution in neonates <2 weeks of age. This case report highlights utilization of pharmacogenomic testing and therapeutic drug monitoring as adjuncts to guide therapy for a newborn with refractory supraventricular tachycardia, who experienced flecainide toxicity immediately post-partum and was trialed unsuccessfully on multiple alternative antiarrhythmics without rhythm control.

A Novel UHPLC-MS/MS Method for the Quantification of Seven Opioids in Different Human Tissues

BACKGROUND

Opioids are considered the cornerstone of pain management: they show good efficacy as a first-line therapy for moderate to severe cancer pain. Since pharmacokinetic/pharmacodynamic information about the tissue-specific effect and toxicity of opioids is still scarce, their quantification in post-mortem autoptic specimens could give interesting insights.

METHODS

We describe an ultra-high-performance liquid chromatography coupled with tandem mass spectrometry method for the simultaneous quantification of methadone, morphine, oxycodone, hydrocodone, oxymorphone, hydromorphone and fentanyl in several tissues: liver, brain, kidney, abdominal adipose tissue, lung and blood plasma. The presented method has been applied on 28 autoptic samples from different organs obtained from four deceased PLWH who used opioids for palliative care during terminal disease.

RESULTS

Sample preparation was based on tissue weighing, disruption, sonication with drug extraction medium and a protein precipitation protocol. The extracts were then dried, reconstituted and injected onto the LX50 QSight 220 (Perkin Elmer, Milan, Italy) system. Separation was obtained by a 7 min gradient run at 40 °C with a Kinetex Biphenyl 2.6 µm, 2.1 × 100 mm. Concerning the analyzed samples, higher opioids concentrations were observed in tissues than in plasma. Particularly, O-MOR and O-COD showed higher concentrations in kidney and liver than other tissues (>15-20 times greater) and blood plasma (>100 times greater).

CONCLUSIONS

Results in terms of linearity, accuracy, precision, recovery and matrix effect fitted the recommendations of FDA and EMA guidelines, and the sensitivity was high enough to allow successful application on human autoptic specimens from an ethically approved clinical study, confirming its eligibility for post-mortem pharmacological/toxicological studies.

Recombinant Technologies Facilitate Drug Metabolism, Pharmacokinetics, and General Biomedical Research

The development of safe and effective medications requires a profound understanding of their pharmacokinetic (PK) and pharmacodynamic properties. PK studies have been built through investigation of enzymes and transporters that drive drug absorption, distribution, metabolism, and excretion (ADME). Like many other disciplines, the study of ADME gene products and their functions has been revolutionized through the invention and widespread adoption of recombinant DNA technologies. Recombinant DNA technologies use expression vectors such as plasmids to achieve heterologous expression of a desired transgene in a specified host organism. This has enabled the purification of recombinant ADME gene products for functional and structural characterization, allowing investigators to elucidate their roles in drug metabolism and disposition. This strategy has also been used to offer recombinant or bioengineered RNA (BioRNA) agents to investigate the posttranscriptional regulation of ADME genes. Conventional research with small noncoding RNAs such as microRNAs (miRNAs) and small interfering RNAs has been dependent on synthetic RNA analogs that are known to carry a range of chemical modifications expected to improve stability and PK properties. Indeed, a novel transfer RNA fused pre-miRNA carrier-based bioengineering platform technology has been established to offer consistent and high-yield production of unparalleled BioRNA molecules from Escherichia coli fermentation. These BioRNAs are produced and processed inside living cells to better recapitulate the properties of natural RNAs, representing superior research tools to investigate regulatory mechanisms behind ADME. SIGNIFICANCE STATEMENT: This review article summarizes recombinant DNA technologies that have been an incredible boon in the study of drug metabolism and PK, providing investigators with powerful tools to express nearly any ADME gene products for functional and structural studies. It further overviews novel recombinant RNA technologies and discusses the utilities of bioengineered RNA agents for the investigation of ADME gene regulation and general biomedical research.

The pharmacogenomic landscape of an Indigenous Australian population

Background: Population genomic studies of individuals of Indigenous ancestry have been extremely limited comprising <0.5% of participants in international genetic databases and genome-wide association studies, contributing to a "genomic gap" that limits their access to personalised medicine. While Indigenous Australians face a high burden of chronic disease and associated medication exposure, corresponding genomic and drug safety datasets are sorely lacking. Methods: To address this, we conducted a pharmacogenomic study of almost 500 individuals from a founder Indigenous Tiwi population. Whole genome sequencing was performed using short-read Illumina Novaseq6000 technology. We characterised the pharmacogenomics (PGx) landscape of this population by analysing sequencing results and associated pharmacological treatment data. Results: We observed that every individual in the cohort carry at least one actionable genotype and 77% of them carry at least three clinically actionable genotypes across 19 pharmacogenes. Overall, 41% of the Tiwi cohort were predicted to exhibit impaired CYP2D6 metabolism, with this frequency being much higher than that for other global populations. Over half of the population predicted an impaired CYP2C9, CYP2C19, and CYP2B6 metabolism with implications for the processing of commonly used analgesics, statins, anticoagulants, antiretrovirals, antidepressants, and antipsychotics. Moreover, we identified 31 potentially actionable novel variants within Very Important Pharmacogenes (VIPs), five of which were common among the Tiwi. We further detected important clinical implications for the drugs involved with cancer pharmacogenomics such as thiopurines and tamoxifen, immunosuppressants like tacrolimus and certain antivirals used in the hepatitis C treatment due to potential differences in their metabolic processing. Conclusion: The pharmacogenomic profiles generated in our study demonstrate the utility of pre-emptive PGx testing and have the potential to help guide the development and application of precision therapeutic strategies tailored to Tiwi Indigenous patients. Our research provides valuable insights on pre-emptive PGx testing and the feasibility of its use in ancestrally diverse populations, emphasizing the need for increased diversity and inclusivity in PGx investigations.

Effect of pharmacogenomics testing guiding on clinical outcomes in major depressive disorder: a systematic review and meta-analysis of RCT

BACKGROUND

Pharmacogenomic testing guided treatment have been developed to guide drug selection or conversion in major depressive disorder patients. Whether patients benefit from pharmacogenetic testing remains unclear. We aim to evaluates the effect of pharmacogenomic testing guiding on clinical outcomes of major depressive disorder.

METHODS

Pubmed, Embase, and Cochrane Library of Clinical Trials were searched from inception until August 2022. Key terms included pharmacogenomic and antidepressive. Odds ratios (RR) with 95% confidence intervals (95%CIs) were calculated using fixed-effects model for low or moderate heterogeneity or random-effects model for high heterogeneity.

RESULTS

Eleven studies (5347 patients) were included. Compared with usual group, pharmacogenomic testing guided group was associated with an increased response rate at week 8 (OR 1.32, 95%CI 1.15-1.53, 8 studies, 4328 participants) and week 12 (OR 1.36, 95%CI 1.15-1.62, 4 studies, 2814 participants). Similarly, guided group was associated with an increased rate of remission at week 8 (OR 1.58, 95%CI 1.31-1.92, 8 studies, 3971 participants) and week 12 (OR 2.23, 95%CI 1.23-4.04, 5 studies, 2664 participants). However, no significant differences were found between the two groups in response rate at week 4 (OR 1.12, 95%CI 0.89-1.41, 2 studies, 2261 participants) and week 24 (OR 1.16, 95%CI 0.96-1.41, 2 studies, 2252 participants), and remission rate at week 4 (OR 1.26, 95%CI 0.93-1.72, 2 studies, 2261 participants) and week 24 (OR 1.06, 95%CI 0.83-1.34, 2 studies, 2252 participants). Medication congruence in 30 days was significantly reduced in the pharmacogenomic guided group compared with the usual care group (OR 2.07, 95%CI 1.69-2.54, 3 studies, 2862 participants). We found significant differences between subgroups of target population in response and remission rate.

CONCLUSION

Patients with major depressive disorder may benefit from pharmacogenomic testing guided treatment by achieving target response and remission rates more quickly.

Allelic and Genotype Frequencies of CYP2B6∗2 (64C > T) and CYP2B6∗3 (777C > A) in Three Dominant Ethnicities of the Iranian Population

Background

Cytochrome P450 complex plays a key role in drug metabolism. CYP2B6 has an essential part in Cytochrome P450 complex metabolism. This study aims to determine the allelic distribution of CYP2B6^∗2 and CYP2B6^∗3 in three main Iranian ethnicities: Fars, Turk, and Kurd.

Methods

The study was conducted on 174 unrelated healthy volunteers from three main Iranian ethnicities. After DNA extraction from peripheral blood samples, genotyping of CYP2B6^∗2 and ^∗3 was performed using tetra ARMS and ARMS PCR, respectively.

Results

The average age of 174 cases was 40.69 ± 11.87 (mean ± SD) and 39.06 ± 11.63 (mean ± SD) for males and females. In the CYP2B6^∗2 variant, the genotyping frequency of wild type (C/C), heterozygous (C/T), and homozygous mutant (T/T) was 8.7%, 86%, and 5.2%, respectively. The CYP2B6^∗2 (c.64C > T) allele frequency was 48.2% (95% CI: (37.8-58.6)). In the CYP2B6^∗3 variant, the frequency of wild type (C/C), heterozygous (C/T), and homozygous mutant (T/T) was 75.3%, 11%, and 13.6%, respectively. The CYP2B6^∗3 (c.777C > A) allelic frequency was 19.1% (95% CI: (17.5-20.7)).

Conclusion

Allelic distribution in three main Iranian ethnicities, i.e., Turk, Kurd, and Fars, is remarkably higher than that in other populations, even that in Southern Iran. High frequencies of CYP2B6^∗2 and ^∗3 in the Iranian population highly affect drug responsiveness. Understanding such variability could help to increase drug efficacy and reduce its toxicity.

Pharmacokinetics of Tamoxifen and Its Major Metabolites and the Effect of the African Ancestry Specific CYP2D6*17 Variant on the Formation of the Active Metabolite, Endoxifen

Tamoxifen (TAM) is widely used in the treatment of hormone receptor-positive breast cancer. TAM is metabolized into the active secondary metabolite endoxifen (ENDO), primarily by CYP2D6. We aimed to investigate the effects of an African-specific CYP2D6 variant allele, CYP2D6*17, on the pharmacokinetics (PK) of TAM and its active metabolites in 42 healthy black Zimbabweans. Subjects were grouped based on CYP2D6 genotypes as CYP2D6*1/*1 or *1/*2 or *2/*2 (CYP2D6*1 or *2), CYP2D6*1/*17 or 2*/*17, and CYP2D6*17/*17. PK parameters for TAM and three metabolites were determined. The pharmacokinetics of ENDO showed statistically significant differences among the three groups. The mean ENDO AUC_0-∞ in CYP2D6*17/*17 subjects was 452.01 (196.94) h·*ng/mL, and the AUC0-∞ in CYP2D6*1/*17 subjects was 889.74 h·ng/mL, which was 5-fold and 2.8-fold lower than in CYP2D6*1 or *2 subjects, respectively. Individuals who were heterozygous or homozygous for CYP2D6*17 alleles showed a 2- and 5-fold decrease in Cmax, respectively, compared to the CYP2D6*1 or *2 genotype. CYP2D6*17 gene carriers have significantly lower ENDO exposure levels than CYP2D6*1 or *2 gene carriers. Pharmacokinetic parameters of TAM and the two primary metabolites, N-desmethyl tamoxifen (NDT) and 4-hydroxy tamoxifen (4OHT), did not show any significant difference in the three genotype groups. The African-specific CYP2D6*17 variant had effects on ENDO exposure levels that could potentially have clinical implications for patients homozygous for this variant.

EA, Udupa KS, Gota V, Mallayasamy S. Factors Influencing Pharmacokinetics of Tamoxifen in Breast Cancer Patients: A Systematic Review of Population Pharmacokinetic Models

BACKGROUND

Tamoxifen is useful in managing breast cancer and it is reported to have significant variability in its pharmacokinetics. This review aimed to summarize reported population pharmacokinetics studies of tamoxifen and to identify the factors affecting the pharmacokinetics of tamoxifen in adult breast cancer patients.

METHOD

A systematic search was undertaken in Scopus, Web of Science, and PubMed for papers published in the English language from inception to 20 August 2022. Studies were included in the review if the population pharmacokinetic modeling was based on non-linear mixed-effects modeling with a parametric approach for tamoxifen in breast cancer patients.

RESULTS

After initial selection, 671 records were taken for screening. A total of five studies were selected from Scopus, Web of Science, PubMed, and by manual searching. The majority of the studies were two-compartment models with first-order absorption and elimination to describe tamoxifen and its metabolites' disposition. The CYP2D6 phenotype and CYP3A4 genotype were the main covariates that affected the metabolism of tamoxifen and its metabolites. Other factors influencing the drug's pharmacokinetics included age, co-medication, BMI, medication adherence, CYP2B6, and CYP2C19 genotype.

CONCLUSION

The disposition of tamoxifen and its metabolites varies primarily due to the CYP2D6 phenotype and CYP3A4 genotype. However, other factors, such as anthropometric characteristics and menopausal status, should also be addressed when accounting for this variability. All these studies should be externally evaluated to assess their applicability in different populations and to use model-informed dosing in the clinical setting.

Delineating gene-environment effects using virtual twins of patients treated with clozapine

Studies that focus on individual covariates, while ignoring their interactions, may not be adequate for model-informed precision dosing (MIPD) in any given patient. Genetic variations that influence protein synthesis should be studied in conjunction with environmental covariates, such as cigarette smoking. The aim of this study was to build virtual twins (VTs) of real patients receiving clozapine with interacting covariates related to genetics and environment and to delineate the impact of interacting covariates on predicted clozapine plasma concentrations. Clozapine-treated patients with schizophrenia (N = 42) with observed clozapine plasma concentrations, demographic, environmental, and genotype data were used to construct VTs in Simcyp. The effect of increased covariate virtualization was assessed by performing simulations under three conditions: "low" (demographic), "medium" (demographic and environmental interaction), and "high" (demographic and environmental/genotype interaction) covariate virtualization. Increasing covariate virtualization with interaction improved the coefficient of variation (R² ) from 0.07 in the low model to 0.391 and 0.368 in the medium and high models, respectively. Whereas R² was similar between the medium and high models, the high covariate virtualization model had improved accuracy, with systematic bias of predicted clozapine plasma concentration improving from -138.48 ng/ml to -74.65 ng/ml. A high level of covariate virtualization (demographic, environmental, and genotype) may be required for MIPD using VTs.

The use of Gammarus pulex as a model organism for ecotoxicological assessment of ibuprofen and propranolol at environmental relevant concentrations

The aim of this study is to assess the toxicity of ibuprofen (IBU) and propranolol (PRO) drugs usingGammarus pulex as a model organism. Firstly, the 96 h LC₅₀ values of IBU and PRO were determined and then three sublethal concentrations of the drugs were exposed to G. pulex. The activities of superoxide dismutase (SOD), catalase (CAT) and acetylcholinesterase (AChE) were evaluated. SOD activity decreased in G. pulex exposed to IBU and PRO compared to control. In all groups exposed to IBU, CAT activity increased at different concentrations at 24 and 96 h. In the groups exposed to different PRO concentrations, CAT activities increased after 24 h compared to the control group (p < 0.05). AChE activities increased in all application groups exposed to IBU for 96 hours (p < 0.05). In conclusion, exposure to IBU and PRO resulted in increased oxidative damage. PRO has been found to cause neurotoxicity.

The impact of CYP2D6*41 on CYP2D6 enzyme activity using phenotyping methods in urine, plasma, and saliva

Aims: The CYP2D6*41 variant is the second or third frequent reduced function allele in Chinese with a frequency of around 3–4%, while it is the major reduced function allele in Indians, Saudi Arabians and Caucasians with frequencies of around 10–20%. The present study was designed to explore the impact of CYP2D6*41 on the metabolic activity of CYP2D6 using phenotyping methods in urine, plasma, and saliva.

Methods: We used dextromethorphan as the probe drug to analyze the phenotypes of 87 subjects with CYP2D6*1/*1 (n = 22), CYP2D6*1/*2 (n = 33), CYP2D6*2/*2 (n = 4), CYP2D6*1/*41 (n = 5), CYP2D6*2/*41 (n = 3), CYP2D6*10/*41 (n = 16), and CYP2D6*5/*41 (n = 4) for CYP2D6. The ratio of parent drug to metabolite in 3 h saliva, 3 h plasma, and in 0–3 h urine was considered the metabolic ratio (MR).

Results: The CYP2D6*41 allele had substantial impact on the metabolic activity of CYP2D6 regardless of the urinary, plasma, or salivary phenotyping method used. In subjects with CYP2D6*1(or *2)/*1(or *2), *1 (or *2)/*41, *10/*41 and *5/*41 (all p < 0.001), the salivary, plasma, or urinary MR value increased. The MRs in saliva, plasma, and urine displayed high correlations.

Conclusion: The activity score system or the consensus activity score system, instead of the traditional phenotype classification, could predict the CYP2D6 enzyme activity more accurately. CYP2D6*41 had similar or more impact on the CYP2D6 enzyme activity as compared with CYP2D6*10. Assigning *41 a score of 0.5 and assigning *10 a score of 0.25 according to the consensus AS system should be reconsidered.

Isolation and Identification of 3,4-Seco-Solanidine-3,4-dioic acid (SSDA) as a Urinary Biomarker of Cytochrome P450 2D6 (CYP2D6) Activity

Cytochrome P450 2D6 (CYP2D6), is responsible for the metabolism and elimination of approximately 25% of clinically used drugs, including antidepressants and antipsychotics, and its activity varies considerably on a population basis primary due to genetic variation. CYP2D6 phenotype can be assessed in vivo following administration of an exogenous probe compound, such as dextromethorphan or debrisoquine, but use of a biomarker that does not require administration of an exogenous compound (i.e., drug) has considerable appeal for assessing CYP2D6 activity in vulnerable populations, such as children. The goal of this study was to isolate, purify and identify an "endogenous" urinary biomarker (M1; m/z 444.3102) of CYP2D6 activity reported previously. Several chromatographic separation techniques (reverse phase HPLC, cation exchange and analytical reverse phase UPLC) were used to isolate and purify 96 μg of M1 from 40 L of urine. Subsequently, 1D and 2D NMR, and functional group modification reactions were used to elucidate its structure. Analysis of mass spectrometry and NMR data revealed M1 to have similar spectroscopic features to the nitrogen-containing steroidal alkaloid, solanidine. 2D NMR characterization by HMBC, COSY, TOCSY, and HSQC-TOCSY proved to be invaluable in the structural elucidation of M1; derivatization of M1 revealed the presence of two carboxylic acid moieties. M1 was determined to be a steroidal alkaloid with a solanidine backbone that had undergone C-C bond scission to yield 3,4-seco-solanidine-3,4-dioic acid (SSDA). SSDA may have value as a dietary biomarker of CYP2D6 activity in populations where potato consumption is common. Significance Statement Endogenous biomarkers of processes involved in drug disposition and response may allow improved individualization of drug treatment, especially in vulnerable populations, such as children. Given that several CYP2D6 substrates are commonly used in pediatrics and the ubiquitous nature of potato consumption in western diets, SSDA has considerable appeal as non-invasive biomarker of CYP2D6 activity to guide treatment with CYP2D6 substrates in children and adults.

Physiologically-Based Pharmacokinetic Models of CYP2D6 Substrate and Inhibitors Nebivolol, Cinacalcet and Mirabegron to Simulate Drug-Drug Interactions

BACKGROUND AND OBJECTIVES

Index substrates and inhibitors to investigate the role of the polymorphic enzyme, cytochrome P450 (CYP) 2D6, in the metabolism of new compounds have been proposed by regulatory agencies. This work describes the development and verification of physiologically-based pharmacokinetic (PBPK) models for the CYP2D6-sensitive substrate, nebivolol and the index CYP2D6 inhibitors, mirabegron and cinacalcet.

METHODS

PBPK models for nebivolol, mirabegron and cinacalcet were developed using in vitro and clinical data. The performance of the PBPK models was verified by comparing the simulated results against reported human systemic exposure and clinical drug-drug interactions (DDIs) studies.

RESULTS

The exposure of nebivolol, cinacalcet and mirabegron predicted by the PBPK models was verified against pharmacokinetic data from 13, 3 and 9 clinical studies, respectively. For nebivolol, the predicted mean maximum plasma concentration (C_max) and area under the plasma concentration-time (AUC) values in CYP2D6 extensive metaboliser subjects were within 0.9- to 1.49-fold of the observed values. In poor metaboliser CYP2D6 subjects, the predicted C_max and AUC values were within 0.41- to 0.81-fold of observed values. For cinacalcet, the predicted C_max and AUC values were within 0.97- to 1.32-fold of the observed data. For mirabegron, the predicted AUC values across all the studies investigated were within 0.71- to 1.88-fold of observed values. The PBPK model-predicted DDIs were in good agreement (within 2-fold) with observed DDIs in all verification studies (n = 8) assessed. The overall precision was 1.26 and 1.21 for C_max and the AUC ratio, respectively.

CONCLUSIONS

The developed PBPK models can be used to assess the DDI potential liability of new chemical entities that are substrates or inhibitors of CYP2D6.

The Effect of CYP2D6 Phenotypes on the Pharmacokinetics of Propafenone: A Systematic Review and Meta-Analysis

Propafenone (PPF) is a class 1C antiarrhythmic agent mainly metabolized by cytochrome (CYP) 2D6, CYP1A2, and CYP3A4. Previous studies have shown that CYP2D6 polymorphism influences the pharmacokinetics (PK) of PPF. However, the small sample sizes of PK studies can lead to less precise estimates of the PK parameters. Thus, this meta-analysis was performed to merge all current PK studies of PPF to determine the effects of the CYP2D6 phenotype more accurately on the PPF PK profile. We searched electronic databases for published studies to investigate the association between the PPF PK and CYP2D6 phenotype. Four PK-related outcomes were included: area under the time–concentration curve (AUC), maximum concentration (Cmax), apparent clearance (CL/F), and half-life (t1/2). A total of five studies were included in this meta-analysis (n = 56). Analyses were performed to compare PK parameters between poor metabolizers (PMs) versus extensive metabolizers (EMs). PPF has a non-linear pharmacokinetics; therefore, analyses were performed according to dose (300 mg and 400 mg). At 300 mg, the AUC mean (95% CI), Cmax, and t1/2 of PPF in PMs were 15.9 (12.5–19.2) µg·h/mL, 1.10 (0.796–1.40) µg/mL, and 12.8 (11.3–14.3) h, respectively; these values were 2.4-, 11.2-, and 4.7-fold higher than those in the EM group, respectively. At 400 mg, a comparison was performed between S- and R-enantiomers. The CL/F was approximately 1.4-fold higher for the R-form compared with the S-form, which was a significant difference. This study demonstrated that CYP2D6 metabolizer status could significantly affect the PPF PK profile. Adjusting the dose of PPF according to CYP2D6 phenotype would help to avoid adverse effects and ensure treatment efficacy.

Pharmacogenomic Profile of Amazonian Amerindians

Given the role of pharmacogenomics in the large variability observed in drug efficacy/safety, an assessment about the pharmacogenomic profile of patients prior to drug prescription or dose adjustment is paramount to improve adherence to treatment and prevent adverse drug reaction events. A population commonly underrepresented in pharmacogenomic studies is the Native American populations, which have a unique genetic profile due to a long process of geographic isolation and other genetic and evolutionary processes. Here, we describe the pharmacogenetic variability of Native American populations regarding 160 pharmacogenes involved in absorption, distribution, metabolism, and excretion processes and biological pathways of different therapies. Data were obtained through complete exome sequencing of individuals from 12 different Amerindian groups of the Brazilian Amazon. The study reports a total of 3311 variants; of this, 167 are exclusive to Amerindian populations, and 1183 are located in coding regions. Among these new variants, we found non-synonymous coding variants in the DPYD and the IFNL4 genes and variants with high allelic frequencies in intronic regions of the MTHFR, TYMS, GSTT1, and CYP2D6 genes. Additionally, 332 variants with either high or moderate (disruptive or non-disruptive impact in protein effectiveness, respectively) significance were found with a minimum of 1% frequency in the Amazonian Amerindian population. The data reported here serve as scientific basis for future design of specific treatment protocols for Amazonian Amerindian populations as well as for populations admixed with them, such as the Northern Brazilian population.

A view of response and resistance to atomoxetine treatment in children with ADHD: effects of CYP2C19 polymorphisms and BDNF levels

OBJECTIVE

Although several genes have previously been studied about the treatment of Attention Deficit Hyperactivity Disorder (ADHD), the number of studies investigating the effects of genes on atomoxetine (ATX) treatment is very limited. In this study, we aimed to investigate the effect of CYP2C19 polymorphisms, which have a role in ATX biotransformation, on the treatment response and also to assess whether there is a relationship between BDNF and treatment response in children and adolescents with ADHD.

METHODS

One hundred children with ADHD and 100 healthy controls (HCs) were included in this study. The treatment response was assessed 2 months after the start of the ATX treatment. DNA samples from peripheral venous blood were replicated using PCR and analyzed using the ILLUMINA next-generation sequencing method. The resulting fastqs were analyzed using Basespace's Variant Interpreter Program. Plasma BDNF levels were evaluated with ELISA kits.

RESULTS

Treatment response was found to be lower in both heterozygous and homozygous carriers of the c.681G > A (CYP2C19*2) polymorphism. When the BDNF level was compared, it was found to be significantly higher in the ADHD group compared to HCs. Also, BDNF has a stronger predictive value for assessing resistance to ATX treatment.

CONCLUSIONS

To our knowledge, this is the first study to assess the effects of CYP2C19 polymorphisms and BDNF levels together on ATX treatment in children. Further studies with an extensive population are needed to better understand the effects of CYP2C19 polymorphisms on treatment and side effects, as well as the effects of BDNF levels.

Therapeutic drug monitoring in children and adolescents with schizophrenia and other psychotic disorders using risperidone

Risperidone is commonly used to treat different psychiatric disorders worldwide. Knowledge on dose–concentration relationships of risperidone treatment in children and adolescents with schizophrenia or other psychotic disorders is, however, scarce and no age-specific therapeutic ranges have been established yet. Multicenter data of a therapeutic drug monitoring service were analyzed to evaluate the relationship between risperidone dose and serum concentration of the active moiety (risperidone (RIS) plus its main metabolite 9-hydroxyrisperidone (9-OH-RIS)) in children and adolescents with psychotic disorders. Patient characteristics, doses, serum concentrations and therapeutic outcomes were assessed by standardized measures. The study also aimed to evaluate whether the therapeutic reference range for adults (20–60 ng/ml) is applicable for minors. In the 64 patients (aged 11–18 years) included, a positive correlation between daily dose and the active moiety (RIS_am) concentration was found (r_s = 0.49, p = 0.001) with variation in dose explaining 24% (r_s² = 0.240) of the variability in serum concentrations. While the RIS_am concentration showed no difference, RIS as well 9-OH-RIS concentrations and the parent to metabolite ratio varied significantly in patients with co-medication of a CYP2D6 inhibitor. Patients with extrapyramidal symptoms (EPS) had on average higher RISam concentrations than patients without (p = 0.05). Considering EPS, the upper threshold of the therapeutic range of RIS_am was determined to be 33 ng/ml. A rough estimation method also indicated a possibly decreased lower limit of the preliminary therapeutic range in minors compared to adults. These preliminary data may contribute to the definition of a therapeutic window in children and adolescents with schizophrenic disorders treated with risperidone. TDM is recommended in this vulnerable population to prevent concentration-related adverse drug reactions.

Exploring the usefulness of plasma level determination and pharmacogenetics for patients treated with clozapine

Aims: The aims of the present study were to assess the variance of plasma clozapine (CLZ) levels and to identify the influence of sociodemographic and pharmacogenetic factors on it and to introduce these tools in a clinical setting. Patients & methods: CLZ concentration was measured and genetic variants of CLZ pharmacokinetic and pharmacodynamic factors were assessed in 23 patients with psychotic disorders. Results: A significant association between mean concentration/dose ratio (C/D) and smoking status, age and weight were found. There was a significant difference in mean plasma CLZ levels and gender. The rs762551 AA genotype in smokers had a significantly lower C/D. Conclusion: In addition to classical factors, monitoring of plasma concentrations together with pharmacogenetics led to greater individualization of treatment.

PacBio long-read amplicon sequencing enables scalable high-resolution population allele typing of the complex CYP2D6 locus

The CYP2D6 enzyme is estimated to metabolize 25% of commonly used pharmaceuticals and is of intense pharmacogenetic interest due to the polymorphic nature of the CYP2D6 gene. Accurate allele typing of CYP2D6 has proved challenging due to frequent copy number variants (CNVs) and paralogous pseudogenes. SNP-arrays, qPCR and short-read sequencing have been employed to interrogate CYP2D6, however these technologies are unable to capture longer range information. Long-read sequencing using the PacBio Single Molecule Real Time (SMRT) sequencing platform has yielded promising results for CYP2D6 allele typing. However, previous studies have been limited in scale and have employed nascent data processing pipelines. We present a robust data processing pipeline "PLASTER" for accurate allele typing of SMRT sequenced amplicons. We demonstrate the pipeline by typing CYP2D6 alleles in a large cohort of 377 Solomon Islanders. This pharmacogenetic method will improve drug safety and efficacy through screening prior to drug administration.

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

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

Individual Pharmacotherapy Management (IPM) - I: a group-matched retrospective controlled clinical study on prevention of complicating delirium in the elderly trauma patients and identification of associated factors

Background

Delirium is one of the most frequent complications in hospitalized elderly patients with additional costs such as prolongation of hospital stays and institutionalization, with risk of reduced functional recovery, long-term cognitive impairment, and increased morbidity and mortality. We analyzed the effect of individual pharmacotherapy management (IPM) in the University Hospital Halle in geriatric trauma patients on complicating delirium and aimed to identify associated factors.

Methods

In a retrospective controlled clinical study of 404 hospitalized trauma patients ≥70 years we compared the IPM intervention group (IG) with a control group (CG) before IPM implementation. Delirium was recorded from the hospital discharge letter. The medication review and data records included baseline data, all medications, diagnoses, electrocardiogram (ECG), laboratory and vital parameters during hospitalization. The IPM internist and the senior trauma physician guaranteed personnel and structural continuity in the implementation of the interdisciplinary patient rounds.

Results

There was a highly matched congruence between CG and IG in terms of age, gender, residency, BMI, most diagnoses, and injury patterns to compare the two groups. The total number of medications per patient was 11.1 ± 4.9 (CG) versus 10.4 ± 3.6 (IG). Our targeted IPM focus on 6 frontline aspects with reduction of antipsychotics, anticholinergic burden, benzodiazepines, serotonergic opioids, elimination of pharmacokinetic and pharmacodynamic drug interactions and overdosage reduced complicating delirium from 5% to almost zero at 0.5%. The association of IPM with a significant 10-fold reduction, OR = 0.09 [95% CI 0.01–0.7], in univariable regression, maintained of clinical relevance in multivariable regression OR = 0.1 [95% CI 0.01–1.1]. Factors most strongly associated with complicating delirium in univariable regression were cognitive dysfunction, nursing home residency, muscle relaxants, antiparkinsonian agents, xanthines, transient disorientation documented in the fall risk scale, antibiotic-requiring infections, antifungals, antipsychotics, and intensive care stay, the two latter maintaining significance in multivariable regression.

Conclusions

IPM is associated with a highly effective prevention of complicating delirium in the elderly trauma patients. For patient safety it should be integrated as an essential preventative contribution. The associated factors help identify patients at risk.

Effects of CYP2D6 Genetic Polymorphism and Drug Interaction on the Metabolism of Dacomitinib

We aim to study the effects of CYP2D6 variants and drug-drug interaction on the metabolism of dacomitinib. CYP2D6 variants were incubated with 25-1000 μM dacomitinib for 40 min at 37 °C, and the reaction was terminated by cooling to -80 °C immediately. For an in vivo experiment, 18 male Sprague-Dawley rats were randomly divided into three groups (n = 6): a single dose of 5 mg/kg dacomitinib (group A), a single dose of 6 mg/kg trazodone (group B), and a combined group (group C). Processed samples were analyzed by ultraperformance liquid chromatography-tandem mass spectrometry (UPLC-MS/MS.) The relative clearance of dacomitinib was reduced for most of the variants. Moreover, the inhibitory potency of classic CYP inhibitors on dacomitinib metabolism was significantly different among the main subtypes of CYP2D6. Interestingly, compared with gefitinib, even the same CYP2D6 variants showed significant differences in metabolic activity, suggesting that the activity of CYP2D6 has strong variability. In addition, the interaction between trazodone and dacomitinib was determined both in vitro and in vivo. When dacomitinib was given in combination with trazodone, the blood exposure to these two drugs increased remarkably. The mechanistic study revealed that the interaction followed the noncompetitive inhibition. We demonstrated that the activity of CYP2D6 variants to metabolize dacomitinib was significantly reduced. In combination with the CYP2D6 inhibitor, the degree of activity inhibition of different variants obviously differed. When trazodone and dacomitinib were used in combination, the body exposure to the two drugs increased significantly. This study provides data for the precise use of dacomitinib in clinical settings.

Machine Learning Prediction of the Three Main Input Parameters of a Simplified Physiologically Based Pharmacokinetic Model Subsequently Used to Generate Time-Dependent Plasma Concentration Data in Humans after Oral Doses of 212 Disparate Chemicals

Physiologically based pharmacokinetic (PBPK) modeling has the potential to play significant roles in estimating internal chemical exposures. The three major PBPK model input parameters (i.e., absorption rate constants, volumes of the systemic circulation, and hepatic intrinsic clearances) were generated in silico for 212 chemicals using machine learning algorithms. These input parameters were calculated based on sets of between 17 and 65 chemical properties that were generated by in silico prediction tools before being processed by machine learning algorithms. The resulting simplified PBPK models were used to estimate plasma concentrations after virtual oral administrations in humans. The estimated absorption rate constants, volumes of the systemic circulation, and hepatic intrinsic clearance values for the 212 test compounds determined traditionally (i.e., based on fitting to measured concentration profiles) and newly estimated had correlation coefficients of 0.65, 0.68, and 0.77 (p < 0.01, n = 212), respectively. When human plasma concentrations were modeled using traditionally determined input parameters and again using in silico estimated input parameters, the two sets of maximum plasma concentrations (r = 0.85, p < 0.01, n = 212) and areas under the curve (r = 0.80, p < 0.01, n = 212) were correlated. Virtual chemical exposure levels in liver and kidney were also estimated using these simplified PBPK models along with human plasma levels. These results indicate that the PBPK model input parameters for humans of a diverse set of compounds can be reliability estimated using chemical descriptors calculated using in silico tools.

Generation of Caco-2 cells stably expressing CYP3A4·POR·UGT1A1 and CYP3A4·POR·UGT1A1*6 using a PITCh system

The small intestine plays a critical role in the absorption and metabolism of orally administered drugs. Therefore, a model capable of evaluating drug absorption and metabolism in the small intestine would be useful for drug discovery. Patients with genotype UGT1A1*6 (exon 1, 211G > A) treated with the antineoplastic drug SN-38 have been reported to exhibit decreased glucuronide conjugation and increased incidence of intestinal toxicity and its severe side effects, including severe diarrhea. To ensure the safety of drugs, we must develop a drug metabolism and toxicity evaluation model which considers UGT1A1*6. In this study, we generated CYP3A4·POR·UGT1A1 KI- and CYP3A4·POR·UGT1A1*6 KI-Caco-2 cells for pharmaceutical research using a PITCh system. The CYP3A4·POR·UGT1A1 KI-Caco-2 cells were shown to express functional CYP3A4 and UGT1A1. The CYP3A4·POR·UGT1A1*6 KI-Caco-2 cells were sensitive to SN-38-induced intestinal toxicity. We thus succeeded in generating CYP3A4·POR·UGT1A1 KI- and CYP3A4·POR·UGT1A1*6 KI-Caco-2 cells, which can be used in pharmaceutical research. We also developed an intestinal epithelial cell model of patients with UGT1A1*6 and showed that it was useful as a tool for drug discovery.

Tools for optimising pharmacotherapy in psychiatry (therapeutic drug monitoring, molecular brain imaging and pharmacogenetic tests): focus on antidepressants

Objectives: More than 40 drugs are available to treat affective disorders. Individual selection of the optimal drug and dose is required to attain the highest possible efficacy and acceptable tolerability for every patient.Methods: This review, which includes more than 500 articles selected by 30 experts, combines relevant knowledge on studies investigating the pharmacokinetics, pharmacodynamics and pharmacogenetics of 33 antidepressant drugs and of 4 drugs approved for augmentation in cases of insufficient response to antidepressant monotherapy. Such studies typically measure drug concentrations in blood (i.e. therapeutic drug monitoring) and genotype relevant genetic polymorphisms of enzymes, transporters or receptors involved in drug metabolism or mechanism of action. Imaging studies, primarily positron emission tomography that relates drug concentrations in blood and radioligand binding, are considered to quantify target structure occupancy by the antidepressant drugs in vivo. Results: Evidence is given that in vivo imaging, therapeutic drug monitoring and genotyping and/or phenotyping of drug metabolising enzymes should be an integral part in the development of any new antidepressant drug.Conclusions: To guide antidepressant drug therapy in everyday practice, there are multiple indications such as uncertain adherence, polypharmacy, nonresponse and/or adverse reactions under therapeutically recommended doses, where therapeutic drug monitoring and cytochrome P450 genotyping and/or phenotyping should be applied as valid tools of precision medicine.

A Drug-Drug Interaction Study Evaluating the Effect of Givosiran, a Small Interfering Ribonucleic Acid, on Cytochrome P450 Activity in the Liver

Givosiran (trade name GIVLAARI) is a small interfering ribonucleic acid that targets hepatic delta-aminolevulinic acid synthase 1 (ALAS1) messenger RNA for degradation through RNA interference (RNAi) that has been approved for the treatment of acute hepatic porphyria (AHP). RNAi therapeutics, such as givosiran, have a low liability for drug-drug interactions (DDIs) because they are not metabolized by cytochrome 450 (CYP) enzymes, and do not directly inhibit or induce CYP enzymes in the liver. The pharmacodynamic effect of givosiran (lowering of hepatic ALAS1, the first and rate limiting enzyme in the heme biosynthesis pathway) presents a unique scenario where givosiran could potentially impact heme-dependent activities in the liver, such as CYP enzyme activity. This study assessed the impact of givosiran on the pharmacokinetics of substrates of 5 major CYP450 enzymes in subjects with acute intermittent porphyria (AIP), the most common type of AHP, by using the validated "Inje cocktail," comprised of caffeine (CYP1A2), losartan (CYP2C9), omeprazole (CYP2C19), dextromethorphan (CYP2D6), and midazolam (CYP3A4). We show that givosiran treatment had a differential inhibitory effect on CYP450 enzymes in the liver, resulting in a moderate reduction in activity of CYP1A2 and CYP2D6, a minor effect on CYP3A4 and CYP2C19, and a similar weak effect on CYP2C9. To date, this is the first study evaluating the DDI for an oligonucleotide therapeutic and highlights an atypical drug interaction due to the pharmacological effect of givosiran. The results of this study suggest that givosiran does not have a large effect on heme-dependent CYP enzyme activity in the liver.

Frequency of Important CYP450 Enzyme Gene Polymorphisms in the Iranian Population in Comparison with Other Major Populations: A Comprehensive Review of the Human Data

Genetic polymorphisms in cytochrome P450 genes can cause alteration in metabolic activity of clinically important medicines. Thus, single nucleotide variants (SNVs) and copy number variations (CNVs) in CYP genes are leading factors of drug pharmacokinetics and toxicity and form pharmacogenetics biomarkers for drug dosing, efficacy, and safety. The distribution of cytochrome P450 alleles differs significantly between populations with important implications for personalized drug therapy and healthcare programs. To provide a meta-analysis of CYP allele polymorphisms with clinical importance, we brought together whole-genome and exome sequencing data from 800 unrelated individuals of Iranian population (100 subjects from 8 major ethnics of Iran) and 63,269 unrelated individuals of five major human populations (EUR, AMR, AFR, EAS and SAS). By integrating these datasets with population-specific linkage information, we evolved the frequencies of 140 CYP haplotypes related to 9 important CYP450 isoenzymes (CYP1A2, CYP2B6, CYP2C8, CYP2C9, CYP2C19, CYP2D6, CYP2E1, CYP3A4 and CYP3A5) giving a large resource for major genetic determinants of drug metabolism. Furthermore, we evaluated the more frequent Iranian alleles and compared the dataset with the Caucasian race. Finally, the similarity of the Iranian population SNVs with other populations was investigated.

Effects of age and gender on the cytochrome P450 2D6 activity in a Korean population

WHAT IS KNOWN AND OBJECTIVE

Age and gender have been reported to play a crucial role in modulating the disposition of pharmacological agents, and to influence the activities of cytochrome P450 (CYP) 2D6, a drug-metabolizing enzyme involved in the disposition of clinically used drugs. In the present study, the effects of age and gender on the CYP2D6 activity were evaluated using dextromethorphan as a probe drug in humans.

METHODS

Healthy young (20 < age < 30 years, n = 60) and old age (age >60 years, n = 60) subjects were enrolled and were given 15 mg dextromethorphan orally. Blood samples were collected before and 3 h after medication. Dextromethorphan and its metabolite dextrorphan were measured using HPLC-fluorescence, and dextromethorphan metabolic ratio (MR, log [dextromethorphan/dextrorphan]) was used to evaluate the CYP2D6 activity.

RESULTS AND DISCUSSION

Mean (±SD) dextromethorphan MR was -2.42 ± 0.46 for the young male group, -2.28 ± 0.56 for the young female group, -2.46 ± 0.55 for the older male group and -2.34 ± 0.65 for the old female group. Based on our findings, the effects of age and gender on CYP2D6 activity were not statistically significant.

WHAT IS NEW AND CONCLUSION

The results of the present study indicate that age and gender play a minor role in the modulation of CYP2D6 activity in the Korean population.

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.

Han Chinese specific cytochrome P450 polymorphisms and their impact on the metabolism of anti-hypertensive drugs with adrenoreceptor blocking properties

Introduction: Cytochrome P450 (CYP) is a monooxygenase superfamily mediating the elimination of anti-hypertensive drugs. Polymorphisms of CYP would lead to differential drug efficacy. Building relationships between genotype and phenotype will benefit individual medical treatment of hypertension.Areas covered: The review systematically summarizes the polymorphisms of four CYPs (CYP2C9, CYP2C19, CYP2D6, and CYP3A4) concentrated distributed in the Han Chinese population. Moreover, the activity of variants on metabolizing anti-hypertensive drugs are reviewed, especially drugs with adrenoceptor blocking properties, as well as their clinical relevancies.Expert opinion: The polymorphisms of CYP can cause stratification in drug exposure of antihypertensive drugs. Although the clinical relevance has been built partially, the translational medicine still lacks reliable data support. Furthermore, the studies have demonstrated that even the same CYP variant will exhibit different catalytic capability for different drugs, which is another obstacle to hinder its application. With the deepening of multiomics research and structural biology, nucleotide polymorphisms can be combined with transcriptome, proteome, metabolome and molecular structure analyses to study the susceptibility to hypertension and drug efficacy. A complete data chain would be further estabolished by combining studies of pharmacokinetics-pharmacodynamics, which can effectively promote the precise application of anti-hypertensive drugs.

Population Pharmacokinetics of Primaquine in the Korean Population

While primaquine has long been used for malaria treatment, treatment failure is common. This study aims to develop a population pharmacokinetic model of primaquine and its metabolite, carboxyprimaquine, and examine factors influencing pharmacokinetic variability. The data was obtained from a clinical study in 24 Korean subjects randomly assigned to normal and obese groups. The participants received primaquine 15 mg daily for 4 days and blood samples were collected at day 4. Pharmacokinetic modeling was performed with NONMEM and using simulations; the influences of doses and covariates on drug exposure were examined. A minimal physiology-based pharmacokinetic model connected with a liver compartment comprehensively described the data, with CYP450 mediated clearance being positively correlated with the body weight and CYP2D6 activity score (p < 0.05). In the simulation, while the weight-normalized area under drug concentration for primaquine in the obese group decreased by 29% at the current recommended dose of 15 mg, it became similar to the normal weight group at a weight-normalized dose of 3.5 mg/kg. This study has demonstrated that the body weight and CYP2D6 activity score significantly influence the pharmacokinetics of primaquine. The developed model is expected to be used as a basis for optimal malaria treatment in Korean patients.

Cytochrome P450 2D6 genotype-phenotype characterization through population pharmacokinetic modeling of tedatioxetine

The cytochrome P450 (CYP) 2D6 enzyme exhibits large interindividual differences in metabolic activity. Patients are commonly assigned a CYP2D6 phenotype based on their CYP2D6 genotype, but there is a lack of consensus on how to translate genotypes into phenotypes, causing inconsistency in genotype‐based dose recommendations. The aim of this study was to quantify and compare the impact of different CYP2D6 genotypes and alleles on CYP2D6 metabolism using a large clinical data set. A population pharmacokinetic (popPK) model of tedatioxetine and its CYP2D6‐dependent metabolite was developed based on pharmacokinetic data from 578 subjects. The CYP2D6‐mediated metabolism was quantified for each subject based on estimates from the final popPK model, and CYP2D6 activity scores were calculated for each allele using multiple linear regression. The activity scores estimated for the decreased function alleles were 0.46 (CYP2D6*9), 0.34 (CYP2D6*10), 0.01 (CYP2D6*17), 0.65 (CYP2D6*29), and 0.21 (CYP2D6*41). The CYP2D6*17 and CYP2D6*41 alleles were thus associated with the lowest CYP2D6 activity, although only the difference to the CYP2D6*9 allele was shown to be statistically significant (p = 0.02 and p = 0.05, respectively). The study provides new in vivo evidence of the enzyme function of different CYP2D6 genotypes and alleles. Our findings suggest that the activity score assigned to CYP2D6*41 should be revisited, whereas CYP2D6*17 appears to exhibit substrate‐specific behavior. Further studies are needed to confirm the findings and to improve the understanding of CYP2D6 genotype–phenotype relationships across substrates.

Genetic variability of CYP2D6, CYP3A4 and CYP3A5 among the Egyptian population

Aim: This study investigated major allelic variants of CYP2D6, CYP3A4 and CYP3A5 in Egyptians, an Arabic population for which there is little information regarding these important pharmacogenes. Patients & methods: CYP2D6*2, *4, *5, *10, *41 and gene copy number variation, as well as CYP3A4*22 and CYP3A5*3 were determined with commercially available TaqMan assays in 145 healthy study participants. Results: The CYP2D6 alleles identified suggest that the prevalence of poor metabolizers is low as none were found among the 145 subjects investigated. The frequency for CYP3A5 nonexpressers was 74.5% and the CYP3A4*22 allele frequency was low at 2.0%. Conclusion: These preliminary findings indicate that pharmacogene variation in Egyptians is different from those of other Middle Eastern/Arabic populations and warrants further investigation.

Pharmacokinetics of duloxetine self-administered in overdose with quetiapine and other antipsychotic drugs in a Japanese patient admitted to hospital

Background

Combinations of antidepressant duloxetine (at doses of 40–60 mg/day) and other antipsychotics are frequently used in clinical treatment; however, several fatal and nonfatal cases of duloxetine overdose have been documented. We experienced a patient who had taken an overdose of duloxetine (780 mg) in combination with other drugs in a suicide attempt.

Case presentation

The patient was a 37-year-old man (body weight, 64 kg) with a history of gender identity disorder and depression. He intentionally took an overdose of duloxetine in combination with three other antipsychotic drugs (18 mg flunitrazepam, 850 mg quetiapine, and 1100 mg trazodone) and was emergently admitted to Kyoto Medical Center. The patient’s plasma concentration of duloxetine during ambulance transport was 57 ng/ml, and the level was still as high as 126 ng/mL at 32 h after administration. Duloxetine disappeared most slowly from plasma, in contrast to quetiapine, which was the fastest to clear among the four medicines determined in this patient. The observed concentrations of duloxetine in this overdose patient were generally within the 95% confidence intervals of the plasma concentration curves predicted using a physiologically based pharmacokinetic (PBPK) model.

Conclusion

Even if more than 1 h (the generally recommended period) has passed after administration of duloxetine in such overdose cases, gastric lavage and/or administration of activated charcoal may be effective in clinical practice up to 6 h because of the typically slow elimination behavior illustrated by the PBPK model. Pharmacokinetic profiles visualized using PBPK modeling can inform treatment decisions in cases of drug overdose for medicines such as duloxetine in emergency clinical practice.

Application of Model Informed Precision Dosing to Address the Impact of Pregnancy Stage and CYP2D6 Phenotype on Foetal Morphine Exposure

Guidance regarding the effect of codeine and its metabolites on foetal development is limited by small studies and inconsistent findings. The primary objective was to use physiologically based pharmacokinetic modelling to investigate the impact of gestational stage and maternal CYP2D6 phenotype on foetal morphine exposure following codeine administration. Full body physiologically based pharmacokinetic models were developed and verified for codeine and morphine using Simcyp (version 19.1). The impact of gestational age and maternal CYP2D6 phenotype on foetal and maternal morphine and codeine exposure following oral codeine administration was modelled in a cohort of 250 pregnant females and foetuses at gestational weeks 0 (mothers only), 6, 12, 24 and 36. Consistent with the known effect on codeine metabolism, a clinically meaningful (> 1.65-fold) increase in foetal morphine AUC was observed in the CYP2D6 UM phenotype cohort compared to the CYP2D6 EM and PM phenotype cohorts. The mean (95% CI) foetal morphine AUC in the CYP2D6 UM cohort of 0.988 (0.902 to 1.073) ng/mL.h was 1.8-fold higher than the CYP2D6 EM cohort of 0.546 (0.492 to 0.600) ng/mL.h (p < 0.001). Despite a 2.8-fold increase in maternal CYP2D6 protein abundance between gestational weeks 6 and 36, the mean foetal morphine AUC in the CYP2D6 EM and UM phenotype cohorts reduced by 1.55- and 1.75-fold, respectively, over this period. Maternal CYP2D6 phenotype is a significant determinant of foetal morphine AUC. Simulations suggest that the greatest risk with respect to foetal morphine exposure is during the first trimester of pregnancy, particularly in CYP2D6 UM phenotype mothers.

Could Cytochrome P450 2D6, 3A4 and 3A5 Polymorphisms Explain the Variability in Clinical Response to Clomiphene Citrate of Anovulatory PCOS Women?

INTRODUCTION

Cytochrome P450 2D6, 3A4 and 3A5 are involved in the metabolism of many drugs. These enzymes have a genetic polymorphism responsible for different metabolic phenotypes. They play a role in the metabolism of clomiphene citrate (CC), which is used to induce ovulation. Response to CC treatment is variable, and no predictive factors have thus far been identified.

OBJECTIVE

To study a possible link between the cytochrome P450 2D6, 3A4 and 3A5 polymorphisms and clinical response to CC.

STUDY DESIGN

Seventy-seven women with anovulatory Polycystic Ovarian Syndrome (PCOS) treated with CC were included which determined their cytochrome P450 2D6, 3A4 and 3A5 genotypes and used the results to predict ovarian response to this drug. Predicted responses based on the cytochrome genotypes were compared with the observed clinical responses using the calculation of a weighted Kappa coefficient.

MAIN OUTCOME MEASURES

Number of dominant follicles assessed by ultrasound at the end of the follicular phase and confirmation of ovulation by blood progesterone assay in the luteal phase.

RESULTS

Concordance between the predicted and observed responses for the combination of the three cytochromes was 36.71%, with a negative Kappa coefficient (K = -0.0240), which corresponds to a major disagreement. Similarly, for predictions based on the cytochrome P450 2D6 genotype alone, only 39.24% of predictions were verified (coefficient K = -0.0609).

CONCLUSION

The genetic polymorphism of cytochromes P450 2D6, 3A4 and 3A5 does not appear to influence clinical response to CC used to induce ovulation in anovulatory PCOS women.

Sources of Interindividual Variability

The efficacy, safety, and tolerability of drugs are dependent on numerous factors that influence their disposition. A dose that is efficacious and safe for one individual may result in sub-therapeutic or toxic blood concentrations in others. A significant source of this variability in drug response is drug metabolism, where differences in presystemic and systemic biotransformation efficiency result in variable degrees of systemic exposure (e.g., AUC, C_max, and/or C_min) following administration of a fixed dose.Interindividual differences in drug biotransformation have been studied extensively. It is recognized that both intrinsic factors (e.g., genetics, age, sex, and disease states) and extrinsic factors (e.g., diet , chemical exposures from the environment, and the microbiome) play a significant role. For drug-metabolizing enzymes, genetic variation can result in the complete absence or enhanced expression of a functional enzyme. In addition, upregulation and downregulation of gene expression, in response to an altered cellular environment, can achieve the same range of metabolic function (phenotype), but often in a less predictable and time-dependent manner. Understanding the mechanistic basis for variability in drug disposition and response is essential if we are to move beyond the era of empirical, trial-and-error dose selection and into an age of personalized medicine that will improve outcomes in maintaining health and treating disease.

Impact of CYP2D6 on serum concentrations of flupentixol, haloperidol, perphenazine and zuclopenthixol

AIMS

To investigate the impact of cytochrome P450 2D6 (CYP2D6) on dose-adjusted serum concentrations of flupentixol, haloperidol, perphenazine and zuclopenthixol in a therapeutic drug monitoring (TDM) cohort of psychiatric patients. We also studied the functional impact of CYP2D6*41 on dose-adjusted serum concentrations in the perphenazine-treated patients.

METHODS

Serum concentrations of flupentixol, haloperidol, perphenazine and zuclopenthixol from CYP-genotyped patients were extracted retrospectively from a routine TDM database in the period March 2005 to May 2019. Samples were divided into three CYP2D6 phenotype subgroups according to genotype; normal metabolizers (NMs), intermediate metabolizers (IMs) and poor metabolizers (PMs). The effect of CYP2D6 phenotype on dose-adjusted serum concentrations of the four antipsychotics was evaluated by multivariable mixed model analyses.

RESULTS

Mean dose-adjusted serum concentrations of perphenazine (564 samples) were 3.9-fold and 1.6-fold higher in CYP2D6 PMs and IMs, respectively, compared with NMs (P < .001 and P < .01). For zuclopenthixol (658 samples), mean dose-adjusted serum concentrations were about 1.5-fold and 1.3-fold higher in CYP2D6 PMs and IMs, respectively, compared with NMs (P < .01 and P < .001). CYP2D6 was of minor or no importance to haloperidol (320 samples) and flupentixol (115 samples). In our data material, the genotype CYP2D6 *1/*41 appears to have a similar impact on dose-adjusted serum concentrations of perphenazine as *1/null (null = variant allele encoding no enzyme function).

CONCLUSIONS

This study shows that CYP2D6 is important for the metabolism of perphenazine and zuclopenthixol, but not for haloperidol and flupentixol. The CYP2D6*41 allele appears to have a reduced function close to nonfunctional variant alleles.

CYP2D6 Predicts Plasma Donepezil Concentrations in a Cohort of Thai Patients with Mild to Moderate Dementia

Purpose

Donepezil, a drug frequently used to treat dementia, is mainly metabolized by cytochrome P450 2D6 (CYP2D6). This study investigated the relationships between CYP2D6 genotype and activity scores as well as predicted phenotype of plasma donepezil concentrations in 86 Thai dementia participants.

Materials and Methods

CYP2D6 was genotyped using bead-chip technology (Luminex xTAG^® v.3). Steady-state trough plasma donepezil concentrations were measured using high-performance liquid chromatography.

Results

Sixteen genotypes were found but the most frequent genotypes detected among our participants were CYP2D6*10/*10 (27.9%) and *1/*10 (26.7%). One-third of the participants had an activity score of 1.25 which predicted that they were normal metabolizers. The overall median (interquartile range) of plasma donepezil concentration was 51.20 (32.59–87.24) ng/mL. Normal metabolizers (NMs) had lower plasma donepezil concentrations compared to intermediate metabolizers (IMs) (41.15 (28.44–67.65) ng/mL vs 61.95 (35.25–97.00) ng/mL). Multivariate analysis showed that CYP2D6 activity score (r² = 0.50) and the predicted phenotype (independent of dose) could predict the plasma donepezil concentration (r² = 0.49).

Conclusion

Plasma donepezil concentration in NMs was lower compared to IMs. Additional studies with larger sample size and use of next-generation sequencing as well as its outcomes are warranted to confirm the benefit of using pharmacogenetic-guided treatment for donepezil.

Transfer learning enables prediction of CYP2D6 haplotype function

Cytochrome P450 2D6 (CYP2D6) is a highly polymorphic gene whose protein product metabolizes more than 20% of clinically used drugs. Genetic variations in CYP2D6 are responsible for interindividual heterogeneity in drug response that can lead to drug toxicity and ineffective treatment, making CYP2D6 one of the most important pharmacogenes. Prediction of CYP2D6 phenotype relies on curation of literature-derived functional studies to assign a functional status to CYP2D6 haplotypes. As the number of large-scale sequencing efforts grows, new haplotypes continue to be discovered, and assignment of function is challenging to maintain. To address this challenge, we have trained a convolutional neural network to predict functional status of CYP2D6 haplotypes, called Hubble.2D6. Hubble.2D6 predicts haplotype function from sequence data and was trained using two pre-training steps with a combination of real and simulated data. We find that Hubble.2D6 predicts CYP2D6 haplotype functional status with 88% accuracy in a held-out test set and explains 47.5% of the variance in in vitro functional data among star alleles with unknown function. Hubble.2D6 may be a useful tool for assigning function to haplotypes with uncurated function, and used for screening individuals who are at risk of being poor metabolizers.