Investigation by site-directed mutagenesis of the role of cytochrome P450 2B4 non-active-site residues in protein-ligand interactions based on crystal structures of the ligand-bound enzyme

Recombinant Technologies Facilitate Drug Metabolism, Pharmacokinetics, and General Biomedical Research

The development of safe and effective medications requires a profound understanding of their pharmacokinetic (PK) and pharmacodynamic properties. PK studies have been built through investigation of enzymes and transporters that drive drug absorption, distribution, metabolism, and excretion (ADME). Like many other disciplines, the study of ADME gene products and their functions has been revolutionized through the invention and widespread adoption of recombinant DNA technologies. Recombinant DNA technologies use expression vectors such as plasmids to achieve heterologous expression of a desired transgene in a specified host organism. This has enabled the purification of recombinant ADME gene products for functional and structural characterization, allowing investigators to elucidate their roles in drug metabolism and disposition. This strategy has also been used to offer recombinant or bioengineered RNA (BioRNA) agents to investigate the posttranscriptional regulation of ADME genes. Conventional research with small noncoding RNAs such as microRNAs (miRNAs) and small interfering RNAs has been dependent on synthetic RNA analogs that are known to carry a range of chemical modifications expected to improve stability and PK properties. Indeed, a novel transfer RNA fused pre-miRNA carrier-based bioengineering platform technology has been established to offer consistent and high-yield production of unparalleled BioRNA molecules from Escherichia coli fermentation. These BioRNAs are produced and processed inside living cells to better recapitulate the properties of natural RNAs, representing superior research tools to investigate regulatory mechanisms behind ADME. SIGNIFICANCE STATEMENT: This review article summarizes recombinant DNA technologies that have been an incredible boon in the study of drug metabolism and PK, providing investigators with powerful tools to express nearly any ADME gene products for functional and structural studies. It further overviews novel recombinant RNA technologies and discusses the utilities of bioengineered RNA agents for the investigation of ADME gene regulation and general biomedical research.

Tyrosine Nitration Contributes to Nitric Oxide-Stimulated Degradation of CYP2B6

Human cytochrome P450 (P450) CYP2B6 undergoes nitric oxide (NO)-dependent proteasomal degradation in response to the NO donor dipropylenetriamine NONOate (DPTA) and biologic NO in HeLa and HuH7 cell lines. CYP2B6 is also downregulated by NO in primary human hepatocytes. We hypothesized that NO or derivative reactive nitrogen species may generate adducts of tyrosine and/or cysteine residues, causing CYP2B6 downregulation, and selected Tyr and Cys residues for mutation based on predicted solvent accessibility. CYP2B6V5-Y317A, -Y380A, and -Y190A mutant proteins expressed in HuH7 cells were less sensitive than wild-type (WT) enzyme to degradation evoked by DPTA, suggesting that these tyrosines are targets for NO-dependent downregulation. The Y317A or Y380A mutants did not show increases in high molecular mass (HMM) species after treatment with DPTA or bortezomib + DPTA, in contrast to the WT enzyme. Carbon monoxide-releasing molecule 2 treatment caused rapid suppression of 2B6 enzyme activity, significant HMM species generation, and ubiquitination of CYP2B6 protein but did not stimulate CYP2B6 degradation. The CYP2B6 inhibitor 4-(4-chlorophenyl)imidazole blocked NO-dependent CYP2B6 degradation, suggesting that NO access to the active site is important. Molecular dynamics simulations predicted that tyrosine nitrations of CYP2B6 would cause significant destabilizing perturbations of secondary structure and remove correlated motions likely required for enzyme function. We propose that cumulative nitrations of Y190, Y317, and Y380 by reactive nitrogen species cause destabilization of CYP2B6, which may act synergistically with heme nitrosylation to target the enzyme for degradation. SIGNIFICANCE STATEMENT: This work provides novel insight into the mechanisms by which nitric oxide, which is produced in hepatocytes in response to inflammation, triggers the ubiquitin-dependent proteasomal degradation of the cytochrome P450 (P450) enzyme CYP2B6. Our data demonstrate that both nitration of specific tyrosine residues and interaction of nitric oxide (NO) with the P450 heme are necessary for NO to trigger ubiquitination and protein degradation.

Dissecting the Structural Plasticity and Dynamics of Cytochrome P450 2B4 by Molecular Dynamics Simulations

