A Negative Cooperativity Mechanism of Human CYP2E1 Inferred from Molecular Dynamics Simulations and Free Energy Calculations

Human cytochrome P450 2E1 (CYP2E1) participates in the metabolism of over 2% of all the oral drugs. A hallmark peculiar feature of this enzyme is that it exhibits a pronounced negative cooperativity in substrate binding. However the mechanism by which the negative cooperativity occurs is unclear. Here, we performed molecular dynamics simulations and free energy calculations on human CYP2E1 to examine the structural differences between the substrate-free and the enzymes with one and two aniline molecules bound. Our results indicate that although the effector substrate does not bind in the active site cavity, it still can directly interact with the active site residues of human CYP2E1. The interaction of the effector substrate with the active site leads to a reorientation of active site residues, which thereby weakens the interactions of the active substrate with this site. We also identify a conserved residue T303 that plays a crucial role in the negative cooperative binding on the short-range effects. This residue is a key factor in the positioning of substrates and in proton delivery to the active site. Additionally, a long-range effect of the effector substrate is identified in which F478 is proposed to play a key role. As located in the interface between the active and effector sites, this residue structurally links the active and effector sites and is found to play a significant role in affecting substrate access and ligand positioning within the active site. In the negative cooperative binding, this residue can decrease the interactions of the active substrate with the active site by π-π stacking which then lowers the hydroxylation activity for the active substrate. These findings are in agreement with previous experimental observations and thus provide detailed atomistic insight into the poorly understood mechanism of the negative cooperativity in human CYP2E1.

[1]  M. J. Coon,et al.  Epoxidation of olefins by cytochrome P450: evidence from site-specific mutagenesis for hydroperoxo-iron as an electrophilic oxidant. , 1998, Proceedings of the National Academy of Sciences of the United States of America.

[2]  C. Ioannides,et al.  Sex differences in the diabetes-induced modulation of rat hepatic cytochrome P450 proteins. , 1993, Biochemical pharmacology.

[3]  T. Goosen,et al.  Differential effects of naturally occurring isothiocyanates on the activities of cytochrome P450 2E1 and the mutant P450 2E1 T303A. , 2001, Archives of biochemistry and biophysics.

[4]  D. Koop.  Hydroxylation of p-nitrophenol by rabbit ethanol-inducible cytochrome P-450 isozyme 3a. , 1986, Molecular pharmacology.

[5]  J. Liu,et al.  Analysis of four residues within substrate recognition site 4 of human cytochrome P450 3A4: role in steroid hydroxylase activity and alpha-naphthoflavone stimulation. , 1998, Archives of biochemistry and biophysics.

[6]  G. Miller.  Advances in the interpretation and prediction of CYP2E1 metabolism from a biochemical perspective , 2008, Expert opinion on drug metabolism & toxicology.

[7]  Ruth Nussinov,et al.  Structural dynamics of the cooperative binding of organic molecules in the human cytochrome P450 3A4. , 2007, Journal of the American Chemical Society.

[8]  Danielson Pb,et al.  The cytochrome P450 superfamily: biochemistry, evolution and drug metabolism in humans. , 2002 .

[9]  M. Komori,et al.  Replacement of Thr-303 of P450 2E1 with serine modifies the regioselectivity of its fatty acid hydroxylase activity. , 1993, Journal of biochemistry.

[10]  G A Petsko,et al.  Aromatic-aromatic interaction: a mechanism of protein structure stabilization. , 1985, Science.

[11]  P. Danielson,et al.  The cytochrome P450 superfamily: biochemistry, evolution and drug metabolism in humans. , 2002, Current drug metabolism.

[12]  G. McGaughey,et al.  pi-Stacking interactions. Alive and well in proteins. , 1998, The Journal of biological chemistry.

[13]  David S. Goodsell,et al.  A semiempirical free energy force field with charge‐based desolvation , 2007, J. Comput. Chem..

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

[15]  C. Lieber,et al.  Induction of cytochrome P450IIE1 in the obese overfed rat. , 1991, Molecular pharmacology.

[16]  F. Guengerich,et al.  Catalytic selectivity of human cytochrome P450 enzymes: relevance to drug metabolism and toxicity. , 1994, Toxicology letters.

[17]  Sason Shaik,et al.  The elusive oxidant species of cytochrome P450 enzymes: characterization by combined quantum mechanical/molecular mechanical (QM/MM) calculations. , 2002, Journal of the American Chemical Society.

[18]  D. Trafalis,et al.  CYP2E1 and risk of chemically mediated cancers , 2010, Expert opinion on drug metabolism & toxicology.

[19]  J. Halpert,et al.  Analysis of Differential Substrate Selectivities of CYP2B6 and CYP2E1 by Site-Directed Mutagenesis and Molecular Modeling , 2003, Journal of Pharmacology and Experimental Therapeutics.

[20]  Kuo-Chen Chou,et al.  Insight into the molecular switch mechanism of human Rab5a from molecular dynamics simulations. , 2009, Biochemical and biophysical research communications.

[21]  Tao Zhang,et al.  Advances in human cytochrome p450 and personalized medicine. , 2011, Current drug metabolism.

[22]  W. Atkins,et al.  Multiple-ligand binding in CYP2A6: probing mechanisms of cytochrome P450 cooperativity by assessing substrate dynamics. , 2008, Biochemistry.

