Molecular cloning and heterologous expression of progesterone 5beta-reductase from Digitalis lanata Ehrh

Addressing the Evolution of Cardenolide Formation in Iridoid-Synthesizing Plants: Site-Directed Mutagenesis of PRISEs (Progesterone-5β-Reductase/Iridoid Synthase-like Enzymes) of Plantago Species

Enzymes capable of processing a variety of compounds enable plants to adapt to diverse environmental conditions. PRISEs (progesterone-5β-reductase/iridoid synthase-like enzymes), examples of such substrate-promiscuous enzymes, are involved in iridoid and cardenolide pathways and demonstrate notable substrate promiscuity by reducing the activated C=C double bonds of plant-borne and exogenous 1,4-enones. In this study, we identified PRISE genes in Plantago media (PmdP5βR1) and Plantago lanceolata (PlP5βR1), and the corresponding enzymes were determined to share a sequence identity of 95%. Despite the high sequence identity, recombinant expressed PmdP5βR1 was 70 times more efficient than PlP5βR1 for converting progesterone. In order to investigate the underlying reasons for this significant discrepancy, we focused on specific residues located near the substrate-binding pocket and adjacent to the conserved phenylalanine "clamp". This clamp describes two phenylalanines influencing substrate preferences by facilitating the binding of smaller substrates, such as 2-cyclohexen-1-one, while hindering larger ones, such as progesterone. Using structural analysis based on templates PDB ID: 5MLH and 6GSD from PRISE of Plantago major, along with in silico docking, we identified positions 156 and 346 as hot spots. In PlP5βR1 amino acid residues, A156 and F346 seem to be responsible for the diminished ability to reduce progesterone. Moreover, the double mutant PlP5βR_F156L_A346L, which contains the corresponding amino acids from PmdP5βR1, showed a 15-fold increase in progesterone 5β-reduction. Notably, this modification did not significantly alter the enzyme's ability to convert other substrates, such as 8-oxogeranial, 2-cyclohexen-1-one, and methyl vinyl ketone. Hence, a rational enzyme design by reducing the number of hotspots selectively, specifically improved the substrate preference of PlP5βR1 for progesterone.

Progesterone Metabolism in Digitalis and Other Plants-60 Years of Research and Recent Results

5β-Cardenolides are pharmaceutically important metabolites from the specialized metabolism of Digitalis lanata. They were used over decades to treat cardiac insufficiency and supraventricular tachycardia. Since the 1960s, plant scientists have known that progesterone is an essential precursor of cardenolide formation. Therefore, biosynthesis of plant progesterone was mainly analyzed in species of the cardenolide-containing genus Digitalis during the following decades. Today, Digitalis enzymes catalyzing the main steps of progesterone biosynthesis are known. Most of them are found in a broad range of organisms. This review will summarize the findings of 60 years of research on plant progesterone metabolism with particular focus on the recent results in Digitalis lanata and other plants.

Four enzymes control natural variation in the steroid core of Erysimum cardenolides

Plants commonly produce families of structurally related metabolites with similar defensive functions. This apparent redundancy raises the question of underlying molecular mechanisms and adaptive benefits of such chemical variation. Cardenolides, a class defensive compounds found in the wallflower genus Erysimum (L., Brassicaceae) and scattered across other plant families, show substantial structural variation, with glycosylation and hydroxylation being common modifications of a steroid core, which itself may vary in terms of stereochemistry and saturation. Through a combination of chemical mutagenesis and analysis of gene coexpression networks, we identified four enzymes involved in cardenolide biosynthesis in Erysimum that work together to determine stereochemistry at carbon 5 of the steroid core: Ec3βHSD, a 3β-hydroxysteroid dehydrogenase, Ec3KSI, a ketosteroid isomerase, EcP5βR2, a progesterone 5β-reductase, and EcDET2, a steroid 5α-reductase. We biochemically characterized the activity of these enzymes in vitro and generated CRISPR/Cas9 knockout lines to confirm activity in vivo. Cardenolide biosynthesis was not eliminated in any of the knockouts. Instead, mutant plants accumulated cardenolides with altered saturation and stereochemistry of the steroid core. Furthermore, we found variation in carbon 5 configuration among the cardenolides of 44 species of Erysimum, where the occurrence of some 5β-cardenolides is associated with the expression and sequence of P5βR2. This may have allowed Erysimum species to fine-tune their defensive profiles to target specific herbivore populations over the course of evolution.

SIGNIFICANCE STATEMENT

Plants use an array of toxic compounds to defend themselves from attack against insects and other herbivores. One mechanism through which plants may evolve more toxic compounds is through modifications to the structure of compounds they already produce. In this study, we show how plants in the wallflower genus Erysimum use four enzymes to fine-tune the structure of toxic metabolites called cardenolides. Natural variation in the sequence and expression of a single enzyme called progesterone 5β-reductase 2 partly explains the variation in cardenolides observed across the Erysimum genus. These alterations to cardenolide structure over the course of evolution suggests that there may be context-dependent benefits to Erysimum to invest in one cardenolide variant over another.

Citronellol biosynthesis in pelargonium is a multistep pathway involving progesterone 5β-reductase and/or iridoid synthase-like enzymes

