Specific residues and conformational plasticity define the substrate specificity of short-chain dehydrogenases/reductases.

Short-chain dehydrogenases/reductases (SDRs) are one of the most prevalent enzyme families distributed among the sequenced microorganisms. Despite the presence of a conserved catalytic tetrad and high structural similarity, these enzymes exhibit different substrate specificities. The insufficient knowledge regarding the amino acids underlying substrate specificity hinders the understanding of the SDRs' roles in diverse and significant biological processes. Here, we performed bioinformatic analysis, molecular modeling, and mutagenesis studies to identify the key residues that regulate the substrate specificities of two homologous microbial SDRs (i.e., DesE and KduD). Further, we investigated the impact of altering the physicochemical properties of these amino acids on enzyme activity. Interestingly, molecular dynamics simulations also suggest a critical role of enzyme conformational flexibility in substrate recognition and catalysis. Overall, our findings improve the understanding of microbial SDR substrate specificity and shed light on future rational design of more efficient and effective biocatalysts.

Exploring the catalytic diversity of two short-chain dehydrogenases/reductases from Stachybotrys chartarum.

Short-chain dehydrogenase/reductases (SDRs) belong to the NAD(P)(H)-dependent oxidoreductase superfamily, which have various functions of catalyzing oxidation/reduction reactions and have been generally used as powerful biocatalysts in the production of pharmaceuticals. In this study, ScSDR1 and ScSDR2, two new SDRs have been identified and characterized from Stachybotrys chartarum 3.5365. Substrate scope investigation revealed that both of the enzymes possessed the ability to oxidize ß-OH to ketone specifically, and exhibited substrate promiscuity and high stereo-selectivity for efficiently catalyzing the structurally different prochiral ketones to chiral alcohols. These findings not only suggest that ScSDR1 and ScSDR2 might be potent synthetic tools in drug research and development, but also provide good examples for further engineered enzymes with higher efficiency and stereo-selectivity.

Quantitative Analysis of mRNA and Protein Expression Levels of Aldo-Keto Reductase and Short-Chain Dehydrogenase/Reductase Isoforms in the Human Intestine.

Enzymes catalyzing the reduction reaction of xenobiotics are mainly members of the aldo-keto reductase (AKR) and short-chain dehydrogenase/reductase (SDR) superfamilies. The intestine, together with the liver, is responsible for first-pass effects and is an organ that determines the bioavailability of orally administered drugs. In this study, we evaluated the mRNA and protein expression levels of 12 AKR isoforms (AKR1A1, AKR1B1, AKR1B10, AKR1B15, AKR1C1, AKR1C2, AKR1C3, AKR1C4, AKR1D1, AKR1E2, AKR7A2, and AKR7A3) and 7 SDR isoforms (CBR1, CBR3, CBR4, DCXR, DHRS4, HSD11B1, and HSD17B12) in each region of the human intestine using next-generation sequencing and data-independent acquisition proteomics. At both the mRNA and protein levels, most AKR isoforms were highly expressed in the upper regions of the intestine, namely the duodenum and jejunum, and then declined toward the rectum. Among the members in the SDR superfamily, CBR1 and DHRS4 were highly expressed in the upper regions, whereas the expression levels of the other isoforms were almost uniform in all regions. Significant positive correlations between mRNA and protein levels were observed in AKR1A1, AKR1B1, AKR1B10, AKR1C3, AKR7A2, AKR7A3, CBR1, and CBR3. The mRNA level of AKR1B10 was highest, followed by AKR7A3 and CBR1, each accounting for more than 10% of the sum of all AKR and SDR levels in the small intestine. This expression profile in the human intestine was greatly different from that in the human liver, where AKR1C isoforms are predominantly expressed. SIGNIFICANCE STATEMENT: In this study comprehensively determined the mRNA and protein expression profiles of aldo-keto reductase (AKR) and short-chain dehydrogenase/reductase isoforms involved in xenobiotic metabolism in the human intestine and found that most of them are highly expressed in the upper region, where AKR1B10, AKR7A3, and CBR1 are predominantly expressed. Since the intestine is significantly involved in the metabolism of orally administered drugs, the information provided here is valuable for pharmacokinetic studies in drug development.