[23]  Ryan M. Laddusaw,et al.  CYP2E1 Substrate Inhibition , 2008, Journal of Biological Chemistry.

[24]  Barry Honig,et al.  Extending the Applicability of the Nonlinear Poisson−Boltzmann Equation: Multiple Dielectric Constants and Multivalent Ions† , 2001 .

[25]  A. Rettie,et al.  Evaluation of atypical cytochrome P450 kinetics with two-substrate models: evidence that multiple substrates can simultaneously bind to cytochrome P450 active sites. , 1998, Biochemistry.

[26]  K. Korzekwa,et al.  Activation of CYP3A4: evidence for the simultaneous binding of two substrates in a cytochrome P450 active site. , 1994, Biochemistry.

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

[28]  Advances in the interpretation and prediction of CYP2E1 metabolism from a biochemical perspective. , 2008, Expert opinion on drug metabolism & toxicology.

[29]  Richard M. Jackson,et al.  Q-SiteFinder: an energy-based method for the prediction of protein-ligand binding sites , 2005, Bioinform..

[30]  A. Nanji,et al.  Markedly enhanced cytochrome P450 2E1 induction and lipid peroxidation is associated with severe liver injury in fish oil-ethanol-fed rats. , 1994, Alcoholism, clinical and experimental research.

[31]  D. Lewis,et al.  Cytochromes P 450 in the Bioactivation of Chemicals , 2004 .

[32]  Kuo-Chen Chou,et al.  Insights from Modeling the 3D Structure of New Delhi Metallo-β-Lactamse and Its Binding Interactions with Antibiotic Drugs , 2011, PloS one.

[33]  Sason Shaik,et al.  Mechanism of oxidation reactions catalyzed by cytochrome p450 enzymes. , 2004, Chemical reviews.

[34]  David S. Goodsell,et al.  Automated docking using a Lamarckian genetic algorithm and an empirical binding free energy function , 1998, J. Comput. Chem..

[35]  D S Goodsell,et al.  Automated docking of flexible ligands: Applications of autodock , 1996, Journal of molecular recognition : JMR.

[36]  D. Lewis,et al.  Cytochromes P450 in the bioactivation of chemicals. , 2004, Current topics in medicinal chemistry.

[37]  Kuo-Chen Chou,et al.  Molecular modeling of cytochrome P450 and drug metabolism. , 2010, Current drug metabolism.

[38]  Ying Wang,et al.  Molecular Dynamics Studies on T1 Lipase: Insight into a Double-Flap Mechanism , 2010, J. Chem. Inf. Model..

[39]  Kuo-Chen Chou,et al.  An Allosteric Mechanism Inferred from Molecular Dynamics Simulations on Phospholamban Pentamer in Lipid Membranes , 2011, PloS one.

[40]  Kuo-Chen Chou,et al.  Insights from studying the mutation-induced allostery in the M2 proton channel by molecular dynamics. , 2010, Protein engineering, design & selection : PEDS.

[41]  M. Ingelman-Sundberg,et al.  Ethanol-, fasting-, and acetone-inducible cytochromes P-450 in rat liver: regulation and characteristics of enzymes belonging to the IIB and IIE gene subfamilies. , 1988, Biochemistry.

[42]  H Koga,et al.  Uncoupling of the cytochrome P-450cam monooxygenase reaction by a single mutation, threonine-252 to alanine or valine: possible role of the hydroxy amino acid in oxygen activation. , 1989, Proceedings of the National Academy of Sciences of the United States of America.

[43]  J. Ponder,et al.  Force fields for protein simulations. , 2003, Advances in protein chemistry.

[44]  Kuo-Chen Chou,et al.  Insights from investigating the interactions of adamantane-based drugs with the M2 proton channel from the H1N1 swine virus. , 2009, Biochemical and biophysical research communications.

[45]  K. Battaile,et al.  Human Cytochrome P450 2E1 Structures with Fatty Acid Analogs Reveal a Previously Unobserved Binding Mode* , 2010, The Journal of Biological Chemistry.

[46]  Kuo-Chen Chou,et al.  Structure of cytochrome p450s and personalized drug. , 2009, Current medicinal chemistry.

[47]  Anna L. Blobaum,et al.  Novel Reversible Inactivation of Cytochrome P450 2E1 T303A by tert-Butyl Acetylene: The Role of Threonine 303 in Proton Delivery to the Active Site of Cytochrome P450 2E1 , 2004, Journal of Pharmacology and Experimental Therapeutics.

[48]  M. Komori,et al.  Different mechanisms of regioselection of fatty acid hydroxylation by laurate (omega-1)-hydroxylating P450s, P450 2C2 and P450 2E1. , 1994, Journal of Biochemistry (Tokyo).

[49]  E. Scott,et al.  Structures of Human Cytochrome P-450 2E1 , 2008, Journal of Biological Chemistry.

[50]  M. J. Coon,et al.  Peroxo-iron and oxenoid-iron species as alternative oxygenating agents in cytochrome P450-catalyzed reactions: switching by threonine-302 to alanine mutagenesis of cytochrome P450 2B4. , 1996, Proceedings of the National Academy of Sciences of the United States of America.

[51]  E. Alexov,et al.  Combining conformational flexibility and continuum electrostatics for calculating pK(a)s in proteins. , 2002, Biophysical journal.