Structural Basis for Cyclic Terpene Biosynthesis by Tobacco 5-Epi-Aristolochene Synthase

See allHide authors and affiliations

Science  19 Sep 1997:
Vol. 277, Issue 5333, pp. 1815-1820
DOI: 10.1126/science.277.5333.1815


Terpene cyclases catalyze the synthesis of cyclic terpenes with 10-, 15-, and 20-carbon acyclic isoprenoid diphosphates as substrates. Plants have been a source of these natural products by providing a homologous set of terpene synthases. The crystal structures of 5-epi-aristolochene synthase, a sesquiterpene cyclase from tobacco, alone and complexed separately with two farnesyl diphosphate analogs were analyzed. These structures reveal an unexpected enzymatic mechanism for the synthesis of the bicyclic product, 5-epi-aristolochene, and provide a basis for understanding the stereochemical selectivity displayed by other cyclases in the biosynthesis of pharmacologically important cyclic terpenes. As such, these structures provide templates for the engineering of novel terpene cyclases.

Terpene cyclases control the synthesis of cyclic terpenoids including flavors and fragrances such as menthol and camphor, plant defense chemicals like capsidiol and lubimin (1), and more common compounds like steroids and lipid-soluble vitamins. Several cyclic terpenoids have pharmacological activity; for example, limonene can inhibit tumorigenesis induced in mice by particular carcinogens (2), and the diterpenoid taxol has antitumor activity (3). Numerous terpene cyclases from plant and microbial sources have been characterized (4, 5). Although the plant cyclases exhibit a significant degree of similarity in amino acid sequence, very little similarity is observed between the bacterial, fungal, and plant terpene cyclases (6). These soluble enzymes convert the acyclic isoprenoid diphosphates geranyl diphosphate (GPP, 10 carbon), farnesyl diphosphate (FPP, 15 carbon), and geranylgeranyl diphosphate (GGPP, 20 carbon) into cyclic monoterpenes, sesquiterpenes, and diterpenes, respectively. In most cases, loss of diphosphate (pyrophosphate) from the enzyme-bound acyclic substrate results in an allylic carbocation that electrophilically attacks a double bond further down the terpene chain to effect the first ring closure. Additional rearrangements involving transient carbocations can include proton abstractions, hydride and alkyl migrations, and additional electrophilic attacks.

TEAS (tobacco 5-epi-aristolochene synthase) (7) fromNicotiana tabacum converts farnesyl diphosphate (FPP) to 5-epi-aristolochene (Fig. 1) (8), a precursor of the antifungal phytoalexin capsidiol. TEAS shares 77% amino acid identity with Hyoscyamus muticus vetispiradiene synthase (HVS) (9). Vetispiradiene (Fig. 1) is a precursor to the phytoalexins solavetivone and lubimin. Both enzymes have similar reaction mechanisms (9); in fact, several TEAS-HVS chimeras produce mixtures of the natural reaction products 5-epi-aristolochene and vetispiradiene (10).

Figure 1

(A) Reactions catalyzed by TEAS and HVS. (B) Schematic view of the TEAS·farnesyl hydroxyphosphonate (FHP) complex. Blue rods represent α-helices in the NH2-terminal domain; orange rods represent α-helices in the COOH-terminal domain. Loop regions shown in green are disordered in the native TEAS structure. The NH2-terminal domain forms a twisted α-barrel made up of eight short (10 to 15 residues) helices surrounding a surface cavity filled by ordered water molecules. The COOH-terminal domain consists of a two layered barrel of α-helices surrounding a hydrophobic and aromatic-rich active site pocket. Helix H is disrupted between segments H1 and H2 by Pro430, but its interhelical packing with helix G is accommodated by a corresponding kink in helix G between G1 and G2 (residues 400 to 404). Within this kink, hydrogen bonds between the hydroxyl group of Thr401 and the carbonyl oxygen of residue 398, and between the hydroxyl group of Thr402 and the carbonyl oxygen of residue 399 disrupt the main chain intrahelical hydrogen bonding of helix G. Naming of helices in the COOH-terminal domain corresponds to the convention used for FPP synthase (FPS) (15). (C) Surface electrostatic representation of the hydrophobic active site pocket of TEAS in the closed, FHP-bound structure. The surface is colored by electrostatic potential, with blue corresponding to positive, red to negative, and white to neutral regions. FHP is shown in stick representation with carbons in green. Panel (B) was generated with RIBBONS (24) and panel (C) was prepared with GRASP (24).