Citronellol is a pleasant-smelling compound produced in rose (Rosa spp.) flowers and in the leaves of many aromatic plants, including pelargoniums (Pelargonium spp.). Although geraniol production has been well studied in several plants, citronellol biosynthesis has been documented only in crab-lipped spider orchid (Caladenia plicata) and its mechanism remains open to question in other species. We therefore profiled 10 pelargonium accessions using RNA sequencing and gas chromatography-MS analysis. Three enzymes from the progesterone 5β-reductase and/or iridoid synthase-like enzymes (PRISE) family were characterized in vitroand subsequently identified as citral reductases (named PhCIRs). Transgenic RNAi lines supported a role for PhCIRs in the biosynthesis of citronellol as well as in the production of mint-scented terpenes. Despite their high amino acid sequence identity, the 3 enzymes showed contrasting stereoselectivity, either producing mainly (S)-citronellal or a racemate of both (R)- and (S)-citronellal. Using site-directed mutagenesis, we identified a single amino acid substitution as being primarily responsible for the enzyme's enantioselectivity. Phylogenetic analysis of pelargonium PRISEs revealed 3 clades and 7 groups of orthologs. PRISEs from different groups exhibited differential affinities toward substrates (citral and progesterone) and cofactors (NADH/NADPH), but most were able to reduce both substrates, prompting hypotheses regarding the evolutionary history of PhCIRs. Our results demonstrate that pelargoniums evolved citronellol biosynthesis independently through a 3-step pathway involving PRISE homologs and both citral and citronellal as intermediates. In addition, these enzymes control the enantiomeric ratio of citronellol thanks to small alterations of the catalytic site.

A cytochrome P450 CYP87A4 imparts sterol side-chain cleavage in digoxin biosynthesis

Digoxin extracted from the foxglove plant is a widely prescribed natural product for treating heart failure. It is listed as an essential medicine by the World Health Organization. However, how the foxglove plant synthesizes digoxin is mostly unknown, especially the cytochrome P450 sterol side chain cleaving enzyme (P450scc), which catalyzes the first and rate-limiting step. Here we identify the long-speculated foxglove P450scc through differential transcriptomic analysis. This enzyme converts cholesterol and campesterol to pregnenolone, suggesting that digoxin biosynthesis starts from both sterols, unlike previously reported. Phylogenetic analysis indicates that this enzyme arises from a duplicated cytochrome P450 CYP87A gene and is distinct from the well-characterized mammalian P450scc. Protein structural analysis reveals two amino acids in the active site critical for the foxglove P450scc's sterol cleavage ability. Identifying the foxglove P450scc is a crucial step toward completely elucidating digoxin biosynthesis and expanding the therapeutic applications of digoxin analogs in future work.

Cardenolide Increase in Foxglove after 2,1,3-Benzothiadiazole Treatment Reveals a Potential Link between Cardenolide and Phytosterol Biosynthesis

Cardenolides are steroidal metabolites in Digitalis lanata with potent cardioactive effects on animals. In plants, cardenolides are likely involved in various stress responses. However, the molecular mechanism of cardenolide increase during stresses is mostly unknown. Additionally, cardenolides are proposed to arise from cholesterol, but indirect results show that phytosterols may also be substrates for cardenolide biosynthesis. Here, we show that cardenolides increased after methyl jasmonate (MJ), sorbitol, potassium chloride (KCl) and salicylic acid analog [2,1,3-benzothiadiazole (BTH)] treatments. However, the expression of three known genes for cardenolide biosynthesis did not correlate well with these increases. Specifically, the expression of progesterone-5β-reductases (P5βR and P5βR2) did not correlate with the cardenolide increase. The expression of 3β-hydroxysteroid dehydrogenase (3βHSD) correlated with changes in cardenolide levels only during the BTH treatment. Mining the D. lanata transcriptome identified genes involved in cholesterol and phytosterol biosynthesis: C24 sterol sidechain reductase 1 (SSR1), C4 sterol methyl oxidase 1, and 3 (SMO1 and SMO3). Surprisingly, the expression of all three genes correlated well with the cardenolide increase after the BTH treatment. Phylogenetic analysis showed that SSR1 is likely involved in both cholesterol and phytosterol biosynthesis. In addition, SMO1 is likely specific to phytosterol biosynthesis, and SMO3 is specific to cholesterol biosynthesis. These results suggest that stress-induced increase of cardenolides in foxglove may correlate with cholesterol and phytosterol biosynthesis. In summary, this work shows that cardenolides are important for stress responses in D. lanata and reveals a potential link between phytosterol and cardenolide biosynthesis.

The Role of Plant Progesterone in Regulating Growth, Development, and Biotic/Abiotic Stress Responses

Progesterone is a steroid hormone that performs important functions in mammals. However, studies on its physiological functions in plants have gradually increased in recent years. Therefore, this review summarizes the regulatory functions of progesterone on plant growth and development, as well as its response to stress. Moreover, the plant metabolic processes of progesterone are also discussed. Overall, progesterone is ubiquitous in plants and can regulate numerous plant physiological processes at low concentrations. Since progesterone shares similar characteristics with plant hormones, it is expected to become a candidate for plant hormone. However, most of the current research on progesterone in plants is limited to the physiological level, and more molecular level research is needed to clarify progesterone signaling pathways.

Characterization of stereospecific enoyl reductase ActVI-ORF2 for pyran ring formation in the actinorhodin biosynthesis of Streptomyces coelicolor A3(2)

Actinorhodin (ACT) is a benzoisochromanequinone antibiotic produced by Streptomyces coelicolor A3(2), which has served as a favored model organism for comprehensive studies of antibiotic biosynthesis and its regulation. (S)-DNPA undergoes various modifications as an intermediate in the ACT biosynthetic pathway, including enoyl reduction to DDHK. It has been suggested that actVI-ORF2 encodes an enoyl reductase (ER). However, its function has not been characterized in vitro. In this study, biochemical analysis of recombinant ActVI-ORF2 revealed that (S)-DNPA is converted to DDHK in a stereospecific manner with NADPH acting as a cofactor. (R)-DNPA was also reduced to 3-epi-DDHK with the comparable efficacy as (S)-DNPA, suggesting that the stereospecificity of ActVI-ORF2 was not affected by the stereochemistry at the C-3 of DNPA. ActVI-ORF2 is a new example of a discrete ER, which is distantly related to known ERs according to phylogenetic analysis.