Genetic architecture of asthma in African American patients.

BACKGROUND: Asthma is a chronic inflammatory disorder with a strong genetic inheritance. Although more than 100 loci were reported through the genome-wide association study of European populations, the genetic underpinning of asthma in African American individuals remains largely elusive. OBJECTIVE: We aimed to identify genetic loci associated with asthma in African American individuals. METHODS: Three cohorts were genotyped at the Children's Hospital of Philadelphia by using the Illumina single-nucleotide polymorphism array platform. Genotype imputation was performed by using the Trans-Omics for Precision Medicine (TOPMed) reference panel, which includes whole genome sequencing data from more than 100,000 individuals. A meta-analysis of 3 Children's Hospital of Philadelphia cohorts and 10 Consortium on Asthma among African Ancestry Populations in the Americas cohorts, totaling 19,628 subjects, was conducted to identify genetic loci associated with asthma in African American individuals. RESULTS: Our study identified 12 loci surpassing the classical genome-wide significance threshold (5 × 10-8). Of those loci, 8 reached the stricter significance threshold (3 × 10-8). The 9p24.1 locus (rs10975467 [P = 1.63 × 10-8]) has previously been associated with asthma in European individuals. Six loci are associated with enhancer activities, 2 loci are in DNase I-hypersensitive regions, and all of them are associated with regulatory motifs. Moreover, the locus 11q13.4 (rs7480008) is an expression quantitative trait locus of XRRA1 in lung (P = 9.4 × 10-10), and the locus 13q14.3 (rs1543525) is a splicing quantitative trait locus of DHRS12 in lung (P = 1.1 × 10-13). CONCLUSIONS: Our findings provide candidate genetic loci for therapeutic target identification and prioritization for African populations.

Quantitative Evaluation of the Contribution of Each Aldo-Keto Reductase and Short-Chain Dehydrogenase/Reductase Isoform to Reduction Reactions of Compounds Containing a Ketone Group in the Human Liver.

Enzymes of the aldo-keto reductase (AKR) and short-chain dehydrogenase/reductase superfamilies are involved in the reduction of compounds containing a ketone group. In most cases, multiple isoforms appear to be involved in the reduction of a compound, and the enzyme(s) that are responsible for the reaction in the human liver have not been elucidated. The purpose of this study was to quantitatively evaluate the contribution of each isoform to reduction reactions in the human liver. Recombinant cytosolic isoforms were constructed, i.e., AKR1C1, AKR1C2, AKR1C3, AKR1C4, and carbonyl reductase 1 (CBR1), and a microsomal isoform, 11ß-hydroxysteroid dehydrogenase type 1 (HSD11B1), and their contributions to the reduction of 10 compounds were examined by extrapolating the relative expression of each reductase protein in human liver preparations to recombinant systems quantified by liquid chromatography-mass spectrometry. The reductase activities for acetohexamide, doxorubicin, haloperidol, loxoprofen, naloxone, oxcarbazepine, and pentoxifylline were predominantly catalyzed by cytosolic isoforms, and the sum of the contributions of individual cytosolic reductases was almost 100%. Interestingly, AKR1C3 showed the highest contribution to acetohexamide and loxoprofen reduction, although previous studies have revealed that CBR1 mainly metabolizes them. The reductase activities of bupropion, ketoprofen, and tolperisone were catalyzed by microsomal isoform(s), and the contributions of HSD11B1 were calculated to be 41%, 32%, and 104%, respectively. To our knowledge, this is the first study to quantitatively evaluate the contribution of each reductase to the reduction of drugs in the human liver. SIGNIFICANCE STATEMENT: To our knowledge, this is the first study to determine the contribution of aldo-keto reductase (AKR)-1C1, AKR1C2, AKR1C3, AKR1C4, carbonyl reductase 1, and 11ß-hydroxysteroid dehydrogenase type 1 to drug reductions in the human liver by utilizing the relative expression factor approach. This study found that AKR1C3 contributes to the reduction of compounds at higher-than-expected rates.

