Mechanism of the decrease in catalytic activity of human cytochrome P450 2C9 polymorphic variants investigated by computational analysis

Cytochrome P450 (CYP) is deeply involved in the metabolism of chemicals including pharmaceuticals. Therefore, polymorphisms of this enzyme have been widely studied to avoid unfavorable side effects of drugs in chemotherapy. In this work, we performed computational analysis of the mechanism of the decrease in enzymatic activity for three typical polymorphisms in CYP 2C9 species: *2, *3, and *5. Based on the equilibrated structure obtained by molecular dynamics simulation, the volume of the binding pocket and the fluctuation of amino residues responsible for substrate holding were compared between the wild type and the three variants. Further docking simulation was carried out to evaluate the appropriateness of the binding pocket to accommodate substrate chemicals. Every polymorphic variant was suggested to be inferior to the wild type in enzymatic ability from the structural viewpoint. F‐G helices were obviously displaced outward in CYP2C9*2. Expansion of the binding pocket, especially the space near F′ helix, was remarkable in CYP2C9*3. Disappearance of the hydrogen bond between K helix and β4 loop was observed in CYP2C9*5. The reduction of catalytic activity of those variants can be explained from the deformation of the binding pocket and the consequent change in binding mode of substrate chemicals. The computational approach is effective for predicting the enzymatic activity of polymorphic variants of CYP. This prediction will be helpful for advanced drug design because calculations forecast unexpected change in drug efficacy for individuals. © 2010 Wiley Periodicals, Inc. J Comput Chem, 2010     

Application of dried blood spot sampling combined with LC-MS/MS for genotyping and phenotyping of CYP450 enzymes in healthy volunteers

An early clinical development study (phase I) was conducted to determine the usefulness of dried blood spot (DBS) sampling as an alternative to venous sampling for phenotyping and genotyping of CYP450 enzymes in healthy volunteers. Midazolam (MDZ) was used as a substrate for phenotyping CYP3A4 activity; the concentrations of MDZ and its main metabolite 1'-hydroxymidazolam (1-OH MDZ) were compared between the DBS method from finger punctures, plasma and whole blood (WB), drawn by venipuncture, whereby several methodological parameters were studied (i.e. punch width, amount of dots analyzed and storage time stability). Genotyping between DBS and venous WB samples was compared for CYP2D6 (*3, *4, *6), CYP2C19 (*2, *3), CYP3A4 (*1B) and CYP3A5 (*3C). In addition, the subject's and phlebotomist's satisfaction with venous blood sampling compared with the DBS method was evaluated using a standardized questionnaire. An LC-MS/MS method for the quantification of the MDZ and 1-OH MDZ concentrations in DBS samples was developed and validated in the range of 0.100-100 ng/mL. No compromises were made for the limits of quantification of the DBS-LC-MS/MS method vs the authentic plasma and WB methods.     

一种细胞色素P4502C8基因多态-CYP2C8.3的结构和功能

细胞色素P450超级家族在代谢众多的外源性化学物质方面发挥重要的作用.细胞色素P4502C8是人体肝脏中主要负责代谢抗癌药物紫杉醇的酶,它至少负责代谢5%的临床药物.细胞色素P450 2C8的基因多态性与用药个体化有着密切的关系.CYP2C8.3是常见的P450 2C8的基因多态之一,其发生了双点突变,分别是R139K...     

Structural and functional insights into CYP2C8.3: A genetic polymorph of cytochrome P450 2C8

The cytochrome P450 (CYP) superfamily plays a key role in the oxidative metabolism of a wide range of exogenous chemicals. CYP2C8 is the principal enzyme responsible for the metabolism of the anti-cancer drug paclitaxel in the human liver, and carries out the oxidative metabolism of at least 5% of clinical drugs. Polymorphisms in CYP2C8 have been closely implicated in individualized medication. CYP2C8.3, a common polymorph of CYP2C8 with dual amino acid substitutions R139K and K399R, is found primarily in Caucasians. In this study, CYP2C8.3 and its wild type (WT) CYP2C8 were expressed in E. coli, and their purified proteins were characterized by UV-visible spectroscopy, mass spectrometry, and circular dichroism. Their thermal stability, substrate binding ability, and metabolic activity against paclitaxel were investigated. The electron transfer kinetics during paclitaxel metabolism by WT CYP2C8 or CYP2C8.3 was studied by stopped-flow kinetics. The results revealed that mutations in CYP2C8.3 did not greatly influence the heme active site or protein thermal stability and paclitaxel binding ability, but the metabolic activity against paclitaxel was significantly depressed to just 11% of that of WT CYP2C8. Electron transfer from CYP reductase to CYP2C8.3 was found to be significantly slower than that to WT CYP2C8 during catalysis, and this might be the main reason for the depressed metabolic activity. Since the polymorph CYP2C8.3 is defective in catalyzing substrates of CYP2C8 in vitro, it might be expected to have important clinical and pathophysiological consequences in homozygous individuals, and this study provides valuable information in this aspect.     

