Crystal structure of a NADPH-cytochrome P450 oxidoreductase (CYPOR) and heme oxygenase 1 fusion protein implies a conformational change in CYPOR upon NADPH/NADP+ binding

A novel signature combing cuproptosis- and ferroptosis-related genes in sepsis-induced cardiomyopathy

Objective: Cardiac dysfunction caused by sepsis, usually termed sepsis-induced cardiomyopathy (SIC), is one of the most serious complications of sepsis, and ferroptosis can play a key role in this disease. In this study, we identified key cuproptosis- and ferroptosis-related genes involved in SIC and further explored drug candidates for the treatment of SIC. Methods: The GSE79962 gene expression profile of SIC patients was downloaded from the Gene Expression Omnibus database (GEO). The data was used to identify differentially expressed genes (DEGs) and to perform weighted correlation network analysis (WGCNA). Furthermore, Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses were conducted. Then, gene set enrichment analysis (GSEA) was applied to further analyze pathway regulation, with an adjusted p-value <0.05 and a false discovery rate (FDR) <0.25. Ferroptosis-related genes were obtained from the FerrDb V2 database, and cuproptosis-related genes were obtained from the literature. We constructed a novel signature (CRF) by combing cuproptosis-related genes with ferroptosis-related genes using the STRING website. The SIC hub genes were obtained by overlapping DEGs, WGCNA-based hub genes and CRF genes, and receiver operating characteristic (ROC) curve analysis was used to determine the diagnostic value of hub genes. A transcription factor-microRNA-hub gene network was also constructed based on the miRnet database. Finally, potential therapeutic compounds for SIC were predicted based on the Drug Gene Interaction Database. Results: We identified 173 DEGs in SIC patients. Four hub modules and 411 hub genes were identified by WGCNA. A total of 144 genes were found in the CRF. Then, POR, SLC7A5 and STAT3 were identified as intersecting hub genes and their diagnostic values were confirmed with ROC curves. Drug screening identified 15 candidates for SIC treatment. Conclusion: We revealed that the cuproptosis- and ferroptosis-related genes, POR, SLC7A5 and STAT3, were significantly correlated with SIC and we also predicted therapeutic drugs for these targets. The findings from this study will make contributions to the development of treatments for SIC.

Enzymological and structural characterization of Arabidopsis thaliana heme oxygenase-1

Arabidopsis thaliana heme oxygenase‐1 (AtHO‐1), a metabolic enzyme in the heme degradation pathway, serves as a prototype for study of the bilin‐related functions in plants. Past biological analyses revealed that AtHO‐1 requires ferredoxin‐NADP⁺ reductase (FNR) and ferredoxin for its enzymatic activity. Here, we characterized the binding and degradation of heme by AtHO‐1, and found that ferredoxin is a dispensable component of the reducing system that provides electrons for heme oxidation. Furthermore, we reported the crystal structure of heme‐bound AtHO‐1, which demonstrates both conserved and previously undescribed features of plant heme oxygenases. Finally, the electron transfer pathway from FNR to AtHO‐1 is suggested based on the known structural information.

Conformational Rearrangements in the Redox Cycling of NADPH-Cytochrome P450 Reductase from Sorghum bicolor Explored with FRET and Pressure-Perturbation Spectroscopy

Simple Summary

NADPH-cytochrome P450 reductase (CPR) enzymes are known to undergo an ample conformational transition between the closed and open states in the process of their redox cycling. To explore the conformational landscape of CPR from the potential biofuel crop Sorghum bicolor (SbCPR), we incorporated a FRET donor/acceptor pair into the enzyme and employed rapid scanning stop-flow and pressure perturbation spectroscopy to characterize the equilibrium between its open and closed states at different stages of the redox cycle. Our results suggest the presence of several open conformational sub-states differing in the system volume change associated with the opening transition (ΔV⁰). Although the closed conformation always predominates in the conformational landscape, the population of the open conformations increases by order of magnitude upon the two-electron reduction and the formation of the disemiquinone state of the enzyme. In addition to elucidating the functional choreography of plant CPRs, our study demonstrates the high exploratory potential of a combination of the pressure-perturbation approach with the FRET-based monitoring of protein conformational transitions.

