Dynamic Conformational States Dictate Selectivity toward the Native Substrate in a Substrate-Permissive Acyltransferase

Protein-Chlorogenic Acid Interactions: Mechanisms, Characteristics, and Potential Food Applications

The interactions between proteins and chlorogenic acid (CGA) have gained significant attention in recent years, not only as a promising approach to modify the structural and techno-functional properties of proteins but also to enhance their bioactive potential in food systems. These interactions can be divided into covalent (chemical or irreversible) and non-covalent (physical or reversible) linkages. Mechanistically, CGA forms covalent bonds with nucleophilic amino acid residues of proteins by alkaline, free radical, and enzymatic approaches, leading to changes in protein structure and functionality, such as solubility, emulsification properties, and antioxidant activity. In addition, the protein-CGA complexes can be obtained by hydrogen bonds, hydrophobic and electrostatic interactions, and van der Waals forces, each offering unique advantages and outcomes. This review highlights the mechanism of these interactions and their importance in modifying the structural, functional, nutritional, and physiological attributes of animal- and plant-based proteins. Moreover, the potential applications of these protein-CGA conjugates/complexes are explored in various food systems, such as beverages, films and coatings, emulsion-based delivery systems, and so on. Overall, this literature review provides an in-depth overview of protein-CGA interactions, offering valuable insights for future research to develop novel protein-based food and non-food products with improved nutritional and functional characteristics.

Genome-wide association identifies a BAHD acyltransferase activity that assembles an ester of glucuronosylglycerol and phenylacetic acid

Genome-wide association studies (GWAS) are an effective approach to identify new specialized metabolites and the genes involved in their biosynthesis and regulation. In this study, GWAS of Arabidopsis thaliana soluble leaf and stem metabolites identified alleles of an uncharacterized BAHD-family acyltransferase (AT5G57840) associated with natural variation in three structurally related metabolites. These metabolites were esters of glucuronosylglycerol, with one metabolite containing phenylacetic acid as the acyl component of the ester. Knockout and overexpression of AT5G57840 in Arabidopsis and heterologous overexpression in Nicotiana benthamiana and Escherichia coli demonstrated that it is capable of utilizing phenylacetyl-CoA as an acyl donor and glucuronosylglycerol as an acyl acceptor. We, thus, named the protein Glucuronosylglycerol Ester Synthase (GGES). Additionally, phenylacetyl glucuronosylglycerol increased in Arabidopsis CYP79A2 mutants that overproduce phenylacetic acid and was lost in knockout mutants of UDP-sulfoquinovosyl: diacylglycerol sulfoquinovosyl transferase, an enzyme required for glucuronosylglycerol biosynthesis and associated with glycerolipid metabolism under phosphate-starvation stress. GGES is a member of a well-supported clade of BAHD family acyltransferases that arose by duplication and neofunctionalized during the evolution of the Brassicales within a larger clade that includes HCT as well as enzymes that synthesize other plant-specialized metabolites. Together, this work extends our understanding of the catalytic diversity of BAHD acyltransferases and uncovers a pathway that involves contributions from both phenylalanine and lipid metabolism.

Discovering a mitochondrion-localized BAHD acyltransferase involved in calystegine biosynthesis and engineering the production of 3β-tigloyloxytropane

Solanaceous plants produce tropane alkaloids (TAs) via esterification of 3α- and 3β-tropanol. Although littorine synthase is revealed to be responsible for 3α-tropanol esterification that leads to hyoscyamine biosynthesis, the genes associated with 3β-tropanol esterification are unknown. Here, we report that a BAHD acyltransferase from Atropa belladonna, 3β-tigloyloxytropane synthase (TS), catalyzes 3β-tropanol and tigloyl-CoA to form 3β-tigloyloxytropane, the key intermediate in calystegine biosynthesis and a potential drug for treating neurodegenerative disease. Unlike other cytosolic-localized BAHD acyltransferases, TS is localized to mitochondria. The catalytic mechanism of TS is revealed through molecular docking and site-directed mutagenesis. Subsequently, 3β-tigloyloxytropane is synthesized in tobacco. A bacterial CoA ligase (PcICS) is found to synthesize tigloyl-CoA, an acyl donor for 3β-tigloyloxytropane biosynthesis. By expressing TS mutant and PcICS, engineered Escherichia coli synthesizes 3β-tigloyloxytropane from tiglic acid and 3β-tropanol. This study helps to characterize the enzymology and chemodiversity of TAs and provides an approach for producing 3β-tigloyloxytropane.

Uncovering the Role of Hydroxycinnamoyl Transferase in Boosting Chlorogenic Acid Accumulation in Carthamus tinctorius Cells under Methyl Jasmonate Elicitation

Chlorogenic acids (CGAs) are bioactive compounds widely used in the food, pharmaceutical, and cosmetic industries. Carthamus tinctorius is an important economic crop, and its suspension cells are rich in CGAs. However, little is known about the biosynthesis and regulation of CGAs in Carthamus tinctorius cells. This study first elucidated the regulatory mechanism of CGA biosynthesis in methyl jasmonate (MeJA)-treated Carthamus tinctorius cells and the role of the MeJA-responsive hydroxycinnamoyl transferase (HCT) gene in enhancing their CGA accumulation. Firstly, temporal changes in intracellular metabolites showed that MeJA increased the intracellular CGA content up to 1.61-fold to 100.23 mg·g-1. Meanwhile, 31 primary metabolites showed significant differences, with 6 precursors related to increasing CGA biosynthesis. Secondly, the transcriptome data revealed 3637 new genes previously unannotated in the Carthamus tinctorius genome and 3653 differentially expressed genes. The genes involved in the plant signaling pathway and the biosynthesis of CGAs and their precursors showed a general up-regulation, especially the HCT gene family, which ultimately promoted CGA biosynthesis. Thirdly, the expression of a newly annotated and MeJA-responsive HCT gene (CtHCT, CtNewGene_3476) was demonstrated to be positively correlated with CGA accumulation in the cells, and transient overexpression of CtHCT enhanced CGA accumulation in tobacco. Finally, in vitro catalysis kinetics and molecular docking simulations revealed the ability and mechanism of the CtHCT protein to bind to various substrates and catalyze the formation of four hydroxycinnamic esters, including CGAs. These findings strengthened our understanding of the regulatory mechanism of CGA biosynthesis, thereby providing theoretical support for the efficient production of CGAs.

Analysing a Group of Homologous BAHD Enzymes Provides Insights into the Evolutionary Transition of Rosmarinic Acid Synthases from Hydroxycinnamoyl-CoA:Shikimate/Quinate Hydroxycinnamoyl Transferases

The interplay of various enzymes and compounds gives rise to the intricate secondary metabolic networks observed today. However, the current understanding of their formation and expansion remains limited. BAHD acyltransferases play important roles in the biosynthesis of numerous significant secondary metabolites. In plants, they are widely distributed and exhibit a diverse range of activities. Among them, rosmarinic acid synthase (RAS) and hydroxycinnamoyl-CoA:shikimate/quinate hydroxycinnamoyl transferase (HCT) have gained significant recognition and have been extensively investigated as prominent members of the BAHD acyltransferase family. Here, we conducted a comprehensive study on a unique group of RAS homologous enzymes in Mentha longifolia that display both catalytic activities and molecular features similar to HCT and Lamiaceae RAS. Subsequent phylogenetic and comparative genome analyses revealed their derivation from expansion events within the HCT gene family, indicating their potential as collateral branches along the evolutionary trajectory, leading to Lamiaceae RAS while still retaining certain ancestral vestiges. This discovery provides more detailed insights into the evolution from HCT to RAS. Our collective findings indicate that gene duplication is the driving force behind the observed evolutionary pattern in plant-specialized enzymes, which probably originated from ancestral enzyme promiscuity and were subsequently shaped by principles of biological adaptation.