Insight into G[bond]T mismatch recognition using molecular dynamics with time-averaged restraints derived from NMR spectroscopy

Molecular dynamics (MD) simulations were conducted for a G[bond]T mismatch-containing DNA decamer, d(CCATGCGTGG)(2), and its Watson-Crick parent sequence, d(CCACGCGTGG)(2). Dynamics in unrestrained MD trajectories were in poor agreement with prior (13)C NMR studies. However, the accuracy of the trajectories was improved by the use of time-averaged interatomic distance restraints derived from (1)H NMR. Postprocess smoothing of the trajectories further improved accuracy. Comparison of restrained and smoothed trajectories of the two DNA molecules revealed distinct differences in dynamics. The major groove width of the mismatched oligomer was more variable over the course of the simulation compared to its parent sequence. Greater variability in helical parameters stretch and opening for the mismatches indicated less kinetically stable base pairing. Interbase helical parameters rise, roll, and tilt were also more variable in certain base steps involving mismatched bases. These dynamic differences between normal and G[bond]T mismatched DNA reflect differences in local flexibility that may play a role in mismatch recognition by the MutS. A potential alternate G[bond]T mismatch binding mode for MutS is also proposed.     

Phenotyping of CYP450 in human liver microsomes using the cocktail approach

The cocktail approach is an advantageous strategy used to monitor the activities of several cytochromes P450 (CYPs) in a single test to increase the throughput of in vitro phenotyping studies. In this study, a cocktail mixture was developed with eight CYP-specific probe substrates to simultaneously evaluate the activity of the most important CYPs, namely, CYP1A2, CYP2A6, CYP2B6, CYP2C9, CYP2C19, CYP2D6, CYP2E1, and the CYP3A subfamily. After cocktail incubation in the presence of human liver microsomes (HLMs), the eight selected substrates and their specific metabolites were analyzed by ultra-high-pressure liquid chromatography and electrospray ionization quadrupole time-of-flight mass spectrometry. Qualitative and quantitative data were simultaneously acquired to produce an overview of the extended phase I biotransformation routes for each probe substrate in the HLMs and to generate phenotypic profiles of various HLMs. A comparison of the cocktail strategy with an individual substrate assay for each CYP produced similar results. Moreover, the cocktail was tested on HLMs with different allelic variants and/or in the presence of selective inhibitors. The results were in agreement with the genetic polymorphisms of the CYPs and the expected effect of the alterations. All of these experiments confirmed the reliability of this cocktail assay for phenotyping of the microsomal CYPs.     

Expression and characterization of mutant forms of penicillin acylase from Alcaligenes faecalis

The sequencing of six plasmids carrying a gene of penicillin acylase from Alcaligenes faecalis VKM B1518 (AfPA) revealed the presence of random mutations in the gene; they occurred during a polymerase chain reaction. Six mutant AfPAs and a wild-type enzyme were expressed in E. coli cells. The activity assay of mutant AfPAs in E. coli cells indicated that several amino acid substitutions affect the expression level of the AfPA gene and the rate of cell growth. Four mutant AfPAs were purified; their catalytic properties and thermal stability were studied. It is shown that the amino acid substitutions under study do not affect the catalytic efficiency value. Within the experimental error, the βQ133R and βK184E (the AfPA M2 mutant) substitutions had no effect on the thermal stability of the enzyme; in the case of mutants AfPA M4 (βY90H), M5 (αD132G, βR97C), and M6 (αV5E, αN183S, and βE439G), the inactivation rate constant increased 2.4, 2.75, and 8.3 times, respectively, as compared to that of the wild-type enzyme.     

