Module 6: Secondary Metabolites & Chemical Defense

Trees invest 2–10% of their dry mass in secondary metabolites — compounds with no primary metabolic role but crucial ecological functions: defense against herbivores and pathogens, attraction of pollinators and seed dispersers, allelopathy, and UV screening. The major biosynthetic pathways (phenylpropanoid, terpenoid, alkaloid) are tightly integrated with primary carbon metabolism and regulated by jasmonate, ethylene, and salicylate signals.

6.1 Phenylpropanoid Pathway

The phenylpropanoid pathway converts phenylalanine into an enormous diversity of phenolic compounds: lignin monomers, flavonoids, stilbenes, coumarins, tannins, and hydroxycinnamic acid esters. It is the single largest metabolic sink for fixed carbon in woody tissues, with lignin alone comprising 15–35% of wood dry mass.

PAL: The Committed Step & MIO Mechanism

Phenylalanine ammonia-lyase (PAL) catalyzes the non-oxidative deamination of L-phenylalanine to trans-cinnamic acid:

\[ \text{L-Phe} \xrightarrow{\text{PAL}} \text{trans-cinnamate} + \text{NH}_3 \]

PAL is unique among ammonia-lyases in using a 4-methylideneimidazol-5-one (MIO) prosthetic group formed autocatalytically from the tripeptide Ala-Ser-Gly (residues 142–144 in Arabidopsis PAL). The MIO electrophile abstracts the pro-S proton from the β-carbon of L-Phe via a Friedel-Crafts-type electrophilic substitution, then eliminates ammonia to form the double bond. PAL is a homotetramer (~330 kDa) in most plants; trees typically have 4–8 PAL genes with different induction kinetics.

Core Pathway: C4H → 4CL → CHS/CHI

trans-Cinnamate → p-Coumarate (C4H)

Cinnamate 4-hydroxylase (C4H, CYP73A1) is a cytochrome P450 monooxygenase that inserts an oxygen atom at the para-position, using NADPH and O₂ (reductive O₂ activation via the Fe²⁺-heme center). Product: p-coumaric acid. C4H is ER-membrane anchored and forms a metabolon with PAL and 4CL.

p-Coumarate → CoA Ester (4CL)

4-Coumarate:CoA ligase (4CL) activates p-coumarate (and caffeate, ferulate, sinapate) to CoA thioesters in a 2-step ATP-dependent reaction: first forming an adenylate intermediate, then the thioester. Products are the activated substrates for CHS, HCT, and monolignol biosynthesis.

p-Coumaroyl-CoA → Naringenin Chalcone (CHS)

Chalcone synthase (CHS) catalyzes iterative Claisen condensations: one molecule of p-coumaroyl-CoA + three malonyl-CoA → naringenin chalcone (C15) + 4 CoA + 3 CO₂. A type III PKS (polyketide synthase), CHS uses a Cys-His-Asn catalytic triad and a hydrophobic cavity that determines starter and extension unit selectivity.

Chalcone → Naringenin (CHI)

Chalcone isomerase (CHI) catalyzes the stereospecific intramolecular Michael addition of the 4'-OH to the α,β-unsaturated carbonyl, forming the flavanone ring (ring closure). Product: (2S)-naringenin. Branching from naringenin leads to flavonols, anthocyanins, isoflavonoids, and condensed tannins.

6.2 Monolignol Biosynthesis & Lignin Polymerization

Lignin is a racemic heteropolymer formed by oxidative coupling of three monolignol monomers synthesized via the phenylpropanoid/monolignol pathway. The three monomers differ in the number of methoxyl groups on the aromatic ring:

p-Coumaryl alcohol (H-unit)

Methoxyl groups: 0

Compression wood, grasses

Coniferyl alcohol (G-unit)

Methoxyl groups: 1 (at C3)

Gymnosperms (dominant ~90%), hardwood G/S

Sinapyl alcohol (S-unit)

Methoxyl groups: 2 (at C3 & C5)

Angiosperms (hardwoods, 30–50%)

Monolignols are synthesized from p-coumaroyl-CoA via a series of hydroxylation (CYP84A1/CCoAOMT), methylation, and reduction steps. The final two steps — CAD (cinnamyl alcohol dehydrogenase) and CCR (cinnamoyl-CoA reductase) — use NADPH to produce the alcohol form. Monolignols are exported to the cell wall via ABC transporters and vesicle-mediated exocytosis.

