Involvement of nonlamellar-prone lipids in the stability increase of human cytochrome P450 1A2 in reconstituted membranes

Molecular Dynamics Simulations Reveal Structural Differences among Allelic Variants of Membrane-Anchored Cytochrome P450 2D6

Cytochrome P450 2D6 (CYP2D6) is an enzyme that is involved in the metabolism of roughly 25% of all marketed drugs and therefore belongs to the most important enzymes in drug metabolism. CYP2D6 features a high degree of genetic polymorphism that can significantly affect the metabolic activity of an individual. In extreme cases, structural changes at the level of single amino acids can either increase its enzymatic activity abolishing the drug therapeutic effect or completely disable the enzyme and elevate drug plasma level potentially leading to adverse effects. In this study, starting from the crystal structure, we built a full-length membrane-anchored all-atom model of the wild-type CYP2D6 as well as five of its variants differing in the enzymatic activity. We validated our models with available experimental data and compared their structural properties with molecular dynamics simulations. The main focus of this study was to identify differences that could mechanistically explain the altered activity of the variants and improve our understanding of their functioning. We observed differences in the opening frequencies and minimal diameters of tunnels that connect the buried active site to the surrounding solvent environment. The variants CYP2D6*4 and CYP2D6*10 associated with missing or decreased activity showed less frequent opening of the tunnels compared to the wild-type. Both CYP2D6*10 and CYP2D6*17 showed a deprivation of an important ligand tunnel suggesting a feasible reason for their altered substrate specificity. Next, the altered fold at the N-terminal anchor region and the decreased active site volume caused by the amino acid mutations of the CYP2D6*4 variant offer an explanation for the absence of its metabolic activity. The mutations in CYP2D6*53 contributed to a significant enlargement of an important ligand tunnel and an extension of the active site cavity. This could explain the altered metabolic profile as well as the enhanced metabolic rates of this particular variant supporting its designation as a possible cause for the ultrarapid metabolizer phenotype. We believe these novel structural insights could advance the fields of personalized medicine and enzyme engineering. Furthermore, they could aid in guiding laboratory as well as computational experiments in the future.

Lipid molecules can induce an opening of membrane-facing tunnels in cytochrome P450 1A2

Cytochrome P450 1A2 (P450 1A2, CYP1A2) is a membrane-bound enzyme that oxidizes a broad range of hydrophobic substrates. The structure and dynamics of both the catalytic and trans-membrane (TM) domains of this enzyme in the membrane/water environment were investigated using a multiscale computational approach, including coarse-grained and all-atom molecular dynamics. Starting from the spontaneous self-assembly of the system containing the TM or soluble domain immersed in randomized dilauroyl phosphatidylcholine (DLPC)/water mixture into their respective membrane-bound forms, we reconstituted the membrane-bound structure of the full-length P450 1A2. This structure includes a TM helix that spans the membrane, while being connected to the catalytic domain by a short flexible loop. Furthermore, in this model, the upper part of the TM helix interacts directly with a conserved and highly hydrophobic N-terminal proline-rich segment of the catalytic domain; this segment and the FG loop are immersed in the membrane, whereas the remaining portion of the catalytic domain remains exposed to aqueous solution. The shallow membrane immersion of the catalytic domain induces a depression in the opposite intact layer of the phospholipids. This structural effect may help in stabilizing the position of the TM helix directly beneath the catalytic domain. The partial immersion of the catalytic domain also allows for the enzyme substrates to enter the active site from either aqueous solution or phospholipid environment via several solvent- and membrane-facing tunnels in the full-length P450 1A2. The calculated tunnel dynamics indicated that the opening probability of the membrane-facing tunnels is significantly enhanced when a DLPC molecule spontaneously penetrates into the membrane-facing tunnel 2d. The energetics of the lipid penetration process were assessed by the linear interaction energy (LIE) approximation, and found to be thermodynamically feasible.

Cytochrome P450 Organization and Function Are Modulated by Endoplasmic Reticulum Phospholipid Heterogeneity

Cytochrome P450s (P450s) comprise a superfamily of proteins that catalyze numerous monooxygenase reactions in animals, plants, and bacteria. In eukaryotic organisms, these proteins not only carry out reactions necessary for the metabolism of endogenous compounds, but they are also important in the oxidation of exogenous drugs and other foreign compounds. Eukaryotic P450 system proteins generally reside in membranes, primarily the endoplasmic reticulum or the mitochondrial membrane. These membranes provide a scaffold for the P450 system proteins that facilitate interactions with their redox partners as well as other P450s. This review focuses on the ability of specific lipid components to influence P450 activities, as well as the role of the membrane in P450 function. These studies have shown that P450s and NADPH-cytochrome P450 reductase appear to selectively associate with specific phospholipids and that these lipid-protein interactions influence P450 activities. Finally, because of the heterogeneous nature of the endoplasmic reticulum as well as other biologic membranes, the phospholipids are not arranged randomly but associate to generate lipid microdomains. Together, these characteristics can affect P450 function by 1) altering the conformation of the proteins, 2) influencing the P450 interactions with their redox partners, and 3) affecting the localization of the proteins into specific membrane microdomains.

