Characterization of Interactions Among CYP1A2, CYP2B4, and NADPH-cytochrome P450 Reductase: Identification of Specific Protein Complexes

Identification of the N-terminal residues responsible for the differential microdomain localization of CYP1A1 and CYP1A2

The endoplasmic reticulum (ER) is organized into ordered regions enriched in cholesterol and sphingomyelin, and disordered microdomains characterized by more fluidity. Rabbit CYP1A1 and CYP1A2 localize into disordered and ordered microdomains, respectively. Previously, a CYP1A2 chimera containing the first 109 amino acids of CYP1A1 showed altered microdomain localization. The goal of this study was to identify specific residues responsible for CYP1A microdomain localization. Thus, CYP1A2 chimeras containing substitutions from homologous regions of CYP1A1 were expressed in HEK 293T/17 cells, and the localization was examined after solubilization with Brij 98. A CYP1A2 mutant with the three amino acids from CYP1A1 (VAG) at positions 27-29 of CYP1A2 was generated that showed a distribution pattern similar to those of CYP1A1/1A2 chimeras containing both the first 109 amino acids and the first 31 amino acids of CYP1A1 followed by remaining amino acids of CYP1A2. Similarly, the reciprocal substitution of three amino acids from CYP1A2 (AVR) into CYP1A1 resulted in a partial redistribution of the chimera into ordered microdomains. Molecular dynamic simulations indicate that the positive charges of the CYP1A1 and CYP1A2 linker regions between the N-termini and catalytic domains resulted in different depths of immersion of the N-termini in the membrane. The overlap of the distribution of positively charged residues in CYP1A2 (AVR) and negatively charged phospholipids was higher in the ordered than disordered microdomain. These findings identify three residues in the CYP1A N-terminus as a novel microdomain-targeting motif of the P450s and provide a mechanistic explanation for the differential microdomain localization of CYP1A.

Roles of Individual Human Cytochrome P450 Enzymes in Drug Metabolism

Our knowledge of the roles of individual cytochrome P450 (P450) enzymes in drug metabolism has developed considerably in the past 30 years, and this base has been of considerable use in avoiding serious issues with drug interactions and issues due to variations. Some newer approaches are being considered for "phenotyping" metabolism reactions with new drug candidates. Endogenous biomarkers are being used for noninvasive estimation of levels of individual P450 enzymes. There is also the matter of some remaining "orphan" P450s, which have yet to be assigned reactions. Practical problems that continue in drug development include predicting drug-drug interactions, predicting the effects of polymorphic and other P450 variations, and evaluating interspecies differences in drug metabolism, particularly in the context of "metabolism in safety testing" regulatory issues ["disproportionate (human) metabolites"]. SIGNIFICANCE STATEMENT: Cytochrome P450 enzymes are the major catalysts involved in drug metabolism. The characterization of their individual roles has major implications in drug development and clinical practice.

Functional characterization of CYP1 enzymes: Complex formation, membrane localization and function

CYP1A1, CYP1A2, and CYP1B1 have a high degree of sequence similarity, similar substrate selectivities and induction characteristics. However, experiments suggest that there are significant differences in their quaternary structures and function. The goal of this study was to characterize the CYP1 proteins regarding their ability to form protein-protein complexes, lipid microdomain localization, and ultimately function. This was accomplished by examining (1) substrate metabolism of the CYP1s as a function of NADPH-cytochrome P450 reductase (POR) concentration, and (2) quaternary structure, using bioluminescence resonance energy transfer (BRET). Both CYP1As were able to form BRET-detectable homomeric complexes, which was not observed with CYP1B1. When activities were measured as a function of [POR], CYP1A1 and CYP1B1 showed a hyperbolic response, consistent with mass-action binding; however, CYP1A2 produced a sigmoidal response, suggesting that the homomeric complex affected its function. Differences were observed in their ability to form heteromeric complexes. Whereas CYP1B1 and CYP1A1 formed a complex, neither the CYP1A1/CYP1A2 nor the CYP1B1/CYP1A2 pair formed BRET-detectable complexes. These proteins also differed in their lipid microdomain localization, with CYP1A2 and CYP1B1 residing in ordered membranes, and CYP1A1 in the disordered lipid regions. Taken together, despite their sequence similarities, there are substantial differences in quaternary structures and microdomain localization that can influence enzymatic activities. As these proteins exist in the endoplasmic reticulum with other ER-resident proteins, the P450s need to be considered as part of multi-enzyme systems rather than simply monomeric proteins interacting with their redox partners.

Assembling the P450 puzzle: on the sources of nonadditivity in drug metabolism

