Structural and functional characterization of a cytochrome P450 2B4 F429H mutant with an axial thiolate-histidine hydrogen bond

Mixed H2S and H2O Clusters─New Insights into Dispersion-Dominated Hydrogen Bonding

Here, we report the results of an IR spectroscopy study on heteroclusters of H2S and H2O and several of their isotopomers using mass-selective IR spectroscopy in superfluid helium nanodroplets in the range of 2560-2800 cm-1. Based on DFT calculations on the B3LYP-D3/6-311++G(d,p) level of theory, we were able to assign the experimentally observed O-D stretching bands to heterodimer and heterotrimer clusters. Since no bands of the S-H-bound conformer HSH···OH2 could be observed, we were able to determine the O-H-bound conformer HOH···SH2 to be the global minimum structure. A trapping of the local minima in helium nanodroplets was not observed. This is in line with the weaker hydrogen bond expected for H2S complexes. In these clusters, the interaction energy is expected to be more dominated by dispersion and less dictated by highly directional electrostatic forces.

Effect of CO binding to P450 BM3 F393 mutants on electron density distribution in the heme cofactor

Resonance Raman spectroscopy has been performed on a set of cytochrome P450 BM3 heme domains in which mutation of the highly conserved Phe393 induces significant variation in heme iron reduction potential. In previous work [Chen, Z., Ost, T.W.B., and Schelvis, J.P.M. (2004) Biochemistry 43, 1798-1808], a correlation between heme vinyl conformation and the heme iron reduction potential indicated a steric control by the protein over the distribution of electron density in the reduced heme cofactor. The current study aims to monitor changes in electron density on the ferrous heme cofactor following CO binding. In addition, ferric-NO complexes have been studied to investigate potential changes to the proximal Cys400 thiolate. We find that binding of CO to the ferrous heme domains results in a reorientation of the vinyl groups to a largely out-of-plane conformation, the extent of which correlates with the size of the residue at position 393. We conclude that FeII dπ back bonding to the CO ligand largely takes away the need for conjugation of the vinyl groups with the porphyrin ring to accommodate FeII dπ back bonding to the porphyrin ligand. The ferrous-CO and ferric-NO data are consistent with a small decrease in σ-electron donation from the proximal Cys400 thiolate in the F393A mutant and, to a lesser extent, the F393H mutant, potentially due to a small increase in hydrogen bonding to the proximal ligand. Phe393 seems strategically placed to preserve robust σ-electron donation to the heme iron and to fine-tune its electron density by limiting vinyl group rotation.

Metabolomics analysis of milk thistle lipids to identify drought-tolerant genes

Milk thistle is an oil and medicinal crop known as an alternative oil crop with a high level of unsaturated fatty acids, which makes it a favorable edible oil for use in food production. To evaluate the importance of Milk thistle lipids in drought tolerance, an experiment was performed in field conditions under three different water deficit levels (Field capacity (FC), 70% FC and 40% FC). After harvesting seeds of the plant, their oily and methanolic extracts were isolated, and subsequently, types and amounts of lipids were measured using GC-MS. Genes and enzymes engaged in biosynthesizing of these lipids were identified and their expression in Arabidopsis was investigated under similar conditions. The results showed that content of almost all measured lipids of milk thistle decreased under severe drought stress, but genes (belonged to Arabidopsis), which were involved in their biosynthetic pathway showed different expression patterns. Genes biosynthesizing lipids, which had significant amounts were selected and their gene and metabolic network were established. Two networks were correlated, and for each pathway, their lipids and respective biosynthesizing genes were grouped together. Four up-regulated genes including PXG3, LOX2, CYP710A1, PAL and 4 down-regulated genes including FATA2, CYP86A1, LACS3, PLA2-ALPHA were selected. The expression of these eight genes in milk thistle was similar to Arabidopsis under drought stress. Thus, PXG3, PAL, LOX2 and CYP86A1 genes that increased expression were selected for protein analysis. Due to the lack of protein structure of these genes in the milk thistle, modeling homology was performed for them. The results of molecular docking showed that the four proteins CYP86A1, LOX2, PAL and PXG3 bind to ligands HEM, 11O, ACT and LIG, respectively. HEM ligand was involved in production of secondary metabolites and dehydration tolerance, and HEM binding site remained conserved in various plants. CA ligands were involved in synthesis of cuticles and waxes. Overall, this study confirmed the importance of lipids in drought stress tolerance in milk thistle.