Functional Characterization and Synthetic Application of Is2-SDR, a Novel Thermostable and Promiscuous Ketoreductase from a Hot Spring Metagenome.

In a metagenome mining-based search of novel thermostable hydroxysteroid dehydrogenases (HSDHs), enzymes that are able to selectively oxidize/reduce steroidal compounds, a novel short-chain dehydrogenase/reductase (SDR), named Is2-SDR, was recently discovered. This enzyme, found in an Icelandic hot spring metagenome, shared a high sequence similarity with HSDHs, but, unexpectedly, showed no activity in the oxidation of the tested steroid substrates, e.g., cholic acid. Despite that, Is2-SDR proved to be a very active and versatile ketoreductase, being able to regio- and stereoselectively reduce a diversified panel of carbonylic substrates, including bulky ketones, α- and ß-ketoesters, and α-diketones of pharmaceutical relevance. Further investigations showed that Is2-SDR was indeed active in the regio- and stereoselective reduction of oxidized steroid derivatives, and this outcome was rationalized by docking analysis in the active site model. Moreover, Is2-SDR showed remarkable thermostability, with an apparent melting temperature (TM) around 75 °C, as determined by circular dichroism analysis, and no significant decrease in catalytic activity, even after 5 h at 80 °C. A broad tolerance to both water-miscible and water-immiscible organic solvents was demonstrated as well, thus, confirming the potential of this new biocatalyst for its synthetic application.

Tunable Toxicity of Bufadienolides is Regulated through a Configuration Inversion Catalyzed by a Short-Chain Dehydrogenase/Reductase.

Bufadienolides are toxic components widely found in amphibious toads that exhibit a wide range of biological activities. Guided by UPLC-QTOF-MS analysis, several 3-epi-bufadienolides with unique structures were isolated from the bile of the Asiatic toad, Bufo gargarizans. However, the enzymatic machinery of this epimerization in toads and its significance in chemical ecology remains poorly understood. Herein, we firstly compared the toxicities of two typical bufadienolides, bufalin (featuring a 14ß-hydroxyl) and resibufogenin (containing a 14, 15-epoxy group), with their corresponding 3-epi isomers in a zebrafish model. The results of the toxicology assays showed that the ratio of maximum non-toxic concentrations of these two pairs of compounds are 256 and 96 times, respectively, thereby indicating that 3-hydroxyl epimerization leads to a significant decrease in toxicity. Aiming to investigate the biotransformation of 3-epi bufadienolides in toads, we applied liver lysate to transform bufalin and found that it could stereoselectively catalyze the conversion of bufalin into its 3α-hydroxyl epimer. Following this, we cloned and characterized a short-chain dehydrogenase/reductase, HSE-1, from the toad liver cDNA library and verified its 3(ß→α)-hydroxysteroid epimerization activity. To the best of our knowledge, this is the first hydroxyl epimerase identified from amphibians that regulates the toxicity of animal-derived natural products.

Structure-based rational design of a short-chain dehydrogenase/reductase for improving activity toward mycotoxin patulin.

Patulin is a fatal mycotoxin that is widely detected in drinking water and fruit-derived products contaminated by diverse filamentous fungi. CgSDR from Candida guilliermondii represents the first NADPH-dependent short-chain dehydrogenase/reductase that catalyzes the reduction of patulin to the nontoxic E-ascladiol. To elucidate the catalytic mechanism of CgSDR, we solved its crystal structure in complex with cofactor and substrate. Structural analyses indicate that patulin is situated in a hydrophobic pocket adjacent to the cofactor, with the hemiacetal ring orienting toward the nicotinamide moiety of NADPH. In addition, we conducted structure-guided engineering to modify substrate-binding residue V187 and obtained variant V187F, V187K and V187W, whose catalytic activity was elevated by 3.9-, 2.2- and 1.7-fold, respectively. The crystal structures of CgSDR variants suggest that introducing additional aromatic stacking or hydrogen-bonding interactions to bind the lactone ring of patulin might account for the observed enhanced activity. These results illustrate the catalytic mechanism of SDR-mediated patulin detoxification for the first time and provide the upgraded variants that exhibit tremendous potentials in industrial applications.

