REVIEW: Epistasis and dominance in the emergence of catalytic function as exemplified by the evolution of plant terpene synthases

Decoding Catalysis by Terpene Synthases

The review by Christianson, published in 2017 on the twentieth anniversary of the emergence of the field, summarizes the foundational discoveries and key advances in terpene synthase/cyclase (TS) biocatalysis (Christianson, D. W. Chem Rev2017, 117 (17), 11570-11648. DOI: 10.1021/acs.chemrev.7b00287). Here, we review the TS literature published since then, bringing the field up to date and looking forward to what could be the near future of TS rational design. Many revealing discoveries have been made in recent years, building on the knowledge and fundamental principles uncovered during those initial two decades of study. We use these to explore TS reaction chemistry and see how a combined experimental and computational approach helps to decipher the complexities of TS catalysis. Revealed are a suite of catalytic motifs which control product outcome in TSs, some obvious, some more subtle. We examine each in detail, using the most recent papers and insights to illustrate how exactly this fascinating class of enzymes takes a single acyclic substrate and turns it into the many thousands of complex terpenoids found in Nature. We then explore some of the recent strategies for TS engineering, including machine learning and other data-driven approaches. From this, rational and predictive engineering of TSs, "designer terpene synthases", will begin to emerge as a realistic goal.

Non-enzymatic formation of isoprene and 2-methyl-3-buten-2-ol (2-MBO) by manganese

It has been suggested that isoprene synthesis by isoprene synthase (IspS) proceeds via a substrate-assisted mechanism. The authors observed a non-enzymatic isoprene formation by Mn²⁺, which represents the basis of IspS enzyme reaction. Because IspS and many other terpene synthases require Mn²⁺ metal ions as cofactor, this study characterized the formation reaction for the first time. Metal ions including Mn²⁺ non-enzymatically produced both isoprene and 2-methyl-3-buten-2-ol (2-MBO) from dimethylallyl pyrophosphate (DMADP). Isoprene formation was most enhanced by Fe²⁺ and, to a lesser extent, by Mn²⁺ or Cu²⁺. Ni²⁺, Co²⁺, Mg²⁺, and Ba²⁺ exhibited a low activity to generate both isoprene and 2-MBO. The proportion of isoprene and 2-MBO varied with the Mn²⁺ concentration: isoprene predominated over 2-MBO at a higher Mn²⁺ concentration. Similarly, isoprene formation by Mn²⁺ increased exponentially as temperature increased with predominance of isoprene over 2-MBO at higher temperature. Both isoprene and 2-MBO formation was enhanced by acidic and neutral pH compared to alkaline conditions. Molecular dynamic simulation of DMADP suggested that the formation reaction is initiated by deprotonation of hydrogen on allyl terminal carbon by phosphate oxygen and generates carbocation and allyl anion intermediates. This is followed by quenching to produce isoprene or by hydroxyl addition to form 2-MBO. Thus, this study provided an insight into reaction mechanism of isoprene and 2-MBO biosynthesis and highlighted some parts of isoprene emission from terrestrial plants, which could be formed by non-enzymatic mechanism.

Evolutionary History of Plant Metabolism

Tremendous chemical diversity is the hallmark of plants and is supported by highly complex biochemical machinery. Plant metabolic enzymes originated and were transferred from eukaryotic and prokaryotic ancestors and further diversified by the unprecedented rates of gene duplication and functionalization experienced in land plants. Unlike microbes, which have frequent horizontal gene transfer events and multiple inputs of energy and organic carbon, land plants predominantly rely on organic carbon generated from CO₂ and have experienced very few, if any, gene transfers during their recent evolutionary history. As such, plant metabolic networks have evolved in a stepwise manner and on existing networks under various evolutionary constraints. This review aims to take a broader view of plant metabolic evolution and lay a framework to further explore evolutionary mechanisms of the complex metabolic network. Understanding the underlying metabolic and genetic constraints is also an empirical prerequisite for rational engineering and redesigning of plant metabolic pathways.

Epistasis and intramolecular networks in protein evolution

