Crystal structure of H2O2-dependent cytochrome P450SPalpha with its bound fatty acid substrate: insight into the regioselective hydroxylation of fatty acids at the alpha position

Cytochrome P450_SPα (CYP152B1) isolated from Sphingomonas paucimobilis is the first P450 to be classified as a H₂O₂-dependent P450. P450_SPα hydroxylates fatty acids with high α-regioselectivity. Herein we report the crystal structure of P450_SPα with palmitic acid as a substrate at a resolution of 1.65 Å. The structure revealed that the C_α of the bound palmitic acid in one of the alternative conformations is 4.5 Å from the heme iron. This conformation explains the highly selective α-hydroxylation of fatty acid observed in P450_SPα. Mutations at the active site and the F–G loop of P450_SPα did not impair its regioselectivity. The crystal structures of mutants (L78F and F288G) revealed that the location of the bound palmitic acid was essentially the same as that in the WT, although amino acids at the active site were replaced with the corresponding amino acids of cytochrome P450_BSβ (CYP152A1), which shows β-regioselectivity. This implies that the high regioselectivity of P450_SPα is caused by the orientation of the hydrophobic channel, which is more perpendicular to the heme plane than that of P450_BSβ.     

Computer modeling of cytochrome P450 2E1 three-dimensional structure

A computer model of human cytochrome P450 2E1 (CYP2E1) three-dimensional structure and active site was constructed based on homology with crystallographic coordinates of CYP2C5 and CYP2C9. A high degree of secondary structure homology for human, mouse, rat and rabbit CYP2E1 was demonstrated. The location of heme and the supporting alpha-helices was established. CYP2E1, CYP2C5 and CYP2C9 active sites are distinguished by pocket size and their amino acid residues composition. Key amino acid residues forming the active site channel and substrate-binding cavity are presented. Active site surface area and volume for CYP2E1, CYP2C5 and CYP2C9 were calculated.     

The Role of Ile476 in the Structural Stability and Substrate Binding of Human Cytochrome P450 2C8

The biological function and stability of a cytochrome P450 (CYP) mainly depend on the subtle properties of the residues in the active site cavity, which are generally more divergent among proteins than other parts of the protein. As the most unique member of human CYP2C family, CYP2C8 has an isoleucine (Ile) 476 instead of phenylalanine (Phe) in substrate recognizing site 6 (SRS6). However, the role of Ile476 of CYP2C8 is still unknown. Therefore, six site-directed mutants of CYP2C8 were constructed to better define this. By UV–visible and circular dichroism spectroscopy studies, we studied for the first time the structural stability and all-trans-retinoic acid binding capability of the CYP2C8 variants. We found that the ferric CYP2C8 went through three states during thermal unfolding. Combined with substrate binding studies, our data revealed that residue 476 was involved in contact with substrate and was important for maintaining the thermal stability of CYP2C8.     

Structural analysis of mammalian cytochrome P450 2B4 covalently bound to the mechanism-based inactivator tert-butylphenylacetylene: insight into partial enzymatic activity

A combined structural and computational analysis of rabbit cytochrome P450 2B4 covalently bound to the mechanism-based inactivator tert-butylphenylacetylene (tBPA) has yielded insight into how the enzyme retains partial activity. Since conjugation to tBPA modifies a highly conserved active site residue, the residual activity of tBPA-labeled 2B4 observed in previous studies was puzzling. Here we describe the first crystal structures of a modified mammalian P450, which show an oxygenated metabolite of tBPA conjugated to Thr 302 of helix I. These results are consistent with previous studies that identified Thr 302 as the site of conjugation. In each structure, the core of 2B4 remains unchanged, but the arrangement of plastic regions differs. This results in one structure that is compact and closed. In this conformation, tBPA points toward helix B', making a 31° angle with the heme plane. This conformation is in agreement with previously performed in silico experiments. However, dimerization of 2B4 in the other structure, which is caused by movement of the B/C loop and helices F through G, alters the position of tBPA. In this case, tBPA lies almost parallel to the heme plane due to the presence of helix F' of the opposite monomer entering the active site to stabilize the dimer. However, docking experiments using this open form show that tBPA is able to rotate upward to give testosterone and 7-ethoxy-4-trifluoromethylcoumarin access to the heme, which could explain the previously observed partial activity.     