Introduction of Chiral Centers to α- and/or β-Positions of Carbonyl Groups by Biocatalytic Asymmetric Reduction of α,β-Unsaturated Carbonyl Compounds

Biocatalytic asymmetric reductions of acyclic and cyclic α,β-unsaturated carbonyl compounds are favorable protocols for introduction of chiral centers to α- and/or β-positions of the carbonyl groups. Representative biocatalytic reductions of electron deficient olefins are compiled from a synthetic point of view according to compound types from the papers in 2012 to early 2022. Applications to syntheses of some enantiomericaly enriched perfumery ingredients are presented to show the feasibility of the biocatalytic reductions.

Identification and functional characterization of three iridoid synthases in Gardenia jasminoides

The discovery of three iridoid synthases (GjISY, GjISY2 and GjISY4) from Gardenia jasminoides and their functional characterization increase the understanding of iridoid scaffold/iridoid glycoside biosynthesis in iridoid-producing plants. Iridoids are a class of noncanonical monoterpenes that are found naturally in the plant kingdom mostly as glycosides. Over 40 iridoid glycosides (e.g., geniposide, gardenoside and shanzhiside) have been isolated from Gardenia jasminoides. They have multiple pharmacological properties and health-promoting effects. However, their biosynthetic pathway is poorly understood, and the iridoid synthase (ISY) responsible for the cyclization of the core scaffold remains unclear. In this study, three homologs of ISYs from G. jasminoides (GjISY, GjISY2 and GjISY4) were identified on the basis of transcriptomic data and functionally characterized. The genomic structure and intron-exon arrangement revealed that all three ISYs contained an intron. Biochemical assays indicated that all three recombinant enzymes reduced 8-oxogeranial to nepetalactol and its open forms (iridodials) as the products of the classical CrISY (Catharanthus roseus). In addition, all three enzymes reduced progesterone to 5-β-prognane-3,20-dione. However, only GjISY2 and GjISY4 reduced 2-cyclohexen-1-one to cyclohexanone. Overall, the GjISY2 expression levels in the flowers and fruits were similar to the GjISY and GjISY4 expression levels. By contrast, the GjISY2 expression levels in the upper and lower leaves were substantially higher than the GjISY and GjISY4 expression levels. Among the three, GjISY2 exhibited the highest catalytic efficiency for 8-oxogeranial. GjISY2 might be the major contributor to iridoid biosynthesis in G. jasminoides. Collectively, our results advance the understanding of iridoid scaffold/iridoid glycoside biosynthesis in G. jasminoides and provide a potential target for metabolic engineering and breeding.

Knockout of Arabidopsis thaliana VEP1, Encoding a PRISE (Progesterone 5β-Reductase/Iridoid Synthase-Like Enzyme), Leads to Metabolic Changes in Response to Exogenous Methyl Vinyl Ketone (MVK)

Small or specialized natural products (SNAPs) produced by plants vary greatly in structure and function, leading to selective advantages during evolution. With a limited number of genes available, a high promiscuity of the enzymes involved allows the generation of a broad range of SNAPs in complex metabolic networks. Comparative metabolic studies may help to understand why—or why not—certain SNAPs are produced in plants. Here, we used the wound-induced, vein patterning regulating VEP1 (AtStR1, At4g24220) and its paralogue gene on locus At5g58750 (AtStR2) from Arabidopsis to study this issue. The enzymes encoded by VEP1-like genes were clustered under the term PRISEs (progesterone 5β-reductase/iridoid synthase-like enzymes) as it was previously demonstrated that they are involved in cardenolide and/or iridoid biosynthesis in other plants. In order to further understand the general role of PRISEs and to detect additional more “accidental” roles we herein characterized A. thaliana steroid reductase 1 (AtStR1) and compared it to A. thaliana steroid reductase 2 (AtStR2). We used A. thaliana Col-0 wildtype plants as well as VEP1 knockout mutants and VEP1 knockout mutants overexpressing either AtStR1 or AtStR2 to investigate the effects on vein patterning and on the stress response after treatment with methyl vinyl ketone (MVK). Our results added evidence to the assumption that AtStR1 and AtStR2, as well as PRISEs in general, play specific roles in stress and defense situations and may be responsible for sudden metabolic shifts.

Plastidial Expression of 3β-Hydroxysteroid Dehydrogenase and Progesterone 5β-Reductase Genes Confer Enhanced Salt Tolerance in Tobacco

The short-chain dehydrogenase/reductase (SDR) gene family is widely distributed in all kingdoms of life. The SDR genes, 3β-hydroxysteroid dehydrogenase (3β-HSD) and progesterone 5-β-reductases (P5βR1, P5βR2) play a crucial role in cardenolide biosynthesis pathway in the Digitalis species. However, their role in plant stress, especially in salinity stress management, remains unexplored. In the present study, transplastomic tobacco plants were developed by inserting the 3β-HSD, P5βR1 and P5βR2 genes. The integration of transgenes in plastomes, copy number and transgene expression at transcript and protein level in transplastomic plants were confirmed by PCR, end-to-end PCR, qRT-PCR and Western blot analysis, respectively. Subcellular localization analysis showed that 3β-HSD and P5βR1 are cytoplasmic, and P5βR2 is tonoplast-localized. Transplastomic lines showed enhanced growth in terms of biomass and chlorophyll content compared to wild type (WT) under 300 mM salt stress. Under salt stress, transplastomic lines remained greener without negative impact on shoot or root growth compared to the WT. The salt-tolerant transplastomic lines exhibited enhanced levels of a series of metabolites (sucrose, glutamate, glutamine and proline) under control and NaCl stress. Furthermore, a lower Na+/K+ ratio in transplastomic lines was also observed. The salt tolerance, mediated by plastidial expression of the 3β-HSD, P5βR1 and P5βR2 genes, could be due to the involvement in the upregulation of nitrogen assimilation, osmolytes as well as lower Na+/K+ ratio. Taken together, the plastid-based expression of the SDR genes leading to enhanced salt tolerance, which opens a window for developing saline-tolerant plants via plastid genetic engineering.

