Next Article in Journal
Relative and Absolute Quantification of Aberrant and Normal Splice Variants in HBBIVSI−110 (G > A) β-Thalassemia
Next Article in Special Issue
Dose-Dependent Effects of Resveratrol on Cisplatin-Induced Hearing Loss
Previous Article in Journal
Electrophile-Induced Conformational Switch of the Human TRPA1 Ion Channel Detected by Mass Spectrometry
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 21
Issue 18
10.3390/ijms21186669
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Interaction Modes of Microsomal Cytochrome P450s with Its Reductase and the Role of Substrate Binding

by
Francisco Esteves
1,*,
Philippe Urban
2,
José Rueff
1,
Gilles Truan
2 and
Michel Kranendonk
1,*
1
Center for Toxicogenomics and Human Health (ToxOmics), Genetics, Oncology and Human Toxicology, NOVA Medical School/Faculty of Medical Sciences, Universidade NOVA de Lisboa, 1169-056 Lisbon, Portugal
2
TBI, Université de Toulouse, CNRS, INRAE, INSA, CEDEX 04, 31077 Toulouse, France
*
Authors to whom correspondence should be addressed.
Int. J. Mol. Sci. 2020, 21(18), 6669; https://doi.org/10.3390/ijms21186669
Submission received: 17 July 2020 / Revised: 7 September 2020 / Accepted: 10 September 2020 / Published: 11 September 2020
(This article belongs to the Special Issue Cytochromes P450: Drug Metabolism, Bioactivation and Biodiversity 3.0)
Browse Figures
Review Reports Versions Notes

Abstract

:
The activity of microsomal cytochromes P450 (CYP) is strictly dependent on the supply of electrons provided by NADPH cytochrome P450 oxidoreductase (CPR). The variant nature of the isoform-specific proximal interface of microsomal CYPs implies that the interacting interface between the two proteins is degenerated. Recently, we demonstrated that specific CPR mutations in the FMN-domain (FD) may induce a gain in activity for a specific CYP isoform. In the current report, we confirm the CYP isoform dependence of CPR’s degenerated binding by demonstrating that the effect of four of the formerly studied FD mutants are indeed exclusive of a specific CYP isoform, as verified by cytochrome c inhibition studies. Moreover, the nature of CYP’s substrate seems to have a modulating role in the CPR:CYP interaction. In silico molecular dynamics simulations of the FD evidence that mutations induces very subtle structural alterations, influencing the characteristics of residues formerly implicated in the CPR:CYP interaction or in positioning of the FMN moiety. CPR seems therefore to be able to form effective interaction complexes with its structural diverse partners via a combination of specific structural features of the FD, which are functional in a CYP isoform dependent manner, and dependent on the substrate bound.

1. Introduction

Human mitochondrial and microsomal cytochromes P450 (CYP) form a superfamily of membrane-bound hemeproteins with monooxygenase activity [1,2]. Microsomal CYP isoforms are dependent on NADPH cytochrome P450 oxidoreductase (CPR)—containing flavin adenine dinucleotide (FAD) (reductase) and flavin mononucleotide (FMN) (transporter) domains—for supply of electrons in their catalytic activities to metabolize a myriad of compounds, both endobiotics including vitamins, steroids and hormones, and xenobiotics, such as therapeutic drugs, carcinogens and environmental chemicals [3,4,5].
The study of affinity and binding in the CPR:CYP microsomal monooxygenase system is crucial for a better understanding of the structural features which guide the formation of an efficient enzyme complex, determinant for electron transfer (ET) from CPR to CYP, thus potentially modulating physiologically important microsomal CYPs activities. Although all CYPs are likely to diverge from a common ancestor gene, their interactions with redox partners seem to have evolved differently, depending on the cellular localization and type of substrates to metabolize [2,6]. Mitochondrial CYPs exhibit conserved basic residues at their proximal side for effective interaction with their cognate redox partner: [2Fe-2S] adrenodoxin (Adx) [7]. Contrarily, microsomal CYPs do not possess such signature sequences to guide their interaction with the FMN-domain (FD) of CPR for ET [6]. This suggests that CPR:CYP electron transport system is more versatile and does not depend on uniquely conserved structural features present on CPR–CYPs interaction interface. Instead, it relies on a degenerate interaction between the FD and the proximal side of the structurally diverse microsomal CYP isoforms [6,8,9,10].
The precise mechanism governing this degenerated protein-protein interaction and how CPR can sustain efficient ET with so many structurally diverse partners has been the subject of many studies over the last decades (reviewed in [2,5,10]). The overall steady-state ET from CPR to microsomal CYPs seems to occur via a concerted mechanism involving the FD and FAD-domain within CPR (inter-flavin ET) and the FD and CYP proximal side (inter-protein ET) [3,4,5]. Multiple studies reported by our group and others indicate that effective docking and ET from CPR to its many redox partners are enabled by: (i) The existence of multiple open conformations, the result of CPR’s open-closed dynamics, that can be sampled by its diverse redox partners [11,12,13,14,15]; and (ii) affinity probing of these open conformers, guided by electron receptor specific affinity for the degenerated interface/binding site of CPR [6,8,9,16,17]. Together with the membrane anchoring segment, implicated in the preliminary association of the CPR with CYPs and in substrate entrance/binding to CYP’s substrate-access channel [18,19,20,21], these features are thus required for an effective ET by CPR to multiple CYPs.
However, these “one-serve-all” types of interactions may still be dependent on the CYP isoform. It is likely that multiple CPR conformations and binding modes are probed in a selection process by which each of the structurally diverse CYP isoforms binds to (a) more favorable conformation(s) of the CPR [5]. Additionally, the in vivo existence of a vast excess of CYPs over CPR, implies competition between individual CYP isoforms in binding to CPR. This evidences the role of affinity in the formation of effective transient enzyme complexes, rather than stable complexes [11,22,23].
Through structural and functional studies, specific patches of acidic residues surrounding the flavin isoalloxazine ring and on the binding interface of the FD were identified. These acidic residues drive the effective docking between CPR and its multiple electron acceptors via electrostatic interactions (Supplementary Figure S1) [14,17,24,25,26]. However, the relative importance of these residues for interaction seems to be dependent on the CPR:redox-partner couple [2,5]. Additionally, basic residues on the proximal side of CYPs seem to steer the negatively charged convex binding-interface of the FD to CYP’s concave docking site near the heme moiety [9]. Interestingly, structural information on CYP substrate/ligand complexes and conformational dynamics for substrate access demonstrated that CYP’s active site adapts to substrates of diverse sizes and shapes, which in turn may have implications in the CPR:CYP interaction as the structural rearrangement associated with substrate binding was shown to extends its effects to the proximal side of CYPs [1,2,27,28,29,30,31].
We previously studied the underlying molecular mechanisms on CPR:CYP interaction and ET, by creating and selecting seven mutants in the FD of human CPR, each supporting a gain in activity for a specific CYP isoform [16]. Most of the identified positions co-localize with specific acidic stretches—formerly indicated to be involved in CPR:CYP interactions—and/or with mutations found in naturally occurring CPR variants—previously shown to cause CYP isoform dependent effects [4,22,32,33,34,35]. The present study aims: (i) To explore CPR’s degenerated binding mode by demonstrating that the generated FD mutants of CPR are exclusive of a specific CYP isoform and thus is detrimental to interaction of CPR with other CYP forms; (ii) to determine the role of substrate-binding in CYPs in modulating this specific recognition. From the set of previously identified CYP-specific CPR mutants [16], we selected P117H and G144C (specific for CYP1A2), G175D (specific for CYP2A6) and N151D (specific for CYP3A4) and crossed-combined each of them with the two other CYPs, i.e., the ones not used in their selection. Effects were evaluated using extended enzyme kinetic analysis with different fluorogenic standard probes, cytochrome c (cyt c) competition, and multiple structural (in silico) analyses. The results presented herein further confirm the preferential interaction of CYP isoforms with FD mutants and describe the structural features that promote the CPR–CYP recognition in an isoform specific manner, thus potentially modulating CPR’s promiscuous interactions with structurally diverse CYP isoforms and CYP-substrate binding states.

2. Results

2.1. Characterization of the Membrane Fractions

