Characterization of variant alleles of cytochrome CYP2D6 in a Spanish population

The CYP2D6 gene codes for a P450 monooxygenase which is involved in the biotransformation of a large number of commonly prescribed drugs. Adverse drug effects and therapeutic failure can be related to abnormal CYP2D6 activity. We investigated the allele and genotype frequencies of cytochrome P4502D6 in a Spanish population to predict the prevalence of ultra-rapid and poor metabolizer phenotypes in our population and to design a feasible CYP2D6 genotyping protocol. The study included 105 healthy unrelated Spanish Caucasian volunteers. CYP2D6 genotyping was performed by a combination of long-PCR, direct sequencing and allele-specific real-time PCR. The frequency of the wild-type CYP2D6*1 allele was 31%. The alleles coding for slightly (CYP2D6*2) or moderately (*9 and *10) reduced activity showed frequencies of 40.47, 2.38 and 1.90%, respectively. Frequencies of defective alleles *3, *4, *5 and *6 were 0.95, 13.8, 3.33 and 0.95%, respectively. The defective CYP2D6 alleles *7, *8, *12, *14, *15 and *21 were not found. Duplicated CYP2D6 alleles were detected at a frequency of 4.27%. Our protocol allows the identification of the four inactive CYP2D6 alleles (*3, *4, *5 and *6) and the detection of alleles with CYP2D6 *1, CYP2D6 *2 and CYP2D6*4 gene duplications. Testing for this reduced CYP2D6 allele set would facilitate its use in clinical practice by assisting in the development of individualized pharmacotherapy.     

CYP3A4 and CYP2D6 inhibitory activities of Indonesian medicinal plants

Thirty samples of Indonesian medicinal plants were analyzed for their capacity to inhibit in vitro metabolism by human cytochrome P450 3A4 (CYP3A4) and CYP2D6 with a radiometric assay. The MeOH-soluble fractions of 25 samples, prepared from water extracts, demonstrated inhibitory activity more than 50% on the metabolism mediated by CYP3A4, and 21 samples on the metabolism mediated by CYP2D6. Among the MeOH-soluble fractions, Piper nigrum leaf showed the highest inhibitory activity against CYP3A4 (91.7%), and Punica granatum against CYP2D6 (98.1%). The water extracts of which MeOH-soluble fraction showed inhibitory activity more than 70% were fractionated with EtOAc. From the EtOAc-soluble fractions, Curcuma heyneana (67.0%), Pi. cubeba (75.0%), Pi. nigrum fruit (84.0%), Pi. nigrum leaf (85.8%), and Zingiber aromaticum (75.3%) demonstrated inhibitory activity more than 50% on the metabolism mediated by CYP3A4, but only Pi. nigrum fruit (72.8%) and Pi. nigrum leaf (69.1%) showed strong inhibitory activity against CYP2D6. For samples that showed more than 70% inhibition, their IC(50) values were determined. The most potent inhibitory activity against CYP3A4 (IC(50) value of 25 microg/ml) was found for the extract of Pi. nigrum leaf, while that of Catharanthus roseus showed the most potent inhibitory effect against CYP2D6 (IC(50) value of 11 microg/ml). These results should indicate once more the possibility of potential medicinal plant-drug interactions.     

细胞色素P450 2D6酶缺陷等位基因的分析

细胞色素P450 2D6(CYP2D6)第1 795位胸腺嘧啶核苷缺失造成CYP2D6酶活性缺陷,该等位基因被称为CYP2D6T.对该等位基因的测定有助于准确预测CYP2D6表现型.利用等位基因特异扩增法的基本原理,建立了测定CYP2D6T的方法.经396例测定,证明比利用PCR扩增后再酶切的方法更为快捷、更少污染,为该项测定应用于临床奠定基础.     

Metabolism of vitamin D by human microsomal CYP2R1