Radical Coupling & Lignin Network Formation

Polymerization is driven by class III peroxidases (and laccases) that use H₂O₂ (peroxidases) or O₂ (laccases) to generate monolignol phenoxy radicals. These radicals undergo spontaneous combinatorial coupling (β-O-4, β-β, β-5, 5-5, 4-O-5 linkages). The resulting lignin is non-repetitive, branched, and covalently crosslinks with polysaccharides (cellulose microfibrils via ester linkages to glucuronoxylan).

The most abundant interunit linkage in gymnosperm (G-lignin) is β-aryl ether (β-O-4), comprising ~50% of all bonds. This ether bond is the primary target of lignocellulose biorefinery enzymes (β-etherases from Novosphingobium):

\[ \text{G-unit-}\beta\text{-O-4-G} \xrightarrow{\text{LigE/LigF/LigG}} \text{GGE} \rightarrow \text{monomers (GVL, GPE)} \]

Gymnosperm lignin (G-lignin) is more condensed than angiosperm G/S-lignin because the C5 position is free in G-units, enabling 5-5 and β-5 condensed bonds. S-units (two methoxyl groups) have blocked C3 and C5, forcing primarily linear β-O-4 chains — making hardwood lignin easier to depolymerize.

6.3 Terpenoid Biosynthesis: MVA vs MEP Pathways

All terpenoids are built from two universal 5-carbon precursors: isopentenyl diphosphate (IPP) and dimethylallyl diphosphate (DMAPP). These are synthesized by two independent and compartmentally separated pathways in plants:

MVA Pathway (cytosol) vs MEP Pathway (plastid)MVA Pathway (Cytosol)Acetyl-CoA (x3)HMG-CoAMevalonate (MVA)Mevalonate-5-PMevalonate-5-PPIPP (cytosol)DMAPPGPP (C10), FPP (C15)Sterols, sesquiterpenesAACTHMGSHMGR*MVKPMKMVDIDIFPPSdiverse TPSMEP Pathway (Plastid)Pyruvate + GAP (G3P)DXPMEPCDP-MECDP-MEPMECIPP (plastid)GGPP (C20)Diterpenes, carotenoids, GADXS*DXRCMSCMKMCSHDS,HDRIDI2GGPPSCPS, KS, ...IPP exchange(limited)

The MVA pathway (cytosol/ER) provides IPP/DMAPP for sesquiterpenes (C15), triterpenes (C30), and sterols. The MEP pathway (plastids) supplies IPP for monoterpenes (C10), diterpenes (C20), and tetraterpenes/carotenoids (C40). Limited cross-talk occurs via MEP/IPP exchange across the plastid inner envelope. Asterisks mark rate-limiting enzymes (HMGR and DXS).

MVA Pathway: 6 Steps from Acetyl-CoA

1. AACT (Acetyl-CoA acetyltransferase)2 Acetyl-CoA → Acetoacetyl-CoA + CoA — thiolase condensation
2. HMGS (HMG-CoA synthase)Acetoacetyl-CoA + Acetyl-CoA + H₂O → HMG-CoA + CoA — irreversible Claisen condensation
3. HMGR* (HMG-CoA reductase)HMG-CoA + 2 NADPH → Mevalonate + CoA — RATE LIMITING, target of statins
4. MVK (Mevalonate kinase)Mevalonate + ATP → Mevalonate-5-phosphate
5. PMK (Phosphomevalonate kinase)MVA-5-P + ATP → MVA-5-PP
6. MVD (Mevalonate diphosphate decarboxylase)MVA-5-PP + ATP → IPP + CO₂ + ADP + Pi — decarboxylation

MEP Pathway: 7 Steps from Pyruvate + G3P

1. DXS* (DXP synthase)Pyruvate + GAP → DXP (1-deoxy-D-xylulose-5-phosphate) — TPP-dependent, RATE LIMITING
2. DXR (DXP reductoisomerase)DXP → MEP (2-C-methyl-D-erythritol-4-P) — NADPH-dependent isomerization/reduction, target of fosmidomycin
3. CMS (CDP-ME synthase)MEP + CTP → CDP-ME + PPi — nucleotidylation
4. CMK (CDP-ME kinase)CDP-ME + ATP → CDP-MEP — phosphorylation
5. MCS (MEC synthase)CDP-MEP → MEC (2-C-methyl-D-erythritol-2,4-cyclodiphosphate) + CMP
6. HDS (HMBPP synthase)MEC + NADPH + Fd_red → HMBPP — [4Fe-4S] enzyme
7. HDR (HMBPP reductase)HMBPP + 2 Fd_red → IPP + DMAPP (4:1 ratio) — [4Fe-4S]

