Structure and function of cytochromes P450: a comparative analysis of three crystal structures

The Dark Side of the pollen: BSA-seq identified genomic regions linked to male sterility in globe artichoke

Globe artichoke (Cynara cardunculus var. scolymus; 2n = 2x = 34) is a food crop consumed for its immature flower heads. Traditionally, globe artichoke varietal types are vegetatively propagated. However, seed propagation makes it possible to treat the crop as annual, increasing field uniformity and reducing farmers costs, as well as pathogens diffusion. Despite globe artichoke's significant agricultural value and the critical role of heterosis in the development of superior varieties, the production of hybrids remains challenging without a reliable system for large-scale industrial seed production. Male sterility (MS) presents a promising avenue for overcoming these challenges by simplifying the hybridization process and enabling cost-effective seed production. However, within the Cynara genus, genic male sterility has been linked to three recessive loci in globe artichoke, with no definitive genetic mechanism elucidated to date. A 250 offsprings F2 population, derived from a cross between a MS globe artichoke and a male fertile (MF) cultivated cardoon (C. cardunculus var. altilis) and fitting a monogenic segregation model (3:1), was analyzed through BSA-seq, aiming at the identification of genomic regions/genes affecting male sterility. Four QTL regions were identified on chromosomes 4, 12, and 14. By analyzing the sequence around the highest pick on chromosome 14, a cytochrome P450 (CYP703A2) was identified, carrying a deleterious substitution (R/Q) fixed in the male sterile parent. A single dCAPS marker was developed around this SNP, allowing the discrimination between MS and MF genotypes within the population, suitable for applications in plant breeding programs. A 3D model of the protein was generated by homology modeling, revealing that the mutated amino acid is part of a highly conserved motif crucial for protein folding.

Disentangling hydroxynitrile glucoside biosynthesis in a barley (Hordeum vulgare) metabolon provides access to elite malting barleys for ethyl carbamate-free whisky production

Barley produces several specialized metabolites, including five α-, β-, and γ-hydroxynitrile glucosides (HNGs). In malting barley, presence of the α-HNG epiheterodendrin gives rise to undesired formation of ethyl carbamate in the beverage production, especially after distilling. Metabolite-GWAS identified QTLs and underlying gene candidates possibly involved in the control of the relative and absolute content of HNGs, including an undescribed MATE transporter. By screening 325 genetically diverse barley accessions, we discovered three H. vulgare ssp. spontaneum (wild barley) lines with drastic changes in the relative ratios of the five HNGs. Knock-out (KO)-lines, isolated from the barley FIND-IT resource and each lacking one of the functional HNG biosynthetic genes (CYP79A12, CYP71C103, CYP71C113, CYP71U5, UGT85F22 and UGT85F23) showed unprecedented changes in HNG ratios enabling assignment of specific and mutually dependent catalytic functions to the biosynthetic enzymes involved. The highly similar relative ratios between the five HNGs found across wild and domesticated barley accessions indicate assembly of the HNG biosynthetic enzymes in a metabolon, the functional output of which was reconfigured in the absence of a single protein component. The absence or altered ratios of the five HNGs in the KO-lines did not change susceptibility to the fungal phytopathogen Pyrenophora teres causing net blotch. The study provides a deeper understanding of the organization of HNG biosynthesis in barley and identifies a novel, single gene HNG-0 line in an elite spring barley background for direct use in breeding of malting barley, eliminating HNGs as a source of ethyl carbamate formation in whisky production.

Identification and characterization of camptothecin tailoring enzymes in Nothapodytes tomentosa

Camptothecin is a complex monoterpenoid indole alkaloid with remarkable antitumor activity. Given that two C-10 modified camptothecin derivatives, topotecan and irinotecan, have been approved as potent anticancer agents, there is a critical need for methods to access other aromatic ring-functionalized congeners (e.g., C-9, C-10, etc.). However, contemporary methods for chemical oxidation are generally harsh and low-yielding when applied to the camptothecin scaffold, thereby limiting the development of modified derivatives. Reported herein, we have identified four tailoring enzymes responsible for C-9 modifications of camptothecin from Nothapodytes tomentosa, via metabolomic and transcriptomic analysis. These consist of a cytochrome P450 (NtCPT9H) which catalyzes the regioselective oxidation of camptothecin to 9-hydroxycamptothecin, as well as two methyltransferases (NtOMT1/2, converting 9-hydroxycamptothecin to 9-methoxycamptothecin), and a uridine diphosphate-glycosyltransferase (NtUGT5, decorating 9-hydroxycamptothecin to 9-β-D-glucosyloxycamptothecin). Importantly, the critical residues that contribute to the specific catalytic activity of NtCPT9H have been elucidated through molecular docking and mutagenesis experiments. This work provides a genetic basis for producing camptothecin derivatives through metabolic engineering. This will hasten the discovery of novel C-9 modified camptothecin derivatives, with profound implications for pharmaceutical manufacture.

Genome-wide identification and expression analyses of CYP450 genes in sweet potato (Ipomoea batatas L.)

BACKGROUND

Cytochrome P450 monooxygenases (CYP450s) play a crucial role in various biochemical reactions involved in the synthesis of antioxidants, pigments, structural polymers, and defense-related compounds in plants. As sweet potato (Ipomoea batatas L.) holds significant economic importance, a comprehensive analysis of CYP450 genes in this plant species can offer valuable insights into the evolutionary relationships and functional characteristics of these genes.

RESULTS

In this study, we successfully identified and categorized 95 CYP450 genes from the sweet potato genome into 5 families and 31 subfamilies. The predicted subcellular localization results indicate that CYP450s are distributed in the cell membrane system. The promoter region of the IbCYP450 genes contains various cis-acting elements related to plant hormones and stress responses. In addition, ten conserved motifs (Motif1-Motif10) have been identified in the IbCYP450 family proteins, with 5 genes lacking introns and only one exon. We observed extensive duplication events within the CYP450 gene family, which may account for its expansion. The gene duplication analysis results showed the presence of 15 pairs of genes with tandem repeats. Interaction network analysis reveals that IbCYP450 families can interact with multiple target genes and there are protein-protein interactions within the family. Transcription factor interaction analysis suggests that IbCYP450 families interact with multiple transcription factors. Furthermore, gene expression analysis revealed tissue-specific expression patterns of CYP450 genes in sweet potatoes, as well as their response to abiotic stress and plant hormones. Notably, quantitative real-time polymerase chain reaction (qRT‒PCR) analysis indicated the involvement of CYP450 genes in the defense response against nonbiological stresses in sweet potatoes.

CONCLUSIONS

These findings provide a foundation for further investigations aiming to elucidate the biological functions of CYP450 genes in sweet potatoes.

Structural Characterization and Molecular Dynamics Study of the REPI Fusion Protein from Papaver somniferum L

REPI is a pivotal point enzyme in plant benzylisoquinoline alkaloid metabolism as it promotes the evolution of the biosynthetic branch of morphinan alkaloids. Experimental studies of its activity led to the identification of two modules (DRS and DRR) that catalyze two sequential steps of the epimerization of (S)- to (R)-reticuline. Recently, special attention has been paid to its genetic characterization and evolutionary history, but no structural analyses of the REPI protein have been conducted to date. We present here a computational structural characterization of REPI with heme and NADP cofactors in the apo state and in three complexes with substrate (S)-reticuline in DRS and intermediate 1,2-dehydroreticuline in DRS and in DRR. Since no experimental structure exists for REPI, we used its AlphaFold model as a scaffold to build up these four systems, which were submitted to all-atom molecular dynamics (MD) simulations. A comparison of MD results for the four systems revealed key dynamic changes associated with cofactor and ligand binding and provided a dynamic picture of the evolution of their structures and interactions. We also explored the possible dynamic occurrence of tunnels and electrostatic highways potentially involved in alternative mechanisms for channeling the intermediate from DRS to DRR.

Computational Exploration of Minimum Energy Reaction Pathway of N2O Formation from Intermediate I of P450nor Using an Active Center Model

P450nor is a heme-containing enzyme that catalyzes the conversion of nitric oxide (NO) to nitrous oxide (N2O). Its catalytic mechanism has attracted attention in chemistry, biology, and environmental engineering. The catalytic cycle of P450nor is proposed to consist of three major steps. The reaction mechanism for the last step, N2O generation, remains unknown. In this study, the reaction pathway of the N2O generation from the intermediate I was explored with the B3LYP calculations using an active center model after the examination of the validity of the model. In the validation, we compared the heme distortions between P450nor and other oxidoreductases, suggesting a small effect of protein environment on the N2O generation reaction in P450nor. We then evaluated the electrostatic environment effect of P450nor on the hydride affinity to the active site with quantum mechanics/molecular mechanics (QM/MM) calculations, confirming that the affinity was unchanged with or without the protein environment. The active center model for P450nor showed that the N2O generation process in the enzymatic reaction undergoes a reasonable barrier height without protein environment. Consequently, our findings strongly suggest that the N2O generation reaction from the intermediate I depends sorely on the intrinsic reactivity of the heme cofactor bound on cysteine residue.

