Molecular modeling and identification of substrate binding site of orphan human cytochrome P450 4F22

Cytochrome P450s are superfamily of heme proteins which generally monooxygenate hydrophobic compounds. The human cytochrome P450 4F22 (CYP4F22) was categorized into "orphan" CYPs because of its unknown function. CYP4F22 is a potential drug target for cancer therapy. However, three-dimensional structure, the active site topology and substrate specificity of CYP4F22 remain unclear. In this study, a three-dimensional model of human P450 4F22 was constructed by comparative modeling using Modeller 9v5. The resulting model was refined by energy minimization subjected to the quality assessment from both geometric and energetic aspects and was found to be of reasonable quality. Docking approach was employed to dock arachidonic acid into the active site of CYP4F22 in order to probe the ligand-binding modes. As a result, several key residues were identified to be responsible for the binding of arachidonic acid with CYP4F22. These findings provide useful information for understanding the biological roles of CYP4F22 and structure-based drug design.     

Inhibition of CYP2B4 by the mechanism-based inhibitor 2-ethynylnaphthalene: inhibitory potential of 2EN is dependent on the size of the substrate

2-Ethynylnaphthalene (2EN) is a mechanism-based inhibitor of CYP2B4 with two components to the inhibition, (1) enzyme inactivation, which requires covalent binding of the 2EN metabolite, and (2) reversible inhibition by 2EN itself. Both inhibitory components were examined using several different CYP2B4 substrates. Preincubation of CYP2B4 with 2EN led to a time-dependent inactivation of each of the CYP2B4-dependent activities examined; however, the ability of 2EN to reversibly inhibit CYP2B4 depended on the substrate employed, which is inconsistent with classical inhibition patterns. The degree 2EN's reversible inhibition was shown not to correlate with the substrate affinity for the active site, but with parameters related to the molecular size of the substrate. The results are consistent with 2EN and the smaller substrates simultaneously fitting in the CYP2B4 active site, leading to very little inhibition. Larger substrates exhibited greater degrees of inhibition because of their inability to co-bind with inhibitor in the active site.     

Essential requirements for substrate binding affinity and selectivity toward human CYP2 family enzymes

A detailed analysis of substrate selectivity within the cytochrome P450 2 (CYP2) family is reported. From a consideration of specific interactions between drug substrates for human CYP2 family enzymes and the putative active sites of CYP2A6, CYP2B6, CYP2C8, CYP2C9, CYP2C19, CYP2D6, and CYP2E1, it is likely that the number and disposition of hydrogen bond donor/acceptors and aromatic rings within the various P450 substrate molecules determines their enzyme selectivity and binding affinity, together with directing their preferred routes of metabolism by the CYP2 enzymes concerned. Although many aliphatic residues are present in most P450 active sites, it would appear that their main contribution centers around hydrophobic interactions and desolvation processes accompanying substrate binding. Molecular modeling studies based on the recent CYP2C5 crystal structure appear to show close agreement with site-directed mutagenesis experiments and with information on substrate metabolism and selectivity within the CYP2 family.     

Inhibition of CYP2B4 by 2-ethynylnaphthalene: evidence for the co-binding of substrate and inhibitor within the active site

2-Ethynylnaphthalene (2EN) is an effective mechanism-based inhibitor of CYP2B4. There are two inhibitory components: (1) irreversible inactivation of CYP2B4 (a typical time-dependent inactivation), and (2) a reversible component. The reversible component was unusual in that the degree of inhibition was not simply a characteristic of the enzyme-inhibitor interaction, but dependent on the size of the substrate molecule used to monitor residual activity. The effect of 2EN on the metabolism of seven CYP2B4 substrates showed that it was not an effective reversible inhibitor of substrates containing a single aromatic ring; substrates with two fused rings were competitively inhibited by 2EN; and larger substrates were non-competitively inhibited. Energy-based docking studies demonstrated that, with increasing substrate size, the energy of 2EN and substrate co-binding in the active site became unfavorable precisely at the point where 2EN became a competitive inhibitor. Hierarchical docking revealed potential allosteric inhibition sites separate from the substrate binding site.     

Structure–function analysis of cytochromes P450 2B