Cryo-EM reveals the architecture of the dimeric cytochrome P450 CYP102A1 enzyme and conformational changes required for redox partner recognition

Cytochrome P450 family 102 subfamily A member 1 (CYP102A1) is a self-sufficient flavohemeprotein and a highly active bacterial enzyme capable of fatty acid hydroxylation at a >3,000 min-1 turnover rate. The CYP102A1 architecture has been postulated to be responsible for its extraordinary catalytic prowess. However, the structure of a functional full-length CYP102A1 enzyme remains to be determined. Herein, we used a cryo-EM single-particle approach, revealing that full-length CYP102A1 forms a homodimer in which both the heme and FAD domains contact each other. The FMN domain of one monomer was located close to the heme domain of the other monomer, exhibiting a trans configuration. Moreover, full-length CYP102A1 is highly dynamic, existing in multiple conformational states, including open and closed states. In the closed state, the FMN domain closely contacts the FAD domain, whereas in the open state, one of the FMN domains rotates away from its FAD domain and traverses to the heme domain of the other monomer. This structural arrangement and conformational dynamics may facilitate rapid intraflavin and trans FMN-to-heme electron transfers (ETs). Results with a variant having a 12-amino-acid deletion in the CYP102A1 linker region, connecting the catalytic heme and the diflavin reductase domains, further highlighted the importance of conformational dynamics in the ET process. Cryo-EM revealed that the Δ12 variant homodimer is conformationally more stable and incapable of FMN-to-heme ET. We conclude that closed-to-open alternation is crucial for redox partner recognition and formation of an active ET complex for CYP102A1 catalysis.

Probing protein-protein and protein-substrate interactions in the dynamic membrane-associated ternary complex of cytochromes P450, b5, and reductase

Cytochrome P450 (cytP450) interacts with two redox partners, cytP450 reductase and cytochrome-b5, to metabolize substrates. Using NMR, we reveal changes in the dynamic interplay when all three proteins are incorporated into lipid nanodiscs in the absence and presence of substrates.

Glutamine-451 Confers Sensitivity to Oxidative Inhibition and Heme-Thiolate Sulfenylation of Cytochrome P450 4B1

Human cytochrome P450 (P450) family 4 enzymes are involved in the metabolism of fatty acids and the bioactivation of carcinogenic arylamines and toxic natural products, e.g., 4-ipomeanol. These and other drug-metabolizing P450s are redox sensitive, showing a loss of activity resulting from preincubation with H2O2 and recovery with mild reducing agents [Albertolle, M. W., et al. (2017) J. Biol. Chem. 292, 11230-11242]. The inhibition is due to sulfenylation of the heme-thiolate ligand, as determined by chemopreoteomics and spectroscopy. This phenomenon may have implications for chemical toxicity and observed disease-drug interactions, in which the decreased metabolism of P450 substrates occurs in patients with inflammatory diseases (e.g., influenza and autoimmunity). Human P450 1A2 was determined to be redox insensitive. To determine the mechanism underlying the differential redox sensitivity, molecular dynamics (MD) simulations were employed using the crystal structure of rabbit P450 4B1 (Protein Data Bank entry 5T6Q ). In simulating either the thiolate (Cys-S-) or the sulfenic acid (Cys-SOH) at the heme ligation site, MD revealed Gln-451 in either an "open" or "closed" conformation, respectively, between the cytosol and heme-thiolate cysteine. Mutation to either an isosteric leucine (Q451L) or glutamate (Q451E) abrogated the redox sensitivity, suggesting that this "open" conformation allows for reduction of the sulfenic acid and religation of the thiolate to the heme iron. In summary, MD simulations suggest that Gln-451 in P450 4B1 adopts conformations that may stabilize and protect the heme-thiolate sulfenic acid; mutating this residue destabilizes the interaction, producing a redox insensitive enzyme.