Acyltransferase families that act on thioesters: Sequences, structures, and mechanisms

Acyltransferases (AT) are enzymes that catalyze the transfer of acyl group to a receptor molecule. This review focuses on ATs that act on thioester-containing substrates. Although many ATs can recognize a wide variety of substrates, sequence similarity analysis allowed us to classify the ATs into fifteen distinct families. Each AT family is originated from enzymes experimentally characterized to have AT activity, classified according to sequence similarity, and confirmed with tertiary structure similarity for families that have crystallized structures available. All the sequences and structures of the AT families described here are present in the thioester-active enzyme (ThYme) database. The AT sequences and structures classified into families and available in the ThYme database could contribute to enlightening the understanding acyl transfer to thioester-containing substrates, most commonly coenzyme A, which occur in multiple metabolic pathways, mostly with fatty acids.

Integrated Analysis of Metabolome and Transcriptome Revealed Different Regulatory Networks of Metabolic Flux in Tea Plants [Camellia sinensis (L.) O. Kuntze] with Varied Leaf Colors

Leaf color variations in tea plants were widely considered due to their attractive phenotypes and characteristic flavors. The molecular mechanism of color formation was extensively investigated. But few studies focused on the transformation process of leaf color change. In this study, four strains of 'Baijiguan' F1 half-sib generation with similar genetic backgrounds but different colors were used as materials, including Green (G), Yellow-Green (Y-G), Yellow (Y), and Yellow-Red (Y-R). The results of broadly targeted metabolomics showed that 47 metabolites were differentially accumulated in etiolated leaves (Y-G, Y, and Y-R) as compared with G. Among them, lipids were the main downregulated primary metabolites in etiolated leaves, which were closely linked with the thylakoid membrane and chloroplast structure. Flavones and flavonols were the dominant upregulated secondary metabolites in etiolated leaves, which might be a repair strategy for reducing the negative effects of dysfunctional chloroplasts. Further integrated analysis with the transcriptome indicated different variation mechanisms of leaf phenotype in Y-G, Y, and Y-R. The leaf color formation of Y-G and Y was largely determined by the increased content of eriodictyol-7-O-neohesperidoside and the enhanced activities of its modification process, while the color formation of Y-R depended on the increased contents of apigenin derivates and the vigorous processes of their transportation and transcription factor regulation. The key candidate genes, including UDPG, HCT, CsGSTF1, AN1/CsMYB75, and bHLH62, might play important roles in the flavonoid pathway.

Loop Evolutionary Patterns Shape Catalytic Efficiency of TRI101/201 for Trichothecenes: Insights into Protein-Substrate Interactions

Trichothecenes are highly toxic mycotoxins produced by Fusarium fungi, while TRI101/201 family enzymes play a crucial role in detoxification through acetylation. Studies on the substrate specificity and catalytic kinetics of TRI101/201 have revealed distinct kinetic characteristics, with significant differences observed in catalytic efficiency toward deoxynivalenol, while the catalytic efficiency for T-2 toxin remains relatively consistent. In this study, we used structural bioinformatics analysis and a molecular dynamics simulation workflow to investigate the mechanism underlying the differential catalytic activity of TRI101/201. The findings revealed that the binding stability between trichothecenes and TRI101/201 hinges primarily on a hydrophobic cage structure within the binding site. An intrinsic disordered loop, termed loop cover, defined the evolutionary patterns of the TRI101/201 protein family that are categorized into four subfamilies (V1/V2/V3/M). Furthermore, the unique loop displayed different conformations among these subfamilies' structures, which served to disrupt (V1/V2/V3) or reinforce (M) the hydrophobic cages. The disrupted cages enhanced the water exposure of the hydrophilic moieties of substrates like deoxynivalenol and thereby hindered their binding to the catalytic sites of V-type enzymes. In contrast, this water exposure does not affect substrates like T-2 toxin, which have more hydrophobic substituents, resulting in a comparable catalytic efficiency of both V- and M-type enzymes. Overall, our studies provide theoretical support for understanding the catalytic mechanism of TRI101/201, which shows how an intrinsic disordered loop could impact the protein-ligand binding and suggests a direction for rational protein design in the future.

Family characteristics, phylogenetic reconstruction, and potential applications of the plant BAHD acyltransferase family

The BAHD acyltransferase family is a class of proteins in plants that can acylate a variety of primary and specialized secondary metabolites. The typically acylated products have greatly improved stability, lipid solubility, and bioavailability and thus show significant differences in their physicochemical properties and pharmacological activities. Here, we review the protein structure, catalytic mechanism, and phylogenetic reconstruction of plant BAHD acyltransferases to describe their family characteristics, acylation reactions, and the processes of potential functional differentiation. Moreover, the potential applications of the BAHD family in human activities are discussed from the perspectives of improving the quality of economic plants, enhancing the efficacy of medicinal plants, improving plant biomass for use in biofuel, and promoting stress resistance of land plants. This review provides a reference for the research and production of plant BAHD acyltransferases.

Characterization and structure-based protein engineering of a regiospecific saponin acetyltransferase from Astragalus membranaceus

Acetylation contributes to the bioactivity of numerous medicinally important natural products. However, little is known about the acetylation on sugar moieties. Here we report a saponin acetyltransferase from Astragalus membranaceus. AmAT7-3 is discovered through a stepwise gene mining approach and characterized as the xylose C3'/C4'-O-acetyltransferse of astragaloside IV (1). To elucidate its catalytic mechanism, complex crystal structures of AmAT7-3/1 and AmAT7-3A310G/1 are obtained, which reveal a large active pocket decided by a specific sequence AADAG. Combining with QM/MM computation, the regiospecificity of AmAT7-3 is determined by sugar positioning modulated by surrounding amino acids including #A310 and #L290. Furthermore, a small mutant library is built using semi-rational design, where variants A310G and A310W are found to catalyze specific C3'-O and C4'-O acetylation, respectively. AmAT7-3 and its variants are also employed to acetylate other bioactive saponins. This work expands the understanding of saponin acetyltransferases, and provide efficient catalytic tools for saponin acetylation.

Phylogenomic analyses across land plants reveals motifs and coexpression patterns useful for functional prediction in the BAHD acyltransferase family

The BAHD acyltransferase family is one of the largest enzyme families in flowering plants, containing dozens to hundreds of genes in individual genomes. Highly prevalent in angiosperm genomes, members of this family contribute to several pathways in primary and specialized metabolism. In this study, we performed a phylogenomic analysis of the family using 52 genomes across the plant kingdom to gain deeper insights into its functional evolution and enable function prediction. We found that BAHD expansion in land plants was associated with significant changes in various gene features. Using pre-defined BAHD clades, we identified clade expansions in different plant groups. In some groups, these expansions coincided with the prominence of metabolite classes such as anthocyanins (flowering plants) and hydroxycinnamic acid amides (monocots). Clade-wise motif-enrichment analysis revealed that some clades have novel motifs fixed on either the acceptor or the donor side, potentially reflecting historical routes of functional evolution. Co-expression analysis in rice and Arabidopsis further identified BAHDs with similar expression patterns, however, most co-expressed BAHDs belonged to different clades. Comparing BAHD paralogs, we found that gene expression diverges rapidly after duplication, suggesting that sub/neo-functionalization of duplicate genes occurs quickly via expression diversification. Analyzing co-expression patterns in Arabidopsis in conjunction with orthology-based substrate class predictions and metabolic pathway models led to the recovery of metabolic processes of most of the already-characterized BAHDs as well as definition of novel functional predictions for some uncharacterized BAHDs. Overall, this study provides new insights into the evolution of BAHD acyltransferases and sets up a foundation for their functional characterization.

Emergence of a proton exchange-based isomerization and lactonization mechanism in the plant coumarin synthase COSY

