Pharmacogenomics (PGx) investigates the influence of genetics on drug responses, enabling tailored treatments for personalized healthcare. This study assessed the accuracy of genotyping six genes using whole genome sequencing with four different computational tools and various sequencing depths. The effects of using different reference genomes (GRCh38 and GRCh37) and sequence aligners (BWA-MEM and Bowtie2) were also explored. The results showed generally minor variations in tool performance across most genes; however, more notable discrepancies were observed in the analysis of the complex CYP2D6 gene. Cyrius, a CYP2D6-specific tool, demonstrated the most robust performance, achieving the highest concordance rates for CYP2D6 in all instances, comparable to the consensus approach in most cases. There were rather small differences between the samples with 20× coverage depth and those with higher depth, but the decreased performance was more evident at lower depths, particularly at 5×. Additionally, variations in CYP2D6 results were observed when samples were aligned to different reference genomes using the same method, or to the same genome using different aligners, which led to reporting incorrect rare star alleles in several cases. These findings inform the selection of optimal PGx tools and methodologies as well as suggest that employing a consensus approach with two or more tools might be preferable for certain genes and tool combinations, especially at lower sequencing depths, to ensure accurate results. Additionally, we show how the upstream alignment can affect the performance of tools, an important factor to take into account.

CYP2D6 variants contain various single nucleotide polymorphisms as well as differing levels of metabolic activity. Among these, one of the less active variants CYP2D6*10 (100C > T) is the most prevalent mutation in East Asians, including Japanese. This mutation leads to an amino acid substitution from proline to serine, which reduces the stability of CYP2D6 and consequently decreases its metabolic activity. In this study, we used a genome editing technology called the Precise Integration into Target Chromosome (PITCh) system to stably express six drug-metabolizing enzymes (CYP3A4, POR, uridine diphosphate glucuronosyltransferase 1A1 (UGT1A1), CYP1A2, CYP2C19, CYP2C9, and CYP2D6*10) in HepG2 (CYP2D6*10 KI-HepG2) cells to examine the effect of CYP2D6*10 on drug metabolism prediction. The protein expression levels of CYP2D6 in CYP2D6*10 KI-HepG2 cells were reduced relative to those in the CYP3A4-POR-UGT1A1-CYP1A2-CYP2C19-CYP2C9-CYP2D6 knock-in-HepG2 (CYPs-UGT1A1 KI-HepG2) cells. Consistent with the CYP2D6 protein expression results, CYP2D6 metabolic activity in CYP2D6*10 KI-HepG2 cells was reduced relative to CYPs-UGT1A1 KI-HepG2 cells. We successfully generated CYP2D6*10 KI-HepG2 cells with highly expressed, functional CYP2D6*10, as well as CYP1A2, 2C9, 2C19 and 3A4. CYP2D6*10 KI-HepG2 cells could be an invaluable model for hepatic metabolism and hepatotoxicity studies in East Asians, including Japanese.

Cytochrome P450 2D6 (CYP2D6) plays a crucial role in metabolizing approximately 20% of medications prescribed clinically. This enzyme is encoded by the CYP2D6 gene, known for its extensive polymorphism with over 170 catalogued haplotypes or star alleles, which can have a profound impact on drug efficacy and safety. Despite its importance, a gap exists in the global genomic databases, which are predominantly representative of European ancestries, thereby limiting comprehensive knowledge of CYP2D6 variation in ethnically diverse populations. In an effort to bridge this knowledge gap, we focused on elucidating the CYP2D6 variation landscape within a multi-ethnic Asian cohort, encompassing individuals of Chinese, Malay, and Indian descent. Our study comprised data analysis of 1850 whole genomes from the SG10K_Health dataset using an in-house consensus algorithm, which integrates the capabilities of Cyrius, Aldy, and StellarPGx. This analysis unveiled distinct population-specific star-allele distribution trends, highlighting the unique genetic makeup of the Singaporean population. Significantly, 46% of our cohort harbored actionable CYP2D6 variants-those with direct implications for drug dosing and treatment strategies. Furthermore, we identified 14 potential novel CYP2D6 star-alleles, of which 7 were observed in multiple individuals, suggesting their broader relevance. Overall, our study contributes novel data on CYP2D6 genetic variations specific to the Southeast Asian context. The findings are instrumental for the advancement of pharmacogenomics and personalized medicine, not only in Southeast Asia but also in other regions with comparable genetic diversity.