Substrate mediated redox partner selectivity of cytochrome P450

Investigating the interplay between cytochrome-P450 and its redox partners (CPR and cytochrome-b5) is vital for understanding the metabolism of most hydrophobic drugs. Dynamic structural interactions with the ternary complex, with and without substrates, captured by NMR reveal a gating mechanism for redox partners to promote P450 function.

Use of Phenoxyaniline Analogues To Generate Biochemical Insights into the Interactio n of Polybrominated Diphenyl Ether with CYP2B Enzymes

Human hepatic cytochromes P450 (CYP) are integral to xenobiotic metabolism. CYP2B6 is a major catalyst of biotransformation of environmental toxicants, including polybrominated diphenyl ethers (PBDEs). CYP2B substrates tend to contain halogen atoms, but the biochemical basis for this selectivity and for species specific determinants of metabolism has not been identified. Spectral binding titrations and inhibition studies were performed to investigate interactions of rat CYP2B1, rabbit CYP2B4, and CYP2B6 with a series of phenoxyaniline (POA) congeners that are analogues of PBDEs. For most congeners, there was a <3-fold difference between the spectral binding constants (K_S) and IC₅₀ values. In contrast, large discrepancies between these values were observed for POA and 3-chloro-4-phenoxyaniline. CYP2B1 was the enzyme most sensitive to POA congeners, so the Val-363 residue from that enzyme was introduced into CYP2B4 or CYP2B6. This substitution partially altered the protein-ligand interaction profiles to make them more similar to that of CYP2B1. Addition of cytochrome P450 oxidoreductase (POR) to titrations of CYP2B6 with POA or 2'4'5'TCPOA decreased the affinity of both ligands for the enzyme. Addition of cytochrome b₅ to a recombinant enzyme system containing POR and CYP2B6 increased the POA IC₅₀ value and decreased the 2'4'5'TCPOA IC₅₀ value. Overall, the inconsistency between K_S and IC₅₀ values for POA versus 2'4'5'TCPOA is largely due to the effects of redox partner binding. These results provide insight into the biochemical basis of binding of diphenyl ethers to human CYP2B6 and changes in CYP2B6-mediated metabolism that are dependent on POA congener and redox partner identity.

Role of cytochrome b5 in the modulation of the enzymatic activities of cytochrome P450 17α-hydroxylase/17,20-lyase (P450 17A1)

Cytochrome b5 (cyt b5) is a small hemoprotein that plays a significant role in the modulation of activities of an important steroidogenic enzyme, cytochrome P450 17α-hydroxylase/17,20-lyase (P450 17A1, CYP17A1). Located in the zona fasciculata and zona reticularis of the adrenal cortex and in the gonads, P450 17A1 catalyzes two different reactions in the steroidogenic pathway; the 17α-hydroxylation and 17,20-lyase, in the endoplasmic reticulum of these respective tissues. The activities of P450 17A1 are regulated by cyt b5 that enhances the 17,20-lyase reaction by promoting the coupling of P450 17A1 and cytochrome P450 reductase (CPR), allosterically. Cyt b5 can also act as an electron donor to enhance the 16-ene-synthase activity of human P450 17A1. In this review, we discuss the many roles of cyt b5 and focus on the modulation of CYP17A1 activities by cyt b5 and the mechanisms involved.

Role of the Proximal Cysteine Hydrogen Bonding Interaction in Cytochrome P450 2B4 Studied by Cryoreduction, Electron Paramagnetic Resonance, and Electron-Nuclear Double Resonance Spectroscopy