The plasticity of cytochromes P450 (P450s) is known to contribute significantly to their catalytic capacity of metabolizing various substrates. Although numerous studies have been performed, factors governing the plasticity and dynamics of P450s are still not fully understood. In this study, taking CYP2B4 as an example, we dissect the protein plasticity and dynamics in different environments. CYP2B4 is featured by a high degree of plasticity, which exhibits open, closed, and intermediate states. By analyzing the CYP2B4 crystal structures, we identified the structural features for the closed, open, and intermediate states. Interestingly, formation of the dimer structure was found in the open and intermediate states. The subsequent molecular dynamics (MD) simulations of the open structure in water confirmed the importance of the dimer form in stabilizing the open conformations. MD simulations of the closed and open structures in the membrane environment and the free energies for opening the F-G cassette obtained from the umbrella sampling calculations indicate that the membrane environment is important for stabilizing the F-G cassette. The dynamical network analysis indicates that Asp105 on the B-C loop plays an important role in transiting the structure from the open to the intermediate state. Our results thus unveil the mechanisms of dimer formation and open-to-intermediate transition for CYP2B4 in the water and membrane environments.

Identification of a cytochrome P450 hydroxylase, CYP81E22, as a causative gene for the high sensitivity of soybean to herbicide bentazon

KEY MESSAGE

A frame shift invoked by a single-base deletion in the gene encoding a cytochrome P450 hydroxylase, CYP81E22, causes the loss of bentazon detoxification function in soybean. Bentazon is an effective herbicide in soybean cultivation applied at post-emergence stages for control of several broadleaf weeds. However, some soybean cultivars are highly sensitive to bentazon and are killed upon application. In this study, the gene related to the high sensitivity of soybean cultivars to bentazon was mapped to chromosome 16, and its location was narrowed down to a 257-kb region where three cytochrome P450 genes were located. In these genes, a single-base deletion of cytosine was detected in the coding region of Glyma.16G149300, CYP81E22, at + 1465 bp downstream from the translation start codon, leading to a frame shift in the open reading frame and creating a premature stop codon. This stop codon resulted in the loss of more than half of the P450, and consequently, the remaining molecule failed to form a functioning protein. This single-base deletion was common among the highly sensitive cultivars screened from the soybean mini-core collection and other previously reported highly sensitive cultivars. Furthermore, we screened plant lines from the targeting-induced local lesions in genomes library of the soybean cultivar Enrei based on a modelled 3D structure of CYP81E22. The lines with mutations in Glyma.16G149300 were highly sensitive to bentazon, which provides strong evidence that Glyma.16G149300 is the gene responsible for high sensitivity to bentazon.

Advances in the Understanding of Protein-Protein Interactions in Drug Metabolizing Enzymes through the Use of Biophysical Techniques

In recent years, a growing appreciation has developed for the importance of protein-protein interactions to modulate the function of drug metabolizing enzymes. Accompanied with this appreciation, new methods and technologies have been designed for analyzing protein-protein interactions both in vitro and in vivo. These technologies have been applied to several classes of drug metabolizing enzymes, including: cytochrome P450's (CYPs), monoamine oxidases (MAOs), UDP-glucuronosyltransferases (UGTs), glutathione S-transferases (GSTs), and sulfotransferases (SULTs). In this review, we offer a brief description and assessment of the impact of many of these technologies to the study of protein-protein interactions in drug disposition. The still expanding list of these techniques and assays has the potential to revolutionize our understanding of how these enzymes carry out their important functions in vivo.

Coumarin Derivatives as Substrate Probes of Mammalian Cytochromes P450 2B4 and 2B6: Assessing the Importance of 7-Alkoxy Chain Length, Halogen Substitution, and Non-Active Site Mutations

Using a combined structural and biochemical approach, the functional importance of a recently described peripheral pocket bounded by the E-, F-, G-, and I-helices in CYP2B4 and 2B6 was probed. Three series of 4-substituted-7-alkoxycoumarin derivatives with -H, -CH3, or -CF3 at the 4 position of the coumarin core were used initially to monitor functional differences between CYP2B4 and 2B6. 7-Ethoxy-4-(trifluoromethyl)coumarin (7-EFC) displayed the highest catalytic efficiency among these substrates. Mutants were made to alter side-chain polarity (V/E194Q) or bulk (F/Y244W) to alter access to the peripheral pocket. Modest increases in catalytic efficiency of 7-EFC O-deethylation by the mutants were magnified considerably by chlorination or bromination of the substrate ethoxy chain. A structure of CYP2B6 Y244W in complex with (+)-α-pinene was solved at 2.2 Å and showed no CYMAL-5 in the peripheral pocket. A ligand free structure of CYP2B4 F244W was solved at 3.0 Å with CYMAL-5 in the peripheral pocket. In both instances, comparison of the respective wild-type and mutant CYP2B enzymes revealed that CYMAL-5 occupancy of the peripheral pocket had little effect on the topology of active site residue side-chains, despite the fact that the peripheral pocket and active site are located on opposite sides of the I-helix. Analysis of available CYP2B structures suggest that the effect of the amino acid substitutions within the peripheral pocket derive from altered interactions between the F and G helices.