RNAi-mediated gene knockdown of progesterone 5β-reductases in Digitalis lanata reduces 5β-cardenolide content

Studying RNAi-mediated DlP5βR1 and DlP5βR2 knockdown shoot culture lines of Digitalis lanata, we here provide direct evidence for the participation of PRISEs (progesterone 5β-reductase/iridoid synthase-like enzymes) in 5β-cardenolide formation. Progesterone 5β-reductases (P5βR) are assumed to catalyze the reduction of progesterone to 5β-pregnane-3,20-dione, which is a crucial step in the biosynthesis of the 5β-cardenolides. P5βRs are encoded by VEP1-like genes occurring ubiquitously in embryophytes. P5βRs are substrate-promiscuous enone-1,4-reductases recently termed PRISEs (progesterone 5β-reductase/iridoid synthase-like enzymes). Two PRISE genes, termed DlP5βR1 (AY585867.1) and DlP5βR2 (HM210089.1) were isolated from Digitalis lanata. To give experimental evidence for the participation of PRISEs in 5β-cardenolide formation, we here established several RNAi-mediated DlP5βR1 and DlP5βR2 knockdown shoot culture lines of D. lanata. Cardenolide contents were lower in D. lanata P5βR-RNAi lines than in wild-type shoots. We considered that the gene knockdowns may have had pleiotropic effects such as an increase in glutathione (GSH) which is known to inhibit cardenolide formation. GSH levels and expression of glutathione reductase (GR) were measured. Both were higher in the Dl P5βR-RNAi lines than in the wild-type shoots. Cardenolide biosynthesis was restored by buthionine sulfoximine (BSO) treatment in Dl P5βR2-RNAi lines but not in Dl P5βR1-RNAi lines. Since progesterone is a precursor of cardenolides but can also act as a reactive electrophile species (RES), we here discriminated between these by comparing the effects of progesterone and methyl vinyl ketone, a small RES but not a precursor of cardenolides. To the best of our knowledge, we here demonstrated for the first time that P5βR1 is involved in cardenolide formation. We also provide further evidence that PRISEs are also important for plants dealing with stress by detoxifying reactive electrophile species (RES).

21-Hydroxypregnane 21-O-malonylation, a crucial step in cardenolide biosynthesis, can be achieved by substrate-promiscuous BAHD-type phenolic glucoside malonyltransferases from Arabidopsis thaliana and homolog proteins from Digitalis lanata

Three putative 21-hydroxypregnane 21-O-malonyltransferases (21MaT) from Digitalis lanata were partially purified. Two of them were supposed to be BAHD-type enzymes. We were unable to purify them in quantities necessary for reliable sequencing. We identified two genes in A. thaliana coding for substrate-promiscuous BAHD-type phenolic glucoside malonyltransferases (AtPMaT1, AtPMaT2) and docked various 21-hydroxypregnanes into the substrate-binding site of a homology model built on the BAHD template 2XR7 (NtMaT1 from N. tabacum). Recombinant forms of Atpmat1 and Atpmat2 were expressed in E. coli and the recombinant enzymes characterized with regard to their substrate preferences. They were shown to malonylate various 21-hydroxypregnanes. The Atpmat1 sequence was used to identify candidate genes in Digitalis lanata (Dlmat1 to Dlmat4). Dlmat1 and Dlmat2 were also expressed in E. coli and shown to possess 21-hydroxypregnane 21-O-malonyltransferase activity.

Identification of key sites determining the cofactor specificity and improvement of catalytic activity of a steroid 5β-reductase from Capsella rubella

Progesterone 5β-reductases (P5βRs) are involved in 5β-cardenolide formation by stereo-specific reduction of the △^4,5 double bond of steroid precursors. In this study a steroid 5β-reductase was identified in Capsella rubella (CrSt5βR1) and its function in steroid 5β-reduction was validated experimentally. CrSt5βR1 is capable of enantioselectively reducing the activated CC bond of broad substrates such as steroids and enones by using NADPH as a cofactor and therefore has the potential as a biocatalyst in organic synthesis. However, for industrial purposes the cheaper NADH is the preferred cofactor. By applying rational design based on literature and complementary mutagenesis strategies, we successfully identified two key amino acid residues determining the cofactor specificity of the enzyme. The R63 K mutation enables the enzyme to convert progesterone to 5β-pregnane-3,20-dione with NADH as cofactor, whereas the wild-type CrSt5βR1 is strictly NADPH-dependent. By further introducing the R64H mutation, the double mutant R63K_R64H of CrSt5βR1 was shown to increase enzymatic activity by13.8-fold with NADH as a cofactor and to increase the NADH/NADPH conversion ratio by 10.9-fold over the R63 K single mutant. This finding was successfully applied to change the cofactor specificity and to improve activity of other members of the same enzyme family, AtP5βR and DlP5βR. CrSt5βR1 mutants are expected to have the potential for biotechnological applications in combination with the well-established NADH regeneration systems.

