Conformational Mobility in Cytochrome P450 3A4 Explored by Pressure-Perturbation EPR Spectroscopy

An efficient EPR spin-labeling method enables insights into conformational changes in DNA

Electron paramagnetic resonance (EPR) is a powerful tool for elucidating both static and dynamic conformational alterations in macromolecules. However, to effectively utilize EPR for such investigations, the presence of paramagnetic centers, known as spin labels, is required. The process of spin labeling, particularly for nucleotides, typically demands intricate organic synthesis techniques. In this study, we introduce a unique addition-elimination reaction method with a simple spin-labeling process, facilitating the monitoring of structural changes within nucleotide sequences. Our investigation focuses on three distinct labeling positions with a DNA sequence, allowing the measurement of distance between two spin labels. The experimental mean distances obtained agreed with the calculated distances, underscoring the efficacy of this straightforward spin-labeling approach in studying complex biological processes such as transcription mechanism using EPR measurements.

Ni2+ and Cu2+ complexes of N-(2,6-dichlorophenyl)-N-mesityl formamidine dithiocarbamate structural and functional properties as CYP3A4 potential substrates

Metal compounds continued to attract diverse applications due to their malleability in several capacities. In this study, we present our findings on the crystal structures and functional properties of Ni2+ and Cu2+ complexes of N'-(2,6-dichlorophenyl)-N-mesitylformamidine dithiocarbamate (L) comprising [Ni-(L)2] (1) and [Cu-(L)2] (2) with a four-coordinate metal center. We established the two complex structures through 1H and 13C nuclear magnetic resonance (NMR), elemental, and single-crystal X-ray analysis. The analyses showed that the two complexes are isomorphous, having P21/c as a space group and a unit-cell similarity index (π) of 0.002. The two complexes conform to a distorted square planar geometry around the metal centers. The calculated and experimental data, including bond lengths, angles, and NMR values, are similar. Hirshfeld surface analysis revealed the variational contribution of the different types of intermolecular contacts driven by the crystal lattice of the two solvated complexes. Our knowledge of the potential biological implication of these structures enabled us to probe the compounds as prospective CYP3A4 inhibitors. This approach mimics current trends in pharmaceutical design and biomedicine by incorporating potentially active molecules into various media to predict their biological efficacies. The simulations show appreciable binding of compounds 1 and 2 to CYP3A4 with average interaction energies of -97 and -87 kcal/mol, respectively. The protein attains at least five conformational states in the three studied models using a Gaussian Mixture Model-based clustering and free energy prediction. Electric field analysis shows the crucial residues to substrate binding at the active site, enabling CYP3A4 structure to function prediction. The predicted inhibition with these Ni2+ and Cu2+ complexes indicates that CYP3A4 overexpression in a diseased state like cancer would reduce, thereby increasing the chemotherapeutic compounds' shelf-lives for adsorption. This multidimensional study addresses various aspects of molecular metal electronics, including their application as substrate-mimicking inhibitors. The outcome would enable further research on bio-metal compounds of critical potential.

Nanodisc-embedded cytochrome P450 P3A4 binds diverse ligands by distributing conformational dynamics to its flexible elements

Cytochrome P450 3A4 (CYP3A4) metabolizes a wide range of drugs and toxins. Interactions of CYP3A4 with ligands are difficult to predict due to promiscuity and conformational flexibility. To better understand CYP3A4 conformational responses to ligands we use hydrogen deuterium exchange mass spectrometry (HDX-MS) to investigate the effect of ligands on nanodisc-embedded CYP3A4. For a subset of CYP3A4-ligand complexes, differences in the low-frequency modes derived by principal component analyses of molecular dynamics trajectories mirrored the HDX-MS results. The effects of ligands are distributed to flexible elements of CYP3A4 between stretches of secondary structure. The largest effects occur in the F- and G-helices, where most ligands increase the flexibility of the F-helix and connecting loops and decrease the flexibility of the C-term of the G-helix. Most ligands affect the E-F-G, CD and HI regions of the protein. Ligand-dependent differences are observed in the A"-A' loop, BC region, E-helix, K-β1 region, proximal loop, and C-term loop. Correlated HDX responses were observed in the CD region and the C-term of the G-helix that were most pronounced for Type II ligands. Collectively, the HDX and molecular dynamics results suggest that CYP3A4 accommodates diverse binding partners by propagating local backbone fluctuations from the binding site onto the flexible regions of the enzyme via long-range interactions that are differentially modulated by ligands. In contrast to the paradigm wherein ligands decrease protein dynamics at their binding site, a wide range of ligands modestly increase CYP3A4 dynamics throughout the protein including effects remote from the active site.