In order to better understand terpene cyclase mechanisms, we have determined x-ray crystal structures at 2.2 to 2.8 Å resolution of TEAS alone and in complexes with two different substrate analogs. Recombinant TEAS was expressed in Eschericia coli, purified and crystallized (11). We determined the initial structure for a TEAS-HVS chimeric cyclase (CH3) (10) by multiple isomorphous replacement with anomalous scattering (MIRAS) with data to 2.8 Å resolution (Table 1). After density modification, a nearly complete CH3 model was built and refined against the CH3 diffraction data to 2.8 Å. This model was then used for construction of the initial TEAS model with diffraction data for native TEAS crystals extending to 2.25 Å (Table 1). The TEAS model was refined to an R factor of 19.9% with good stereochemistry (12).

Table 1

Data collection, structure determination, and refinement. Diffraction data for chimeric and derivative crystals were collected on a MacScience imaging plate detector DIP2020K (MacScience Corp.) with double focusing Pt/Ni coated mirrors and CuKα x-rays. All other data sets were collected at the Stanford Synchrotron Radiation Laboratory, beamline 7-1 (λ = 1.08 Å) on a MAR imaging plate system. Crystals were stabilized by soaking in synthetic mothor liquor containing 20% ethylene glycol, and subsequently frozen in a 100 K nitrogen gas stream. Data sets were processed with DENZO (25) and scaled with SCALEPACK (25). The structure was solved by MIRAS; heavy atom sites were located with difference Patterson and difference Fourier maps. Sites were refined and initial phases were calculated with ML-PHARE (26), giving a figure of merit of 0.39. This initial phase set was improved and the phases extended to 2.8 Å by solvent flattening and histogram matching with the program DM (27), assuming a solvent content of 65%. All models were constructed with program O (28) and refined with XPLOR (29).

View this table:

TEAS consists entirely of α-helices and short connecting loops and turns, and is organized into two structural domains (Fig. 1B). The backbone of the NH2-terminal domain (residues 36 to 230) aligns structurally with the catalytic cores of two glycosyl hydrolases: glucoamylase (PDB code 3GLY) from Aspergillus awamori (13), and endoglucanase CelD (PDB code 1CLC) from Clostridium thermocellum. (14). The function of this domain in TEAS is unknown. The COOH-terminal domain aligns structurally with avian FPP synthase (FPS), despite their lack of sequence similarity (15).

The enzyme active site was identified in the hydrophobic pocket of the COOH-terminal domain by the location of two Mg2+coordination sites at its opening, by the location of bound substrate analogs, and by comparison of the structure with the published FPS structures (15). Two Mg2+ ions are coordinated on opposite sides of the entrance to the active site pocket in positions analogous to those in FPS, and constitute a diphosphate binding site. Asp301 coordinates Mga 2+ in the native TEAS structure, and the side chain carboxyl of Glu379 provides a longer range interaction. Asp305 provides an additional coordination bond in the enzyme with substrate analogs bound. Asp301 and Asp305 are part of a -DDXXD- sequence found in terpene cyclases and were previously assumed to constitute a binding site for the required divalent metal (16). Asp301directly contacts Mga 2+, whereas Asp302 demonstrates no direct metal coordination. The side chains of Asp444, Thr448, Glu452, and one water molecule coordinate Mgb 2+. In the native TEAS structure the A-C and J-K loops and the residues NH2-terminal of residue 36 are disordered (Fig.1B). Given their location near the entrance to the active site pocket, the flexibility of these regions results in an open active site that displays a number of solvent accessible hydrophobic residues poised to bind the hydrophobic farnesyl chain of FPP.