Targeting of the highly conserved threonine 302 residue of cytochromes P450 2B family during mechanism-based inactivation by aryl acetylenes

Cytochromes P450 (CYPs or P450s) contain a highly conserved threonine residue in the active site, which is referred to as Thr302 in the amino acid sequence of CYP2B4. Extensive biochemical and crystallographic studies have established that this Thr302 plays a critical role in activating molecular oxygen to generate Compound I, a putative iron(IV)-oxo porphyrin cation radical, that carries out the preliminary oxygenation of CYP substrates. Because of its proximity to the center of the P450 active site, this Thr302 is susceptible to mechanism-based inactivation under certain conditions. In this article, we review recent studies on the mechanism-based inactivation of three mammalian P450s in the 2B family, CYP2B1 (rat), 2B4 (rabbit) and 2B6 (human) by tert-butylphenylacetylene (tBPA). These studies showed that tBPA is a potent mechanism-based inactivator of CYP2B1, 2B4 and 2B6 with high k_inact/K_I ratios (0.23–2.3 min⁻¹ μM⁻¹) and low partition ratios (0–5). Furthermore, mechanistic studies revealed that tBPA inactivates these three CYP2B enzymes through the formation of a single ester adduct with the Thr302 in the active site. These inhibitory properties of tBPA allowed the preparation of a modified CYP2B4 where the Thr302 was covalently and stoichiometrically labeled by a reactive intermediate of tBPA in quantities large enough to permit spectroscopic and crystallographic studies of the consequences of covalent modification of Thr302. Molecular modeling studies revealed a unique binding mode of tBPA in the active site that may shed light on the potency of this inhibition. The results from these studies may serve as a basis for designing more specific and potent inhibitors for P450s by targeting this highly conserved threonine residue which is present in the active sites of most mammalian P450s.     

Impact of incorporating the 2C5 crystal structure into comparative models of cytochrome P450 2D6

Cytochrome P450 2D6 (CYP2D6) metabolizes approximately one third of the drugs in current clinical use. To gain insight into its structure and function, we have produced four different sets of comparative models of 2D6: one based on the structures of P450s from four different microorganisms (P450 terp, P450 eryF, P450 cam, and P450 BM3), another on the only mammalian P450 (2C5) structure available, and the other two based on alternative amino acid sequence alignments of 2D6 with all five of these structures. Principal component analysis suggests that inclusion of the 2C5 crystal structure has a profound effect on the modeling process, altering the general topology of the active site, and that the models produced differ significantly from all of the templates. The four models of 2D6 were also used in conjunction with molecular docking to produce complexes with the substrates codeine and 1-methyl-4-phenyl-1,2,3,6-tetrahydropyridine (MPTP); this identified Glu 216 [in the F-helix; substrate recognition site (SRS) 2] as a key determinant in the binding of the basic moiety of the substrate. Our studies suggest that both Asp 301 and Glu 216 are required for metabolism of basic substrates. Furthermore, they suggest that Asp 301 (I-helix, SRS-4), a residue thought from mutagenesis studies to bind directly to the basic moiety of substrates, may play a key role in positioning the B'-C loop (SRS-1) and that the loss of activity on mutating Asp 301 may therefore be the result of an indirect effect (movement of the B'-C loop) on replacing this residue.     

Structure of the human lung cytochrome P450 2A13

The human lung cytochrome P450 2A13 (CYP2A13) activates the nicotine-derived procarcinogen 4-(methylnitrosamino)-1-(3-pyridyl)-1-butanone (NNK) into DNA-altering compounds that cause lung cancer. Another cytochrome P450, CYP2A6, is also present in human lung, but at much lower levels. Although these two enzymes are 93.5% identical, CYP2A13 metabolizes NNK with much lower K(m) values than does CYP2A6. To investigate the structural differences between these two enzymes the structure of CYP2A13 was determined to 2.35A by x-ray crystallography and compared with structures of CYP2A6. As expected, the overall CYP2A13 and CYP2A6 structures are very similar with an average root mean square deviation of 0.5A for the Calpha atoms. Like CYP2A6, the CYP2A13 active site cavity is small and highly hydrophobic with a cluster of Phe residues composing the active site roof. Active site residue Asn(297) is positioned to hydrogen bond with an adventitious ligand, identified as indole. Amino acid differences between CYP2A6 and CYP2A13 at positions 117, 300, 301, and 208 relate to different orientations of the ligand plane in the two protein structures and may underlie the significant variations observed in binding and catalysis of many CYP2A ligands. In addition, docking studies suggest that residues 365 and 366 may also contribute to differences in NNK metabolism.     