Characterization of the biotransformation pathways of clomiphene, tamoxifen and toremifene as assessed by LC-MS/(MS) following in vitro and excretion studies

The use of selective oestrogen receptor modulators has been prohibited since 2005 by the World Anti-Doping Agency regulations. As they are extensively cleared by hepatic and intestinal metabolism via oxidative and conjugating enzymes, a complete investigation of their biotransformation pathways and kinetics of excretion is essential for the anti-doping laboratories to select the right marker(s) of misuse. This work was designed to characterize the chemical reactions and the metabolizing enzymes involved in the metabolic routes of clomiphene, tamoxifen and toremifene. To determine the biotransformation pathways of the substrates under investigation, urine samples were collected from six subjects (three females and three males) after oral administration of 50 mg of clomiphene citrate or 40 mg of tamoxifen or 60 mg of toremifene, whereas the metabolizing enzymes were characterized in vitro, using expressed cytochrome P450s and uridine diphosphoglucuronosyltransferases. The separation, identification and determination of the compounds formed in the in vivo and in vitro experiments were carried out by liquid chromatography coupled with mass spectrometry techniques using different acquisition modes. Clomiphene, tamoxifen and toremifene were biotransformed to 22, 23 and 18 metabolites respectively, these phase I reactions being catalyzed mainly by CYP3A4 and CYP2D6 isoforms and, to a lesser degree, by CYP3A5, CYP2B6, CYP2C9, CYP2C19 isoforms. The phase I metabolic reactions include hydroxylation in different positions, N-oxidation, dehalogenation, carboxylation, hydrogenation, methoxylation, N-dealkylation and combinations of them. In turn, most of the phase I metabolites underwent conjugation reaction to form the corresponding glucuro-conjugated mainly by UGT1A1, UGT1A3, UGT1A4, UGT2B7, UGT2B15 and UGT2B17 isoenzymes.     

Rapid testing of clopidogrel resistance by genotyping of CYP2C19 and CYP2C9 polymorphisms using denaturing on-chip capillary electrophoresis

Antiplatelet therapy is a cornerstone of cardiovascular treatment in patients with coronary artery disease and after myocardial infarction. Clopidogrel has become a popular antiplatelet agent due to its fast action and low frequency of adverse effects. Kinetics of clopidogrel metabolism is driven by enzymatic activity of the Cytochrome P450 system. Genotyping of CYP2C19 and CYP2C9 polymorphisms allows to identify slow metabolizers showing resistance to clopidogrel therapy. Today, a number of PCR-based techniques for single nucleotide polymorphism genotyping directed at clopidogrel resistance polymorphisms are in use. Here, we describe a new alternative genotyping approach combining the separation power of denaturing capillary electrophoresis with the analysis speed and ease of use of Bioanalyzer chipCE platform. Using an upgraded heater control, we present an optimization for allele separation of CYP2C19 I331V, CYP2C9 R144C, and CYP2C9 I359L polymorphisms employing run temperatures of up to 55°C. We demonstrate rapid and accessible approach to reproducible clopidogrel resistance with feasibility and low cost.     

一种由全氟碘代烷合成全氟羧酸的新方法

本文报道全氟碘代烷在光氧化条件下生成相应酸羧的反应, 并对其机理进行了初步的探讨。以甲醇为溶剂, RFCF2I在Et3N存在下能顺利地进行光氧化反应, 反应完毕, 除去溶剂及Et3N即可分得全氟羧酸(RFCO2H), 反应操作简便, 收率较好, 无RFCF2H生成。     

Effective virtual screening protocol for CYP2C9 ligands using a screening site constructed from flurbiprofen and S-warfarin pockets