SDR enzymes oxidize specific lipidic alkynylcarbinols into cytotoxic protein-reactive species.

Hundreds of cytotoxic natural or synthetic lipidic compounds contain chiral alkynylcarbinol motifs, but the mechanism of action of those potential therapeutic agents remains unknown. Using a genetic screen in haploid human cells, we discovered that the enantiospecific cytotoxicity of numerous terminal alkynylcarbinols, including the highly cytotoxic dialkynylcarbinols, involves a bioactivation by HSD17B11, a short-chain dehydrogenase/reductase (SDR) known to oxidize the C-17 carbinol center of androstan-3-alpha,17-beta-diol to the corresponding ketone. A similar oxidation of dialkynylcarbinols generates dialkynylketones, that we characterize as highly protein-reactive electrophiles. We established that, once bioactivated in cells, the dialkynylcarbinols covalently modify several proteins involved in protein-quality control mechanisms, resulting in their lipoxidation on cysteines and lysines through Michael addition. For some proteins, this triggers their association to cellular membranes and results in endoplasmic reticulum stress, unfolded protein response activation, ubiquitin-proteasome system inhibition and cell death by apoptosis. Finally, as a proof-of-concept, we show that generic lipidic alkynylcarbinols can be devised to be bioactivated by other SDRs, including human RDH11 and HPGD/15-PGDH. Given that the SDR superfamily is one of the largest and most ubiquitous, this unique cytotoxic mechanism-of-action could be widely exploited to treat diseases, in particular cancer, through the design of tailored prodrugs.

Nitric oxide enhances resistance of Pleurotus eryngii to cadmium stress by alleviating oxidative damage and regulating of short-chain dehydrogenase/reductase family.

The function and mechanism of nitric oxide (NO) in regulating Pleurotus eryngii biological response to cadmium (Cd) stress was evaluated by using anti-oxidation and short-chain dehydrogenase/reductase (SDR) family analysis. The fresh biomass of P. eryngii mycelia sharply decreased after treatment with 50 µM Cd; the lipid peroxidation and H2O2 accumulation in P. eryngii were found responsible for it. Proper exogenous supply of NO (150 µM SNP) alleviated the oxidative damage induced by Cd stress in P. eryngii, which reduced the accumulation of thiobarbituric acid reactive substances (TBARS) and H2O2. The activities of antioxidant enzymes (superoxide dismutase, peroxidase) were significantly increased to deal with Cd stress when treated with SNP (150 µM), and the content of proline was also closely related to NO-mediated reduction of Cd toxicity. Moreover, SDR family members were widely involved in the response to Cd stress, especially PleSCH70 gene was observed for the first time in participating in NO-mediated enhancement of Cd tolerance in P. eryngii. Taken together, this study provides new insights in understanding the tolerance mechanisms of P. eryngii to heavy metal and lays a foundation for molecular breeding of P. eryngii to improve its tolerance to environmental stress.

Crystal structure of L-arabinose 1-dehydrogenase as a short-chain reductase/dehydrogenase protein.

l-Arabinose 1-dehydrogenase (AraDH) catalyzes the NAD(P)+-dependent oxidation of l-arabinose to L-arabinono-1,4-lactone in the non-phosphorylative l-arabinose pathway, and is classified into glucose-fructose oxidoreductase and short-chain dehydrogenase/reductase (SDR). We herein report the crystal structure of a SDR-type AraDH (from Herbaspirillum huttiense) for the first time. The interactions between Asp49 and the 2'- and 3'-hydroxyl groups of NAD+ were consistent with strict specificity for NAD+. In a binding model for the substrate, Ser155 and Tyr168, highly conserved in the SDR superfamily, interacted with the C1 and/or C2 hydroxyl(s) of l-arabinose, whereas interactions between Asp107, Arg109, and Gln206 and the C2 and/or C3 hydroxyl(s) were unique to AraDH. Trp200 significantly contributed to the selectivities of the C4 hydroxyl and C6 methyl of substrates.

