An LIR motif in the Rift Valley fever virus NSs protein is critical for the interaction with LC3 family members and inhibition of autophagy.

Rift Valley fever virus (RVFV) is a viral zoonosis that causes severe disease in ruminants and humans. The nonstructural small (NSs) protein is the primary virulence factor of RVFV that suppresses the host's antiviral innate immune response. Bioinformatic analysis and AlphaFold structural modeling identified four putative LC3-interacting regions (LIR) motifs (NSs 1-4) in the RVFV NSs protein, which suggest that NSs interacts with the host LC3-family proteins. Using, isothermal titration calorimetry, X-ray crystallography, co-immunoprecipitation, and co-localization experiments, the C-terminal LIR motif (NSs4) was confirmed to interact with all six human LC3 proteins. Phenylalanine at position 261 (F261) within NSs4 was found to be critical for the interaction of NSs with LC3, retention of LC3 in the nucleus, as well as the inhibition of autophagy in RVFV infected cells. These results provide mechanistic insights into the ability of RVFV to overcome antiviral autophagy through the interaction of NSs with LC3 proteins.

Coevolutionary Information Captures Catalytic Functions and Reveals Divergent Roles of Terpene Synthase Interdomain Connections.

Inferring the historical and biophysical causes of diversity within protein families is a complex puzzle. A key to unraveling this problem is characterizing the rugged topography of sequence-function adaptive landscapes. Using biochemical data from a 29 = 512 combinatorial library of tobacco 5-epi-aristolochene synthase (TEAS) mutants engineered to make the native major product of Egyptian henbane premnaspirodiene synthase (HPS) and a complementary 512 mutant HPS library, we address the question of how product specificity is controlled. These data sets reveal that HPS is far more robust and resistant to mutations than TEAS, where most mutants are promiscuous. We also combine experimental data with a sequence Potts Hamiltonian model and direct coupling analysis to quantify mutant fitness. Our results demonstrate that the Hamiltonian captures variation in product outputs across both libraries, clusters native family members based on their substrate specificities, and exposes the divergent catalytic roles of couplings between the catalytic and noncatalytic domains of TEAS versus HPS. Specifically, we found that the role of the interdomain connectivities in specifying product output is more important in TEAS than connectivities within the catalytic domain. Despite being 75% identical, this property is not shared by HPS, where connectivities within the catalytic domain are more important for specificity. By solving the X-ray crystal structure of HPS, we assessed structural bases for their interdomain network differences. Last, we calculate the product profile Shannon entropies of the two libraries, which showcases that site-site connectivities also play divergent roles in catalytic accuracy.

Terpene synthases and pathways in animals: enzymology and structural evolution in the biosynthesis of volatile infochemicals.

Covering: up to the beginning of 2023Many animals release volatile or semi-volatile terpenes as semiochemicals in intra- and inter-specific interactions. Terpenes are important constituents of pheromones and serve as chemical defenses to ward off predators. Despite the occurrence of terpene specialized metabolites from soft corals to mammals, the biosynthetic origin of these compounds has largely remained obscure. An increasing number of animal genome and transcriptome resources is facilitating the identification of enzymes and pathways that allow animals to produce terpenes independent of their food sources or microbial endosymbionts. Substantial evidence has emerged for the presence of terpene biosynthetic pathways such as in the formation of the iridoid sex pheromone nepetalactone in aphids. In addition, terpene synthase (TPS) enzymes have been discovered that are evolutionary unrelated to canonical plant and microbial TPSs and instead resemble precursor enzymes called isoprenyl diphosphate synthases (IDSs) in central terpene metabolism. Structural modifications of substrate binding motifs in canonical IDS proteins presumably facilitated the transition to TPS function at an early state in insect evolution. Other arthropods such as mites appear to have adopted their TPS genes from microbial sources via horizontal gene transfer. A similar scenario likely occurred in soft corals, where TPS families with closer resemblance to microbial TPSs have been discovered recently. Together, these findings will spur the identification of similar or still unknown enzymes in terpene biosynthesis in other lineages of animals. They will also help develop biotechnological applications for animal derived terpenes of pharmaceutical value or advance sustainable agricultural practices in pest management.

Sparse Epistatic Patterns in the Evolution of Terpene Synthases.