In the last 4 years, breakthroughs were made in the field of P450 2B (CYP2B) structure–function through determination of one ligand-free and two inhibitor-bound X-ray crystal structures of CYP2B4, which revealed many of the structural features required for binding ligands of different size and shape. Large conformational changes of several plastic regions of CYP2B4 can dramatically reshape the active site of the enzyme to fit the size and shape of the bound ligand without perturbing the overall P450 fold. Solution biophysical studies using isothermal titration calorimetry (ITC) have revealed the large difference in the thermodynamic parameters of CYP2B4 in binding inhibitors of different ring chemistry and side chains. Other studies have revealed that the effects of site-specific mutations on steady-state kinetic parameters and mechanism-based inactivation are often substrate dependent. These findings agree with the structural data that the enzymes adopt different conformations to bind various ligands. Thus, the substrate specificity of an individual enzyme is determined not only by active site residues but also non-active site residues that modulate conformational changes that are important for substrate access and rearrangement of the active site to accommodate the bound substrate.     

In vitro metabolism of fipronil by human and rat cytochrome P450 and its interactions with testosterone and diazepam

Fipronil (5-amino-1-[2,6-dichloro-4-(trifluoromethyl)phenyl]-4-[(trifluoromethyl)sulfinyl]-1H-pyrazole-3-carbonitrile) is a highly active, broad spectrum insecticide from the phenyl pyrazole family, which targets the gamma-amino butyric acid (GABA) receptor. Although fipronil is presently widely used as an insecticide and acaricide, little information is available with respect to its metabolic fate and disposition in mammals. This study was designed to investigate the in vitro human metabolism of fipronil and to examine possible metabolic interactions that fipronil may have with other substrates. Fipronil was incubated with human liver microsomes (HLM) and several recombinant cytochrome P450 (CYP) isoforms obtained from BD Biosciences. HPLC was used for metabolite identification and quantification. Fipronil sulfone was the predominant metabolite via CYP oxidation. The K(m) and V(max) values for human liver microsomes are 27.2 microM and 0.11 nmol/mg proteinmin, respectively; for rat liver microsomes (RLM) the K(m) and V(max) are 19.9 microM and 0.39 nmol/mg proteinmin, respectively. CYP3A4 is the major isoform responsible for fipronil oxidation in humans while CYP2C19 is considerably less active. Other human CYP isoforms have minimal or no activity toward fipronil. Co-expression of cytochrome b(5) (b(5)) is essential for CYP3A4 to manifest high activity toward fipronil. Ketoconazole, a specific inhibitor of CYP3A4, inhibits 78% of the HLM activity toward fipronil at a concentration of 2 microM. Oxidative activity toward fipronil in 19 single-donor HLMs correlated well with their ability to oxidize testosterone. The interactions of fipronil and other CYP3A4 substrates, such as testosterone and diazepam, were also investigated. Fipronil metabolism was activated by testosterone in HLM but not in CYP3A4 Supersomes. Testosterone 6beta-hydroxylation in HLM was inhibited by fipronil. Fipronil inhibited diazepam demethylation but had little effect on diazepam hydroxylation. The results suggest that fipronil has the potential to interact with a wide range of xenobiotics or endogenous chemicals that are CYP3A4 substrates and that fipronil may be a useful substrate for the characterization of CYP3A4 in HLM.     

Finding homes for orphan cytochrome P450s: CYP4V2 and CYP4F22 in disease states

The cytochrome P450 (CYP) 4 family of enzymes contains several recently identified membersthat are referred to as “orphan P450s” because their endogenous substrates are unknown.Human CYP4V2 and CYP4F22 are two such orphan P450s that are strongly linked to ocular andskin disease, respectively. Genetic analyses have identified a wide spectrum of mutations in the CYP4V2gene from patients suffering from Bietti’s crystalline corneoretinal dystrophy, and mutations in theCYP4F22 gene have been linked to lamellar ichthyosis. The strong gene–disease associations provideunique opportunities for elucidating the substrate specificity of these orphan P450s and unraveling thebiochemical pathways that may be impacted in patients with CYP4V2 and CYP4F22 functional deficits.     

Heterologous expression of CYP102A5 variant from Bacillus cereus CYPPB-1: Validation of model for predicting drug metabolism of human P450 probe substrates