Crystallographic studies have shown that the F429H mutation of cytochrome P450 2B4 introduces an H-bond between His429 and the proximal thiolate ligand, Cys436, without altering the protein fold but sharply decreases the enzymatic activity and stabilizes the oxyferrous P450 2B4 complex. To characterize the influence of this hydrogen bond on the states of the catalytic cycle, we have used radiolytic cryoreduction combined with electron paramagnetic resonance (EPR) and (electron-nuclear double resonance (ENDOR) spectroscopy to study and compare their characteristics for wild-type (WT) P450 2B4 and the F429H mutant. (i) The addition of an H-bond to the axial Cys436 thiolate significantly changes the EPR signals of both low-spin and high-spin heme-iron(III) and the hyperfine couplings of the heme-pyrrole (14)N but has relatively little effect on the (1)H ENDOR spectra of the water ligand in the six-coordinate low-spin ferriheme state. These changes indicate that the H-bond introduced between His and the proximal cysteine decreases the extent of S → Fe electron donation and weakens the Fe(III)-S bond. (ii) The added H-bond changes the primary product of cryoreduction of the Fe(II) enzyme, which is trapped in the conformation of the parent Fe(II) state. In the wild-type enzyme, the added electron localizes on the porphyrin, generating an S = (3)/2 state with the anion radical exchange-coupled to the Fe(II). In the mutant, it localizes on the iron, generating an S = (1)/2 Fe(I) state. (iii) The additional H-bond has little effect on g values and (1)H-(14)N hyperfine couplings of the cryogenerated, ferric hydroperoxo intermediate but noticeably slows its decay during cryoannealing. (iv) In both the WT and the mutant enzyme, this decay shows a significant solvent kinetic isotope effect, indicating that the decay reflects a proton-assisted conversion to Compound I (Cpd I). (v) We confirm that Cpd I formed during the annealing of the cryogenerated hydroperoxy intermediate and that it is the active hydroxylating species in both WT P450 2B4 and the F429H mutant. (vi) Our data also indicate that the added H-bond of the mutation diminishes the reactivity of Cpd I.

F429 Regulation of Tunnels in Cytochrome P450 2B4: A Top Down Study of Multiple Molecular Dynamics Simulations

The root causes of the outcomes of the single-site mutation in enzymes remain by and large not well understood. This is the case of the F429H mutant of the cytochrome P450 (CYP) 2B4 enzyme where the substitution, on the proximal surface of the active site, of a conserved phenylalanine 429 residue with histidine seems to hamper the formation of the active species, Compound I (porphyrin cation radical-Fe^(IV) = O, Cpd I) from the ferric hydroperoxo (Fe^(III)OOH^-, Cpd 0) precursor. Here we report a study based on extensive molecular dynamic (MD) simulations of 4 CYP-2B4 point mutations compared to the WT enzyme, having the goal of better clarifying the importance of the proximal Phe429 residue on CYP 2B4 catalytic properties. To consolidate the huge amount of data coming from five simulations and extract the most distinct structural features of the five species studied we made an extensive use of cluster analysis. The results show that all studied single polymorphisms of F429, with different side chain properties: i) drastically alter the reservoir of conformations accessible by the protein, perturbing global dynamics ii) expose the thiolate group of residue Cys436 to the solvent, altering the electronic properties of Cpd0 and iii) affect the various ingress and egress channels connecting the distal sites with the bulk environment, altering the reversibility of these channels. In particular, it was observed that the wild type enzyme exhibits unique structural features as compared to all mutant species in terms of weak interactions (hydrogen bonds) that generate a completely different dynamical behavior of the complete system. Albeit not conclusive, the current computational investigation sheds some light on the subtle and critical effects that proximal single-site mutations can exert on the functional mechanisms of human microsomal CYPs which should go rather far beyond local structure characterization.