We explore sequence determinants of enzyme activity and specificity in a major enzyme family of terpene synthases. Most enzymes in this family catalyze reactions that produce cyclic terpenes-complex hydrocarbons widely used by plants and insects in diverse biological processes such as defense, communication, and symbiosis. To analyze the molecular mechanisms of emergence of terpene cyclization, we have carried out in-depth examination of mutational space around (E)-ß-farnesene synthase, an Artemisia annua enzyme which catalyzes production of a linear hydrocarbon chain. Each mutant enzyme in our synthetic libraries was characterized biochemically, and the resulting reaction rate data were used as input to the Michaelis-Menten model of enzyme kinetics, in which free energies were represented as sums of one-amino-acid contributions and two-amino-acid couplings. Our model predicts measured reaction rates with high accuracy and yields free energy landscapes characterized by relatively few coupling terms. As a result, the Michaelis-Menten free energy landscapes have simple, interpretable structure and exhibit little epistasis. We have also developed biophysical fitness models based on the assumption that highly fit enzymes have evolved to maximize the output of correct products, such as cyclic products or a specific product of interest, while minimizing the output of byproducts. This approach results in nonlinear fitness landscapes that are considerably more epistatic. Overall, our experimental and computational framework provides focused characterization of evolutionary emergence of novel enzymatic functions in the context of microevolutionary exploration of sequence space around naturally occurring enzymes.

The Omics Dashboard for interactive exploration of gene-expression data.

The Omics Dashboard is a software tool for interactive exploration and analysis of gene-expression datasets. The Omics Dashboard is organized as a hierarchy of cellular systems. At the highest level of the hierarchy the Dashboard contains graphical panels depicting systems such as biosynthesis, energy metabolism, regulation and central dogma. Each of those panels contains a series of X-Y plots depicting expression levels of subsystems of that panel, e.g. subsystems within the central dogma panel include transcription, translation and protein maturation and folding. The Dashboard presents a visual read-out of the expression status of cellular systems to facilitate a rapid top-down user survey of how all cellular systems are responding to a given stimulus, and to enable the user to quickly view the responses of genes within specific systems of interest. Although the Dashboard is complementary to traditional statistical methods for analysis of gene-expression data, we show how it can detect changes in gene expression that statistical techniques may overlook. We present the capabilities of the Dashboard using two case studies: the analysis of lipid production for the marine alga Thalassiosira pseudonana, and an investigation of a shift from anaerobic to aerobic growth for the bacterium Escherichia coli.

A conserved amino acid residue critical for product and substrate specificity in plant triterpene synthases.

Triterpenes are structurally complex plant natural products with numerous medicinal applications. They are synthesized through an origami-like process that involves cyclization of the linear 30 carbon precursor 2,3-oxidosqualene into different triterpene scaffolds. Here, through a forward genetic screen in planta, we identify a conserved amino acid residue that determines product specificity in triterpene synthases from diverse plant species. Mutation of this residue results in a major change in triterpene cyclization, with production of tetracyclic rather than pentacyclic products. The mutated enzymes also use the more highly oxygenated substrate dioxidosqualene in preference to 2,3-oxidosqualene when expressed in yeast. Our discoveries provide new insights into triterpene cyclization, revealing hidden functional diversity within triterpene synthases. They further open up opportunities to engineer novel oxygenated triterpene scaffolds by manipulating the precursor supply.

Pathway collages: personalized multi-pathway diagrams.

BACKGROUND: Metabolic pathway diagrams are a classical way of visualizing a linked cascade of biochemical reactions. However, to understand some biochemical situations, viewing a single pathway is insufficient, whereas viewing the entire metabolic network results in information overload. How do we enable scientists to rapidly construct personalized multi-pathway diagrams that depict a desired collection of interacting pathways that emphasize particular pathway interactions? RESULTS: We define software for constructing personalized multi-pathway diagrams called pathway-collages using a combination of manual and automatic layouts. The user specifies a set of pathways of interest for the collage from a Pathway/Genome Database. Layouts for the individual pathways are generated by the Pathway Tools software, and are sent to a Javascript Pathway Collage application implemented using Cytoscape.js. That application allows the user to re-position pathways; define connections between pathways; change visual style parameters; and paint metabolomics, gene expression, and reaction flux data onto the collage to obtain a desired multi-pathway diagram. We demonstrate the use of pathway collages in two application areas: a metabolomics study of pathogen drug response, and an Escherichia coli metabolic model. CONCLUSIONS: Pathway collages enable facile construction of personalized multi-pathway diagrams.