PRISEs (progesterone 5β-reductase and/or iridoid synthase-like 1,4-enone reductases): Catalytic and substrate promiscuity allows for realization of multiple pathways in plant metabolism

PRISEs (progesterone 5β-reductase and/or iridoid synthase-like 1,4-enone reductases) are involved in cardenolide and iridoid biosynthesis. We here investigated a PRISE (rAtSt5βR) from Arabidopsis thaliana, a plant producing neither cardenolides nor iridoids. The structure of rAtSt5βR was elucidated with X-ray crystallography and compared to the known structures of PRISEs from Catharanthus roseus (rCrISY) and Digitalis lanata (rDlP5βR). The three enzymes show a high degree of sequence and structure conservation in the active site. Amino acids previously considered to allow discrimination between progesterone 5β-reductase and iridoid synthase were interchanged among rAtSt5βR, rCrISY and rDlP5βR applying site-directed mutagenesis. Structural homologous substitutions had different effects, and changes in progesterone 5β-reductase and iridoid synthase activity were not correlated in all cases. Our results help to explain fortuitous emergence of metabolic pathways and product accumulation. The fact that PRISEs are found ubiquitously in spermatophytes insinuates that PRISEs might have a more general function in plant metabolism such as, for example, the detoxification of reactive carbonyl species.

Short-chain dehydrogenase/reductase governs steroidal specialized metabolites structural diversity and toxicity in the genus Solanum

Thousands of specialized, steroidal metabolites are found in a wide spectrum of plants. These include the steroidal glycoalkaloids (SGAs), produced primarily by most species of the genus Solanum, and metabolites belonging to the steroidal saponins class that are widespread throughout the plant kingdom. SGAs play a protective role in plants and have potent activity in mammals, including antinutritional effects in humans. The presence or absence of the double bond at the C-5,6 position (unsaturated and saturated, respectively) creates vast structural diversity within this metabolite class and determines the degree of SGA toxicity. For many years, the elimination of the double bond from unsaturated SGAs was presumed to occur through a single hydrogenation step. In contrast to this prior assumption, here, we show that the tomato GLYCOALKALOID METABOLISM25 (GAME25), a short-chain dehydrogenase/reductase, catalyzes the first of three prospective reactions required to reduce the C-5,6 double bond in dehydrotomatidine to form tomatidine. The recombinant GAME25 enzyme displayed 3β-hydroxysteroid dehydrogenase/Δ^5,4 isomerase activity not only on diverse steroidal alkaloid aglycone substrates but also on steroidal saponin aglycones. Notably, GAME25 down-regulation rerouted the entire tomato SGA repertoire toward the dehydro-SGAs branch rather than forming the typically abundant saturated α-tomatine derivatives. Overexpressing the tomato GAME25 in the tomato plant resulted in significant accumulation of α-tomatine in ripe fruit, while heterologous expression in cultivated eggplant generated saturated SGAs and atypical saturated steroidal saponin glycosides. This study demonstrates how a single scaffold modification of steroidal metabolites in plants results in extensive structural diversity and modulation of product toxicity.

Transcriptome and Metabolite analysis reveal candidate genes of the cardiac glycoside biosynthetic pathway from Calotropis procera

Calotropis procera is a medicinal plant of immense importance due to its pharmaceutical active components, especially cardiac glycosides (CG). As genomic resources for this plant are limited, the genes involved in CG biosynthetic pathway remain largely unknown till date. Our study on stage and tissue specific metabolite accumulation showed that CG's were maximally accumulated in stems of 3 month old seedlings. De novo transcriptome sequencing of same was done using high throughput Illumina HiSeq platform generating 44074 unigenes with average mean length of 1785 base pair. Around 66.6% of unigenes were annotated by using various public databases and 5324 unigenes showed significant match in the KEGG database involved in 133 different pathways of plant metabolism. Further KEGG analysis resulted in identification of 336 unigenes involved in cardenolide biosynthesis. Tissue specific expression analysis of 30 putative transcripts involved in terpenoid, steroid and cardenolide pathways showed a positive correlation between metabolite and transcript accumulation. Wound stress elevated CG levels as well the levels of the putative transcripts involved in its biosynthetic pathways. This result further validated the involvement of identified transcripts in CGs biosynthesis. The identified transcripts will lay a substantial foundation for further research on metabolic engineering and regulation of cardiac glycosides biosynthesis pathway genes.

Disorder prediction-based construct optimization improves activity and catalytic efficiency of Bacillus naganoensis pullulanase

Pullulanase is a well-known starch-debranching enzyme. However, the production level of pullulanase is yet low in both wide-type strains and heterologous expression systems. We predicted the disorder propensities of Bacillus naganoensis pullulanase (PUL) using the bioinformatics tool, Disorder Prediction Meta-Server. On the basis of disorder prediction, eight constructs, including PULΔN5, PULΔN22, PULΔN45, PULΔN64, PULΔN78 and PULΔN106 by deleting the first 5, 22, 45, 64, 78 and 106 residues from the N-terminus, and PULΔC9 and PULΔC36 by deleting the last 9 and 36 residues from the C-terminus, were cloned into the recombinant expression vector pET-28a-PelB and auto-induced in Escherichia coli BL21 (DE3) cells. All constructs were evaluated in production level, specific activities and kinetic parameters. Both PULΔN5 and PULΔN106 gave higher production levels of protein than the wide type and displayed increased specific activities. Kinetic studies showed that substrate affinities of the mutants were improved in various degrees and the catalytic efficiency of PULΔN5, PULΔN45, PULΔN78, PULΔN106 and PULΔC9 were enhanced. However, the truncated mutations did not change the advantageous properties of the enzyme involving optimum temperature and pH for further application. Therefore, Disorder prediction-based truncation would be helpful to efficiently improve the enzyme activity and catalytic efficiency.