The Activation Parameters of a Cold-Adapted Short Chain Dehydrogenase Are Insensitive to Enzyme Oligomerization.

The structural principles of enzyme cold adaptation are of fundamental interest both for understanding protein evolution and for biotechnological applications. It has become clear in recent years that structural flexibility plays a major role in tuning enzyme activity at low temperatures, which is reflected by characteristic changes in the thermodynamic activation parameters for psychrophilic enzymes, compared to those of mesophilic and thermophilic ones. Hence, increased flexibility of the enzyme surface has been shown to lead to a lower enthalpy and a more negative entropy of activation, which leads to higher activity in the cold. This immediately raises the question of how enzyme oligomerization affects the temperature dependence of catalysis. Here, we address this issue by computer simulations of the catalytic reaction of a cold-adapted bacterial short chain dehydrogenase in different oligomeric states. Reaction free energy profiles are calculated at different temperatures for the tetrameric, dimeric, and monomeric states of the enzyme, and activation parameters are obtained from the corresponding computational Arrhenius plots. The results show that the activation free energy, enthalpy, and entropy are remarkably insensitive to the oligomeric state, leading to the conclusion that assembly of the subunit interfaces does not compromise cold adaptation, even though the mobilities of interfacial residues are indeed affected.

Crystal structure of a short-chain dehydrogenase/reductase from Burkholderia phymatum in complex with NAD.

Burkholderia phymatum is an important symbiotic nitrogen-fixing betaproteobacterium. B. phymatum is beneficial, unlike other Burkholderia species, which cause disease or are potential bioagents. Structural genomics studies at the SSGCID include characterization of the structures of short-chain dehydrogenases/reductases (SDRs) from multiple Burkholderia species. The crystal structure of a short-chain dehydrogenase from B. phymatum (BpSDR) was determined in space group C2221 at a resolution of 1.80 Å. BpSDR shares less than 38% sequence identity with any known structure. The monomer is a prototypical SDR with a well conserved cofactor-binding domain despite its low sequence identity. The substrate-binding cavity is unique and offers insights into possible functions and likely inhibitors of the enzymatic functions of BpSDR.

Crystal structure of a putative short-chain dehydrogenase/reductase from Paraburkholderia xenovorans.

Paraburkholderia xenovorans degrades organic wastes, including polychlorinated biphenyls. The atomic structure of a putative dehydrogenase/reductase (SDR) from P. xenovorans (PxSDR) was determined in space group P21 at a resolution of 1.45 Å. PxSDR shares less than 37% sequence identity with any known structure and assembles as a prototypical SDR tetramer. As expected, there is some conformational flexibility and difference in the substrate-binding cavity, which explains the substrate specificity. Uniquely, the cofactor-binding cavity of PxSDR is not well conserved and differs from those of other SDRs. PxSDR has an additional seven amino acids that form an additional unique loop within the cofactor-binding cavity. Further studies are required to determine how these differences affect the enzymatic functions of the SDR.

Unexpected Role of a Short-Chain Dehydrogenase/Reductase Family Protein in Type II Polyketide Biosynthesis.

We investigated the biosynthetic pathway of type II polyketide murayaquinone. The murayaquinone biosynthetic cluster contains genes for three putative short-chain dehydrogenase/reductase family enzymes including MrqF and MrqH with known functions and MrqM with unclear function. We report the functional characterization of MrqM for its role in murayaquinone biosynthesis. Our gene deletion experiment and structural elucidation of the accumulated intermediates revealed that MrqM is related with the second polyketide ring cyclization, because the inactivation of mrqM resulted in the accumulation of an off-pathway intermediate SEK43 and disrupted the second and third ring cyclization. Site-directed mutagenesis studies showed that two conserved residues in MrqM and homologous proteins Y151 and K155 are essential for its activity. The previously proposed second/third ring cyclase, MrqD, might instead play another important role in the chain releasing step of the murayaquinone biosynthesis.

Dehydrogenase reductase 9 (SDR9C4) and related homologs recognize a broad spectrum of lipid mediator oxylipins as substrates.