Plants contain rapidly evolving specialized enzymes that support the biosynthesis of functionally diverse natural products. In coumarin biosynthesis, a BAHD acyltransferase-family enzyme COSY was recently discovered to accelerate coumarin formation as the only known BAHD enzyme to catalyze an intramolecular acyl transfer reaction. Here we investigate the structural and mechanistic basis for COSY's coumarin synthase activity. Our structural analyses reveal an unconventional active-site configuration adapted to COSY's specialized activity. Through mutagenesis studies and deuterium exchange experiments, we identify a unique proton exchange mechanism at the α-carbon of the o-hydroxylated trans-hydroxycinnamoyl-CoA substrates during the catalytic cycle of COSY. Quantum mechanical cluster modeling and molecular dynamics further support this key mechanism for lowering the activation energy of the rate-limiting trans-to-cis isomerization step in coumarin production. This study unveils an unconventional catalytic mechanism mediated by a BAHD-family enzyme, and sheds light on COSY's evolutionary origin and its recruitment to coumarin biosynthesis in eudicots.

BAHD acyltransferase induced by histone deacetylase inhibitor catalyzes 3-O-hydroxycinnamoylquinic acid formation in bamboo cells

Suspension-cultured cells of a bamboo species (Bambusa multiplex; Bm) produce 3-O-feruloylquinic acid (3-FQA) and 3-O-p-coumaroylquinic acid (3-pCQA) by treatment with the histone deacetylase inhibitor suberoyl bis-hydroxamic acid (SBHA). Acyltransferases catalyzing the formation of 5-O-hydroxycinnamoylquinic acid esters by transesterification from hydroxycinnamoyl-CoAs to the C-5 hydroxy group of quinic acid (hydroxycinnamoyl-CoA:quinate hydroxycinnamoyltransferase, HQT) have been identified in the biosynthesis of chlorogenic acids and monolignols; however, an HQT that catalyzes the acylation of the C-3 hydroxy group of quinic acid has not been identified previously. In the present study, we purified a native HQT from SBHA-treated Bm cells. The purified enzyme preferentially accepted feruloyl-/p-coumaroyl-CoAs as acyl-donors and quinic acid as the acyl-acceptor, and the enzyme specifically formed 3-FQA and 3-pCQA but not 5-O-hydroxycinnamoylquinic acid esters or esters with shikimic acid. A cDNA (BmHQT1) encoding this HQT was isolated. Although BmHQT1 is a phylogenetically unique member of the BAHD acyltransferase superfamily that does not cluster with other HQTs, functional characterization of the recombinant enzyme verified that BmHQT1 catalyzes the regiospecific formation of 3-O-hydroxycinnamoylquinic acid esters. Transcript levels of BmHQT1 markedly increased in Bm cells cultured in the presence of SBHA. Moreover, elevated acetylation levels of histone H3 were observed in the coding region of BmHQT1 in the presence of SBHA, indicating that the induced accumulation of 3-FQA/3-pCQA by SBHA is caused by transcriptional activation of BmHQT1 by the action of SBHA as a histone deacetylase inhibitor. The results demonstrate the utility of HDAC inhibitors for discovery of cryptic secondary metabolites and unknown biosynthetic enzymes.

Orthology-based analysis helps map evolutionary diversification and predict substrate class use of BAHD acyltransferases

Large enzyme families catalyze metabolic diversification by virtue of their ability to use diverse chemical scaffolds. How enzyme families attain such functional diversity is not clear. Furthermore, duplication and promiscuity in such enzyme families limits their functional prediction, which has produced a burgeoning set of incompletely annotated genes in plant genomes. Here, we address these challenges using BAHD acyltransferases as a model. This fast-evolving family expanded drastically in land plants, increasing from one to five copies in algae to approximately 100 copies in diploid angiosperm genomes. Compilation of >160 published activities helped visualize the chemical space occupied by this family and define eight different classes based on structural similarities between acceptor substrates. Using orthologous groups (OGs) across 52 sequenced plant genomes, we developed a method to predict BAHD acceptor substrate class utilization as well as origins of individual BAHD OGs in plant evolution. This method was validated using six novel and 28 previously characterized enzymes and helped improve putative substrate class predictions for BAHDs in the tomato genome. Our results also revealed that while cuticular wax and lignin biosynthetic activities were more ancient, anthocyanin acylation activity was fixed in BAHDs later near the origin of angiosperms. The OG-based analysis enabled identification of signature motifs in anthocyanin-acylating BAHDs, whose importance was validated via molecular dynamic simulations, site-directed mutagenesis and kinetic assays. Our results not only describe how BAHDs contributed to evolution of multiple chemical phenotypes in the plant world but also propose a biocuration-enabled approach for improved functional annotation of plant enzyme families.

Probing the transcriptome of Boehmeria nivea reveals candidate genes associated with the biosynthesis of chlorogenic acid

Boehmeria nivea (L.) Gaudich is used in traditional Chinese medicine. Chlorogenic acids are major medically active components of Boehmeria nivea, which can be used clinically to treat hyperglycemia, pneumonia, and cancer. To identify the genes involved in chlorogenic acid biosynthesis, we analyzed transcriptome data from leaf, root, and stem tissues of Boehmeria nivea using the Illumina Hi-Seq 4000 platform. A total of 146,790 unigenes were obtained from Boehmeria nivea, of which 106,786 were annotated in public databases. In analyses of the KEGG (Kyoto Encyclopedia of Genes and Genome) database, 484 unigenes that encode the five key enzymes involved in chlorogenic acid biosynthesis were identified, and shikimate O-hydroxycinnamoyl transferase was spatially simulated. Some of these key enzyme unigenes expression levels were verified by RT-qPCR (real-time quantitative Polymerase Chain Reaction). Furthermore, multiple genes encoding plant resistance proteins or transcription factors were identified and analyzed. Differentially expressed genes were identified by performing pairwise comparison of genes between tissues. This study increases the number of public transcript datasets of this species and identifies candidate genes related to the biosynthesis of chlorogenic acid, laying a foundation for the further exploration of this pathway in Boehmeria nivea.

The Biological Activity Mechanism of Chlorogenic Acid and Its Applications in Food Industry: A Review

Chlorogenic acid (CGA), also known as coffee tannic acid and 3-caffeoylquinic acid, is a water-soluble polyphenolic phenylacrylate compound produced by plants through the shikimic acid pathway during aerobic respiration. CGA is widely found in higher dicotyledonous plants, ferns, and many Chinese medicine plants, which enjoy the reputation of “plant gold.” We have summarized the biological activities of CGA, which are mainly shown as anti-oxidant, liver and kidney protection, anti-bacterial, anti-tumor, regulation of glucose metabolism and lipid metabolism, anti-inflammatory, protection of the nervous system, and action on blood vessels. We further determined the main applications of CGA in the food industry, including food additives, food storage, food composition modification, food packaging materials, functional food materials, and prebiotics. With a view to the theoretical improvement of CGA, biological activity mechanism, and subsequent development and utilization provide reference and scientific basis.

Substrate promiscuity of acyltransferases contributes to the diversity of hydroxycinnamic acid derivatives in purple coneflower

High pliability and promiscuity are observed widely exist in plant specialized metabolism, especially the hydroxycinnamic acid metabolism. Here, we identified an addition BAHD acyltransferase (EpHMT) that catalyzes phaselic acid biosynthesis and found that the substrate promiscuities of identified BAHD and SCPL acyltransferases are responsible for the diversity of hydroxycinnamic acid derivatives in purple coneflower.

Rerouting of the lignin biosynthetic pathway by inhibition of cytosolic shikimate recycling in transgenic hybrid aspen