Identification and Characterization of the Iridoid Synthase Involved in Oleuropein Biosynthesis in Olive (Olea europaea) Fruits

The secoiridoids are the main class of specialized metabolites present in olive (Olea europaea L.) fruit. In particular, the secoiridoid oleuropein strongly influences olive oil quality because of its bitterness, which is a desirable trait. In addition, oleuropein possesses a wide range of pharmacological properties, including antioxidant, anti-inflammatory, and anti-cancer activities. In accordance, obtaining high oleuropein varieties is a main goal of molecular breeding programs. Here we use a transcriptomic approach to identify candidate genes belonging to the secoiridoid pathway in olive. From these candidates, we have functionally characterized the olive homologue of iridoid synthase (OeISY), an unusual terpene cyclase that couples an NAD (P)H-dependent 1,4-reduction step with a subsequent cyclization, and we provide evidence that OeISY likely generates the monoterpene scaffold of oleuropein in olive fruits. OeISY, the first pathway gene characterized for this type of secoiridoid, is a potential target for breeding programs in a high value secoiridoid-accumulating species.

Structure of iridoid synthase in complex with NADP(+)/8-oxogeranial reveals the structural basis of its substrate specificity

Iridoid synthase (IS), as a vegetal enzyme belonging to the short-chain dehydrogenase/reductase (SDR) superfamily, produces the ring skeletons for downstream alkaloids with various pharmaceutical activities, including the commercially available antineoplastic agents, vinblastine and vincristine. Here, we present the crystal structures of IS in apo state and in complex with NADP(+)/8-oxogeranial, exhibiting an active center that lacks the classical Tyr/Lys/Ser triad spatially conserved in SDRs, with only the catalytically critical function of triad tyrosine remained in Tyr178. In consistent, mutation of Tyr178 to a phenylalanine residue significantly abolished the catalytic activity of IS. Within the substrate binding pocket, the linear-shaped 8-oxogeranial adopts an entirely extended conformation with its two aldehyde ends hydrogen-bonded to Tyr178-OH and Ser349-OH, respectively. In addition, the intermediate carbon chain of bound substrate is harbored by a well-ordered hydrophobic scaffold, involving residues Ile145, Phe149, Leu203, Met213, Phe342, Ile345 and Leu352. Mutagenesis studies showed that both Ser349 and the hydrophobic residues around are determinant to the substrate specificity and, consequently, the catalytic activity of IS. In contrast, the Gly150-Pro160 loop previously proposed as a factor involved in substrate binding might have very limited contribution, because the deletion of residues Ile151-His161 has only slight influence on the catalytic activity. We believe that the present work will help to elucidate the substrate specificity of IS and to integrate its detailed catalytic mechanism.

Steroid 5β-Reductase from Leaves of Vitis vinifera: Molecular Cloning, Expression, and Modeling

A steroid 5β-reductase gene corresponding to the hypothetical protein LOC100247199 from leaves of Vitis vinifera (var. 'Chardonnay') was cloned and overexpressed in Escherichia coli. The recombinant protein showed 5β-reductase activity when progesterone was used as a substrate. The reaction was stereoselective, producing only 5β-products such as 5β-pregnane-3,20-dione. Other small substrates (terpenoids and enones) were also accepted as substrates, indicating the highly promiscuous character of the enzyme class. Our results show that the steroid 5β-reductase gene, encoding an orthologous enzyme described as a key enzyme in cardenolide biosynthesis, is also expressed in leaves of the cardenolide-free plant V. vinifera. We emphasize the fact that, on some occasions, different reductases (e.g., progesterone 5β-reductase and monoterpenoid reductase) can also use molecules that are similar to the final products as a substrate. Therefore, in planta, the different reductases may contribute to the immense number of diverse small natural products finally leading to the flavor of wine.

Steroidogenesis in plants--Biosynthesis and conversions of progesterone and other pregnane derivatives

In plants androstanes, estranes, pregnanes and corticoids have been described. Sometimes 17β-estradiol, androsterone, testosterone or progesterone were summarized as sex hormones. These steroids influence plant development: cell divisions, root and shoot growth, embryo growth, flowering, pollen tube growth and callus proliferation. First reports on the effect of applicated substances and of their endogenous occurrence date from the early twenties of the last century. This caused later on doubts on the identity of the compounds. Best investigated is the effect of progesterone. Main steps of the progesterone biosynthetic pathway have been analyzed in Digitalis. Cholesterol-side-chain-cleavage, pregnenolone and progesterone formation as well as the stereospecific reduction of progesterone are described and the corresponding enzymes are presented. Biosynthesis of androstanes, estranes and corticoids is discussed. Possible progesterone receptors and physiological reactions on progesterone application are reviewed.

Iridoid synthase activity is common among the plant progesterone 5β-reductase family

Catharanthus roseus, the Madagascar periwinkle, synthesizes bioactive monoterpenoid indole alkaloids, including the anti-cancer drugs vinblastine and vincristine. The monoterpenoid branch of the alkaloid pathway leads to the secoiridoid secologanin and involves the enzyme iridoid synthase (IS), a member of the progesterone 5β-reductase (P5βR) family. IS reduces 8-oxogeranial to iridodial. Through transcriptome mining, we show that IS belongs to a family of six C. roseus P5βR genes. Characterization of recombinant CrP5βR proteins demonstrates that all but CrP5βR3 can reduce progesterone and thus can be classified as P5βRs. Three of them, namely CrP5βR1, CrP5βR2, and CrP5βR4, can also reduce 8-oxogeranial, pointing to a possible redundancy with IS (corresponding to CrP5βR5) in secoiridoid synthesis. In-depth functional analysis by subcellular protein localization, gene expression analysis, in situ hybridization, and virus-induced gene silencing indicate that besides IS, CrP5βR4 may also participate in secoiridoid biosynthesis. We cloned a set of P5βR genes from angiosperm plant species not known to produce iridoids and demonstrate that the corresponding recombinant proteins are also capable of using 8-oxogeranial as a substrate. This suggests that IS activity is intrinsic to angiosperm P5βR proteins and has evolved early during evolution.

