Influence of phenylalanine 120 on cytochrome P450 2D6 catalytic selectivity and regiospecificity: crucial role in 7-methoxy-4-(aminomethyl)-coumarin metabolism

Removal of BPA by Pseudomonas asiatica P1: Synergistic response mechanism of toxicity resistance and biodegradation

Bisphenol A (BPA) is a globally concerning toxic pollutant, and microbial degradation is considered an effective method to treat BPA contamination. However, the inherent microbial toxicity of BPA is often overlooked, particularly the microbial mechanisms of resistance and detoxification against BPA. This study found that under the toxic stress of BPA, cbb3-type cytochrome c oxidase (cbb3-Cox) in the cells of Pseudomonas asiatica P1 (P. asiatica P1) was the first to resist the toxicity. Genes such as ccoNOQPG showed significant upregulation with an average log2FC value of 3.56. Subsequently, genes that are related to metal ion binding, transport, and DNA repair were upregulated in the middle to later phase, which enhanced the metabolic functions of the strains and induced strain mutations to assist P. asiatica P1 in resisting the BPA toxicity. Meanwhile, three potential BPA degradation genes were identified, among which sdrP1 was crucial to the BPA degradation and detoxification. After genetic recombination, sdrP1 achieved a degradation rate of 92.52 % for BPA. Furthermore, through various methods such as alkyl interactions, sdrP1 exhibited oxidation and demethylation to form lower toxic intermediate products and complete the biological detoxification of BPA. This study provides a systematic analysis of the toxicity resistance, biodegradation, and detoxification processes in bacterial BPA removal, refines the mechanism of BPA biodegradation and contributes to a more comprehensive and systematic understanding of the overall process of microbial removal of toxic pollutants.

Molecular Mechanisms Underlying Single Nucleotide Polymorphism-Induced Reactivity Decrease in CYP2D6

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.

Coumarin-Synthetic Methodologies, Pharmacology, and Application as Natural Fluorophore

Coumarins are secondary metabolites made up of benzene and α-pyrone rings fused together that can potentially treat various ailments, including cancer, metabolic, and degenerative disorders. Coumarins are a diverse category of both naturally occurring as well as synthesized compounds with numerous biological and therapeutic properties. Coumarins as fluorophores play a key role in fluorescent labeling of biomolecules, metal ion detection, microenvironment polarity detection, and pH detection. This review provides a detailed insight into the characteristics of coumarins as well as their biosynthesis in plants and metabolic pathways. Various synthetic strategies for coumarin core involving both conventional and green methods have been discussed comparing advantages and disadvantages of each method. Conventional methods discussed are Pechmann, Knoevenagel, Perkin, Wittig, Kostanecki, Buchwald-Hartwig, and metal-induced coupling reactions such as Heck and Suzuki, as well as green approaches involving microwave or ultrasound energy. Various pharmacological applications of coumarin derivatives are discussed in detail. The structural features and conditions responsible for influencing the fluorescence of coumarin core are also elaborated.

Elucidating the Mechanism of Metabolism of Cannabichromene by Human Cytochrome P450s

Cannabichromene (CBC) is a nonpsychoactive phytocannabinoid well-known for its wide-ranging health advantages. However, there is limited knowledge regarding its human metabolism following CBC consumption. This research aimed to explore the metabolic pathways of CBC by various human liver cytochrome P450 (CYP) enzymes and support the outcomes using in vivo data from mice. The results unveiled two principal CBC metabolites generated by CYPs: 8'-hydroxy-CBC and 6',7'-epoxy-CBC, along with a minor quantity of 1″-hydroxy-CBC. Notably, among the examined CYPs, CYP2C9 demonstrated the highest efficiency in producing these metabolites. Moreover, through a molecular dynamics simulation spanning 1 μs, it was observed that CBC attains stability at the active site of CYP2J2 by forming hydrogen bonds with I487 and N379, facilitated by water molecules, which specifically promotes the hydroxy metabolite's formation. Additionally, the presence of cytochrome P450 reductase (CPR) amplified CBC's binding affinity to CYPs, particularly with CYP2C8 and CYP3A4. Furthermore, the metabolites derived from CBC reduced cytokine levels, such as IL6 and NO, by approximately 50% in microglia cells. This investigation offers valuable insights into the biotransformation of CBC, underscoring the physiological importance and the potential significance of these metabolites.

"Commandeuring" Xenobiotic Metabolism: Advances in Understanding Xenobiotic Metabolism

The understanding of how exogenous chemicals (xenobiotics) are metabolized, distributed, and eliminated is critical to determine the impact of the chemical and its metabolites to the (human) organism. This is part of the research and educational discipline ADMET (absorption, distribution, metabolism, elimination, and toxicity). Here, we review the work of Jan Commandeur and colleagues who have not only made a significant impact in understanding of phase I and phase II metabolism of several important compounds but also contributed greatly to the development of experimental techniques for the study of xenobiotic metabolism. Jan Commandeur’s work has covered a broad area of research, such as the development of online screening methodologies, the use of a combination of enzyme mutagenesis and molecular modeling for structure–activity relationship (SAR) studies, and the development of novel probe substrates. This work is the bedrock of current activities and brings the field closer to personalized (cohort-based) pharmacology, toxicology, and hazard/risk assessment.

D3PM: a comprehensive database for protein motions ranging from residue to domain

Background