Structure-Function Analysis of Mammalian CYP2B Enzymes Using 7-Substituted Coumarin Derivatives as Probes: Utility of Crystal Structures and Molecular Modeling in Understanding Xenobiotic Metabolism

Crystal structures of CYP2B35 and CYP2B37 from the desert woodrat were solved in complex with 4-(4-chlorophenyl)imidazole (4-CPI). The closed conformation of CYP2B35 contained two molecules of 4-CPI within the active site, whereas the CYP2B37 structure demonstrated an open conformation with three 4-CPI molecules, one within the active site and the other two in the substrate access channel. To probe structure-function relationships of CYP2B35, CYP2B37, and the related CYP2B36, we tested the O-dealkylation of three series of related substrates-namely, 7-alkoxycoumarins, 7-alkoxy-4-(trifluoromethyl)coumarins, and 7-alkoxy-4-methylcoumarins-with a C1-C7 side chain. CYP2B35 showed the highest catalytic efficiency (kcat/KM) with 7-heptoxycoumarin as a substrate, followed by 7-hexoxycoumarin. In contrast, CYP2B37 showed the highest catalytic efficiency with 7-ethoxy-4-(trifluoromethyl)coumarin (7-EFC), followed by 7-methoxy-4-(trifluoromethyl)coumarin (7-MFC). CYP2B35 had no dealkylation activity with 7-MFC or 7-EFC. Furthermore, the new CYP2B-4-CPI-bound structures were used as templates for docking the 7-substituted coumarin derivatives, which revealed orientations consistent with the functional studies. In addition, the observation of multiple -Cl and -NH-π interactions of 4-CPI with the aromatic side chains in the CYP2B35 and CYP2B37 structures provides insight into the influence of such functional groups on CYP2B ligand binding affinity and specificity. To conclude, structural, computational, and functional analysis revealed striking differences between the active sites of CYP2B35 and CYP2B37 that will aid in the elucidation of new structure-activity relationships.

Functional importance of a peripheral pocket in mammalian cytochrome P450 2B enzymes

The functional importance of a peripheral pocket found in previously published X-ray crystal structures of CYP2B4 and CYP2B6 was probed using a biophysical approach. Introduction of tryptophan within the pocket of CYP2B4 at F202 or I241 leads to marked impairment of 7-ethoxy-4-(trifluoromethyl)coumarin (7-EFC) or 7-benzyloxyresorufin O-dealkylation efficiency; a similar substitution at F195, near the surface access to the pocket, does not affect these activities. The analogous CYP2B6 F202W mutant is inactive in the 7-EFC O-dealkylation assay. The stoichiometry of 7-EFC deethylation suggested that the decreased activity of F202W and I241W in CYP2B4 and lack of activity of F202W in CYP2B6 coincided with a sharp increase in the flux of reducing equivalents through the oxidase shunt to produce excess water. The results indicate that the chemical identity of residues within this peripheral pocket, but not at the mouth of the pocket, is important in substrate turnover and redox coupling, likely through effects on active site topology.

Substrate tunnels in enzymes: structure-function relationships and computational methodology

In enzymes, the active site is the location where incoming substrates are chemically converted to products. In some enzymes, this site is deeply buried within the core of the protein, and, in order to access the active site, substrates must pass through the body of the protein via a tunnel. In many systems, these tunnels act as filters and have been found to influence both substrate specificity and catalytic mechanism. Identifying and understanding how these tunnels exert such control has been of growing interest over the past several years because of implications in fields such as protein engineering and drug design. This growing interest has spurred the development of several computational methods to identify and analyze tunnels and how ligands migrate through these tunnels. The goal of this review is to outline how tunnels influence substrate specificity and catalytic efficiency in enzymes with buried active sites and to provide a brief summary of the computational tools used to identify and evaluate these tunnels.