Allostery‐driven changes in dynamics regulate the activation of bacterial copper transcription factor

Metalloregulators bind and respond to metal ions by regulating the transcription of metal homeostasis genes. Copper efflux regulator (CueR) is a copper‐responsive metalloregulator that is found in numerous Gram‐negative bacteria. Upon Cu(I) coordination, CueR initiates transcription by bending the bound DNA promoter regions facilitating interaction with RNA polymerase. The structure of Escherichia coli CueR in presence of DNA and metal ion has been reported using X‐ray crystallography and cryo‐EM, providing information about the mechanism of action. However, the specific role of copper in controlling this transcription mechanism remains elusive. Herein, we use room temperature electron paramagnetic resonance (EPR) experiments to follow allosterically driven dynamical changes in E. coli CueR induced by Cu(I) binding. We suggest that more than one Cu(I) ion binds per CueR monomer, leading to changes in site‐specific dynamics at the Cu(I) binding domain and at the distant DNA binding site. Interestingly, Cu(I) binding leads to an increase in dynamics about 27 Å away at the DNA binding domain. These changes in the dynamics of the DNA binding domain are important for exact coordination with the DNA. Thus, Cu(I) binding is critical to initiate a series of conformational changes that regulate and initiate gene transcription.

High-Pressure ESR Spectroscopy: On the Rotational Motion of Spin Probes in Pressurized Ionic Liquids

We report high-pressure (up to 50 MPa) ESR-spectroscopic investigations on the rotational correlation times of the nitroxide radicals 2,2,6,6-tetramethylpiperidine 1-oxyl (TEMPO), 4-hydroxy-2,2,6,6-tetramethylpiperidine 1-oxyl (TEMPOL), and 4-amino-2,2,6,6-tetramethylpiperidine 1-oxyl (ATEMPO) in the ionic liquids 1-ethyl-3-methylimidazolium tetrafluoroborate (emimBF₄), 1-butyl-3-methylimidazolium hexafluorophosphate (bmimPF₆), 1-butyl-3-methylimidazolium tetrafluoroborate (bmimBF₄), 1-methyl-3-octylimidazolium tetrafluoroborate (omimBF₄), and 1-methyl-3-octylimidazolium hexafluorophosphate (omimPF₆). The activation volumes (38.5–56.6 ų) determined from pressure dependent rotational diffusion coefficients agree well with the pressure dependent viscosities of the ionic liquids. Experimentally, the fractional exponent of the generalized Stokes–Einstein–Debye relation is found to be close to one.

Probing functional interactions between cytochromes P450 with principal component analysis of substrate saturation profiles and targeted proteomics

We investigated the correspondence between drug metabolism routes and the composition of the P450 ensemble in human liver microsomes (HLM). As a probe, we used Coumarin 152 (C152), a fluorogenic substrate metabolized by multiple P450 species. Studying the substrate-saturation profiles (SSP) in seven pooled HLM preparations, we sought to correlate them with the P450 pool's composition characterized by targeted proteomics. This analysis, complemented with the assays with specific inhibitors of CYP3A4 and CYP2C19, the primary C152 metabolizers, demonstrated a significant contrast between different HLM samples. To unveil the source of these differences, we implemented Principal Component Analysis (PCA) of the SSP series obtained with HLM samples with a known composition of the P450 pool. Our analysis revealed that the parameters of C152 metabolism are primarily determined by the content of CYP2A6, CYP2B6, CYP2C8, CYP2E1, and CYP3A5 of those only CYP2B6 and CYP3A5 can metabolize C152. To validate this finding, we studied the effect of enriching HLM with CYP2A6, CYP2E1, and CYP3A5. The incorporation of CYP3A5 into HLM decreases the rate of C152 metabolism while increasing the role of CYP2B6 in its turnover. In contrast, incorporation of CYP2A6 and CYP2E1 reroutes the C152 demethylation towards some P450 enzyme with a moderate affinity to the substrate, most likely CYP3A4. Our results reveal a sharp non-additivity of the individual P450 properties and suggest a pivotal role of P450-P450 interactions in determining drug metabolism routes. This study demonstrates the high potential of our new PCA-based approach in unveiling functional interrelationships between different P450 species.