Lignin is a phenolic polymer deposited in the plant cell wall, and is mainly polymerized from three canonical monomers (monolignols), i.e. p-coumaryl, coniferyl and sinapyl alcohols. After polymerization, these alcohols form different lignin substructures. In dicotyledons, monolignols are biosynthesized from phenylalanine, an aromatic amino acid. Shikimate acts at two positions in the route to the lignin building blocks. It is part of the shikimate pathway that provides the precursor for the biosynthesis of phenylalanine, and is involved in the transesterification of p-coumaroyl-CoA to p-coumaroyl shikimate, one of the key steps in the biosynthesis of coniferyl and sinapyl alcohols. The shikimate residue in p-coumaroyl shikimate is released in later steps, and the resulting shikimate becomes available again for the biosynthesis of new p-coumaroyl shikimate molecules. In this study, we inhibited cytosolic shikimate recycling in transgenic hybrid aspen by accelerated phosphorylation of shikimate in the cytosol through expression of a bacterial shikimate kinase (SK). This expression elicited an increase in p-hydroxyphenyl units of lignin and, by contrast, a decrease in guaiacyl and syringyl units. Transgenic plants with high SK activity produced a lignin content comparable to that in wild-type plants, and had an increased processability via enzymatic saccharification. Although expression of many genes was altered in the transgenic plants, elevated SK activity did not exert a significant effect on the expression of the majority of genes responsible for lignin biosynthesis. The present results indicate that cytosolic shikimate recycling is crucial to the monomeric composition of lignin rather than for lignin content.

Crystal structure of the plant feruloyl-coenzyme A monolignol transferase provides insights into the formation of monolignol ferulate conjugates

Lignin is a highly complex phenolic polymer which is essential for plants, but also makes it difficult for industrial processing. Engineering lignin by introducing relatively labile linkages into the lignin backbone can render it more amenable to chemical depolymerization. It has been reported that introducing a feruloyl-coenzyme A monolignol transferase from Angelica sinensis (AsFMT) into poplar could incorporate monolignol ferulate conjugates (ML-FAs) into lignin polymers, suggesting a promising way to manipulate plants for readily deconstructing. FMT catalyzes a reaction between monolignols and feruloyl-CoA to produce ML-FAs and free CoA-SH. However, the mechanisms of substrate specificity and catalytic process of FMT remains poorly understood. Here we report the structure of AsFMT, which adopts a typical fold of BAHD acyltransferase family. Structural comparisons with other BAHD homologs reveal several unique structural features of AsFMT, different from those of the BAHD homologs. Further molecular docking studies showed that T375 in AsFMT may function as an oxyanion hole to stabilize the reaction intermediate and also proposed a role of H278 in the binding of the nucleophilic hydroxyl group of monolignols. Together, this study provides important structural insights into the reactions catalyzed by AsFMT and will shed light on its future application in lignin engineering.

The BAHD Gene Family in Cacao (Theobroma cacao, Malvaceae): Genome-Wide Identification and Expression Analysis

The benzyl alcohol O-acetyl transferase, anthocyanin O-hydroxycinnamoyl transferase, N-hydroxycinnamoyl anthranilate benzoyl transferase, and deacetylvindoline 4-O-acetyltransferase (BAHD) enzymes play a critical role in regulating plant metabolites and affecting cell stability. In the present study, members of the BAHD gene family were recognized in the genome of Theobroma cacao and characterized using various bioinformatics tools. We found 27 non-redundant putative tcBAHD genes in cacao for the first time. Our findings indicate that tcBAHD genes are diverse based on sequence structure, physiochemical properties, and function. When analyzed with BAHDs of Gossypium raimondii and Corchorus capsularis clustered into four main groups. According to phylogenetic analysis, BAHD genes probably evolved drastically after their divergence. The divergence time of duplication events with purifying selection pressure was predicted to range from 1.82 to 15.50 MYA. Pocket analysis revealed that serine amino acid is more common in the binding site than other residuals, reflecting its key role in regulating the activity of tcBAHDs. Furthermore, cis-acting elements related to the responsiveness of stress and hormone, particularly ABA and MeJA, were frequently observed in the promoter region of tcBAHD genes. RNA-seq analysis further illustrated that tcBAHD13 and tcBAHD26 are involved in response to Phytophthora megakarya fungi. In conclusion, it is likely that evolutionary processes, such as duplication events, have caused high diversity in the structure and function of tcBAHD genes.

Molecular and structural characterization of agmatine coumaroyltransferase in Triticeae, the key regulator of hydroxycinnamic acid amide accumulation

Hydroxycinnamic acid amides (HCAAs) are involved in stress-induced defense in many plant species. Barley accumulates high concentrations of HCAAs irrespective of exogenous stressors, while other major cereals such as wheat and rice accumulate relatively low levels of HCAAs in intact tissues. The primary HCAA species in barley are biosynthesized by agmatine p-coumaroyltransferase (ACT), an N-acyltransferase of the BAHD superfamily. However, the molecular basis underlying barley's uniquely high HCAA accumulation has not been elucidated, and information regarding the structural details of BAHD N-acyltransferases is limited. Hence, we aimed to investigate the ACTs of family Poaceae. We isolated ACT (-like) genes, including those previously undescribed, and investigated their enzymatic and genetic features. All the identified enzymes belonged to clade IVa of the BAHD superfamily. The barley and wheat ACTs were further categorized, based on catalytic properties and primary structures, into ACT1 and ACT2 groups, the encoding loci of which are neighbors on the same chromosome. While all ACTs exhibited similar K_m values for CoA-thioesters (acyl-group donors), members of the ACT1 group showed a distinctly higher affinity for agmatine (acyl-acceptor). Among the ACTs tested, an ACT isozyme in barley (HvACT1-1) showed the highest catalytic efficiency and transcript level, indicating that ACT regulates high-level HCAA accumulation in barley. For further enzymatic characterization of the ACTs, we crystalized wheat ACT2 (TaACT2) and determined its structure at 2.3 Å resolution. Structural alignment of TaACT2 and HvACT1-1 showed that the architectures of the substrate binding pockets were well conserved. However, the structure of a loop located at the entrance to acyl-acceptor binding site may be more flexible in TaACT2, which could be responsible for the lower affinity of TaACT2 to agmatine. Mutations of HvACT1-1 at Glu372 and Asp374 within one of the clade-IV specific motifs facing the deduced acyl-acceptor binding pocket caused significant catalytic deterioration toward agmatine both in K_m and k_cat, suggesting their key roles in acyl acceptor binding by the clade-IV enzymes. This study elucidated the molecular basis of how plants accumulate defensive specialized metabolites and provided insights into developing efficient and eco-friendly agricultural methods.

Function of the HYDROXYCINNAMOYL-CoA:SHIKIMATE HYDROXYCINNAMOYL TRANSFERASE is evolutionarily conserved in embryophytes

The plant phenylpropanoid pathway generates a major class of specialized metabolites and precursors of essential extracellular polymers that initially appeared upon plant terrestrialization. Despite its evolutionary significance, little is known about the complexity and function of this major metabolic pathway in extant bryophytes, which represent the non-vascular stage of embryophyte evolution. Here, we report that the HYDROXYCINNAMOYL-CoA:SHIKIMATE HYDROXYCINNAMOYL TRANSFERASE (HCT) gene, which plays a critical function in the phenylpropanoid pathway during seed plant development, is functionally conserved in Physcomitrium patens (Physcomitrella), in the moss lineage of bryophytes. Phylogenetic analysis indicates that bona fide HCT function emerged in the progenitor of embryophytes. In vitro enzyme assays, moss phenolic pathway reconstitution in yeast and in planta gene inactivation coupled to targeted metabolic profiling, collectively indicate that P. patens HCT (PpHCT), similar to tracheophyte HCT orthologs, uses shikimate as a native acyl acceptor to produce a p-coumaroyl-5-O-shikimate intermediate. Phenotypic and metabolic analyses of loss-of-function mutants show that PpHCT is necessary for the production of caffeate derivatives, including previously reported caffeoyl-threonate esters, and for the formation of an intact cuticle. Deep conservation of HCT function in embryophytes is further suggested by the ability of HCT genes from P. patens and the liverwort Marchantia polymorpha to complement an Arabidopsis thaliana CRISPR/Cas9 hct mutant, and by the presence of phenolic esters of shikimate in representative species of the three bryophyte lineages.