A review of metabolic and enzymatic engineering strategies for designing and optimizing performance of microbial cell factories

Microbial cell factories (MCFs) are of considerable interest to convert low value renewable substrates to biofuels and high value chemicals. This review highlights the progress of computational models for the rational design of an MCF to produce a target bio-commodity. In particular, the rational design of an MCF involves: (i) product selection, (ii) de novo biosynthetic pathway identification (i.e., rational, heterologous, or artificial), (iii) MCF chassis selection, (iv) enzyme engineering of promiscuity to enable the formation of new products, and (v) metabolic engineering to ensure optimal use of the pathway by the MCF host. Computational tools such as (i) de novo biosynthetic pathway builders, (ii) docking, (iii) molecular dynamics (MD) and steered MD (SMD), and (iv) genome-scale metabolic flux modeling all play critical roles in the rational design of an MCF. Genome-scale metabolic flux models are of considerable use to the design process since they can reveal metabolic capabilities of MCF hosts. These can be used for host selection as well as optimizing precursors and cofactors of artificial de novo biosynthetic pathways. In addition, recent advances in genome-scale modeling have enabled the derivation of metabolic engineering strategies, which can be implemented using the genomic tools reviewed here as well.

Human cytochrome P450 17A1 conformational selection: modulation by ligand and cytochrome b5

Crystallographic studies of different membrane cytochrome P450 enzymes have provided examples of distinct structural conformations, suggesting protein flexibility. It has been speculated that conformational selection is an integral component of substrate recognition and access, but direct evidence of such substate interconversion has thus far remained elusive. In the current study, solution NMR revealed multiple and exchanging backbone conformations for certain structural features of the human steroidogenic cytochrome P450 17A1 (CYP17A1). This bifunctional enzyme is responsible for pregnenolone C17 hydroxylation, followed by a 17,20-lyase reaction to produce dehydroepiandrosterone, the key intermediate in human synthesis of androgen and estrogen sex steroids. The distribution of CYP17A1 conformational states was influenced by temperature, binding of these two substrates, and binding of the soluble domain of cytochrome b5 (b5). Notably, titration of b5 to CYP17A1·pregnenolone induced a set of conformational states closely resembling those of CYP17A1·17α-hydroxypregnenolone without b5, providing structural evidence consistent with the reported ability of b5 to selectively enhance 17,20-lyase activity. Solution NMR thus revealed a set of conformations likely to modulate human steroidogenesis by CYP17A1, demonstrating that this approach has the potential to make similar contributions to understanding the functions of other membrane P450 enzymes involved in drug metabolism and disease states.

Functional characterization of cytochromes P450 2B from the desert woodrat Neotoma lepida

Mammalian detoxification processes have been the focus of intense research, but little is known about how wild herbivores process plant secondary compounds, many of which have medicinal value or are drugs. cDNA sequences that code for three enzymes of the cytochrome P450 (CYP) 2B subfamily, here termed 2B35, 2B36, and 2B37 have been recently identified from a wild rodent, the desert woodrat (Malenke et al., 2012). Two variant clones of each enzyme were engineered to increase protein solubility and to facilitate purification, as reported for CYP2B enzymes from multiple species. When expressed in Escherichia coli each of the woodrat proteins gave the characteristic maximum at 450nm in a reduced carbon monoxide difference spectrum but generally expressed at lower levels than rat CYP2B1. Two enzymes, 2B36 and 2B37, showed dealkylation activity with the model substrates 7-ethoxy-4-(trifluoromethyl)coumarin and 7-benzyloxyresorufin, whereas 2B35 was inactive. Binding of the monoterpene (+)-α-pinene produced a Type I shift in the absorbance spectrum of each enzyme. Mutation of 2B37 at residues 114, 262, or 480, key residues governing ligand interactions with other CYP2B enzymes, did not significantly change expression levels or produce the expected functional changes. In summary, two catalytic and one ligand-binding assay are sufficient to distinguish among CYP2B35, 2B36, and 2B37. Differences in functional profiles between 2B36 and 2B37 are partially explained by changes in substrate recognition site residue 114, but not 480. The results advance our understanding of the mechanisms of detoxification in wild mammalian herbivores and highlight the complexity of this system.

A structural snapshot of CYP2B4 in complex with paroxetine provides insights into ligand binding and clusters of conformational states