CYP102A5 variant (ADL27534) from isolated Bacillus cereus CYPPB-1 was heterologously expressed in Escherichia coli Top 10 cells. Comparative sequence analysis of purified CYP102A5 variant with respect to reported CYP102A5 (AAP10153) from Bacillus cereus ATCC 14579 revealed amino acid sequence changes at positions P245S and M318I of heme domain. The binding affinities of 15 selected human P450 probe substrates towards isolated CYP102A5 were analyzed in silico using a homology model together with molecular docking techniques to predict the human drug metabolism. In vitro analysis suggested that the purified CYP102A5 metabolizes typical substrates of human CYP2C9, CYP2D6, CYP2E1, and CYP3A4, such as coumarin, propranolol, aniline, chlorzoxazone, p-nitrophenol, and nifedipine. The calculated K _M values for propranolol, chloroxazone, coumarin, aniline, and 4-nitrophenol were calculated to be 0.962?±?0.041, 1.254?±?0.057, 2.859?±?0.083, 2.732?±?0.106, and 2.528?±?0.11 mM, respectively. Importantly, taking a ChemScore cutoff value of ?31 kJ/mol, substrate binding at active site and in vitro activity as the distinguishing lines between “substrates” and “nonsubstrates” revealed one false-positive and one false-negative results out of the 15 compounds examined. This is the first report on validation of CYP102A family homology model for in silico prediction of human drug metabolism.     

Thermodynamic fidelity of the mammalian cytochrome P450 2B4 active site in binding substrates and inhibitors

Structural plasticity of mammalian cytochromes P450 (CYP) has recently been explored in our laboratory and elsewhere to understand the ligand-binding promiscuity. CYP2B4 exhibits very different conformations and thermodynamic signatures in binding the small inhibitor 4-(4-chlorophenyl)imidazole (4-CPI) versus the large bifonazole. Using four key active-site mutants (F296A, T302A, I363A, and V367L) that are involved in binding one or both inhibitors, we dissected the thermodynamic basis for the ability of CYP2B4 to bind substrates and inhibitors of different sizes and chemistry. In all cases, 1:1 binding stoichiometry was observed. The inhibitors 4-CPI, 1-(4-chlorophenyl)imidazole, and 1-(2-(benzyloxy)ethyl)imidazole bind to the mutants with a free energy difference (ΔΔG) of ∼ 0.5 to 1 kcal/mol compared with the wild type but with a large entropy-enthalpy compensation of up to 50 kcal/mol. The substrate testosterone binds to all four mutants with a ΔΔG of ∼ 0.5 kcal/mol but with as much as 40 kcal/mol of entropy-enthalpy compensation. In contrast, benzphetamine binding to V367L and F296A is accompanied by a ΔΔG of ∼ 1.5 and 3 kcal/mol, respectively. F296A, I363A, and V367L exhibit very different benzphetamine metabolite profiles, indicating the different substrate-binding orientations in the active site of each mutant. Overall, the findings indicate that malleability of the active site allows mammalian P450s to exhibit a high degree of thermodynamic fidelity in ligand binding.     

Acceptor specificities and selective inhibition of recombinant human Gal- and GlcNAc-transferases that synthesize core structures 1, 2, 3 and 4 of O-glycans

Background

Modifications of proteins by O-glycosylation determine many of the properties and functions of proteins. We wish to understand the mechanisms of O-glycosylation and develop inhibitors that could affect glycoprotein functions and alter cellular behavior.

Methods

We expressed recombinant soluble human Gal- and GlcNAc-transferases that synthesize the O-glycan cores 1 to 4 and are critical for the overall structures of O-glycans. We determined the properties and substrate specificities of these enzymes using synthetic acceptor substrate analogs. Compounds that were inactive as substrates were tested as inhibitors.

Results

Enzymes significantly differed in their recognition of the sugar moieties and aglycone groups of substrates. Core 1 synthase was active with glycopeptide substrates but GlcNAc-transferases preferred substrates with hydrophobic aglycone groups. Chemical modifications of the acceptors shed light on enzyme–substrate interactions. Core 1 synthase was weakly inhibited by its substrate analog benzyl 2-butanamido-2-deoxy-α-d-galactoside while two of the three GlcNAc-transferases were selectively and potently inhibited by bis-imidazolium salts which are not substrate analogs.

Conclusions

This work delineates the distinct specificities and properties of the enzymes that synthesize the common O-glycan core structures 1 to 4. New inhibitors were found that could selectively inhibit the synthesis of cores 1, 2 and 3 but not core 4.

General significance

These studies help our understanding of the mechanisms of action of enzymes critical for O-glycosylation. The results may be useful for the re-engineering of O-glycosylation to determine the roles of O-glycans and the enzymes critical for O-glycosylation, and for biotechnology with potential therapeutic applications.     

Homology modeling and substrate binding study of human CYP2C18 and CYP2C19 enzymes