Discovery of germacrene A synthases in Barnadesia spinosa: The first committed step in sesquiterpene lactone biosynthesis in the basal member of the Asteraceae.

The Andes-endemic Barnadesioideae lineage is the oldest surviving and phylogenetically basal subfamily of the Asteraceae (Compositae), a prolific group of flowering plants with world-wide distribution (∼24,000 species) marked by a rich diversity of sesquiterpene lactones (STLs). Intriguingly, there is no evidence that members of the Barnadesioideae produce STLs, specialized metabolites thought to have contributed to the adaptive success of the Asteraceae family outside South America. The biosynthesis of STLs requires the intimate expression and functional integration of germacrene A synthase (GAS) and germacrene A oxidase (GAO) to sequentially cyclize and oxidize farnesyl diphosphate into the advanced intermediate germacrene A acid leading to diverse STLs. Our previous discovery of GAO activity conserved across all major subfamilies of Asteraceae, including the phylogenetically basal lineage of Barnadesioideae, prompted further investigation of the presence of the gateway GAS in Barnadesioideae. Herein we isolated two terpene synthases (BsGAS1/BsGAS2) from the basal Barnadesia spinosa (Barnadesioideae) that displayed robust GAS activity when reconstituted in yeast and characterized in vitro. Despite the apparent lack of STLs in the Barnadesioideae, this work unambiguously confirms the presence of GAS in the basal genera of the Asteraceae. Phylogenetic analysis reveals that the two BsGASs fall into two distinct clades of the Asteraceae's GASs, and BsGAS1 clade is only retained in the evolutionary closer Cichorioideae subfamily, implicating BsGAS2 is likely the ancestral base of most GASs found in the lineages outside the Barnadesioideae. Taken together, these results show the enzymatic capacities of GAS and GAO emerged prior to the subsequent radiation of STL-producing Asteraceae subfamilies.

The Comparative Genome Dashboard.

The Comparative Genome Dashboard is a web-based software tool for interactive exploration of the similarities and differences in gene functions between organisms. It provides a high-level graphical survey of cellular functions, and enables the user to drill down to examine subsystems of interest in greater detail. At its highest level the Comparative Dashboard contains panels for cellular systems such as biosynthesis, energy metabolism, transport, and response to stimulus. Each panel contains a set of bar graphs that plot the numbers of compounds or gene products for each organism across a set of subsystems of that panel. Users can interactively drill down to focus on subsystems of interest and see grids of compounds produced or consumed by each organism, specific GO term assignments, pathway diagrams, and links to more detailed comparison pages. For example, the dashboard enables users to compare the cofactors that a set of organisms can synthesize, the metal ions that they are able to transport, their DNA damage repair capabilities, their biofilm-formation genes, and their viral response proteins. The dashboard enables users to quickly perform comprehensive comparisons at varying levels of detail.

The Genome Explorer genome browser.

Are two adjacent genes in the same operon? What are the order and spacing between several transcription factor binding sites? Genome browsers are software data visualization and exploration tools that enable biologists to answer questions such as these. In this paper, we report on a major update to our browser, Genome Explorer, that provides nearly instantaneous scaling and traversing of a genome, enabling users to quickly and easily zoom into an area of interest. The user can rapidly move between scales that depict the entire genome, individual genes, and the sequence; Genome Explorer presents the most relevant detail and context for each scale. By downloading the data for the entire genome to the user's web browser and dynamically generating visualizations locally, we enable fine control of zoom and pan functions and real-time redrawing of the visualization, resulting in smoother and more intuitive exploration of a genome than is possible with other browsers. Further, genome features are presented together, in-line, using familiar graphical depictions. In contrast, many other browsers depict genome features using data tracks, which have low information density and can visually obscure the relative positions of features. Genome Explorer diagrams have a high information density that provides larger amounts of genome context and sequence information to be presented in a given-sized monitor than for tracks-based browsers. Genome Explorer provides optional data tracks for the analysis of large-scale data sets and a unique comparative mode that aligns genomes at orthologous genes with synchronized zooming. IMPORTANCE: Genome browsers provide graphical depictions of genome information to speed the uptake of complex genome data by scientists. They provide search operations to help scientists find information and zoom operations to enable scientists to view genome features at different resolutions. We introduce the Genome Explorer browser, which provides extremely fast zooming and panning of genome visualizations and displays with high information density.