A Minireview on Temperature Dependent Protein Conformational Sampling

In this minireview we discuss the role of the more subtle conformational change-protein conformational sampling and connect it to the classic relationship of protein structure and function. The theory of pre-existing functional states of protein are discussed in context of alternate protein conformational sampling. Last, we discuss how temperature, ligand binding and mutations affect the protein conformational sampling mode which is linked to the protein function regulation. The review includes several protein systems that showed temperature dependent protein conformational sampling. We also specifically included two enzyme systems, thermophilic alcohol dehydrogenase (ht-ADH) and thermolysin which we previously studied when discussing temperature dependent protein conformational sampling.

Anthocyanin 5,3'-aromatic acyltransferase from Gentiana triflora, a structural insight into biosynthesis of a blue anthocyanin

The acylation of anthocyanins contributes to their structural diversity. Aromatic acylation is responsible for the blue color of anthocyanins and certain flowers. Aromatic acyltransferase from Gentiana triflora Pall. (Gentianaceae) (Gt5,3'AT) catalyzes the acylation of glucosyl moieties at the 5 and 3' positions of anthocyanins. Anthocyanin acyltransferase transfers an acyl group to a single position, such that Gt5,3'AT possesses a unique enzymatic activity. Structural investigation of this aromatic acyl group transfer is fundamental to understand the molecular mechanism of the acylation of double positions. In this study, structural analyses of Gt5,3'AT were conducted to identify the underlying mechanism. The crystal structure indicated that Gt5,3'AT shares structural similarities with other BAHD family enzymes, consisting of N and C terminal lobes. Structural comparison revealed that acyl group preference (aromatic or aliphatic) for the enzymes was determined by four amino acid positions, which are well conserved in aromatic and aliphatic CoA-binding acyltransferases. Although a complex structure with anthocyanins was not obtained, the binding of delphinidin 3,5,3'-triglucoside to Gt5,3'AT was investigated by evaluating the molecular dynamics. The simulation indicated that acyl transfer by Gt5,3'AT preferentially occurs at the 5-position rather than at the 3'-position, with interacting amino acids that are mainly located in the C-terminal lobe. Subsequent assays of chimeric enzymes (exchange of the N-terminal lobe and the C-terminal lobe between Gt5,3'AT and lisianthus anthocyanin 5AT) demonstrated that acyl transfer selectivity may be caused by the C-terminal lobe.

In vitro inhibition of shikimate hydroxycinnamoyltransferase by acibenzolar acid, the first metabolite of the plant defence inducer acibenzolar-S-methyl

Acibenzolar acid, the first metabolite formed in planta from the defence inducer acibenzolar-S-methyl (ASM), has been shown to be an inhibitor of the enzyme shikimate hydroxycinnamoyltransferase (HST), extracted from grapevine or tobacco cell suspension cultures. Using a purified recombinant Arabidopsis thaliana HST, the inhibition was found to be competitive, acibenzolar acid binding reversibly to the shikimate binding site of the HST:p-coumaroyl-CoA complex, with a Ki value of 250 μM. The other hydroxycinnamoyltransferases tested in the course of this study, using either hydroxypalmitic acid, putrescine, tyramine, or quinic acid as acyl acceptors were not, or only slightly, inhibited by acibenzolar acid. To understand the specificity of the interaction of acibenzolar acid with HST, we analyzed the structure-activity relationship of a series of benzoic or acibenzolar acid analogues, tested either as AtHST substrates or as inhibitors. This analysis confirmed previously published data on the substrate flexibility of HST and demonstrated that both the carboxyl group and the thiadiazole moiety of acibenzolar acid are playing an important role in the interaction with the shikimate binding site. Acibenzolar acid, which cannot form an ester bond with p-coumaric acid, was however a less potent inhibitor than protocatechuic or 3-hydroxybenzoic acids, which are used as acyl acceptors by HST. Our results show that the interaction of acibenzolar acid with HST, which is probably directly linked to the substrate promiscuity of HST, is unlikely to play a direct role in the defence-inducing properties of ASM in plants.

Versatility in acyltransferase activity completes chicoric acid biosynthesis in purple coneflower

Purple coneflower (Echinacea purpurea (L.) Moench) is a popular native North American herbal plant. Its major bioactive compound, chicoric acid, is reported to have various potential physiological functions, but little is known about its biosynthesis. Here, taking an activity-guided approach, we identify two cytosolic BAHD acyltransferases that form two intermediates, caftaric acid and chlorogenic acid. Surprisingly, a unique serine carboxypeptidase-like acyltransferase uses chlorogenic acid as its acyl donor and caftaric acid as its acyl acceptor to produce chicoric acid in vacuoles, which has evolved its acyl donor specificity from the better-known 1-O-β-D-glucose esters typical for this specific type of acyltransferase to chlorogenic acid. This unusual pathway seems unique to Echinacea species suggesting convergent evolution of chicoric acid biosynthesis. Using these identified acyltransferases, we have reconstituted chicoric acid biosynthesis in tobacco. Our results emphasize the flexibility of acyltransferases and their roles in the evolution of specialized metabolism in plants.

Characterization and functional analysis of the Hydroxycinnamoyl-CoA: shikimate hydroxycinnamoyl transferase (HCT) gene family in poplar

Hydroxycinnamoyl-CoA: shikimate hydroxycinnamoyl transferase (HCT) divides the mass flux to H, G and S units in monolignol biosynthesis and affects lignin content. Ten HCT homologs were identified in the Populus trichocarpa (Torr. & Gray) genome. Both genome duplication and tandem duplication resulted in the expansion of HCT orthologs in Populus. Comprehensive analysis including motif analysis, phylogenetic analysis, expression profiles and co-expression analysis revealed the divergence and putative function of these candidate PoptrHCTs. PoptrHCT1 and 2 were identified as likely involved in lignin biosynthesis. PoptrHCT9 and 10- are likely to be involved in plant development and the response to cold stress. Similar functional divergence was also identified in Populus tomentosa Carr. Enzymatic assay of PtoHCT1 showed that PtoHCT1 was able to synthesize caffeoyl shikimate using caffeoyl-CoA and shikimic acid as substrates.

Targeting hydroxycinnamoyl CoA: shikimate hydroxycinnamoyl transferase for lignin modification in Brachypodium distachyon

BACKGROUND

Hydroxycinnamoyl CoA: shikimate hydroxycinnamoyl transferase (HCT) is a central enzyme of the so-called "esters" pathway to monolignols. As originally envisioned, HCT functions twice in this pathway, to form coumaroyl shikimate and then, in the "reverse" direction, to convert caffeoyl shikimate to caffeoyl CoA. The discovery of a caffeoyl shikimate esterase (CSE) that forms caffeic acid directly from caffeoyl shikimate calls into question the need for the reverse HCT reaction in lignin biosynthesis. Loss of function of HCT gives severe growth phenotypes in several dicot plants, but less so in some monocots, questioning whether this enzyme, and therefore the shikimate shunt, plays the same role in both monocots and dicots. The model grass Brachypodium distachyon has two HCT genes, but lacks a classical CSE gene. This study was therefore conducted to evaluate the utility of HCT as a target for lignin modification in a species with an "incomplete" shikimate shunt.

RESULTS

The kinetic properties of recombinant B. distachyon HCTs were compared with those from Arabidopsis thaliana, Medicago truncatula, and Panicum virgatum (switchgrass) for both the forward and reverse reactions. Along with two M. truncatula HCTs, B. distachyon HCT2 had the least kinetically unfavorable reverse HCT reaction, and this enzyme is induced when HCT1 is down-regulated. Down regulation of B. distachyon HCT1, or co-down-regulation of HCT1 and HCT2, by RNA interference led to reduced lignin levels, with only modest changes in lignin composition and molecular weight.