An effective virtual screening protocol was developed against an extended active site of CYP2C9, which was derived from X-ray structures complexed with flubiprofen and S-warfarin. Virtual screening has been effectively supported by our structure-based pharmacophore model. Importance of hot residues identified by mutation data and structural analysis was first estimated in an enrichment study. Key role of Arg108 and Phe114 in ligand binding was also underlined. Our screening protocol successfully identified 76% of known CYP2C9 ligands in the top 1% of the ranked database resulting 76-fold enrichment relative to random situation. Relevance of the protocol was further confirmed in selectivity studies, when 89% of CYP2C9 ligands were retrieved from a mixture of CYP2C9 and CYP2C8 ligands, while only 22% of CYP2C8 ligands were found applying the structure-based pharmacophore constraints. Moderate discrimination of CYP2C9 ligands from CYP2C18 and CYP2C19 ligands could also be achieved extending the application domain of our virtual screening protocol for the entire CYP2C family. Our findings further demonstrate the existence of an active site comprising of at least two binding pockets and strengthens the need of involvement of protein flexibility in virtual screening.     

Flexibility of human cytochromes P450: molecular dynamics reveals differences between CYPs 3A4, 2C9, and 2A6, which correlate with their substrate preferences

Molecular dynamics (MD) simulations at normal and high temperature were used to study the flexibility and malleability of three microsomal cytochromes P450 (CYPs): CYP3A4, CYP2C9, and CYP2A6. Comparison of B-factors (describing the atomic fluctuations) between X-ray and MD data shows that the X-ray B-factors are significantly lower in the regions where the crystal contacts occur than for other regions. Consequently, the conclusions about CYP flexibility based solely on the X-ray data might be misleading. Comparison of flexibility patterns of the three CYPs enabled common features and variations in flexibility and malleability of the studied CYPs to be identified. The previously described pattern of flexibility in topological elements of microsomal CYPs (a rigid heme binding core, a malleable distal side and intermediately flexible proximal side) was confirmed. These topological features provide an important combination of high stereo- and regio-specificity (mediated by the relative rigidity in the neighborhood of the heme), together with high substrate promiscuity due to the more flexible active site and the malleability of the distal side. The data acquired here show that the malleability of the three studied CYPs correlates with their substrate specificity: CYP2A6 has a narrow substrate range and is the most rigid, CYP3A4 is the most promiscuous CYP known and is the most malleable, and CYP2C9 is intermediate in terms of both its substrate specificity and malleability. Thus, the malleability of CYPs is probably a major determinant of their substrate specificity.     

Molecular Docking and Molecular Dynamics Study on the Effect of ERCC1 Deleterious Polymorphisms in ERCC1-XPF Heterodimer

Excision repair cross complementation group 1 (ERCC1) is an important protein in the nucleotide excision repair (NER) pathway, which is responsible for removing DNA adducts induced by platinum based compounds. The heterodimer ERCC1-XPF is one of two endonucleases required for NER. Genetic variations or polymorphisms in ERCC1 gene alter DNA repair capacity. Reduced DNA repair (NER) capacity may result in tumors and enhances cisplatin chemotherapy in cancer patients, which functions by causing DNA damage. Therefore, ERCC1 variants have the potential to be used as a strong candidate biomarker in cancer treatments. In this study we identified five variants V116M, R156Q, A199T, S267P, and R322C of ERCC1 gene as highly deleterious. Further structural and functional analysis has been conducted for ERCC1 protein in the presence of three variants V116M, R156Q, and A199T. Occurrence of theses variations adversely affected the regular interaction between ERCC1 and XPF protein. Analysis of 20 ns molecular dynamics simulation trajectories reveals that the predicted deleterious variants altered the ERCC1-XPF complex stability, flexibility, and surface area. Notably, the number of hydrogen bonds in ERCC1-XPF mutant complexes decreased in the molecular dynamic simulation periods. Overall, this study explores the link between the ERCC1 deleterious variants and cisplatin chemotherapy for various cancers with the help of molecular docking and molecular dynamic approaches.     

Synthesis and Crystal Structure of Ethyl2—（2—Amino—3—ethoxycarbonyl—4H—benzopyran—4—yl）cyanoacetate

1 INTRODUCTION Inorganic solid supports as catalysts resulting in higher selectivity, milder conditions and easier work-up have been reported as useful catalysts for many reactions [1~3]. Recently, we have reported the Knoevenagel condensation catalyzed by KF-Al2O3[4]. In this paper, we discussed the crystal structure of the title compound synthesized by the reaction of salicylaldehyde and ethyl cyanoacetate in DMF using the catalyst KF-Al2O3 at room temperature KF-Al2O3. In or…     

Molecular dynamics simulation of folding of a short helical peptide with many charged residues