The Comparative Genome Dashboard.

The Comparative Genome Dashboard is a web-based software tool for interactive exploration of the similarities and differences in gene functions between organisms. It provides a high-level graphical survey of cellular functions, and enables the user to drill down to examine subsystems of interest in greater detail. At its highest level the Comparative Dashboard contains panels for cellular systems such as biosynthesis, energy metabolism, transport, and response to stimulus. Each panel contains a set of bar graphs that plot the numbers of compounds or gene products for each organism across a set of subsystems of that panel. Users can interactively drill down to focus on subsystems of interest and see grids of compounds produced or consumed by each organism, specific GO term assignments, pathway diagrams, and links to more detailed comparison pages. For example, the dashboard enables users to compare the cofactors that a set of organisms can synthesize, the metal ions that they are able to transport, their DNA damage repair capabilities, their biofilm-formation genes, and their viral response proteins. The dashboard enables users to quickly perform comprehensive comparisons at varying levels of detail.

Emergence of terpene chemical communication in insects: Evolutionary recruitment of isoprenoid metabolism.

Insects have evolved a chemical communication system using terpenoids, a structurally diverse class of specialized metabolites, previously thought to be exclusively produced by plants and microbes. Gene discovery, bioinformatics, and biochemical characterization of multiple insect terpene synthases (TPSs) revealed that isopentenyl diphosphate synthases (IDS), enzymes from primary isoprenoid metabolism, are their likely evolutionary progenitors. However, the mutations underlying the emergence of the TPS function remain a mystery. To address this gap, we present the first structural and mechanistic model for the evolutionary emergence of TPS function in insects. Through identifying key mechanistic differences between IDS and TPS enzymes, we hypothesize that the loss of isopentenyl diphosphate (IPP) binding motifs strongly correlates with the gain of the TPS function. Based on this premise, we have elaborated the first explicit structural definition of isopentenyl diphosphate-binding motifs (IBMs) and used the IBM definitions to examine previously characterized insect IDSs and TPSs and to predict the functions of as yet uncharacterized insect IDSs. Consistent with our hypothesis, we observed a clear pattern of disruptive substitutions to IBMs in characterized insect TPSs. In contrast, insect IDSs maintain essential consensus residues for binding IPP. Extending our analysis, we constructed the most comprehensive phylogeny of insect IDS sequences (430 full length sequences from eight insect orders) and used IBMs to predict the function of TPSs. Based on our analysis, we infer multiple, independent TPS emergence events across the class of insects, paving the way for future gene discovery efforts.

Ancient origin and conserved gene function in terpene pheromone and defense evolution of stink bugs and hemipteran insects.

Insects use diverse arrays of small molecules such as metabolites of the large class of terpenes for intra- and inter-specific communication and defense. These molecules are synthesized by specialized metabolic pathways; however, the origin of enzymes involved in terpene biosynthesis and their evolution in insect genomes is still poorly understood. We addressed this question by investigating the evolution of isoprenyl diphosphate synthase (IDS)-like genes with terpene synthase (TPS) function in the family of stink bugs (Pentatomidae) within the large order of piercing-sucking Hemipteran insects. Stink bugs include species of global pest status, many of which emit structurally related 15-carbon sesquiterpenes as sex or aggregation pheromones. We provide evidence for the emergence of IDS-type TPS enzymes at the onset of pentatomid evolution over 100 million years ago, coinciding with the evolution of flowering plants. Stink bugs of different geographical origin maintain small IDS-type families with genes of conserved TPS function, which stands in contrast to the diversification of TPS genes in plants. Expanded gene mining and phylogenetic analysis in other hemipteran insects further provides evidence for an ancient emergence of IDS-like genes under presumed selection for terpene-mediated chemical interactions, and this process occurred independently from a similar evolution of IDS-type TPS genes in beetles. Our findings further suggest differences in TPS diversification in insects and plants in conjunction with different modes of gene functionalization in chemical interactions.

Lettuce costunolide synthase (CYP71BL2) and its homolog (CYP71BL1) from sunflower catalyze distinct regio- and stereoselective hydroxylations in sesquiterpene lactone metabolism.