Evidence for communality in the primary determinants of CYP74 catalysis and of structural similarities between CYP74 and classical mammalian P450 enzymes

In silico structural analysis of CYP74C3, a membrane-associated P450 enzyme from the plant Medicago truncatula (barrel medic) with hydroperoxide lyase (HPL) specificity, showed that it had strong similarities to the structural folds of the classical microsomal P450 enzyme from rabbits (CYP2C5). It was not only the secondary structure predictions that supported the analysis but site directed mutagenesis of the substrate interacting residues was also consistent with it. This led us to develop a substrate-binding model of CYP74C3 which predicted three amino acid residues, N285, F287, and G288 located in the putative I-helix and distal haem pocket of CYP74C3 to be in close proximity to the preferred substrate 13-HPOTE. These residues were judged to be in equivalent positions to those identified in SRS-4 of CYP2C5. Significance of the residues and their relevance to the model were further assessed by site directed mutagenesis of the three residues followed by EPR spectroscopic and detailed kinetic investigations of the mutated proteins in the presence and absence of detergent. Although point mutation of the residues had no effect on the haem content of the mutated proteins, significant effects on the spin state equilibrium of the haem iron were noted. Kinetic effects of the mutations, which were investigated using three different substrates, were dramatic in nature. In the presence of detergent with the preferred substrate (13-HPOTE), the catalytic center activities and substrate binding affinities of the mutant proteins were reduced by a factor of 8-32 and 4-12, respectively, compared with wild-type--a two orders of magnitude reduction in catalytic efficiencies. We believe this is the first report where primary determinants of catalysis for any CYP74 enzyme, which are fully consistent with our model, have been identified. Our working model predicts that N285 is close enough to suggest that a hydrogen bond with the peroxy group of the enzyme substrate 13-HPOTE is warranted, whereas significance of F287 may arise from a strong hydrophobic interaction between the alkyl group(s) of the substrate and the phenyl ring of F287. We believe that G288 is crucial because of its size. Any other residue with a relatively bulky side chain will hinder the access of substrate to the active site. The effects of the mutations suggests that subtle protein conformational changes in the putative substrate-binding pocket regulate the formation of a fully active monomer-micelle complex with low-spin haem iron and that structural communication exists between the substrate- and micelle-binding sites of CYP74C3. Conservation in CYP74 sequence alignments suggests that N285, F287, and G288 in CYP74C3 and the equivalent residues at positions in other CYP74 enzymes are likely to be critical to catalysis. To support this we show that G324 in CYP74D4 (Arabidopsis AOS), equivalent to G288 in CYP74C3, is a primary determinant of positional specificity. We suggest that the overall structure of CYP74 enzymes is likely to be very similar to those described for classical P450 monooxygenase enzymes.     

Structural Analysis of the Streptomyces avermitilis CYP107W1-Oligomycin A Complex and Role of the Tryptophan 178 Residue

CYP107W1 from Streptomyces avermitilis is a cytochrome P450 enzyme involved in the biosynthesis of macrolide oligomycin A. A previous study reported that CYP107W1 regioselectively hydroxylated C12 of oligomycin C to produce oligomycin A, and the crystal structure of ligand free CYP107W1 was determined. Here, we analyzed the structural properties of the CYP107W1-oligomycin A complex and characterized the functional role of the Trp178 residue in CYP107W1. The crystal structure of the CYP107W1 complex with oligomycin A was determined at a resolution of 2.6 Å. Oligomycin A is bound in the substrate access channel on the upper side of the prosthetic heme mainly by hydrophobic interactions. In particular, the Trp178 residue in the active site intercalates into the large macrolide ring, thereby guiding the substrate into the correct binding orientation for a productive P450 reaction. A Trp178 to Gly mutation resulted in the distortion of binding titration spectra with oligomycin A, whereas binding spectra with azoles were not affected. The Gly178 mutant’s catalytic turnover number for the 12-hydroxylation reaction of oligomycin C was highly reduced. These results indicate that Trp178, located in the open pocket of the active site, may be a critical residue for the productive binding conformation of large macrolide substrates.     