In order to pinpoint the contributions of particular side chains to substrate and product specificity and to elucidate the stereochemical control exerted by TEAS during farnesyl cyclization, we examined the complexes between TEAS and two FPP analogs: a farnesyl hydroxyphosphonate (FHP), and a trifluoro-farnesyl diphosphate (F3-FPP). Complexes with each were obtained both by incubation of TEAS with analog prior to crystallization and by soaking native crystals in solutions containing the analogs.

The native TEAS model, lacking the J-K loop and the residues NH2-terminal of residue 36, was refined against data from the TEAS·FHP complex (17). Difference electron density maps calculated with phases derived from this partial model revealed new electron density for FHP, and the previously disordered NH2-terminal segment, A-C loop, and J-K loop (Fig. 2, top panel). The ordering of these loops upon FHP binding results in a closed, solvent-inaccessible active site pocket, as expected for binding the hydrophobic substrate and for protecting reactive carbocation intermediates from attack by water (Fig. 1C). As the J-K loop (residues 521 to 534) becomes ordered, it forms a lid that clamps down over the active site entrance in the presence of FHP. This clamp places Tyr527 next to Trp273, forming an extended aromatic box deep within the active site pocket. The A-C loop, which contains Arg264 and Arg266, translates inward toward the active site on FHP binding, positioning the side chain of Arg264 in close proximity to the C1 hydroxyl group of FHP. Arg266 hydrogen bonds with residues on both the J-K loop and the NH2-terminal segment in the FHP complex. Both Arg264 and Arg266 are conserved among many terpene cyclases. On FHP binding, a third Mg2+ site, designated Mgc 2+, forms. The single phosphonate group of FHP binds between Mgb 2+ and Mgc 2+. Much of the FHP farnesyl tail exhibits strong electron density as it extends back into the hydrophobic active site.

Figure 2

Stereo views of the TEAS active site. Blue dashed lines are hydrogen bonds; green dashed lines are coordination bonds. The top panel illustrates the TEAS·FHP complex. FHP is emphasized with dark gray bonds. The lower panel illustrates the TEAS·F3FPP complex. The diphosphate moiety is emphasized with dark bonds. The rest of the F3-FPP molecule, which exhibits weak electron density, is indicated with dashed bonds. Figure was prepared with MOLSCRIPT (24) and RASTER3D (24).

In the TEAS·F3-FPP complex (18), we observe a well-ordered diphosphate binding pocket; the A-C loop and the NH2-terminal segment exhibit well-defined electron density, the A-C loop translates toward the active site, and there is strong electron density for the F3-FPP diphosphate moiety. The hydrophobic pocket, however, remains flexible; the J-K loop and the F3-FPP farnesyl chain are disordered (Fig.2, bottom panel).

The positions of the two analogs in the TEAS active site suggest a binding mode for the natural substrate FPP that is consistent with the known product specificity of TEAS. Both the phosphonate of FHP and the β-phosphate of F3-FPP reside in nearly identical positions between Mgb 2+ and Mgc 2+; the C1 hydroxyl of FHP also lies near the position of the α-phosphate of F3-FPP, near Mgc 2+. Utilizing these structural constraints, we modeled the binding mode for FPP (Fig.3A). This model of TEAS complexed with its natural substrate defines our proposed catalytic mechanisms for 5-epi-aristolochene biosynthesis in TEAS, and vetispiradiene biosynthesis in the related enzyme HVS, that are consistent with the cyclization pathways proposed on the basis of biochemical data (9).

Figure 3

Proposed catalytic mechanisms of TEAS and HVS based on the TEAS crystal structures. (A) through (F) are reaction steps common to both TEAS and HVS. (G1 ) through (I1 ) are specific to TEAS, and (G2 ) through (I2 ) are specific to HVS. Carbocations are yellow, protons are light blue, and aromatic quadrupoles are indicated with purple lines. Intermediate structures are numbered according to the acyclic precursor, FPP (A). Figure was generated with Chem3D (CambridgeSoft).