The activation of vitamin D requires 25-hydroxylation in the liver and 1alpha-hydroxylation in the kidney. However, it remains unclear which enzyme is relevant to vitamin D 25-hydroxylation. Recently, human CYP2R1 has been reported to be a potential candidate for a hepatic vitamin D 25-hydroxylase. Thus, vitamin D metabolism by CYP2R1 was compared with human mitochondrial CYP27A1, which used to be considered a physiologically important vitamin D(3) 25-hydroxylase. A clear difference was observed between CYP2R1 and CYP27A1 in the metabolism of vitamin D(2). CYP2R1 hydroxylated vitamin D(2) at the C-25 position while CYP27A1 hydroxylated it at positions C-24 and C-27. The K(m) and k(cat) values for the CYP2R1-dependent 25-hydroxylation activity toward vitamin D(3) were 0.45microM and 0.97min(-1), respectively. The k(cat)/K(m) value of CYP2R1 was 26-fold higher than that of CYP27A1. These results strongly suggest that CYP2R1 plays a physiologically important role in the vitamin D 25-hydroxylation in humans.     

Interactions studies of CYP2D6 with quercetin and hyperoside by spectral analysis and molecular dynamics simulations

Some ingredients from herbal medicine can significantly affect the activity of CYP2D6, thus leading to serious interactions between herbs and drugs. Quercetin and hyperoside are active ingredients widely found in vegetables, fruits, and herbal medicines. Quercetin and hyperoside have many biological activities. In this work, the characteristic bindings of CYP2D6 with quercetin/hyperoside are revealed by multi-spectroscopy analysis, molecular docking, and molecular dynamics simulations. The fluorescence of CYP2D6 is statically quenched by quercetin and hyperoside. The binding constant (K_a) values of CYP2D6–quercetin/hyperoside range from 10⁴ L mol⁻¹, which indicates that these two flavonoids bind moderately to CYP2D6. Meanwhile, quercetin has a stronger quenching ability to CYP2D6 than that of hyperoside. The secondary structure of CYP2D6 is obviously changed by binding with quercetin/hyperoside. The docking results reveal that the quercetin/hyperoside enters the active site of CYP2D6 near heme and binds to CYP2D6 by hydrogen bonds and van der Waals forces. The molecular dynamics simulation results indicate that the binding of quercetin/hyperoside can stabilize the two complexes, enhance the flexibility of CYP2D6 backbone atoms, and make a more unfolded and looser structure of CYP2D6.     

Pharmacophore-directed Homology Modeling and Molecular Dynamics Simulation of G Protein-coupled Receptor： Study of Possible Binding Modes of 5-HT2c Receptor Agonists

A new pharmacophore-based modeling procedure, including homology modeling, pharmacophore study, flexible molecular docking, and long-time molecular dynamics (MD) simulations, was employed to construct the structure of the human 5-HT_(2C) receptor and determine the characteristics of binding modes of 5-HT_(2C) receptor agonists. An agonist-receptor complex has been constructed based on homology modeling and a pharmacophore hypothesis model based on some high active compounds. Then MD simulations of the ligand-receptor complex in an explicit membrane environment were carried out. The conformation of the 5- HT_(2C) receptor during MD simulation was explored, and the stable binding modes of the studied agonist were determined. Flexible molecular docking of several structurally diverse agonists of the human 5-HT_(2C) receptor was carried out, and the general binding modes of these agonists were investigated. According to the models presented in this work and the results of Flexi-Dock, the involvement of the amino acid residues Asp134, Ser138, Ash210, Asn331, Tyr358, Ile131, Ser132, Val135, Thr139, Ile189, Val202, Val208, Leu209, Phe214, Val215, Gly218, Ser219, Phe223, Trp324, Phe327, and Phe328 in agonist recognition was studied. The obtained binding modes of the human 5-HT_(2C) receptor agonists have good agreement with the site-directed mutagenesis data and other studies.     

Probing ligand binding modes of human cytochrome P450 2J2 by homology modeling, molecular dynamics simulation, and flexible molecular docking