Cytochrome P450 enzymes including CYP2C19 and CYP2D6 are important for antidepressant metabolism and polymorphisms of these genes have been determined to predict metabolite levels. Nonetheless, more evidence is needed to understand the impact of genetic variations on antidepressant response. In this study, individual clinical and genetic data from 13 studies of European and East Asian ancestry populations were collected. The antidepressant response was clinically assessed as remission and percentage improvement. Imputed genotype was used to translate genetic polymorphisms to metabolic phenotypes (poor, intermediate, normal, and rapid+ultrarapid) of CYP2C19 and CYP2D6. CYP2D6 structural variants cannot be imputed from genotype data, limiting the determination of metabolic phenotypes, and precluding testing for association with response. The association of CYP2C19 metabolic phenotypes with treatment response was examined using normal metabolizers as the reference. Among 5843 depression patients, a higher remission rate was found in CYP2C19 poor metabolizers compared to normal metabolizers at nominal significance but did not survive after multiple testing correction (OR = 1.46, 95% CI [1.03, 2.06], p = 0.033, heterogeneity I2 = 0%, subgroup difference p = 0.72). No metabolic phenotype was associated with percentage improvement from baseline. After stratifying by antidepressants primarily metabolized by CYP2C19, no association was found between metabolic phenotypes and antidepressant response. Metabolic phenotypes showed differences in frequency, but not effect, between European- and East Asian-ancestry studies. In conclusion, metabolic phenotypes imputed from genetic variants using genotype were not associated with antidepressant response. CYP2C19 poor metabolizers could potentially contribute to antidepressant efficacy with more evidence needed. Sequencing and targeted pharmacogenetic testing, alongside information on side effects, antidepressant dosage, depression measures, and diverse ancestry studies, would more fully capture the influence of metabolic phenotypes.

Human cytochrome P450 (CYP) proteins metabolize 75% of small-molecule pharmaceuticals, which makes structure-based modeling of CYP metabolism and inhibition, bolstered by the current availability of X-ray crystal structures of CYP globular catalytic domains, an attractive prospect. Accounting for this broad metabolic capacity is a combination of the existence of multiple different CYP proteins and the capacity of a single CYP protein to metabolize multiple different small molecules. It is thought that structural plasticity and flexibility contribute to this latter property; therefore, incorporating diverse conformational states of a particular CYP is likely an important consideration in structure-based CYP metabolism and inhibition modeling. While all-atom explicit-solvent molecular dynamics simulations can be used to generate conformational ensembles under biologically relevant conditions, existing CYP crystal structures are of the globular domain only, whereas human CYPs contain N-terminal transmembrane and linker peptides that anchor the globular catalytic domain to the endoplasmic reticulum. To determine whether this can cause significant differences in the sampled binding site conformations, microsecond scale all-atom explicit-solvent molecular dynamics simulations of the CYP2D6 globular domain in an aqueous environment were compared with those of the full-length protein anchored in a POPC lipid bilayer. While bilayer-anchoring damped some structural fluctuations in the globular domain relative to the aqueous simulations, none of the affected residues included binding site pocket residues. Furthermore, clustering of molecular dynamics snapshots based on either pairwise binding site pocket RMSD or volume differences demonstrated a lack of separation of snapshots from the two simulation conditions into different clusters. These results suggest the substantially simpler and computationally cheaper aqueous simulation approach can be used to generate a relevant conformational ensemble of the CYP2D6 binding site for structure-based metabolism and inhibition modeling.

Cytochrome P450 2D6 (CYP2D6) is one of the most important enzymes involved in drug metabolism. Genetic polymorphism can influence drug metabolism by CYP2D6 such that a therapy is seriously affected by under- or overdosing of drugs. However, a general explanation at the atomistic level for poor activity is missing so far. Here we show for the 20 most common single nucleotide polymorphisms (SNPs) of CYP2D6 that poor metabolism is driven by four mechanisms. We found in extensive all-atom molecular dynamics simulations that the rigidity of the I-helix (central helix), distance between central phenylalanines (stabilizing bound substrate), availability of basic residues on the surface of CYP2D6 (binding of cytochrome P450 reductase), and position of arginine 132 (electron transfer to heme) are essential for an extensive function of the enzyme. These results were applied to SNPs with unknown effects, and potential SNPs that may lead to poor drug metabolism were identified. The revealed molecular mechanisms might be important for other drug-metabolizing cytochrome P450 enzymes.

Schizophrenia spectrum disorders (SSD) are a group of diseases characterized by one or more abnormal features in perception, thought processing and behavior. Patients suffering from SSD are at risk of developing life-threatening complications. Pharmacogenetic studies have shown promising results on personalized treatment of psychosis. In the current study, 103 patients diagnosed with SSD treated with risperidone as antipsychotic monotherapy were enrolled. Socio-demographics and clinical data were recorded, and laboratory tests and genotyping standard procedure for cytochrome P450 (CYP) 2D6*4 were performed. Patients were evaluated by the Positive and Negative Syndrome Scale (PANSS) on admission and at discharge. Based on the reduction in the PANSS total score, subjects were divided into non-responders, partial responders and full responders. Only 11 subjects had a full response to risperidone (10.67%), 53 subjects (51.45%) had a partial response, and 39 participants (37.86%) were non-responders. Patients at first episode psychosis showed significantly higher levels of blood glucose and prolactin levels, while chronic patients showed significantly higher LDL levels. Adverse drug reactions (ADR) such as tremor and stiffness significantly correlated with genetic phenotypes (p = 0.0145). While CYP2D6 showed no impact on treatment response, ADR were significantly more frequent among poor and intermediate metabolizers.