The Localization of Cytochrome P450s CYP1A1 and CYP1A2 into Different Lipid Microdomains Is Governed by Their N-terminal and Internal Protein Regions

In cellular membranes, different lipid species are heterogeneously distributed forming domains with different characteristics. Ordered domains are tightly packed with cholesterol, sphingomyelin, and saturated fatty acids, whereas disordered domains contain high levels of unsaturated fatty acids. Our laboratory has shown that membrane heterogeneity affects the organization of cytochrome P450s and their cognate redox partner, the cytochrome P450 reductase (CPR). Despite the high degree of sequence similarity, CYP1A1 was found to localize to disordered regions, whereas CYP1A2 resided in ordered domains. We hypothesized that regions of amino acid sequence variability may contain signal motifs that direct CYP1A proteins into ordered or disordered domains. Thus, chimeric constructs of CYP1A1 and CYP1A2 were created, and their localization was tested in HEK293T cells. CYP1A2, containing the N-terminal regions from CYP1A1, no longer localized in ordered domains, whereas the N terminus of CYP1A2 partially directed CYP1A1 into ordered regions. In addition, intact CYP1A2 containing a 206-302-residue peptide segment of CYP1A1 had less affinity to bind to ordered microdomains. After expression, the catalytic activity of CYP1A2 was higher than that of the CYP1A1-CYP1A2 chimera containing the N-terminal end of CYP1A1 with subsaturating CPR concentrations, but it was approximately equal with excess CPR suggesting that the localization of the CYP1A enzyme in ordered domains favored its interaction with CPR. These data demonstrate that both the N-terminal end and an internal region of CYP1A2 play roles in targeting CYP1A2 to ordered domains, and domain localization may influence P450 function under conditions that resemble those found in vivo.

Relationship between CYP1A2 localization and lipid microdomain formation as a function of lipid composition

Cytochrome P450 (P450) function requires the interaction of P450 and NADPH-cytochrome P450 reductase (CPR) in membranes, and is frequently studied using reconstituted systems composed solely of phosphatidylcholine. There is increasing evidence that other endoplasmic reticulum (ER) lipids can affect P450 structure, activity, and interactions with CPR. Some of these lipid effects have been attributed to the formation of organized liquid-ordered (l(o)) domains. The goal of this study was to determine if l(o) domains were formed in P450 reconstituted systems mimicking the ER membrane. CYP1A2, when incorporated in "ER-like" lipid vesicles, displayed detergent insolubility after treatment with Brij 98 and centrifugation in a sucrose gradient. Lipid probes were employed to identify domain formation in both ER-like vesicles and model membranes known to form l(o) domains. Changes in fluorescence resonance energy transfer (FRET) using an established donor/acceptor FRET pair in both ER-like and model l(o)-forming systems demonstrated the coexistence of l(o)- and liquid-disordered domains as a function of cholesterol and sphingomyelin content. Similarly, 6-dodecanoyl-2-dimethylaminonaphthalene (laurdan), a probe that reports on membrane organization, showed that cholesterol and sphingomyelin increased membrane order. Finally, brominated-phosphatidylcholine allowed for monitoring of the location of both CPR and CYP1A2 within the l(o) regions of ER-like systems. Taken together, the results demonstrate that ER-like vesicles generate microdomains, and both CYP1A2 and CPR predominantly localize into l(o) membrane regions. Probe fluorescent responses suggest that lipid microdomains form in these vesicles whether or not enzymes are included in the reconstituted systems. Thus, it does not appear that the proteins are critical for stabilizing l(o) domains.

Membrane phospholipid augments cytochrome P4501a enzymatic activity by modulating structural conformation during detoxification of xenobiotics