Knowledge of protein motions is significant to understand its functions. While currently available databases for protein motions are mostly focused on overall domain motions, little attention is paid on local residue motions. Albeit with relatively small scale, the local residue motions, especially those residues in binding pockets, may play crucial roles in protein functioning and ligands binding.

Results

A comprehensive protein motion database, namely D3PM, was constructed in this study to facilitate the analysis of protein motions. The protein motions in the D3PM range from overall structural changes of macromolecule to local flip motions of binding pocket residues. Currently, the D3PM has collected 7679 proteins with overall motions and 3513 proteins with pocket residue motions. The motion patterns are classified into 4 types of overall structural changes and 5 types of pocket residue motions. Impressively, we found that less than 15% of protein pairs have obvious overall conformational adaptations induced by ligand binding, while more than 50% of protein pairs have significant structural changes in ligand binding sites, indicating that ligand-induced conformational changes are drastic and mainly confined around ligand binding sites. Based on the residue preference in binding pocket, we classified amino acids into “pocketphilic” and “pocketphobic” residues, which should be helpful for pocket prediction and drug design.

Conclusion

D3PM is a comprehensive database about protein motions ranging from residue to domain, which should be useful for exploring diverse protein motions and for understanding protein function and drug design. The D3PM is available on www.d3pharma.com/D3PM/index.php.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-022-04595-0.

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.

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

Metabolic phenotype can be affected by multiple factors, including allelic variation and interactions with inhibitors. Human CYP2D6 is responsible for approximately 20% of cytochrome P450-mediated drug metabolism but consists of more than 100 known variants; several variants are commonly found in the population, whereas others are quite rare. Four CYP2D6 allelic variants-three with a series of mutations distal to the active site (*34, *17-2, *17-3) and one ultra-metabolizer with mutations near the active site (*53), along with reference *1 and an active site mutant of *1 (Thr309Ala)-were expressed, purified, and studied for interactions with the typical substrates dextromethorphan and bufuralol and the inactivator SCH 66712. We found that *34, *17-2, and *17-3 displayed reduced enzyme activity and NADPH coupling while producing the same metabolites as *1, suggesting a possible role for Arg296 in NADPH coupling. A higher-activity variant, *53, displayed similar NADPH coupling to *1 but was less susceptible to inactivation by SCH 66712. The Thr309Ala mutant showed similar activity to that of *1 but with greatly reduced NADPH coupling. Overall, these results suggest that kinetic and metabolic analysis of individual CYP2D6 variants is required to understand their possible contributions to variable drug response and the complexity of personalized medicine.

Engineering a self-sufficient Mycobacterium tuberculosis CYP130 by gene fusion with the reductase-domain of CYP102A1 from Bacillus megaterium

CYP130 belongs to the subset of cytochrome P450s from Mycobacterium tuberculosis (Mtb) that have been structurally characterized. Despite several efforts for its functional characterization, CYP130 is still considered an orphan enzyme for which no endogenous or exogenous substrate has been identified. In addition, functional redox-partners for CYP130 have not been clearly established yet, hampering the elucidation of its physiological role. In the present study, a catalytically active fusion protein involving CYP130 and the NADPH reductase-domain of CYP102A1 from Bacillus megaterium was created. By screening a panel of known substrates of human P450s, dextromethorphan N-demethylation was identified as a reaction catalyzed by CYP130. The fusion enzyme showed higher catalytic activity, when compared to CYP130 reconstituted with a selection of non-native redox-partners. Molecular dynamics simulation studies based on the crystal structure of CYP130 revealed two primary docking poses of dextromethorphan within the active site consistent with the experimentally observed N-demethylation reaction during the entire molecular dynamics simulation. The dextromethorphan N-demethylation reaction was strongly inhibited by azole-drugs and maybe applied to identify mechanism-based inhibitors of CYP130. Furthermore, the present active CYP130-fusion protein may facilitate the identification of endogenous substrates from Mtb.

Functional characterization of wild-type and 49 CYP2D6 allelic variants for N-desmethyltamoxifen 4-hydroxylation activity

Genetic variations in cytochrome P450 2D6 (CYP2D6) contribute to interindividual variability in the metabolism of clinically used drugs, e.g., tamoxifen. CYP2D6 is genetically polymorphic and is associated with large interindividual variations in therapeutic efficacy and drug toxicity. In this study, we performed an in vitro analysis of 50 allelic variants of CYP2D6 proteins. Wild-type CYP2D6.1 and 49 variants were transiently expressed in COS-7 cells, and the enzymatic activities of the CYP2D6 variants were characterized using N-desmethyltamoxifen as a substrate. The kinetic parameters K(m), V(max), and intrinsic clearance (V(max)/K(m)) of N-desmethyltamoxifen 4-hydroxylation were determined. Among the 50 CYP2D6 variants, the kinetic parameters for N-desmethyltamoxifen 4-hydroxylation were determined for 20 CYP2D6 variants. On the other hand, the kinetic parameters of 30 CYP2D6 variants could not be determined because the amount of metabolite produced was at or below the detection limit at the lower substrate concentrations. Among them, 8 variants, i.e., CYP2D6.2, .9, .26, .28, .32, .43, .45, and .70, showed decreased intrinsic clearance at <50% of CYP2D6.1. The comprehensive in vitro assessment of CYP2D6 variants provides novel insights into allele-specific activity towards tamoxifen and may be valuable when interpreting in vivo studies.

Functional characterization of a first avian cytochrome P450 of the CYP2D subfamily (CYP2D49)

The CYP2D family members are instrumental in the metabolism of 20-25% of commonly prescribed drugs. Although many CYP2D isoforms have been well characterized in other animal models, research concerning the chicken CYP2Ds is limited. In this study, a cDNA encoding a novel CYP2D enzyme (CYP2D49) was cloned from the chicken liver for the first time. The CYP2D49 cDNA contained an open reading frame of 502 amino acids that shared 52%-57% identities with other CYP2Ds. The gene structure and neighboring genes of CYP2D49 are conserved and similar to those of human CYP2D6. Additionally, similar to human CYP2D6, CYP2D49 is un-inducible in the liver and expressed predominantly in the liver, kidney and small intestine, with detectable levels in several other tissues. Metabolic assays of the CYP2D49 protein heterologously expressed in E. coli and Hela cells indicated that CYP2D49 metabolized the human CYP2D6 substrate, bufuralol, but not debrisoquine. Moreover, quinidine, a potent inhibitor of human CYP2D6, only inhibited the bufuralol 1'-hydroxylation activity of CYP2D49 to a negligible degree. All these results indicated that CYP2D49 had functional characteristics similar to those of human CYP2D6 but measurably differed in the debrisoquine 4'-hydroxylation and quinidine inhibitory profile. Further structure-function investigations that employed site-directed mutagenesis and circular dichroism spectroscopy identified the importance of Val-126, Glu-222, Asp-306, Phe-486 and Phe-488 in keeping the enzymatic activity of CYP2D49 toward bufuralol as well as the importance of Asp-306, Phe-486 and Phe-488 in maintaining the conformation of CYP2D49 protein. The current study is only the first step in characterizing the metabolic mechanism of CYP2D49; further studies are still required.

Molecular docking and enzyme kinetic studies of dihydrotanshinone on metabolism of a model CYP2D6 probe substrate in human liver microsomes

The effects of Danshen and its active components (tanshinone I, tanshinone IIA, dihydrotanshinone and cryptotanshinone) on CYP2D6 activity was investigated by measuring the metabolism of a model CYP2D6 probe substrate, dextromethorphan to dextrorphan in human pooled liver microsomes. The ethanolic extract of crude Danshen (6.25-100 μg/ml) decreased dextromethorphan O-demethylation in vitro (IC(50)=23.3 μg/ml) and the water extract of crude Danshen (0.0625-1 mg/ml) showed no inhibition. A commercially available Danshen pill (31.25-500 μg/ml) also decreased CYP2D6 activity (IC(50)=265.8 μg/ml). Among the tanshinones, only dihydrotanshinone significantly inhibited CYP2D6 activity (IC(50)=35.4 μM), compared to quinidine, a specific CYP2D6 inhibitor (IC(50)=0.9 μM). Crytotanshinone, tanshinone I and tanshinone IIA produced weak inhibition, with IC(20) of 40.8 μM, 16.5 μM and 61.4 μM, respectively. Water soluble components such as salvianolic acid B and danshensu did not affect CYP2D6-mediated metabolism. Enzyme kinetics studies showed that inhibition of CYP2D6 activity by the ethanolic extract of crude Danshen and dihydrotanshinone was concentration-dependent, with K(i) values of 4.23 μg/ml and 2.53 μM, respectively, compared to quinidine, K(i)=0.41 μM. Molecular docking study confirmed that dihydrotanshinone and tanshinone I interacted with the Phe120 amino acid residue in the active cavity of CYP2D6 through Pi-Pi interaction, but did not interact with Glu216 and Asp301, the key residues for substrate binding. The logarithm of free binding energy of dihydrotanshinone (-7.6 kcal/mol) to Phe120 was comparable to quinidine (-7.0 kcal/mol) but greater than tanshinone I (-5.4 kcal/mol), indicating dihydrotanshinone has similar affinity to quinidine in binding to the catalytic site on CYP2D6.

Using a homology model of cytochrome P450 2D6 to predict substrate site of metabolism

CYP2D6 is an important enzyme that is involved in first pass metabolism and is responsible for metabolizing ~25% of currently marketed drugs. A homology model of CYP2D6 was built using X-ray structures of ligand-bound CYP2C5 complexes as templates. This homology model was used in docking studies to rationalize and predict the site of metabolism of known CYP2D6 substrates. While the homology model was generally found to be in good agreement with the recently solved apo (ligand-free) X-ray structure of CYP2D6, significant differences between the structures were observed in the B' and F-G helical region. These structural differences are similar to those observed between ligand-free and ligand-bound structures of other CYPs and suggest that these conformational changes result from induced-fit adaptations upon ligand binding. By docking to the homology model using Glide, it was possible to identify the correct site of metabolism for a set of 16 CYP2D6 substrates 85% of the time when the 5 top scoring poses were examined. On the other hand, docking to the apo CYP2D6 X-ray structure led to a loss in accuracy in predicting the sites of metabolism for many of the CYP2D6 substrates considered in this study. These results demonstrate the importance of describing substrate-induced conformational changes that occur upon binding. The best results were obtained using Glide SP with van der Waals scaling set to 0.8 for both the receptor and ligand atoms. A discussion of putative binding modes that explain the distribution of metabolic sites for substrates, as well as a relationship between the number of metabolic sites and substrate size, are also presented. In addition, analysis of these binding modes enabled us to rationalize the typical hydroxylation and O-demethylation reactions catalyzed by CYP2D6 as well as the less common N-dealkylation.

New insights into the structural characteristics and functional relevance of the human cytochrome P450 2D6 enzyme

To date, the crystal structures of at least 12 human CYPs (1A2, 2A6, 2A13, 2C8, 2C9, 2D6, 2E1, 2R1, 3A4, 7A1, 8A1, and 46A1) have been determined. CYP2D6 accounts for only a small percentage of all hepatic CYPs (< 2%), but it metabolizes approximately 25% of clinically used drugs with significant polymorphisms. CYP2D6 also metabolizes procarcinogens and neurotoxins, such as 1-methyl-4-phenyl-1,2,3,6-tetrahydropyridine, 1,2,3,4-tetrahydroquinoline, and indolealkylamines. Moreover, the enzyme utilizes hydroxytryptamines and neurosteroids as endogenous substrates. Typical CYP2D6 substrates are usually lipophilic bases with an aromatic ring and a nitrogen atom, which can be protonated at physiological pH. Substrate binding is generally followed by oxidation (5-7 A) from the proposed nitrogen-Asp301 interaction. A number of homology models have been constructed to explore the structural features of CYP2D6, while antibody studies also provide useful structural information. Site-directed mutagenesis studies have demonstrated that Glu216, Asp301, Phe120, Phe481, and Phe483 play important roles in determining the binding of ligands to CYP2D6. The structure of human CYP2D6 has been recently determined and shows the characteristic CYP fold observed for other members of the CYP superfamily. The lengths and orientations of the individual secondary structural elements in the CYP2D6 structure are similar to those seen in other human CYP2 members, such as CYP2C9 and 2C8. The 2D6 structure has a well-defined active-site cavity located above the heme group with a volume of approximately 540 A(3), which is larger than equivalent cavities in CYP2A6 (260 A(3)), 1A2 (375 A(3)), and 2E1 (190 A(3)), but smaller than those in CYP3A4 (1385 A(3)) and 2C8 (1438 A(3)). Further studies are required to delineate the molecular mechanisms involved in CYP2D6 ligand interactions and their implications for drug development and clinical practice.

Cytochrome P450 125 (CYP125) catalyses C26-hydroxylation to initiate sterol side-chain degradation in Rhodococcus jostii RHA1

The cyp125 gene of Rhodococcus jostii RHA1 was previously found to be highly upregulated during growth on cholesterol and the orthologue in Mycobacterium tuberculosis (rv3545c) has been implicated in pathogenesis. Here we show that cyp125 is essential for R. jostii RHA1 to grow on 3-hydroxysterols such as cholesterol, but not on 3-oxo sterol derivatives, and that CYP125 performs an obligate first step in cholesterol degradation. The involvement of cyp125 in sterol side-chain degradation was confirmed by disrupting the homologous gene in Rhodococcus rhodochrous RG32, a strain that selectively degrades the cholesterol side-chain. The RG32 Omega cyp125 mutant failed to transform the side-chain of cholesterol, but degraded that of 5-cholestene-26-oic acid-3beta-ol, a cholesterol catabolite. Spectral analysis revealed that while purified ferric CYP125(RHA1) was < 10% in the low-spin state, cholesterol (K(D)(app) = 0.20 +/- 0.08 microM), 5 alpha-cholestanol (K(D)(app) = 0.15 +/- 0.03 microM) and 4-cholestene-3-one (K(D)(app) = 0.20 +/- 0.03 microM) further reduced the low spin character of the haem iron consistent with substrate binding. Our data indicate that CYP125 is involved in steroid C26-carboxylic acid formation, catalysing the oxidation of C26 either to the corresponding carboxylic acid or to an intermediate state.

Polymorphism of human cytochrome P450 enzymes and its clinical impact

Pharmacogenetics is the study of how interindividual variations in the DNA sequence of specific genes affect drug response. This article highlights current pharmacogenetic knowledge on important human drug-metabolizing cytochrome P450s (CYPs) to understand the large interindividual variability in drug clearance and responses in clinical practice. The human CYP superfamily contains 57 functional genes and 58 pseudogenes, with members of the 1, 2, and 3 families playing an important role in the metabolism of therapeutic drugs, other xenobiotics, and some endogenous compounds. Polymorphisms in the CYP family may have had the most impact on the fate of therapeutic drugs. CYP2D6, 2C19, and 2C9 polymorphisms account for the most frequent variations in phase I metabolism of drugs, since almost 80% of drugs in use today are metabolized by these enzymes. Approximately 5-14% of Caucasians, 0-5% Africans, and 0-1% of Asians lack CYP2D6 activity, and these individuals are known as poor metabolizers. CYP2C9 is another clinically significant enzyme that demonstrates multiple genetic variants with a potentially functional impact on the efficacy and adverse effects of drugs that are mainly eliminated by this enzyme. Studies into the CYP2C9 polymorphism have highlighted the importance of the CYP2C9*2 and *3 alleles. Extensive polymorphism also occurs in other CYP genes, such as CYP1A1, 2A6, 2A13, 2C8, 3A4, and 3A5. Since several of these CYPs (e.g., CYP1A1 and 1A2) play a role in the bioactivation of many procarcinogens, polymorphisms of these enzymes may contribute to the variable susceptibility to carcinogenesis. The distribution of the common variant alleles of CYP genes varies among different ethnic populations. Pharmacogenetics has the potential to achieve optimal quality use of medicines, and to improve the efficacy and safety of both prospective and currently available drugs. Further studies are warranted to explore the gene-dose, gene-concentration, and gene-response relationships for these important drug-metabolizing CYPs.

