Proteomics, modeling, and fluorescence assays delineate cytochrome b5 residues involved in binding and stimulation of cytochrome P450 17A1 17,20-lyase.

Cytochrome b5 (b5) is known to stimulate some catalytic activities of cytochrome P450 (P450, CYP) enzymes, although mechanisms still need to be defined. The reactions most strongly enhanced by b5 are the 17,20-lyase reactions of P450 17A1 involved in steroid biosynthesis. We had previously used a fluorescently labeled human b5 variant (Alexa 488-T70C-b5) to characterize human P450 17A1-b5 interactions, but subsequent proteomic analyses indicated that lysines in b5 were also modified with Alexa 488 maleimide in addition to Cys-70, due to disulfide dimerization of the T70C mutant. A series of b5 variants were constructed with Cys replacements for the identified lysine residues and labeled with the dye. Fluorescence attenuation and the function of b5 in the steroid lyase reaction depended on the modified position. Apo-b5 (devoid of heme group) studies revealed the lack of involvement of the b5 heme in the fluorescence attenuation. A structural model of b5 with P450 17A1 was predicted using AlphaFold-Multimer algorithms/Rosetta docking, based upon the individual structures, which predicted several new contacts not previously reported, that is, interactions of b5 Glu-48:17A1 Arg-347, b5 Glu-49:17A1 Arg-449, b5 Asp-65:17A1 Arg-126, b5 Asp-65:17A1 Arg-125, and b5 Glu-61:17A1 Lys-91. Fluorescence polarization assays with two modified b5 variants yielded Kd values (for b5-P450 17A1) of 120 to 380 nM, the best estimate of binding affinity. We conclude that both monomeric and dimeric b5 can bind to P450 17A1 and stimulate activity. Results with the mutants indicate that several Lys residues in b5 are sensitive to the interaction with P450 17A1, including Lys-88 and Lys-91.

Function and regulation of a steroidogenic CYP450 enzyme in the mitochondrion of Toxoplasma gondii.

As an obligate intracellular parasite, Toxoplasma gondii must import essential nutrients from the host cell into the parasitophorous vacuole. We previously reported that the parasite scavenges cholesterol from host endocytic organelles for incorporation into membranes and storage as cholesteryl esters in lipid droplets. In this study, we have investigated whether Toxoplasma utilizes cholesterol as a precursor for the synthesis of metabolites, such as steroids. In mammalian cells, steroidogenesis occurs in mitochondria and involves membrane-bound type I cytochrome P450 oxidases that are activated through interaction with heme-binding proteins containing a cytochrome b5 domain, such as members of the membrane-associated progesterone receptor (MAPR) family. Our LC-MS targeted lipidomics detect selective classes of hormone steroids in Toxoplasma, with a predominance for anti-inflammatory hydroxypregnenolone species, deoxycorticosterone and dehydroepiandrosterone. The genome of Toxoplasma contains homologs encoding a single type I CYP450 enzyme (we named TgCYP450mt) and a single MAPR (we named TgMAPR). We showed that TgMAPR is a hemoprotein with conserved residues in a heme-binding cytochrome b5 domain. Both TgCYP450 and TgMAPR localize to the mitochondrion and show interactions in in situ proximity ligation assays. Genetic ablation of cyp450mt is not tolerated by Toxoplasma; we therefore engineered a conditional knockout strain and showed that iΔTgCYP450mt parasites exhibit growth impairment in cultured cells. Parasite strains deficient for mapr could be generated; however, ΔTgMAPR parasites suffer from poor global fitness, loss of plasma membrane integrity, aberrant mitochondrial cristae, and an abnormally long S-phase in their cell cycle. Compared to wild-type parasites, iΔTgCYP450mt and ΔTgMAPR lost virulence in mice and metabolomics studies reveal that both mutants have reduced levels of steroids. These observations point to a steroidogenic pathway operational in the mitochondrion of a protozoan that involves an evolutionary conserved TgCYP450mt enzyme and its binding partner TgMAPR.

The association between CYB5A gene rs1790834 polymorphism and clinical improvement after 6 months of leflunomide treatment in women with rheumatoid arthritis.