Proteins are molecular machines composed of complex, highly connected amino acid networks. Their functional optimization requires the reorganization of these intramolecular networks by evolution. In this review, we discuss the mechanisms by which epistasis, that is, the dependence of the effect of a mutation on the genetic background, rewires intramolecular interactions to alter protein function. Deciphering the biophysical basis of epistasis is crucial to our understanding of evolutionary dynamics and the elucidation of sequence-structure-function relationships. We featured recent studies that provide insights into the molecular mechanisms giving rise to epistasis, particularly at the structural level. These studies illustrate the convoluted and fascinating nature of the intramolecular networks co-opted by epistasis during the evolution of protein function.

Enzyme Evolution: An Epistatic Ratchet versus a Smooth Reversible Transition

Evolutionary trajectories are deemed largely irreversible. In a newly diverged protein, reversion of mutations that led to the functional switch typically results in loss of both the new and the ancestral functions. Nonetheless, evolutionary transitions where reversions are viable have also been described. The structural and mechanistic causes of reversion compatibility versus incompatibility therefore remain unclear. We examined two laboratory evolution trajectories of mammalian paraoxonase-1, a lactonase with promiscuous organophosphate hydrolase (OPH) activity. Both trajectories began with the same active-site mutant, His115Trp, which lost the native lactonase activity and acquired higher OPH activity. A neo-functionalization trajectory amplified the promiscuous OPH activity, whereas the re-functionalization trajectory restored the native activity, thus generating a new lactonase that lacks His115. The His115 revertants of these trajectories indicated opposite trends. Revertants of the neo-functionalization trajectory lost both the evolved OPH and the original lactonase activity. Revertants of the trajectory that restored the original lactonase function were, however, fully active. Crystal structures and molecular simulations show that in the newly diverged OPH, the reverted His115 and other catalytic residues are displaced, thus causing loss of both the original and the new activity. In contrast, in the re-functionalization trajectory, reversion compatibility of the original lactonase activity derives from mechanistic versatility whereby multiple residues can fulfill the same task. This versatility enables unique sequence-reversible compositions that are inaccessible when the active site was repurposed toward a new function.

Complex Patterns of Cannabinoid Alkyl Side-Chain Inheritance in Cannabis

The cannabinoid alkyl side-chain represents an important pharmacophore, where genetic targeting of alkyl homologs has the potential to provide enhanced forms of Cannabis for biopharmaceutical manufacture. Delta(9)-tetrahydrocannabinolic acid (THCA) and cannabidiolic acid (CBDA) synthase genes govern dicyclic (CBDA) and tricyclic (THCA) cannabinoid composition. However, the inheritance of alkyl side-chain length has not been resolved, and few studies have investigated the contributions and interactions between cannabinoid synthesis pathway loci. To examine the inheritance of chemical phenotype (chemotype), THCAS and CBDAS genotypes were scored and alkyl cannabinoid segregation analysed in 210 F₂ progeny derived from a cross between two Cannabis chemotypes divergent for alkyl and cyclic cannabinoids. Inheritance patterns of F₂ progeny were non-Gaussian and deviated from Mendelian expectations. However, discrete alkyl cannabinoid segregation patterns consistent with digenic as well as epistatic modes of inheritance were observed among F₂ THCAS and CBDAS genotypes. These results suggest linkage between cannabinoid pathway loci and highlight the need for further detailed characterisation of cannabinoid inheritance to facilitate metabolic engineering of chemically elite germplasm.

Understanding enzyme function evolution from a computational perspective

In this review, we will explore recent computational approaches to understand enzyme evolution from the perspective of protein structure, dynamics and promiscuity. We will present quantitative methods to measure the size of evolutionary steps within a structural domain, allowing the correlation between change in substrate and domain structure to be assessed, and giving insights into the evolvability of different domains in terms of the number, types and sizes of evolutionary steps observed. These approaches will help to understand the evolution of new catalytic and non-catalytic functionality in response to environmental demands, showing potential to guide de novoenzyme design and directed evolution experiments.