The cytochrome P450 (CYP) family of enzymes plays a central role in the metabolism of many drugs. CYP genes are highly polymorphic, which is known to affect protein levels, but for some low frequent CYP genotypes the correlation between genotype and CYP protein expression is less established. In this study, we determined the CYP2C9, CYP2C19, CYP2D6, and CYP3A5 genotypes of 250 Danish individuals included in a postmortem study. For 116 of the individuals, the hepatic CYP protein levels were investigated by a proteomics approach. Overall, we found the postmortem genetic and proteomic data to be in agreement with those of other studies performed on fresh hepatic tissue, showing the usability of postmortem hepatic tissue for this type of investigation. For less investigated genotypes, we could corroborate previously found results: 1) statistically significantly lower levels of hepatic CYP2C9 protein in individuals carrying the CYP2C9*3 variant compared with individuals with two wild type (wt) alleles; 2) comparable levels of CYP2C19 in CYP2C19*2/*17 and CYP2C19*1/*2 individuals; 3) reduced CYP2D6 protein levels in heterozygous individuals with the CYP2D6*3, CYP2D6*4, and CYP2D6*5 gene deletion variants; and 4) significantly lower levels of CYP3A5 protein in CYP3A5*3 homozygous individuals compared with individuals who were heterozygous for the CYP3A5*3 allele or homozygous individuals for the wt alleles. In conclusion, the use of postmortem tissue significantly increases the access to human specimens for research purposes, and postmortem proteomics can be used to investigate the link between CYP genotypes and hepatic protein expression. SIGNIFICANCE STATEMENT: In tissue samples from a large postmortem cohort (n = 250) we determined the CYP2C9, CYP2C19, CYP2D6, and CYP3A5 genotypes. Hepatic CYP protein levels were investigated in 116 individuals using a proteomics approach. For common genotypes, we found results similar to previous knowledge, pointing toward the usability of postmortem tissue. For the less investigated genotypes, we were able to corroborate genotype/protein expression correlations. It is a novel approach to use a large postmortem cohort to investigate genetic/protein expression correlations.

BACKGROUND: Several opioids have pharmacogenetic and drug-drug interactions which may compromise their analgesic effectiveness, but are not routinely implemented into supportive pain management. We hypothesized that CYP2D6 phenotypes and concomitant use of CYP2D6 substrates or inhibitors would correlate with opioid analgesic outcomes.
METHODS: An observational cross-sectional study was conducted with 263 adult chronic non cancer pain (CNCP) patients from a real-world pain unit under long-term CYP2D6-related opioid treatment (tramadol, hydromorphone, tapentadol or oxycodone). Metabolizer phenotype (ultrarapid [UM], normal [NM], intermediate [IM] or poor [PM]) was determined by the CYP2D6 genotype. The socio-demographic (sex, age, employment status), clinical (pain intensity and relief, neuropathic component, quality of life, disability, anxiety and depression), pharmacological (opioid doses and concomitant pharmacotherapy) and safety (adverse events) variables were recorded.
RESULTS: The whole population (66 % female, 65 (14) years old, 70 % retired and 63 % attended for low back pain) were classified as PM (5 %), IM (32 %), NM (56 %) and UM (6 %). Multiple linear and logistic regressions showed higher pain intensity and neuropathic component at younger ages when using any CYP2D6 substrate (p = 0.022) or inhibitor (p = 0.030) drug, respectively, with poorer pain relief when CYP2D6 inhibitors (p=0.030) were present.
CONCLUSIONS: The concomitant use of CYP2D6 substrates or inhibitors during opioid therapy for CNCP may result in lack of analgesic effectiveness. This aspect could be relevant for pharmacological decision making during CNCP management.

OBJECTIVE: Primaquine is essential for the radical cure of Plasmodium vivax malaria and must be metabolized into its bioactive metabolites. Accordingly, polymorphisms in primaquine-metabolizing enzymes can impact the treatment efficacy. This pioneering study explores the influence of monoamine oxidase-A (MAO-A) on primaquine metabolism and its impact on malaria relapses.
METHODS: Samples from 205 patients with P. vivax malaria were retrospectively analysed by genotyping polymorphisms in MAO-A and cytochrome P450 2D6 (CYP2D6) genes. We measured the primaquine and carboxyprimaquine blood levels in 100 subjects for whom blood samples were available on the third day of treatment. We also examined the relationship between the enzyme variants and P. vivax malaria relapses in a group of subjects with well-documented relapses.
RESULTS: The median carboxyprimaquine level was significantly reduced in individuals carrying low-expression MAO-A alleles plus impaired CYP2D6. In addition, this group experienced significantly more P. vivax relapses. The low-expression MAO-A status was not associated with malaria relapses when CYP2D6 had normal activity. This suggests that the putative carboxyprimaquine contribution is irrelevant when the CYP2D6 pathway is fully active.
CONCLUSIONS: We found evidence that the low-expression MAO-A variants can potentiate the negative impact of impaired CYP2D6 activity, resulting in lower levels of carboxyprimaquine metabolite and multiple relapses. The findings support the hypothesis that carboxyprimaquine may be further metabolized through CYP-mediated pathways generating bioactive metabolites that act against the parasite.