Homology model of human retinoic acid metabolising enzyme cytochrome P450 26A1 (CYP26A1): active site architecture and ligand binding

Homology models of cytochrome P450 RA1 (CYP26A1) were constructed using three human P450 structures, CYP2C8, CYP2C9 and CYP3A4 as templates for the model building. Using MOE software the lowest energy CYP26A1 model was then assessed for stereochemical quality and side chain environment. Further active site optimisation of the CYP26A1 model built using the CYP3A4 template was performed by molecular dynamics to generate a final CYP26A1 model. The natural substrate, all-trans-retinoic acid (atRA), and inhibitor R 15866, were docked into the model allowing further validation of the active site architecture. Using the docking studies structurally and functionally important residues were identified with subsequent characterisation of secondary structure. Multiple hydrophobic interactions, including the side chains of TRP112, PHE299, PHE222, PHE84, PHE374 and PRO371, are important for binding of atRA and R115866. Additional hydrogen bonding interactions were noted as follows: atRA-- C==O of the atRA carboxylate group and ARG86; R115866--benzothiazole nitrogen and the backbone NH of SER115.     

Crystal structures of substrate-free and retinoic acid-bound cyanobacterial cytochrome P450 CYP120A1

The crystal structures of substrate-free and all-trans-retinoic acid-bound CYP120A1 from Synechocystis sp. PCC 6803 were determined at 2.4 and 2.1 A resolution, respectively, representing the first structural characterization of a cyanobacterial P450. Features of CYP120A1 not observed in other P450 structures include an aromatic ladder flanking the channel leading to the active site and a triple-glycine motif within SRS5. Using spectroscopic methods, CYP120A1 is shown to bind 13-cis-retinoic acid, 9-cis-retinoic acid, and retinal with high affinity and dissociation constants of less than 1 microM. Metabolism of retinoic acid by CYP120A1 suggests that CYP120A1 hydroxylates a variety of retinoid derivatives in vivo. On the basis of the retinoic acid-bound CYP120A1 crystal structure, we propose that either carbon 2 or the methyl groups (C16 or C17) of the beta-ionone ring are modified by CYP120A1.     

The numerous effector functions of Nef.

P450BM-3, a catalytically self-sufficient, soluble bacterial P450, contains on the same polypeptide a heme domain and a reductase domain. P450BM-3 catalyzes the oxidation of short- and long-chain, saturated and unsaturated fatty acids. The three-dimensional structure of the heme domain both in the absence and in the presence of fatty acid substrates has been determined; however, the fatty acid in the substrate-bound form is not adequately close to the heme iron to permit a prediction regarding the stereoselectivity of oxidation. In the case of long-chain fatty acids, the products can also serve as substrate and be metabolized several times. In the current study, we have determined the absolute configuration of the three primary products of palmitic acid hydroxylation (15-, 14-, and 13-OH palmitic acid). While the 15- and 14-hydroxy compounds are produced in a highly stereoselective manner (98% R, 2% S), the 13-hydroxy is a mixture of 72% R and 28% S. We have also examined the binding of these three hydroxy acids to P450BM-3 and shown that only two of them (14-OH and 13-OH palmitic acid) can bind to and be further metabolized by P450BM-3. The results indicate that in contrast to the flexibility of palmitoleic acid bound to the oxidized enzyme, palmitic acid is rigidly bound in the active site during catalytic turnover.     

Crystal structure of CYP199A2, a para-substituted benzoic acid oxidizing cytochrome P450 from Rhodopseudomonas palustris