CONCLUSIONS

Down-regulation of HCT1, or co-down-regulation of both HCT genes, in B. distachyon results in less extensive changes in lignin content/composition and cell wall structure than observed following HCT down-regulation in dicots, with little negative impact on biomass yield. Nevertheless, HCT down-regulation leads to significant improvements in biomass saccharification efficiency, making this gene a preferred target for biotechnological improvement of grasses for bioprocessing.

The evolution and engineering of enzyme activity through tuning conformational landscapes

Proteins are dynamic molecules whose structures consist of an ensemble of conformational states. Dynamics contribute to protein function and a link to protein evolution has begun to emerge. This increased appreciation for the evolutionary impact of conformational sampling has grown from our developing structural biology capabilities and the exploration of directed evolution approaches, which have allowed evolutionary trajectories to be mapped. Recent studies have provided empirical examples of how proteins can evolve via conformational landscape alterations. Moreover, minor conformational substates have been shown to be involved in the emergence of new enzyme functions as they can become enriched through evolution. The role of remote mutations in stabilizing new active site geometries has also granted insight into the molecular basis underpinning poorly understood epistatic effects that guide protein evolution. Finally, we discuss how the growth of our understanding of remote mutations is beginning to refine our approach to engineering enzymes.

Arabidopsis ECERIFERUM2-LIKEs Are Mediators of Condensing Enzyme Function

Condensing enzymes catalyze the committed reaction of fatty acid elongation and determine the chain length of fatty acids accepted and produced by the elongation complex. While necessary for the elongation of very-long-chain fatty acids (VLCFAs), identified plant condensing enzymes cannot efficiently produce VLCFAs longer than 28 carbons, which are precursors for the most abundant cuticular waxes of most plant species that have been surveyed. The eceriferum2 (cer2) mutant of Arabidopsis thaliana has a severe wax-deficient phenotype and specifically lacks waxes longer than 28 carbons, but the CER2 protein does not share sequence similarity with condensing enzymes. Instead, CER2 is homologous to BAHD acyltransferases. Heterologous expression in yeast previously demonstrated that CER2, and a small clade of BAHD acyltransferases with high sequence identity to CER2, can extend the chain-length specificity of the condensing enzyme CER6. This biochemical function is distinct from that of the broader BAHD acyltransferase family. The product specificity and physiological functions of individual CER2-LIKE proteins are unique. Here, we demonstrate that CER2 physically interacts with the fatty acid elongase. We cloned chimeric CER2-LIKE proteins and expressed these in yeast cells to identify the features that define the substrate specificities of CER2-LIKEs. We generated homology-based structural models to compare CER2-LIKEs and BAHD acyltransferases. In addition, based on the current phylogenetic analysis of the CER2-LIKE clade, we describe two further Arabidopsis CER2-LIKE genes, CER2-LIKE3 and CER2-LIKE4. We used yeast expression and mutant analysis to characterize these genes. Collectively, these results expand our knowledge of the functions of CER2-LIKEs, the BAHD acyltransferase family and cuticular wax metabolism.

Genome-wide analysis of general phenylpropanoid and monolignol-specific metabolism genes in sugarcane

Lignin is the main component of secondary cell walls and is essential for plant development and defense. However, lignin is recognized as a major recalcitrant factor for efficiency of industrial biomass processing. Genes involved in general phenylpropanoid and monolignol-specific metabolism in sugarcane have been previously analyzed at the transcriptomic level. Nevertheless, the number of genes identified in this species is still very low. The recently released sugarcane genome sequence has allowed the genome-wide characterization of the 11 gene families involved in the monolignol biosynthesis branch of the phenylpropanoid pathway. After an exhaustive analysis of sugarcane genomes, 438 haplotypes derived from 175 candidate genes from Saccharum spontaneum and 144 from Saccharum hybrid R570 were identified as associated with this biosynthetic route. The phylogenetic analyses, combined with the search for protein conserved residues involved in the catalytic activity of the encoded enzymes, were employed to identify the family members potentially involved in developmental lignification. Accordingly, 15 candidates were identified as bona fide lignin biosynthesis genes: PTAL1, PAL2, C4H4, 4CL1, HCT1, HCT2, C3'H1, C3'H2, CCoAOMT1, COMT1, F5H1, CCR1, CCR2, CAD2, and CAD7. For this core set of lignin biosynthetic genes, we searched for the chromosomal location, the gene expression pattern, the promoter cis-acting elements, and microRNA targets. Altogether, our results present a comprehensive characterization of sugarcane general phenylpropanoid and monolignol-specific genes, providing the basis for further functional studies focusing on lignin biosynthesis manipulation and biotechnological strategies to improve sugarcane biomass utilization.

NO GAMETOPHORES 2 Is a Novel Regulator of the 2D to 3D Growth Transition in the Moss Physcomitrella patens

The colonization of land by plants was one of the most transformative events in the history of life on Earth. The transition from water, which coincided with and was likely facilitated by the evolution of three-dimensional (3D) growth, enabled the generation of morphological diversity on land. In many plants, the transition from two-dimensional (2D) to 3D growth occurs during embryo development. However, in the early divergent moss Physcomitrella patens, 3D growth is preceded by an extended filamentous phase that can be maintained indefinitely. Here, we describe the identification of the cytokinin-responsive NO GAMETOPHORES 2 (PpNOG2) gene, which encodes a shikimate o-hydroxycinnamoyltransferase. In mutants lacking PpNOG2 function, transcript levels of CLAVATA and SCARECROW genes are significantly reduced, excessive gametophore initial cells are produced, and buds undergo premature developmental arrest. Mutants also exhibit misregulation of auxin-responsive genes. Our results suggest that PpNOG2 functions in the ascorbic acid pathway leading to cuticle formation and that NOG2-related genes were co-opted into the lignin biosynthesis pathway after the divergence of bryophytes and vascular plants. We present a revised model of 3D growth in which PpNOG2 comprises part of a feedback mechanism that is required for the modulation of gametophore initial cell frequency. We also propose that the 2D to 3D growth transition in P. patens is underpinned by complex auxin-cytokinin crosstalk that is regulated, at least in part, by changes in flavonoid metabolism.

Crystal structure of barley agmatine coumaroyltransferase, an N-acyltransferase from the BAHD superfamily

The enzymes of the BAHD superfamily, a large group of acyl-CoA-dependent acyltransferases in plants, are involved in the biosynthesis of diverse secondary metabolites. While the structures of several O-acyltransferases from the BAHD superfamily, such as hydroxycinnamoyl-CoA shikimate hydroxycinnamoyl transferase, have been elucidated, no structural information on N-acyltransferases is available. Hordeum vulgare agmatine coumaroyltransferase (HvACT) is an N-acyltransferase from the BAHD superfamily and is one of the most important enzymes in the secondary metabolism of barley. Here, an apo-form structure of HvACT is reported as the first structure of an N-acyltransferase from the BAHD superfamily. HvACT crystals diffracted to 1.8 Å resolution and belonged to the monoclinic space group P2₁, with unit-cell parameters a = 57.6, b = 59.5, c = 73.6 Å, α = 90, β = 91.3 , γ = 90°. Like other known BAHD superfamily structures, HvACT contains two domains that adopt a two-layer αβ-sandwich architecture and a solvent-exposed channel that penetrates the enzyme core.

Functional characterization and catalytic activity improvement of BAHD acyltransferase from Celastrus angulatus Maxim