Prenyl Transferases: Assembling the Carbon Skeleton

IPP and DMAPP are condensed by prenyl transferases in a head-to-tail manner via ionization of the allylic diphosphate (DMAPP) and electrophilic addition to the double bond of IPP (SN1-like mechanism):

\[ DMAPP_{C5} + IPP_{C5} \xrightarrow{GPPS} GPP_{C10} \xrightarrow{+IPP, FPPS} FPP_{C15} \xrightarrow{+IPP, GGPPS} GGPP_{C20} \]

GPP (C10)

via GPPS

Monoterpenes: limonene, α/β-pinene, myrcene, linalool

FPP (C15)

via FPPS

Sesquiterpenes: farnesol, α-humulene, β-caryophyllene, abscisic acid

GGPP (C20)

via GGPPS

Diterpenes: taxol, gibberellins, abietate, copalate (resin acids)

Squalene (C30)

via SQS (2× FPP)

Triterpenes: sterols, triterpenoid saponins, betulin

6.4 Conifer Oleoresin: Monoterpene Synthases & Traumatic Resin Ducts

Conifers produce oleoresin — a mixture of monoterpenes (volatile "turpentine" fraction) and diterpene resin acids (non-volatile "rosin" fraction) — as a highly effective defense against bark beetles and associated fungal pathogens. The physical properties of oleoresin (viscosity, crystallization kinetics) are as important as its toxicity: resin pressurizes the resin canal system (often 3–20 bar), and when a beetle attacks, the resin physically pitches out the beetle and its fungal inoculant.

Monoterpene Synthases: Multi-product Enzymes

Monoterpene synthases (MTSs) convert GPP to cyclic or acyclic monoterpenes. Many are multi-product enzymes; for example, γ-terpinene synthase from Citrus sinensisgenerates 22 distinct monoterpene products from a single GPP substrate. The reaction begins with ionization of GPP diphosphate by a Mg²⁺ cluster (DDxxD motif), forming the geranyl carbocation. A conformational change induced by closure of the active site cap then governs the specific cyclization cascade.

(−)-α-Pinene synthase

Product: α-Pinene (bicyclic)

Major component of Norway spruce oleoresin; also aggregation kairomone for some bark beetles

(−)-β-Pinene synthase

Product: β-Pinene

Bicyclic, also β-phellandrene; widespread in Pinus spp.; strong antifungal activity

Limonene synthase

Product: (+)/(−)-Limonene

Monocyclic; major Abies and Citrus component; potent insect deterrent

Diterpene Resin Acids: Abietic & Pimaric

Diterpene resin acids (DRAs) are the dominant non-volatile components of conifer resin (40–90% of rosin in Pinus). Their biosynthesis begins with GGPP → copalyl diphosphate (CPP) via copalyl diphosphate synthase (CPS, a class II diterpene synthase), followed by cyclization to different diterpene skeletons via class I diterpene synthases (KS, isopimaradiene synthase, palustradiene synthase), and then three oxidation steps by CYP720B cytochrome P450s to yield the carboxylic acid.

Abietic acid (the most abundant DRA in many species) has antimicrobial activity against bark beetle-associated fungi (Ophiostoma spp.) with MIC values of ~50–200 μg/mL. Its mechanism involves disruption of fungal plasma membrane ergosterol — structurally analogous to plant sterol disruption by saponins.

Traumatic Resin Duct Formation

Conifers maintain constitutive resin ducts in their axial and radial xylem. However, following bark beetle attack, mechanical wounding, or pathogen inoculation, they can produce massive induced traumatic resin ducts (TRDs) — large, resin-filled cavities surrounded by epithelial cells that secrete oleoresin. TRD formation is driven by a jasmonate + ethylene synergistic signal:

Wounding/beetle attack → OPDA/JA-Ile accumulation (within hours) + ET burst (peak ~24 h)
→ COI1-JAZ co-receptor activation → MYC2/MYC3 TF activation
→ Upregulation of TPS (terpene synthases), CYP720B (acid oxidation), GGPPS
→ Resin biosynthesis ×5–20 above constitutive levels
→ Cambial cell proliferation → traumatic resin duct formation (3–7 days post-attack)