Cytochrome P450 (P450) 2J2 catalyzes epoxidation of arachidonic acid to eicosatrienoic acids, which are related to a variety of diseases such as coronary artery disease, hypertension, and carcinogenesis. Recent experimental data also suggest that P450 2J2 could be a novel biomarker and a potential target for cancer therapy. However, the active site topology and substrate specificity of this enzyme remain unclear. In this study, a three-dimensional model of human P450 2J2 was first constructed on the basis of the crystal structure of human P450 2C9 in complex with a substrate using homology modeling method, and refined by molecular dynamics simulation. Flexible docking approaches were then employed to dock four ligands into the active site of P450 2J2 in order to probe the ligand-binding modes. By analyzing the results, active site architecture and certain key residues responsible for substrate specificity were identified on the enzyme, which might be very helpful for understanding the enzyme's biological role and providing insights for designing novel inhibitors of P450 2J2.     

Determination of risperidone and 9-hydroxyrisperidone in human plasma by liquid chromatography: application to the evaluation of CYP2D6 drug interactions

A high-pressure liquid chromatography with ultra-violet detection method for the simultaneous determination of risperidone and 9-hydroxyrisperidone in plasma after liquid-liquid extraction has been developed. The limit of quantitation was 5 nmol/L, and the inter-day coefficient of variation was less than 8% for both compounds. The mean recoveries of risperidone and 9-hydroxyrisperidone added to plasma were 96.8 and 99.4%, with an intra-day coefficient of variation of under 5 and 6%, respectively. Studies of analytical interference showed that the most commonly co-administered antidepressants and benzodiazepines did not interfere. The method was used for the determination of the plasma concentrations of a schizophrenic patient treated daily with an oral dose of 4.5 mg risperidone. The patient suffered severe extrapyramidal side-effects after adding risperidone to his previous medication of haloperidol and levomepromazine. The risperidone plasma concentration was well above the average (182 nmol/L), which suggests that a pharmacokinetic interaction occurred, presumably due to inhibition of the enzyme CYP2D6.     

Refinement of 3D models of horseradish peroxidase isoenzyme C: Predictions of 2D NMR assignments and substrate binding sites

In this study, two alternative three-dimensional (3D) models of horseradish peroxidase (HRP-C)—differing mainly in the structure of a long untemplated insertion—were refined, systematically assessed, and used to make predictions that can both guide and be tested by future experimental studies. A key first step in the model-building process was a procedure for multiple sequence alignment based on structurally conserved regions and key conserved residues, including those side chains providing ligands to the two Ca²⁺ binding sites. The model refinements reported here include (1) optimization of side-chain conformations; (3) addition of structural waters using a template-independent procedure; (2) structural refinement of the untemplated 34 amino acid insertion located between the F and G helices, using both energy criteria and NMR data; (4) unconstrained energy optimization of the refined models. Using these procedures, two refined structures of HRP-C were obtained, differing mainly in the conformation of this long insertion. The presence of residues in this insertion that could potentially interact with bound substrates suggests a functional role that may be related to the general ability of class III peroxidases to form stable 1:1 complexes with a variety of substrates. The structural validity of the models was systematically assessed by a variety of criteria. Most notably, the ProsaII z scores and Profiles 3D scores of the two HRP-C models indicated that they are significantly better than would be obtained by simple amino acid replacement, using any of the known structures as a template. These two 3D HRP-C models, were then used to predict candidate residues for the assignment of NOESY cross-peaks previously noted in 2D-NMR studies. Specifically, the residues known as Ile X, Phe A, Phe B, aliphatic residue Q, and Ile T. Candidate substrate binding sites were also identified and compared with experimentally based predictions. This work is timely because new X-ray structures are anticipated that will facilitate the validation of these procedures. © 1996 Wiley-Liss, Inc.     

Hybrid homology modeling and mutational analysis of cytochrome P450C24A1 (CYP24A1) of the Vitamin D pathway: insights into substrate specificity and membrane bound structure-function