Metabolism ofN-substituted 7-methoxy-4-(aminomethyl) -coumarins by cytochrome P450 2D6 mutants and the indication of additional substrate interaction points

Previous studies have shown the critical roles residues F120 and F483 play in the oxidative metabolism of 7-methoxy-4-(aminomethyl)-coumarin (MAMC) by cytochrome P450 2D6 (CYP2D6). In the present study, a series of N-alkyl-7-methoxy-4-(aminomethyl)-coumarins (MAMC analogues) were used as substrates for the F120A and F483A mutants in order to probe the CYP2D6 active site. The F120A and F483A mutants of CYP2D6 displayed significant activity towards the MAMC analogues. Automated docking studies of the MAMC analogues in a CYP2D6 homology model suggested a distal hydrophobic active site binding cleft for the substrate N-alkyl chains, consisting of the residues L213 and V308.

Molecular modeling of human cytochrome P450 2W1 and its interactions with substrates

The human cytochrome P450 2W1 (CYP2W1) was categorized into the so-called "orphan" CYPs because of its unknown enzymatic function. However, recent studies showed that the recombinant CYP2W1 exhibited broad catalytic activity towards several chemicals. Furthermore, this enzyme was selectively expressed in some forms of cancers, whereas a very low expression was found in human normal issues. These render CYP2W1 as a potential drug target for cancer therapy. At present, however, little information is available on the active site topology and the substrate binding modes of CYP2W1. In this study, the three-dimensional model of CYP2W1 was constructed using the homology modeling method. Two known substrates, benzphetamine and indole, were then docked into the active site, and refined by molecular dynamics simulations. The interaction energy between the substrates and the enzyme was calculated and analyzed by using the MM-GBSA method. The results indicated that the constructed CYP2W1 model can account for the regioselectivity of this enzyme towards the known substrates and van der Waals interactions were the driving force for the substrate binding. Several key residues were identified to be responsible for the binding of indole and benzphetamine with CYP2W1. These findings provide useful information for the detailed characterization of the biological roles of CYP2W1 and structure-based drug design of this enzyme.

Identification of selectivity-determining residues in cytochrome P450 monooxygenases: a systematic analysis of the substrate recognition site 5

The large and diverse family of cytochrome P450 monooxygenases (CYPs) was systematically analyzed to identify selectivity- and specificity-determining residues in the substrate recognition site 5, which is located in close vicinity to the heme center. A positively charged heme-interacting residue was identified in the structures of 29 monooxygenases and in 97.7% of the 6379 CYP sequences investigated here. This heme-interacting residue restricts the conformation of the substrate recognition site 5 and is preferentially located at position 10 or 11 after the conserved ExxR motif (in 94.4% of the sequences), in 3.3% of the sequences at position 9 or 12. As a result, a classification by the position of the heme-interacting residue allows to predict residues that are closest to the heme center and restrict its accessibility. In 98.4% of all CYP sequences a preferentially hydrophobic residue is located at position 5 after the ExxR motif that is predicted to point close to the heme center. Replacing this residue by hydrophobic residues of different size has been shown to change substrate specificity and regioselectivity for CYPs of different superfamilies. Twenty-seven percent of all CYPs are predicted to contain a second selectivity-determining residue at position 9 after the ExxR motif that can be identified by the pattern EXXR-X(7)-{P}-x-P-[HKR].

Cytochrome P450 bio-affinity detection coupled to gradient HPLC: on-line screening of affinities to cytochrome P4501A2 and 2D6

Here we describe novel on-line human CYP1A2 and CYP2D6 Enzyme Affinity Detection (EAD) systems coupled to gradient HPLC. The use of the systems lies in the detection of individual inhibitory ligands in mixtures (e.g. metabolic mixtures or herbal extracts) towards two relevant drug metabolizing human CYPs. The systems can rapidly detect individual compounds in mixtures with affinities to CYP1A2 or 2D6. The HPLC-EAD systems were first evaluated and validated in flow injection analysis mode. IC50 values of known ligands for both CYPs, tested both in flow injection and in HPLC mode, were well comparable with those measured in microplate reader formats. Both EAD systems were also connected to gradient HPLC and used to screen known compound mixtures for the presence of CYP1A2 and 2D6 inhibitors. Finally, the on-line CYP2D6 EAD system was used to screen for the inhibitory activities of stereoisomers of a mixture of five methylenedioxy-alkylamphetamines (XTC analogs) on a chiral analytical column.

Free energies of binding of R- and S-propranolol to wild-type and F483A mutant cytochrome P450 2D6 from molecular dynamics simulations

Detailed molecular dynamics (MD) simulations have been performed to reproduce and rationalize the experimental finding that the F483A mutant of CYP2D6 has lower affinity for R-propranolol than for S-propranolol. Wild-type (WT) CYP2D6 does not show this stereospecificity. Four different approaches to calculate the free energy differences have been investigated and were compared to the experimental binding data. From the differences between calculations based on forward and backward processes and the closure of thermodynamic cycles, it was clear that not all simulations converged sufficiently. The approach that calculates the free energies of exchanging R-propranolol with S-propranolol in the F483A mutant relative to the exchange free energy in WT CYP2D6 accurately reproduced the experimental binding data. Careful inspection of the end-points of the MD simulations involved in this approach, allowed for a molecular interpretation of the observed differences.