The characterization of contents of heterologous expressed proteins found in the membrane fractions (isolated using the bacterial CPR/CYP co-expression system, BTC) [22,36] containing the wildtype CPR (CPRwt) and mutant CPRs (CPRmut) combined with CYP1A2, 2A6 or 3A4 is shown in Table 1. Co-expression of CPRmut with CYP1A2, 2A6, or 3A4 was achieved with a stoichiometry approximating the one obtained with CPRwt for each CYP, demonstrating no significant differences (p > 0.05). As such, any difference detected in the CYP activities can be ascribed to differences in the CPRmut:CYP coupling. Indeed, the CPR:CYP stoichiometries measured in membrane fractions were similar to those described previously for the BTC system in our former studies [11,16,33,36,37,38,39], and comparable with the ranges of those reported for human liver microsomes, which range from 1:2 to 1:13 [40,41]. These ratios evidence a submolar presence of CPR vs. CYP, which is recapitulated in the BTC system, as already demonstrated in our former studies [11,22,36,39].

2.2. Detailed Enzyme Activity

2.2.1. Reaction Velocities of the Different CPRmut:CYP Combinations

Comparison of the kinetic parameters for each CPRmut/CYP couple are presented in Figure 1 and detailed numbers in Supplementary Table S1. We confirm the significantly higher relative turnovers of these mutants in combination with specific CYPs, namely CPRP117H:CYP1A2 and CPRG144C:CYP1A2 [ethoxyresorufin (EthR) O-deethylation, EROD], CPRG175D:CYP2A6 (coumarin 7-hydroxylation, C7H) and CPRN151D:CYP3A4 [dibenzylfluorescein (DBF) O-debenzylation, DBODF], as was observed in our previous study [16]. Except for CPRG175D:CYP3A4 (increased kcat, 1.22 X, p < 0.005) (Figure 1C), all the new CPRmut:CYP combinations demonstrated equal (CPRP117H:CYP2A6; CPRN151D:CYP2A6) or significantly lower activities (CPRG175D:CYP1A2; CPRN151D:CYP1A2; CPRG144C:CYP2A6; CPRP117H:CYP3A4; CPRG144C:CYP3A4) when compared to the corresponding CPRwt:CYP couples. According to this data, one single mutation in the FD is able to induce increased ET efficiency for one specific CPR:CYP combination, while maintaining or decreasing this activity for other couples. In few cases, the apparent affinity constant (KM) changed with the FD mutant CPRs. The apparent affinity constant of the CYP enzyme complex is a combination of multiple factors such as substrate entrance, binding and product egress but also ET rates. We thus hypothesize that mutations in the FD, by either affecting the binding of CPR to CYPs or directly modifying ET, may affect the apparent affinity constants.

2.2.2. Reaction Velocities of CPRmut:CYP1A2 Using Different Substrates

Microsomal CYP turn-overs can be modulated by many parameters and, amongst those, the efficiency of the complex formation between CPR and CYP and the substrate nature and presence in the active site are known to be critical [5,42]. We thus compared the effects of the four FD mutations of CPR with CYP1A2 using two other substrates (Figure 2 and Supplementary Table S2). Globally, the effects of the mutations are substrate dependent: The relative kcat profiles observed for the CPRmut:CYP1A2 in [3-cyano-7-ethoxycoumarin (CEC) O-dealkylation, CECOD] (Figure 2A) and DBODF (Figure 2B) are distinct from the ones of EROD. Contrarily to the activity profiles observed for EROD (Figure 1A, i.e., increased CYP1A2 activity with P117H and G144C; p < 0.005), the kcat of CYP1A2 generated with these CPR mutants demonstrated CPRwt-like profiles for CECOD (Figure 2A) and DBODF (Figure 2B). The G175D mutant (selected previously for increased CYP2A6-C7H activity) showed a decrease in CYP1A2-EROD (p < 0.05) and -CECOD (p < 0.05) activities, when compared with the CPRwt. Inversely, this mutant increased considerably the CYP1A2-DBODF activity (p < 0.005). A decrease in CYP1A2-EROD (p < 0.005) and -CECOD (p < 0.005) activities and a CPRwt-like CYP1A2-DBODF activity were observed for the N151D mutant (selected previously for increased CYP3A4-DBODF activity).
When comparing the relative kcat profiles of the different reactions, DBODF, a typical CYP3A4 substrate, showed turnover profiles (Figure 2B) distinct from the other substrates studied. This may be explained by DBF’s less favorable chemical structure, when compared to the other, more typical CYP1A2 substrates studied. P117H, G144C and N151D demonstrated relative DBODF turnovers very similar to the ones of CPRwt:CYP1A2, while for G175D we observed a completely different pattern. For this specific mutant, relative kcat values increased dramatically (1.83 X, p < 0.005). This was not observed, either for the other two substrates or by any of the other CPRmut:CYP1A2 combinations. Additionally, the CPRmut:CYP1A2 reaction velocities profiles for the CYP1A2 mediated de-alkylation of methoxyresorufine (MROD) (Supplementary Figure S2) were identical to the ones observed for the EROD reaction. These results are not surprising because substrates (typical CYP1A2 substrates) are similar and deviate only by one CH2 group.

2.3. Cytochrome c Competition with CPRmut/CYP

In two recent studies, we applied cytochrome b5 (CYB5) to probe affinity differences between CPR mutants and CYP [11,16]. Although informative, the inhibition data were difficult to interpret due to a dual effect by CYB5, namely an initial stimulating effect on CYP activity (depending on the CYP-isoform) and subsequent inhibitory effects at higher CYB5 levels when CPR binding is prevented [16]. An alternative competitive inhibitor is cyt c. Structural and functional studies of the FD binding site of rat CPR evidenced that the surface patches involved in protein:protein contacts with both CYP and cyt c are probably overlapping (detailed in Supplementary Table S3, and Figures S1 and S3) [5,14,17,24,25,43,44,45,46,47,48]. Therefore, cyt c competition experiments may shed light on differences in affinity between the various FD mutants when binding to different CYPs, without any confounding stimulatory effects when CYB5 is applied for this purpose. Relative velocities (kobs mutant/kobs wild-type) were determined for the four FD mutants with various cyt c/CPR ratios for CYP1A2, 2A6, or 3A4-mediated reactions (Figure 3). Globally, cyt c seems to have an inhibitory effect on CYP dependent reactions, suggesting direct competition between cyt c and CYP for CPR binding. The relative velocity traces of the CPRmut:CYP couples CPRP117H:CYP1A2, CPRG144C:CYP1A2, CPRG175D:CYP2A6 and CPRN151D:CYP3A4 show less inhibited activities than those from the corresponding CPRwt/CYP couples. Except for CPRN151D:CYP2A6 (Figure 3B), which showed a pronounced inhibition by cyt c, the profiles of all other CPRmut:CYP couples were comparable to the ones observed for the CPRwt:CYP. One-way analysis of variance demonstrated significant differences (p < 0.05) between the distribution profile of the means of the activity traces of CPR:CYP1A2 (Figure 3A) and CPR:CYP2A6 (Figure 3B) couples.
These results confirm that the competition exerted by cyt c with CYP for CPR-binding (i.e., ability to inhibit CPR:CYP interaction and ET) depends on the strength of interaction between CPR and a particular CYP, thus reinforcing the notion that the four mutants have indeed increased affinities for the CYP with which they were selected. It is also most probable that the level of ionic interactions in the CPR:CYP complex are CYP isoform dependent, which our former study demonstrated [11,16] and current data confirmed.

2.4. Structural Analysis of CPR Mutants