Insights into oral bioavailability enhancement of therapeutic herbal constituents by cytochrome P450 3A inhibition

Herbal plants typically have complex compositions and diverse mechanisms. Among them, bioactive constituents with relatively high exposure in vivo are likely to exhibit therapeutic efficacy. On the other hand, their bioavailability may be influenced by the synergistic effects of different bioactive components. Cytochrome P450 3A (CYP3A) is one of the most abundant CYP enzymes, responsible for the metabolism of 50% of approved drugs. In recent years, many therapeutic herbal constituents have been identified as CYP3A substrates. It is more evident that CYP3A inhibition derived from the herbal formula plays a critical role in improving the oral bioavailability of therapeutic constituents. CYP3A inhibition may be the mechanism of the synergism of herbal formula. In this review, we explored the multiplicity of CYP3A, summarized herbal monomers with CYP3A inhibitory effects, and evaluated herb-mediated CYP3A inhibition, thereby providing new insights into the mechanisms of CYP3A inhibition-mediated oral herb bioavailability.

Effects of alcohol-induced increase in CYP2E1 content in human liver microsomes on the activity and cooperativity of CYP3A4

We investigate the effect of the alcohol-induced increase in the content of CYP2E1 in human liver microsomes (HLM) on the function of CYP3A4. Membrane incorporation of the purified CYP2E1 into HLM considerably increases the rate of metabolism of 7-benzyloxyquinoline (BQ) and attenuates the homotropic cooperativity observed with this CYP3A4-specific substrate. It also eliminates the activating effect of α-naphthoflavone (ANF) seen in some HLM samples. To probe the physiological relevance of these effects, we compared three pooled preparations of HLM from normal donors (HLM-N) with a pooled preparation from ten heavy alcohol consumers (HLM-A). The composition of the P450 pool in all samples was characterized by the mass-spectrometric determination of 11 cytochrome P450 species. The fractional content of CYP2E1 in HLM-A was from 2.0 to 3.4 times higher than in HLM-N. In contrast, the content of CYP3A4 in HLM-A was the lowest among all samples. Despite that, HLM-A exhibited a much higher metabolism rate and a lower homotropic cooperativity with BQ, similar to CYP2E1-enriched HLM-N. To substantiate the involvement of interactions between CYP2E1 and CYP3A4 in these effects, we probed hetero-association of these proteins in CYP3A4-containing Supersomes™ with a technique employing CYP2E1 labeled with BODIPY-618 maleimide. These experiments evinced the interactions between the two enzymes and revealed an inhibitory effect of ANF on their association. Our results demonstrate that the functional properties of CYP3A4 are fundamentally dependent on the composition of the cytochrome P450 ensemble and suggest a possible impact of chronic alcohol exposure on the pharmacokinetics of drugs metabolized by CYP3A4.

Pressure tolerance of deep-sea enzymes can be evolved through increasing volume changes in protein transitions: a study with lactate dehydrogenases from abyssal and hadal fishes

We explore the principles of pressure tolerance in enzymes of deep-sea fishes using lactate dehydrogenases (LDH) as a case study. We compared the effects of pressure on the activities of LDH from hadal snailfishes Notoliparis kermadecensis and Pseudoliparis swirei with those from a shallow-adapted Liparis florae and an abyssal grenadier Coryphaenoides armatus. We then quantified the LDH content in muscle homogenates using mass-spectrometric determination of the LDH-specific conserved peptide LNLVQR. Existing theory suggests that adaptation to high pressure requires a decrease in volume changes in enzymatic catalysis. Accordingly, evolved pressure tolerance must be accompanied with an important reduction in the volume change associated with pressure-promoted alteration of enzymatic activity ( Δ V PP ∘ ). Our results suggest an important revision to this paradigm. Here, we describe an opposite effect of pressure adaptation-a substantial increase in the absolute value of Δ V PP ∘ in deep-living species compared to shallow-water counterparts. With this change, the enzyme activities in abyssal and hadal species do not substantially decrease their activity with pressure increasing up to 1-2 kbar, well beyond full-ocean depth pressures. In contrast, the activity of the enzyme from the tidepool snailfish, L. florae, decreases nearly linearly from 1 to 2500 bar. The increased tolerance of LDH activity to pressure comes at the expense of decreased catalytic efficiency, which is compensated with increased enzyme contents in high-pressure-adapted species. The newly discovered strategy is presumably used when the enzyme mechanism involves the formation of potentially unstable excited transient states associated with substantial changes in enzyme-solvent interactions.