Abstract

NADPH-cytochrome P450 reductase (CPR) from Sorghum bicolor (SbCPR) serves as an electron donor for cytochrome P450 essential for monolignol and lignin production in this biofuel crop. The CPR enzymes undergo an ample conformational transition between the closed and open states in their functioning. This transition is triggered by electron transfer between the FAD and FMN and provides access of the partner protein to the electron-donating FMN domain. To characterize the electron transfer mechanisms in the monolignol biosynthetic pathway better, we explore the conformational transitions in SbCPR with rapid scanning stop-flow and pressure-perturbation spectroscopy. We used FRET between a pair of donor and acceptor probes incorporated into the FAD and FMN domains of SbCPR, respectively, to characterize the equilibrium between the open and closed states and explore its modulation in connection with the redox state of the enzyme. We demonstrate that, although the closed conformation always predominates in the conformational landscape, the population of open state increases by order of magnitude upon the formation of the disemiquinone state. Our results are consistent with several open conformation sub-states differing in the volume change (ΔV⁰) of the opening transition. While the ΔV⁰ characteristic of the oxidized enzyme is as large as −88 mL/mol, the interaction of the enzyme with the nucleotide cofactor and the formation of the double-semiquinone state of CPR decrease this value to −34 and −18 mL/mol, respectively. This observation suggests that the interdomain electron transfer in CPR increases protein hydration, while promoting more open conformation. In addition to elucidating the functional choreography of plant CPRs, our study demonstrates the high exploratory potential of a combination of the pressure-perturbation approach with the FRET-based monitoring of protein conformational transitions.

Identification of Crucial Genes and Infiltrating Immune Cells Underlying Sepsis-Induced Cardiomyopathy via Weighted Gene Co-Expression Network Analysis

Sepsis-induced cardiomyopathy (SIC), with a possibly reversible cardiac dysfunction, is a potential complication of septic shock. Despite quite a few mechanisms including the inflammatory mediator, exosomes, and mitochondrial dysfunction, having been confirmed in the existing research studies we still find it obscure about the overall situation of gene co-expression that how they can affect the pathological process of SIC. Thus, we intended to find out the crucial hub genes, biological signaling pathways, and infiltration of immunocytes underlying SIC. It was weighted gene co-expression network analysis that worked as our major method on the ground of the gene expression profiles: hearts of those who died from sepsis were compared to hearts donated by non-failing humans which could not be transplanted for technical reasons (GSE79962). The top 25 percent of variant genes were abstracted to identify 10 co-expression modules. In these modules, brown and green modules showed the strongest negative and positive correlation with SIC, which were primarily enriched in the bioenergy metabolism, immunoreaction, and cell death. Next, nine genes (LRRC39, COQ10A, FSD2, PPP1R3A, TNFRSF11B, IL1RAP, DGKD, POR, and THBS1) including two downregulated and seven upregulated genes which were chosen as hub genes that meant the expressive level of which was higher than the counterparts in control groups. Then, the gene set enrichment analysis (GSEA) demonstrated a close relationship of hub genes to the cardiac metabolism and the necroptosis and apoptosis of cells in SIC. Concerning immune cells infiltration, a higher level of neutrophils and B cells native and a lower level of mast cells resting and plasma cells had been observed in patients with SIC. In general, nine candidate biomarkers were authenticated as a reliable signature for deeper exploration of basic and clinical research studies on SIC.

Crystal structure of higher plant heme oxygenase-1 and its mechanism of interaction with ferredoxin