In order to interpret the biochemical properties of the mutants with structural modifications, we performed molecular dynamics simulations. Four independent simulations were run for the wild-type and each of the mutants, in water, during 80 ns. A first analysis of root mean square deviations (RMSD) along the progress of the simulations revealed a common pattern: A rise in the RMSD in the first 20 ns, corresponding to an increase in the energy as the starting coordinates were directly issued from a crystal structure and then a stabilization of the RMSD corresponding to a stable energy state, with rather small variations between the different simulations (Supplementary Figures S4 and S5). Comparisons of the average structures from the simulations with the FMN domain from the crystal structure reveals minor deviations which mostly affect helix A of the FMN domain. This helix, located at the N-terminus, is probably more subjected to translation movements compared to the rest of the structure. In one of the simulations (P117H, simulation 1), the time to reach the energy equilibrium is longer (40 ns instead of ~20 ns). In the average structure from this simulation, the D helix is slightly tilted toward the solvent. The average RMSD value for N151D in simulation 3 is also slightly higher, which might reflect a larger deviation from the starting structure. However, for both forms and simulations, the distribution of RMSD after the equilibrium is similar to the others, implying that the equilibrium, when reached, was the same for both particular simulations. In any case, the strongest deviations in RMSDs do not cluster (albeit possibly for P117H simulation) close to the mutated positions (Supplementary Figure S5). Therefore, the various selected mutations do not, in an 80 ns simulation, perturb the overall structure or folding of the FD, but probably have rather specific and localized effects. We thus looked more closely into the different simulations in order to reveal these localized effects and link them to known properties for interactions between CPR and CYPs.
We first examined the effects of P117H and N151D mutations (Figure 4) on the pKa of acidic residues around both positions. For P117H, the presence of a histidine at the position 117 leads to a drastic change of the pKa of D116 and E118, but not E119 (Figure 4A). This change can be explained by the short distances between the ND1 (tele or τ in IUPAC numbering) atom of H117 and the oxygen atoms of the carboxylic functions of either D116 in simulation 1 and 4 (Supplementary Figure S6A) or E118 in simulation 2, 3, and 4 (Supplementary Figure S6B), leading to potential hydrogen bonds in both simulations (Supplementary Figure S6C). This difference between the two simulations is interesting as it reveals the two possible states of the sidechain of H117 which may be interchanged in longer simulations times. In any case, these changes propagate to the electrostatic potential which is less negatively charged in the H117 simulations than the one present in the wild-type simulations (Figure 4B). On the contrary, the presence of an aspartate residue instead of an asparagine at position 151 does not seem to modify the pKa of the surrounding acidic residues (Figure 4C). However, the pKa of D151, while being one pH unit higher than the theoretical value of aspartate (Figure 4C), introduces a supplementary negative charge in this area (Supplementary Figure S7). This is also reflected on the electrostatic potential (Figure 4D), where, contrarily to P117H mutant (Figure 4B), the added negative charge leads to a quite visible increase of the negative potential in the same region. Overall, both mutations affect the electrostatic potential of the exact same area, although in a complete opposite way, decreasing the charge density for P117H and increasing the negative potential for N151D. Interestingly, D116 and E118 were formerly identified to be involved in hydrogen bonds with R440 in a docking model of the human FD of CPR with CYP2D6 [49,50]. Furthermore, a mutation of this arginine residue into histidine leads to a 5- to 100-fold decrease of CYP2D6 activity [49].
For G114C mutation, the situation is quite different as the glycine residue at this position in CPR is well conserved amongst species. As seen in Figure 5A, the cysteine at this position extends toward the surface but does not seem to perturb the positions of the residues surrounding the FMN moiety except for E145 for which the sidechain moves backward to the FMN. Analysis of the pKa of the ionizable residues close to the FMN (E145, E182 and H183) reveals no significant modifications compared to the wild-type (Supplementary Figure S8). A more detailed analysis of the positioning of the FMN in this pocket (Supplementary Figure S9) shows slight modifications of certain distances: Longer average distances between the OH atom of Y143 or the N atom of C144 (vs. G144) and the C7M or C8M atoms of FMN in the G144C mutant compared to the wild-type, longer distances between the CG2 of V191 and C7M, C8M, N3, and N5 atoms of FMN. We can thus conclude that the mutation G144C modifies slightly the positioning of the FMN moiety, which, in return, may impact the redox potential and electron transfer properties of the cofactor. However, the most drastic change can be seen on the electrostatic potential of the surface around the FMN (Figure 5B) which is much less negative than in the wild-type, due to probably both the movement of the E145 sidechain and the presence of the sulfur atom of C144 around. Such a change can also affect the redox potential of the FMN cofactor but also the positioning of the whole FD within the complex with CYPs.
Quite surprisingly, for the G175D mutant (Figure 5C,D), the addition of an acidic residue on β-strand 3 [45] does not modify drastically the structure of this region. The D175 residue is fully embedded in the structure and its modelled pKa value (~9.0) is more than 5 pH units higher than the theoretical value for aspartate residues (3.8), a quite logical value considering the total absence of water molecules surrounding this aspartate. Consequently, the presence of an aspartate does not modify the electrostatic potential around the mutated position. However, the OD1 atom of D175 is at ~3.1 Å of the N atom of T142 strongly suggesting that a hydrogen bond between β-strand 4 and α-helix E may be created. Analysis of the various φ and ψ angles for β-strands 3 and 4 shows little effects of the mutation on the overall secondary structures around the mutation position (Supplementary Figure S10). The average φ angle value of D175 is smaller (Supplementary Figures S10, blue arrows and S11), demonstrating a small but visible distortion of the β-strand 3 (Figure 5D, green for WT, magenta for G175D). This is also reflected in the average φ values for A141 and M140 (Supplementary Figure S10, red arrows) where the changes of the average φ angle value are a bit smaller or bigger (respectively) than for D175. Furthermore, for A141, the deviations around the average φ values are bigger than for the wild-type, suggesting more flexibility at the end of the β-strand 4 (Supplementary Figure S11). These changes impacted the residues surrounding the FMN moiety: The average distances from the OH atom from Y143 and N atom of G144 to the C7M and C8M atoms of FMN were increased (Supplementary Figure S12). Overall, the effects introduced by 3 of the analysed mutations (P117H, G144C, N151D) lead to visible electrostatic surface changes, while, for the G175D mutation, no changes in the electrostatic could be detected but possible changes in the mobility of the residues around the FMN moiety are likely to occur.

3. Discussion

