Dissecting the Cytochrome P450 1A2- and 3A4-Mediated Metabolism of Aflatoxin B1 in Ligand and Protein Contributions.

Aflatoxin B1 (AFB1) is a chemically intriguing compound because it has several potential sites of metabolism (SOMs), although only some of them are observed experimentally. Cytochrome P450 (CYP) 3A4 and 1A2 are the major isoforms involved in its metabolism. Here, we systematically investigate reactivity and accessibility of all possible SOMs in these two CYPs to elucidate AFB1 metabolism. DFT calculations were used to determine activation energies for each possible reaction. Aliphatic hydroxylation on position 9A and 3α are energetically favored, whereas position 9 is the preferred site for epoxidation. Docking studies, molecular dynamics (MD) simulations, and free energy (MM/GBSA) calculations were applied to elucidate the accessibility of each SOM. The most stable binding modes in CYP3A4 favor the formation of the 3α-hydroxylated and 8,9-exo-epoxide metabolites. Conversion of the methoxy group is also sterically possible, but not observed experimentally due to its low reactivity. In the CYP1A2 active site, AFB1 cannot orient position 3 towards the catalytic center, whereas the 8,9-exo-epoxide and 9A-hydroxylated metabolites are formed from the most stable and the 8,9-endo-epoxide from a less stable binding mode, respectively. The results agree with experimental data and suggest that both reactivity and the shape of the enzyme active site need to be considered to understand the distribution of SOMs and to improve current SOM prediction methods.

[1]  Giulio Rastelli,et al.  Application of a post-docking procedure based on MM-PBSA and MM-GBSA on single and multiple protein conformations. , 2012, European journal of medicinal chemistry.

[2]  Scott P. Brown,et al.  Large-scale application of high-throughput molecular mechanics with Poisson-Boltzmann surface area for routine physics-based scoring of protein-ligand complexes. , 2009, Journal of medicinal chemistry.

[3]  G. Ciccotti,et al.  Numerical Integration of the Cartesian Equations of Motion of a System with Constraints: Molecular Dynamics of n-Alkanes , 1977 .

[4]  Sason Shaik,et al.  A density functional study of the factors that influence the regioselectivity of toluene hydroxylation by cytochrome P450 enzymes , 2007 .

[5]  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.

[6]  H. Berendsen,et al.  Molecular dynamics with coupling to an external bath , 1984 .

[7]  T. Poulos,et al.  Dissecting cytochrome P450 3A4-ligand interactions using ritonavir analogues. , 2013, Biochemistry.

[8]  A. Klamt,et al.  COSMO : a new approach to dielectric screening in solvents with explicit expressions for the screening energy and its gradient , 1993 .

[9]  Yasuo Shinohara,et al.  Flexible structure of cytochrome P450: promiscuity of ligand binding in the CYP3A4 heme pocket. , 2009, Anticancer research.

[10]  Chris Oostenbrink,et al.  Cooperative binding of aflatoxin B1 by cytochrome P450 3A4: a computational study. , 2014, Chemical research in toxicology.

[11]  Parr,et al.  Development of the Colle-Salvetti correlation-energy formula into a functional of the electron density. , 1988, Physical review. B, Condensed matter.

[12]  Stewart B Kirton,et al.  Prediction of binding modes for ligands in the cytochromes P450 and other heme‐containing proteins , 2005, Proteins.

[13]  Charlotte M. Deane,et al.  Rapid, Precise, and Reproducible Prediction of Peptide-MHC Binding Affinities from Molecular Dynamics That Correlate Well with Experiment. , 2015, Journal of chemical theory and computation.

[14]  Adrian J Mulholland,et al.  Understanding the determinants of selectivity in drug metabolism through modeling of dextromethorphan oxidation by cytochrome P450 , 2011, Proceedings of the National Academy of Sciences.

[15]  T. Poulos,et al.  Structural and Mechanistic Insights into the Interaction of Cytochrome P4503A4 with Bromoergocryptine, a Type I Ligand* , 2011, The Journal of Biological Chemistry.

[16]  Ulf Ryde,et al.  The reaction mechanism of iron and manganese superoxide dismutases studied by theoretical calculations , 2006, J. Comput. Chem..

[17]  Hans W. Horn,et al.  Fully optimized contracted Gaussian basis sets for atoms Li to Kr , 1992 .

[18]  David E. Gloriam,et al.  SMARTCyp: A 2D Method for Prediction of Cytochrome P450-Mediated Drug Metabolism. , 2010, ACS medicinal chemistry letters.

[19]  F. Guengerich,et al.  2. Cytochromes P450 , 2012 .

[20]  A. Becke.  A New Mixing of Hartree-Fock and Local Density-Functional Theories , 1993 .

[21]  P Willett,et al.  Development and validation of a genetic algorithm for flexible docking. , 1997, Journal of molecular biology.

[22]  Lars Olsen,et al.  Prediction of activation energies for hydrogen abstraction by cytochrome p450. , 2006, Journal of medicinal chemistry.

[23]  Christine M. Bathelt,et al.  Aromatic hydroxylation by cytochrome P450: model calculations of mechanism and substituent effects. , 2003, Journal of the American Chemical Society.

[24]  Chris Oostenbrink,et al.  Catalytic site prediction and virtual screening of cytochrome P450 2D6 substrates by consideration of water and rescoring in automated docking. , 2006, Journal of medicinal chemistry.

[25]  S. Genheden,et al.  The MM/PBSA and MM/GBSA methods to estimate ligand-binding affinities , 2015, Expert opinion on drug discovery.

[26]  E. Gallagher,et al.  The kinetics of aflatoxin B1 oxidation by human cDNA-expressed and human liver microsomal cytochromes P450 1A2 and 3A4. , 1996, Toxicology and applied pharmacology.

[27]  L. Olsen,et al.  Use of density functional theory in drug metabolism studies , 2014, Expert opinion on drug metabolism & toxicology.

[28]  D. Kroetz,et al.  Structure-function relationships of human liver cytochromes P450 3A: aflatoxin B1 metabolism as a probe. , 1998, Biochemistry.

[29]  A. Becke.  Density-functional thermochemistry. III. The role of exact exchange , 1993 .

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

[31]  Sason Shaik,et al.  A predictive pattern of computed barriers for C-h hydroxylation by compound I of cytochrome p450. , 2004, Journal of the American Chemical Society.

[32]  U. Zanger.  Introduction to Drug Metabolism , 2012 .

[33]  Tingjun Hou,et al.  Assessing the performance of the molecular mechanics/Poisson Boltzmann surface area and molecular mechanics/generalized Born surface area methods. II. The accuracy of ranking poses generated from docking , 2011, J. Comput. Chem..

[34]  Sason Shaik,et al.  The "Rebound Controversy": An Overview and Theoretical Modeling of the Rebound Step in C-H Hydroxylation by Cytochrome P450 , 2004 .

[35]  Zhen Zhang,et al.  Aflatoxin B1-induced hepatocellular carcinoma in developing countries: Geographical distribution, mechanism of action and prevention , 2013, Oncology letters.

[36]  Chris Oostenbrink,et al.  Prediction of cytochrome P450 mediated metabolism. , 2015, Advanced drug delivery reviews.

[37]  Thomas E. Cheatham,et al.  Quantum mechanically derived AMBER‐compatible heme parameters for various states of the cytochrome P450 catalytic cycle , 2012, J. Comput. Chem..

[38]  Andreas Bender,et al.  Computational Prediction of Metabolism: Sites, Products, SAR, P450 Enzyme Dynamics, and Mechanisms , 2012, J. Chem. Inf. Model..

[39]  Holger Gohlke,et al.  The Amber biomolecular simulation programs , 2005, J. Comput. Chem..

[40]  T. Sjögren,et al.  Structural basis for ligand promiscuity in cytochrome P450 3A4 , 2006, Proceedings of the National Academy of Sciences.

[41]  David E. Gloriam,et al.  The SMARTCyp cytochrome P450 metabolism prediction server , 2010, Bioinform..

[42]  C David Stout,et al.  Adaptations for the Oxidation of Polycyclic Aromatic Hydrocarbons Exhibited by the Structure of Human P450 1A2*♦ , 2007, Journal of Biological Chemistry.

[43]  S. Shaik,et al.  Multi-state epoxidation of ethene by cytochrome P450: a quantum chemical study. , 2001, Journal of the American Chemical Society.

[44]  Youyong Li,et al.  Assessing the performance of MM/PBSA and MM/GBSA methods. 4. Accuracies of MM/PBSA and MM/GBSA methodologies evaluated by various simulation protocols using PDBbind data set. , 2014, Physical chemistry chemical physics : PCCP.

[45]  Peter A. Kollman,et al.  AMBER, a package of computer programs for applying molecular mechanics, normal mode analysis, molecular dynamics and free energy calculations to simulate the structural and energetic properties of molecules , 1995 .

[46]  Lars Olsen,et al.  Prediction of activation energies for aromatic oxidation by cytochrome P450. , 2008, The journal of physical chemistry. A.

[47]  Christine M. Bathelt,et al.  Quantum Mechanics/Molecular Mechanics Modeling of Regioselectivity of Drug Metabolism in Cytochrome P450 2C9 , 2013, Journal of the American Chemical Society.