Heme oxygenase (HO) converts heme to carbon monoxide, biliverdin, and free iron, products that are essential in cellular redox signaling and iron recycling. In higher plants, HO is also involved in the biosynthesis of photoreceptor pigment precursors. Despite many common enzymatic reactions, the amino acid sequence identity between plant-type and other HOs is exceptionally low (∼19.5%), and amino acids that are catalytically important in mammalian HO are not conserved in plant-type HOs. Structural characterization of plant-type HO is limited to spectroscopic characterization by electron spin resonance, and it remains unclear how the structure of plant-type HO differs from that of other HOs. Here, we have solved the crystal structure of Glycine max (soybean) HO-1 (GmHO-1) at a resolution of 1.06 Å and carried out the isothermal titration calorimetry measurements and NMR spectroscopic studies of its interaction with ferredoxin, the plant-specific electron donor. The high-resolution X-ray structure of GmHO-1 reveals several novel structural components: an additional irregularly structured region, a new water tunnel from the active site to the surface, and a hydrogen-bonding network unique to plant-type HOs. Structurally important features in other HOs, such as His ligation to the bound heme, are conserved in GmHO-1. Based on combined data from X-ray crystallography, isothermal titration calorimetry, and NMR measurements, we propose the evolutionary fine-tuning of plant-type HOs for ferredoxin dependency in order to allow adaptation to dynamic pH changes on the stroma side of the thylakoid membrane in chloroplast without losing enzymatic activity under conditions of fluctuating light.

Conformational Equilibrium of NADPH-Cytochrome P450 Oxidoreductase Is Essential for Heme Oxygenase Reaction

Heme oxygenase (HO) catalyzes heme degradation using electrons supplied by NADPH-cytochrome P450 oxidoreductase (CPR). Electrons from NADPH flow first to FAD, then to FMN, and finally to the heme in the redox partner. Previous biophysical analyses suggest the presence of a dynamic equilibrium between the open and the closed forms of CPR. We previously demonstrated that the open-form stabilized CPR (ΔTGEE) is tightly bound to heme-HO-1, whereas the reduction in heme-HO-1 coupled with ΔTGEE is considerably slow because the distance between FAD and FMN in ΔTGEE is inappropriate for electron transfer from FAD to FMN. Here, we characterized the enzymatic activity and the reduction kinetics of HO-1 using the closed-form stabilized CPR (147CC514). Additionally, we analyzed the interaction between 147CC514 and heme-HO-1 by analytical ultracentrifugation. The results indicate that the interaction between 147CC514 and heme-HO-1 is considerably weak, and the enzymatic activity of 147CC514 is markedly weaker than that of CPR. Further, using cryo-electron microscopy, we confirmed that the crystal structure of ΔTGEE in complex with heme-HO-1 is similar to the relatively low-resolution structure of CPR complexed with heme-HO-1 in solution. We conclude that the "open-close" transition of CPR is indispensable for electron transfer from CPR to heme-HO-1.

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

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

Coupling of Redox and Structural States in Cytochrome P450 Reductase Studied by Molecular Dynamics Simulation

Cytochrome P450 reductase (CPR) is the key protein that regulates the electron transfer from NADPH to various heme-containing monooxygenases. CPR has two flavin-containing domains: one with flavin adenine dinucleotide (FAD), called FAD domain, and the other with flavin mononucleotide (FMN), called FMN domain. It is considered that the electron transfer occurs via FAD and FMN (NADPH → FAD → FMN → monooxygenase) and is regulated by an interdomain open-close motion. It is generally thought that the structural state is coupled with the redox state, which, however, has not yet been firmly established. In this report, we studied the coupling of the redox and the structural states by full-scale molecular dynamics (MD) simulation of CPR (total 86.4 μs). Our MD result showed that while CPR predominantly adopts the closed state both in the oxidized and reduced states, it exhibits a tendency to open in the reduced state. We also found a correlation between the FAD-FMN distance and the predicted FMN-monooxygenase distance, which is embedded in the equilibrium thermal fluctuation of CPR. Based on these results, a physical mechanism for the electron transfer by CPR is discussed.