An X-ray crystal structure of CYP2B4 in complex with the drug paroxetine [(3S,4R)-3-[(2H-1,3-benzodioxol-5-yloxy)methyl]-4-(4-fluorophenyl)piperidine] was solved at 2.14 Å resolution. The structure revealed a conformation intermediate to that of the recently solved complex with amlodipine and that of the more compact complex with 4-(4-chlorophenyl)imidazole in terms of the placement of the F-G cassette. Moreover, comparison of the new structure with 15 previously solved structures of CYP2B4 revealed some new insights into the determinants of active-site size and shape. The 2B4-paroxetine structure is nearly superimposable on a previously solved closed structure in a ligand-free state. Despite the overall conformational similarity among multiple closed structures, the active-site cavity volume of the paroxetine complex is enlarged. Further analysis of the accessible space and binding pocket near the heme reveals a new subchamber that resulted from the movement of secondary structural elements and rearrangements of active-site side chains. Overall, the results from the comparison of all 16 structures of CYP2B4 demonstrate a cluster of protein conformations that were observed in the presence or absence of various ligands.

The key role of a non-active-site residue Met148 on the catalytic efficiency of meta-cleavage product hydrolase BphD

meta-Cleavage product (MCP) hydrolases (EC 3.7.1.9) can catalyze a specific C-C bond fission during the microbial aerobic degradation of aromatics. The previous studies on structure-function relationship of MCP hydrolases mainly focus on the active site residues by site-directed mutagenesis. However, the information about the role of the non-active-site residues is still unclear. In this study, a non-active-site residue Met148 of MCP hydrolase BphD was selected as the mutagenesis site according to the sequence alignments, structure superimpose and the tunnel analysis, which underwent the saturation mutagenesis resulting 19 mutants. The catalytic efficiencies of the mutants on 6-oxo-6-phenylhexa-2,4-dienoic acid (HOPDA) were all decreased compared with the wild-type one except for the M148D mutant. Especially, the M148P mutant exhibited 290-fold lower k cat/K m than that of the wild-type BphD. Transient kinetic analyses of M148P showed the reciprocal relaxation time corresponded to C-C bond cleavage and product release steps (9.6 s(-1)) was 4.08-fold lower than BphD WT (39.2 s(-1)). Tunnel cluster analysis of BphD WT, M148P and M148W demonstrated that only the bulky Trp148 could block tunnel T2 in the BphD WT, but it exhibited slight effects on the catalytic efficiency (0.94-fold of BphD WT). Therefore, product release was not the main reason for the efficiency decrease of M148P. On the other hand, molecular dynamics simulations on the BphD WT and BphD M148P in complex with HOPDA indicated that the dramatic decrease of the catalytic efficiencies of BphD M148P should be due to the unproductive binding of HOPDA. The study demonstrated the catalytic efficiency of MCP hydrolase can be engineered by modification of non-active site residue.

Steroid-based facial amphiphiles for stabilization and crystallization of membrane proteins

Amphiphile selection is a critical step for structural studies of membrane proteins (MPs). We have developed a family of steroid-based facial amphiphiles (FAs) that are structurally distinct from conventional detergents and previously developed FAs. The unique FAs stabilize MPs and form relatively small protein-detergent complexes (PDCs), a property considered favorable for MP crystallization. We attempted to crystallize several MPs belonging to different protein families, including the human gap junction channel protein connexin 26, the ATP binding cassette transporter MsbA, the seven-transmembrane G protein-coupled receptor-like bacteriorhodopsin, and cytochrome P450s (peripheral MPs). Using FAs alone or mixed with other detergents or lipids, we obtained 3D crystals of the above proteins suitable for X-ray crystallographic analysis. The fact that FAs enhance MP crystallizability compared with traditional detergents can be attributed to several properties, including increased protein stability, formation of small PDCs, decreased PDC surface flexibility, and potential to mediate crystal lattice contacts.

Distal effect of amino acid substitutions in CYP2C9 polymorphic variants causes differences in interatomic interactions against (S)-warfarin