In our proposed mechanism, as diphosphate is generated (Fig. 3A), its additional negative charge would be offset by interactions with the three Mg2+ ions, Arg264, and Arg441. This concentrated region of positive charge may serve to direct the diphosphate away from the hydrophobic pocket where a number of highly reactive carbocations form throughout the reaction. The released diphosphate, if not constrained, would likely capture the newly formed allylic carbocation at C1 or C3, thereby regenerating FPP (C1 attack) or its tertiary allylic isomer nerolidyl diphosphate (C3 attack). The allylic carbocation of the intermediate is positioned near the main chain carbonyl oxygens of residues 401 and 402 at the kink in helix G. These carbonyls are not involved in hydrogen bonding within the helix; instead, their dipole moments are directed toward the positive charge delocalized over C1, C2, and C3. Movement of the allylic carbocation deeper into the active site would improve this dipolar stabilization by positioning the carbocation closer to the peptide carbonyls of residues 401 and 402 as well as closer to the dipole of Thr403 hydroxyl group (Fig. 3B). Just as the Mg2+ and arginines directed the diphosphate away from the active site, these dipoles would direct the cationic end of the farnesyl chain into the hydrophobic active site, poised for attack on C10.

Once C1 has been positioned near the π orbitals of the C10–C11 bond, electrophilic attack at C10 would create a C1–C10 bond and a tertiary carbocation on C11 (Fig. 3C). The quadrupole of Tyr527 is nicely positioned to stabilize the positive charge on C11. In turn, the newly formed carbocation at C11 will substantially increase the acidity of the C12 and C13 methyl protons. As modeled, the carboxyl group of Asp525 would abstract a proton from C13 (cis methyl group) leading to the formation of a neutral germacrene intermediate (Fig. 3D).

Subsequently, the reaction requires proton addition at C6. In the closed conformation observed in the TEAS·FHP complex, an Asp444-Tyr520-Asp525 triad organizes on the interior surface of the active site as the J/K loop, containing Tyr520 and Asp525, clamps over the active site. A carboxyl oxygen from each Asp side chain is within hydrogen bonding distance of Tyr520's hydroxyl group (Fig. 2, top panel). In the proposed TEAS·germacrene complex, the Tyr520 hydroxyl group also resides near C6 of germacrene (Fig. 3D). In our proposed reaction scheme, the carboxyl group of Asp444 would remove the hydroxyl proton from Tyr520. At the same time, the phenolic oxygen at Tyr520 would accept a proton from Asp525. Tyr520 then could donate this newly positioned proton to the double bond of germacrene at C6 (Fig.3E). In a concerted fashion, Tyr520 would re-accept the proton currently on the Asp444 carboxyl group. Concomitant with the C6 protonation by Tyr520, germacrene would undergo a second ring closure as the π orbitals on C2 and C7 line up for the formation of the C2–C7 sigma bond. The resulting bicyclic eudesmane carbocation intermediate (Fig. 3F) with the positive charge on C3 would be stabilized by the summed dipoles of the main chain carbonyls of residues 401 and 402 and the Thr403 hydroxyl group.

In order to ensure either 5-epi-aristolochene (TEAS) or vetispiradiene (HVS) synthesis, both enzymes subsequently stabilize a carbocation centered on C7. In both cases, the C2 hydride first migrates to the planar C3 carbocation along the top of the bicyclic ring system leaving a tertiary carbocation with planar geometry at C2 (Fig. 3F). C7 abuts the aromatic face of Trp273, and stabilization afforded by this aromatic quadrupole would selectively position the electrophilic center at C7. Shifting the carbocation to C7 from C2 occurs via two different routes that distinguish the TEAS catalyzed synthesis of 5-epi-aristolochene (Fig. 3, G1 to I1) and the HVS catalyzed synthesis of vetispiradiene (Fig. 3, G2 to I2). In TEAS, the eudesmane carbocation intermediate with the positive charge centered at C2 would exist in a conformation that orients the C7–C14 sigma bond parallel to the empty p orbital on C2. This conformer would facilitate the migration of the C14 methyl group from C7 to C2 (Fig. 3G1). The conformation of the resulting carbocation with a positive charge positioned on C7 directs a proton on C8 toward the indole ring of Trp273. The presence of the carbocation center on C7 greatly increases the acidity of the proton at C8 which may now be removed by Trp273 (Fig.3H1), giving rise to a positive arenium ion (TrpH+ at residue 273) and the final reaction product, 5-epi-aristolochene (Fig. 3I1).