INTRODUCTION/OBJECTIVES: Rheumatoid arthritis (RA) is a disease affecting around 1% of the population in developed countries and can be treated with leflunomide. The higher prevalence of RA among women and numerous previous studies suggested the crucial role of sex hormones. Cytochrome CYB5A regulates the synthesis of androgens. Therefore, the aim of this study was to determine the association between common CYB5A gene polymorphism and the response to leflunomide in women with RA. METHODS: This study included 111 patients. All of them received oral leflunomide monotherapy at a dose of 20 mg daily. Women were genotyped for the presence of CYB5A rs1790834 polymorphism and evaluated monthly for 6 months following the initiation of treatment. RESULTS: After 6 months of therapy, patients with the GG genotype had higher DAS28 values and less improvement in DAS28 compared to patients with the GA and AA genotypes (p = 0.04). No statistically significant differences were found in relation to other disease activity parameters. CONCLUSIONS: The results of the current study suggest a possible association of the CYB5A rs1790834 polymorphism with some disease activity parameters in RA patients treated with leflunomide during the initial therapy period. However, confirmation of the effect of this polymorphism on the efficacy of leflunomide treatment requires further studies. Key Points • Leflunomide is the synthetic disease-modifying anti-rheumatic drug used in the therapy of rheumatoid arthritis. • CYB5A rs1790834 gene polymorphism may influence the clinical improvement after 6 months of leflunomide treatment in women with rheumatoid arthritis.

Inflammatory phenotype, clinicopathologic variables, and effects of long-term methylene blue in dogs with hereditary methemoglobinemia caused by cytochrome b5 reductase deficiency.

OBJECTIVE: To determine whether dogs with cytochrome b5 reductase (CYB5R) deficiency have a constitutive proinflammatory phenotype, characterize hematologic and serum chemistry results, and describe changes in methemoglobin (MetHb) levels and serum C-reactive protein (CRP) concentrations after long-term per os (PO) methylene blue (MB) therapy. ANIMALS: 21 client-owned dogs (CYB5R deficient, n = 10; healthy controls, 11). PROCEDURES: In this prospective, case-control study, methemoglobin levels were measured using a blood gas analyzer with co-oximetry. Plasma tumor necrosis factor-α (TNF-α), interleukin-6 (IL-6), and interleukin-10 (IL-10) concentrations were measured using a canine-specific multiplex bead-based assay. Serum CRP concentrations were measured with a canine-specific commercial ELISA kit. Serum CRP concentration and MetHb levels were measured in 6 dogs with CYB5R deficiency after ≥ 60 days of PO MB therapy. RESULTS: As expected, MetHb levels were higher in dogs with CYB5R deficiency compared to controls (P < .001). Plasma TNF-α, IL-6, IL-10, and serum CRP concentrations were no different between CYB5R-deficient and control dogs. Dogs with CYB5R deficiency had lower absolute lymphocyte (P = .005) and eosinophil counts (P = .04) and higher alanine transaminase (P = .04) and alkaline phosphatase activity (P = .02) than controls, but these changes were not clinically relevant. Methemoglobin levels decreased after PO MB therapy (P = .03). CLINICAL RELEVANCE: These results suggest that otherwise healthy dogs with CYB5R deficiency do not have a constitutive proinflammatory phenotype and clinically relevant abnormalities in hematologic and serum chemistry panels are not expected. Dogs with decreased quality of life attributed to methemoglobinemia from CYB5R deficiency might benefit from PO MB therapy.

Structural characterisation of a MAPR-related archaeal cytochrome b5M protein.

We recently reported that the membrane-associated progesterone receptor (MAPR) protein family (mammalian members: PGRMC1, PGRMC2, NEUFC and NENF) originated from a new class of prokaryotic cytochrome b5 (cytb5 ) domain proteins, called cytb5M (MAPR-like). Relative to classical cytb5 proteins, MAPR and ctyb5M proteins shared unique sequence elements and a distinct heme-binding orientation at an approximately 90° rotation relative to classical cytb5 , as demonstrated in the archetypal crystal structure of a cytb5M protein (PDB accession number 6NZX). Here, we present the crystal structure of an archaeal cytb5M domain (Methanococcoides burtonii WP_011499504.1, PDB:6VZ6). It exhibits similar heme binding to the 6NZX cytb5M , supporting the deduction that MAPR-like heme orientation was inherited from the prokaryotic ancestor of the original eukaryotic MAPR gene.

Regulating Effect of Cytochrome b5 Overexpression on Human Breast Cancer Cells.