Sesquiterpene lactones (STLs) are terpenoid natural products possessing the γ-lactone, well known for their diverse biological and medicinal activities. The occurrence of STLs is sporadic in nature, but most STLs have been isolated from plants in the Asteraceae family. Despite the implication of the γ-lactone group in many reported bioactivities of STLs, the biosynthetic origins of the γ-lactone ring remains elusive. Germacrene A acid (GAA) has been suggested as a central precursor of diverse STLs. The regioselective (C6 or C8) and stereoselective (α or ß) hydroxylation on a carbon of GAA adjacent to its carboxylic acid at C12 is responsible for the γ-lactone formation. Here, we report two cytochrome P450 monooxygenases (P450s) capable of catalyzing 6α- and 8ß-hydroxylation of GAA from lettuce and sunflower, respectively. To identify these P450s, sunflower trichomes were isolated to generate a trichome-specific transcript library, from which 10 P450 clones were retrieved. Expression of these clones in a yeast strain metabolically engineered to synthesize substrate GAA identified a P450 catalyzing 8ß-hydroxylation of GAA, but the STL was not formed by spontaneous lactonization. Subsequently, we identified the closest homolog of the GAA 8ß-hydroxylase from lettuce and discovered 6α-hydroxylation of GAA by the recombinant enzyme. The resulting 6α-hydroxy-GAA spontaneously undergoes a lactonization to yield the simplest form of STL, costunolide. Furthermore, we demonstrate the milligram per liter scale de novo synthesis of costunolide using the lettuce P450 in an engineered yeast strain, an important advance that will enable exploitation of STLs. Evolution and homology models of these two P450s are discussed.

Artemisia annua L. plants lacking Bornyl diPhosphate Synthase reallocate carbon from monoterpenes to sesquiterpenes except artemisinin.

The monoterpene camphor is produced in glandular secretory trichomes of the medicinal plant Artemisia annua, which also produces the antimalarial drug artemisinin. We have found that, depending on growth conditions, camphor can accumulate at levels ranging from 1- 10% leaf dry weight (LDW) in the Artemis F1 hybrid, which has been developed for commercial production of artemisinin at up to 1% LDW. We discovered that a camphor null (camphor-0) phenotype segregates in the progeny of self-pollinated Artemis material. Camphor-0 plants also show reduced levels of other less abundant monoterpenes and increased levels of the sesquiterpene precursor farnesyl pyrophosphate plus sesquiterpenes, including enzymatically derived artemisinin pathway intermediates but not artemisinin. One possible explanation for this is that high camphor concentrations in the glandular secretory trichomes play an important role in generating the hydrophobic conditions required for the non-enzymatic conversion of dihydroartemisinic acid tertiary hydroperoxide to artemisinin. We established that the camphor-0 phenotype associates with a genomic deletion that results in loss of a Bornyl diPhosphate Synthase (AaBPS) gene candidate. Functional characterization of the corresponding enzyme in vitro confirmed it can catalyze the first committed step in not only camphor biosynthesis but also in a number of other monoterpenes, accounting for over 60% of total volatiles in A. annua leaves. This in vitro analysis is consistent with loss of monoterpenes in camphor-0 plants. The AaBPS promoter drives high reporter gene expression in A. annua glandular secretory trichomes of juvenile leaves with expression shifting to non-glandular trichomes in mature leaves, which is consistent with AaBPS transcript abundance.

Physical constraints on sesquiterpene diversity arising from cyclization of the eudesm-5-yl carbocation.

The biogenic origins of complex cyclic terpenes derive from the interplay of enzymes and the intrinsic reactivity of carbocation species at major branch-points along intramolecular cyclization pathways to ultimately determine the distribution of terpene skeletal types in nature. Solanaceous plants biosynthesize chemical defense compounds, largely derived from the eremophilane and spirovetivane-type sesquiterpenes. These hydrocarbon skeletons share a common biogenic origin, stemming from alternative Wagner-Meerwein rearrangements of the eudesm-5-yl carbocation during the cyclization of farnesyl pyrophosphate (FPP) catalyzed by sesquiterpene synthases. While the spirojatamane skeleton shares the same carbocation intermediate, this class of sesquiterpenes has not been reported in the Solanaceae and is exceedingly rare in nature. To investigate the physical basis for alternative rearrangements of the eudesm-5-yl carbocation, we carried out quantum mechanics (QM) analyses to calculate the allowable conformations, energies, and transition states linking conformers of the eudesm-5-yl carbocation to the eremophilene, spirovetivane, and spirojatamane skeletons. Additionally, we conducted parallel investigations on simplified decalin carbocation systems to examine the contribution of ring substituents to allowable conformations and rearrangement pathways. Our study reveals that ring substituents expand the conformational space accessible to the eudesm-5-yl carbocation while sterically blocking rearrangements in certain contexts. From our analysis, we define a conformational threshold for each possible rearrangement based on dihedral angles describing transition state geometry. Further, our calculations indicate that methylene migration rearrangements leading to spiro compounds are thermodynamically dominant in the eudesm-5-yl and simpler decalin carabocation systems. Interestingly, the theoretical abundance of sesquiterpene skeletal types arising from the intrinsic reactivity of the eudesm-5-yl carbocation stands in sharp contrast to their currently known natural abundance. The implications of these results for the catalytic tragectories catalyzed by sesquiterpene synthases are discussed.