Functional Characterization of 29 Cytochrome P450 4F2 Variants Identified in a Population of 8380 Japanese Subjects and Assessment of Arachidonic Acid ω-Hydroxylation

Cytochrome P450 4F2 (CYP4F2) is an enzyme that is involved in the metabolism of arachidonic acid (AA), vitamin E and K, and xenobiotics including drugs. CYP4F2*3 polymorphism (rs2108622; c.1297G>A; p.Val433Met) has been associated with hypertension, ischemic stroke, and variation in the effectiveness of the anticoagulant drug warfarin. In this study, we characterized wild-type CYP4F2 and 28 CYP4F2 variants, including a Val433Met substitution, detected in 8380 Japanese subjects. The CYP4F2 variants were heterologously expressed in 293FT cells to measure the concentrations of CYP4F2 variant holoenzymes using carbon monoxide-reduced difference spectroscopy, where the wild type and 18 holoenzyme variants showed a peak at 450 nm. Kinetic parameters [Vmax , substrate concentration producing half of Vmax (S50 ), and intrinsic clearance (CL int ) as Vmax /S50 ] of AA ω-hydroxylation were determined for the wild type and 21 variants with enzyme activity. Compared with the wild type, two variants showed significantly decreased CL int values for AA ω-hydroxylation. The values for seven variants could not be determined because no enzymatic activity was detected at the highest substrate concentration used. Three-dimensional structural modeling was performed to determine the reason for reduced enzymatic activity of the CYP4F2 variants. Our findings contribute to a better understanding of CYP4F2 variant-associated diseases and possible future therapeutic strategies. SIGNIFICANCE STATEMENT: CYP4F2 is involved in the metabolism of arachidonic acid and vitamin K, and CYP4F2*3 polymorphisms have been associated with hypertension and variation in the effectiveness of the anticoagulant drug warfarin. This study presents a functional analysis of 28 CYP4F2 variants identified in Japanese subjects, demonstrating that seven gene polymorphisms cause loss of CYP4F2 function, and proposes structural changes that lead to altered function.

CYP26B1-related disorder: expanding the ends of the spectrum through clinical and molecular evidence

CYP26B1 metabolizes retinoic acid in the developing embryo to regulate its levels. A limited number of individuals with pathogenic variants in CYP26B1 have been documented with a varied phenotypic spectrum, spanning from a severe manifestation involving skull anomalies, craniosynostosis, encephalocele, radio-humeral fusion, oligodactyly, and a narrow thorax, to a milder presentation characterized by craniosynostosis, restricted radio-humeral joint mobility, hearing loss, and intellectual disability. Here, we report two families with CYP26B1-related phenotypes and describe the data obtained from functional studies of the variants. Exome and Sanger sequencing were used for variant identification in family 1 and family 2, respectively. Family 1 reflects a mild phenotype, which includes craniofacial dysmorphism with brachycephaly (without craniosynostosis), arachnodactyly, reduced radioulnar joint movement, conductive hearing loss, learning disability-and compound heterozygous CYP26B1 variants: (p.[(Pro118Leu)];[(Arg234Gln)]) were found. In family 2, a stillborn fetus presented a lethal phenotype with spina bifida occulta, hydrocephalus, poor skeletal mineralization, synostosis, limb defects, and a synonymous homozygous variant in CYP26B1: c.1083C > A. A minigene assay revealed that the synonymous variant created a new splice site, removing part of exon 5 (p.Val361_Asp382del). Enzymatic activity was assessed using a luciferase assay, demonstrating a notable reduction in exogenous retinoic acid metabolism for the variant p.Val361_Asp382del. (~ 3.5 × decrease compared to wild-type); comparatively, the variants p.(Pro118Leu) and p.(Arg234Gln) demonstrated a partial loss of metabolism (1.7× and 2.3× reduction, respectively). A proximity-dependent biotin identification assay reaffirmed previously reported ER-resident protein interactions. Additional work into these interactions is critical to determine if CYP26B1 is involved with other biological events on the ER. Immunofluorescence assay suggests that mutant CYP26B1 is still localized in the endoplasmic reticulum. These results indicate that novel pathogenic variants in CYP26B1 result in varying levels of enzymatic activity that impact retinoic acid metabolism and relate to the distinct phenotypes observed.

Two CYP93A enzymes play a dual role in isoflavonoid biosynthesis in Glycine max L

Glycine max L. is rich in isoflavonoids with diverse biological activities. However, isoflavonoid biosynthetic pathway is not fully elucidated in soybean. In the present study, we investigated characteristics of all the thirteen CYP93 subfamily members, and found GmCYP93A1, GmCYP93A2, and GmCYP93A3 are closely clustered, preferentially expressed in roots, and highly inducible by elicitor. When expressed in yeast, GmCYP93A1 was active towards liquiritigenin, naringenin, and 3,9-dihydroxyptercarpan, GmCYP93A2 towards 3,9-dihydroxyptercarpan with strict substrate specificity, whereas GmCYP93A3 did not show any activity towards all the tested substrates. Both GmCYP93A1 and GmCYP93A2 could catalyze 3,9-dihydroxyptercarpan into daidzein and glycinol, with both hydroxylation and aryl migration activity. Site-directed mutagenesis assays revealed that mutation in Thr446 to Ser446 in heme-binding domain increased the enzyme activity of GmCYP93A1 towards 3,9-dihydroxyptercarpan, which highlights its key amino acid residues as shown with its molecular docking with 3,9-dihydroxyptercarpan and HEM. Overexpression of GmCYP93A1 and GmCYP93A2 in the soybean hairy roots reduced the content of daidzein, whereas knockdown of these two genes increased genistein content, indicating changes in expression level of GmCYP93A1 and GmCYP93A2 altered isoflavonoid flux in soybean. Our studies on the activity of GmCYP93A1 and GmCYP93A2 enriched diverse functions of CYP93 subfamily in soybean isoflavonoid pathway, which is valuable for further understanding and bioengineering of isoflavonoid pathway in soybean.

Biological response and cell death signaling pathways modulated by tetrahydroisoquinoline-based aldoximes in human cells

The uncharged 3-hydroxy-2-pyridine aldoximes with protonatable tertiary amines are studied as antidotes in toxic organophosphates (OP) poisoning. Due to some of their specific structural features, we hypothesize that these compounds could exert diverse biological activity beyond their main scope of application. To examine this further, we performed an extensive cell-based assessment to determine their effects on human cells (SH-SY5Y, HEK293, HepG2, HK-2, myoblasts and myotubes) and possible mechanism of action. As our results indicated, aldoxime having a piperidine moiety did not induce significant toxicity up to 300µM within 24hours, while those with a tetrahydroisoquinoline moiety, in the same concentration range, showed time-dependent effects and stimulated mitochondria-mediated activation of the intrinsic apoptosis pathway through ERK1/2 and p38-MAPK signaling and subsequent activation of initiator caspase 9 and executive caspase 3 accompanied with DNA damage as observed already after 4hour exposure. Mitochondria and fatty acid metabolism were also likely targets of 3-hydroxy-2-pyridine aldoximes with tetrahydroisoquinoline moiety, due to increased phosphorylation of acetyl-CoA carboxylase. In silico analysis predicted kinases as their most probable target class, while pharmacophores modeling additionally predicted the inhibition of a cytochrome P450cam. Overall, if the absence of significant toxicity for piperidine bearing aldoxime highlights the potential of its further studies in medical counter-measures, the observed biological activity of aldoximes with tetrahydroisoquinoline moiety could be indicative for future design of compounds either in a negative context in OP antidotes design, or in a positive one for design of compounds for the treatment of other phenomena like cell proliferating malignancies.

Genome-Wide Identification and Expression Profiling of Cytochrome P450 Monooxygenase Superfamily in Foxtail Millet