Cytochrome P450 is a superfamily of membrane-bound hemoprotein that gets involved with the degradation of xenobiotics and internal metabolites. Accumulated body of evidence indicates that phospholipids play a crucial role in determining the enzymatic activity of cytochrome P450 in the microenvironment by modulating its structure during detoxification; however, the structure-function relationship of cytochrome P4501A, a family of enzymes responsible for degrading lipophilic aromatic hydrocarbons, is still not well defined. Inducibility of cytochrome P4501A in cultured catfish hepatocytes in response to carbofuran, a widely used pesticide around the world, was studied earlier in our laboratory. In this present investigation, we observed that treating catfish with carbofuran augmented total phospholipid in the liver. We examined the role of phospholipid on the of cytochrome P4501A-marker enzyme which is known as ethoxyresorufin-O-deethylase (EROD) in the context of structure and function. We purified the carbofuran-induced cytochrome P4501A protein from catfish liver. Subsequently, we examined the enzymatic activity of purified P4501A protein in the presence of phospholipid, and studied how the structure of purified protein was influenced in the phospholipid environment. Membrane phospholipid appeared to accelerate the enzymatic activity of EROD by changing its structural conformation and thus controlling the detoxification of xenobiotics. Our study revealed the missing link of how the cytochrome P450 restores its enzymatic activity by changing its structural conformation in the phospholipid microenvironment.

Effects of phospholipids on the functional regulation of tBID in membranes

The functional interplay between tBID and phospholipids was investigated in this study. The binding of tBID to model membranes was increased by an incorporation of phosphatidylserine (PS) into the liposomes. Using limited proteolysis and mass spectrometry, two peptide regions, which correspond to Ser(100)-Arg(114) and His(89)-Arg(114) in BID, revealed the specific PS-binding site. tBID also decreased the light scattering values of PS-containing liposomes and increased the leakage of fluorescent dye encapsulated in vesicles, which suggest that tBID reduces membrane integrity by fragmentation. The membrane fragmentation by tBID was also observed using confocal and transmission electron microscopy. The activity of tBID paralleled results that were obtained with cardiolipin (CL)-containing membranes. However, other anionic phospholipids had little effect. CL- and PS-induced conformational changes of tBID were observed by circular dichroism and intrinsic fluorescence. CL and PS also stimulated the insertion of BID into lipid monolayers. tBID stimulated the leakage of Ca(2+) from purified microsomes and mitochondria in a protein concentration-dependent manner. In contrast, BID showed significantly reduced effects when compared to tBID in all of the experiments performed. These results suggest that tBID specifically interacts with PS as well as CL and decreases membrane integrity without the aid of other pro-apoptotic proteins.

Simultaneous measurement of CYP1A2 activity, regioselectivity, and coupling: Implications for environmental sensitivity of enzyme-substrate binding

The cytochrome P450 (CYP) reaction mechanism often yields a broad array of coupled and uncoupled products from a single substrate. While it is well known that reaction conditions can drastically affect the rate of P450 catalysis, their effects on regioselectivity and coupling are not well characterized. To investigate such effects, the CYP1A2 oxidation of 7-ethoxymethoxy-3-cyanocoumarin (EOMCC) was examined as a function of buffer type, buffer concentration, pH, and temperature. A high-throughput, optical method was developed to simultaneously measure the rate of substrate depletion, NADPH depletion, and generation of the O-dealkylated product. Increasing the phosphate buffer concentration and temperature increased both the NADPH and EOMCC depletion rates by 6-fold, whereas coupling was constant at 7.9% and the regioselectivity of O-dealkylation to other coupled pathways was constant at 21.7%. Varying the buffer type and pH increased NADPH depletion by 2.5-fold and EOMCC depletion by 3.5-fold; however, neither coupling nor regioselectivity was constant, with variations of 14.4% and 21.6%, respectively. Because the enzyme-substrate binding interaction is a primary determinant of both coupling and regioselectivity, it is reasonable to conclude that ionic strength, as varied by the buffer concentration, and temperature alter the rate without affecting binding while buffer type and pH alter both.

Conformational change of Escherichia coli signal recognition particle Ffh is affected by the functionality of signal peptides of ribose-binding protein

We examined the effects of synthetic signal peptides, wild-type (WT) and export-defective mutant (MT) of ribose-binding protein, on the conformational changes of signal recognition particle 54 homologue (Ffh) in Escherichia coli. Upon interaction of Ffh with WT peptide, the intrinsic Tyr fluorescence, the transition temperature of thermal unfolding, and the GTPase activity of Ffh decreased in a peptide concentration-dependent manner, while the emission intensity of 8-anilinonaphthalene-1-sulfonic acid increased. In contrast, the secondary structure of the protein was not affected. Additionally, polarization of fluorescein-labeled WT increased upon association with Ffh. These results suggest that WT peptide induces the unfolded states of Ffh. The WT-mediated conformational change of Ffh was also revealed to be important in the interaction between SecA and Ffh. However, MT had marginal effect on these conformational changes suggesting that the in vivo functionality of signal peptide is important in the interaction with Ffh and concomitant structural change of the protein.

Physical incorporation of NADPH-cytochrome P450 reductase and cytochrome P450 into phospholipid vesicles using glycocholate and Bio-Beads