CYP199A2, a cytochrome P450 enzyme from Rhodopseudomonas palustris, oxidatively demethylates 4-methoxybenzoic acid to 4-hydroxybenzoic acid. 4-Ethylbenzoic acid is converted to a mixture of predominantly 4-(1-hydroxyethyl)-benzoic acid and 4-vinylbenzoic acid, the latter being a rare example of CC bond dehydrogenation of an unbranched alkyl group. The crystal structure of CYP199A2 has been determined at 2.0-Å resolution. The enzyme has the common P450 fold, but the B′ helix is missing and the G helix is broken into two (G and G′) by a kink at Pro204. Helices G and G′ are bent back from the extended BC loop and the I helix to open up a clearly defined substrate access channel. Channel openings in this region of the P450 fold are rare in bacterial P450 enzymes but more common in eukaryotic P450 enzymes. The channel is hydrophobic except for the basic residue Arg246 at the entrance, which probably plays a role in the specificity of this enzyme for charged benzoates over neutral phenols and benzenes. The substrate binding pocket is hydrophobic, with Ser97 and Ser247 being the only polar residues. Computer docking of 4-ethylbenzoic acid into the active site suggests that the substrate carboxylate oxygens interact with Ser97 and Ser247, and the β-methyl group is located over the heme iron by Phe185, the side chain of which is only 6.35 Å above the iron in the native structure. This binding orientation is consistent with the observed product profile of exclusive attack at the para substituent. Putidaredoxin of the CYP101A1 system from Pseudomonas putida supports substrate oxidation by CYP199A2 at ∼6% of the activity of the physiological ferredoxin. Comparison of the heme proximal faces of CYP199A2 and CYP101A1 suggests that charge reversal surrounding the surface residue Leu369 in CYP199A2 may be a significant factor in this low cross-activity.     

Modelling of three-dimensional structures of cytochromes P450 11B1 and 11B2.

The final steps of the biosynthesis of glucocorticoids and mineralocorticoids in the adrenal cortex require the action of two different cytochromes P450--CYP11B1 and CYP11B2. Homology modelling of the three-dimensional structures of these cytochromes was performed based on crystallographic coordinates of two bacterial P450s, CYP102 (P450BM-3) and CYP108 (P450terp). Principal attention was given to the modelling of the active sites and a comparison of the active site structures of CYP11B1 and CYP11B2 was performed. It can be demonstrated that key residue contacts within the active site appear to depend on the orientation of the heme. The obtained 3D structures of CYP11B1 and CYP11B2 were used for investigation of structure-function relationships of these enzymes. Previously obtained results on naturally occurring mutants and on mutants obtained by site-directed mutagenesis are discussed.     

Atomic structure of Mycobacterium tuberculosis CYP121 to 1.06 A reveals novel features of cytochrome P450

The first structure of a P450 to an atomic resolution of 1.06 A has been solved for CYP121 from Mycobacterium tuberculosis. A comparison with P450 EryF (CYP107A1) reveals a remarkable overall similarity in fold with major differences residing in active site structural elements. The high resolution obtained allows visualization of several unusual aspects. The heme cofactor is bound in two distinct conformations while being notably kinked in one pyrrole group due to close interaction with the proline residue (Pro(346)) immediately following the heme iron-ligating cysteine (Cys(345)). The active site is remarkably rigid in comparison with the remainder of the structure, notwithstanding the large cavity volume of 1350 A(3). The region immediately surrounding the distal water ligand is remarkable in several aspects. Unlike other bacterial P450s, the I helix shows no deformation, similar to mammalian CYP2C5. In addition, the positively charged Arg(386) is located immediately above the heme plane, dominating the local structure. Putative proton relay pathways from protein surface to heme (converging at Ser(279)) are identified. Most interestingly, the electron density indicates weak binding of a dioxygen molecule to the P450. This structure provides a basis for rational design of putative antimycobacterial agents.     

Homology model of 1alpha,25-dihydroxyvitamin D3 24-hydroxylase cytochrome P450 24A1 (CYP24A1): active site architecture and ligand binding

Homology models of cytochrome P450 24A1 (CYP24A1) were constructed using three human P450 structures, CYP2C8, CYP2C9 and CYP3A4 as templates for the model building. Using molecular operating environment (MOE) software the lowest energy CYP24A1 model was then assessed for stereochemical quality and side chain environment. Further active site optimisation of the CYP24A1 model built using the CYP3A4 template was performed by molecular dynamics to generate a final CYP24A1 model. The natural substrate, 1,25-dihydroxyvitamin D(3) (calcitriol) and the CYP24 inhibitor (R)-N-(2-(1H-imidazol-1-yl)-2-phenylethyl)-4'-chlorobiphenyl-4-carboxamide ((R)-VID-400) were docked into the model allowing further validation of the active site architecture. Using the docking studies structurally and functionally important residues were identified with subsequent characterisation of secondary structure.     

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.     

Human Cytochrome P450 2E1 Structures with Fatty Acid Analogs Reveal a Previously Unobserved Binding Mode