It is well established that the variable binding-site architecture and composition of the P450 metabolizing heme proteins are major modulators of substrate and product specificity. Even the three closely related human liver isozymes, CYP2C9, CYP2C18, and CYP2C19, do not share all substrates and do not always lead to the same preferred hydroxylation products. The lack of knowledge of their three-dimensional (3D) structures has hindered efforts to understand the differences in their specificities. Building on previous work for the CYP2C9 enzyme, 3D models of CYP2C18 and 2C19 have been constructed and validated by computational methods developed and tested in our laboratory. They were used to characterize explicit enzyme-substrate complexes using the isoform-specific substrates progesterone and (S)-mephenytoin for 2C19 and 2-[2,3-dichloro-4-(3-hydroxypropyloxy)benzoyl]thiophene for 2C18. The results allowed both common and unique binding-site residues to be identified in each model. The calculated preferred hydroxylation site was obtained for each substrate and was found to be consistent with experimental observation. Comparisons were made among the 2C9, 2C18, and 2C19 model binding sites to investigate the subtle differences among them. These models can be used as structure-based guides for mutagenesis studies and screening of potential pharmaceuticals or toxins.     

Distinct binding of cholesterol and 5beta-cholestane-3alpha,7alpha,12alpha-triol to cytochrome P450 27A1: evidence from modeling and site-directed mutagenesis studies

Cytochrome P450 27A1 (P450 27A1 or CYP27A1) is an important enzyme that participates in different pathways of cholesterol degradation as well as in the activation of vitamin D(3). Several approaches were utilized to investigate how two physiological substrates, cholesterol and 5beta-cholestane-3alpha,7alpha,12alpha-triol, interact with CYP27A1. The enzyme active site was first probed spectrally by assessing binding of the two substrates and five substrate analogues followed by computer modeling and site-directed mutagenesis. The computer models suggest that the spatial positions and orientations of cholesterol and 5beta-cholestane-3alpha,7alpha,12alpha-triol are different in the enzyme active site. As a result, some of the active site residues interact with both substrates, although they are situated differently relative to each steroid, and some residues bind only one substrate. Mutation of the overlapping substrate-contact residues (W100, H103, T110, M301C, V367, I481, and V482) affected CYP27A1 binding and enzyme activity in a substrate-dependent manner and allowed identification of several important side chains. T110 is proposed to interact with the 12alpha-hydroxyl of 5beta-cholestane-3alpha,7alpha,12alpha-triol, whereas V367 seems to be crucial for correct positioning of the cholesterol C26 methyl group and for regioselective hydroxylation of this substrate. Distinct binding of the CYP27A1 substrates may provide insight into why phenotypic manifestations of cerebrotendinous xanthomatosis, a disease associated with CYP27A1 deficiency, are so diverse.     

Substrate preferences and catalytic parameters determined by structural characteristics of sterol 14alpha-demethylase (CYP51) from Leishmania infantum

Leishmaniasis is a major health problem that affects populations of ～90 countries worldwide, with no vaccine and only a few moderately effective drugs. Here we report the structure/function characterization of sterol 14α-demethylase (CYP51) from Leishmania infantum. The enzyme catalyzes removal of the 14α-methyl group from sterol precursors. The reaction is essential for membrane biogenesis and therefore has great potential to become a target for antileishmanial chemotherapy. Although L. infantum CYP51 prefers C4-monomethylated sterol substrates such as C4-norlanosterol and obtusifoliol (V(max) of ～10 and 8 min(-1), respectively), it is also found to 14α-demethylate C4-dimethylated lanosterol (V(max) = 0.9 min(-1)) and C4-desmethylated 14α-methylzymosterol (V(max) = 1.9 min(-1)). Binding parameters with six sterols were tested, with K(d) values ranging from 0.25 to 1.4 μM. Thus, L. infantum CYP51 is the first example of a plant-like sterol 14α-demethylase, where requirements toward the composition of the C4 atom substituents are not strict, indicative of possible branching in the postsqualene portion of sterol biosynthesis in the parasite. Comparative analysis of three CYP51 substrate binding cavities (Trypanosoma brucei, Trypanosoma cruzi, and L. infantum) suggests that substrate preferences of plant- and fungal-like protozoan CYP51s largely depend on the differences in the enzyme active site topology. These minor structural differences are also likely to underlie CYP51 catalytic rates and drug susceptibility and can be used to design potent and specific inhibitors.     

Characterization of differences in substrate specificity among CYP1A1, CYP1A2 and CYP1B1: an integrated approach employing molecular docking and molecular dynamics simulations