There is an increasing number of indications of an oversimplification in the premise that the cumulative properties of the human drug-metabolizing ensemble represent a simple aggregate of the properties of the constituting enzymes. Recent studies of the functional effects of hetero-association of multiple cytochrome P450 species and their interactions with metabolically related enzymes revealed a tight integration in the drug-metabolizing ensemble. In our opinion, the sources of interindividual variability in drug metabolism can be elucidated only when considering this ensemble as a multienzyme system, the functional parameters of which are determined by interactions between its constituents. In this article, we present a conceptual model providing a mechanistic explanation for the functional effects of the interactions between multiple P450 species and propose a clue to understanding the nonadditive behavior of the drug-metabolizing ensemble.

Identification of the contact region responsible for the formation of the homomeric CYP1A2•CYP1A2 complex

Previous studies showed that cytochrome P450 1A2 (CYP1A2) forms a homomeric complex that influences its metabolic characteristics. Specifically, CYP1A2 activity exhibits a sigmoidal response as a function of NADPH-cytochrome P450 reductase (POR) concentration and is consistent with an inhibitory CYP1A2•CYP1A2 complex that is disrupted by increasing [POR] (Reed et al. (2012) Biochem. J. 446, 489-497). The goal of this study was to identify the CYP1A2 contact regions involved in homomeric complex formation. Examination of X-ray structure of CYP1A2 implicated the proximal face in homomeric complex formation. Consequently, the involvement of residues L91-K106 (P1 region) located on the proximal face of CYP1A2 was investigated. This region was replaced with the homologous region of CYP2B4 (T81-S96) and the protein was expressed in HEK293T/17 cells. Complex formation and its disruption was observed using bioluminescence resonance energy transfer (BRET). The P1-CYP1A2 (CYP1A2 with the modified P1 region) exhibited a decreased BRET signal as compared with wild-type CYP1A2 (WT-CYP1A2). On further examination, P1-CYP1A2 was much less effective at disrupting the CYP1A2•CYP1A2 homomeric complex, when compared with WT-CYP1A2, thereby demonstrating impaired binding of P1-CYP1A2 to WT-CYP1A2 protein. In contrast, the P1 substitution did not affect its ability to form a heteromeric complex with CYP2B4. P1-CYP1A2 also showed decreased activity as compared with WT-CYP1A2, which was consistent with a decrease in the ability of P1-CYP1A2 to associate with WT-POR, again implicating the P1 region in POR binding. These results indicate that the contact region responsible for the CYP1A2•CYP1A2 homomeric complex resides in the proximal region of the protein.

Heteromeric complex formation between human cytochrome P450 CYP1A1 and heme oxygenase-1

P450 and heme oxygenase-1 (HO-1) receive their necessary electrons by interaction with the NADPH-cytochrome P450 reductase (POR). As the POR concentration is limiting when compared with P450 and HO-1, they must effectively compete for POR to function. In addition to these functionally required protein-protein interactions, HO-1 forms homomeric complexes, and several P450s have been shown to form complexes with themselves and with other P450s, raising the question, 'How are the HO-1 and P450 systems organized in the endoplasmic reticulum?' Recently, CYP1A2 was shown to associate with HO-1 affecting the function of both proteins. The goal of this study was to determine if CYP1A1 formed complexes with HO-1 in a similar manner. Complex formation among POR, HO-1, and CYP1A1 was measured using bioluminescence resonance energy transfer, with results showing HO-1 and CYP1A1 form a stable complex that was further stabilized in the presence of POR. The POR•CYP1A1 complex was readily disrupted by the addition of HO-1. CYP1A1 also was able to affect the POR•HO-1 complex, although the effect was smaller. This interaction between CYP1A1 and HO-1 also affected function, where the presence of CYP1A1 inhibited HO-1-mediated bilirubin formation by increasing the KmPOR•HO-1 without affecting the Vmaxapp. In like manner, HO-1 inhibited CYP1A1-mediated 7-ethoxyresorufin dealkylation by increasing the KmPOR•CYP1A1. Based on the mathematical simulation, the results could not be explained by a model where CYP1A1 and HO-1 simply compete for POR, and are consistent with the formation of a stable CYP1A1•HO-1 complex that affected the functional characteristics of both moieties.

Differential effects on human cytochromes P450 by CRISPR/Cas9-induced genetic knockout of cytochrome P450 reductase and cytochrome b5 in HepaRG cells