In vetispiradiene synthesis by HVS, an energetically allowed 60° rotation around the C2–C7 bond of the eudesmane carbocation intermediate with the positive charge centered on C2 would result in a conformation in which the C7–C8 bond is now parallel to the empty p orbital of C2 (Fig. 3G2). In this conformation, the methylene group C8 migrates from C7 to C2, resulting in the collapse of the original six-membered ring to a five-membered ring. The resulting carbocation intermediate directs a proton of C6 toward the indole ring of Trp273. As described for TEAS, the greatly enhanced acidity of this proton could allow facile deprotonation by Trp273 (Fig. 3H2), giving rise to a positive arenium ion (TrpH+residue 273), and the reaction product vetispiradiene (Fig. 3I2). In both enzymes, the active site would again become solvent-accessible following product dissociation as the A-C and J-K loops and residues at the NH2-terminus become mobile. The TrpH+cation would then lose a proton to the surrounding aqueous environment.

Trp273 is postulated to act as a proton acceptor during the final deprotonation step for two reasons. First, the size and positioning of the indole ring allows it to accept an acidic proton from either C8 (Fig. 3H1) or C6 (Fig. 3H2) as would be required in the chimeric TEAS or HVS enzymes that catalyze the synthesis of both 5-epi-aristolochene and vetispiradiene (10). No other proton acceptors are positioned to abstract a proton from C6 or C8 without a substantial reorientation of one or more of the cyclic intermediates. Mutation of Trp273 to either nonaromatic or aromatic residues results in cessation of epi-aristolochene or vetispiradiene synthesis (19). Second, the greatly reduced nucleophilicity of the indole ring, compared to general bases such as histidine, aspartate, or glutamate, greatly diminishes the likelihood of inadvertent enzyme alkylation by the highly reactive carbocation intermediate. While the enzymatic roles proposed here for Trp273 are novel, there is chemical precedence for both quadrupole mediated cation interactions (20) and for electrophilic aromatic substitutions in indole ring systems (21). In the solvent shielded active site of TEAS or HVS, the arenium cation formed upon protonation of the indole ring of Trp273 would be more energetically stable than the eudesmane carbocations, due to the presence of the nitrogen atom and the extended conjugation of the indole ring.

Modeling the structure of HVS based on our crystallographic model of TEAS suggests that the residues which directly line the active site are conserved between the two enzymes; product selectivity, therefore, likely derives from the amino acid differences in layers surrounding the active site. This pattern of conservation suggests that the catalytic residues required to form the two products are the same, and that specificity depends on having a particular active site conformation determined by the surrounding layers. This model is consistent with the observation that TEAS and HVS chimeric enzymes can produce mixtures of 5-epi-aristolochene and vetispiradiene. These multiproduct chimeric enzymes, in which the active site topology is influenced by a combination of interactions from the two wild-type enzymes, still catalyze terpene cyclization, but have lost the ability to select a single conformation of the eudesmane intermediate that would give rise to a single product.

Sequence alignments between TEAS and other plant cyclases viewed in the context of the TEAS structure suggest that all plant cyclases share similar three-dimensional structures. For example, the large sequence insertions in taxadiene synthase (TS) (22) with respect to TEAS occur at the NH2-terminus, or within solvent exposed loops in the NH2-terminal domain. TS converts geranylgeranyl diphosphate to the tricyclic taxol precursor, taxadiene. Modeling the TS structure on the basis of the three-dimensional fold of TEAS suggests that the TS active site is more spacious than that of TEAS, consistent with the larger substrate of TS. The TS active site is also rich in aromatic residues, including Phe602 (Leu290 in TEAS), Trp753(Cys440 in TEAS), Phe835 (Thr519 in TEAS), and Tyr841 (Tyr527 in TEAS). Given the likely involvement of aromatic quadrupoles in carbocation stabilization, this extensive aromatic surface may reflect the greater spatial distribution of carbocations on various intermediates in the proposed TS reaction mechanism (23). In closing, the high resolution structures and comparative models that we describe provide a structural foundation for efforts to create novel terpene cyclases.

  • * To whom correspondence should be addressed. E-mail: noel{at}


View Abstract

Stay Connected to Science

Navigate This Article