Imbalance in the cellular redox system is thought to be associated with the induction and progression of breast cancers, and heme proteins may regulate the redox balance. Cytochrome b5 (Cyt b5) is a small mitochondrial heme protein. Its function and regulating mechanism in breast cancer remain unknown. In this study, we elucidated the level of endogenous oxidative stress in breast cancer cells, MCF-7 cells (hormone receptor-positive cells) and MDA-MB-231 cells (triple-negative cells), and investigated the difference in Cyt b5 content. Based on the low content of Cyt b5 in MDA-MB-231 cells, the overexpression of Cyt b5 was found to regulate the oxidative stress and apoptosis cascades, including ERK1/2 and Akt signaling pathways. The overexpressed Cyt b5 MDA-MB-231 cells were shown to exhibit decreased oxidative stress, less phosphorylation of ERK1/2 and Akt, and less cleavage of caspases 3 and 9 upon treatment with H2O2, as compared to those of normal MDA-MB-231 cells. Moreover, the overexpressed Cyt b5 most likely functioned by interacting with its protein partner, Cyt c, as suggested by co-immunoprecipitation studies. These results indicated that Cyt b5 has different effects on breast cancer cells of different phenotypes, which provides useful information for understanding the multiple roles of Cyt b5 and provides clues for clinical treatment.

Significant genes in response to low temperature in Fenneropenaeus chinensis screened through multiple transcriptome group comparisons.

Fenneropenaeus chinensis is a migratory marine species with a suitable growth at 18-30°C. To prolong breeding season and reduce mortality in winter, breeding new shrimp varieties with cold tolerance is essential. Genes upregulated and highly expressed at low temperature are reasonable candidate genetic markers for the breeding of cold tolerant strain variants. This study screened genes with these features by comparing multiple low-vs. normal-temperature transcriptome groups. The results showed that nine genes were upregulated and highly expressed at low temperature in more than seven of the nine comparison groups. Six of them were identified as genes encoding transcription factor ATF2, RNA recognition motif domain-containing protein, cytochrome b5-like protein, troponin C, tubulin alpha-1, and 18S/5.8S/28S rRNA, respectively. Cold-inducible upregulations of ATF2, cytochrome b5, and rRNAs were novel findings in this study. The other three novel genes were predicted to encode a membrane-bound extracellular protein and two lncRNAs. Four of the screened genes were verified by real time RT-PCR, and their expression levels were consistent with the sequencing results, demonstrating the accuracy of the transcriptome sequencing data. Function analysis showed that ATF2 might be the master transcription factor regulating the expressions of proteins involved in cellular responses to cold. The other genes played a role in events such as enhancing translation, increasing energy, inhibiting apoptosis, and preserving cell integrity. The expression features of these nine genes suggested that they were of great significance to the cold tolerance of shrimp.

Biochemical Investigation of Membrane-Bound Cytochrome b5 and the Catalytic Domain of Cytochrome b5 Reductase from Arabidopsis thaliana.

The endoplasmic reticulum (ER) membrane of plant cells contains several enzymes responsible for the biosynthesis of a diverse range of molecules essential for plant growth and holds potential for industrial applications. Many of these enzymes are dependent on electron transfer proteins to sustain their catalytic cycles. In plants, two crucial ER-bound electron transfer proteins are cytochrome b5 and cytochrome b5 reductase, which catalyze the stepwise transfer of electrons from NADH to redox enzymes such as fatty acid desaturases, cytochrome P450s, and plant aldehyde decarbonylase. Despite the high significance of plant cytochrome b5 and cytochrome b5 reductase, they have eluded detailed characterization to date. Here, we overexpressed the full-length membrane-bound cytochrome b5 isoform B from the model plant Arabidopsis thaliana in Escherichia coli, purified the protein employing detergents as well as styrene-maleic acid (SMA) copolymers, and biochemically characterized the protein. The SMA-encapsulated cytochrome b5 exhibits a discoidal shape and the characteristic features of the active heme-bound state. We also overexpressed and purified the soluble domain of cytochrome b5 reductase from A. thaliana, establishing its activity, stability, and kinetic parameters. Further, we demonstrated that the plant cytochrome b5, purified in detergents and styrene maleic acid lipid particles (SMALPs), readily accepts electrons from the cognate plant cytochrome b5 reductase and distant electron mediators such as plant NADPH-cytochrome P450 oxidoreductase and cyanobacterial NADPH-ferredoxin reductase. We also measured the kinetic parameters of cytochrome b5 reductase for cytochrome b5. Our studies are the first to report the purification and detailed biochemical characterization of the plant cytochrome b5 and cytochrome b5 reductase from the bacterial overexpression system.

