COMPARISON OF METHODS FOR THE PREDICTION OF THE METABOLIC SITES FOR CYP3A4-MEDIATED METABOLIC REACTIONS

Predictions of the metabolic sites for new chemical entities, synthesized or only virtual, are important in the early phase of drug discovery to guide chemistry efforts in the synthesis of new compounds with reduced metabolic liability. This information can now be obtained from in silico predictions, and therefore, a thorough and unbiased evaluation of the computational techniques available is needed. Several computational methods to predict the metabolic hot spots are emerging. In this study, metabolite identification using MetaSite and a docking methodology, GLUE, were compared. Moreover, the published CYP3A4 crystal structure and computed CYP3A4 homology models were compared for their usefulness in predicting metabolic sites. A total of 227 known CYP3A4 substrates reported to have one or more metabolites adding up to 325 metabolic pathways were analyzed. Distance-based fingerprints and four-point pharmacophore derived from GRID molecular interaction fields were used to characterize the substrate and protein in MetaSite and the docking methodology, respectively. The CYP3A4 crystal structure and homology model with the reactivity factor enabled achieved a similar prediction success (78%) using the MetaSite method. The docking method had a relatively lower prediction success (∼57% for the homology model), although it still may provide useful insights for interactions between ligand and protein, especially for uncommon reactions. The MetaSite methodology is automated, rapid, and has relatively accurate predictions compared with the docking methodology used in this study.

[1]  F. Guengerich,et al.  Cytochrome P-450 3A4: regulation and role in drug metabolism. , 1999, Annual review of pharmacology and toxicology.

[2]  B. Ring,et al.  Comparative metabolic capabilities of CYP3A4, CYP3A5, and CYP3A7. , 2002, Drug metabolism and disposition: the biological fate of chemicals.

[3]  G. Cruciani,et al.  MetaSite: understanding metabolism in human cytochromes from the perspective of the chemist. , 2005, Journal of medicinal chemistry.

[4]  D. Goodsell,et al.  Automated docking to multiple target structures: Incorporation of protein mobility and structural water heterogeneity in AutoDock , 2002, Proteins.

[5]  Francesca Fanelli,et al.  Theoretical investigation of substrate specificity for cytochromes P450 IA2, P450 IID6 and P450 IIIA4 , 2000, J. Comput. Aided Mol. Des..

[6]  Gabriele Cruciani,et al.  Predicting drug metabolism: a site of metabolism prediction tool applied to the cytochrome P450 2C9. , 2003, Journal of medicinal chemistry.

[7]  Stewart B Kirton,et al.  Impact of incorporating the 2C5 crystal structure into comparative models of cytochrome P450 2D6 , 2002, Proteins.

[8]  J. Halpert,et al.  Analysis of homotropic and heterotropic cooperativity of diazepam oxidation by CYP3A4 using site-directed mutagenesis and kinetic modeling. , 2003, Archives of biochemistry and biophysics.

[9]  P. Goodford A computational procedure for determining energetically favorable binding sites on biologically important macromolecules. , 1985, Journal of medicinal chemistry.

[10]  R. Sheridan,et al.  A model for predicting likely sites of CYP3A4-mediated metabolism on drug-like molecules. , 2003, Journal of medicinal chemistry.

[11]  Eric F. Johnson,et al.  Structure of Mammalian Cytochrome P450 2B4 Complexed with 4-(4-Chlorophenyl)imidazole at 1.9-Å Resolution , 2004, Journal of Biological Chemistry.

[12]  Eric F. Johnson,et al.  The Structure of Human Microsomal Cytochrome P450 3A4 Determined by X-ray Crystallography to 2.05-Å Resolution* , 2004, Journal of Biological Chemistry.

[13]  G. Cruciani,et al.  An Innovative Application of the “Flexible” GRID/PCA Computational Method: Study of Differences in Selectivity between PGAs from Escherichia coli and a Providentia rettgeri Mutant , 2004, Biotechnology progress.

[14]  J Berendzen,et al.  The catalytic pathway of cytochrome p450cam at atomic resolution. , 2000, Science.

[15]  J. Halpert,et al.  Phenylalanine and tryptophan scanning mutagenesis of CYP3A4 substrate recognition site residues and effect on substrate oxidation and cooperativity. , 2001, Biochemistry.

[16]  A. Alex,et al.  A novel approach to predicting P450 mediated drug metabolism. CYP2D6 catalyzed N-dealkylation reactions and qualitative metabolite predictions using a combined protein and pharmacophore model for CYP2D6. , 1999, Journal of medicinal chemistry.

[17]  Eric F. Johnson,et al.  Structure of mammalian cytochrome P450 2B4 , 2000 .

[18]  Slobodan Petar Rendic,et al.  Human cytochrome P450 enzymes: a status report summarizing their reactions, substrates, inducers, and inhibitors. , 1997, Drug metabolism reviews.

[19]  Wei Tang,et al.  Heterotropic cooperativity of cytochrome P450 3A4 and potential drug-drug interactions. , 2001, Current drug metabolism.

[20]  M Pastor,et al.  A novel strategy for improving ligand selectivity in receptor-based drug design. , 1995, Journal of medicinal chemistry.

[21]  Ramon Carbó-Dorca,et al.  Fitted electronic density functions from H to Rn for use in quantum similarity measures: cis-diamminedichloroplatinum(II) complex as an application example , 1999, J. Comput. Chem..

[22]  Jose Cosme,et al.  Crystal Structures of Human Cytochrome P450 3A4 Bound to Metyrapone and Progesterone , 2004, Science.

[23]  Chris de Graaf,et al.  Binding mode prediction of cytochrome p450 and thymidine kinase protein-ligand complexes by consideration of water and rescoring in automated docking. , 2005, Journal of medicinal chemistry.

[24]  C David Stout,et al.  Structure of mammalian cytochrome P450 2C5 complexed with diclofenac at 2.1 A resolution: evidence for an induced fit model of substrate binding. , 2003, Biochemistry.

[25]  A. Galetin,et al.  Multisite kinetic analysis of interactions between prototypical CYP3A4 subgroup substrates: midazolam, testosterone, and nifedipine. , 2003, Drug metabolism and disposition: the biological fate of chemicals.