In a previous study from our laboratory (Drug Metab Dispos 34: 660-666, 2006), we found several limitations with published methods (cholate gel filtration and cholate dialysis) for the incorporation of cytochromes P450 and P450 reductase into phospholipid vesicles. We found that a significant proportion of reductase was not incorporated in the vesicles when the amount of reductase was equal to or greater than that of CYP2B4 in the systems reconstituted with phosphatidylcholine. Furthermore, implementation of these methods compromised the ability of the CYP2B4 to form a ferrous carbon monoxy complex. In the current study, a comparison of results using the detergent-dialysis method with five similar detergents having the "bile salt" ring structure showed that glycocholate results in the greatest incorporation of reductase and the least loss in the ferrous carbon monoxy CYP2B4 complex. The method is further improved by using Bio-Beads SM-2 to remove detergent instead of the lengthy dialysis procedure or size exclusion chromatography that significantly dilutes the protein and lipid concentrations of the preparation. The method is shown to be applicable over a range of lipid/CYP2B4 ratios, and by using assay methods for total lipid, reductase, and CYP2B4, this improved reconstitution method resulted in increased incorporation efficiencies while minimizing the protein degradation inherent with these procedures.

Heterologous expression and characterization of wild-type human cytochrome P450 1A2 without conventional N-terminal modification in Escherichia coli

In this study, wild-type human CYP1A2 without the conventional N-terminal modification (second codon GCT) or the truncation of the N-terminal hydrophobic region was functionally expressed in Escherichia coli. Its enzymatic properties were compared with N-terminally modified CYP1A2. Although modified CYP1A2 is almost all high-spin, some wild-type CYP1A2 shifted to low-spin. Spectral binding titrations with several ligands could be performed with wild-type enzyme, but not with modified enzyme. Kinetic parameters for several substrates were similar for the two CYP1A2 enzymes. However, the oxidation rates of phenacetin by modified enzyme were approximately 2-fold higher than those by wild-type enzyme. The intermolecular isotope effects were approximately 2 for phenacetin O-deethylation catalyzed by both enzymes. However, the wild-type enzyme, but not the modified enzyme, increased C-hydroxylation when O-deethylation rates were lowered by deuterium substitution. Molecular switching indicates that phenacetin rotates within the active site of wild-type enzyme and suggests a looser conformation in the active site of the wild-type enzyme than of the modified enzyme. These results reveal that the overall enzymatic properties of wild-type CYP1A2 enzyme are quite similar to those of modified CYP1A2, although its active site environment seems to differ from that of the modified enzyme.

Lateral segregation of anionic phospholipids in model membranes induced by cytochrome P450 2B1: bi-directional coupling between CYP2B1 and anionic phospholipid

The lateral segregation of anionic phospholipids phosphatidic acid (PA), phosphatidylinositol (PI), and phosphatidylserine (PS) was detected after addition of cytochrome P450 2B1 (CYP2B1). The tendency of lipid clustering was highly dependent on the type of anionic phospholipids examined. PA was the most highly clustered while PI and PS clustered to a lesser degree. Moreover, liposomes containing anionic phospholipids form anionic phospholipid-rich microdomains in the presence of CYP2B1. Anionic phospholipids (mostly notably PA) also increased the ability of CYP2B1 to bind to lipid monolayers. In addition to the ability of CYP2B1 to modulate the physical properties of the membrane, the membrane itself can have reciprocal effects on the activity and conformation of CYP2B1. The catalytic activity of CYP2B1 increased as a function of anionic phospholipid concentration and in the presence of 10 mol% PA, the activity increased by 85%. These results suggest a bi-directional coupling between the CYP2B1 and anionic phospholipids.

Production of polyclonal antibodies against peptide antigens using polystyrene beads as a carrier

We describe a method for producing polyclonal antibody (PAb) against peptide antigen (the epitope region of cytochromes P450 1A2 and 3A4, HA tag, FLAG tag, and c-myc tag) using a polymer bead as a carrier. Carboxylated polystyrene beads (0.05, 0.5, 1, and 2 microm diam) were conjugated with the peptide, which is known as an epitope region, using a chemical cross-linker and injected into rabbits with adjuvant. By immunizing the animals intradermally twice at two week intervals, an immunoblot assay with anti-sera showed the efficient generation of the PAb.

Diacylglycerols, multivalent membrane modulators

Diacylglycerols are second messengers confined to biomembranes and, although relatively simple molecules from the structural point of view, they are able of triggering a surprisingly wide range of biological responses. Diacylglycerols are recognized by a well conserved protein motif, such as the C1 domain. This domain was observed for the first time in protein kinases C but is now known to be present in many other proteins. The effect of diacylglycerols is not limited to binding to C1 domains and they are able to alter the biophysical properties of biomembranes and hence modulate the activity of membrane associated proteins and also facilitate some processes like membrane fusion.