Abstract
In our previous kinetics studies the natural products oroxylin and wogonin were shown to have strong biological affinity for, and inhibitory effects against, human cytochrome P450 1A2, with IC50 values of 579 and 248 nM, respectively; this might lead to the occurrence of drug–drug interactions when co-administered clinically. However, their inhibitory mechanisms against 1A2 remain elusive. In this study, molecular docking and molecular dynamics simulations were performed to better understand the molecular basis of their inhibitory mechanisms towards 1A2. Structural analysis revealed that oroxylin has a different binding pattern from wogonin and another very strongly binding inhibitor α-naphthoflavone (ANF, IC50 = 49 nM). The O7 atom of oroxylin forms hydrogen bonds with the OD1/OD2 atoms of Asp313, which is not observed in the 1A2–wogonin complex. Because of energetically unfavorable repulsions with the methoxy group at the 6 position of the oroxylin ring, significant conformational changes were observed for the sidechain of Thr118 in the MD simulated model. As a result, the larger and much more open binding-site architecture of the 1A2–oroxylin complex may account for its weaker inhibitory effect relative to the 1A2–ANF complex. Energy analysis indicated that oroxylin has a less negative predicted binding free energy of −19.8 kcal/mol than wogonin (−21.1 kcal/mol), which is consistent with our experimental assays. Additionally, our energy results suggest that van der Waals/hydrophobic and hydrogen-bonding interactions are important in the inhibitory mechanisms of oroxylin whereas the former is the underlying force responsible for strong inhibition by ANF and wogonin.
![](http://media.springernature.com/m312/springer-static/image/art%3A10.1007%2Fs00249-011-0785-1/MediaObjects/249_2011_785_Fig1_HTML.gif)
![](http://media.springernature.com/m312/springer-static/image/art%3A10.1007%2Fs00249-011-0785-1/MediaObjects/249_2011_785_Fig2_HTML.gif)
![](http://media.springernature.com/m312/springer-static/image/art%3A10.1007%2Fs00249-011-0785-1/MediaObjects/249_2011_785_Fig3_HTML.gif)
![](http://media.springernature.com/m312/springer-static/image/art%3A10.1007%2Fs00249-011-0785-1/MediaObjects/249_2011_785_Fig4_HTML.gif)
![](http://media.springernature.com/m312/springer-static/image/art%3A10.1007%2Fs00249-011-0785-1/MediaObjects/249_2011_785_Fig5_HTML.gif)
![](http://media.springernature.com/m312/springer-static/image/art%3A10.1007%2Fs00249-011-0785-1/MediaObjects/249_2011_785_Fig6_HTML.gif)
Similar content being viewed by others
References
Arimoto R (2006) Computational models for predicting interaction with cytochrome P450 enzyme. Curr Top Med Chem 6:1909–1918
Babu TH, Rao VRS, Tiwari AK, Babu KS, Srinivas PV, Ali AZ, Rao JM (2008) Synthesis and biological evaluation of novel 8-aminomethylated oroxylin A analogues as alpha-glucosidase inhibitors. Bioorg Med Chem Lett 18:1659–1662
Berendsen HJC, Postma JPM, van Gunsteren WF, DiNola A, Haak JR (1984) Molecular dynamics with coupling to an external bath. J Chem Phys 81:3684–3690
Birt DF, Hendrich S, Wang W (2001) Dietary agents in cancer prevention: flavonoids and isoflavonoids. Pharmacol Ther 90:157–177
Butler MA, Iwasaki M, Guengerich FP, Kadlubar FF (1989) Human cytochrome P-450PA (P-4501A2), the phenacetin O-deethylase, is primarily responsible for the hepatic 3-demethylation of caffeine and N-oxidation of carcinogenic arylamines. Proc Natl Acad Sci USA 86:7696–7700
Case DA, Darden TA, Cheatham TE III et al (2008) AMBER 10. University of California, San Francisco
Cho US, Park EY, Dong MS, Park BS, Kim K, Kim KH (2003) Tight-binding inhibition by alpha-naphthoflavone of human cytochrome P450 1A2. Biochim Biophys Acta 1648:195–202
Darden T, York D, Pedersen L (1993) Particle mesh Ewald: an N log(N) method for Ewald sums in large systems. J Chem Phys 98:10089–10092
Eaton DL, Gallagher EP, Bammler TK, Kunze KL (1995) Role of cytochrome P450 1A2 in chemical carcinogenesis: implications for human variability in expression and enzyme activity. Pharmacogenet Genom 5:259–274
Faber MS, Jetter A, Fuhr U (2005) Assessment of 1A2 activity in clinical practice: why, how, and when? Basic Clin Pharmacol Toxicol 97:125–134
Fang L, Zhang H, Cui W, Ji MJ (2008) Studies of the mechanism of selectivity of protein tyrosine phosphatase 1B (PTP1B) bidentate inhibitors using molecular dynamics simulations and free energy calculations. J Chem Info Model 48:2030–2041
Fogolari F, Brigo A, Molinari H (2003) Protocol for MM/PBSA molecular dynamics simulations of proteins. Biophys J 85:159–166
Frisch MJ, Trucks GW, Schlegel HB et al (2004) Gaussian 03, Revision E01. Gaussian, Inc, Pittsburgh PA
Genheden S, Ryde U (2010) How to obtain statistically converged MM/GBSA results. J Comput Chem 31:837–846
Harris DL, Park JY, Gruenke L, Waskell L (2004) Theoretical study of the ligand-CYP2B4 complexes: effect of structure on binding free energies and heme spin state. Proteins 55:895–914
He F, Bi HC, Xie ZY, Zuo Z, Li JK, Li X, Zhao LZ, Chen X, Huang M (2007) Rapid determination of six metabolites from multiple cytochrome P450 probe substrates in human liver microsome by liquid chromatography/mass spectrometry: application to high-throughput inhibition screening of terpenoids. Rapid Commun Mass Spectrom 21:635–643
He L, He F, Bi H, Li J, Zeng S, Luo HB, Huang M (2010) Isoform-selective inhibition of chrysin towards human cytochrome P450 1A2. Kinetics analysis, molecular docking, and molecular dynamics simulations. Bioorg Med Chem Lett 20:6008–6012
Heim KE, Tagliaferro AR, Bobilya DJ (2002) Flavonoid antioxidants: chemistry, metabolism and structure-activity relationships. J Nutr Biochem 13:572–584
Hong CC, Tang BK, Hammond GL, Tritchler D, Yaffe M, Boyd NF (2004) Cytochrome P450 1A2 (CYP1A2) activity and risk factors for breast cancer: a cross-sectional study. Breast Cancer Res 6:R352–R365
Ioria F, da Fonsecab R, Joao Ramosb M, Menziani MC (2005) Theoretical quantitative structure—activity relationships of flavone ligands interacting with cytochrome P450 1A1 and 1A2 isozymes. Bioorg Med Chem 13:4366–4374
Jain AN (2007) Surflex-Dock 2.1: robust performance from ligand energetic modeling, ring flexibility, and knowledge-base search. J Comput Aid Mol Des 21:281–306
Kemp CA, Flanagan JU, van Eldik AJ, Maréchal JD, Wolf CR, Roberts GCK, Paine MJI, Sutcliffe MJ (2004) Validation of model of cytochrome P450 2D6: an in silico tool for predicting metabolism and inhibition. J Med Chem 47:5340–5346
Kunze LK, William FT (1993) Isoform-selective mechanism-based inhibition of human cytochrome P450 1A2 by furafylline. Chem Res Toxicol 6:649–656
Labute P (2008) The generalized Born/volume integral implicit solvent model: estimation of the free energy of hydration using London dispersion instead of atomic surface area. J Comput Chem 29(10):1693–1698
Li HB, Chen F (2005) Isolation and puritication of baicalein, wogonin and oroxylin a from the medicinal plant Scutellaria baicalensis by high-speed counter-current chromatography. J Chromatogr A 1074:107–110
Li JK, He F, Bi HC, Zuo Z, Liu BD, Luo HB, Huang M (2008) Inhibition of human cytochrome P450 CYP 1A2 by flavonoids: a quantitative structure-activity relationship study. Acta Pharm Sin 43:1198–1204
Lin CC, Shieh DE (1996) The anti-inflammatory activity of Scutellaria rivularis extracts and its active components, baicalin, baicalein, and wogonin. Am J Chin Med 24:31–36
Liu M, Yuan M, Luo M, Bu X, Luo HB, Hu X (2010) Binding of curcumin with glyoxalase I: molecular docking, molecular dynamics simulations, and kinetics analysis. Biophys Chem 147:28–34
Li-Weber M (2009) New therapeutic aspects of flavones: the anticancer properties of Scutellaria and its main active constituents wogonin, baicalein and baicalin. Cancer Treat Rev 35:57–68
Miyamoto S, Kollman PA (1992) Settle: an analytical version of the SHAKE and RATTLE algorithm for rigid water models. J Comput Chem 13:8952–8962
MOE 2008.10 (2008) Chemical Computing Group Inc., Montreal, Quebec, Canada
Moon YJ, Wang X, Morris ME (2006) Dietary flavonoids: effects on xenobiotic and carcinogen metabolism. Toxicol In Vitro 20:187–210
Neugebauer RC, Uchiechowska U, Meier R et al (2008) Structure-activity studies on splitomicin derivatives as sirtuin inhibitors and computational prediction of binding mode. J Med Chem 51:1203–1213
Park H, Lee S, Suh J (2005) Structural and dynamical basis of broad substrate specificity, catalytic mechanism, and inhibition of cytochrome P450 3A4. J Am Chem Soc 127:13634–13642
Rastelli G, Del RA, Degliesposti G, Sgobba M (2010) Fast and accurate predictions of binding free energies using MM-PBSA and MM-GBSA. J Comput Chem 31:797–810
Sanner MF, Olson AJ, Spehner JC (1996) Reduced surface: an efficient way to compute molecular surfaces. Biopolymers 38:305–320
Sansen S, Yano JK, Reynald RL, Schoch GA, Griffin KJ, Stout CD, Johnson EF (2007) Adaptations for the oxidation of polycyclic aromatic hydrocarbons exhibited by the structure of human P450 1A2. J Biol Chem 282:14348–14355
Si DY, Wang Y, Zhou YH, Guo Y, Wang J, Zhou H, Li ZS, Fawcett JP (2009) Mechanism of CYP2C9 inhibition by flavones and flavonols. Drug Metab Disp 37:629–634
Still WC, Tempczyk A, Hawley RC, Hendrickson T (1990) Semianalytical treatment of solvation for molecular mechanics and dynamics. J Am Chem Soc 112:6127–6129
Wang J, Wolf RM, Caldwell JW, Kollman PA, Case DA (2004) Development and testing of a general AMBER force field. J Comput Chem 25:1157–1174
Wang J, Wang W, Kollman PA, Case DA (2006) Automatic atom type and bond type perception in molecular mechanical calculations. J Mol Graph Model 25:247–260
Wolf CR, Smith G, Smith RL (2000) Science, medicine and the future pharmacogenetics. Br Med J 320:987–990
Zeng J, Li W, Zhao Y, Liu G, Tang Y, Jiang H (2008) Insights into ligand selectivity in estrogen receptor isoforms: molecular dynamics simulations and binding free energy calculations. J Phys Chem B 112:2719–2726
Acknowledgments
We are grateful for financial support from the Natural Science Foundation of China (21103234), the Natural Science Foundation of Guangdong Province (S2011030003190), the Science Foundation of the Department of Education in Guangdong Province (CXZD1006), the Science Foundation of Guangzhou City (2010Y1-C531), and Fundamental Research Funds for the Central Universities (10ykjc20).
Author information
Authors and Affiliations
Corresponding author
Additional information
Yong-Xian Shao and Peng Zhao contributed equally.
Electronic supplementary material
Below is the link to the electronic supplementary material.
Rights and permissions
About this article
Cite this article
Shao, YX., Zhao, P., Li, Z. et al. The molecular basis for the inhibition of human cytochrome P450 1A2 by oroxylin and wogonin. Eur Biophys J 41, 297–306 (2012). https://doi.org/10.1007/s00249-011-0785-1
Received:
Revised:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s00249-011-0785-1