An alternative CYB5A transcript is expressed in aneuploid ALL and enriched in relapse.

BACKGROUND: B-cell precursor acute lymphoblastic leukemia (BCP-ALL) is a genetically heterogenous malignancy with poor prognosis in relapsed adult patients. The genetic basis for relapse in aneuploid subtypes such as near haploid (NH) and high hyperdiploid (HeH) BCP-ALL is only poorly understood. Pathogenic genetic alterations remain to be identified. To this end, we investigated the dynamics of genetic alterations in a matched initial diagnosis-relapse (ID-REL) BCP-ALL cohort. Here, we firstly report the identification of the novel genetic alteration CYB5Aalt, an alternative transcript of CYB5A, in two independent cohorts. METHODS: We identified CYB5alt in the RNAseq-analysis of a matched ID-REL BCP-ALL cohort with 50 patients and quantified its expression in various molecular BCP-ALL subtypes. Findings were validated in an independent cohort of 140 first diagnosis samples from adult BCP-ALL patients. Derived from patient material, the alternative open reading frame of CYB5Aalt was cloned (pCYB5Aalt) and pCYB5Aalt or the empty vector were stably overexpressed in NALM-6 cells. RNA sequencing was performed of pCYB5Aalt clones and empty vector controls followed by differential expression analysis, gene set enrichment analysis and complementing cell death and viability assays to determine functional implications of CYB5Aalt. RESULTS: RNAseq data analysis revealed non-canonical exon usage of CYB5Aalt starting from a previously undescribed transcription start site. CYB5Aalt expression was increased in relapsed BCP-ALL and its occurrence was specific towards the shared gene expression cluster of NH and HeH BCP-ALL in independent cohorts. Overexpression of pCYB5Aalt in NALM-6 cells induced a distinct transcriptional program compared to empty vector controls with downregulation of pathways related to reported functions of CYB5A wildtype. Interestingly, CYB5A wildtype expression was decreased in CYB5Aalt samples in silico and in vitro. Additionally, pCYB5Aalt NALM-6 elicited a more resistant drug response. CONCLUSIONS: Across all age groups, CYB5Aalt was the most frequent secondary genetic event in relapsed NH and HeH BCP-ALL. In addition to its high subgroup specificity, CYB5Aalt is a novel candidate to be potentially implicated in therapy resistance in NH and HeH BCP-ALL. This is underlined by overexpressing CYB5Aalt providing first evidence for a functional role in BCL2-mediated apoptosis.

Functional variants in cytochrome b5 type A (CYB5A) are enriched in Southwest American Indian individuals and associate with obesity.

OBJECTIVE: This study aimed to identify genetic variants enriched in Southwest American Indian (SWAI) individuals that associate with BMI. METHODS: Whole genome sequencing data (n = 296) were used to identify potentially functional variants that are common in SWAI individuals (minor allele frequency ≥10%) but rare in other ethnic groups (minor allele frequency < 0.1%). Enriched variants were tested for association with BMI in 5,870 SWAI individuals. One variant was studied using a luciferase reporter, and haplotypes that included this variant were analyzed for association with various measures of obesity (n = 917-5,870), 24-hour energy expenditure (24-h EE; n = 419), and skeletal muscle biopsy expression data (n = 207). RESULTS: A 5' untranslated region variant in cytochrome b5 type A (CYB5A), rs548402150, met the enrichment criteria and associated with increased BMI (ß = 2%, p = 0.004). Functionally, rs548402150 decreased luciferase expression by 30% (p = 0.003) and correlated with decreased skeletal muscle CYB5A expression (ß = -0.5 SD, p = 0.0008). Combining rs548402150 with two splicing quantitative trait loci in CYB5A identified a haplotype carried almost exclusively in SWAI individuals that associated with increased BMI (ß = 3%, p = 0.0003) and decreased CYB5A expression, whereas the most common haplotype in all ethnic groups associated with lower BMI and percentage of body fatness, increased 24-h EE, and increased CYB5A expression. CONCLUSIONS: Further studies on the effects of CYB5A on 24-h EE and BMI may provide insights into obesity-related physiology.

Identification and functional validation of novel pharmacogenomic variants using a next-generation sequencing-based approach for clinical pharmacogenomics.