The cytochrome P450 monooxygenases (CYP450) are the largest enzyme family in plant metabolism and widely involved in the biosynthesis of primary and secondary metabolites. Foxtail millet (Setaria italica (L.) P. Beauv) can respond to abiotic stress through a highly complex polygene regulatory network, in which the SiCYP450 family is also involved. Although the CYP450 superfamily has been systematically studied in a few species, the research on the CYP450 superfamily in foxtail millet has not been completed. In this study, three hundred and thirty-one SiCYP450 genes were identified in the foxtail millet genome by bioinformatics methods, which were divided into four groups, including forty-six subgroups. One hundred and sixteen genes were distributed in thirty-three tandem duplicated gene clusters. Chromosome mapping showed that SiCYP450 was distributed on seven chromosomes. In the SiCYP450 family of foxtail millet, 20 conserved motifs were identified. Cis-acting elements in the promoter region of SiCYP450 genes showed that hormone response elements were found in all SiCYP450 genes. Of the three hundred and thirty-one SiCYP450 genes, nine genes were colinear with the Arabidopsis thaliana genes. Two hundred SiCYP450 genes were colinear with the Setaria viridis genes, including two hundred and forty-five gene duplication events. The expression profiles of SiCYP450 genes in different organs and developmental stages showed that SiCYP450 was preferentially expressed in specific tissues, and many tissue-specific genes were identified, such as SiCYP75B6, SiCYP96A7, SiCYP71A55, SiCYP71A61, and SiCYP71A62 in the root, SiCYP78A1 and SiCYP94D9 in leaves, and SiCYP78A6 in the ear. The RT-PCR data showed that SiCYP450 could respond to abiotic stresses, ABA, and herbicides in foxtail millet. Among them, the expression levels of SiCYP709B4, SiCYP71A11, SiCYP71A14, SiCYP78A1, SiCYP94C3, and SiCYP94C4 were significantly increased under the treatment of mesotrione, florasulam, nicosulfuron, fluroxypyr, and sethoxydim, indicating that the same gene might respond to multiple herbicides. The results of this study will help reveal the biological functions of the SiCYP450 family in development regulation and stress response and provide a basis for molecular breeding of foxtail millet.

Molecular cloning, expression analysis and functional characterization of chicken cytochrome P450 27A1: A novel mitochondrial vitamin D3 25-hydroxylase

Vitamin D₃ is hydroxylated by cytochrome P450 (CYP) before exerting biological effects. The chicken CYP involved in vitamin D₃ 25-hydroxylation has yet to be cloned, and little is known about its functional characteristics, tissue distribution, and cellular expression. We identified a novel, full-length CYP27A1 gene cloned from chicken hepatocyte cDNA that encodes a putative protein of 518 amino acids. Swiss modeling revealed that chicken CYP27A1 has a classic open-fold form. Multisequence homology alignment determined that CYP27A1 contains conserved motifs for substrate recognition and binding. Quantitative real-time PCR analysis in 2-mo-old Partridge Shank broilers demonstrated that CYP27A1 mRNA levels were highest in the liver, followed by the thigh muscles, the breast muscles, and kidneys. The transcripts of CYP27A1 in breast muscles were significantly higher in males than in females. A subcellular localization analysis demonstrated that CYP27A1 was mainly expressed in the mitochondria. In vitro enzyme assays suggested that recombinant CYP27A1 hydroxylates vitamin D₃ at the C-25 position to form 25-hydroxyvitamin D₃ (25(OH)D₃). The K_m and V_max values for CYP27A1-dependent vitamin D₃ 25-hydroxylation were estimated to be 4.929 μM and 0.389 mol min^-1 mg^-1 protein, respectively. In summary, these results suggest that CYP27A1 encodes a mitochondrial CYP that plays an important physiologic role in the 25-hydroxylation of vitamin D₃ in chickens, providing novel insights into vitamin D₃ metabolism in this species.

Engineering Nicotiana tabacum trichomes for triterpenic acid production

In this work, we aimed at implementing the biosynthesis of triterpenic acids in Nicotiana tabacum glandular trichomes. Although endogenous genes coding for enzymes involved in such biosynthetic pathway are found in the Nicotiana tabacum genome, implementing such pathway specifically in glandular trichomes required to boost endogenous enzymatic activities. Five transgenes coding for a farnesyl-diphosphate synthase, a squalene synthase, a squalene epoxidase, a beta-amyrin synthase and a beta-amyrin 28-monooxygenase were introduced in N.tabacum, their expression being driven by pMALD1, a trichome-specific transcriptional promoter. This study aimed at testing whether sinking isoprenoid precursors localized in plastids, by exploiting potential cross-talks allowing the exchange of terpenoid pools from the chloroplast to the cytosol, could be a way to improve overall yield. By analyzing metabolites extracted from entire leaves, a low amount of ursolic acid was detected in plants expressing the five transgenes. Our study shows that the terpene biosynthetic pathway could be, in part, redirected in N.tabacum glandular trichomes with no deleterious phenotype at the whole plant level (chlorosis, dwarfism,…). In light of our results, possible ways to improve the final yield are discussed.

Genes encoding cytochrome P450 monooxygenases and glutathione S-transferases associated with herbicide resistance evolved before the origin of land plants

Cytochrome P450 (CYP) monooxygenases and glutathione S-transferases (GST) are enzymes that catalyse chemical modifications of a range of organic compounds. Herbicide resistance has been associated with higher levels of CYP and GST gene expression in some herbicide-resistant weed populations compared to sensitive populations of the same species. By comparing the protein sequences of 9 representative species of the Archaeplastida-the lineage which includes red algae, glaucophyte algae, chlorophyte algae, and streptophytes-and generating phylogenetic trees, we identified the CYP and GST proteins that existed in the common ancestor of the Archaeplastida. All CYP clans and all but one land plant GST classes present in land plants evolved before the divergence of streptophyte algae and land plants from their last common ancestor. We also demonstrate that there are more genes encoding CYP and GST proteins in land plants than in algae. The larger numbers of genes among land plants largely results from gene duplications in CYP clans 71, 72, and 85 and in the GST phi and tau classes [1,2]. Enzymes that either metabolise herbicides or confer herbicide resistance belong to CYP clans 71 and 72 and the GST phi and tau classes. Most CYP proteins that have been shown to confer herbicide resistance are members of the CYP81 family from clan 71. These results demonstrate that the clan and class diversity in extant plant CYP and GST proteins had evolved before the divergence of land plants and streptophyte algae from a last common ancestor estimated to be between 515 and 474 million years ago. Then, early in embryophyte evolution during the Palaeozoic, gene duplication in four of the twelve CYP clans, and in two of the fourteen GST classes, led to the large numbers of CYP and GST proteins found in extant land plants. It is among the genes of CYP clans 71 and 72 and GST classes phi and tau that alleles conferring herbicide resistance evolved in the last fifty years.

Genetics of tolerance in honeybees to the neonicotinoid clothianidin

The effects of neonicotinoid insecticides (NNIs) on honeybee health are intensely debated, with numerous studies showing negative effects of exposure, while others report no such effects. We carried out experiments to study the genetic and molecular basis of NNI tolerance in honeybees, which may underlie the discrepancies observed in the literature. We discovered that worker survival post-exposure to an acute oral dose of clothianidin is heritable (H ² = 37.8%). Tolerance to clothianidin was not associated with differences in the expression of detoxification enzymes in our experiments. Instead, mutations in the primary neonicotinoid detoxification genes CYP9Q1 and CYP9Q3 were strongly associated with worker survival post-clothianidin exposure. In some instances, the strong association between CYP9Q haplotypes and worker survival was associated with the protein's predicted binding affinity for clothianidin. Our findings have implications regarding future toxicological studies utilizing honeybees as a model pollinator.

Genome-wide analysis of cytochrome P450 genes in Citrus clementina and characterization of a CYP gene encoding flavonoid 3'-hydroxylase