Cytochrome P450C24A1 (CYP24A1), a peripheral inner mitochondrial membrane hemoprotein and candidate oncogene, regulates the side-chain metabolism and biological function of vitamin D and many of its related analog drugs. Rational mutational analysis of rat CYP24A1 based on hybrid (2C5/BM-3) homology modeling and affinity labeling studies clarified the role of key domains (N-terminus, A', A, and F-helices, beta3a strand, and beta5 hairpin) in substrate binding and catalysis. The scope of our study was limited by an inability to purify stable mutant enzyme targeting soluble domains (B', G, and I-helices) and suggested greater conformational flexibility among CYP24A1's membrane-associated domains. The most notable mutants developed by modeling were V391T and I500A, which displayed defective-binding function and profound metabolic defects for 25-hydroxylated vitamin D3 substrates similar to a non-functional F-helix mutant (F249T) that we previously reported. Val-391 (beta3a strand) and Ile-500 (beta5 hairpin) are modeled to interact with Phe-249 (F-helix) in a hydrophobic cluster that directs substrate-binding events through interactions with the vitamin D cis-triene moiety. Prior affinity labeling studies identified an amino-terminal residue (Ser-57) as a putative active-site residue that interacts with the 3beta-OH group of the vitamin D A-ring. Studies with 3-epi and 3-deoxy-1,25(OH)2D3 analogs confirmed interactions between the 3beta-OH group and Ser-57 effect substrate recognition and trafficking while establishing that the trans conformation of A-ring hydroxyl groups (1alpha and 3beta) is obligate for high-affinity binding to rat CYP24A1. Our work suggests that CYP24A1's amphipathic nature allows for monotopic membrane insertion, whereby a pw2d-like substrate access channel is formed to shuttle secosteroid substrate from the membrane to the active-site. We hypothesize that CYP24A1 has evolved a unique amino-terminal membrane-binding motif that contributes to substrate specificity and docking through coordinated interactions with the vitamin D A-ring.     

Regional and cellular expression of CYP2D6 in human brain: higher levels in alcoholics