Nonadditivity in human microsomal drug metabolism revealed in a study with coumarin 152, a polyspecific cytochrome P450 substrate

We closely characterized 7-Dimethylamino-4-trifluromethylcoumarin (Coumarin 152, C152), a substrate metabolized by multiple P450 species, to establish a new fluorogenic probe for the studies of functional integration in the cytochrome P450 ensemble. Scanning fluorescence spectroscopy and LC/MS-MS were used to characterize the products of N-demethylation of C152 and optimize their fluorometric detection. The metabolism of C152 by the individual P450 species was characterized using the microsomes containing cDNA-expressed enzymes. C152 metabolism in human liver microsomes (HLM) was studied in a preparation with quantified content of eleven P450 species. C152 is metabolized by CYP2B6, CYP3A4, CYP3A5, CYP2C19, CYP1A2, CYP2C9, and CYP2C8 listed in the order of decreasing turnover. The affinities exhibited by CYP3A5, CYP2C9, and CYP2C8 were lower than those characteristic to the other enzymes. The presumption of additivity suggests the participation of CYP3A4, CYP2B6, and CYP2C19 to be 84, 8, and 0.2%, respectively. Contrary to this prediction, inhibitory analysis identified CYP2C19 as the principal C152-metabolizing enzyme. We thoroughly characterize C152 for the studies of drug metabolism in HLM and demonstrate the limitations of the proportional projection approach by providing an example, where the involvement of individual P450 species cannot be predicted from their content.

Structure and Function of the Cytochrome P450 Monooxygenase Cinnamate 4-hydroxylase from Sorghum bicolor

Cinnamate 4-hydroxylase (C4H; CYP73A) is a cytochrome P450 monooxygenase associated externally with the endoplasmic reticulum of plant cells. The enzyme uses NADPH-cytochrome P450 reductase as a donor of electrons and hydroxylates cinnamic acid to form 4-coumaric acid in phenylpropanoid metabolism. In order to better understand the structure and function of this unique class of plant P450 enzymes, we have characterized the enzyme C4H1 from lignifying tissues of sorghum (Sorghum bicolor), encoded by Sobic.002G126600 Here we report the 1.7 Å resolution crystal structure of CYP73A33. The obtained structural information, along with the results of the steady-state kinetic analysis and the absorption spectroscopy titration, displays a high degree of similarity of the structural and functional features of C4H to those of other P450 proteins. Our data also suggest the presence of a putative allosteric substrate-binding site in a hydrophobic pocket on the enzyme surface. In addition, comparing the newly resolved structure with those of well-investigated cytochromes P450 from mammals and bacteria enabled us to identify those residues of critical functional importance and revealed a unique sequence signature that is potentially responsible for substrate specificity and catalytic selectivity of C4H.

Structural and dynamic origins of ESR lineshapes in spin-labeled GB1 domain: the insights from spin dynamics simulations based on long MD trajectories

Site-directed spin labeling (SDSL) ESR is a valuable tool to probe protein systems that are not amenable to characterization by x-ray crystallography, NMR or EM. While general principles that govern the shape of SDSL ESR spectra are known, its precise relationship with protein structure and dynamics is still not fully understood. To address this problem, we designed seven variants of GB1 domain bearing R1 spin label and recorded the corresponding MD trajectories (combined length 180 μs). The MD data were subsequently used to calculate time evolution of the relevant spin density matrix and thus predict the ESR spectra. The simulated spectra proved to be in good agreement with the experiment. Further analysis confirmed that the spectral shape primarily reflects the degree of steric confinement of the R1 tag and, for the well-folded protein such as GB1, offers little information on local backbone dynamics. The rotameric preferences of R1 side chain are determined by the type of the secondary structure at the attachment site. The rotameric jumps involving dihedral angles χ1 and χ2 are sufficiently fast to directly influence the ESR lineshapes. However, the jumps involving multiple dihedral angles tend to occur in (anti)correlated manner, causing smaller-than-expected movements of the R1 proxyl ring. Of interest, ESR spectra of GB1 domain with solvent-exposed spin label can be accurately reproduced by means of Redfield theory. In particular, the asymmetric character of the spectra is attributable to Redfield-type cross-correlations. We envisage that the current MD-based, experimentally validated approach should lead to a more definitive, accurate picture of SDSL ESR experiments.