Truncation of N-terminal regions of Digitalis lanata progesterone 5β-reductase alters catalytic efficiency and substrate preference

N-Terminal truncated forms of progesterone 5β-reductase (P5βR) were synthesized taking a full-length cDNA encoding for Digitalis lanata P5βR with a hexa-histidine tag attached at the C-terminus (rDlP5βRc) as the starting point. Four pETite-c-His/DlP5βR constructs coding for P5βR derivatives truncated in the N-terminal region, termed rDlP5βRcn-10, rDlP5βRcn-20, rDlP5βRcn-30, and rDlP5βRcn-40 were obtained by site-directed mutagenesis. The cDNAs coding for full-length rDlP5βRc, rDlP5βRcn-10 and rDlP5βRcn-20 were over-expressed in Escherichia coli and the respective enzymes were soluble and catalytically active (progesterone and 2-cyclohexen-1-one as substrates). GST-tagged recombinant DlP5βR (rDlP5βR-GST) and rDlP5βR-GSTr, with the GST-tag removed by protease treatment were produced as well and served as controls. The Km values and substrate preferences considerably differed between the various DlP5βR derivatives. As for the C-terminal His-tagged rDlP5βR the catalytic efficiency for progesterone was highest for the full-length rDlP5βRc whereas the N-terminal truncated forms preferred 2-cyclohexen-1-one as the substrate. Affinity tags and artifacts resulting from the cloning strategy used may alter substrate specificity. Therefore enzyme properties determined with recombinant proteins should not be used to infer in vivo scenarios and should be considered for each particular case.

Medicinal plants: a public resource for metabolomics and hypothesis development

Specialized compounds from photosynthetic organisms serve as rich resources for drug development. From aspirin to atropine, plant-derived natural products have had a profound impact on human health. Technological advances provide new opportunities to access these natural products in a metabolic context. Here, we describe a database and platform for storing, visualizing and statistically analyzing metabolomics data from fourteen medicinal plant species. The metabolomes and associated transcriptomes (RNAseq) for each plant species, gathered from up to twenty tissue/organ samples that have experienced varied growth conditions and developmental histories, were analyzed in parallel. Three case studies illustrate different ways that the data can be integrally used to generate testable hypotheses concerning the biochemistry, phylogeny and natural product diversity of medicinal plants. Deep metabolomics analysis of Camptotheca acuminata exemplifies how such data can be used to inform metabolic understanding of natural product chemical diversity and begin to formulate hypotheses about their biogenesis. Metabolomics data from Prunella vulgaris, a species that contains a wide range ofantioxidant, antiviral, tumoricidal and anti-inflammatory constituents, provide a case study of obtaining biosystematic and developmental fingerprint information from metabolite accumulation data in a little studied species. Digitalis purpurea, well known as a source of cardiac glycosides, is used to illustrate how integrating metabolomics and transcriptomics data can lead to identification of candidate genes encoding biosynthetic enzymes in the cardiac glycoside pathway. Medicinal Plant Metabolomics Resource (MPM) [1] provides a framework for generating experimentally testable hypotheses about the metabolic networks that lead to the generation of specialized compounds, identifying genes that control their biosynthesis and establishing a basis for modeling metabolism in less studied species. The database is publicly available and can be used by researchers in medicine and plant biology.

Vein Patterning 1-encoded progesterone 5β-reductase: activity-guided improvement of catalytic efficiency

Progesterone 5β-reductases (P5βR; EC 1.3.99.6) encoded by Vein Patterning 1 (VEP1) genes are capable of reducing the CC double-bond of a variety of enones enantioselectively. Sequence and activity data of orthologous P5βRs were used to define a set of residues possibly responsible for the large differences in enzyme activity seen between rAtSt5βR and rDlP5βR, recombinant forms of P5βRs from Arabidopsis thaliana and Digitalis lanata, respectively. Tyrosine-156, asparagine-205 and serine-248 were identified as hot spots in the rDlP5βR responsible for its low catalytic efficiency. These positions were individually substituted for amino acids found in the strong rAtSt5βR in the corresponding sites. Kinetic constants were determined for rDlP5βR and its mutants as well as for rAtSt5βR using progesterone and 2-cyclohexen-1-one as substrates. Enzyme mutants in which asparagine-205 was substituted for methionine or alanine showed considerably lower km and higher K(cat)/k(m) values than the wild-type DlP5βR, approaching the catalytic efficiency of strong P5βRs. The introduced mutations not only lead to an improved capability to reduce progesterone but also to altered substrate preference. Our findings provided structural insights into the differences seen among the natural P5βRs with regard to their substrate preferences and catalytic efficiencies.

The Vein Patterning 1 (VEP1) gene family laterally spread through an ecological network