A BAHD terpene alcohol acyltransferase, CaAT20, was identified from Celastrus angulatus Maxim, expressed in E. coli and functionally characterized. S405A mutant of CaAT20 increased the enzyme activity. Acylation is a diversely physiological process in the biosynthesis of plant secondary metabolites. Plant BAHD acyltransferases play an important role in the modification of volatile esters with biological activities. In this research, a BAHD acyltransferase (CaAT20) was identified from Celastrus angulatus Maxim and the function of this enzyme was characterized. CaAT20 could convert geraniol to geranyl esters by using benzoyl-CoA and acetyl-CoA as the acyl donors respectively. Furthermore, the catalytic activity of CaAT20 for benzoyl-CoA was higher than that of acetyl-CoA. Site-directed mutation of CaAT20 was carried out based on the results of molecular simulation. In vitro site-directed mutant S405A of CaAT20 increased the volume of binding cavity so as to facilitate the entry of geraniol, indicating a more efficient acylation for geraniol and benzoyl-CoA. Our research provides new insight for the catalytic functions of CaAT20.

The physical basis and practical consequences of biological promiscuity

Proteins interact with metabolites, nucleic acids, and other proteins to orchestrate the myriad catalytic, structural and regulatory functions that support life from the simplest microbes to the most complex multicellular organisms. These molecular interactions are often exquisitely specific, but never perfectly so. Adventitious "promiscuous" interactions are ubiquitous due to the thousands of macromolecules and small molecules crowded together in cells. Such interactions may perturb protein function at the molecular level, but as long as they do not compromise organismal fitness, they will not be removed by natural selection. Although promiscuous interactions are physiologically irrelevant, they are important because they can provide a vast reservoir of potential functions that can provide the starting point for evolution of new functions, both in nature and in the laboratory.

Identification of genes from the general phenylpropanoid and monolignol-specific metabolism in two sugarcane lignin-contrasting genotypes

The phenylpropanoid pathway is an important route of secondary metabolism involved in the synthesis of different phenolic compounds such as phenylpropenes, anthocyanins, stilbenoids, flavonoids, and monolignols. The flux toward monolignol biosynthesis through the phenylpropanoid pathway is controlled by specific genes from at least ten families. Lignin polymer is one of the major components of the plant cell wall and is mainly responsible for recalcitrance to saccharification in ethanol production from lignocellulosic biomass. Here, we identified and characterized sugarcane candidate genes from the general phenylpropanoid and monolignol-specific metabolism through a search of the sugarcane EST databases, phylogenetic analysis, a search for conserved amino acid residues important for enzymatic function, and analysis of expression patterns during culm development in two lignin-contrasting genotypes. Of these genes, 15 were cloned and, when available, their loci were identified using the recently released sugarcane genomes from Saccharum hybrid R570 and Saccharum spontaneum cultivars. Our analysis points out that ShPAL1, ShPAL2, ShC4H4, Sh4CL1, ShHCT1, ShC3H1, ShC3H2, ShCCoAOMT1, ShCOMT1, ShF5H1, ShCCR1, ShCAD2, and ShCAD7 are strong candidates to be bona fide lignin biosynthesis genes. Together, the results provide information about the candidate genes involved in monolignol biosynthesis in sugarcane and may provide useful information for further molecular genetic studies in sugarcane.

Structural and Biochemical Insights Into Two BAHD Acyltransferases (AtSHT and AtSDT) Involved in Phenolamide Biosynthesis

Phenolamides represent one of the largest classes of plant-specialized secondary metabolites and function in diverse physiological processes, including defense responses and development. The biosynthesis of phenolamides requires the BAHD-family acyltransferases, which transfer acyl-groups from different acyl-donors specifically to amines, the acyl-group acceptors. However, the mechanisms of substrate specificity and multisite-acylation of the BAHD-family acyltransferases remain poorly understood. In this study, we provide a structural and biochemical analysis of AtSHT and AtSDT, two representative BAHD-family members that catalyze the multisite acylation of spermidine but show different product profiles. By determining the structures of AtSHT and AtSDT and using structure-based mutagenesis, we identified the residues important for substrate recognition in AtSHT and AtSDT and hypothesized that the acyl acceptor spermidine might adopt a free-rotating conformation in AtSHT, which can undergo mono-, di-, or tri-acylation; while the spermidine molecule in AtSDT might adopt a linear conformation, which only allows mono- or di-acylation to take place. In addition, through sequence similarity network (SSN) and structural modeling analysis, we successfully predicted and verified the functions of two uncharacterized Arabidopsis BAHD acyltransferases, OAO95042.1 and NP_190301.2, which use putrescine as the main acyl-acceptor. Our work provides not only an excellent starting point for understanding multisite acylation in BAHD-family enzymes, but also a feasible methodology for predicting possible acyl acceptor specificity of uncharacterized BAHD-family acyltransferases.

PBS3 and EPS1 Complete Salicylic Acid Biosynthesis from Isochorismate in Arabidopsis

Salicylic acid (SA) is an important phytohormone mediating both local and systemic defense responses in plants. Despite over half a century of research, how plants biosynthesize SA remains unresolved. In Arabidopsis, a major part of SA is derived from isochorismate, a key intermediate produced by the isochorismate synthase, which is reminiscent of SA biosynthesis in bacteria. Whereas bacteria employ an isochorismate pyruvate lyase (IPL) that catalyzes the turnover of isochorismate to pyruvate and SA, plants do not contain an IPL ortholog and generate SA from isochorismate through an unknown mechanism. Combining genetic and biochemical approaches, we delineated the SA biosynthetic pathway downstream of isochorismate in Arabidopsis. We found that PBS3, a GH3 acyl adenylase-family enzyme important for SA accumulation, catalyzes ATP- and Mg²⁺-dependent conjugation of L-glutamate primarily to the 8-carboxyl of isochorismate and yields the key SA biosynthetic intermediate, isochorismoyl-glutamate A. Moreover, we discovered that EPS1, a BAHD acyltransferase-family protein with a previously implicated role in SA accumulation upon pathogen attack, harbors a noncanonical active site and an unprecedented isochorismoyl-glutamate A pyruvoyl-glutamate lyase activity that produces SA from the isochorismoyl-glutamate A substrate. Together, PBS3 and EPS1 form a two-step metabolic pathway to produce SA from isochorismate in Arabidopsis, which is distinct from how SA is biosynthesized in bacteria. This study closes a major knowledge gap in plant SA metabolism and would help develop new strategies for engineering disease resistance in crop plants.

Independent evolution of rosmarinic acid biosynthesis in two sister families under the Lamiids clade of flowering plants

As a means to maintain their sessile lifestyle amid challenging environments, plants produce an enormous diversity of compounds as chemical defenses against biotic and abiotic insults. The underpinning metabolic pathways that support the biosynthesis of these specialized chemicals in divergent plant species provide a rich arena for understanding the molecular evolution of complex metabolic traits. Rosmarinic acid (RA) is a phenolic natural product first discovered in plants of the mint family (Lamiaceae) and is recognized for its wide range of medicinal properties and potential applications in human dietary and medical interventions. Interestingly, the RA chemotype is present sporadically in multiple taxa of flowering plants as well as some hornworts and ferns, prompting the question whether its biosynthesis arose independently across different lineages. Here we report the elucidation of the RA biosynthetic pathway in Phacelia campanularia (desert bells). This species represents the borage family (Boraginaceae), an RA-producing family closely related to the Lamiaceae within the Lamiids clade. Using a multi-omics approach in combination with functional characterization of candidate genes both in vitro and in vivo, we found that RA biosynthesis in P. campanularia involves specific activities of a BAHD acyltransferase and two cytochrome P450 hydroxylases. Further phylogenetic and comparative structure-function analyses of the P. campanularia RA biosynthetic enzymes clearly indicate that RA biosynthesis has evolved independently at least twice in the Lamiids, an exemplary case of chemotypic convergence through disparate evolutionary trajectories.

COSY catalyses trans-cis isomerization and lactonization in the biosynthesis of coumarins