A molecular dynamics simulation of the folding of conantokin-T (con-T), a short helical peptide with 5 helical turns of 21 amino acids with 10 charged residues, was carried out to examine folding pathways for this peptide and to predict the folding rate. In the 18 trajectories run at 300 K, 16 trajectories folded, with an averaged folding time of approximately 50 ns. Two trajectories did not fold in up to 200 ns simulation. The folded structure in folded trajectories is in good agreement with experimental structure. An analysis of the trajectories showed that, at the beginning of a few nanoseconds, helix formation started from residues 5-9 with assistance of a hydrophobic clustering involving Tyr5, Met8, and Leu9. The peptide formed a U-shape mainly due to charge-charge interactions between charged residues at the N- and C-terminus segments. In the next approximately 10 ns, several nonnative charge-charge interactions were broken and nonnative Gla10-Lys18 (this denotes a salt bridge between Gal10 and Lys18) and/or Gla10-Lys19 interactions appeared more frequently in this folding step and the peptide became a fishhook J-shape. From this structure, the peptide folded to the folded state in 7 of all 16 folded trajectories in approximately 15 ns. Alternatively, in approximately 30 ns, the con-T went to a conformation in an L-shape with 4 helical turns and a kink at the Arg13 and Gla14 segment in the other 9 trajectories. Con-T in the L-shape then required another approximately 15 ns to fold into the folded state. In addition, in overall folding times, the former 7 trajectories folded faster with the total folding times all shorter than 45 ns, while the latter 9 trajectories folded at a time longer than 45 ns, resulting in an average folding time of approximately 50 ns. Two major folding intermediates found in 2 nonfolded trajectories are stabilized by charge clusters of 5 and 6 charged residues, respectively. With inclusion of friction and solvent-solvent interactions, which were ignored in the present GB/SA solvation model, the folding time obtained above should be multiplied by a factor of 1.25-1.7 according to a previous, similar simulation study. This results in a folding time of 65-105 ns, slightly shorter than the folding time of 127 ns for an alanine-based peptide of the same length. This suggests that the energy barrier of folding for this type of peptides with many charged residues is slightly lower than alanine-based helical peptides by less than 1 kcal/mol.     

Rational reprogramming of the R2 subunit of Escherichia coli ribonucleotide reductase into a self-hydroxylating monooxygenase.

The outcome of O2 activation at the diiron(II) cluster in the R2 subunit of Escherichia coli (class I) ribonucleotide reductase has been rationally altered from the normal tyrosyl radical (Y122*) production to self-hydroxylation of a phenylalanine side-chain by two amino acid substitutions that leave intact the (histidine)2-(carboxylate)4 ligand set characteristic of the diiron-carboxylate family. Iron ligand Asp (D) 84 was replaced with Glu (E), the amino acid found in the cognate position of the structurally similar diiron-carboxylate protein, methane monooxygenase hydroxylase (MMOH). We previously showed that this substitution allows accumulation of a mu-1,2-peroxodiiron(III) intermediate, which does not accumulate in the wild-type (wt) protein and is probably a structural homologue of intermediate P (H(peroxo)) in O2 activation by MMOH. In addition, the near-surface residue Trp (W) 48 was replaced with Phe (F), blocking transfer of the "extra" electron that occurs in wt R2 during formation of the formally Fe(III)Fe(IV) cluster X. Decay of the mu-1,2-peroxodiiron(III) complex in R2-W48F/D84E gives an initial brown product, which contains very little Y122* and which converts very slowly (t1/2 approximately 7 h) upon incubation at 0 degrees C to an intensely purple final product. X-ray crystallographic analysis of the purple product indicates that F208 has undergone epsilon-hydroxylation and the resulting phenol has shifted significantly to become a ligand to Fe2 of the diiron cluster. Resonance Raman (RR) spectra of the purple product generated with 16O2 or 18O2 show appropriate isotopic sensitivity in bands assigned to O-phenyl and Fe-O-phenyl vibrational modes, confirming that the oxygen of the Fe(III)-phenolate species is derived from O2. Chemical analysis, experiments involving interception of the hydroxylating intermediate with exogenous reductant, and M?ssbauer and EXAFS characterization of the brown and purple species establish that F208 hydroxylation occurs during decay of the peroxo complex and formation of the initial brown product. The slow transition to the purple Fe(III)-phenolate species is ascribed to a ligand rearrangement in which mu-O2- is lost and the F208-derived phenolate coordinates. The reprogramming to F208 monooxygenase requires both amino acid substitutions, as very little epsilon-hydroxyphenylalanine is formed and pathways leading to Y122* formation predominate in both R2-D84E and R2-W48F.     