Bioactive oxylipins play multiple roles during inflammation and in the immune response, with termination of their actions partly dependent on the activity of yet-to-be characterized dehydrogenases. Here, we report that human microsomal dehydrogenase reductase 9 (DHRS9, also known as SDR9C4 of the short-chain dehydrogenase/reductase (SDR) superfamily) exhibits a robust oxidative activity toward oxylipins with hydroxyl groups located at carbons C9 and C13 of octadecanoids, C12 and C15 carbons of eicosanoids, and C14 carbon of docosanoids. DHRS9/SDR9C4 is also active toward lipid inflammatory mediator dihydroxylated Leukotriene B4 and proresolving mediators such as tri-hydroxylated Resolvin D1 and Lipoxin A4, although notably, with lack of activity on the 15-hydroxyl of prostaglandins. We also found that the SDR enzymes phylogenetically related to DHRS9, i.e., human SDR9C8 (or retinol dehydrogenase 16), the rat SDR9C family member known as retinol dehydrogenase 7, and the mouse ortholog of human DHRS9 display similar activity toward oxylipin substrates. Mice deficient in DHRS9 protein are viable, fertile, and display no apparent phenotype under normal conditions. However, the oxidative activity of microsomal membranes from the skin, lung, and trachea of Dhrs9-/- mice toward 1 µM Leukotriene B4 is 1.7- to 6-fold lower than that of microsomes from wild-type littermates. In addition, the oxidative activity toward 1 µM Resolvin D1 is reduced by about 2.5-fold with DHRS9-null microsomes from the skin and trachea. These results strongly suggest that DHRS9 might play an important role in the metabolism of a wide range of bioactive oxylipins in vivo.

Inactive pseudoenzyme subunits in heterotetrameric BbsCD, a novel short-chain alcohol dehydrogenase involved in anaerobic toluene degradation.

Anaerobic toluene degradation proceeds by fumarate addition to produce (R)-benzylsuccinate as first intermediate, which is further degraded via ß-oxidation by five enzymes encoded in the conserved bbs operon. This study characterizes two enzymes of this pathway, (E)-benzylidenesuccinyl-CoA hydratase (BbsH), and (S,R)-2-(α-hydroxybenzyl)succinyl-CoA dehydrogenase (BbsCD) from Thauera aromatica. BbsH, a member of the enoyl-CoA hydratase family, converts (E)-benzylidenesuccinyl-CoA to 2-(α-hydroxybenzyl)succinyl-CoA and was subsequently used in a coupled enzyme assay with BbsCD, which belongs to the short-chain dehydrogenases/reductase (SDR) family. The BbsCD crystal structure shows a C2-symmetric heterotetramer consisting of BbsC2 and BbsD2 dimers. BbsD subunits are catalytically active and capable of binding NAD+ and substrate, whereas BbsC subunits represent built-in pseudoenzyme moieties lacking all motifs of the SDR family required for substrate binding or catalysis. Molecular modeling studies predict that the active site of BbsD is specific for conversion of the (S,R)-diastereomer of 2-(α-hydroxybenzyl)succinyl-CoA to (S)-2-benzoylsuccinyl-CoA by hydride transfer to the re-face of nicotinamide adenine dinucleotide (NAD)+ . Furthermore, BbsC subunits are not engaged in substrate binding and merely serve as scaffold for the BbsD dimer. BbsCD represents a novel clade of related enzymes within the SDR family, which adopt a heterotetrameric architecture and catalyze the ß-oxidation of aromatic succinate adducts.

SDR7-6, a short-chain alcohol dehydrogenase/reductase family protein, regulates light-dependent cell death and defence responses in rice.