HepaRG cells are increasingly accepted as model for human drug metabolism and other hepatic functions. We used lentiviral transduction of undifferentiated HepaRG cells to deliver Cas9 and two alternative sgRNAs targeted at NADPH:cytochrome P450 oxidoreductase (POR), the obligate electron donor for microsomal cytochromes P450 (CYP). Cas9-expressing HepaRG^VC (vector control) cells were phenotypically similar to wild type HepaRG cells and could be differentiated into hepatocyte-like cells by DMSO. Genetic POR-knockout resulted in phenotypic POR knockdown of up to 90% at mRNA, protein, and activity levels. LC–MS/MS measurement of seven CYP-activities showed differential effects of POR-knockdown with CYP2C8 being least and CYP2C9 being most affected. Further studies on cytochrome b5 (CYB5), an alternative NADH-dependent electron donor indicated particularly strong support of CYP2C8-dependent amodiaquine N-deethylation by CYB5 and this was confirmed by genetic CYB5 single- and POR/CYB5 double-knockout. POR-knockdown also affected CYP expression on mRNA and protein level, with CYP1A2 being induced severalfold, while CYP2C9 was strongly downregulated. In summary our results show that POR/NADPH- and CYB5/NADH-electron transport systems influence human drug metabolizing CYPs differentially and differently than mouse Cyps. Our Cas9-expressing HepaRG^VC cells should be suitable to study the influence of diverse genes on drug metabolism and other hepatic functions.

Heme oxygenase-1 affects cytochrome P450 function through the formation of heteromeric complexes: Interactions between CYP1A2 and heme oxygenase-1

Heme oxygenase 1 (HO-1) and the cytochromes P450 (P450s) are endoplasmic reticulum-bound enzymes that rely on the same protein, NADPH-cytochrome P450 reductase (POR), to provide the electrons necessary for substrate metabolism. Although the HO-1 and P450 systems are interconnected owing to their common electron donor, they generally have been studied separately. As the expressions of both HO-1 and P450s are affected by xenobiotic exposure, changes in HO-1 expression can potentially affect P450 function and, conversely, changes in P450 expression can influence HO-1. The goal of this study was to examine interactions between the P450 and HO-1 systems. Using bioluminescence resonance energy transfer (BRET), HO-1 formed HO-1•P450 complexes with CYP1A2, CYP1A1, and CYP2D6, but not all P450s. Studies then focused on the HO-1-CYP1A2 interaction. CYP1A2 formed a physical complex with HO-1 that was stable in the presence of POR. As expected, both HO-1 and CYP1A2 formed BRET-detectable complexes with POR. The POR•CYP1A2 complex was readily disrupted by the addition of HO-1, whereas the POR•HO-1 complex was not significantly affected by the addition of CYP1A2. Interestingly, enzyme activities did not follow this pattern. BRET data suggested substantial inhibition of CYP1A2-mediated 7-ethoxyresorufin de-ethylation in the presence of HO-1, whereas its activity was actually stimulated at subsaturating POR. In contrast, HO-1-mediated heme metabolism was inhibited at subsaturating POR. These results indicate that HO-1 and CYP1A2 form a stable complex and have mutual effects on the catalytic behavior of both proteins that cannot be explained by a simple competition for POR.

Differential effects of variations in human P450 oxidoreductase on the aromatase activity of CYP19A1 polymorphisms R264C and R264H

Aromatase (CYP19A1) converts androgens into estrogens and is required for female sexual development and growth and development in both sexes. CYP19A1 is a member of cytochrome P450 family of heme-thiolate monooxygenases located in the endoplasmic reticulum and depends on reducing equivalents from the reduced nicotinamide adenine dinucleotide phosphate via the cytochrome P450 oxidoreductase coded by POR. Both the CYP19A1 and POR genes are highly polymorphic, and mutations in both these genes are linked to disorders of steroid biosynthesis. We have previously shown that R264C and R264H mutations in CYP19A1, as well as mutations in POR, reduce CYP19A1 activity. The R264C is a common polymorphic variant of CYP19A1, with high frequency in Asian and African populations. Polymorphic alleles of POR are found in all populations studied so far and, therefore, may influence activities of CYP19A1 allelic variants. So far, the effects of variations in POR on enzymatic activities of allelic variants of CYP19A1 or any other steroid metabolizing cytochrome P450 proteins have not been studied. Here we are reporting the effects of three POR variants on the aromatase activities of two CYP19A1 variants, R264C, and R264H. We used bacterially expressed and purified preparations of WT and variant forms of CYP19A1 and POR and constructed liposomes with embedded CYP19A1 and POR proteins and assayed the CYP19A1 activities using radiolabeled androstenedione as a substrate. With the WT-POR as a redox partner, the R264C-CYP19A1 showed only 15% of aromatase activity, but the R264H had 87% of aromatase activity compared to WT-CYP19A1. With P284L-POR as a redox partner, R264C-CYP19A1 lost all activity but retained 6.7% of activity when P284T-POR was used as a redox partner. The R264H-CYP19A1 showed low activities with both the POR-P284 L as well as the POR-P284 T. When the POR-Y607C was used as a redox partner, the R264C-CYP19A1 retained approximately 5% of CYP19A1 activity. Remarkably, The R264H-CYP19A1 had more than three-fold higher activity compared to WT-CYP19A1 when the POR-Y607C was used as the redox partner, pointing toward a beneficial effect. The slight increase in activity of R264C-CYP19A1 with the P284T-POR and the three-fold increase in activity of the R264H-CYP19A1 with the Y607C-POR point toward a conformational effect and role of protein-protein interaction governed by the R264C and R264H substitutions in the CYP19A1 as well as P284 L, P284 T and Y607C variants of POR. These studies demonstrate that the allelic variants of P450 when present with a variant form of POR may show different activities, and combined effects of variations in the P450 enzymes as well as in the POR should be considered when genetic data are available. Recent trends in the whole-exome and whole-genome sequencing as diagnostic tools will permit combined evaluation of variations in multiple genes that are interdependent and may guide treatment options by adjusting therapeutic interventions based on laboratory analysis.