Heterotropic and homotropic cooperativity by a drug-metabolising mutant of cytochrome P450 BM3

Recently, we described a triple mutant of the bacterial cytochrome P450 BM3 as the first mutant with affinity for drug-like compounds. In this paper, we show that this mutant, but not wild-type BM3, is able to metabolise testosterone and several drug-like molecules such as amodiaquine, dextromethorphan, acetaminophen, and 3,4-methylenedioxymethylamphetamine that are known substrates of human P450s. Interestingly, the metabolism of 3,4-methylenedioxymethylamphetamine and acetaminophen could be stimulated up to 70-fold by the addition of caffeine, a known activator of rat P450 3A2. With testosterone metabolism, homotropic cooperativity was observed. This shows that heterotropic and homotropic cooperativity, known to occur in the P450 3A family, can also take place in BM3. BM3 therefore can be used as a model system to study atypical kinetics in mammalian P450s. Second, this study shows that BM3 can be engineered to a drug-metabolising enzyme, making it a promising candidate to use as biocatalyst in drug discovery and synthesis.

Catalytic site prediction and virtual screening of cytochrome P450 2D6 substrates by consideration of water and rescoring in automated docking

Automated docking strategies successfully applied to binding mode predictions of ligands in Cyt P450 crystal structures in an earlier study (de Graaf et al. J. Med. Chem. 2005, 7, 2308-2318) were used for the catalytic site prediction (CSP) of 65 substrates in a CYP2D6 homology model. The consideration of water molecules at predicted positions in the active site and the rescoring of pooled docking poses from four different docking programs (AutoDock, FlexX, GOLD-Goldscore, and GOLD-Chemscore) with the SCORE scoring function enabled the successful prediction of experimentally reported sites of catalysis of more than 80% of the substrates. Three docking algorithms (FlexX, GOLD-Goldscore, and GOLD-Chemscore) were subsequently used in combination with six scoring functions (Chemscore, DOCK, FlexX, GOLD, PMF, and SCORE) to assess the ability of docking-based virtual screening methods to prioritize known CYP2D6 substrates seeded into a drug-like chemical database (in the absence and presence of active-site water molecules). Finally, the optimal docking strategy in terms of virtual screening accuracy, GOLD-Chemscore with the consideration of active-site water (60% of known substrates recovered in the top 5% of the ranked drug-like database), was verified experimentally; it was successfully used to identify high-affinity CYP2D6 ligands among a larger proprietary database.

IN SILICO PREDICTION OF DRUG BINDING TO CYP2D6: IDENTIFICATION OF A NEW METABOLITE OF METOCLOPRAMIDE

Patients with cancer often take many different classes of drugs to treat the effects of their malignancy and the side effects of treatment, as well as their comorbidities. The potential for drug-drug interactions that may affect the efficacy of anticancer treatment is high, and a major source of such interactions is competition for the drug-metabolizing enzymes, cytochromes P450 (P450s). We have examined a series of 20 drugs commonly prescribed to cancer patients to look for potential interactions via CYP2D6. We used a homology model of CYP2D6, together with molecular docking techniques, to perform an in silico screen for binding to CYP2D6. Experimental IC50 values were determined for these compounds and compared with the model predictions to reveal a correlation with a regression coefficient of r2= 0.61. Importantly, the docked conformation of the commonly prescribed antiemetic metoclopramide predicted a new site of metabolism that was further investigated through in vitro analysis with recombinant CYP2D6. An aromatic N-hydroxy metabolite of metoclopramide, consistent with predictions from our modeling studies, was identified by high-performance liquid chromatography/mass spectrometry. This metabolite was found to represent a major product of metabolism in human liver microsomes, and CYP2D6 was identified as the main P450 isoform responsible for catalyzing its formation. In view of the prevalence of interindividual variation in the CYP2D6 genotype and phenotype, we suggest that those experiencing adverse reactions with metoclopramide, e.g., extrapyramidal syndrome, are likely to have a particular CYP2D6 genotype/phenotype. This warrants further investigation.

Functional interaction of nitrogenous organic bases with cytochrome P450: A critical assessment and update of substrate features and predicted key active-site elements steering the access, binding, and orientation of amines

The widespread use of nitrogenous organic bases as environmental chemicals, food additives, and clinically important drugs necessitates precise knowledge about the molecular principles governing biotransformation of this category of substrates. In this regard, analysis of the topological background of complex formation between amines and P450s, acting as major catalysts in C- and N-oxidative attack, is of paramount importance. Thus, progress in collaborative investigations, combining physico-chemical techniques with chemical-modification as well as genetic engineering experiments, enables substantiation of hypothetical work resulting from the design of pharmacophores or homology modelling of P450s. Based on a general, CYP2D6-related construct, the majority of prospective amine-docking residues was found to cluster near the distal heme face in the six known SRSs, made up by the highly variant helices B', F and G as well as the N-terminal portion of helix C and certain beta-structures. Most of the contact sites examined show a frequency of conservation < 20%, hinting at the requirement of some degree of conformational versatility, while a limited number of amino acids exhibiting a higher level of conservation reside close to the heme core. Some key determinants may have a dual role in amine binding and/or maintenance of protein integrity. Importantly, a series of non-SRS elements are likely to be operative via long-range effects. While hydrophobic mechanisms appear to dominate orientation of the nitrogenous compounds toward the iron-oxene species, polar residues seem to foster binding events through H-bonding or salt-bridge formation. Careful uncovering of structure-function relationships in amine-enzyme association together with recently developed unsupervised machine learning approaches will be helpful in both tailoring of novel amine-type drugs and early elimination of potentially toxic or mutagenic candidates. Also, chimeragenesis might serve in the construction of more efficient P450s for activation of amine drugs and/or bioremediation.