Lesion mimic mutants resembling the hypersensitive response without pathogen attack are an ideal material to understand programmed cell death, the defence response, and the cross-talk between defence response and development in plants. In this study, mic, a lesion mimic mutant from cultivar Yunyin treated with ethyl methanesulphonate (EMS), was screened. By map-based cloning, a short-chain alcohol dehydrogenase/reductase with an atypical active site HxxxK was isolated and designated as SDR7-6. It functions as a homomultimer in rice and is localized at the endoplasmic reticulum. The lesion mimic phenotype of the mutant is light-dependent. The mutant displayed an increased resistance response to bacterial blight, but reduced resistance to rice blast disease. The mutant and knockout lines showed increased reactive oxygen species, jasmonic acid content, antioxidant enzyme activity, and expression of pathogenicity-related genes, while chlorophyll content was significantly reduced. The knockout lines showed significant reduction in grain size, seed setting rate, 1000-grain weight, grain weight per plant, panicle length, and plant height. SDR7-6 is a new lesion mimic gene that encodes a short-chain alcohol dehydrogenase with atypical catalytic site. Disruption of SDR7-6 led to cell death and had adverse effects on multiple agricultural characters. SDR7-6 may act at the interface of the two defence pathways of bacterial blight and rice blast disease in rice.

The biosynthesis of thymol, carvacrol, and thymohydroquinone in Lamiaceae proceeds via cytochrome P450s and a short-chain dehydrogenase.

Thymol and carvacrol are phenolic monoterpenes found in thyme, oregano, and several other species of the Lamiaceae. Long valued for their smell and taste, these substances also have antibacterial and anti-spasmolytic properties. They are also suggested to be precursors of thymohydroquinone and thymoquinone, monoterpenes with anti-inflammatory, antioxidant, and antitumor activities. Thymol and carvacrol biosynthesis has been proposed to proceed by the cyclization of geranyl diphosphate to γ-terpinene, followed by a series of oxidations via p-cymene. Here, we show that γ-terpinene is oxidized by cytochrome P450 monooxygenases (P450s) of the CYP71D subfamily to produce unstable cyclohexadienol intermediates, which are then dehydrogenated by a short-chain dehydrogenase/reductase (SDR) to the corresponding ketones. The subsequent formation of the aromatic compounds occurs via keto-enol tautomerisms. Combining these enzymes with γ-terpinene in in vitro assays or in vivo in Nicotiana benthamiana yielded thymol and carvacrol as products. In the absence of the SDRs, only p-cymene was formed by rearrangement of the cyclohexadienol intermediates. The nature of these unstable intermediates was inferred from reactions with the γ-terpinene isomer limonene and by analogy to reactions catalyzed by related enzymes. We also identified and characterized two P450s of the CYP76S and CYP736A subfamilies that catalyze the hydroxylation of thymol and carvacrol to thymohydroquinone when heterologously expressed in yeast and N. benthamiana Our findings alter previous views of thymol and carvacrol formation, identify the enzymes involved in the biosynthesis of these phenolic monoterpenes and thymohydroquinone in the Lamiaceae, and provide targets for metabolic engineering of high-value terpenes in plants.

Genome-Wide Identification and Characterization of Short-Chain Dehydrogenase/Reductase (SDR) Gene Family in Medicago truncatula.

Short-chain dehydrogenase/reductase (SDR) belongs to the NAD(P)(H)-dependent oxidoreductase superfamily. Limited investigations reveal that SDRs participate in diverse metabolisms. A genome-wide identification of the SDR gene family in M. truncatula was conducted. A total of 213 MtSDR genes were identified, and they were distributed on all chromosomes unevenly. MtSDR proteins were categorized into seven subgroups based on phylogenetic analysis and three types including 'classic', 'extended', and 'atypical', depending on the cofactor-binding site and active site. Analysis of the data from M. truncatula Gene Expression Atlas (MtGEA) showed that above half of MtSDRs were expressed in at least one organ, and lots of MtSDRs had a preference in a tissue-specific expression. The cis-acting element responsive to plant hormones (salicylic acid, ABA, auxin, MeJA, and gibberellin) and stresses were found in the promoter of some MtSDRs. Many genes of MtSDR7C,MtSDR65C, MtSDR110C, MtSDR114C, and MtSDR108E families were responsive to drought, salt, and cold. The study provides useful information for further investigation on biological functions of MtSDRs, especially in abiotic stress adaptation, in the future.