Inter-individual variability in pharmacokinetics and drug response is heavily influenced by single-nucleotide variants (SNVs) and copy-number variations (CNVs) in genes with importance for drug disposition. Nowadays, a plethora of studies implement next generation sequencing to capture rare and novel pharmacogenomic (PGx) variants that influence drug response. To address these issues, we present a comprehensive end-to-end analysis workflow, beginning from targeted PGx panel re-sequencing to in silico analysis pipelines and in vitro validation assays. Specifically, we show that novel pharmacogenetic missense variants that are predicted or putatively predicted to be functionally deleterious, significantly alter protein activity levels of CYP2D6 and CYP2C19 proteins. We further demonstrate that variant priorization pipelines tailored with functional in vitro validation assays provide supporting evidence for the deleterious effect of novel PGx variants. The proposed workflow could provide the basis for integrating next-generation sequencing for PGx testing into routine clinical practice.

Mechanism of electron transfers mediated by cytochromes c and b5 in mitochondria and endoplasmic reticulum: classical and murburn perspectives.

We explore the mechanism of electron transfers mediated by cytochrome c, a soluble protein involved in mitochondrial oxidative phosphorylation and cytochrome b5, a microsomal membrane protein acting as a redox aide in xenobiotic metabolism. We found minimal conservation in the sequence and surface amino acid residues of cytochrome c/b5 proteins among divergent species. Therefore, we question the evolutionary logic for electron transfer (ET) occurring through affinity binding via recognition of specific surface residues/topography. Also, analysis of putative protein-protein interactions in the crystal structures of these proteins and their redox partners did not point to any specific interaction logic. A comparison of the kinetic and thermodynamic constants of wildtype vs. mutants did not provide strong evidence to support the binding-based ET paradigm, but indicated support for diffusible reactive species (DRS)-mediated process. Topographically divergent cytochromes from one species have been substituted for reaction with proteins from other species, implying the involvement of non-specific interactions. We provide a viable alternative (murburn concept) to classical protein-protein binding-based long range ET mechanism. To account for the promiscuity of interactions and solvent-accessible hemes, we propose that the two proteins act as non- specific redox capacitors, mediating one-electron redox equilibriums involving DRS and unbound ions.Communicated by Ramaswamy H. Sarma.

Morphologic and Molecular Characterization of Adrenals and Adrenal Rest Affected by Congenital Adrenal Hyperplasia.

Introduction: Adrenocortical hyperplasia and adrenal rest tumor (ART) formation are common in congenital adrenal hyperplasia (CAH). Although driven by excessive corticotropin, much is unknown regarding the morphology and transformation of these tissues. Our study objective was to characterize CAH-affected adrenals and ART and compare with control adrenal and gonadal tissues. Patients/Methods: CAH adrenals, ART and control tissues were analyzed by histology, immunohistochemistry, and transcriptome sequencing. We investigated protein expression of the ACTH receptor (MC2R), steroidogenic (CYP11B2, CYP11B1, CYB5A) and immune (CD20, CD3, CD68) biomarkers, and delta-like 1 homolog (DLK1), a membrane bound protein broadly expressed in fetal and many endocrine cells. RNA was isolated and gene expression was analyzed by RNA sequencing (RNA-seq) followed by principle component, and unsupervised clustering analyses. Results: Based on immunohistochemistry, CAH adrenals and ART demonstrated increased zona reticularis (ZR)-like CYB5A expression, compared to CYP11B1, and CYP11B2, markers of zona fasciculata and zona glomerulosa respectively. CYP11B2 was mostly absent in CAH adrenals and absent in ART. DLK1 was present in CAH adrenal, ART, and also control adrenal and testis, but was absent in control ovary. Increased expression of adrenocortical marker MC2R, was observed in CAH adrenals compared to control adrenal. Unlike control tissues, significant nodular lymphocytic infiltration was observed in CAH adrenals and ART, with CD20 (B-cell), CD3 (T-cell) and CD68 (macrophage/monocyte) markers of inflammation. RNA-seq data revealed co-expression of adrenal MC2R, and testis-specific INSL3, HSD17B3 in testicular ART indicating the presence of both gonadal and adrenal features, and high expression of DLK1 in ART, CAH adrenals and control adrenal. Principal component analysis indicated that the ART transcriptome was more similar to CAH adrenals and least similar to control testis tissue. Conclusions: CAH-affected adrenal glands and ART have similar expression profiles and morphology, demonstrating increased CYB5A with ZR characteristics and lymphocytic infiltration, suggesting a common origin that is similarly affected by the abnormal hormonal milieu. Immune system modulators may play a role in tumor formation of CAH.