Bisabolyl-derived sesquiterpenes from tobacco 5-epi-aristolochene synthase-catalyzed cyclization of (2Z,6E)-farnesyl diphosphate.

We report the structures and stereochemistry of seven bisabolyl-derived sesquiterpenes arising from an unprecedented 1,6-cyclization (cisoid pathway) efficiently catalyzed by tobacco 5-epi-aristolochene synthase (TEAS). The use of (2Z,6E)-farnesyl diphosphate as an alternate substrate for recombinant TEAS resulted in a robust enzymatic cyclization to an array of products derived exclusively (>/=99.5%) from the cisoid pathway, whereas these same products account for ca. 2.5% of the total hydrocarbons obtained using (2E,6E)-farnesyl diphosphate. Chromatographic fractionations of extracts from preparative incubations with the 2Z,6E substrate afforded, in addition to the acyclic allylic alcohols (2Z,6E)-farnesol (6.7%) and nerolidol (3.6%), five cyclic sesquiterpene hydrocarbons and two cyclic sesquiterpene alcohols: (+)-2-epi-prezizaene (44%), (-)-alpha-cedrene (21.5%), (R)-(-)-beta-curcumene (15.5%), alpha-acoradiene (3.9%), 4-epi-alpha-acoradiene (1.3%), and equal amounts of alpha-bisabolol (1.8%) and epi-alpha-bisalolol (1.8%). The structures, stereochemistry, and enantiopurities were established by comprehensive spectroscopic analyses, optical rotations, chemical correlations with known sesquiterpenes, comparisons with literature data, and GC analyses. The major product, (+)-2-epi-prezizaene, is structurally related to the naturally occurring tricyclic alcohol, jinkohol (2-epi-prezizaan-7beta-ol). Cisoid cyclization pathways are proposed by which all five sesquiterpene hydrocarbons are derived from a common (7R)-beta-bisabolyl(+)/pyrophosphate(-) ion pair intermediate. The implications of the "cisoid" catalytic activity of TEAS are discussed.

Evolving biosynthetic tangos negotiate mechanistic landscapes.

The dependence of polyketide synthase and terpene cyclase mechanistic adaptation on the chemistry of their oligomeric substrates illuminates a convergent evolutionary strategy for shaping cyclization in these otherwise disparate reactions. Evolution of these enzyme families relies on rhythmic tangos, in which the enzymes and substrates together determine product outcome by negotiating decision networks governing intrinsic and induced chemical reactivities.

Quantitative exploration of the catalytic landscape separating divergent plant sesquiterpene synthases.

Throughout molecular evolution, organisms create assorted chemicals in response to varying ecological niches. Catalytic landscapes underlie metabolic evolution, wherein mutational steps alter the biosynthetic properties of enzymes. Here we report the first systematic quantitative characterization of the catalytic landscape underlying the evolution of sesquiterpene chemical diversity. On the basis of our previous discovery of a set of nine naturally occurring amino acid substitutions that functionally interconverted orthologous sesquiterpene synthases from Nicotiana tabacum and Hyoscyamus muticus, we created a library of all possible residue combinations (2(9) = 512) in the N. tabacum enzyme. The product spectra of 418 active enzymes revealed a rugged landscape where several minimal combinations of the nine mutations encode convergent solutions to the interconversions of parental activities. Quantitative comparisons indicated context dependence for mutational effects--epistasis--in product specificity and promiscuity. These results provide a measure of the mutational accessibility of phenotypic variability in a diverging lineage of terpene synthases.