Coumarins, also known as 1,2-benzopyrones, comprise a large class of secondary metabolites that are ubiquitously found throughout the plant kingdom. In many plant species, coumarins are particularly important for iron acquisition and plant defence. Here, we show that COUMARIN SYNTHASE (COSY) is a key enzyme in the biosynthesis of coumarins. Arabidopsis thaliana cosy mutants have strongly reduced levels of coumarin and accumulate o-hydroxyphenylpropanoids instead. Accordingly, cosy mutants have reduced iron content and show growth defects when grown under conditions in which there is a limited availability of iron. Recombinant COSY is able to produce umbelliferone, esculetin and scopoletin from their respective o-hydroxycinnamoyl-CoA thioesters by two reaction steps-a trans-cis isomerization followed by a lactonization. This conversion happens partially spontaneously and is catalysed by light, which explains why the need for an enzyme for this conversion has been overlooked. The combined results show that COSY has an essential function in the biosynthesis of coumarins in organs that are shielded from light, such as roots. These findings provide routes to improving coumarin production in crops or by microbial fermentation.

Structural and dynamic basis of substrate permissiveness in hydroxycinnamoyltransferase (HCT)

Substrate permissiveness has long been regarded as the raw materials for the evolution of new enzymatic functions. In land plants, hydroxycinnamoyltransferase (HCT) is an essential enzyme of the phenylpropanoid metabolism. Although essential enzymes are normally associated with high substrate specificity, HCT can utilize a variety of non-native substrates. To examine the structural and dynamic basis of substrate permissiveness in this enzyme, we report the crystal structure of HCT from Selaginella moellendorffii and molecular dynamics (MD) simulations performed on five orthologous HCTs from several major lineages of land plants. Through altogether 17-μs MD simulations, we demonstrate the prevalent swing motion of an arginine handle on a submicrosecond timescale across all five HCTs, which plays a key role in native substrate recognition by these intrinsically promiscuous enzymes. Our simulations further reveal how a non-native substrate of HCT engages a binding site different from that of the native substrate and diffuses to reach the catalytic center and its co-substrate. By numerically solving the Smoluchowski equation, we show that the presence of such an alternative binding site, even when it is distant from the catalytic center, always increases the reaction rate of a given substrate. However, this increase is only significant for enzyme-substrate reactions heavily influenced by diffusion. In these cases, binding non-native substrates ‘off-center’ provides an effective rationale to develop substrate permissiveness while maintaining the native functions of promiscuous enzymes.

Examples abound of enzymes that can process substrates other than their native ones. However, the structural and dynamic basis of this promiscuity remains to be fully understood. In this work, we examine HCT, an intrinsically promiscuous acyltransferase with conserved function in all land plants. We uncover the sub-microsecond swing motion of a key arginine residue facilitating the recognition of both native and non-native substrates of HCT. We also quantify the impact of an off-center binding site on the non-native reaction rate. Although our calculations were inspired by HCT, the results apply in general, i.e., for enzymes heavily influenced by diffusion, binding non-native substrates ‘off-center’, even with rather weak affinity, can accelerate non-native reactions to appreciable levels.

Isolation and functional characterization of hydroxycinnamoyltransferases from the liverworts Plagiochasma appendiculatum and Marchantia paleacea

Hydroxycinnamoyl-CoA shikimate/quinate hydroxycinnamoyl transferase (HCT, EC: 2.3.1.133) is a key metabolic entry point for the synthesis of monolignols in vascular plants; however, little is known about HCT in liverworts. Here, the isolation and characterization of HCTs encoded by the two liverwort species, Plagiochasma appendiculatum and Marchantia paleacea, are described. The sequences of the two enzymes harbor features typical of BAHD family members, except for the presence of a stretch of >100 residues that are not represented in higher plant HCTs. When truncated versions of both genes, which were constructed to clarify the significance of these extra residues, were investigated, it became apparent that the full-length and the truncated gene products shared similar catalytic activity and recognized the same substrates in vitro. They also functioned equivalently in vivo either when transiently expressed in tobacco to cause a higher total production of CGA (5-CQA) and 4-CQA or stably expressed in liverworts to accumulate the lignin-like contents. A structural model of MpHCT suggests that its active site bind to its substrate similar to that of Arabidopsis thaliana HCT. While truncated forms of HCT were deposited in the nucleocytoplasm, the full-length versions occurred exclusively in the cytoplasm. The conclusion is that liverworts produce bona fide HCTs that represent a point of departure in studying the evolution of lignin synthesis in plants.

Structural and Biochemical Studies of a Biocatalyst for the Enzymatic Production of Wax Esters

Wax esters are high-value products whose enzymatic synthesis is of increasing biotechnological interest. The fabrication of cell factories that mass-produce wax esters may provide a facile route towards a sustainable, and environment-friendly approach to a large-scale process for this commodity chemical. An expedient route for wax-ester biocatalysis may be facilitated by the action of enzymes termed wax ester synthases/diacylglycerol acyltransferases (WS/DGAT), which produce wax esters using fatty acids and alcohols as a precursor. In this work, we report the structure for a member of the WS/DGAT superfamily. The structural data in conjunction with bioinformatics and mutational analyses allowed us to identify the substrate binding pockets, and residues that may be important for catalysis. Using this information as a guide, we generated a mutant with preference towards shorter acyl-substrates. This study demonstrates the efficacy of a structure-guided engineering effort towards a WS/DGAT variant with preference towards wax esters of desired lengths.

Promiscuity, impersonation and accommodation: evolution of plant specialized metabolism

Specialized metabolic enzymes and metabolite diversity evolve through a variety of mechanisms including promiscuity, changes in substrate specificity, modifications of gene expression and gene duplication. For example, gene duplication and substrate binding site changes led to the evolution of the glucosinolate biosynthetic enzyme, AtIPMDH1, from a Leu biosynthetic enzyme. BAHD acyltransferases illustrate how enzymatic promiscuity leads to metabolite diversity. The examples 4-coumarate:CoA ligase and aromatic acid methyltransferases illustrate how promiscuity can potentiate the evolution of these specialized metabolic enzymes.

The relationship between folding and activity in UreG, an intrinsically disordered enzyme

A growing body of literature on intrinsically disordered proteins (IDPs) led scientists to rethink the structure-function paradigm of protein folding. Enzymes are often considered an exception to the rule of intrinsic disorder (ID), believed to require a unique structure for catalysis. However, recent studies revealed the presence of disorder in several functional native enzymes. In the present work, we address the importance of dynamics for catalysis, by investigating the relationship between folding and activity in Sporosarcina pasteurii UreG (SpUreG), a P-loop GTPase and the first discovered native ID enzyme, involved in the maturation of the nickel-containing urease. The effect of denaturants and osmolytes on protein structure and activity was analyzed using circular dichroism (CD), Site-Directed Spin Labeling (SDSL) coupled to EPR spectroscopy, and enzymatic assays. Our data show that SpUreG needs a "flexibility window" to be catalytically competent, with both too low and too high mobility being detrimental for its activity.

A phenol-enriched cuticle is ancestral to lignin evolution in land plants

Lignin, one of the most abundant biopolymers on Earth, derives from the plant phenolic metabolism. It appeared upon terrestrialization and is thought critical for plant colonization of land. Early diverging land plants do not form lignin, but already have elements of its biosynthetic machinery. Here we delete in a moss the P450 oxygenase that defines the entry point in angiosperm lignin metabolism, and find that its pre-lignin pathway is essential for development. This pathway does not involve biochemical regulation via shikimate coupling, but instead is coupled with ascorbate catabolism, and controls the synthesis of the moss cuticle, which prevents desiccation and organ fusion. These cuticles share common features with lignin, cutin and suberin, and may represent the extant representative of a common ancestor. Our results demonstrate a critical role for the ancestral phenolic metabolism in moss erect growth and cuticle permeability, consistent with importance in plant adaptation to terrestrial conditions.

The phenolic polymer lignin is thought to have contributed to adaptation of early land plants to terrestrial environments. Here Renault et al. show that moss, which does not produce lignin, contains an ancestral phenolic metabolism pathway that produces a phenol-enriched cuticle and prevents desiccation.