CYP3A4 and microRNA-122 are involved in the apoptosis of HepG2 cells induced by ILs 1-decyl-3-methylimidazolium bromide

Ionic liquids (ILs) as green alternatives for volatile organic solvents are increasingly used in commercial applications. It is necessary to explore the cytotoxic mechanism of ILs to reduce the risk to human health. For this purpose, cell viability, apoptosis, cytochrome P450 3A4 (CYP3A4), glucose transporter type 2 (GLUT2), and microRNA-122 (miR-122) gene expression in HepG2 cells was evaluated after IL exposure. The results showed that ILs reduced the viability of HepG2 cells through apoptotic cell death. Moreover, ILs markedly upregulated the transcription and protein levels of CYP3A4, but did not affect the expression of GLUT2 in either messenger RNA level or protein level. Finally, ILs increased the expression of miR-122 and inhibition of miR-122 with miR-122 inhibitor blocked ILs-induced apoptosis in HepG2 cells. This finding may contribute to an increased understanding of the in vitro molecular toxicity mechanism of ILs to further understand IL-related human health risks.

Neonatal cytochrome P450 CYP3A7: A comprehensive review of its role in development, disease, and xenobiotic metabolism

The human cytochrome P450 CYP3A7, once thought to be an enzyme exclusive to fetal livers, has more recently been identified in neonates and developing infants as old as 24 months post-gestational age. CYP3A7 has been demonstrated to metabolize two endogenous compounds that are known to be important in the growth and development of the fetus and neonate, namely dehydroepiandrosterone sulfate (DHEA-S) and all-trans retinoic acid (atRA). In addition, it is also known to metabolize a variety of drugs and xenobiotics, albeit generally to a lesser extent relative to CYP3A4/5. CYP3A7 is an important component in the development and protection of the fetal liver and additionally plays a role in certain disease states, such as cancer and adrenal hyperplasia. Ultimately, a full understanding of the expression, regulation, and metabolic properties of CYP3A7 is needed to provide neonates with appropriate individualized pharmacotherapy. This article summarizes the current state of knowledge of CYP3A7, including its discovery, distribution, alleles, RNA splicing, expression and regulation, metabolic properties, substrates, and inhibitors.

Effects of polymorphic variation on the thermostability of heterogenous populations of CYP3A4 and CYP2C9 enzymes in solution

The effect of non-synonymous single nucleotide polymorphisms (SNPs) on cytochrome P450 (CYP450) drug metabolism is currently poorly understood due to the large number of polymorphisms, the diversity of potential substrates and the complexity of CYP450 function. Previously we carried out in silico studies to explore the effect of SNPs on CYP450 function, using in silico calculations to predict the effect of mutations on protein stability. Here we have determined the effect of eight CYP3A4 and seven CYP2C9 SNPs on the thermostability of proteins in solution to test these predictions. Thermostability assays revealed distinct CYP450 sub-populations with only 65–70% of wild-type CYP3A4 and CYP2C9 susceptible to rapid heat-induced P450 to P420 conversion. CYP3A4 mutations G56D, P218R, S222P, I223R, L373F and M445T and CYP2C9 mutations V76M, I359L and I359T were destabilising, increasing the proportion of protein sensitive to the rapid heat-induced P450 to P420 conversion and/or reducing the half-life of this conversion. CYP2C9 Q214L was the only stabilising mutation. These results corresponded well with the in silico protein stability calculations, confirming the value of these predictions and together suggest that the changes in thermostability result from destabilisation/stabilisation of the protein fold, changes in the haem-binding environment or effects on oligomer formation/conformation.

Steroid bioconjugation to a CYP3A4 allosteric site and its effect on substrate binding and coupling efficiency

Human cytochrome P450 3A4 (CYP3A4) is an important drug metabolizing enzyme involved in a number of drug-drug and food-drug interactions. As such, much effort has been devoted into investigating its mechanism of interaction with ligands. CYP3A4 has one of the highest levels of substrate promiscuity for an enzyme, and can even bind multiple ligands simultaneously. The location and orientation of these ligands depend on the chemical structure and stoichiometry, and are generally poorly understood. In the case of the steroid testosterone, up to three copies of the molecule can associate with the enzyme at once, likely two in the active site and one at a postulated allosteric site. Recently, we demonstrated that steroid bioconjugation at the allosteric site results in an increase in activity of CYP3A4 toward testosterone and 7-benzyloxy-4-trifluoromethylcoumarin oxidation. Here, using the established bioconjugation methodology, we show how steroid bioconjugation at the allosteric site affects the heme spin state, the binding affinity (K_S) of CYP3A4 for testosterone, as well as the enzyme coupling efficiency.