ET between CPR and its partners are mediated by weak and transient complexes that are highly dynamic [5,42] and can be described as an ensemble of various protein orientations rapidly interchanging between themselves [44]. Formation of these complexes is mainly guided by long-range electrostatic interactions between complementary charges on the interacting surfaces [5,14,24,25,31,43,45,46].
Several studies questioned the role of specific residues of the FD of CPR in the complex formation, either using site-specific modifications [3,14,17] or from analysis of naturally occurring variants of CPR (reviewed in [4,32]). In a previous study [16], we demonstrated that selected FD mutations can increase individual CYP isoforms activities. In the present study, we questioned the specificity of four (P117H, G144C, N151D, and G175D) of the previously generated mutations toward CYPs by analyzing the effect of each individual mutation on the three other CYPs and using different substrates. As expected, the majority of those mutations had either no or detrimental effects when tested with other CYP isoforms and can thus be considered CYP specific. It is known that CPR can interact with distinct CYPs and promote ET to several isoforms simultaneously. Our data underpin the notion that these “degenerated” interactions can be modulated by specific positions of the FD, increasing the interactions of CPR with a specific CYP while, at the same time, decreasing the formation of interacting complexes with other CYPs.
However, surprisingly, in the case of the CYP1A2, mutations that had negative effects with EROD and CECOD had neutral (N151D) or enhancing effects (G175D) when DBODF was used as substrate. This implies that the structural determinants present in the proximal side of CYPs can be modulated by the type of substrate, which, in return, may change CPR:CYP interactions, potentially annihilating the specific effects of the FD specific mutations. Moreover, the G175D mutant, selected previously for increased CYP2A6-C7H activity, also improved significantly CYP3A4-DBODF activities. The positive and pleiotropic effect of the G175D mutation was not observed with any of the other mutants but occurred for two of the substrates studied (coumarin and DBF). Therefore, binding of these substrates in different CYP isoforms may provoke similar structural alterations which then lead to better complex formation with the G175D FD mutant.
The inhibitory effect of cyt c on CYP dependent reactions confirmed the existence of coinciding conserved patches of acidic residues of CPR in its interaction with CYP and cyt c. The requirement for relatively high cyt c/CPR ratios to observe this effect indicates that the affinity of CPR for CYP is higher than that of CPR for cyt c. Our data also demonstrated that the amplitude of the cyt c effect is distinct for each CPR variant studied. Furthermore, the inhibition by cyt c is less effective with all CPRmut previously reported to support a gain in activity for the CYP they were isolated when compared with CPRwt. Overall, our results confirm that the strengths of interactions between CYP isoforms and CPR are probably different, as suggested in our previous studies [11,16,18].
Structural effects of the four mutations were analyzed via molecular dynamics simulations. None of them affected drastically the stability or the secondary structures of the FD, in line with the facts that the effects of the mutations were quite specific and the effects measured rather small. Both P117H and N151D affect the electrostatic potential of the surface of the FD, albeit with opposite effects: While the histidine residue in P117H, by being involved in hydrogen bonding with either D116 or E117, slightly neutralizes the negatively charged surface, the presence of an aspartate in N151D, by adding a supplementary negative charge, increases the negative potential of the surface. These conserved acidic patches (D116, E118 and E119) were formerly suggested to play a role in CPR:CYP interactions [14,24,25,43,45,46], while the natural occurring variant A115V was shown to have CYP-isoform dependent effects as well [35,51]. N151D is also in the vicinity of the naturally occurring Q153R variant, a mutation linked with severely impaired steroidogenic activity, which has been shown to have CYP isoform dependent effects as well [33,35,51]. Molecular docking also demonstrated the existence of direct charge interactions between D116 and E118 of CPR with R440 of CYP2D6 [49,50]. Consequently, we hypothesize that the electrostatic potential around D116, E118, and E119 may control the formation of specific contacts between CYPs and CPR and thus, mutations that change the charge distribution modify the binding strength and impacts specifically some CYPs, but not necessarily the others.
In position 144 lays the highly conserved glycine position (100%, among 1221 CPR sequences [16]. This residue is in the vicinity of the conserved acidic patch D147/150, adjacent to the conserved tyrosine residue Y143 considered to be directly involved in FMN binding [26,52]. Additionally, G144 is located close to T142, a position found mutated in a natural occurring variant (T142A), reported to lead to a loss of activity in several CYPs, albeit in an isoform dependent manner [32,51]. Interestingly, while a cysteine residue at this position does not promote major structural rearrangements, the presence of an ionizable sulfur atom induces changes in the electrostatic potential around the FMN moiety. While we do not rule out that such changes also affect the redox potential of the FMN cofactor, the less negative electrostatic potential will also certainly modify the fine positioning of the FD into the CYP binding site prior to ET.
Finally, for the G175D mutant, the aspartate residue, being completely hindered within the structure, is almost non ionizable. However, the aspartate presence induces a small bending of β-strand 3, allowing for smaller φ angle values for A141 (but with larger standard deviations around the average value), impacting the fine positioning of the FMN cofactor.

4. Materials and Methods

4.1. Reagents

l-Arginine, thiamine, chloramphenicol, ampicillin, kanamycin sulfate, isopropyl β-D-thiogalactoside (IPTG) (dioxane-free), δ-aminolevulinic acid, cytochrome c (horse heart), glucose 6-phosphate, glucose 6-phosphate dehydrogenase, nicotinamide adenine dinucleotide phosphate (NADP+ and NADPH), resorufin, 7-hydroxy coumarin (7CH), 3-cyano-7-ethoxycoumarin (CEC), and fluorescein were obtained from Sigma-Aldrich (St. Louis, MO, USA). LB Broth, bacto tryptone and bacto peptone were purchased from BD Biosciences (San Jose, CA, USA). Bacto yeast extract was obtained from Formedium (Norwich, UK). Ethoxyresorufin and coumarin were obtained from BD Biosciences (San Jose, CA, USA) and dibenzylfluorescein from Santa Cruz Biotechnology (Santa Cruz, CA, USA). A polyclonal antibody from rabbit serum raised against recombinant human CPR obtained from Genetex (Irvine, CA, USA) was used for immune-detection of the membrane-bound CPR. All other chemicals and solvents were of the highest grade commercially available.

4.2. Membrane Fractions Preparation

The Escherichia coli cell model BTC was used for the heterologous co-expression of membrane-bound human CYP isoforms (CYP1A2, 2A6 or 3A4) together with human full-length CPR, using a biplasmid co-expression system, as previously reported [11,16,18,22,33,36,37]. CPR FD mutants P117H, G144C, G175D, and N151D, harbored in the pLCM_POR plasmid were used to transform E. coli PD301, already containing the expression vector (pCWori) for human CYP1A2, 2A6 or 3A4 (expressed in N-terminal modified forms), creating the following BTC strains: CPRP117H/CYP2A6, CPRP117H/CYP3A4, CPRG144C/CYP2A6, CPRG144C/CYP3A4, CPRG175D/CYP1A2, CPRG175D/CYP3A4, CPRN151D/CYP1A2, and CPRN151D/CYP2A6. Strains were cultured for heterologous expression of the human proteins [33,37,38], and membrane fractions of the different strains were prepared and characterized for total protein content (Bradford assay) as well as CYP- (CO-difference spectrophotometry) and CPR- (immune-detection by Western blotting) contents, as described previously [11,22,26,33,36,37,38].

4.3. CYP Enzyme Activity Assays

The catalytic activities of CYP1A2 (EROD), 2A6 (C7H), or 3A4 (DBODF), sustained by FD mutants were evaluated as reported previously, using 0–5 µM EthR (CYP1A2), 0–20 µM coumarin (CYP2A6), and 0–10 µM DBF (CYP3A4) [11,16]. When testing the effects of FD mutants in CYP1A2-activities (CECOD and DBODF), the concentrations ranges were 0–50 µM CEC and 0–37.5 µM DBF. Initial velocities were measured in triplicate and plots of velocity traces vs. substrates concentrations could be fitted according to the Michaelis–Menten equation (r2 ≥ 0.95) to determine steady-state kinetic parameters, using GraphPad Prism 5.01 Software (La Jolla, CA, USA) [22,33,38].

4.4. Cytochrome c Competition Assays

The effect of cyt c on CYP1A2, 2A6, or 3A4 activities, sustained by CPRwt and CPRmut, was assessed using the same enzyme assay conditions as described above, except substrate concentrations were hold constant (5 µM EthR, 20 µM coumarin or 7.5 µM DBF), and a gradient of cyt c (0–40 µM for CYP1A2 and 3A4, 0–50 µM for CYP2A6) was applied in 100 mM potassium phosphate buffer (pH 7.6), and NADPH regenerating system (NADPH 200 µM, glucose 6-phosphate 500 µM and glucose 6-phosphate dehydrogenase 40 U L−1, final concentrations). Velocities were measured in triplicate in 96-well format using multi-mode microtiter plate reader (SpectraMax®i3x, Molecular Devices, USA; SoftMax Pro 2.0) [11,16]. All reaction velocities (kobs) were measured at least in triplicate [pmol of fluorescent product formed/(pmol of CYP per minute)].

4.5. Structural Analysis of CPR Mutants

The structural analysis of the four mutants (P117H, G144C, N151D and G175D) was performed using the CPR FD only, obtained from the crystal structure of human CPRwt (PDB 5FA6) as the starting material. Mutations were constructed with the YASARA software (www.yasara.org). Fully hydrated molecular dynamics simulations were performed for 80 ns using a step size of 100 ps with the Yasara software suite using the AMBER14 force field at constant pressure (water density at 0.997 g/mL) and temperature (298 K). Parameters were set as follows: Cuboid cell extending 10 nm on each side of the protein, 0.9% NaCl and pH = 7.4. The Bio3D suite in R was used to calculate the coordinates of the average structure from all snapshots, the root mean square deviations (RMSD), the various reported distances and φ or ψ angles. Averages and standard deviations of the data was calculated with R. pKa of ionizable residues were calculated from the average structures using the AMBER force-field and the propka option of the pdb2pqr software at a pH of 7.5. Electrostatic potentials were calculated on the average structures with the integrated module of the Pymol software (www.pymol.org) using default parameters (non-linear Poisson-Boltzmann equation, protein dielectric 2.0, solvent dielectric 78.0, temperature 310 K, ion charge +1 of radius 2.0 Å and −1 of radius 1.8 both at 0.15 M) and plotted on the solvent excluded surfaces (Connolly surfaces). Figures were prepared with the Pymol software and R ggplot function. All modelling data is available upon request.

4.6. Statistical Analysis

Variance in data was analyzed through one-way ANOVA with Bonferroni’s multiple comparison test. The unpaired Student t-test was performed for calculation of the two-tailed p value. The analysis was performed with 95% confidence interval using the GraphPad Prism 5.01 Software. The significance level considered in all the statistical tests was 0.05.

5. Conclusions

Our data confirm and add additional evidences that CPR–CYP interactions are degenerated. While the purpose of this degeneracy is evident as only one donor serves multiple partners (and sometimes at the same time), it does not preclude the possibility of a gain of specificity, as demonstrated by the effects of our selected FD mutations. It is noteworthy that both surfaces (CPR and CYP) can be modified in order to get stronger interactions: Via mutations on the CPR side and via mutation and/or substrate binding on the CYP side. This latter effect may well explain how, in situation when specific CYP isoforms are metabolizing actively foreign or endogenous substrates, CPR could actually distinguish between non active and active CYPs, thus rendering CYP metabolism more efficient without any other transcriptional or posttranscriptional control, an elegant and efficient way to adjust CYP oxidative activities.

Supplementary Materials

Supplementary Materials can be found at https://www.mdpi.com/1422-0067/21/18/6669/s1. Table S1: Kinetic steady-state parameters of CYP1A2, 2A6 and 3A4-activities; Table S2: Kinetic steady-state parameters of CYP1A2-activities with CEC and DBF; Table S3: Summary of determinant CPR residues directly involved in hydrogen bonds/salt bridges on both FD:cyt c and FD:CYP binding sites; Figure S1: Protein sequence the human CPR FD; Figure S2: Relative turnover rates (kcat) (x fold) of the CPRmut:CYP1A2 MROD activities (using 0–2 μM MthR), normalized by the kcat demonstrated by CPRwt/CYP (determined in membrane fractions) (technical replicates N = 3); Figure S3: Structural models of the FD binding site of rat CPR; Figure S4: RMSD of the 80 ns simulations; Figure S5: Comparison of the deviations between each average structure obtained from simulation and the FMN domain from the crystal structure of human CPR (PDB 5FA6); Figure S6: Hydrogen bonds between H117 and D116 or E118 in the P117H mutant; Figure S7: Acidic residues surrounding the aspartate 151 in the N151 D mutant; Figure S8: Variation of pKa of the surrounding charged residues for G144C mutant; Figure S9: Distances between the OH atom of Y143 or Y181, N atom of C144 or G144, CG2 atom of V191 and either the C7M, C8M, N3 or N5 atoms of FMN in WT or G144C mutant; Figure S10: Analyses of the various φ and ψ angles for β-strands 3 and 4; Figure S11: Analysis of the φ angles along the four simulations for residues 175 and 141 in the WT and G175 FD mutant; Figure S12: Distances between the OH atom of Y143 or Y181, N atom of C144 or G144, CG2 atom of V191 and either the C7M, C8M, N3 or N5 atoms of FMN in WT or G175D mutant.

Author Contributions

Conceptualization, F.E., G.T. and M.K.; methodology, F.E. and G.T.; formal analysis, F.E., G.T. and M.K.; writing—original draft preparation, F.E., G.T. and M.K.; writing—review and editing, all authors; project administration and supervision, G.T. and M.K.; funding acquisition, M.K. and G.T. All authors have read and agreed to the published version of the manuscript.

Funding

This research was in part funded by a joint project, funded by the Portuguese Fundação para a Ciência e a Tecnologia [Grant FCT-ANR/BEX-BCM/0002/2013 as well as Grant UID/BIM/0009/2016 of the Research Center for Toxicogenomics and Human Health (ToXomics)] and the L’Agence Nationale de la Recherche (Grant ANR-13-ISV5-0001).

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

C7HCoumarin 7-hydroxylation
CEC3-cyano-7-ethoxycoumarin
CECOD3-cyano-7-ethoxycoumarin O-dealkylation
CPRNADPH cytochrome P450 oxidoreductase
CYPCytochrome P450
Cytochrome cCyt c
DBFDibenzylfluorescein
DBODFDibenzylfluorescein O-debenzylation
ERODEthoxyresorufin O-deethylation
ETElectron transfer
EthREthoxyresorufin
FADFlavin adenine dinucleotide
FDFMN binding domain
FMNFlavin mononucleotide

References

  1. Guengerich, F.P. Human Cytochrome P450 Enzymes. In Cytochrome P450: Structure, Mechanism, and Biochemistry, 4th ed.; Ortiz de Montellano, P.R., Ed.; Springer International Publishing: Cham, Switzerland, 2015; pp. 523–785. ISBN 978-3-319-12108-6. [Google Scholar]
  2. Hlavica, P. Mechanistic basis of electron transfer to cytochromes P450 by natural redox partners and artificial donor constructs. Adv. Exp. Med. Biol. 2015, 851, 247–297. [Google Scholar] [CrossRef] [PubMed]
  3. Strobel, H.W.; Hodgson, A.V.; Shen, S. NADPH Cytochrome P450 Reductase and Its Structural and Functional Domains. In Cytochrome P450: Structure, Mechanism, and Biochemistry, 2nd ed.; Ortiz de Montellano, P.R.O., Ed.; Springer: Boston, MA, USA, 1995; pp. 225–244. ISBN 978-1-4757-2391-5. [Google Scholar]
  4. Pandey, A.V.; Flück, C.E. NADPH P450 oxidoreductase: Structure, function, and pathology of diseases. Pharmacol. Ther. 2013, 138, 229–254. [Google Scholar] [CrossRef] [PubMed]
  5. Waskell, L.; Kim, J.-J.P. Electron Transfer Partners of Cytochrome P450. In Cytochrome P450: Structure, Mechanism, and Biochemistry, 4th ed.; Ortiz de Montellano, P.R., Ed.; Springer International Publishing: Cham, Switzerland, 2015; pp. 33–68. ISBN 978-3-319-12108-6. [Google Scholar]
  6. Yablokov, E.O.; Sushko, T.A.; Ershov, P.V.; Florinskaya, A.V.; Gnedenko, O.V.; Shkel, T.V.; Grabovec, I.P.; Strushkevich, N.V.; Kaluzhskiy, L.A.; Usanov, S.A.; et al. A large-scale comparative analysis of affinity, thermodynamics and functional characteristics of interactions of twelve cytochrome P450 isoforms and their redox partners. Biochimie 2019, 162, 156–166. [Google Scholar] [CrossRef] [PubMed]
  7. Pikuleva, I.A.; Cao, C.; Waterman, M.R. An additional electrostatic interaction between adrenodoxin and P450c27 (CYP27A1) results in tighter binding than between adrenodoxin and p450scc (CYP11A1). J. Biol. Chem. 1999, 274, 2045–2052. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Gutierrez, A.; Paine, M.; Wolf, C.R.; Scrutton, N.S.; Roberts, G.C.K. Relaxation kinetics of cytochrome P450 reductase: Internal electron transfer is limited by conformational change and regulated by coenzyme binding. Biochemistry 2002, 41, 4626–4637. [Google Scholar] [CrossRef]
  9. Hlavica, P.; Lehnerer, M. Oxidative biotransformation of fatty acids by cytochromes P450: Predicted key structural elements orchestrating substrate specificity, regioselectivity and catalytic efficiency. Curr. Drug Metab. 2010, 11, 85–104. [Google Scholar] [CrossRef]
  10. Kandel, S.E.; Lampe, J.N. Role of protein-protein interactions in cytochrome P450-mediated drug metabolism and toxicity. Chem. Res. Toxicol. 2014, 27, 1474–1486. [Google Scholar] [CrossRef] [Green Version]
  11. Campelo, D.; Esteves, F.; Brito Palma, B.; Costa Gomes, B.; Rueff, J.; Lautier, T.; Urban, P.; Truan, G.; Kranendonk, M. Probing the Role of the Hinge Segment of Cytochrome P450 Oxidoreductase in the Interaction with Cytochrome P450. Int. J. Mol. Sci. 2018, 19, 3914. [Google Scholar] [CrossRef] [Green Version]
  12. Aigrain, L.; Pompon, D.; Moréra, S.; Truan, G. Structure of the open conformation of a functional chimeric NADPH cytochrome P450 reductase. EMBO Rep. 2009, 10, 742–747. [Google Scholar] [CrossRef] [Green Version]
  13. Ellis, J.; Gutierrez, A.; Barsukov, I.L.; Huang, W.-C.; Grossmann, J.G.; Roberts, G.C.K. Domain motion in cytochrome P450 reductase: Conformational equilibria revealed by NMR and small-angle x-ray scattering. J. Biol. Chem. 2009, 284, 36628–36637. [Google Scholar] [CrossRef] [Green Version]
  14. Hamdane, D.; Xia, C.; Im, S.-C.; Zhang, H.; Kim, J.-J.P.; Waskell, L. Structure and function of an NADPH-cytochrome P450 oxidoreductase in an open conformation capable of reducing cytochrome P450. J. Biol. Chem. 2009, 284, 11374–11384. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Frances, O.; Fatemi, F.; Pompon, D.; Guittet, E.; Sizun, C.; Pérez, J.; Lescop, E.; Truan, G. A well-balanced preexisting equilibrium governs electron flux efficiency of a multidomain diflavin reductase. Biophys. J. 2015, 108, 1527–1536. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Esteves, F.; Campelo, D.; Gomes, B.C.; Urban, P.; Bozonnet, S.; Lautier, T.; Rueff, J.; Truan, G.; Kranendonk, M. The Role of the FMN-Domain of Human Cytochrome P450 Oxidoreductase in Its Promiscuous Interactions with Structurally Diverse Redox Partners. Front. Pharmacol. 2020, 11, 299. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Shen, A.L.; Kasper, C.B. Role of acidic residues in the interaction of NADPH-cytochrome P450 oxidoreductase with cytochrome P450 and cytochrome c. J. Biol. Chem. 1995, 270, 27475–27480. [Google Scholar] [CrossRef] [Green Version]
  18. Campelo, D.; Lautier, T.; Urban, P.; Esteves, F.; Bozonnet, S.; Truan, G.; Kranendonk, M. The Hinge Segment of Human NADPH-Cytochrome P450 Reductase in Conformational Switching: The Critical Role of Ionic Strength. Front. Pharmacol. 2017, 8, 755. [Google Scholar] [CrossRef]
  19. Barnaba, C.; Gentry, K.; Sumangala, N.; Ramamoorthy, A. The catalytic function of cytochrome P450 is entwined with its membrane-bound nature. F1000Research 2017, 6, 662. [Google Scholar] [CrossRef]
  20. Barnaba, C.; Ramamoorthy, A. Picturing the Membrane-assisted Choreography of Cytochrome P450 with Lipid Nanodiscs. ChemPhysChem 2018, 19, 2603–2613. [Google Scholar] [CrossRef]
  21. Denisov, I.G.; Sligar, S.G. Nanodiscs in Membrane Biochemistry and Biophysics. Chem. Rev. 2017, 117, 4669–4713. [Google Scholar] [CrossRef]
  22. Kranendonk, M.; Marohnic, C.C.; Panda, S.P.; Duarte, M.P.; Oliveira, J.S.; Masters, B.S.S.; Rueff, J. Impairment of human CYP1A2-mediated xenobiotic metabolism by Antley-Bixler syndrome variants of cytochrome P450 oxidoreductase. Arch. Biochem. Biophys. 2008, 475, 93–99. [Google Scholar] [CrossRef] [Green Version]
  23. Backes, W.L.; Kelley, R.W. Organization of multiple cytochrome P450s with NADPH-cytochrome P450 reductase in membranes. Pharmacol. Ther. 2003, 98, 221–233. [Google Scholar] [CrossRef]
  24. Shen, S.; Strobel, H.W. Probing the putative cytochrome P450- and cytochrome c-binding sites on NADPH-cytochrome P450 reductase by anti-peptide antibodies. Biochemistry 1994, 33, 8807–8812. [Google Scholar] [CrossRef] [PubMed]
  25. Zhao, Q.; Modi, S.; Smith, G.; Paine, M.; McDonagh, P.D.; Wolf, C.R.; Tew, D.; Lian, L.Y.; Roberts, G.C.; Driessen, H.P. Crystal structure of the FMN-binding domain of human cytochrome P450 reductase at 1.93 A resolution. Protein Sci. 1999, 8, 298–306. [Google Scholar] [CrossRef] [Green Version]
  26. Marohnic, C.C.; Panda, S.P.; McCammon, K.; Rueff, J.; Masters, B.S.S.; Kranendonk, M. Human cytochrome P450 oxidoreductase deficiency caused by the Y181D mutation: Molecular consequences and rescue of defect. Drug Metab. Dispos. Biol. Fate Chem. 2010, 38, 332–340. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Reed, J.R.; Backes, W.L. The functional effects of physical interactions involving cytochromes P450: Putative mechanisms of action and the extent of these effects in biological membranes. Drug Metab. Rev. 2016, 48, 453–469. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Nair, P.C.; McKinnon, R.A.; Miners, J.O. Cytochrome P450 structure-function: Insights from molecular dynamics simulations. Drug Metab. Rev. 2016, 48, 434–452. [Google Scholar] [CrossRef]
  29. Gentry, K.A.; Anantharamaiah, G.M.; Ramamoorthy, A. Probing protein-protein and protein-substrate interactions in the dynamic membrane-associated ternary complex of cytochromes P450, b5, and reductase. Chem. Commun. 2019, 55, 13422–13425. [Google Scholar] [CrossRef] [PubMed]
  30. Liou, S.-H.; Chuo, S.-W.; Qiu, Y.; Wang, L.-P.; Goodin, D.B. Linkage between Proximal and Distal Movements of P450cam Induced by Putidaredoxin. Biochemistry 2020, 59, 2012–2021. [Google Scholar] [CrossRef] [PubMed]
  31. Šrejber, M.; Navrátilová, V.; Paloncýová, M.; Bazgier, V.; Berka, K.; Anzenbacher, P.; Otyepka, M. Membrane-attached mammalian cytochromes P450: An overview of the membrane’s effects on structure, drug binding, and interactions with redox partners. J. Inorg. Biochem. 2018, 183, 117–136. [Google Scholar] [CrossRef]
  32. Burkhard, F.Z.; Parween, S.; Udhane, S.S.; Flück, C.E.; Pandey, A.V. P450 Oxidoreductase deficiency: Analysis of mutations and polymorphisms. J. Steroid Biochem. Mol. Biol. 2017, 165, 38–50. [Google Scholar] [CrossRef] [Green Version]
  33. Moutinho, D.; Marohnic, C.C.; Panda, S.P.; Rueff, J.; Masters, B.S.; Kranendonk, M. Altered human CYP3A4 activity caused by Antley-Bixler syndrome-related variants of NADPH-cytochrome P450 oxidoreductase measured in a robust in vitro system. Drug Metab. Dispos. Biol. Fate Chem. 2012, 40, 754–760. [Google Scholar] [CrossRef] [Green Version]
  34. Huang, N.; Agrawal, V.; Giacomini, K.M.; Miller, W.L. Genetics of P450 oxidoreductase: Sequence variation in 842 individuals of four ethnicities and activities of 15 missense mutations. Proc. Natl. Acad. Sci. USA 2008, 105, 1733–1738. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Udhane, S.S.; Parween, S.; Kagawa, N.; Pandey, A.V. Altered CYP19A1 and CYP3A4 Activities Due to Mutations A115V, T142A, Q153R and P284L in the Human P450 Oxidoreductase. Front. Pharmacol. 2017, 8, 580. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Duarte, M.P.; Palma, B.B.; Gilep, A.A.; Laires, A.; Oliveira, J.S.; Usanov, S.A.; Rueff, J.; Kranendonk, M. The stimulatory role of human cytochrome b5 in the bioactivation activities of human CYP1A2, 2A6 and 2E1: A new cell expression system to study cytochrome P450-mediated biotransformation (a corrigendum report on Duarte et al. (2005) Mutagenesis 20, 93–100). Mutagenesis 2007, 22, 75–81. [Google Scholar] [CrossRef] [Green Version]
  37. Palma, B.B.; Silva E Sousa, M.; Urban, P.; Rueff, J.; Kranendonk, M. Functional characterization of eight human CYP1A2 variants: The role of cytochrome b5. Pharmacogenet. Genom. 2013, 23, 41–52. [Google Scholar] [CrossRef] [PubMed]
  38. Esteves, F.; Campelo, D.; Urban, P.; Bozonnet, S.; Lautier, T.; Rueff, J.; Truan, G.; Kranendonk, M. Human cytochrome P450 expression in bacteria: Whole-cell high-throughput activity assay for CYP1A2, 2A6 and 3A4. Biochem. Pharmacol. 2018, 158, 134–140. [Google Scholar] [CrossRef] [PubMed]
  39. Duarte, M.P.; Palma, B.B.; Laires, A.; Oliveira, J.S.; Rueff, J.; Kranendonk, M. Escherichia coli BTC, a human cytochrome P450 competent tester strain with a high sensitivity towards alkylating agents: Involvement of alkyltransferases in the repair of DNA damage induced by aromatic amines. Mutagenesis 2005, 20, 199–208. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  40. Paine, M.F.; Khalighi, M.; Fisher, J.M.; Shen, D.D.; Kunze, K.L.; Marsh, C.L.; Perkins, J.D.; Thummel, K.E. Characterization of interintestinal and intraintestinal variations in human CYP3A-dependent metabolism. J. Pharmacol. Exp. Ther. 1997, 283, 1552–1562. [Google Scholar]
  41. Venkatakrishnan, K.; von Moltke, L.L.; Court, M.H.; Harmatz, J.S.; Crespi, C.L.; Greenblatt, D.J. Comparison between cytochrome P450 (CYP) content and relative activity approaches to scaling from cDNA-expressed CYPs to human liver microsomes: Ratios of accessory proteins as sources of discrepancies between the approaches. Drug Metab. Dispos. 2000, 28, 1493–1504. [Google Scholar]
  42. Scott, E.E.; Wolf, C.R.; Otyepka, M.; Humphreys, S.C.; Reed, J.R.; Henderson, C.J.; McLaughlin, L.A.; Paloncýová, M.; Navrátilová, V.; Berka, K.; et al. The Role of Protein-Protein and Protein-Membrane Interactions on P450 Function. Drug Metab. Dispos. Biol. Fate Chem. 2016, 44, 576–590. [Google Scholar] [CrossRef] [Green Version]
  43. Jang, H.-H.; Jamakhandi, A.P.; Sullivan, S.Z.; Yun, C.-H.; Hollenberg, P.F.; Miller, G.P. Beta sheet 2-alpha helix C loop of cytochrome P450 reductase serves as a docking site for redox partners. Biochim. Biophys. Acta 2010, 1804, 1285–1293. [Google Scholar] [CrossRef] [Green Version]
  44. Huang, R.; Zhang, M.; Rwere, F.; Waskell, L.; Ramamoorthy, A. Kinetic and structural characterization of the interaction between the FMN binding domain of cytochrome P450 reductase and cytochrome c. J. Biol. Chem. 2015, 290, 4843–4855. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Wang, M.; Roberts, D.L.; Paschke, R.; Shea, T.M.; Masters, B.S.; Kim, J.J. Three-dimensional structure of NADPH-cytochrome P450 reductase: Prototype for FMN- and FAD-containing enzymes. Proc. Natl. Acad. Sci. USA 1997, 94, 8411–8416. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Hong, Y.; Li, H.; Yuan, Y.-C.; Chen, S. Sequence-function correlation of aromatase and its interaction with reductase. J. Steroid Biochem. Mol. Biol. 2010, 118, 203–206. [Google Scholar] [CrossRef] [Green Version]
  47. Nicolo, C.; Flück, C.E.; Mullis, P.E.; Pandey, A.V. Restoration of mutant cytochrome P450 reductase activity by external flavin. Mol. Cell. Endocrinol. 2010, 321, 245–252. [Google Scholar] [CrossRef] [PubMed]
  48. Nisimoto, Y. Localization of cytochrome c-binding domain on NADPH-cytochrome P-450 reductase. J. Biol. Chem. 1986, 261, 14232–14239. [Google Scholar] [PubMed]
  49. Allorge, D.; Bréant, D.; Harlow, J.; Chowdry, J.; Lo-Guidice, J.-M.; Chevalier, D.; Cauffiez, C.; Lhermitte, M.; Blaney, F.E.; Tucker, G.T.; et al. Functional analysis of CYP2D6.31 variant: Homology modeling suggests possible disruption of redox partner interaction by Arg440His substitution. Proteins 2005, 59, 339–346. [Google Scholar] [CrossRef]
  50. Sündermann, A.; Oostenbrink, C. Molecular dynamics simulations give insight into the conformational change, complex formation, and electron transfer pathway for cytochrome P450 reductase. Protein Sci. Publ. Protein Soc. 2013, 22, 1183–1195. [Google Scholar] [CrossRef] [Green Version]
  51. Huang, N.; Pandey, A.V.; Agrawal, V.; Reardon, W.; Lapunzina, P.D.; Mowat, D.; Jabs, E.W.; Van Vliet, G.; Sack, J.; Flück, C.E.; et al. Diversity and function of mutations in p450 oxidoreductase in patients with Antley-Bixler syndrome and disordered steroidogenesis. Am. J. Hum. Genet. 2005, 76, 729–749. [Google Scholar] [CrossRef] [Green Version]
  52. Xia, C.; Panda, S.P.; Marohnic, C.C.; Martásek, P.; Masters, B.S.; Kim, J.-J.P. Structural basis for human NADPH-cytochrome P450 oxidoreductase deficiency. Proc. Natl. Acad. Sci. USA 2011, 108, 13486–13491. [Google Scholar] [CrossRef] [Green Version]
Ijms 21 06669 g001 550
Figure 1. Relative turnover rates (kcat) (x fold) of the CPRmut:CYP1A2 EROD (A), CPRmut:CYP2A6 C7H (B), and CPRmut:CYP3A4 DBODF (C) activities, normalized by the kcat demonstrated by CPRwt/CYP (determined in membrane fractions; technical replicates N = 3). CPR-FMN-domain wildtype (CPRwt:CYP, black stripes). CPR-FMN-domain mutants previously demonstrated to support a gain in CYP1A2-mediated EROD- (P117H and G144C: green), CYP2A6-mediated C7H- (G175D: blue), and CYP3A4-mediated DBODF-activity (N151D: red). kcat values of the CPRmut/CYP were compared with the ones of the CPRwt/CYP applying the unpaired t test (*** p < 0.0005; ** p < 0.005; * p < 0.05).
Figure 1. Relative turnover rates (kcat) (x fold) of the CPRmut:CYP1A2 EROD (A), CPRmut:CYP2A6 C7H (B), and CPRmut:CYP3A4 DBODF (C) activities, normalized by the kcat demonstrated by CPRwt/CYP (determined in membrane fractions; technical replicates N = 3). CPR-FMN-domain wildtype (CPRwt:CYP, black stripes). CPR-FMN-domain mutants previously demonstrated to support a gain in CYP1A2-mediated EROD- (P117H and G144C: green), CYP2A6-mediated C7H- (G175D: blue), and CYP3A4-mediated DBODF-activity (N151D: red). kcat values of the CPRmut/CYP were compared with the ones of the CPRwt/CYP applying the unpaired t test (*** p < 0.0005; ** p < 0.005; * p < 0.05).
Ijms 21 06669 g001
Ijms 21 06669 g002 550
Figure 2. Relative turnover rates (kcat) (x fold) of the CPRmut:CYP1A2 CECOD (A), and DBODF (B) activities, normalized by the kcat demonstrated by CPRwt/CYP (determined in membrane fractions) (technical replicates N = 3). CPR-FMN-domain wildtype (CPRwt:CYP, black stripes). CPR-FMN-domain mutants previously demonstrated to support a gain in CYP1A2-mediated EROD- (P117H and G144C, green), CYP2A6-mediated C7H- (G175D, blue), and CYP3A4-mediated DBODF-activity (N151D, red). kcat values of the CPRmut/CYP1A2 were compared with the ones of the CPRwt/CYP1A2 applying the unpaired t test (** p < 0.005; * p < 0.05).
Figure 2. Relative turnover rates (kcat) (x fold) of the CPRmut:CYP1A2 CECOD (A), and DBODF (B) activities, normalized by the kcat demonstrated by CPRwt/CYP (determined in membrane fractions) (technical replicates N = 3). CPR-FMN-domain wildtype (CPRwt:CYP, black stripes). CPR-FMN-domain mutants previously demonstrated to support a gain in CYP1A2-mediated EROD- (P117H and G144C, green), CYP2A6-mediated C7H- (G175D, blue), and CYP3A4-mediated DBODF-activity (N151D, red). kcat values of the CPRmut/CYP1A2 were compared with the ones of the CPRwt/CYP1A2 applying the unpaired t test (** p < 0.005; * p < 0.05).
Ijms 21 06669 g002
Ijms 21 06669 g003a 550Ijms 21 06669 g003b 550
Figure 3. Relative CYP activities of the CPR-FMN-domain wildtype and four mutants plotted in function of the cyt c/CPR ratio. (A) CYP1A2 (EROD), (B) CYP2A6 (C7H), and (C) CYP3A4 (DBODF). Observed rate constants sustained by the CPR (kobs) were normalized by the CYP activity without cyt c [kobs (0 nM cyt c)]. FMN wildtype (CPRwt:CYP, black lines). FMN mutants previously demonstrated to support a gain in CYP1A2-mediated EROD- (P117H and G144C, green lines), CYP2A6-mediated C7H- (G175D, blue lines), and CYP3A4-mediated DBODF-activity (N151D, red lines). Normalized CYP activities represent the average of the three replicates and the error bars the standard deviation. Significant differences between the means of the activity profiles of CPR:CYP1A2 EROD (p = 0.015) (A) and CPR:2A6 C7H (p = 0.049) (B) were observed (ANOVA, one-way analysis of variance). (D) CPRwt supporting CYP1A2-mediated EROD- (green line), CYP2A6-mediated C7H- (blue line), and CYP3A4-mediated DBODF-activity (red line).
Figure 3. Relative CYP activities of the CPR-FMN-domain wildtype and four mutants plotted in function of the cyt c/CPR ratio. (A) CYP1A2 (EROD), (B) CYP2A6 (C7H), and (C) CYP3A4 (DBODF). Observed rate constants sustained by the CPR (kobs) were normalized by the CYP activity without cyt c [kobs (0 nM cyt c)]. FMN wildtype (CPRwt:CYP, black lines). FMN mutants previously demonstrated to support a gain in CYP1A2-mediated EROD- (P117H and G144C, green lines), CYP2A6-mediated C7H- (G175D, blue lines), and CYP3A4-mediated DBODF-activity (N151D, red lines). Normalized CYP activities represent the average of the three replicates and the error bars the standard deviation. Significant differences between the means of the activity profiles of CPR:CYP1A2 EROD (p = 0.015) (A) and CPR:2A6 C7H (p = 0.049) (B) were observed (ANOVA, one-way analysis of variance). (D) CPRwt supporting CYP1A2-mediated EROD- (green line), CYP2A6-mediated C7H- (blue line), and CYP3A4-mediated DBODF-activity (red line).
Ijms 21 06669 g003aIjms 21 06669 g003b
Ijms 21 06669 g004 550
Figure 4. Electrostatics effects of the P117H and N151D mutations. Panels (A,C) present the variations of pKa of the surrounding aspartates and glutamates of P117H and N151D, respectively. pKas values were calculated on the average structures generated with four independent molecular dynamics simulations for wild-type (WT), P117H (P117H), or N151D (N151D) with the propka options of the pdb2pqr software, using a pH = 7.5 value and the AMBER force field. Theoretical pKa values of aspartate (3.8) or glutamate (4.5) were then subtracted from the simulated ones. Panels (B,D) present a mapping of the electrostatic potential onto the solvent excluded surface of P117H (Panel (B)) or N151D (Panel (D)) mutants compared to the wildtype FD. In Panel (B), the arrow represents the position of the H117 residue (hidden by the surface). In Panel (D), the magenta surface corresponds to the D151 residue. In Panels (B,D), the light orange oval shows the locations of the electrostatic potential and surface changes. Electrostatic potentials were calculated on the average structure of simulation 1 with the APBS module of Pymol using defaults parameters. The color scale is indicated below each surface and is in kT/e.
Figure 4. Electrostatics effects of the P117H and N151D mutations. Panels (A,C) present the variations of pKa of the surrounding aspartates and glutamates of P117H and N151D, respectively. pKas values were calculated on the average structures generated with four independent molecular dynamics simulations for wild-type (WT), P117H (P117H), or N151D (N151D) with the propka options of the pdb2pqr software, using a pH = 7.5 value and the AMBER force field. Theoretical pKa values of aspartate (3.8) or glutamate (4.5) were then subtracted from the simulated ones. Panels (B,D) present a mapping of the electrostatic potential onto the solvent excluded surface of P117H (Panel (B)) or N151D (Panel (D)) mutants compared to the wildtype FD. In Panel (B), the arrow represents the position of the H117 residue (hidden by the surface). In Panel (D), the magenta surface corresponds to the D151 residue. In Panels (B,D), the light orange oval shows the locations of the electrostatic potential and surface changes. Electrostatic potentials were calculated on the average structure of simulation 1 with the APBS module of Pymol using defaults parameters. The color scale is indicated below each surface and is in kT/e.
Ijms 21 06669 g004
Ijms 21 06669 g005 550
Figure 5. Effects of the G144C and G175D mutations. Panel (A) presents the wild-type (green, simulation 1) and G144C (magenta, simulation 1) average structures. Panel (B) presents a mapping of the electrostatic potential onto the solvent excluded surface of wild-type (left) or G144C mutant (right), the orange color corresponds to the surface of the flavin mononucleotide (FMN) residue. Panel (C,D) presents detailed view of the wild-type (green, simulation 1) and G175D (magenta, simulation 1) average structures.
Figure 5. Effects of the G144C and G175D mutations. Panel (A) presents the wild-type (green, simulation 1) and G144C (magenta, simulation 1) average structures. Panel (B) presents a mapping of the electrostatic potential onto the solvent excluded surface of wild-type (left) or G144C mutant (right), the orange color corresponds to the surface of the flavin mononucleotide (FMN) residue. Panel (C,D) presents detailed view of the wild-type (green, simulation 1) and G175D (magenta, simulation 1) average structures.
Ijms 21 06669 g005
Table
Table 1. Cytochromes P450 (CYP) and cytochrome P450 oxidoreductase (CPR) contents of bacterial CPR/CYP co-expression system (BTC) membrane fractions.
Table 1. Cytochromes P450 (CYP) and cytochrome P450 oxidoreductase (CPR) contents of bacterial CPR/CYP co-expression system (BTC) membrane fractions.
Membrane FractionsProtein Contents
Cyp
Isoform		CPR
Form		CYPCPRCPR/CYP
(pmol/mg Protein)Ratios a
CYP1A2wt a54 ± 14.1 ± 1.51:13
P117H b148 ± 112.2 ± 0.51:12
G144C b124 ± 214.1 ± 0.21:9
G175D379 ± 226.5 ± 0.41:14
N151D446 ± 427.7 ± 0.51:16
CPRnull b91 ± 6NANA
CYP2A6wt b139 ± 110.5 ± 1.31:13
P117H205 ± 316.5 ± 0.91:12
G144C315 ± 828.3 ± 0.81:11
G175D b112 ± 611.3 ± 0.51:10
N151D221 ± 119.7 ± 0.81:11
CPRnull b178 ± 4NANA
CYP3A4wt a83 ± 319.8 ± 0.21:4
P117H230 ± 935.5 ± 0.81:6
G144C189 ± 1126.5 ± 0.31:7
G175D143 ± 531.5 ± 0.51:5
N151D b45 ± 112.7 ± 0.61:4
CPRnull b85 ± 3NANA
CYP and CPR contents are mean ± sd (technical replicates N = 3). Color code indicates the CYP isoform with which the CPR-FMN-domain mutant was isolated (CYP1A2 green; CYP2A6: blue; CYP3A4: red). (NA) not applicable. a Co-expression of CPRmut with CYPs was achieved with stoichiometry’s approximating those obtained with CPRwt. CPR/CYP ratios of the CPRmut were not significantly different from the ones of the corresponding CPRwt (p > 0.05), with each tested CYP. b Membrane fractions and values from batches used in our former studies [11,16].

Share and Cite

MDPI and ACS Style

Esteves, F.; Urban, P.; Rueff, J.; Truan, G.; Kranendonk, M. Interaction Modes of Microsomal Cytochrome P450s with Its Reductase and the Role of Substrate Binding. Int. J. Mol. Sci. 2020, 21, 6669. https://doi.org/10.3390/ijms21186669

AMA Style

Esteves F, Urban P, Rueff J, Truan G, Kranendonk M. Interaction Modes of Microsomal Cytochrome P450s with Its Reductase and the Role of Substrate Binding. International Journal of Molecular Sciences. 2020; 21(18):6669. https://doi.org/10.3390/ijms21186669

Chicago/Turabian Style

Esteves, Francisco, Philippe Urban, José Rueff, Gilles Truan, and Michel Kranendonk. 2020. "Interaction Modes of Microsomal Cytochrome P450s with Its Reductase and the Role of Substrate Binding" International Journal of Molecular Sciences 21, no. 18: 6669. https://doi.org/10.3390/ijms21186669

APA Style

Esteves, F., Urban, P., Rueff, J., Truan, G., & Kranendonk, M. (2020). Interaction Modes of Microsomal Cytochrome P450s with Its Reductase and the Role of Substrate Binding. International Journal of Molecular Sciences, 21(18), 6669. https://doi.org/10.3390/ijms21186669

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop