Comparative structural modeling of six old yellow enzymes (OYEs) from the necrotrophic fungus Ascochyta rabiei: insight into novel OYE classes with differences in cofactor binding, organization of active site residues and stereopreferences

Unlocking the function promiscuity of old yellow enzyme to catalyze asymmetric Morita-Baylis-Hillman reaction

Exploring the promiscuity of native enzymes presents a promising strategy for expanding their synthetic applications, particularly for catalyzing challenging reactions in non-native contexts. In this study, we explore the promiscuous potential of old yellow enzymes (OYEs) to facilitate the Morita-Baylis-Hillman reaction (MBH reaction), leveraging substrate similarities between MBH reaction and reduction reaction. Using mass spectrometry and spectroscopic techniques, we confirm promiscuity of GkOYE in both MBH and reduction reactions. By blocking H- and H+ transfer pathways, we engineer GkOYE.8, which loses its reduction ability but enhances its MBH activity. The structural basis of MBH reaction catalyzed by GkOYE.8 is obtained through mutation studies and kinetic simulations. Furthermore, enantiocomplementary mutants GkOYE.11 and GkOYE.13 are obtained by directed evolution, exhibiting the ability to accept various aromatic aldehydes and alkenes as substrates. This study demonstrates the potential of leveraging substrate similarities to unlock enzyme functionalities, enabling the catalysis of new-to-nature reactions.

Loop 6 and the β-hairpin flap are structural hotspots that determine cofactor specificity in the FMN-dependent family of ene-reductases

Flavin mononucleotide (FMN)-dependent ene-reductases constitute a large family of oxidoreductases that catalyze the enantiospecific reduction of carbon-carbon double bonds. The reducing equivalents required for substrate reduction are obtained from reduced nicotinamide by hydride transfer. Most ene-reductases significantly prefer, or exclusively accept, either NADPH or NADH. Despite their usefulness in biocatalytic applications, the structural determinants for cofactor preference remain elusive. We employed the NADPH-preferring 12-oxophytodienoic acid reductase 3 from Solanum lycopersicum (SlOPR3) as a model enzyme of the ene-reductase family and applied computational and structural methods to investigate the binding specificity of the reducing coenzymes. Initial docking results indicated that the arginine triad R283, R343, and R366 residing on and close to a critical loop at the active site (loop 6) are the main contributors to NADPH binding. In contrast, NADH binds unfavorably in the opposite direction toward the β-hairpin flap within a largely hydrophobic region. Notably, the crystal structures of SlOPR3 in complex with either NADPH4 or NADH4 corroborated these different binding modes. Molecular dynamics simulations confirmed NADH binding near the β-hairpin flap and provided structural explanations for the low binding affinity of NADH to SlOPR3. We postulate that cofactor specificity is determined by the arginine triad/loop 6 and the residue(s) controlling access to a hydrophobic cleft formed by the β-hairpin flap. Thus, NADPH preference depends on a properly positioned arginine triad, whereas granting access to the hydrophobic cleft at the β-hairpin flap favors NADH binding.

Regiodivergent Radical Termination for Intermolecular Biocatalytic C-C Bond Formation

Radical hydrofunctionalizations of electronically unbiased dienes are challenging to render regioselective, because the products are nearly identical in energy. Here, we report two engineered FMN-dependent "ene"-reductases (EREDs) that catalyze regiodivergent hydroalkylations of cyclic and linear dienes. While previous studies focused exclusively on the stereoselectivity of alkene hydroalkylation, this work highlights that EREDs can control the regioselectivity of hydrogen atom transfer, providing a method for selectively preparing constitutional isomers that would be challenging to prepare using traditional synthetic methods. Engineering the ERED from Gluconabacter sp. (GluER) furnished a variant that favors the γ,δ-unsaturated ketone, while an engineered variant from a commercial ERED panel favors the δ,ε-unsaturated ketone. The effect of beneficial mutations has been investigated using substrate docking studies and the mechanism probed by isotope labeling experiments. A variety of α-bromo ketones can be coupled with cyclic and linear dienes. These interesting building blocks can also be further modified to generate difficult-to-access heterocyclic compounds.

The Family Keeps on Growing: Four Novel Fungal OYEs Characterized

Aiming at expanding the portfolio of Old Yellow Enzymes (OYEs), which have been systematically studied to be employed in the chemical and pharmaceutical industries as useful biocatalysts, we decided to explore the immense reservoir of filamentous fungi. We drew from the genome of the two Ascomycetes Aspergillus niger and Botryotinia fuckeliana four new members of the OYE superfamily belonging to the classical and thermophilic-like subfamilies. The two BfOYEs show wider substrate spectra than the AnOYE homologues, which appear as more specialized biocatalysts. According to their mesophilic origins, the new enzymes neither show high thermostability nor extreme pH optimums. The crystal structures of BfOYE4 and AnOYE8 have been determined, revealing the conserved features of the thermophilic-like subclass as well as unique properties, such as a peculiar N-terminal loop involved in dimer surface interactions. For the classical representatives BfOYE1 and AnOYE2, model structures were built and analyzed, showing surprisingly wide open access to the active site cavities due to a shorter β6-loop and a disordered capping subdomain.

Identification of four ene reductases and their preliminary exploration in the asymmetric synthesis of (R)-dihydrocarvone and (R)-profen derivatives

The ene reductases (ERs) from the old yellow enzymes (OYEs) family have the ability to reduce activated alkenes to generate up to two stereocenters, therefore they have been received extensive attention as powerful biocatalysts. In this study, through gene mining, four ERs were identified from the genomes of Ensifer adhaerens, Pseudomonas fluorescens, and Pseudomonas veronil. The biocatalytic properties of these four ERs were identified, and their applications in the synthesis process of dihydrocarvone and profen derivatives were further evaluated. Among them, three ERs (EaER2, PvER1, and PvER2) belonging to the classic OYEs showed the best catalytic activity at 30 °C and pH 7.0 (100 mM potassium phosphate buffer) and the PfER2, which belongs to the thermophilic-like OYEs exhibited the best catalytic at 40 °C and pH 7.0 (100 mM potassium phosphate buffer). When exploring the influence of organic solvents on the catalytic efficiency, it was found that the four ERs were more sensitive to toluene and had tolerance to several other selected organic solvents. In addition, EaER2, PfER2, PvER1 and PvER2 showed excellent catalytic activity toward carvone, and the stereoselectivity of PvER2 toward carvone could reach up to 88.7 % de. EaER2 and PfER2 can catalyze the synthesis of a variety of profen derivatives with a stereoselectivity over 99 % ee. Moreover, through homology modeling and molecular docking, we preliminarily explained the mechanism of catalytic activity and stereoselectivity of the four ERs, which provided a solid base on the rational design of their stereo-preference in the future. The discovery of EaER2, PfER2, PvER1, and PvER2 provides four new enzyme sources for the study of the OYEs family and enriches the biocatalytic toolbox of ERs. Our exploration of the enzymatic properties of these four ERs will provide the sufficient data basis for future research and industrialization progress.

The Interaction of Porcine Dihydropyrimidine Dehydrogenase with the Chemotherapy Sensitizer: 5-Ethynyluracil

Dihydropyrimidine dehydrogenase (DPD) is a complex enzyme that reduces the 5,6-vinylic bond of pyrimidines, uracil, and thymine. 5-Fluorouracil (5FU) is also a substrate for DPD and a common chemotherapeutic agent used to treat numerous cancers. The reduction of 5FU to 5-fluoro-5,6-dihydrouracil negates its toxicity and efficacy. Patients with high DPD activity levels typically have poor outcomes when treated with 5FU. DPD is thus a central mitigating factor in the treatment of a variety of cancers. 5-Ethynyluracil (5EU) covalently inactivates DPD by cross-linking with the active-site general acid cysteine in the pyrimidine binding site. This reaction is dependent on the simultaneous binding of 5EU and nicotinamide adenine dinucleotide phosphate (NADPH). This ternary complex induces DPD to become activated by taking up two electrons from the NADPH. The covalent inactivation of DPD by 5EU occurs concomitantly with this reductive activation with a rate constant of ∼0.2 s^-1. This k_inact value is correlated with the rate of reduction of one of the two flavin cofactors and the localization of a mobile loop in the pyrimidine active site that places the cysteine that serves as the general acid in catalysis proximal to the 5EU ethynyl group. Efficient cross-linking is reliant on enzyme activation, but this process appears to also have a conformational aspect in that nonreductive NADPH analogues can also induce a partial inactivation. Cross-linking then renders DPD inactive by severing the proton-coupled electron transfer mechanism that transmits electrons 56 Å across the protein.

Quaternary Charge-Transfer Complex Enables Photoenzymatic Intermolecular Hydroalkylation of Olefins

Intermolecular C-C bond-forming reactions are underdeveloped transformations in the field of biocatalysis. Here we report a photoenzymatic intermolecular hydroalkylation of olefins catalyzed by flavin-dependent 'ene'-reductases. Radical initiation occurs via photoexcitation of a rare high-order enzyme-templated charge-transfer complex that forms between an alkene, α-chloroamide, and flavin hydroquinone. This unique mechanism ensures that radical formation only occurs when both substrates are present within the protein active site. This active site can control the radical terminating hydrogen atom transfer, enabling the synthesis of enantioenriched γ-stereogenic amides. This work highlights the potential for photoenzymatic catalysis to enable new biocatalytic transformations via previously unknown electron transfer mechanisms.

Old yellow enzymes: structures and structure-guided engineering for stereocomplementary bioreduction

Since the first discovery of old yellow enzyme 1 (OYE1) from Saccharomyces pastorianus in 1932, biocatalytic asymmetric reduction of activated alkenes by OYEs has become a valuable reaction in organic synthesis. To access stereocomplementary C=C-bond bioreduction, the mining of novel OYEs and especially the protein engineering of existing OYEs have been performed, which successfully achieved the stereocomplementary reduction in several cases and further raise the potential of applications. In this review, we analyzed the structures, active sites, and substrate recognition of OYEs, which are the bases for their substrate specificity and stereospecificity. Sequence similarity network of OYEs superfamily was also constructed to investigate the scope of characterized OYEs. The structure-guided engineering to switch the stereoselectivity of OYEs and thus access stereocomplementary bioreduction over the last decade (2009-2020) was then reviewed and discussed, which might give new insights into the mining and engineering of related biocatalysts. KEY POINTS: • The sequence similarity network of OYEs superfamily was constructed and annotated. • The structures and active sites of OYEs from different classes were compared. • "Left/right" binding mode was used to explain the stereopreferences of OYEs. • Structure-guided engineering of OYEs to switch their stereoselectivity was reviewed.

Transient-State Analysis of Human Isocitrate Dehydrogenase I: Accounting for the Interconversion of Active and Non-Active Conformational States

Human isocitrate dehydrogenase 1 (HsICDH1) is a cytoplasmic homodimeric Mg(II)-dependent enzyme that converts d-isocitrate (D-ICT) and NADP⁺ to α-ketoglutarate (AKG), CO₂, and NADPH. The active sites are formed at the subunit interface and incorporate residues from both protomers. The turnover number titrates hyperbolically from 17.5 s^-1 to a minimum of 7 s^-1 with an increasing enzyme concentration. As isolated, the enzyme adopts an inactive open conformation and binds NADPH tightly. The open conformation displaces three of the eight residues that bind D-ICT and Mg(II). Enzyme activation occurs with the addition of Mg(II) or D-ICT with a rate constant of 0.12 s^-1. The addition of both Mg(II) and D-ICT activates the enzyme with a rate constant of 0.6 s^-1 and displaces half of the bound NADPH. This indicates that HsICDH1 may have a half-site mechanism in which the active sites alternate in catalysis. The X-ray crystal structure of the half-site activated complex reveals asymmetry in the homodimer with a single NADPH bound. The structure also indicates a pseudotetramer interface that impedes the egress of NADPH consistent with the suppression of the turnover number at high enzyme concentrations. When the half-site activated form of the enzyme is reacted with NADP⁺, NADPH forms with a rate constant of 204 s^-1 followed by a shift in the NADPH absorption spectrum with a rate constant of 28 s^-1. These data indicate the accumulation of two intermediate states. Once D-ICT is exhausted, HsICDH1 relaxes to the inactive open state with a rate constant of ∼3 s^-1.

Daidzein reductase of Eggerthella sp. YY7918, its octameric subunit structure containing FMN/FAD/4Fe-4S, and its enantioselective production of R-dihydroisoflavones

S-Equol is a metabolite of daidzein, a type of soy isoflavone, and three reductases are involved in the conversion of daidzein by specific intestinal bacteria. S-Equol is thought to prevent hormone-dependent diseases. We previously identified the equol producing gene cluster (eqlABC) of Eggerthella sp. YY7918. Daidzein reductase (DZNR), encoded by eqlA, catalyzes the reduction of daidzein to dihydrodaidzein (the first step of equol synthesis), which was confirmed using a recombinant enzyme produced in Escherichia coli. Here, we purified recombinant DZNR to homogeneity and analyzed its enzymological properties. DZNR contained FMN, FAD, and one 4Fe-4S cluster per 70-kDa subunit as enzymatic cofactors. DZNR reduced the CC bond between the C-2 and C-3 positions of daidzein, genistein, glycitein, and formononetin in the presence of NADPH. R-Dihydrodaidzein and R-dihydrogenistein were highly stereo-selectively produced from daidzein and genistein. The K_m and k_cat for daidzein were 11.9 μM and 6.7 s^-1, and these values for genistein were 74.1 μM and 28.3 s^-1, respectively. This enzyme showed similar kinetic parameters and wide substrate specificity for isoflavone molecules. Thus, this enzyme appears to be an isoflavone reductase. Gel filtration chromatography and chemical cross-linking analysis of the active form of DZNR suggested that the enzyme consists of an octameric subunit structure. We confirmed this by small-angle X-ray scattering and transmission electron microscopy at a magnification of ×200,000. DZNR formed a globular four-petal cloverleaf structure with a central vertical hole. The maximum particle size was 173 Å.

Old Yellow Enzyme homologues in Mucor circinelloides: expression profile and biotransformation

The reduction of C=C double bond, a key reaction in organic synthesis, is mostly achieved by traditional chemical methods. Therefore, the search for enzymes capable of performing this reaction is rapidly increasing. Old Yellow Enzymes (OYEs) are flavin-dependent oxidoreductases, initially isolated from Saccharomyces pastorianus. In this study, the presence and activation of putative OYE enzymes was investigated in the filamentous fungus Mucor circinelloides, which was previously found to mediate C=C reduction. Following an in silico approach, using S. pastorianus OYE1 amminoacidic sequence as template, ten putative genes were identified in the genome of M. circinelloides. A phylogenetic analysis revealed a high homology of McOYE1-9 with OYE1-like proteins while McOYE10 showed similarity with thermophilic-like OYEs. The activation of mcoyes was evaluated during the transformation of three different model substrates. Cyclohexenone, α-methylcinnamaldehyde and methyl cinnamate were completely reduced in few hours and the induction of gene expression, assessed by qRT-PCR, was generally fast, suggesting a substrate-dependent activation. Eight genes were activated in the tested conditions suggesting that they may encode for active OYEs. Their expression over time correlated with C=C double bond reduction.

Transcription Factor Repertoire of Necrotrophic Fungal Phytopathogen Ascochyta rabiei: Predominance of MYB Transcription Factors As Potential Regulators of Secretome

Transcription factors (TFs) are the key players in gene expression and their study is highly significant for shedding light on the molecular mechanisms and evolutionary history of organisms. During host-pathogen interaction, extensive reprogramming of gene expression facilitated by TFs is likely to occur in both host and pathogen. To date, the knowledge about TF repertoire in filamentous fungi is in infancy. The necrotrophic fungus Ascochyta rabiei, that causes destructive Ascochyta blight (AB) disease of chickpea (Cicer arietinum), demands more comprehensive study for better understanding of Ascochyta-legume pathosystem. In the present study, we performed the genome-wide identification and analysis of TFs in A. rabiei. Taking advantage of A. rabiei genome sequence, we used a bioinformatic approach to predict the TF repertoire of A. rabiei. For identification and classification of A. rabiei TFs, we designed a comprehensive pipeline using a combination of BLAST and InterProScan software. A total of 381 A. rabiei TFs were predicted and divided into 32 fungal specific families of TFs. The gene structure, domain organization and phylogenetic analysis of abundant families of A. rabiei TFs were also carried out. Comparative study of A. rabiei TFs with that of other necrotrophic, biotrophic, hemibiotrophic, symbiotic, and saprotrophic fungi was performed. It suggested presence of both conserved as well as unique features among them. Moreover, cis-acting elements on promoter sequences of earlier predicted A. rabiei secretome were also identified. With the help of published A. rabiei transcriptome data, the differential expression of TF and secretory protein coding genes was analyzed. Furthermore, comprehensive expression analysis of few selected A. rabiei TFs using quantitative real-time polymerase chain reaction revealed variety of expression patterns during host colonization. These genes were expressed in at least one of the time points tested post infection. Overall, this study illustrates the first genome-wide identification and analysis of TF repertoire of A. rabiei. This work would provide the basis for further studies to dissect role of TFs in the molecular mechanisms during A. rabiei-chickpea interactions.

Structural insights into the ene-reductase synthesis of profens

Reduction of double bonds of α,β-unsaturated carboxylic acids and esters by ene-reductases remains challenging and it typically requires activation by a second electron-withdrawing moiety, such as a halide or second carboxylate group. We showed that profen precursors, 2-arylpropenoic acids and their esters, were efficiently reduced by Old Yellow Enzymes (OYEs). The XenA and GYE enzymes showed activity towards acids, while a wider range of enzymes were active towards the equivalent methyl esters. Comparative co-crystal structural analysis of profen-bound OYEs highlighted key interactions important in determining substrate binding in a catalytically active conformation. The general utility of ene reductases for the synthesis of (R)-profens was established and this work will now drive future mutagenesis studies to screen for the production of pharmaceutically-active (S)-profens.

Draft genome sequencing and secretome analysis of fungal phytopathogen Ascochyta rabiei provides insight into the necrotrophic effector repertoire

Constant evolutionary pressure acting on pathogens refines their molecular strategies to attain successful pathogenesis. Recent studies have shown that pathogenicity mechanisms of necrotrophic fungi are far more intricate than earlier evaluated. However, only a few studies have explored necrotrophic fungal pathogens. Ascochyta rabiei is a necrotrophic fungus that causes devastating blight disease of chickpea (Cicer arietinum). Here, we report a 34.6 megabase draft genome assembly of A. rabiei. The genome assembly covered more than 99% of the gene space and 4,259 simple sequence repeats were identified in the assembly. A total of 10,596 high confidence protein-coding genes were predicted which includes a large and diverse inventory of secretory proteins, transporters and primary and secondary metabolism enzymes reflecting the necrotrophic lifestyle of A. rabiei. A wide range of genes encoding carbohydrate-active enzymes capable for degradation of complex polysaccharides were also identified. Comprehensive analysis predicted a set of 758 secretory proteins including both classical and non-classical secreted proteins. Several of these predicted secretory proteins showed high cysteine content and numerous tandem repeats. Together, our analyses would broadly expand our knowledge and offer insights into the pathogenesis and necrotrophic lifestyle of fungal phytopathogens.

Genome Wide Identification of LIM Genes in Cicer arietinum and Response of Ca-2LIMs in Development, Hormone and Pathogenic Stress

The eukaryotic lineage-specific LIM protein (LIN11, ISL1, and MEC3) family play pivotal role in modulation of actin dynamics and transcriptional regulation. The systematic investigation of this family has not been carried in detail and rare in legumes. Current study involves the mining of Cicer arietinum genome for the genes coding for LIM domain proteins and displayed significant homology with LIM genes of other species. The analysis led to the identification of 15 members, which were positioned on chickpea chromosomes. The phylogenetic and motif analysis suggested their categorization into two sub-families i.e., Ca-2LIMs and Ca-DA1/DAR, which comprised of nine and six candidates, respectively. Further sub-categories of Ca-2LIMs were recognised as αLIM, βLIM, δLIM and γLIM. The LIM genes within their sub-families displayed conserved genomic and motif organization. The expression pattern of Ca-2LIMs across developmental and reproductive tissues demonstrated strong correlation with established consensus. The Ca-2LIM belongs to PLIM and GLIM (XLIM) was found highly expressed in floral tissue. Others showed ubiquitous expression pattern with their dominance in stem. Under hormonal and pathogenic conditions these LIMs were found to up-regulate during salicylic acid, abscisic acid and Ascochyta rabiei treatment or infection; and down-regulated in response to jasmonic acid treatment. The findings of this work, particularly in terms of modulation of LIM genes under biotic stress will open up the way to further explore and establish the role of chickpea LIMs in plant defense response.

Characterization and proteomic analysis of the Pseudomonas sp. HK-6 xenB knockout mutant under RDX (hexahydro-1,3,5-trinitro-1,3,5-triazine) stress

Pseudomonas sp. HK-6 is able to utilize RDX (hexahydro-1,3,5-trinitro-1,3,5-triazine) as its sole nitrogen source. The role of the xenB gene, encoding xenobiotic reductase B, was investigated using HK-6 xenB knockout mutants. The xenB mutant degraded RDX to a level that was 10-fold less than that obtained with the wild-type HK-6 strain. After 60 days of culture with 25 or 50 μM RDX, no residual RDX was detected in the supernatants of the wild-type aerobically grown cultures, whereas approximately 90 % of the RDX remained in the xenB mutant cultures. The xenB mutant bacteria exhibited a 10(2)-10(4)-fold decrease in survival rate compared to the wild-type. The expression of DnaK and GroEL proteins, two typical stress shock proteins (SSPs), in the xenB mutant increased after immediate exposure to RDX, yet dramatically decreased after 4 h of exposure. In addition, DnaK and GroEL were more highly expressed in the cultures with 25 μM RDX in the medium but showed low expression in the cultures with 50 or 75 μM RDX. The expression levels of the dnaK and groEL genes measured by RT-qPCR were also much lower in the xenB genetic background. Analyses of the proteomes of the HK-6 and xenB mutant cells grown under conditions of RDX stress showed increased induction of several proteins, such as Alg8, alginate biosynthesis sensor histidine kinase, and OprH in the xenB mutants when compared to wild-type. However, many proteins, including two SSPs (DnaK and GroEL) and proteins involved in metabolism, exhibited lower expression levels in the xenB mutant than in the wild-type HK-6 strain. The xenB knockout mutation leads to reduced RDX degradation ability, which renders the mutant more sensitive to RDX stress and results in a lower survival rate and an altered proteomic profile under RDX stress.