Human microsomal cytochrome P450 (CYP) 2E1 is widely known for its ability to oxidize >70 different, mostly compact, low molecular weight drugs and other xenobiotic compounds. In addition CYP2E1 oxidizes much larger C9–C20 fatty acids that can serve as endogenous signaling molecules. Previously structures of CYP2E1 with small molecules revealed a small, compact CYP2E1 active site, which would be insufficient to accommodate medium and long chain fatty acids without conformational changes in the protein. In the current work we have determined how CYP2E1 can accommodate a series of fatty acid analogs by cocrystallizing CYP2E1 with ω-imidazolyl-octanoic fatty acid, ω-imidazolyl-decanoic fatty acid, and ω-imidazolyl-dodecanoic fatty acid. In each structure direct coordination of the imidazole nitrogen to the heme iron mimics the position required for native fatty acid substrates to yield the ω-1 hydroxylated metabolites that predominate experimentally. In each case rotation of a single Phe²⁹⁸ side chain merges the active site with an adjacent void, significantly altering the active site size and topology to accommodate fatty acids. The binding of these fatty acid ligands is directly opposite the channel to the protein surface and the binding observed for fatty acids in the bacterial cytochrome P450 BM3 (CYP102A1) from Bacillus megaterium. Instead of the BM3-like binding mode in the CYP2E1 channel, these structures reveal interactions between the fatty acid carboxylates and several residues in the F, G, and B′ helices at successive distances from the active site.     

Comparative modeling and molecular docking of orphan human CYP4V2 protein with fatty acid substrates: Insights into substrate specificity

Cytochromes P450 (CYPs) are a super family of heme-containing enzymes well-known for their monooxgenase reaction. There are 57 CYP isoenzymes found in human which exhibit specific physiological functions. Thirteen members of this super family are classified as "orphan" CYP because of their unknown enzymatic functions. CYP4V2 is found to be a potential drug target for Bietti crystalline corneoretinal dystrophy (BCD). However, three-dimensional structure, the active site topology and substrate binding modes of CYP4V2 remain unclear. In this study, the three-dimensional model of CYP4V2 was constructed using the homology modeling method. Four possible fatty acid substrates namely, caprylic, lauric, myrisitc and palmitic acids were optimized and evaluated for drug likeness using Lipinski's rule of five. Further, these substrates were docked into active sites of CYP4V2 and several key residues responsible for substrate binding were identified. These findings will be helpful for the structure-based drug design and detailed characterization of the biological roles of CYP4V2.     

Inhibitory Potency of 8-Methoxypsoralen on Cytochrome P450 2A6 (CYP2A6) Allelic Variants CYP2A6*15, CYP2A6*16, CYP2A6*21 and CYP2A6*22: Differential Susceptibility Due to Different Sequence Locations of the Mutations

Human cytochrome P450 2A6 (CYP2A6) is a highly polymorphic isoform of CYP2A subfamily. Our previous kinetic study on four CYP2A6 allelic variants (CYP2A6*15, CYP2A6*16, CYP2A6*21 and CYP2A6*22) have unveiled the functional significance of sequence mutations in these variants on coumarin 7-hydroxylation activity. In the present study, we further explored the ability of a typical CYP2A6 inhibitor, 8-methoxypsoralen (8-MOP), in inhibition of these alleles and we hypothesized that translational mutations in these variants are likely to give impact on 8-MOP inhibitory potency. The CYP2A6 variant and the wild type proteins were subjected to 8-MOP inhibition to yield IC₅₀ values. In general, a similar trend of change in the IC₅₀ and K_m values was noted among the four mutants towards coumarin oxidation. With the exception of CYP2A6*16, differences in IC₅₀ values were highly significant which implied compromised interaction of the mutants with 8-MOP. Molecular models of CYP2A6 were subsequently constructed and ligand-docking experiments were performed to rationalize experimental data. Our docking study has shown that mutations have induced enlargement of the active site volume in all mutants with the exception of CYP2A6*16. Furthermore, loss of hydrogen bond between 8-MOP and active site residue Asn297 was evidenced in all mutants. Our data indicate that the structural changes elicited by the sequence mutations could affect 8-MOP binding to yield differential enzymatic activities in the mutant CYP2A6 proteins.