High-pressure EPR spectroscopy studies of the E. coli lipopolysaccharide transport proteins LptA and LptC

The use of pressure is an advantageous approach to the study of protein structure and dynamics because it can shift the equilibrium populations of protein conformations toward higher energy states that are not of sufficient population to be observable at atmospheric pressure. Recently, the Hubbell group at the University of California, Los Angeles, reintroduced the application of high pressure to the study of proteins by electron paramagnetic resonance (EPR) spectroscopy. This methodology is possible using X-band EPR spectroscopy due to advances in pressure intensifiers, sample cells, and resonators. In addition to the commercial availability of the pressure generation and sample cells by Pressure Biosciences Inc., a five-loop-four-gap resonator required for the initial high pressure EPR spectroscopy experiments by the Hubbell group, and those reported here, was designed by James S. Hyde and built and modified at the National Biomedical EPR Center. With these technological advances, we determined the effect of pressure on the essential periplasmic lipopolysaccharide (LPS) transport protein from Escherichia coli, LptA, and one of its binding partners, LptC. LptA unfolds from the N-terminus to the C-terminus, binding of LPS does not appreciably stabilize the protein under pressure, and monomeric LptA unfolds somewhat more readily than oligomeric LptA upon pressurization to 2 kbar. LptC exhibits a fold and relative lack of stability upon LPS binding similar to LptA, yet adopts an altered, likely monomeric, folded conformation under pressure with only its C-terminus unraveling. The pressure-induced changes likely correlate with functional changes associated with binding and transport of LPS.

High-Level Production and Properties of the Cysteine-Depleted Cytochrome P450 3A4

Human drug-metabolizing cytochrome P450 3A4 (CYP3A4) is a dynamic enzyme with a large and highly malleable active site that can fit structurally diverse compounds. Despite extensive investigations, structure-function relationships and conformational dynamics in CYP3A4 are not fully understood. This study was undertaken to engineer a well-expressed and functionally active cysteine-depleted CYP3A4 that can be used in biochemical and biophysical studies. cDNA codon optimization and screening mutagenesis were utilized to boost the level of bacterial expression of CYP3A4 and identify the least harmful substitutions for all six non-heme-ligating cysteines. The C58A/C64M/C98A/C239T/C377A/C468S (Cys-less) mutant was found to be expressed as highly as the optimized wild-type (opt-WT) CYP3A4. The high-resolution X-ray structures of opt-WT and Cys-less CYP3A4 revealed that gene optimization leads to a different folding in the Phe108 and Phe189 regions and promotes binding of the active site glycerol that interlocks Ser119 and Arg212, critical for ligand association, and the hydrophobic cluster adjacent to Phe108. Crowding and decreased flexibility of the active site, as well as structural alterations observed at the C64M, C239T, and C468S mutational sites, might be responsible for the distinct ligand binding behavior of opt-WT and Cys-less CYP3A4. Nonetheless, the Cys-less mutant could be used for structure-function investigations because it orients bromoergocryptine and ritonavir (a high-affinity substrate and a high-potency inhibitor, respectively) like the WT and has a higher activity toward 7-benzyloxy(4-trifluoromethyl)coumarin.

Physical Studies of P450-P450 Interactions: Predicting Quaternary Structures of P450 Complexes in Membranes from Their X-ray Crystal Structures

Cytochrome P450 enzymes, which catalyze oxygenation reactions of both exogenous and endogenous chemicals, are membrane bound proteins that require interaction with their redox partners in order to function. Those responsible for drug and foreign compound metabolism are localized primarily in the endoplasmic reticulum of liver, lung, intestine, and other tissues. More recently, the potential for P450 enzymes to exist as supramolecular complexes has been shown by the demonstration of both homomeric and heteromeric complexes. The P450 units in these complexes are heterogeneous with respect to their distribution and function, and the interaction of different P450s can influence P450-specific metabolism. The goal of this review is to examine the evidence supporting the existence of physical complexes among P450 enzymes. Additionally, the review examines the crystal lattices of different P450 enzymes derived from X-ray diffraction data to make assumptions regarding possible quaternary structures in membranes and in turn, to predict how the quaternary structures could influence metabolism and explain the functional effects of specific P450-P450 interactions.