Molecular Dynamics Simulations Reveal Structural Differences among Allelic Variants of Membrane-Anchored Cytochrome P450 2D6

Cytochrome P450 2D6 (CYP2D6) is an enzyme that is involved in the metabolism of roughly 25% of all marketed drugs and therefore belongs to the most important enzymes in drug metabolism. CYP2D6 features a high degree of genetic polymorphism that can significantly affect the metabolic activity of an individual. In extreme cases, structural changes at the level of single amino acids can either increase its enzymatic activity abolishing the drug therapeutic effect or completely disable the enzyme and elevate drug plasma level potentially leading to adverse effects. In this study, starting from the crystal structure, we built a full-length membrane-anchored all-atom model of the wild-type CYP2D6 as well as five of its variants differing in the enzymatic activity. We validated our models with available experimental data and compared their structural properties with molecular dynamics simulations. The main focus of this study was to identify differences that could mechanistically explain the altered activity of the variants and improve our understanding of their functioning. We observed differences in the opening frequencies and minimal diameters of tunnels that connect the buried active site to the surrounding solvent environment. The variants CYP2D6*4 and CYP2D6*10 associated with missing or decreased activity showed less frequent opening of the tunnels compared to the wild-type. Both CYP2D6*10 and CYP2D6*17 showed a deprivation of an important ligand tunnel suggesting a feasible reason for their altered substrate specificity. Next, the altered fold at the N-terminal anchor region and the decreased active site volume caused by the amino acid mutations of the CYP2D6*4 variant offer an explanation for the absence of its metabolic activity. The mutations in CYP2D6*53 contributed to a significant enlargement of an important ligand tunnel and an extension of the active site cavity. This could explain the altered metabolic profile as well as the enhanced metabolic rates of this particular variant supporting its designation as a possible cause for the ultrarapid metabolizer phenotype. We believe these novel structural insights could advance the fields of personalized medicine and enzyme engineering. Furthermore, they could aid in guiding laboratory as well as computational experiments in the future.

Functional characterization of 50 CYP2D6 allelic variants by assessing primaquine 5-hydroxylation

Cytochrome P450 2D6 (CYP2D6) is responsible for the metabolic activation of primaquine, an antimalarial drug. CYP2D6 is genetically polymorphic, and these polymorphisms are associated with interindividual variations observed in the therapeutic efficacy of primaquine. To further understand this association, we performed in vitro enzymatic analyses of the wild-type CYP2D6.1 and 49 CYP2D6 allelic variants, which were expressed in 293FT cells, using primaquine as a substrate. The concentrations of CYP2D6 variant holoenzymes were measured by using carbon monoxide (CO)-reduced difference spectroscopy, and the wild type and 27 variants showed a peak at 450 nm. The kinetic parameters K_m, V_max, and intrinsic clearance (V_max/K_m) of primaquine 5-hydroxylation were characterized. The kinetic parameters of the wild type and 16 variants were measured, but the values for the remaining 33 variants could not be determined because of low metabolite concentrations. Among the variants, six (i.e., CYP2D6.17, .18, .35, .39, .53, and .70) showed significantly reduced intrinsic clearance compared with that of CYP2D6.1. Three-dimensional structural modeling analysis was performed to elucidate the mechanism of changes in the kinetics of CYP2D6 variants. Our findings provide insights into the allele-specific activity of CYP2D6 for primaquine, which could be clinically useful for malaria treatment and eradication efforts.

Oxidative metabolism of 5-methoxy-N,N-diisopropyltryptamine (Foxy) by human liver microsomes and recombinant cytochrome P450 enzymes

In vitro quantitative studies of the oxidative metabolism of (5-methoxy-N,N-diisopropyltryptamine, 5-MeO-DIPT, Foxy) were performed using human liver microsomal fractions and recombinant CYP enzymes and synthetic 5-MeO-DIPT metabolites. 5-MeO-DIPT was mainly oxidized to O-demethylated (5-OH-DIPT) and N-deisopropylated (5-MeO-IPT) metabolites in pooled human liver microsomes. In kinetic studies, 5-MeO-DIPT O-demethylation showed monophasic kinetics, whereas its N-deisopropylation showed triphasic kinetics. Among six recombinant CYP enzymes (CYP1A2, CYP2C8, CYP2C9, CYP2C19, CYP2D6 and CYP3A4) expressed in yeast or insect cells, only CYP2D6 exhibited 5-MeO-DIPT O-demethylase activity, while CYP1A2, CYP2C8, CYP2C9, CYP2C19 and CYP3A4 showed 5-MeO-DIPT N-deisopropylase activities. The apparent Km value of CYP2D6 was close to that for 5-MeO-DIPT O-demethylation, and the Km values of other CYP enzymes were similar to those of the low-Km (CYP2C19), intermediate-Km (CYP1A2, CYP2C8 and CYP3A4) and high-Km phases (CYP2C9), respectively, for N-deisopropylation in human liver microsomes. In inhibition studies, quinidine (1 microM), an inhibitor of CYP2D6, almost completely inhibited human liver microsomal 5-MeO-DIPT O-demethylation at a substrate concentration of 10 microM. Furafylline, a CYP1A2 inhibitor, quercetin, a CYP2C8 inhibitor, sulfaphenazole, a CYP2C9 inhibitor and ketoconazole, a CYP3A4 inihibitor (5 microM each) suppressed about 60%, 45%, 15% and 40%, respectively, of 5-MeO-DIPT N-deisopropylation at 50 microM substrate. In contrast, omeprazole (10 microM), a CYP2C19 inhibitor, suppressed only 10% of N-deisopropylation by human liver microsomes, whereas at the same concentration the inhibitor suppressed the reaction by recombinant CYP2C19 almost completely. These results indicate that CYP2D6 is the major 5-MeO-DIPT O-demethylase, and CYP1A2, CYP2C8 and CYP3A4 are the major 5-MeO-DIPT N-deisopropylase enzymes in the human liver.