Recent trends in new drug discovery of anticancer drugs have made oncologists more aware of the fact that the new drug discovery must target the developing mechanism of tumorigenesis to improve the therapeutic efficacy of antineoplastic drugs. The drugs designed are expected to have high affinity towards the novel targets selectively. Current research highlights overexpression of CYP450s, particularly cytochrome P450 1A1 (CYP1A1), in tumour cells, representing a novel target for anticancer therapy. However, the CYP1 family is identified as posing significant problems in selectivity of anticancer molecules towards CYP1A1. Three members have been identified in the human CYP1 family: CYP1A1, CYP1A2 and CYP1B1. Although sequences of the three isoform have high sequence identity, they have distinct substrate specificities. The understanding of macromolecular features that govern substrate specificity is required to understand the interplay between the protein function and dynamics, design novel antitumour compounds that could be specifically metabolized by only CYP1A1 to mediate their antitumour activity and elucidate the reasons for differences in substrate specificity profile among the three proteins. In the present study, we employed a combination of computational methodologies: molecular docking and molecular dynamics simulations. We utilized eight substrates for elucidating the difference in substrate specificity of the three isoforms. Lastly, we conclude that the substrate specificity of a particular substrate depends upon the type of the active site residues, the dynamic motions in the protein structure upon ligand binding and the physico‐chemical characteristics of a particular ligand. Copyright © 2016 John Wiley & Sons, Ltd.     

The characterisation of two members of the cytochrome P450 CYP150 family: CYP150A5 and CYP150A6 from Mycobacterium marinum

BackgroundActinobacteria, including the Mycobacteria, have a large component of cytochrome P450 family monooxygenases. This includes Mycobacterium tuberculosis, M. ulcerans and M. marinum, and M. vanbaalenii. These enzymes can abstract CH bonds and have important roles in natural product biosynthesis.MethodsTwo members of the bacterial CYP150 family, CYP150A5 and CYP150A6 from M. marinum, were produced, purified and characterised. The potential substrate ranges of both enzymes were analysed and the monooxygenase activity of CYP150A5 was reconstituted using a physiological electron transfer partner system. CYP150A6 was structurally characterised by X-ray crystallography.ResultsCYP150A5 was shown to bind various norisoprenoids and terpenoids. It could regioselectively hydroxylate β-ionol. The X-ray crystal structure of substrate-free CYP150A6 was solved to 1.5 Å. This displayed an open conformation with short F and G helices, an unresolved F-G loop region and exposed active site pocket. The active site residues could be identified and important variations were found among the CYP150A enzymes. Haem-binding azole inhibitors were identified for both enzymes.ConclusionsThe structure of CYP150A6 will facilitate the identification of physiological substrates and the design of better inhibitors for members of this P450 family. Based on the observed differences in substrate binding preference and sequence variations among the active site residues, their roles are predicted to be different.General significanceMultiple CYP150 family members were found in many bacteria and are prevalent in the Mycobacteria including several human pathogens. Inhibition and structural data are reported here for these enzymes for the first time.     

Analysis of human cytochrome P450 2C8 substrate specificity using a substrate pharmacophore and site-directed mutants

The structural determinants of substrate specificity of human liver cytochrome P450 2C8 (CYP2C8) were investigated using site-directed mutants chosen on the basis of a preliminary substrate pharmacophore and a three-dimensional (3D) model. Analysis of the structural features common to CYP2C8 substrates exhibiting a micromolar K(m) led to a substrate pharmacophore in which the site of oxidation by CYP2C8 is 12.9, 8.6, 4.4, and 3.9 A from features that could establish ionic or hydrogen bonds, and hydrophobic interactions with protein amino acid residues. Comparison of this pharmacophore with a 3D model of CYP2C8 constructed using the X-ray structure of CYP2C5 suggested potential CYP2C8 amino acid residues that could be involved in substrate recognition. Twenty CYP2C8 site-directed mutants were constructed and expressed in yeast to compare their catalytic activities using five CYP2C8 substrates that exhibit different structures and sizes [paclitaxel, fluvastatin, retinoic acid, a sulfaphenazole derivative (DMZ), and diclofenac]. Mutation of arginine 241 had marked effects on the hydroxylation of anionic substrates of CYP2C8 such as retinoic acid and fluvastatin. Serine 100 appears to be involved in hydrogen bonding interactions with a polar site of the CYP2C8 substrate pharmacophore, as shown by the 3-4-fold increase in the K(m) of paclitaxel and DMZ hydroxylation after the S100A mutation. Residues 114, 201, and 205 are predicted to be in close contact with substrates, and their mutations lead either to favorable hydrophobic interactions or to steric clashes with substrates. For instance, the S114F mutant was unable to catalyze the 6alpha-hydroxylation of paclitaxel. The S114F and F205A mutants were the best catalysts for retinoic acid and paclitaxel (or fluvastatin) hydroxylation, respectively, with k(cat)/K(m) values 5 and 2.1 (or 2.4) times higher, respectively, than those found for CYP2C8. Preliminary experiments of docking of the substrate into the experimentally determined X-ray structure of substrate-free CYP2C8, which became available quite recently [Schoch, G. A., et al. (2004) J. Biol. Chem. 279, 9497], were consistent with key roles for S100, S114, and F205 residues in substrate binding. The results suggest that the effects of mutation of arginine 241 on anionic substrate hydroxylation could be indirect and result from alterations of the packing of helix G with helix B'.     