Cytochrome P450 (CYP) 2D6 is expressed in liver, brain and other extrahepatic tissues where it metabolizes a range of centrally acting drugs and toxins. As ethanol can induce CYP2D in rat brain, we hypothesized that CYP2D6 expression is higher in brains of human alcoholics. We examined regional and cellular expression of CYP2D6 mRNA and protein by RT-PCR, Southern blotting, slot blotting, immunoblotting and immunocytochemistry. A significant correlation was found between mean mRNA and CYP2D6 protein levels across 13 brain regions. Higher expression was detected in 13 brain regions of alcoholics (n = 8) compared to nonalcoholics (n = 5) (anovap < 0.0001). In hippocampus this was localized in CA1-3 pyramidal cells and dentate gyrus granular neurons. In cerebellum this was localized in Purkinje cells and their dendrites. Both of these brain regions, and these same cell-types, are known to be susceptible to alcohol damage. For one case, a poor metabolizer (CYP2D6*4/*4), there was no detectable CYP2D6 protein, confirming the specificity of the antibody used. These data suggest that in alcoholics elevated brain CYP2D6 expression may contribute to altered sensitivity to centrally acting drugs and to the mediation of neurotoxic and behavioral effects of alcohol.     

Comparative modeling and molecular docking of orphan human CYP4V2 protein with fatty acid substrates: Insights into substrate specificity

Cytochromes P450 (CYPs) are a super family of heme-containing enzymes well-known for their monooxgenase reaction. There are 57 CYP isoenzymes found in human which exhibit specific physiological functions. Thirteen members of this super family are classified as "orphan" CYP because of their unknown enzymatic functions. CYP4V2 is found to be a potential drug target for Bietti crystalline corneoretinal dystrophy (BCD). However, three-dimensional structure, the active site topology and substrate binding modes of CYP4V2 remain unclear. In this study, the three-dimensional model of CYP4V2 was constructed using the homology modeling method. Four possible fatty acid substrates namely, caprylic, lauric, myrisitc and palmitic acids were optimized and evaluated for drug likeness using Lipinski's rule of five. Further, these substrates were docked into active sites of CYP4V2 and several key residues responsible for substrate binding were identified. These findings will be helpful for the structure-based drug design and detailed characterization of the biological roles of CYP4V2.     

小菜蛾中肠氨肽酶基因PxAPN5的克隆、原核表达及蛋白质同源建模分析

【目的】克隆和表达小菜蛾Plutella xylostella氨肽酶基因,并进行基因序列分析和同源建模分析。【方法】以小菜蛾中肠cDNA为模板克隆分析氨肽酶基因序列, 原核表达氨肽酶蛋白并进行酶活性测定, 应用配体印迹分析氨肽酶与Cry2Ab蛋白的结合, 通过蛋白质建模对突变位点进行分析。【结果】从小菜蛾中肠cDNA 扩增出氨肽酶基因, 该基因全长2 853 bp, 编码950个氨基酸, 预测蛋白分子量为107.3871 kDa, 等电点为5.24; 进化树分析显示, 克隆得到的氨肽酶基因属于APN家族5, 将其命名为PxAPN5（GenBank登录号: KM034756）。PxAPN5蛋白具有鳞翅目昆虫氨肽酶蛋白的保守性特征, 即含有N-糖基化位点、O-糖基化位点和GPI锚定位点, 具有“HEXXH”锌蛋白酶结构域和C端跨膜区域。在大肠杆菌Escherichia coli中原核表达PxAPN5, 表达产物经SDS-PAGE电泳, 在110 kDa附近出现特异性条带; 酶活性测试显示菌体破碎上清液具有氨肽酶活性, 比活力为1 047.2 U/g。配体印迹结果显示表达的PxAPN5能与Cry2Ab蛋白特异性结合。多序列比对结果表明, 与其他已报道的小菜蛾氨肽酶相比, PxAPN5氨基酸序列有3个保守性位点发生了突变,并通过蛋白质建模的方式表征突变位点。【结论】本研究成功克隆和表达了具有氨肽酶活性的小菜蛾氨肽酶, 并发现其能与Cry2Ab蛋白特异性结合; 通过蛋白质建模对氨肽酶突变位点的特征及功能进行了预测。 这些结果对小菜蛾氨肽酶的功能性研究提供了理论基础。     

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.     

Prediction of possible sites for posttranslational modifications in human gamma crystallins: effect of glycation on the structure of human gamma-B-crystallin as analyzed by molecular modeling

Crystallins are recognized as one of the long-lived proteins of lens tissue that might serve as the target for several posttranslational modifications leading to cataract development. We have studied several such sites present in the human gamma-crystallins based either on PROSITE pattern search results or earlier experimental evidences. Their probabilities were examined on the basis of the database analysis of the gamma-crystallin sequences and on their specific locations in the constructed homology models. An N-glycosylation site in human gammaD-crystallin and several phosphorylation sites in all four human gamma-crystallins were predicted by the PROSITE search. Some of these sites were found to be strongly conserved in the gamma-crystallin sequences from different sources. An extensive analysis of these sites was performed to predict their probabilities as potential sites for protein modifications. Glycation studies were performed separately by attaching sugars to the human gammaB-crystallin model, and the effect of binding was analyzed. The studies showed that the major effect of alphaD-glucose (alphaD-G) and alphaD-glucose-6-phosphate (alphaD-G6P) binding was the disruption of charges not only at the surface but also within the molecule. Only a minor alteration in the distances of sulfhydryl groups of cysteines and on their positions in the three-dimensional models were observed, leading us to assume that glycation alone is not responsible for intra- and intermolecular disulfide bond formation.     

Construction of molecular assemblies via docking: modeling of tetramers with D2 symmetry

Comparative modeling methods are commonly used to construct models of homologous proteins or oligomers. However, comparative modeling may be inapplicable when the number of subunits in a modeled oligomer is different than in the modeling template. Thus, a dimer cannot be a template for a tetramer because a new monomer-monomer interface must be predicted. We present in this study a new prediction approach, which combines protein-protein docking with either of two tetramer-forming algorithms designed to predict the structures of tetramers with D2 symmetry. Both algorithms impose symmetry constraints. However, one of them requires identification of two of the C2 dimers within the tetramer in the docking step, whereas the other, less demanding algorithm, requires identification of only one such dimer. Starting from the structure of one subunit, the procedures successfully reconstructed 16 known D2 tetramers, which crystallize with either a monomer, a dimer or a tetramer in the asymmetric unit. In some cases we obtained clusters of native-like tetramers that differ in the relative rotation of the two identical dimers within the tetramer. The predicted structural pliability for concanavalin-A, phosphofructokinase, and fructose-1,6-bisphosphatase agrees semiquantitatively with the observed differences between the several experimental structures of these tetramers. Hence, our procedure identifies a structural soft-mode that allows regulation via relative rigid-body movements of the dimers within these tetramers. The algorithm also predicted three nearly correct tetramers from model structures of single subunits, which were constructed by comparative modeling from subunits of homologous tetrameric, dimeric, or hexameric systems.     

Locomotor activity is regulated by D2-like receptors in Drosophila: an anatomic and functional analysis

In mammals, dopamine 2-like receptors are expressed in distinct pathways within the central nervous system, as well as in peripheral tissues. Selected neuronal D2-like receptors play a critical role in modulating locomotor activity and, as such, represent an important therapeutic target (e.g. in Parkinson's disease). Previous studies have established that proteins required for dopamine (DA) neurotransmission are highly conserved between mammals and the fruit fly Drosophila melanogaster. These include a fly dopamine 2-like receptor (DD2R; Hearn et al. PNAS 2002 99(22):14554) that has structural and pharmacologic similarity to the human D2-like (D2R). In the current study, we define the spatial expression pattern of DD2R, and functionally characterize flies with reduced DD2 receptor levels. We show that DD2R is expressed in the larval and adult nervous systems, in cell groups that include the Ap-let cohort of peptidergic neurons, as well as in peripheral tissues including the gut and Malpighian tubules. To examine DD2R function in vivo, we generated RNA-interference (RNAi) flies with reduced DD2R expression. Behavioral analysis revealed that these flies show significantly decreased locomotor activity, similar to the phenotype observed in mammals with reduced D2R expression. The fly RNAi phenotype can be rescued by administration of the DD2R synthetic agonist bromocriptine, indicating specificity for the RNAi effect. These results suggest Drosophila as a useful system for future studies aimed at identifying modifiers of dopaminergic signaling/locomotor function.     

Generation of a homology model for the human cytochrome P450, CYP24A1, and the testing of putative substrate binding residues by site-directed mutagenesis and enzyme activity studies

A systematic analysis of conserved H-bonding patterns and tertiary structural motifs from 13 crystal structures was used to create a homology model for the human multicatalytic cytochrome P450, CYP24A1, involved in catabolism of 1alpha,25-dihydroxyvitamin D3. The substrate was docked in the active site and used to identify potential substrate contact residues in the B' helix, B'/C loop, F-helix and the beta-5 hairpin. Seven CYP24A1 mutants were created and studied by mammalian cell transfection and CYP24A1 activity assay. Mutants showed reduced metabolic rates and altered metabolite patterns compared to wild-type. We conclude that: Ile-131 positions substrate via A-ring and cis-triene contacts; Trp-134 and Gly-499 are determinants of substrate access; Leu-148 contacts the substrate side-chain; Met-246 is important in mediating regioselectivity. Our findings validate the new model of CYP24A1, which can now be used to predict structural modifications for rational vitamin D drug design.     

Influence of pKa on the biotransformation of indene H1-antihistamines by CYP2D6

Structure-activity relationship studies were conducted to reduce CYP2D6-mediated metabolism in a series of indene H₁-antihistamines. Reductions in pK_a via incorporation of a β-fluoro substituent or a heteroaryl moiety were shown to reduce contributions to metabolism through this pathway. Several compounds, including 8l, 8o, and 12f were identified with promising primary in vitro profiles and reduced biotransformation via CYP2D6.