Rapid genotyping of CYP2D6, CYP2C19 and TPMT polymorphisms by primer extension reaction in a dipstick format

In recent years an increasing amount of interest has been directed at the study and routine testing of polymorphisms responsible for variations in drug metabolism. Most of the current methods involve either time-consuming electrophoresis steps or specialized and expensive equipment. In this context, we have developed a rapid, simple and robust method for genotyping of CYP2D6*3, CYP2D6*4, CYP2C19*2, CYP2C19*3 and TPMT*2 single nucleotide polymorphisms (SNP). Genomic DNA is isolated from whole blood and the segments that span the SNP of interest are amplified by PCR. The products are subjected directly (without purification) to two primer extension (PEXT) reactions (three cycles each) using normal and mutant primers in the presence of biotin-dUTP. The PEXT primers contain a (dA)₃₀ segment at the 5′ end. The PEXT products are detected visually by a dry-reagent dipstick-type assay in which the biotinylated extension products are captured from immobilized streptavidin on the test zone of the strip and detected by hybridization with oligo(dT)-functionalized gold nanoparticles. Patient samples (76 variants in total) were genotyped and the results were fully concordant with those obtained by direct DNA sequencing.     

以Schiff碱N-丁酰基水杨酰肼为配体的三核镍(Ⅱ)配合物的合成与结构

The trinuclear nickel (Ⅱ) complex Ni3(C11H11N2O3)2(C5H5N)4 was prepared by the reaction of Ni(OAc)2 @ 4H2Owith N-butylsalicylhydrazide Schiff base ligand and characterized by X-ray crystallography. The crystal belongs totriclinic, Mr=930.91, space group P1 with cell parameters a =9. 8489(2)A, b = 12. 3110(2) A, c = 18. 4035(3) A, α=71.353(2)°, β=76.638(2)°, γ=84.815(2)°, V=2056.72(6)A3, Z=2, Dc=1.503g@cm-3,μ(MoKα) = 1. 417mm-1, F(000) = 964, R = 0. 0317, wR = 0. 0868. A total of 5393 independent reflectionswere collected, of which 4448 reflections with I ≥ 2 σ(I) were observed. There are two centrosymmetrical trin-angle of 180°. The interatomic distances between the nickel atoms on the two sides are 9. 2030(8) A and 9. 1876(9) A for the two molecules, respectively. The central nickel atom adopts an axially elongated octahedral geometrywhereas the nickel atoms on the two sides have square-planar coordination environment. CCDC: 194083.     

The role of plastic beta-hairpin and weak hydrophobic core in the stability and unfolding of a full sequence design protein

In this study, the thermal stability of a designed alpha/beta protein FSD (full sequence design) was studied by explicit solvent simulations at three moderate temperatures, 273 K, 300 K, and 330 K. The average properties of the ten trajectories at each temperature were analyzed. The thermal unfolding, as judged by backbone root-mean-square deviation and percentage of native contacts, was displayed with increased sampling outside of the native basin as the temperature was raised. The positional fluctuation of the hairpin residues was significantly higher than that of the helix residues at all three temperatures. The hairpin segment displayed certain plasticity even at 273 K. Apart from the terminal residues, the highest fluctuation was shown in the turn residues 7-9. Secondary structure analysis manifested the structural heterogeneity of the hairpin segment. It was also revealed by the simulation that the hydrophobic core was vulnerable to thermal denaturation. Consistent with the experiment, the I7Y mutation in the double mutant FSD-EY (FSD with mutations Q1E and I7Y) dramatically increased the protein stability in the simulation, suggesting that the plasticity of the hairpin can be partially compensated by a stronger hydrophobic core. As for the unfolding pathway, the breathing of the hydrophobic core and the separation of the two secondary structure elements (alpha helix and beta hairpin) was the initiation step of the unfolding. The loss of global contacts from the separation further destabilized the hairpin structure and also led to the unwinding of the helix.