Binding of bufuralol, dextromethorphan, and 3,4-methylenedioxymethylamphetamine to wild-type and F120A mutant cytochrome P450 2D6 studied by resonance Raman spectroscopy

Cytochrome P450 2D6 (CYP2D6) is one of the most important drug-metabolizing enzymes in humans. Resonance Raman data, reported for the first time for CYP2D6, show that the CYP2D6 heme is found to be in a six-coordinated low-spin state in the absence of substrates, and it is perturbed to different extents by bufuralol, dextromethorphan, and 3,4-methylenedioxymethylamphetamine (MDMA). Dextromethorphan and MDMA induce in CYP2D6 a significant amount of five-coordinated high-spin heme species and reduce the polarity of its heme-pocket, whereas bufuralol does not. Spectra of the F120A mutant CYP2D6 suggest that Phe120 is involved in substrate-binding of dextromethorphan and MDMA, being responsible for the spectral differences observed between these two compounds and bufuralol. These differences could be explained postulating a different substrate mobility for each compound in the CYP2D6 active site, consistently with the role previously suggested for Phe120 in binding dextromethorphan and MDMA.

Topological role of cytochrome P450 2D6 active site residues

Recent reports have identified Phe120, Asp301, Thr309, and Glu216 as important residues in cytochrome P450 2D6 (CYP2D6) substrate binding and catalysis. Complementary homology models have located these amino acids within the binding pocket of CYP2D6 and in the present study we have used aryldiazenes to test these models and gain further insight in the role these amino acids have in maintaining the integrity of the active site cavity. When Phe120 was replaced to alanine, there was a significant increase in probe migration to pyrrole nitrogens C and D, in agreement with homology models which have located the phenyl side-chain of Phe120 above these two pyrrole rings. No changes in topology were observed with the D301Q mutant, supporting claims that in this mutant the electrostatic interactions with the B/C-loop are largely maintained and the loop retains its native orientation. The T309V mutation resulted in significant topological alteration suggesting that, in addition to its potential role in dioxygen activation, Thr309 plays an important structural role within the active site crevice. Replacement of Ile106 with Glu, engineered to cause electrostatic repulsion with Glu216, had a profound topological effect in the higher region within the active site cavity and impaired the catalytic activity towards CYP2D6 probe substrates.

Crystal Structure of Human Cytochrome P450 2D6

Cytochrome P450 2D6 is a heme-containing enzyme that is responsible for the metabolism of at least 20% of known drugs. Substrates of 2D6 typically contain a basic nitrogen and a planar aromatic ring. The crystal structure of human 2D6 has been solved and refined to 3.0A resolution. The structure shows the characteristic P450 fold as seen in other members of the family, with the lengths and orientations of the individual secondary structural elements being very similar to those seen in 2C9. There are, however, several important differences, the most notable involving the F helix, the F-G loop, the B'helix, beta sheet 4, and part of beta sheet 1, all of which are situated on the distal face of the protein. The 2D6 structure has a well defined active site cavity above the heme group, containing many important residues that have been implicated in substrate recognition and binding, including Asp-301, Glu-216, Phe-483, and Phe-120. The crystal structure helps to explain how Asp-301, Glu-216, and Phe-483 can act as substrate binding residues and suggests that the role of Phe-120 is to control the orientation of the aromatic ring found in most substrates with respect to the heme. The structure has been compared with published homology models and has been used to explain much of the reported site-directed mutagenesis data and help understand the metabolism of several compounds.

Binding of 7-methoxy-4-(aminomethyl)-coumarin to wild-type and W128F mutant cytochrome P450 2D6 studied by time-resolved fluorescence spectroscopy

Enzyme structure and dynamics may play a main role in substrate binding and the subsequent steps in the CYP (cytochrome P450) catalytic cycle. In the present study, changes in the structure of human CYP2D6 upon binding of the substrate are studied using steady-state and time-resolved fluorescence methods, focusing not only on the emission of the tryptophan residues, but also on emission of the substrate. As a substrate, MAMC [7-methoxy-4-(aminomethyl)-coumarin] was selected, a compound exhibiting native fluorescence. As well as the wild-type, the W128F (Trp128-->Phe) mutant of CYP2D6 was studied. After binding, a variety of energy transfer possibilities exist, and molecular dynamics simulations were performed to calculate distances and relative orientations of donors and acceptors. Energy transfer from Trp128 to haem appeared to be important; its emission was related to the shortest of the three average tryptophan fluorescence lifetimes observed for CYP2D6. MAMC to haem energy transfer was very efficient as well: when bound in the active site, the emission of MAMC was fully quenched. Steady-state anisotropy revealed that besides the MAMC in the active site, another 2.4% of MAMC was bound outside of the active site to wild-type CYP2D6. The tryptophan residues in CYP2D6 appeared to be less accessible for the external quenchers iodide and acrylamide in presence of MAMC, indicating a tightening of the enzyme structure upon substrate binding. However, the changes in the overall enzyme structure were not very large, since the emission characteristics of the enzyme were not very different in the presence of MAMC.

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.

