CYP2D6 Allelic Variants *34, *17-2, *17-3, and *53 and a Thr309Ala Mutant Display Altered Kinetics and NADPH Coupling in Metabolism of Bufuralol and Dextromethorphan and Altered Susceptibility to Inactivation by SCH 66712

The Identification of Novel CYP2D6 Variants in US Hmong: Results From Genome Sequencing and Clinical Genotyping

Objective: Hmong individuals represent a unique East Asian subpopulation in whom limited information concerning pharmacogenetic variation exists. The objectives of this study were to comprehensively characterize the highly polymorphic CYP2D6 gene in Hmong, estimate allele and phenotype frequencies and to compare results between two testing platforms. Methods: DNA from 48 self-identified Hmong participants were sequenced using a targeted next-generation sequencing (NGS) panel. Star allele calls were made using Astrolabe, manual inspection of NGS variant calls and confirmatory Sanger sequencing. Structural variation was determined by long-range (XL)-PCR and digital droplet PCR (ddPCR). The consensus diplotypes were subsequently translated into phenotype utilizing the activity score system. Clinical grade pharmacogenetic testing was obtained for 12 of the 48 samples enabling an assessment of concordance between the consensus calls and those determined by clinical testing platforms. Results: A total of 13 CYP2D6 alleles were identified. The most common alleles were CYP2D6*10 and its structural arrangements (37.5%, 36/96) and the *5 gene deletion (13.5%, 13/96). Three novel suballeles (*10.007, *36.004, and *75.002) were also identified. Phenotype frequencies were as follows: ultrarapid metabolizers (4.2%, 2/48), normal metabolizers (41.7%, 20/48) and intermediate metabolizers (52.1%, 25/48); none of the 48 participants were predicted to be poor metabolizers. Concordance of diplotype and phenotype calls between the consensus and clinical testing were 66.7 and 50%, respectively. Conclusion: Our study to explore CYP2D6 genotypes in the Hmong population suggests that this subpopulation is unique regarding CYP2D6 allelic variants; also, a higher portion of Hmong participants (50%) are predicted to have an intermediate metabolizer phenotype for CYP2D6 compared to other East Asians which range between 27 and 44%. Results from different testing methods varied considerably. These preliminary findings underscore the importance of thoroughly interrogating unique subpopulations to accurately predict a patient's CYP2D6 metabolizer status.

Tryptophan-75 Is a Low-Energy Channel-Gating Residue that Facilitates Substrate Egress/Access in Cytochrome P450 2D6

CYP2D6 is a major drug metabolizing enzyme with a buried active site. Channels leading to the active site from various enzyme surfaces are believed to facilitate ligand egress and access to the active site. The present study used molecular dynamics (MD) and in vitro studies with CYP2D6*1 and a Trp75-to-Ala mutant to examine channel gating in CYP2D6 by Trp75. MD simulations measured energy landscapes of Trp75 conformations and simulated substrate passage within channel 2b using bufuralol as a model substrate. Trp75 alternated between multiple stable states that supported substrate transport along channel 2b with low-energy barriers between states (∼ -1 kcal/mol). Trp75 conformations were stabilized primarily by hydrogen bonding between Trp75 and Glu222, Asn226, Ala225, or Gln72. Energy barriers were low between Trp75 conformations, allowing Trp75 to easily move between various conformations over time and to function in both binding to and moving substrates in the 2b channel of CYP2D6. Michaelis-Menten kinetic studies completed with purified enzyme in a reconstituted system showed overall reduced enzyme efficiency for metabolism of bufuralol and dextromethorphan by the Trp75Ala mutant compared with CYP2D6*1. In stopped-flow measurements, k off for dextromethorphan was decreased in the absence of Trp75. Our results support a role for Trp75 in substrate shuttling to the active site of CYP2D6. SIGNIFICANCE STATEMENT: Using combined molecular dynamics and in vitro assays, this study shows for the first time a role for Trp75 as a channel entrance gating residue in the mechanism of substrate binding/unbinding in CYP2D6. Energy landscapes derived from molecular dynamics were used to quantitate the strength of gating, and kinetics assays showed the impact on enzyme efficiency and k off of a Trp75Ala mutation.

Spontaneous Ligand Access Events to Membrane-Bound Cytochrome P450 2D6 Sampled at Atomic Resolution

The membrane-anchored enzyme Cytochrome P450 2D6 (CYP2D6) is involved in the metabolism of around 25% of marketed drugs and its metabolic performance shows a high interindividual variation. While it was suggested that ligands access the buried active site of the enzyme from the membrane, no proof from unbiased simulations has been provided to support this hypothesis. Laboratory experiments fail to capture the access process which is suspected to influence binding kinetics. Here, we applied unbiased molecular dynamics (MD) simulations to investigate the access of ligands to wild-type CYP2D6, as well as the allelic variant CYP2D6*53. In multiple simulations, substrates accessed the active site of the enzyme from the protein-membrane interface to ultimately adopt a conformation that would allow a metabolic reaction. We propose the necessary steps for ligand access and the results suggest that the increased metabolic activity of CYP2D6*53 might be caused by a facilitated ligand uptake.

Rolapitant Is a Reversible Inhibitor of CYP2D6

Rolapitant [(Varubi), 5S,8S)-8-[[(1R)-1-[3,5 bis(trifluoromethyl phenyl]ethoxy]methyl]-8-phenyl-1,7-diazaspiro[4.5]decan-2-one] is a high-affinity NK1 receptor antagonist that was approved in September 2015 as a treatment for nausea and vomiting caused by chemotherapy. In vivo rolapitant moderately inhibits CYP2D6 for at least 7 days after one 180 mg dose. Due to the long inhibition time, we investigated rolapitant as a possible mechanism-based inactivator of CYP2D6. Rolapitant docked in the active site of CYP2D6 and displayed type I binding to CYP2D6 with a K _s value of 1.2 ± 0.4 µM. However, in NADPH-, time-, and concentration-dependent assays of CYP2D6 activity, no evidence for mechanism-based inactivation and no metabolites of rolapitant were observed. Stopped-flow binding studies yielded a k_on /k_off (K _d) value of 6.2 µM. The IC₅₀ value for rolapitant inhibition of CYP2D6 activity was 24 µM, suggesting that inhibition is not due to tight binding of rolapitant to CYP2D6. By Lineweaver-Burk analysis, rolapitant behaved as a mixed, reversible inhibitor. The K _i values of 20 and 34 µM were determined by Dixon analysis, with bufuralol and dextromethorphan as reporter substrates, respectively, and drug-drug interaction modeling did not predict the reported in vivo inhibition. The interaction of rolapitant with CYP2D6 was also examined in 1 microsecond molecular dynamics simulations. Rolapitant adopted multiple low-energy binding conformations near the active site, but at distances not consistent with metabolism. Given these findings, we do not see evidence that rolapitant is a mechanism-based inactivator. Moreover, the reversible inhibition of CYP2D6 by rolapitant may not fully account for the moderate inhibition described in vivo.