The Biosynthesis of Enzymatically Oxidized Lipids

Enzymatically oxidized lipids are a specific group of biomolecules that function as key signaling mediators and hormones, regulating various cellular and physiological processes from metabolism and cell death to inflammation and the immune response. They are broadly categorized as either polyunsaturated fatty acid (PUFA) containing (free acid oxygenated PUFA "oxylipins", endocannabinoids, oxidized phospholipids) or cholesterol derivatives (oxysterols, steroid hormones, and bile acids). Their biosynthesis is accomplished by families of enzymes that include lipoxygenases (LOX), cyclooxygenases (COX), cytochrome P450s (CYP), and aldo-keto reductases (AKR). In contrast, non-enzymatically oxidized lipids are produced by uncontrolled oxidation and are broadly considered to be harmful. Here, we provide an overview of the biochemistry and enzymology of LOXs, COXs, CYPs, and AKRs in humans. Next, we present biosynthetic pathways for oxylipins, oxidized phospholipids, oxysterols, bile acids and steroid hormones. Last, we address gaps in knowledge and suggest directions for future work.

Proteomic response of gill microsomes of Crassostrea brasiliana exposed to diesel fuel water-accommodated fraction

Diesel fuel water-accommodated fraction (diesel-WAF) is a complex mixture of organic compounds that may cause harmful effects to marine invertebrates. Expression of microsomal proteins can be changed by oil exposure, causing functional alterations in endoplasmic reticulum (ER). The aim of this study was to investigate changes in protein expression signatures in microsomes of oysterl Crassostrea brasiliana (=C.gasar) gill after exposure to 10% diesel-WAF for 24 and 72 h. Protein expression signatures of gills of oysters exposed to diesel-WAF were compared to those of unexposed oysters using two-dimensional electrophoresis (2-DE) to identify differentially expressed proteins. A total of 458 protein spots with molecular weights between 30-75 kDa were detected by 2-DE in six replicates of exposed oyster proteomes compared to unexposed ones. Fourteen differentially expressed proteins (six up-regulated and eight down-regulated) were identified. They are: proteins related to xenobiotic biotransformation (cytochrome P450 6 A, NADPH-cytochrome P450 reductase); cytoskeleton (α-tubulin, β-tubulin, gelsolin); processing and degradation of proteins pathways (thioredoxin domain-containing protein E3 ubiquitin-protein ligase MIB2); involved in the biosynthesis of glycolipids and glycoproteins (beta-1,3-galactosyltransferase 1); associated with stress responses (glutamate receptor 4 and 14-3-3 protein zeta, corticotropin-releasing factor-binding protein); plasmalogen biosynthesis (fatty acyl-CoA reductase 1), and sodium-and chloride-dependent glycine transporter 2 and glyoxylate reductase/hydroxypyruvate reductase. Different patterns of protein responses were observed between 24 and 72 h-exposed groups. Expression pattern of microsomal proteins provided a first insight on the potential diesel-WAF effects at protein level in microsomal fraction of oyster gills and indicated new potential biomarkers of exposure and effect. The present work can be a basis for future ecotoxicological studies in oysters aiming to elucidate the molecular mechanisms behind diesel-WAF toxicity and for environmental monitoring programs.