Lateral gene transfer (LGT) is a major evolutionary mechanism in prokaryotes. Knowledge about LGT— particularly, multicellular— eukaryotes has only recently started to accumulate. A widespread assumption sees the gene as the unit of LGT, largely because little is yet known about how LGT chances are affected by structural/functional features at the subgenic level. Here we trace the evolutionary trajectory of VEin Patterning 1, a novel gene family known to be essential for plant development and defense. At the subgenic level VEP1 encodes a dinucleotide-binding Rossmann-fold domain, in common with members of the short-chain dehydrogenase/reductase (SDR) protein family. We found: i) VEP1 likely originated in an aerobic, mesophilic and chemoorganotrophic α-proteobacterium, and was laterally propagated through nets of ecological interactions, including multiple LGTs between phylogenetically distant green plant/fungi-associated bacteria, and five independent LGTs to eukaryotes. Of these latest five transfers, three are ancient LGTs, implicating an ancestral fungus, the last common ancestor of land plants and an ancestral trebouxiophyte green alga, and two are recent LGTs to modern embryophytes. ii) VEP1's rampant LGT behavior was enabled by the robustness and broad utility of the dinucleotide-binding Rossmann-fold, which provided a platform for the evolution of two unprecedented departures from the canonical SDR catalytic triad. iii) The fate of VEP1 in eukaryotes has been different in different lineages, being ubiquitous and highly conserved in land plants, whereas fungi underwent multiple losses. And iv) VEP1-harboring bacteria include non-phytopathogenic and phytopathogenic symbionts which are non-randomly distributed with respect to the type of harbored VEP1 gene. Our findings suggest that VEP1 may have been instrumental for the evolutionary transition of green plants to land, and point to a LGT-mediated ‘Trojan Horse’ mechanism for the evolution of bacterial pathogenesis against plants. VEP1 may serve as tool for revealing microbial interactions in plant/fungi-associated environments.

Occurrence of progesterone and related animal steroids in two higher plants

Previously, the presence of a wide variety of chemically diverse steroids has been identified in both flora and fauna. Despite the relatively small differences in chemical structures and large differences in physiological function of steroids, new discoveries indicate that plants and animals are more closely related than previously thought. In this regard, the present study gathers supporting evidence for shared phylogenetic roots of structurally similar steroids produced by these two eukaryotic taxa. Definitive proof for the presence of progesterone (3) in a vascular plant, Juglans regia, is provided. Additional evidence is gleaned from the characterization of five new plant steroids from Adonis aleppica: three 3-O-sulfated pregnenolones (6a/ b, 7), a sulfated H-5beta cardenolide, strophanthidin-3-O-sulfate (8), and spirophanthigenin (10), a novel C-18 oxygenated spirocyclic derivative of strophanthidin. The ab initio isolation and structure elucidation (NMR, MS) of these genuine minor plant steroids offers information on preparative metabolomic profiling at the ppm level and provides striking evidence for the conserved structural space of pregnanes and its congeners across the phylogenetic tree.

Synthesis of a putative substrate for malonyl-coenzyme A: 21-hydroxypregnane 21-O-malonyltransferase and development of an HPLC method for the quantification of the enzyme reaction

The butenolide ring is the main common characteristic of all cardenolides. Its formation is supposed to be initiated by the transfer of a malonyl moiety from malonyl-coenzyme A to an appropriate 21-hydroxypregnane. A new, reliable, fast and sensitive method to determine malonyl-coenzyme A: 21-hydroxypregnane 21-O-malonyltransferase activity had to be developed since previous attempts employing HPLC, TLC or GC did not prove successful. A surrogate substrate was synthesized containing a side chain resembling the sugar side chain attached to C-3 of putative cardenolide precursors and containing a chromophor allowing UV detection. 3beta-benzoyloxy-5beta-pregnane-14beta,21-dihydroxy-20-one and its 21-O-malonylated derivative were synthesized, the latter being the expected product of the enzyme reaction. The new substrate was well accepted by the enzyme. An HPLC method has been established to detect and quantify 3beta-benzoyloxy-5beta-pregnane-14beta,21-dihydroxy-20-one and its 21-O-malonylated derivative, 3beta-benzoyloxy-5beta-pregnane-14beta-hydroxy-20-one 21-O-malonylhemiester. The method was validated.

Progesterone: its occurrence in plants and involvement in plant growth

Progesterone is a mammalian gonadal hormone. In the current study, we identified and quantified progesterone in a range of higher plants by using GC-MS and examined its effects on the vegetative growth of plants. The growth of Arabidopsis (Arabidopsis thaliana) seedlings was promoted by progesterone at low concentrations but suppressed at higher concentrations under both light and dark growth conditions. The growth of the gibberellin-deficient mutant lh of pea (Pisum sativum) was also promoted by progesterone. An earlier study demonstrated that progesterone binds to MEMBRANE STEROID BINDING PROTEIN 1 (MSBP1) of Arabidopsis. In this work, we cloned the homologous genes of Arabidopsis, MSBP2 and STEROID BINDING PROTEIN (SBP), as well as of rice (Oryza sativa), OsMSBP1, OsMSBP2 and OsSBP and examined their expression in plant tissues. All of these genes, except OsMSBP1, were expressed abundantly in plant tissues. The roles of progesterone in plant growth were also discussed.

Crystallization and preliminary crystallographic analysis of selenomethionine-labelled progesterone 5beta-reductase from Digitalis lanata Ehrh

Progesterone 5beta-reductase (5beta-POR) catalyzes the reduction of progesterone to 5beta-pregnane-3,20-dione and is the first stereospecific enzyme in the putative biosynthetic pathway of Digitalis cardenolides. Selenomethionine-derivatized 5beta-POR from D. lanata was successfully overproduced and crystallized. The crystals belong to space group P4(3)2(1)2, with unit-cell parameters a = 71.73, c = 186.64 A. A MAD data set collected at 2.7 A resolution allowed the identification of six out of eight possible Se-atom positions. A first inspection of the MAD-phased electron-density map shows that 5beta-POR is a Rossmann-type reductase and the quality of the map is such that it is anticipated that a complete atomic model of 5beta-POR will readily be built.