A shuffled CYP2C library with a high degree of structural integrity and functional versatility

Cytochrome P450 (CYP) enzymes involved in mammalian xenobiotic metabolism are attractive targets for the engineering of biocatalysts since they have broad and overlapping substrate and reaction substrate specificities. In this report, a library of chimeric mutants was prepared from CYP2C8, CYP2C9, CYP2C18 and CYP2C19 by DNA family shuffling. Twelve randomly selected clones were fully sequenced and showed 9 ± 2 crossovers and 1.5 ± 0.5 spontaneous mutations per ∼1.5 kbp open reading frame. CYP hemoprotein expression was observed in 50% (microaerobic culture) to 54% (aerobic culture) of clones. The functional diversity of the library was assessed using three luminogenic substrates, diclofenac and indole as probe substrates. A random sample of 26 clones revealed two clones with activity towards luciferin ME, one towards luciferin H and five towards diclofenac 4′-hydroxylation. One mutant showed activity towards all three substrates. Of 96 clones screened on solid media, one showed elevated indigo production compared to the parental forms. Turnover rates for luciferin ME and H metabolism by CYP2C9 and mutants were at least one order of magnitude higher in experiments with membranes compared to whole cells, consistent with impaired product egress from cells. Apparent K_m values were increased in whole cell incubations with luciferin H suggesting impaired access of the substrate to the active site of the enzymes in whole cells. Finally screening with a panel of CYP2C ligands using CYP2C9 or active mutants revealed different patterns of inhibition and heteroactivation of metabolism of luciferin analogs.     

The Tomato yellow leaf curl virus (TYLCV) V2 protein interacts with the host papain-like cysteine protease CYP1

The V2 protein of Tomato yellow leaf curl geminivirus (TYLCV) is an RNA-silencing suppressor that counteracts the innate immune response of the host plant. However, this anti-host defense function of V2 may include targeting of other defensive mechanisms of the plant. Specifically, we show that V2 recognizes and directly binds the tomato CYP1 protein, a member of the family of papain-like cysteine proteases which are involved in plant defense against diverse pathogens. This binding occurred both in vitro and in vivo, within living plant cells. The V2 binding site within mCYP1 was identified in the direct proximity to the papain-like cysteine protease active site.     

π–π Stacking mediated drug–drug interactions in human CYP2E1

Because of having many low molecular mass substrates, CYP2E1 is of particular interests to the pharmaceutical industry. Many evidences showed that this enzyme can adopt multiple substrates to significantly reduce the oxidation rate of the substrates. The detailed mechanism for this observation is still unclear. In the current study, we employed GPU‐accelerated molecular dynamics simulations to study the multiple‐binding mode of human CYP2E1, with an aim of offering a mechanistic explanation for the unexplained multiple‐substrate binding. Our results showed that Thr303 and Phe478 were key factors for the substrate recognition and multiple‐substrate binding. The former can form a significant hydrogen bond to recognize and position the substrate in the productive binding orientation in the active site. The latter acted as a mediator for the substrate communications via π–π stacking interactions. In the multiple‐binding mode, the aforementioned π–π stacking interactions formed by the aromatic rings of both substrates and Phe478 drove the first substrate far away from the catalytic center, orienting in an additional binding position and going against the substrate metabolism. All these findings could give atomic insights into the detailed mechanism for the multiple‐substrate binding in human CYP2E1, providing useful information for the drug metabolism mechanism and personalized use of clinical drugs. Proteins 2013; © 2012 Wiley Periodicals, Inc.