A genome-wide association study reveals that the cytochrome b5 involved in seed reserve mobilization during seed germination in rice.

KEY MESSAGE: A candidate gene cytochrome b5 for the major QTL qSRMP9 for rice seed reserve mobilization was validated during seed germination using a genome-wide association study approach. Seed reserve mobilization plays important roles in the early seedling growth in rice. However, the genetic basis underlying this process is poorly understood. In this study, the genetic architecture of variation in seed reserve mobilization during seed germination was studied using a genome-wide association study approach in rice. Three quantitative trait loci (QTL) including qSRMP6, qSRMP9, and qSRMP12 for seed reserve mobilization percentage were identified. In which, the candidate gene cytochrome b5 (OsCyb5) for the major QTL qSRMP9 was validated. Disruption of this gene in Oscyb5 mutants reduced the seed reserve mobilization and seedling growth compared with wild-type (WT) in rice. There were no significant differences of grain size, starch, protein and total soluble sugar content in the mature grains between Oscyb5 mutants and WT. However, the α-amylase activity in the germinating seeds of Oscyb5 mutants was significantly decreased compared to that of WT, and then, the starch and sugar mobilization and the glucose accumulation during seed germination were significantly decreased in Oscyb5 mutants. Two elite haplotypes of OsCyb5 associated with the higher seed reserve mobilization percentage and its elite single nucleotide polymorphism variations were mainly existed in the INDICA and AUS accessions. The natural variation of OsCyb5 contributing to seed reserve mobilization might be useful for the future rice breeding.

Lack of association between CYB5A gene rs1790834 polymorphism and the response to leflunomide in women with rheumatoid arthritis.

AIM: Leflunomide is a disease-modifying antirheumatic drug used in therapy for rheumatoid arthritis (RA). Previous studies indicated that oestrogens and androgens may affect the response to leflunomide in RA patients. The synthesis of androgens is regulated by cytochrome CYB5A. The aim of this study was to examine the association between the CYB5A gene rs1790834 polymorphism and the response to leflunomide in women with RA. METHODS: The study included 111 women diagnosed with RA. Leflunomide was administered in monotherapy at a dose of 20 mg/day. All patients underwent a monthly evaluation for 12 months after the initiation of treatment with leflunomide. RESULTS: After 12 months of therapy, the changes in individual disease activity parameters, such as: DAS28, ESR, CRP and VAS, were not statistically significantly different between rs1790834 genotypes in the Kruskal-Wallis test. CONCLUSIONS: The results of our study suggest lack of statistically significant association between the CYB5A gene rs1790834 polymorphism and the response to leflunomide in women with RA.

Genome Mining and Evolutionary Analysis Reveal Diverse Type III Polyketide Synthase Pathways in Cyanobacteria.

Cyanobacteria are prolific producers of natural products, including polyketides and hybrid compounds thereof. Type III polyketide synthases (PKSs) are of particular interest, due to their wide substrate specificity and simple reaction mechanism, compared with both type I and type II PKSs. Surprisingly, only two type III PKS products, hierridins, and (7.7)paracyclophanes, have been isolated from cyanobacteria. Here, we report the mining of 517 cyanobacterial genomes for type III PKS biosynthesis gene clusters. Approximately 17% of the genomes analyzed encoded one or more type III PKSs. Together with already characterized type III PKSs, the phylogeny of this group of enzymes was investigated. Our analysis showed that type III PKSs in cyanobacteria evolved into three major lineages, including enzymes associated with 1) (7.7)paracyclophane-like biosynthesis gene clusters, 2) hierridin-like biosynthesis gene clusters, and 3) cytochrome b5 genes. The evolutionary history of these enzymes is complex, with some sequences partitioning primarily according to speciation and others putatively according to their reaction type. Protein modeling showed that cyanobacterial type III PKSs generally have a smaller active site cavity (mean = 109.035 Å3) compared with enzymes from other organisms. The size of the active site did not correlate well with substrate size, however, the "Gatekeeper" amino acid residues within the active site were strongly correlated to enzyme phylogeny. Our study provides unprecedented insight into the distribution, diversity, and molecular evolution of cyanobacterial type III PKSs, which could facilitate the discovery, characterization, and exploitation of novel enzymes, biochemical pathways, and specialized metabolites from this biosynthetically talented clade of microorganisms.

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

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