6.5 Volatile Organic Compound (VOC) Emissions: The Guenther Algorithm

Forests emit approximately 500 Tg C yr⁻¹ as biogenic VOCs (BVOCs), with isoprene alone accounting for ~70% of total BVOC flux. These emissions influence atmospheric chemistry (ozone and SOA formation, OH radical scavenging) and regulate leaf temperature through evaporative cooling.

Isoprene Synthase & the Thermal Protection Hypothesis

Isoprene is synthesized in chloroplasts from DMAPP by isoprene synthase (IspS), a class I terpene synthase that is among the fastest known (kcat ~0.6 s⁻¹). IspS has an unusually high thermal optimum (~46°C) and its expression is tightly regulated by temperature, light, and CO₂. One proposed function is thermal protection: isoprene partitions into membrane lipid bilayers and quenches thermally generated reactive oxygen species, stabilizing thylakoid membranes during transient heat pulses.

Guenther Emission Algorithm

The Guenther (1993, 1995) algorithm describes isoprene emission rate as:

\[ E = E_s \cdot C_T \cdot C_L \]

Temperature function:

\[ C_T = \frac{\exp\!\left(\frac{C_{T1}(T - T_s)}{RT T_s}\right)}{1 + \exp\!\left(\frac{C_{T2}(T - T_m)}{RT T_s}\right)} \]

where \( C_{T1} = 95000 \) J mol⁻¹, \( C_{T2} = 230000 \) J mol⁻¹, \( T_s = 303 \) K, \( T_m = 314 \) K.

Light function:

\[ C_L = \frac{\alpha C_{L1} \cdot PAR}{\sqrt{1 + \alpha^2 \cdot PAR^2}} \]

where \( \alpha = 0.0027 \), \( C_{L1} = 1.066 \), PAR in μmol m⁻² s⁻¹.

The isoprene emission function shows a characteristic optimum near 40°C with a steep decline above 42°C (denominator becomes large as \( T \rightarrow T_m \)). Monoterpene emissions follow a simpler exponential temperature dependence (beta model):\( E_{mono} = E_s \cdot e^{\beta(T-T_s)} \) with \( \beta \approx 0.09 \) K⁻¹, because stored monoterpene pools drive emissions without requiring active photosynthate.

6.6 Bark Tannins, Stilbenes & Wood Durability

Condensed Tannins (Proanthocyanidins)

Condensed tannins (CTs) are oligomeric/polymeric flavan-3-ol polymers linked by C–C (4→8 or 4→6) or C–O–C (4→8 ether) bonds. Monomers include catechin (procyanidin B precursor) and afzelechin (propelargonidin). CTs bind and precipitate proteins (herbivore digestive enzymes, collagen), chelate metals (Fe, Cu), and are powerful antioxidants. In Scots pine bark, CTs can reach 10–25% dry mass.

Biosynthesis from naringenin: CHI → eriodictyol → leucoanthocyanidins (via F3H, F3'H, DFR) → flavan-3-ols (via ANR/LAR) → condensation by TT12/MATE transporter-directed vacuolar polymerization.

Stilbenes: Resveratrol in Pine Heartwood

Stilbenes are C14 polyphenols synthesized by stilbene synthase (STS), a type III PKS that competes with CHS for p-coumaroyl-CoA + malonyl-CoA substrates. The product pinosylvin (found in pine heartwood) and resveratrol (mainly in grapevine) are antifungal through membrane disruption and inhibition of cytochrome bc₁ complex (Complex III of the respiratory chain) in fungi.

The deposition of pinosylvin, pinosylvin methyl ether, and diterpene resin acids during heartwood formation (the programmed cell death event in aging sapwood) creates the characteristic durability of species like Pinus sylvestris and P. palustris (longleaf pine).

Python: Terpenoid Emission Modeling (Guenther Algorithm)

This simulation implements the full Guenther algorithm for isoprene (temperature + light response) and the beta-model for monoterpene emissions. It shows: (1) temperature response for both compound classes at standard PAR, (2) isoprene light response at 3 temperatures, and (3) the terpenoid chemical diversity by carbon class, illustrating why sesquiterpenes (C15) dominate known terpenoid diversity.

Python
script.py116 lines

Click Run to execute the Python code

Code will be executed with Python 3 on the server