Role of the conserved threonine 309 in mechanism of oxidation by cytochrome P450 2D6

Based on sequence alignments and homology modeling, threonine 309 in cytochrome P450 2D6 (CYP2D6) is proposed to be the conserved I-helix threonine, which is supposed to be involved in dioxygen activation by CYPs. The T309V mutant of CYP2D6 displayed a strong shift from O-dealkylation to N-dealkylation reactions in oxidation of dextromethorphan and 3,4-methylenedioxymethylamphetamine. This may be explained by an elevated ratio of hydroperoxo-iron to oxenoid-iron of the oxygenating species. In consistence, using cumene hydroperoxide, which directly forms the oxenoid-iron, the T309V mutant again selectively catalyzed the O-dealkylation reactions. The changed ratio of oxygenating species can also explain the decreased activity and changed regioselectivity that were observed in 7-methoxy-4-(aminomethyl)-coumarin and bufuralol oxidation, respectively, by the T309V mutant. Interestingly, the T309V mutant always showed a significantly increased, up to 75-fold, higher activity compared to that of the wild-type when using cumene hydroperoxide. These results indicate that T309 in CYP2D6 is involved in maintaining the balance of multiple oxygenating species and thus influences substrate and regioselectivity.

Why Is Quinidine an Inhibitor of Cytochrome P450 2D6?

We have previously shown that Phe(120), Glu(216), and Asp(301) in the active site of cytochrome P450 2D6 (CYP2D6) play a key role in substrate recognition by this important drug-metabolizing enzyme (Paine, M. J., McLaughlin, L. A., Flanagan, J. U., Kemp, C. A., Sutcliffe, M. J., Roberts, G. C., and Wolf, C. R. (2003) J. Biol. Chem. 278, 4021-4027 and Flanagan, J. U., Maréchal, J.-D., Ward, R., Kemp, C. A., McLaughlin, L. A., Sutcliffe, M. J., Roberts, G. C., Paine, M. J., and Wolf, C. R. (2004) Biochem. J. 380, 353-360). We have now examined the effect of mutations of these residues on interactions of the enzyme with the prototypical CYP2D6 inhibitor, quinidine. Abolition of the negative charge at either or both residues 216 and 301 decreased quinidine inhibition of bufuralol 1'-hydroxylation and dextromethorphan O-demethylation by at least 100-fold. The apparent dissociation constants (K(d)) for quinidine binding to the wild-type enzyme and the E216D and D301E mutants were 0.25-0.50 microm. The amide substitution of Glu(216) or Asp(301) resulted in 30-64-fold increases in the K(d) for quinidine. The double mutant E216Q/D301Q showed the largest decrease in quinidine affinity, with a K(d) of 65 microm. Alanine substitution of Phe(120), Phe(481),or Phe(483) had only a minor effect on the inhibition of bufuralol 1'-hydroxylation and dextromethorphan O-demethylation and on binding. In contrast to the wild-type enzyme, a number of the mutants studied were found to be able to metabolize quinidine. E216F produced O-demethylated quinidine, and F120A and E216Q/D301Q produced both O-demethylated quinidine and 3-hydroxyquinidine metabolites. Homology modeling and molecular docking were used to predict the modes of quinidine binding to the wild-type and mutant enzymes; these were able to rationalize the experimental observations.

The role of phenylalanine 483 in cytochrome P450 2D6 is strongly substrate dependent

The polymorphic cytochrome P450 2D6 (CYP2D6) is involved in the metabolism of 30% of the drugs currently prescribed, and is thus clinically relevant. Typical CYP2D6 substrates generally contain a basic nitrogen atom and an aromatic moiety adjacent to the site of metabolism. Recently, we demonstrated the importance of active site residue F120 in substrate binding and catalysis in CYP2D6. On the basis of protein homology models, it is claimed that another active site phenylalanine, F483, may also play an important role in the interaction with the aromatic moiety of CYP2D6 substrates. Experimental data to support this hypothesis, however, is not yet available. In fact, in the only study performed, mutation of F483 to isoleucine or tryptophan did not affect the 1'-hydroxylation of bufuralol at all [Smith G, Modi S, Pillai I, Lian LY, Sutcliffe MJ, Pritchard MP, et al., Determinants of the substrate specificity of human cytochrome P-450 CYP2D6: design and construction of a mutant with testosterone hydroxylase activity. Biochem J 1998;331:783-92]. In the present study, the role of F483 in ligand binding and metabolism by CYP2D6 was examined experimentally using site-directed mutagenesis. Replacement of F483 by alanine resulted in a 30-fold lower V(max) for bufuralol 1'-hydroxylation, while the K(m) was hardly affected. The V(max) for 3,4-methylenedioxy-methylamphetamine O-demethylenation on the other hand decreased only two-fold, whereas the effect on the K(m) was much larger. For dextromethorphan, in addition to dextrorphan (O-demethylation) and 3-methoxymorphinan (N-demethylation), two other metabolites were formed that could not be detected for the wild-type. The substrate 7-methoxy-4-(aminomethyl)-coumarin was not metabolised at all by CYP2D6[F483A], a phenomenon that was reported also for CYP2D6[F120A]. The presented data show that next to F120, residue F483 plays a very important role in the metabolism of typical CYP2D6 substrates. The influence of F483 on metabolism was found to be strongly substrate-dependent.