Exome sequencing of a Pakistani family with spastic paraplegia identified an 18 bp deletion in the cytochrome B5 domain of FA2H.

Hereditary spastic paraplegias (HSPs) are a diverse class of neurodegenerative disorders that mainly affect the corticospinal tract of the body and result in various clinical conditions such as lower limb spasticity and muscle weakness in the lower extremities. Worldwide, more than 70 chromosomal loci/genes have been reported to be associated with HSPs, out of which, six genes viz., ATL1, FA2H, GJC2, AP4E1, ALDH18A1 and ATP13A2 have been mapped in Pakistani families. In the present genetic study, we report on a large consanguineous Pakistani family with a complex form of HSP segregating with a 18 bp deletion in the first exon of the Fatty Acid 2-Hydroxylase (FA2H) gene (NM_024306.5:c.159_176del). The identified in-frame deletion results in loss of six amino acids (p.Arg53_Ile58del) within the cytochrome B5 domain of the protein. FA2H is required for alpha-hydroxylation of free fatty acids to form alpha-hydroxylated sphingolipids. Its cytochrome b5-like heme-binding domain, which spans from residues 15 to 85, imparts the redox activity to FA2H. This mutation has previously been reported in a Pakistani family presenting with a similar form of complex HSP. Together with our findings the pathogenic role of the observed variant is further supported. Mutation studies on additional Pakistani families for FA2H will further elucidate its mutational spectrum, which may help in developing a prenatal diagnostic test for Khyber Pakhtunkhwa resident Pakistani families.

Modulation of CYP2C9 activity and hydrogen peroxide production by cytochrome b5.

Cytochromes P450 (CYP) play a major role in drug detoxification, and cytochrome b5 (cyt b5) stimulates the catalytic cycle of mono-oxygenation and detoxification reactions. Collateral reactions of this catalytic cycle can lead to a significant production of toxic reactive oxygen species (ROS). One of the most abundant CYP isoforms in the human liver is CYP2C9, which catalyzes the metabolic degradation of several drugs including nonsteroidal anti-inflammatory drugs. We studied modulation by microsomal membrane-bound and soluble cyt b5 of the hydroxylation of salicylic acid to gentisic acid and ROS release by CYP2C9 activity in human liver microsomes (HLMs) and by CYP2C9 baculosomes. CYP2C9 accounts for nearly 75% of salicylic acid hydroxylation in HLMs at concentrations reached after usual aspirin doses. The anti-cyt b5 antibody SC9513 largely inhibits the rate of salicylic acid hydroxylation by CYP2C9 in HLMs and CYP2C9 baculosomes, increasing the KM approximately threefold. Besides, soluble human recombinant cyt b5 stimulates the Vmax nearly twofold while it decreases nearly threefold the Km value in CYP2C9 baculosomes. Regarding NADPH-dependent ROS production, soluble recombinant cyt b5 is a potent inhibitor both in HLMs and in CYP2C9 baculosomes, with inhibition constants of 1.04 ± 0.25 and 0.53 ± 0.06 µM cyt b5, respectively. This study indicates that variability in cyt b5 might be a major factor underlying interindividual variability in the metabolism of CYP2C9 substrates.

Heterologous expression of high-activity cytochrome P450 in mammalian cells.

The evaluation of Cytochrome P450 (CYP) enzymatic activity is essential to estimate drug pharmacokinetics. Numerous CYP allelic variants have been identified; the functional characterisation of these variants is required for their application in precision medicine. Results from heterologous expression systems using mammalian cells can be integrated in in vivo studies; however, other systems such as E. coli, bacteria, yeast, and baculoviruses are generally used owing to the difficulty in expressing high CYP levels in mammalian cells. Here, by optimising transfection and supplementing conditions, we developed a heterologous expression system using 293FT cells to evaluate the enzymatic activities of three CYP isoforms (CYP1A2, CYP2C9, and CYP3A4). Moreover, we established co-expression with cytochrome P450 oxidoreductase and cytochrome b5. This expression system would be a potential complementary or beneficial alternative approach for the pharmacokinetic evaluation of clinically used and developing drugs in vitro.