Cytochrome P450 2C9 (CYP2C9) is crucial in excretion of commonly prescribed drugs. However, changes in metabolic activity caused by CYP2C9 polymorphisms inevitably result in adverse drug effects. CYP2C9*2 and *3 are prevalent in Caucasian populations whereas CYP2C9*13 is remarkable in Asian populations. Single amino acid substitutions caused by these mutations are located outside catalytic cavity but affect kinetic activities of mutants compared to wild-type enzyme. To relate distal effects of these mutations and defective drug metabolisms, simulations of CYP2C9 binding to anti-coagulant (S)-warfarin were performed as a system model. Representative (S)-warfarin-bound forms of wild-type and mutants were sorted and assessed through knowledge-based scoring function. Interatomic interactions towards (S)-warfarin were predicted to be less favorable in mutant structures in correlation with larger distance between hydroxylation site of (S)-warfarin and reactive oxyferryl heme than wild-type structure. Using computational approach could delineate complication of CYP polymorphism in management of drug therapy.

Conformational adaptation of human cytochrome P450 2B6 and rabbit cytochrome P450 2B4 revealed upon binding multiple amlodipine molecules

Structures of human cytochrome P450 2B6 and rabbit cytochrome P450 2B4 in complex with two molecules of the calcium channel blocker amlodipine have been determined by X-ray crystallography. The presence of two drug molecules suggests clear substrate access channels in each P450. According to a previously established nomenclature, amlodipine molecules were trapped in access pathway 2f in P450 2B6 and in pathway 2a or 2f in P450 2B4. These pathways overlap for part of the length and then diverge as they extend toward the protein surface. A previously described solvent channel was also found in each enzyme. The results indicate that key residues located on the surface and at the entrance of the substrate access channels in each of these P450s may play a crucial role in guiding substrate entry. In addition, the region of P450 2B6 and 2B4 involving helices B', F, F', and G' and part of helix G is substantially more open in the amlodipine complexes than in the corresponding 4-(4-chlorophenyl)imidazole complexes. The increased active site volume observed results from the major retraction of helices F, F', and B' and the β4 sheet region located close to the binding cavity to accommodate amlodipine. These structures demonstrate novel insight into distinct conformational states not observed with previous P450 2B structures and provide clear evidence of the substrate access channels in two drug-metabolizing P450s. In addition, the structures exhibit the versatility that can be exploited via in silico studies with other P450 2B6 ligands as large as raloxifene and itraconazole.

Bioinformatics and variability in drug response: a protein structural perspective

Marketed drugs frequently perform worse in clinical practice than in the clinical trials on which their approval is based. Many therapeutic compounds are ineffective for a large subpopulation of patients to whom they are prescribed; worse, a significant fraction of patients experience adverse effects more severe than anticipated. The unacceptable risk-benefit profile for many drugs mandates a paradigm shift towards personalized medicine. However, prior to adoption of patient-specific approaches, it is useful to understand the molecular details underlying variable drug response among diverse patient populations. Over the past decade, progress in structural genomics led to an explosion of available three-dimensional structures of drug target proteins while efforts in pharmacogenetics offered insights into polymorphisms correlated with differential therapeutic outcomes. Together these advances provide the opportunity to examine how altered protein structures arising from genetic differences affect protein-drug interactions and, ultimately, drug response. In this review, we first summarize structural characteristics of protein targets and common mechanisms of drug interactions. Next, we describe the impact of coding mutations on protein structures and drug response. Finally, we highlight tools for analysing protein structures and protein-drug interactions and discuss their application for understanding altered drug responses associated with protein structural variants.

A novel type of allosteric regulation: functional cooperativity in monomeric proteins

Cooperative functional properties and allosteric regulation in cytochromes P450 play an important role in xenobiotic metabolism and define one of the main mechanisms of drug-drug interactions. Recent experimental results suggest that ability to bind simultaneously two or more small organic molecules can be the essential feature of cytochrome P450 fold, and often results in rich and complex pattern of allosteric behavior. Manifestations of non-Michaelis kinetics include homotropic and heterotropic activation and inhibition effects depending on the stoichiometric ratios of substrate and effector, changes in the regio- and stereospecificity of catalytic transformations, and often give rise to the clinically important drug-drug interactions. In addition, functional response of P450 systems is modulated by the presence of specific and non-specific effector molecules, metal ions, membrane incorporation, formation of homo- and hetero-oligomers, and interactions with the protein redox partners. In this article we briefly overview the main factors contributing to the allosteric effects in cytochromes P450 with the main focus on the sources of cooperative behavior in xenobiotic metabolizing monomeric heme enzymes with their conformational flexibility and extremely broad substrate specificity. The novel mechanism of functional cooperativity in P450 enzymes does not require substantial binding cooperativity, rather it implies the presence of one or more binding sites with higher affinity than the single catalytically active site in the vicinity of the heme iron.