Cytochrome P450s (CYPs) are the largest family of enzymes in plant and play multifarious roles in development and defense but the available information about the CYP superfamily in citrus is very limited. Here we provide a comprehensive genome-wide analysis of the CYP superfamily in Citrus clementina genome, identifying 301 CYP genes grouped into ten clans and 49 families. The characteristics of both gene structures and motif compositions strongly supported the reliability of the phylogenetic relationship. Duplication analysis indicated that tandem duplication was the major driving force of expansion for this superfamily. Promoter analysis revealed numerous cis-acting elements related to various responsiveness. RNA-seq data elucidated their expression patterns in citrus fruit peel both during development and in response to UV-B. Furthermore, we characterize a UV-B-induced CYP gene (Ciclev10019637m, designated CitF3'H) as a flavonoid 3'-hydroxylase for the first time. CitF3'H catalyzed numerous flavonoids and favored naringenin in yeast assays. Virus-induced silencing of CitF3'H in citrus seedlings significantly reduced the levels of 3'-hydroxylated flavonoids and their derivatives. These results together with the endoplasmic reticulum-localization of CitF3'H in plant suggest that this enzyme is responsible for the biosynthesis of 3'-hydroxylated flavonoids in citrus. Taken together, our findings provide extensive information about the CYP superfamily in citrus and contribute to further functional verification.

Spotlight on CYP4B1

The mammalian cytochrome P450 monooxygenase CYP4B1 can bioactivate a wide range of xenobiotics, such as its defining/hallmark substrate 4-ipomeanol leading to tissue-specific toxicities. Similar to other members of the CYP4 family, CYP4B1 has the ability to hydroxylate fatty acids and fatty alcohols. Structural insights into the enigmatic role of CYP4B1 with functions in both, xenobiotic and endobiotic metabolism, as well as its unusual heme-binding characteristics are now possible by the recently solved crystal structures of native rabbit CYP4B1 and the p.E310A variant. Importantly, CYP4B1 does not play a major role in hepatic P450-catalyzed phase I drug metabolism due to its predominant extra-hepatic expression, mainly in the lung. In addition, no catalytic activity of human CYP4B1 has been observed owing to a unique substitution of an evolutionary strongly conserved proline 427 to serine. Nevertheless, association of CYP4B1 expression patterns with various cancers and potential roles in cancer development have been reported for the human enzyme. This review will summarize the current status of CYP4B1 research with a spotlight on its roles in the metabolism of endogenous and exogenous compounds, structural properties, and cancer association, as well as its potential application in suicide gene approaches for targeted cancer therapy.

New dual functional CYP450 gene involves in isoflavone biosynthesis in Glycine max L

Glycine max L. accumulates a large amount of isoflavonoid compounds, which is beneficial for plant defense, plant-microbe symbiotic interactions, and human health. Several CYP450 subfamily genes are involved in the flavonoid biosynthetic pathway in plants. In the present study, we found 24 CYP82 subfamily genes were differentially expressed in various tissues of soybean, in Phytophthora sojae-infected soybean varieties and in soybean hairy roots treated with cell wall glucan elicitor. Six of them (GmCYP82A2, GmCYP82A3, GmCYP82A4, GmCYP82A23, GmCYP82C20 and GmCYP82D26) were co-expressed with other known isoflavonoid pathway genes in soybean. Their enzymatic activity in yeast feeding assays showed that only GmCYP82D26 was able to convert naringenin to daidzein with both aryl migration and dehydration function. When GmCYP82D26 was over-expressed in soybean hairy roots, the contents of the two major isoflavonoid aglycones in soybean (daidzein and genistein) were reduced, but total flavonoids were not affected. When GmCYP82D26 was suppressed by RNAi in the hairy roots, daidzein content was decreased but genistein content was increased, with unchanged total flavonoid content. GmCYP82D26 was found to be localized in the endoplasmic reticulum at subcellular level when transiently expressed in tobacco leaf epidermis. GmCYP82D26 gene was preferentially expressed in roots, with low expression level in other tissues in soybean. Homology modeling and molecular docking showed that GmCYP82D26 could form hydrogen bond with both HEM and naringenin at C5-OH and C4 carbonyl. All these results indicated that GmCYP82D26 possesses new and dual enzymatic activity, which bridges the two branches (daidzein and genistein branch) of isoflavonoid pathway in soybean.

Functional characterization of Cinnamate 4-hydroxylase gene family in soybean (Glycine max)

Cinnamate 4-hydroxylase (C4H) is the first key cytochrome P450 monooxygenase (P450) enzyme in the phenylpropanoid pathway. It belongs to the CYP73 family of P450 superfamily, and catalyzes the conversion of trans-cinnamic acid to p-coumaric acid. Since p-coumaric acid serves as the precursor for the synthesis of a wide variety of metabolites involved in plant development and stress resistance, alteration in the expression of soybean C4H genes is expected to affect the downstream metabolite levels, and its ability to respond to stress. In this study, we identified four C4H genes in the soybean genome that are distributed into both class I and class II CYP73 family. GmC4H2, GmC4H14 and GmC4H20 displayed tissue- and developmental stage-specific gene expression patterns with their transcript accumulation at the highest level in root tissues. GmC4H10 appears to be a pseudogene as its transcript was not detected in any soybean tissues. Furthermore, protein homology modelling revealed substrate docking only for GmC4H2, GmC4H14 and GmC4H20. To demonstrate the function of GmC4Hs, we modified a cloning vector for the heterologous expression of P450s in yeast, and used it for microsomal protein production and enzyme assay. Our results confirmed that GmC4H2, GmC4H14 and GmC4H20 contain the ability to hydroxylate trans-cinnamic acid with varying efficiencies.

Structural comparison of the cytochrome P450 enzymes CYP106A1 and CYP106A2 provides insight into their differences in steroid conversion

Understanding the structural basis of the selectivity of steroid hydroxylation requires detailed structural and functional investigations on various steroid hydroxylases with different selectivities, such as the bacterial cytochrome P450 enzymes. Here, the crystal structure of the cytochrome P450 CYP106A1 from Priestia megaterium was solved. CYP106A1 exhibits a rare additional structural motif of a cytochrome P450, a sixth β-sheet. The protein was found in different unusual conformations corresponding to both open and closed forms even when crystallized without any known substrate. The structural comparison of CYP106A1 with the previously investigated CYP106A2, including docking studies for both isoforms with the substrate cortisol, reveals a completely different orientation of the steroid molecule in the active sites. This distinction convincingly explains the experimentally observed differences in substrate conversion and product formation by the two enzymes.

Potential of cytochrome P450, a family of xenobiotic metabolizing enzymes, in cancer therapy

SIGNIFICANCE

Targeted cancer therapy with minimal off-target consequences has shown promise for some cancer types. Although cytochrome P450 (CYP) consists of 18 families, CYP1-4 families play key role in metabolizing xenobiotics and cancer drugs. This eventually affects the process of carcinogenesis, treatment outcome, and cancer drug resistance. Differential overexpression of CYPs in transformed cells, together with phenotypic alterations in tumors, presents a potential for therapeutic intervention.

RECENT ADVANCES

Recent advances in molecular tools and information technology have helped utilize CYPs as cancer targets. The precise expression in various tumors, X-ray crystal structures, improved understanding of the structure-activity relationship, and new approaches in the development of prodrugs have supported the ongoing efforts to develop CYPs-based drugs with a better therapeutic index.

CRITICAL ISSUES

Narrow therapeutic index, off-target effects, drug resistance, and tumor heterogeneity limit the benefits of CYP-based conventional cancer therapies. In this review, we address the CYP1-4 families as druggable targets in cancer. An emphasis is given to the CYP expression, function, and the possible mechanisms that drive expression and activity in normal and transformed tissues. The strategies that inhibit or activate CYPs for therapeutic benefits are also discussed.

FUTURE DIRECTIONS

Efforts are needed to develop more selective tools that will help comprehend molecular and metabolic alterations in tumor tissues with biological end-points in relation to CYPs. This will eventually translate to developing more specific CYP inhibitors/inducers.

Transgenic Expression of Haemonchus contortus Cytochrome P450 Hco-cyp-13A11 Decreases Susceptibility to Particular but Not All Macrocyclic Lactones in the Model Organism Caenorhabditis elegans

The number of reported macrocyclic lactones (ML) resistance cases across all livestock hosts is steadily increasing. Different studies in the parasitic nematode Haemonchus contortus assume the participation of cytochrome P450s (Cyps) enzymes in ML resistance. Still, functional data about their individual contribution to resistance or substrate specificity is missing. Via microinjection, transgenic Caenorhabditis elegans expressing HCON_00141052 (transgene-Hco-cyp-13A11) from extrachromosomal arrays were generated. After 24 h of exposure to different concentrations of ivermectin (IVM), ivermectin aglycone (IVMa), selamectin (SEL), doramectin (DRM), eprinomectin (EPR), and moxidectin (MOX), motility assays were performed to determine the impact of the H. contortus Cyp to the susceptibility of the worms against each ML. While transgene-Hco-cyp-13A11 significantly decreased susceptibility to IVM (four-fold), IVMa (2-fold), and SEL (3-fold), a slight effect for DRM and no effect for MOX, and EPR was observed. This substrate specificity of Hco-cyp-13A11 could not be explained by molecular modeling and docking studies. Hco-Cyp-13A11 molecular models were obtained for alleles from isolates with different resistance statuses. Although 14 amino acid polymorphisms were detected, none was resistance specific. In conclusion, Hco-cyp-13A11 decreased IVM, IVMa, and SEL susceptibility to a different extent, but its potential impact on ML resistance is not driven by polymorphisms.

Biosynthesis of Argolaphos Illuminates the Unusual Biochemical Origins of Aminomethylphosphonate and Nε-Hydroxyarginine Containing Natural Products

Phosphonate natural products have a history of successful application in medicine and biotechnology due to their ability to inhibit essential cellular pathways. This has inspired efforts to discover phosphonate natural products by prioritizing microbial strains whose genomes encode uncharacterized biosynthetic gene clusters (BGCs). Thus, success in genome mining is dependent on establishing the fundamental principles underlying the biosynthesis of inhibitory chemical moieties to facilitate accurate prediction of BGCs and the bioactivities of their products. Here, we report the complete biosynthetic pathway for the argolaphos phosphonopeptides. We uncovered the biochemical origins of aminomethylphosphonate (AMPn) and N^ε-hydroxyarginine, two noncanonical amino acids integral to the antimicrobial function of argolaphos. Critical to this pathway were dehydrogenase and transaminase enzymes dedicated to the conversion of hydroxymethylphosphonate to AMPn. The interconnected activities of both enzymes provided a solution to overcome unfavorable energetics, empower cofactor regeneration, and mediate intermediate toxicity during these transformations. Sequential ligation of l-arginine and l-valine was afforded by two GCN5-related N-acetyltransferases in a tRNA-dependent manner. AglA was revealed to be an unusual heme-dependent monooxygenase that hydroxylated the N^ε position of AMPn-Arg. As the first biochemically characterized member of the YqcI/YcgG protein family, AglA enlightens the potential functions of this elusive group, which remains biochemically distinct from the well-established P450 monooxygenases. The widespread distribution of AMPn and YqcI/YcgG genes among actinobacterial genomes suggests their involvement in diverse metabolic pathways and cellular functions. Our findings illuminate new paradigms in natural product biosynthesis and realize a significant trove of AmPn and N^ε-hydroxyarginine natural products that await discovery.

Second Sphere Effects on Oxygen Reduction and Peroxide Activation by Mononuclear Iron Porphyrins and Related Systems

Activation and reduction of O₂ and H₂O₂ by synthetic and biosynthetic iron porphyrin models have proved to be a versatile platform for evaluating second-sphere effects deemed important in naturally occurring heme active sites. Advances in synthetic techniques have made it possible to install different functional groups around the porphyrin ligand, recreating artificial analogues of the proximal and distal sites encountered in the heme proteins. Using judicious choices of these substituents, several of the elegant second-sphere effects that are proposed to be important in the reactivity of key heme proteins have been evaluated under controlled environments, adding fundamental insight into the roles played by these weak interactions in nature. This review presents a detailed description of these efforts and how these have not only demystified these second-sphere effects but also how the knowledge obtained resulted in functional mimics of these heme enzymes.

Genome-Wide Identification and Expression Profiles of 13 Key Structural Gene Families Involved in the Biosynthesis of Rice Flavonoid Scaffolds

Flavonoids are a class of key polyphenolic secondary metabolites with broad functions in plants, including stress defense, growth, development and reproduction. Oryza sativa L. (rice) is a well-known model plant for monocots, with a wide range of flavonoids, but the key flavonoid biosynthesis-related genes and their molecular features in rice have not been comprehensively and systematically characterized. Here, we identified 85 key structural gene candidates associated with flavonoid biosynthesis in the rice genome. They belong to 13 families potentially encoding chalcone synthase (CHS), chalcone isomerase (CHI), flavanone 3-hydroxylase (F3H), flavonol synthase (FLS), leucoanthocyanidin dioxygenase (LDOX), anthocyanidin synthase (ANS), flavone synthase II (FNSII), flavanone 2-hydroxylase (F2H), flavonoid 3′-hydroxylase (F3′H), flavonoid 3′,5′-hydroxylase (F3′5′H), dihydroflavonol 4-reductase (DFR), anthocyanidin reductase (ANR) and leucoanthocyanidin reductase (LAR). Through structural features, motif analyses and phylogenetic relationships, these gene families were further grouped into five distinct lineages and were examined for conservation and divergence. Subsequently, 22 duplication events were identified out of a total of 85 genes, among which seven pairs were derived from segmental duplication events and 15 pairs were from tandem duplications, demonstrating that segmental and tandem duplication events play important roles in the expansion of key flavonoid biosynthesis-related genes in rice. Furthermore, these 85 genes showed spatial and temporal regulation in a tissue-specific manner and differentially responded to abiotic stress (including six hormones and cold and salt treatments). RNA-Seq, microarray analysis and qRT-PCR indicated that these genes might be involved in abiotic stress response, plant growth and development. Our results provide a valuable basis for further functional analysis of the genes involved in the flavonoid biosynthesis pathway in rice.

Overexpression of the Panax ginseng CYP703 Alters Cutin Composition of Reproductive Tissues in Arabidopsis

Cytochrome P450 (CYP) catalyzes a wide variety of monooxygenation reactions in plant primary and secondary metabolisms. Land plants contain CYP703, belonging to the CYP71 clan, which catalyzes the biochemical pathway of fatty acid hydroxylation, especially in male reproductive tissues. Korean/Asian ginseng (Panax ginseng Meyer) has been regarded as one of important medicinal plant for a long time, however the molecular mechanism is less known on its development. In this study, we identified and characterized a CYP703A gene in P. ginseng (PgCYP703A4), regarding reproductive development. PgCYP703A4 shared a high-sequence identity (81–83%) with predicted amino acid as CYP703 in Dancus carota, Pistacia vera, and Camellia sinensis as well as 76% of amino acid sequence identity with reported CYP703 in Arabidopsis thaliana and 75% with Oryza sativa. Amino acid alignment and phylogenetic comparison of P. ginseng with higher plants and known A. thaliana members clearly distinguish the CYP703 members, each containing the AATDTS oxygen binding motif and PERH as a clade signature. The expression of PgCYP704B1 was only detected in P. ginseng flower buds, particularly in meiotic cells and the tapetum layer of developing anther, indicating the conserved role on male reproduction with At- and Os- CYP703. To acquire the clue of function, we transformed the PgCYP703A4 in A. thaliana. Independent overexpressing lines (PgCYP703A4ox) increased silique size and seed number, and altered the contents of fatty acids composition of cutin monomer in the siliques. Our results indicate that PgCYP703A4 is involved in fatty acid hydroxylation which affects cutin production and fruit size.

Cyanogenesis in the Sorghum Genus: From Genotype to Phenotype

Domestication has resulted in a loss of genetic diversity in our major food crops, leading to susceptibility to biotic and abiotic stresses linked with climate change. Crop wild relatives (CWR) may provide a source of novel genes potentially important for re-gaining climate resilience. Sorghum bicolor is an important cereal crop with wild relatives that are endemic to Australia. Sorghum bicolor is cyanogenic, but the cyanogenic status of wild Sorghum species is not well known. In this study, leaves of wild species endemic in Australia are screened for the presence of the cyanogenic glucoside dhurrin. The direct measurement of dhurrin content and the potential for dhurrin-derived HCN release (HCNp) showed that all the tested Australian wild species were essentially phenotypically acyanogenic. The unexpected low dhurrin content may reflect the variable and generally nutrient-poor environments in which they are growing in nature. Genome sequencing of six CWR and PCR amplification of the CYP79A1 gene from additional species showed that a high conservation of key amino acids is required for correct protein function and dhurrin synthesis, pointing to the transcriptional regulation of the cyanogenic phenotype in wild sorghum as previously shown in elite sorghum.

CYP11B1 variants influence skeletal maturation via alternative splicing

We performed genome-wide association study meta-analysis to identify genetic determinants of skeletal age (SA) deviating in multiple growth disorders. The joint meta-analysis (N = 4557) in two multiethnic cohorts of school-aged children identified one locus, CYP11B1 (expression confined to the adrenal gland), robustly associated with SA (rs6471570-A; β = 0.14; P = 6.2 × 10^-12). rs6410 (a synonymous variant in the first exon of CYP11B1 in high LD with rs6471570), was prioritized for functional follow-up being second most significant and the one closest to the first intron-exon boundary. In 208 adrenal RNA-seq samples from GTEx, C-allele of rs6410 was associated with intron 3 retention (P = 8.11 × 10^-40), exon 4 inclusion (P = 4.29 × 10^-34), and decreased exon 3 and 5 splicing (P = 7.85 × 10^-43), replicated using RT-PCR in 15 adrenal samples. As CYP11B1 encodes 11-β-hydroxylase, involved in adrenal glucocorticoid and mineralocorticoid biosynthesis, our findings highlight the role of adrenal steroidogenesis in SA in healthy children, suggesting alternative splicing as a likely underlying mechanism.

A Promiscuous Bacterial P450: The Unparalleled Diversity of BM3 in Pharmaceutical Metabolism

CYP102A1 (BM3) is a catalytically self-sufficient flavocytochrome fusion protein isolated from Bacillus megaterium, which displays similar metabolic capabilities to many drug-metabolizing human P450 isoforms. BM3's high catalytic efficiency, ease of production and malleable active site makes the enzyme a desirable tool in the production of small molecule metabolites, especially for compounds that exhibit drug-like chemical properties. The engineering of select key residues within the BM3 active site vastly expands the catalytic repertoire, generating variants which can perform a range of modifications. This provides an attractive alternative route to the production of valuable compounds that are often laborious to synthesize via traditional organic means. Extensive studies have been conducted with the aim of engineering BM3 to expand metabolite production towards a comprehensive range of drug-like compounds, with many key examples found both in the literature and in the wider industrial bioproduction setting of desirable oxy-metabolite production by both wild-type BM3 and related variants. This review covers the past and current research on the engineering of BM3 to produce drug metabolites and highlights its crucial role in the future of biosynthetic pharmaceutical production.

Harnessing P450 Enzyme for Biotechnology and Synthetic Biology

Cytochrome P450 enzymes (P450s, CYPs) catalyze the oxidative transformation of a wide range of organic substrates. Their functions are crucial to xenobiotic metabolism and steroid transformation in humans and other organisms. The enzymes are promising for synthetic biology applications but limited by several drawbacks including low turnover rates, poor stability, the dependance of expensive cofactors and redox partners, and the narrow substrate scope. To conquer these obstacles, emerging strategies including substrate engineering, usage of decoy and decoy-based small molecules auxiliaries, designing of artificial enzyme cascades and the incorporation of materials have been explored based on the unique properties of P450s. These strategies can be applied to a wide range of P450s and can be combined with protein engineering to improve the enzymatic activities. This minireview will focus on some recent developments of these strategies which have been used to leverage P450 catalysis. Remaining challenges and future opportunities will also be discussed.

Engineering Sequence and Selectivity of Late-Stage C-H Oxidation in the MycG Iterative Cytochrome P450

MycG is a multifunctional P450 monooxygenase that catalyzes sequential hydroxylation and epoxidation or a single epoxidation in mycinamicin biosynthesis. In the mycinamicin-producing strain Micromonospora griseorubida A11725, very low-level accumulation of mycinamicin V generated by the initial C-14 allylic hydroxylation of MycG is observed due to its subsequent epoxidation to generate mycinamicin II, the terminal metabolite in this pathway. Herein, we investigated whether MycG can be engineered for production of the mycinamicin II intermediate as the predominant metabolite. Thus, mycG was subject to random mutagenesis and screening was conducted in Escherichia coli whole-cell assays. This enabled efficient identification of amino acid residues involved in reaction profile alterations, which included MycG R111Q/V358L, W44R, and V135G/E355K with enhanced monohydroxylation to accumulate mycinamicin V. The MycG V135G/E355K mutant generated 40-fold higher levels of mycinamicin V compared to wild-type M. griseorubida A11725. In addition, the E355K mutation showed improved ability to catalyze sequential hydroxylation and epoxidation with minimal mono-epoxidation product mycinamicin I compared to the wild-type enzyme. These approaches demonstrate the ability to selectively coordinate the catalytic activity of multifunctional P450s and efficiently produce the desired compounds.

Genome-Wide Identification and Analyses of Drought/Salt-Responsive Cytochrome P450 Genes in Medicago truncatula

Cytochrome P450 monooxygenases (P450s) catalyze a great number of biochemical reactions and play vital roles in plant growth, development and secondary metabolism. As yet, the genome-scale investigation on P450s is still lacking in the model legume Medicago truncatula. In particular, whether and how many MtP450s are involved in drought and salt stresses for Medicago growth, development and yield remain unclear. In this study, a total of 346 MtP450 genes were identified and classified into 10 clans containing 48 families. Among them, sixty-one MtP450 genes pairs are tandem duplication events and 10 MtP450 genes are segmental duplication events. MtP450 genes within one family exhibit high conservation and specificity in intron–exon structure. Meanwhile, many Mt450 genes displayed tissue-specific expression pattern in various tissues. Specifically, the expression pattern of 204 Mt450 genes under drought/NaCl treatments were analyzed by using the weighted correlation network analysis (WGCNA). Among them, eight genes (CYP72A59v1, CYP74B4, CYP71AU56, CYP81E9, CYP71A31, CYP704G6, CYP76Y14, and CYP78A126), and six genes (CYP83D3, CYP76F70, CYP72A66, CYP76E1, CYP74C12, and CYP94A52) were found to be hub genes under drought/NaCl treatments, respectively. The expression levels of these selected hub genes could be induced, respectively, by drought/NaCl treatments, as validated by qPCR analyses, and most of these genes are involved in the secondary metabolism and fatty acid pathways. The genome-wide identification and co-expression analyses of M. truncatulaP450 superfamily genes established a gene atlas for a deep and systematic investigation of P450 genes in M. truncatula, and the selected drought-/salt-responsive genes could be utilized for further functional characterization and molecular breeding for resistance in legume crops.

Massively parallel characterization of CYP2C9 variant enzyme activity and abundance

CYP2C9 encodes a cytochrome P450 enzyme responsible for metabolizing up to 15% of small molecule drugs, and CYP2C9 variants can alter the safety and efficacy of these therapeutics. In particular, the anti-coagulant warfarin is prescribed to over 15 million people annually and polymorphisms in CYP2C9 can affect individual drug response and lead to an increased risk of hemorrhage. We developed click-seq, a pooled yeast-based activity assay, to test thousands of variants. Using click-seq, we measured the activity of 6,142 missense variants in yeast. We also measured the steady-state cellular abundance of 6,370 missense variants in a human cell line by using variant abundance by massively parallel sequencing (VAMP-seq). These data revealed that almost two-thirds of CYP2C9 variants showed decreased activity and that protein abundance accounted for half of the variation in CYP2C9 function. We also measured activity scores for 319 previously unannotated human variants, many of which may have clinical relevance.

Novel mutations of the CYP17A1 gene in four Chinese 46,XX cases with partial 17a-hydroxylase/17,20-lyase deficiency

The condition 17a-Hydroxylase/17,20-lyase deficiency (17-OHD) is a rare kind of congenital adrenal hyperplasia (CAH) characterized by failure to synthetize cortisol, adrenal androgens and gonadal steroids. Partial deficiency is much rarer, presenting with subtler symptoms. In this study, we summarized the clinical characteristics and identified the underlying gene mutation in four Chinese 46,XX patients with partial 17-OHD. Mutational analysis of the CYP17A1 gene was performed by polymerase chain reaction (PCR) and Sanger sequencing. Clinical and hormonal findings in these patients were consistent with typical manifestations of partial 17-OHD. All patients were found to have a compound heterozygous mutation of the CYP17A1 gene, with five mutations identified. Among them, c.887 T > C(p. I296T), c.1019G > A(p. R340H) and c.1346G > A(p. R449H) were novel missense mutations. In conclusion, we identified three novel missense mutations of the CYP17A1 gene from four patients with partial 17-OHD deficiency. Genotype-phenotype correlation analysis revealed that these novel mutations can lead to partial 17-OHD. Our findings thus provide novel insight into the clinical evaluations and molecular basis of 17-OHD.

Calculation of CYP450 protein-ligand binding and dissociation free energy paths

The function of an enzyme depends on its dynamic structure, and the catalytic mechanism has long been an active focus of research. The principle for interpreting protein selectivity and fidelity stems from optimization of the active site upon protein-substrate complexation, i.e., a lock-and-key configuration, on which most protein-substrate molecule binding recognition, and hence drug discovery, relies. Yet another thought has been to incorporate the protein folding interior tunnels for stereo- and regio-selectivity along the protein-substrate or protein-ligand/inhibitor binding process. Free energy calculations provide valuable information for molecular recognition and protein-ligand binding dynamics and kinetics. In this study, we focused on the kinetics of cytochrome P450 proteins (CYP450s) and the protein interior tunnel structure-dynamics relationship in terms of the substrate binding and leaving mechanism. A case in point is given by the prostaglandin H₂ (PGH₂) homologous isomerase of prostacyclin synthase. To calculate the reactant and product traversing the tunnels to and from the heme site, the free energy paths and tunnel potentials of mean force are constructed from steered molecular dynamics simulations and adaptive basing force umbrella sampling simulations. We explore the binding tunnels and critical residue lining characteristics for the ligand traverse and the underlying mechanism of CYP450 activity. Our theoretical analysis provides insights into the decisive role of the substrate tunnel binding process of the CYP450 mechanism and may be useful in drug design and protein engineering contexts.

Uncovering of cytochrome P450 anatomy by SecStrAnnotator

Protein structural families are groups of homologous proteins defined by the organization of secondary structure elements (SSEs). Nowadays, many families contain vast numbers of structures, and the SSEs can help to orient within them. Communities around specific protein families have even developed specialized SSE annotations, always assigning the same name to the equivalent SSEs in homologous proteins. A detailed analysis of the groups of equivalent SSEs provides an overview of the studied family and enriches the analysis of any particular protein at hand. We developed a workflow for the analysis of the secondary structure anatomy of a protein family. We applied this analysis to the model family of cytochromes P450 (CYPs)-a family of important biotransformation enzymes with a community-wide used SSE annotation. We report the occurrence, typical length and amino acid sequence for the equivalent SSE groups, the conservation/variability of these properties and relationship to the substrate recognition sites. We also suggest a generic residue numbering scheme for the CYP family. Comparing the bacterial and eukaryotic part of the family highlights the significant differences and reveals a well-known anomalous group of bacterial CYPs with some typically eukaryotic features. Our workflow for SSE annotation for CYP and other families can be freely used at address https://sestra.ncbr.muni.cz .

A Pharmacogenetic Study of CYP2C19 in Acute Coronary Syndrome Patients of Colombian Origin Reveals New Polymorphisms Potentially Related to Clopidogrel Therapy

Clopidogrel, an oral platelet P2Y₁₂ receptor blocker, is used in the treatment of acute coronary syndrome. Interindividual variability in treatment response and the occurrence of adverse effects has been attributed to genetic variants in CYP2C19. The analysis of relevant pharmacogenes in ethnically heterogeneous and poorly studied populations contributes to the implementation of personalized medicine. We analyzed the coding and regulatory regions of CYP2C19 in 166 patients with acute coronary syndrome (ACS) treated with clopidogrel. The allele frequencies of CYP2C19 alleles *1, *2, *4, *17, *27 and *33 alleles were 86.1%, 7.2%, 0.3%, 10.2%, 0.3% and 0.3%, respectively. A new potentially pathogenic mutation (p.L15H) and five intronic variants with potential splicing effects were detected. In 14.4% of the patients, a new haplotype in strong linkage disequilibrium was identified. The clinical outcome indicated that 13.5% of the patients presented adverse drugs reactions with a predominance of bleeding while 25% of these patients were carriers of at least one polymorphic allele. We propose that new regulatory single-nucleotide variants (SNVs) might potentially influence the response to clopidogrel in Colombian individuals.

The cytochrome P450 (CYP) superfamily in cnidarians

The cytochrome P450 (CYP) superfamily is a diverse and important enzyme family, playing a central role in chemical defense and in synthesis and metabolism of major biological signaling molecules. The CYPomes of four cnidarian genomes (Hydra vulgaris, Acropora digitifera, Aurelia aurita, Nematostella vectensis) were annotated; phylogenetic analyses determined the evolutionary relationships amongst the sequences and with existing metazoan CYPs. 155 functional CYPs were identified and 90 fragments. Genes were from 24 new CYP families and several new subfamilies; genes were in 9 of the 12 established metazoan CYP clans. All species had large expansions of clan 2 diversity, with H. vulgaris having reduced diversity for both clan 3 and mitochondrial clan. We identified potential candidates for xenobiotic metabolism and steroidogenesis. That each genome contained multiple, novel CYP families may reflect the large evolutionary distance within the cnidarians, unique physiology in the cnidarian classes, and/or different ecology of the individual species.

PCPD: Plant cytochrome P450 database and web-based tools for structural construction and ligand docking

Plant cytochrome P450s play key roles in the diversification and functional modification of plant natural products. Although over 200,000 plant P450 gene sequences have been recorded, only seven crystalized P450 genes severely hampered the functional characterization, gene mining and engineering of important P450s. Here, we combined Rosetta homologous modeling and MD-based refinement to construct a high-resolution P450 structure prediction process (PCPCM), which was applied to 181 plant P450s with identified functions. Furthermore, we constructed a ligand docking process (PCPLD) that can be applied for plant P450s virtual screening. 10 examples of virtual screening indicated the process can reduce about 80% screening space for next experimental verification. Finally, we constructed a plant P450 database (PCPD: http://p450.biodesign.ac.cn/), which includes the sequences, structures and functions of the 181 plant P450s, and a web service based on PCPCM and PCPLD. Our study not only developed methods for the P450-specific structure analysis, but also introduced a universal approach that can assist the mining and functional analysis of P450 enzymes.

Pervasive cooperative mutational effects on multiple catalytic enzyme traits emerge via long-range conformational dynamics

Multidimensional fitness landscapes provide insights into the molecular basis of laboratory and natural evolution. To date, such efforts usually focus on limited protein families and a single enzyme trait, with little concern about the relationship between protein epistasis and conformational dynamics. Here, we report a multiparametric fitness landscape for a cytochrome P450 monooxygenase that was engineered for the regio- and stereoselective hydroxylation of a steroid. We develop a computational program to automatically quantify non-additive effects among all possible mutational pathways, finding pervasive cooperative signs and magnitude epistasis on multiple catalytic traits. By using quantum mechanics and molecular dynamics simulations, we show that these effects are modulated by long-range interactions in loops, helices and β-strands that gate the substrate access channel allowing for optimal catalysis. Our work highlights the importance of conformational dynamics on epistasis in an enzyme involved in secondary metabolism and offers insights for engineering P450s.

Connecting conformational dynamics and epistasis has so far been limited to a few proteins and a single fitness trait. Here, the authors provide evidence of positive epistasis on multiple catalytic traits in the evolution and dynamics of engineered cytochrome P450 monooxygenase, offering insights for in silico protein design.

Identification of Wild-Type CYP321A2 and Comparison of Allelochemical-Induced Expression Profiles of CYP321A2 with Its Paralog CYP321A1 in Helicoverpa zea

One possible way to overcome the diversity of toxic plant allelochemicals idiosyncratically distributed among potential host plants is to have more counterdefense genes via gene duplication or fewer gene losses. Cytochrome P450 is the most important gene family responsible for detoxification of the diversity of plant allelochemicals. We have recently reported the identification and cloning of the transposon (HzSINE1)-disrupted non-functional CYP321A2, a duplicated paralog of the xenobiotic-metabolizing P450 CYP321A1 from a laboratory colony of Helicoverpa zea. Here we report the identification of the wild-type intact allele of CYP321A2 from another H. zea colony. This CYP321A2 allele encodes a deduced protein of 498 amino acids and has the P450 signature motifs. Quantitative RT-PCR experiments showed that this CYP321A2 allele was highly expressed in midgut and fat body and achieved the highest expression level in the developmental stage of 5th and 3rd instar larvae. CYP321A2 and CYP321A1 were constitutively expressed in low levels but can be differentially and significantly induced by a range of the plant allelochemicals and plant signal molecules, among which xanthotoxin, flavone, and coumarin were the most prominent inducers of CYP321A2 both in midgut and fat body, whereas flavone, coumarin, and indole-3-carbinol were the prominent inducers of CYP321A1 in midgut and fat body. Moreover, xanthotoxin- and flavone-responsive regulatory elements of CYP321A1 were also detected in the promoter region of CYP321A2. Our results enrich the P450 inventory by identifying an allelochemical broadly induced CYP321A2, a paralog of CYP321A1 in H. zea. Our data also suggest that the CYP321A2/CYP321A1 paralogs are a pair of duplicated genes of multigene families and CYP321A2 could potentially be involved in the detoxification of plant allelochemicals and adaptation of H. zea to its chemical environment.

Molecular and Structural Evolution of Cytochrome P450 Aromatase

Aromatase is the cytochrome P450 enzyme converting androgens into estrogen in the last phase of steroidogenesis. As estrogens are crucial in reproductive biology, aromatase is found in vertebrates and the invertebrates of the genus Branchiostoma, where it carries out the aromatization reaction of the A-ring of androgens that produces estrogens. Here, we investigate the molecular evolution of this unique and highly substrate-selective enzyme by means of structural, sequence alignment, and homology modeling, shedding light on its key role in species conservation. The alignments led to the identification of a core structure that, together with key and unique amino acids located in the active site and the substrate recognition sites, has been well conserved during evolution. Structural analysis shows what their roles are and the reason why they have been preserved. Moreover, the residues involved in the interaction with the redox partner and some phosphorylation sites appeared late during evolution. These data reveal how highly substrate-selective cytochrome P450 has evolved, indicating that the driving forces for evolution have been the optimization of the interaction with the redox partner and the introduction of phosphorylation sites that give the possibility of modulating its activity in a rapid way.

Functional Characterization of 40 CYP3A4 Variants by Assessing Midazolam 1'-Hydroxylation and Testosterone 6β-Hydroxylation

CYP3A4 is among the most abundant liver and intestinal drug-metabolizing cytochrome P450 enzymes, contributing to the metabolism of more than 30% of clinically used drugs. Therefore, interindividual variability in CYP3A4 activity is a frequent cause of reduced drug efficacy and adverse effects. In this study, we characterized wild-type CYP3A4 and 40 CYP3A4 variants, including 11 new variants, detected among 4773 Japanese individuals by assessing CYP3A4 enzymatic activities for two representative substrates (midazolam and testosterone). The reduced carbon monoxide-difference spectra of wild-type CYP3A4 and 31 CYP3A4 variants produced with our established mammalian cell expression system were determined by measuring the increase in maximum absorption at 450 nm after carbon monoxide treatment. The kinetic parameters of midazolam and testosterone hydroxylation by wild-type CYP3A4 and 29 CYP3A4 variants (K _m , k _cat , and catalytic efficiency) were determined, and the causes of their kinetic differences were evaluated by three-dimensional structural modeling. Our findings offer insight into the mechanism underlying interindividual differences in CYP3A4-dependent drug metabolism. Moreover, our results provide guidance for improving drug administration protocols by considering the information on CYP3A4 genetic polymorphisms. SIGNIFICANCE STATEMENT: CYP3A4 metabolizes more than 30% of clinically used drugs. Interindividual differences in drug efficacy and adverse-effect rates have been linked to ethnicity-specific differences in CYP3A4 gene variants in Asian populations, including Japanese individuals, indicating the presence of CYP3A4 polymorphisms resulting in the increased expression of loss-of-function variants. This study detected alterations in CYP3A4 activity due to amino acid substitutions by assessing the enzymatic activities of coding variants for two representative CYP3A4 substrates.

Advanced Understanding of the Electron Transfer Pathway of Cytochrome P450s

Cytochrome P450s are heme-thiolate enzymes that participate in carbon source assimilation, natural compound biosynthesis and xenobiotic metabolism in all kingdoms of life. P450s can catalyze various reactions by using a wide range of organic compounds, thus exhibiting great potential in biotechnological applications. The catalytic reactions of P450s are driven by electron equivalents that are sourced from pyridine nucleotides and delivered by cognate or matching redox partners (RPs). The electron transfer (ET) route from RPs to P450s involves one or more redox center-containing domains. As the rate of ET is one of the main determinants of P450 efficacy, an in-depth understanding of the P450 ET pathway should increase our knowledge of these important enzymes and benefit their further applications. Here, the various P450 RP systems along with current understanding of their ET routes will be reviewed. Notably, state-of-the-art structural studies of the two main types of self-sufficient P450 will also be summarized.

Cloning and characterization of estrogen hydroxylase (cyp1a1 and cyp1b1) genes in the stinging catfish Heteropneustes fossilis and induction of mRNA expression during final oocyte maturation

Estrogen hydroxylases (EHs) are cytochrome P450 Family 1 (Cyp1, Clan 2) proteins involved in estrogen hydroxylations at 2-, 4- or 16- carbon positions to form catecholestrogens. EHs are encoded by CYP1A1, CYP1A2 and CYP1B1 in mammals. In the catfish Heteropneustes fossilis, cyp1a1 and cyp1b1 cDNAs were cloned and characterized from liver and ovary. The cyp1a1 cDNA is 2071 bp long and codes for a 518 amino acids (aa) long protein. The cloned cyp1b1 cDNA is 1927 bp long and codes for a 509 residue protein. The deduced proteins clustered distinctly into teleost Cyp1a1 and Cyp1b1 clades, distinct from the tetrapod clusters and featured common function domains and homology with other teleost proteins. In the qPCR assay, the transcripts were the most abundant in the liver, followed by brain and ovary, and moderate in gill, kidney and muscle. Evidence was presented to show the involvement of the genes in reproduction. Expression of brain and ovarian transcripts showed significant seasonal variations with the highest abundance in the spawning phase. In situ hybridization showed the transcripts in the follicular layer (theca and granulosa) of the ovarian follicles. Periovulatory changes in the expression cyp1a1 and cyp1b1 were obtained during final oocyte maturation (FOM) and ovulation induced by human chorionic gonadotropin (hCG), both in vivo and in vitro, and by 2-hydroxyestradiol-17β (catecholestrogen) in vitro. In the brain, the transcript levels increased with time but in the ovary, the increase was maximal at 16 h and decreased at 24 h. The periovulatory activation of the cyp1 genes was reported in this study and discussed on the basis of complex regulation of FOM and ovulation.

A Review of the Important Role of CYP2D6 in Pharmacogenomics

Cytochrome P450 2D6 (CYP2D6) is a critical pharmacogene involved in the metabolism of ~20% of commonly used drugs across a broad spectrum of medical disciplines including psychiatry, pain management, oncology and cardiology. Nevertheless, CYP2D6 is highly polymorphic with single-nucleotide polymorphisms, small insertions/deletions and larger structural variants including multiplications, deletions, tandem arrangements, and hybridisations with non-functional CYP2D7 pseudogenes. The frequency of these variants differs across populations, and they significantly influence the drug-metabolising enzymatic function of CYP2D6. Importantly, altered CYP2D6 function has been associated with both adverse drug reactions and reduced drug efficacy, and there is growing recognition of the clinical and economic burdens associated with suboptimal drug utilisation. To date, pharmacogenomic clinical guidelines for at least 48 CYP2D6-substrate drugs have been developed by prominent pharmacogenomics societies, which contain therapeutic recommendations based on CYP2D6-predicted categories of metaboliser phenotype. Novel algorithms to interpret CYP2D6 function from sequencing data that consider structural variants, and machine learning approaches to characterise the functional impact of novel variants, are being developed. However, CYP2D6 genotyping is yet to be implemented broadly into clinical practice, and so further effort and initiatives are required to overcome the implementation challenges and deliver the potential benefits to the bedside.

Development of a method to determine cytochrome P450 1A2, 2C9, 2D6 and 3A4 activity sheep hepatic microsomes

INTRODUCTION

Ex vivo studies of human fetal hepatic drug metabolism are uncommon as it requires access to functional liver tissue and therefore raises practical and ethical concerns. Large animal models provide an alternative opportunity to study changes in cytochrome P450 (CYP) activity in the mother and fetus during pregnancy. We aimed to develop methods to determine the activity of CYP1A2, CYP2C9, CYP2D6 and CYP3A4 in sheep hepatic microsomes.

METHODS

We identified optimal conditions to determine the activity of CYP1A2 (using the probe drug phenacetin), CYP2C9 (diclofenac), CYP2D6 (dextromethorphan) and CYP3A4 (midazolam) by varying techniques for microsome extraction, probe drug concentration, incubation time and microsome concentration. The specificity of each probe drug was assessed by determining the rate of metabolism when specific CYP enzyme inhibitors were included in the reaction.

RESULTS

The optimum incubation time and probe drug concentration was six hours with 5 μM phenacetin (CYP1A2), four hours with 10 μM diclofenac (CYP2C9), 30 min with 1 μM of midazolam (CYP3A4) and 10 min with 1 μM dextromethorphan (CYP2D6). For both CYP2D6 and CYP3A4 reactions required 20 μg of microsomal protein, whereas for CYP1A2 and CYP2C9, reactions required 40 μg of microsomal protein. Metabolism of phenacetin, dextromethorphan and midazolam was reduced by specific enzyme inhibitors, but the specific CYP2C9 inhibitor sulfaphenazole did not substantially inhibit diclofenac metabolism.

DISCUSSION

This study identifies the optimal conditions for determining CYP activity in maternal sheep hepatic microsomes. In doing so, we have developed a standardised protocol for assessment of microsomal activity of CYP3A4, CYP1A2 and CYP2D6, but we were unable to optimise conditions for assessment of CYP2C9. This approach can be applied to investigate the impact of pregnancy complications on maternal and fetal hepatic drug metabolism.

A Novel Thermostable Cytochrome P450 from Sequence-Based Metagenomics of Binh Chau Hot Spring as a Promising Catalyst for Testosterone Conversion

Biotechnological applications of cytochromes P450 show difficulties, such as low activity, thermal and/or solvent instability, narrow substrate specificity and redox partner dependence. In an attempt to overcome these limitations, an exploitation of novel thermophilic P450 enzymes from nature via uncultured approaches is desirable due to their great advantages that can resolve nearly all mentioned impediments. From the metagenomics library of the Binh Chau hot spring, an open reading frame (ORF) encoding a thermostable cytochrome P450—designated as P450-T3—which shared 66.6% amino acid sequence identity with CYP109C2 of Sorangium cellulosum So ce56 was selected for further identification and characterization. The ORF was synthesized artificially and heterologously expressed in Escherichia coli C43(DE3) using the pET17b system. The purified enzyme had a molecular weight of approximately 43 kDa. The melting temperature of the purified enzyme was 76.2 °C and its apparent half-life at 60 °C was 38.7 min. Redox partner screening revealed that P450-T3 was reduced well by the mammalian AdR-Adx4-108 and the yeast Arh1-Etp1 redox partners. Lauric acid, palmitic acid, embelin, retinoic acid (all-trans) and retinoic acid (13-cis) demonstrated binding to P450-T3. Interestingly, P450-T3 also bound and converted testosterone. Overall, P450-T3 might become a good candidate for biocatalytic applications on a larger scale.

Structural basis for plant lutein biosynthesis from α-carotene

Two cytochrome P450 enzymes, CYP97A3 and CYP97C1, catalyze hydroxylations of the β- and ε-rings of α-carotene to produce lutein. Chirality is introduced at the C-3 atom of both rings, and the reactions are both pro-3R-stereospecific. We determined the crystal structures of CYP97A3 in substrate-free and complex forms with a nonnatural substrate and the structure of CYP97C1 in a detergent-bound form. The structures of CYP97A3 in different states show the substrate channel and the structure of CYP97C1 bound with octylthioglucoside confirms the binding site for the carotenoid substrate. Biochemical assays confirm that the ferredoxin-NADP⁺ reductase (FNR)-ferredoxin pair is used as the redox partner. Details of the pro-3R stereospecificity are revealed in the retinal-bound CYP97A3 structure. Further analysis indicates that the CYP97B clan bears similarity to the β-ring-specific CYP97A clan. Overall, our research describes the molecular basis for the last steps of lutein biosynthesis.