Roles of phenylalanine at position 120 and glutamic acid at position 222 in the oxidation of chiral substrates by cytochrome P450 2D6

The roles of Phe-120 and Glu-222 in the oxidation of chiral substrates bunitrolol (BTL) and bufuralol (BF) by CYP2D6 are discussed. Wild-type CYP2D6 (CYP2D6-WT) oxidized BTL to 4-hydroxybunitrolol (4-OH-BTL) with substrate enantioselectivity of (R)-(+)-BTL > (S)-(-)-BTL. The same enzyme converted BF into 1''-hydroxybufuralol with substrate enantioselectivity of (R)-BF >> (S)-BF and metabolite diastereoselectivity of (1''R)-OH < (1''S)-OH. The substitution of Phe-120 by alanine markedly increased the apparent K(m) and V(max) values for enantiomeric BTL 4-hydroxylation by CYP2D6. In contrast, the same substitution caused an increase only in V(max) values of (S)-BF 1''-hydroxylation without changing apparent K(m) values, while kinetic parameters (K(m) and V(max) values) for (R)-BF 1''-hydroxylation remained unchanged. Furthermore, the substitution of Glu-222 as well as Glu-216 by alanine remarkably decreased both the apparent K(m) and V(max) values without changing substrate enantioselectivity or metabolite diastereoselectivity. A computer-assisted simulation study using energy minimization and molecular dynamics techniques indicated that the hydrophobic interaction of an aromatic moiety of the substrate with Phe-120 and the ionic interaction of a basic nitrogen atom of the substrate with Glu-222 in combination with Glu-216 play important roles in the binding of BF and BTL by CYP2D6 and the orientation of these substrates in the active-site cavity. This modeling yielded a convincing explanation for the reversal of substrate enantioselectivity in BTL 4-hydroxylation between CYP2D6-WT and CYP2D6-V374M having methionine in place of Val-374, which supports the validity of this modeling.

The effect of dimethyl sulfoxide on the function of cytochrome P450 2D6 in HepG2 cells upon the co-expression with NADPH-cytochrome P450 reductase

HepG2 cells, a human hepatoma cell line, stably expressing NADPH-cytochrome P450 reductase (OR) and/or cytochrome P450 2D6 wild-type (CYP2D6-WT) or its variants (Pro34Ser, Gly42Arg, Arg296Cys and Ser486Thr) were established in the present study. The cultivation of HepG2 cells expressing CYP2D6-WT in the culture medium containing dimethyl sulfoxide (DMSO, 0.1% of final concentration) markedly increased the bufuralol (BF) 1''-hydroxylase activity compared with that of control cells when cultivated without DMSO. A similar effect was also observed in HepG2 cells stably expressing CYP2D6 and OR. The addition of hemin in place of DMSO to the culture medium resulted in no increase in the enzyme activity. Western blot analysis revealed that the levels of CYP2D6 protein were similar between DMSO-treated and non-treated HepG2 cells regardless of OR expression. Spectrophotometric analysis of reduced carbon monoxide-difference spectra of HepG2 cells expressing CYP2D6-WT and/or OR demonstrated that the addition of DMSO increased the peak height of functional CYP2D6 at 450 nm. These results suggest that the increase in CYP2D6 activity is attributable to the radical-scavenging effect of DMSO. The HepG2 cell lines stably expressing OR and CYP2D6 or its variants in combination with DMSO treatment may be useful for screening the cytotoxicity of chemical compounds which undergo oxidation by CYP2D6.

Change in enantioselectivity in bufuralol 1?-hydroxylation by the substitution of phenylalanine-120 by alanine in cytochrome P450 2D6

The functional roles of phenylalanine at position 120 in drug oxidation by cytochrome P450 2D6 (CYP2D6) were examined using a yeast cell expression system and bufuralol (BF) enantiomers as a chiral substrate. Two mutated cDNAs, one encoding a CYP2D6 mutant having alanine instead of Phe-120 (F120A) and another encoding a mutant having alanine instead of Glu-222 (E222A), were prepared by site-directed mutagenesis and transformed into yeast cells via pGYRI vectors. The enantiomeric BF 1''-hydroxylase activities of the mutants were compared with those of the wild type. When enantiomeric BF 1''-hydroxylase activities at a substrate concentration of 100 microM were compared, the CYP2D6 wild type showed substrate enantioselectivity of (R-BF >> S-BF) and the F120A mutant exhibited substrate enantioselectivity of (R-BF < or = S-BF), whereas the product diastereoselectivity of (1''R-OH-BF << 1''-S-OH-BF) was similar between the wild type and the mutant. The activities of the other mutant (E222A) were much lower than those of the wild type and the F120A mutant, while its substrate enantioselectivity and product diastereoselectivity were the same as those of the wild type. The kinetics demonstrated that apparent K(m) values were similar among the recombinant enzymes, and V(max) values clearly reflected the selectivity described above. These results indicate that Phe-120 has a key role in the enantioselective BF 1''-hydroxylation by CYP2D6.