Part V — Chapter 14

Alkaloid Biosynthesis

Nitrogen-containing secondary metabolites derived from amino acid precursors — from the opium poppy’s morphine to the coffee bean’s caffeine, alkaloids represent some of the most pharmacologically potent natural products.

Alkaloid Classes from Amino Acid Precursors

Tryptophan-Derived: Indole Alkaloids

The monoterpene indole alkaloids (MIAs) are the largest and most structurally diverse group, with over 3,000 known compounds. They are produced by the Catharanthus roseus (periwinkle),Strychnos, Rauvolfia, and related plants.

Strictosidine — the Universal MIA Intermediate

Strictosidine synthase (STR) catalyses the Pictet-Speraserase condensation of tryptamine (from Trp decarboxylation) with secologanin (a secoiridoid terpene from the MEP pathway) to form strictosidine — the key intermediate for virtually all MIAs.

\( \text{tryptamine} + \text{secologanin} \xrightarrow{\text{STR}} \text{strictosidine} + \text{H}_2\text{O} \)

Pathway to Vincristine (Anti-cancer)

Strictosidine → catharanthine + tabersonine (via ~10 enzymatic steps) → vindoline (formed in vacuoles via AVLBS, T16H, etc.) → vinblastine (coupling by peroxidase in vacuole) → vincristine (N-oxidation). The pathway spans 4 subcellular compartments (cytosol, ER, vacuole, chloroplast).

Tyrosine-Derived: Benzylisoquinoline Alkaloids (BIAs)

BIAs are produced by the Papaveraceae (opium poppy), Berberidaceae, and related families. The pathway begins with dopamine and 4-hydroxyphenylacetaldehyde (both from Tyr) condensing to form (S)-norcoclaurine, the key intermediate.

Norcoclaurine Synthase (NCS)

Pictet-Speraserase condensation of dopamine + 4-HPAA → (S)-norcoclaurine. Committed step to all BIAs.

\( \text{dopamine} + \text{4-HPAA} \xrightarrow{\text{NCS}} \text{(S)-norcoclaurine} \)

Morphine Branch

Norcoclaurine → reticuline (6-OMT, CNMT, CYP80B1) → salutaridine (SalSyn CYP450) → thebaine → codeinone → codeine → morphine (T6ODM, COR, CODM).

~17 enzymatic steps from Tyr to morphine.

Berberine Branch

Reticuline → scoulerine (SMT) → tetrahydrocolumbamine → canadine → berberine (BBE, ChBER). Yellow antimicrobial pigment; model for metabolic engineering.

Ornithine-Derived: Tropane Alkaloids

Tropane alkaloids are produced primarily by Solanaceae (Atropa belladonna,Datura, Hyoscyamus, Erythroxylum coca). The bicyclic tropane ring is assembled from ornithine via putrescine, N-methylputrescine, 4-methylaminobutanal, and the N-methyl-D1-pyrrolinium cation.

Key biosynthetic step: Tropinone formation

\( \text{N-methyl-}\Delta^1\text{-pyrrolinium}^+ + \text{acetyl-CoA} \xrightarrow{\text{PYKS/CS}} \text{tropinone} \)

Two tropinone reductases (TRI, TRII) determine the stereospecific branch: TRI makes tropine (for atropine/scopolamine); TRII makes pseudotropine.

Atropine (muscarinic antagonist)

Tropine + tropic acid (from phenylalanine) → hyoscyamine (lit by hyoscyamine 6beta-hydroxylase, H6H) → scopolamine. Blocks acetylcholine at muscarinic receptors; pupil dilation, decreased secretions.

Cocaine (stimulant)

Produced by Erythroxylum coca. Pseudotropine branch: acetyl-CoA condensation, ester formation with methylecgonine, then cocaine synthase.

Purine Alkaloids — Caffeine

Caffeine is biosynthesised from xanthosine (released from nucleotide catabolism) by 3 sequential N-methylations using SAM, catalysed by two key enzymes: XMT and DXMT.

\( \text{xanthosine} \xrightarrow{\text{XMT}} \text{7-MX} \xrightarrow{\text{MXMT}} \text{theobromine} \xrightarrow{\text{DXMT}} \text{caffeine} \)

Caffeine evolved independently 5+ times in plant evolution (coffee, tea, cacao, guarana, Ilex). Functions as a feeding deterrent and allelopathic agent.

Python: Caffeine Accumulation in Coffee Bean Development

Model caffeine accumulation kinetics using Michaelis-Menten equations for each methyltransferase step (XMT, MXMT, DXMT) during coffee bean maturation.

Python

script.py128 lines

import numpy as np
import matplotlib
matplotlib.use('Agg')
import matplotlib.pyplot as plt

# Model caffeine accumulation in coffee bean development
# Caffeine biosynthesis from xanthosine: 4 methylation steps
# Xanthosine -> 7-methylxanthosine -> 7-methylxanthine -> theobromine -> caffeine
# Each step uses S-adenosylmethionine (SAM) as methyl donor

# Enzyme kinetics: Michaelis-Menten for each methyltransferase step
Km_values = {
    'XMT (xanthosine->7-methylxanthosine)': 0.15,   # mM
    'MXMT (7-methylxanthine->theobromine)': 0.08,
    'DXMT (theobromine->caffeine)': 0.12,
}
Vmax_values = {
    'XMT (xanthosine->7-methylxanthosine)': 2.5,    # nmol/min/mg protein
    'MXMT (7-methylxanthine->theobromine)': 1.8,
    'DXMT (theobromine->caffeine)': 3.2,
}

# Coffee bean development stages (weeks after flowering)
weeks = np.linspace(0, 32, 200)

# Xanthosine pool: rises then falls as precursors are committed
xanthosine = 2.0 * np.exp(-(weeks - 8)**2 / 30) + 0.3 * np.exp(-weeks / 20)

# Enzyme activity profiles (increase then plateau as bean matures)
def enzyme_activity(weeks, peak_week, width, max_act):
    return max_act / (1 + np.exp(-(weeks - peak_week) / width))

xmt_activity = enzyme_activity(weeks, 6, 2, 1.0)
mxmt_activity = enzyme_activity(weeks, 10, 3, 0.85)
dxmt_activity = enzyme_activity(weeks, 14, 3, 1.2)

# Caffeine accumulation model (simplified ODE integration)
dt = weeks[1] - weeks[0]
caffeine = np.zeros(len(weeks))
theobromine = np.zeros(len(weeks))
intermediate = np.zeros(len(weeks))  # 7-methylxanthine

caffeine[0] = 0.0
theobromine[0] = 0.0
intermediate[0] = 0.0

for i in range(1, len(weeks)):
    sub_xmt = max(xanthosine[i], 0)
    sub_mxmt = max(intermediate[i-1], 0)
    sub_dxmt = max(theobromine[i-1], 0)

Km_xmt = Km_values['XMT (xanthosine->7-methylxanthosine)']
    Km_mxmt = Km_values['MXMT (7-methylxanthine->theobromine)']
    Km_dxmt = Km_values['DXMT (theobromine->caffeine)']

v_xmt = Vmax_values['XMT (xanthosine->7-methylxanthosine)'] * xmt_activity[i] * sub_xmt / (Km_xmt + sub_xmt)
    v_mxmt = Vmax_values['MXMT (7-methylxanthine->theobromine)'] * mxmt_activity[i] * sub_mxmt / (Km_mxmt + sub_mxmt)
    v_dxmt = Vmax_values['DXMT (theobromine->caffeine)'] * dxmt_activity[i] * sub_dxmt / (Km_dxmt + sub_dxmt)

intermediate[i] = intermediate[i-1] + (v_xmt - v_mxmt) * dt * 0.1
    theobromine[i] = theobromine[i-1] + (v_mxmt - v_dxmt) * dt * 0.1
    caffeine[i] = caffeine[i-1] + v_dxmt * dt * 0.1

# Clamp to non-negative
    intermediate[i] = max(0, intermediate[i])
    theobromine[i] = max(0, theobromine[i])

fig, axes = plt.subplots(2, 2, figsize=(13, 9))
fig.patch.set_facecolor('#0f172a')
for ax in axes.flat:
    ax.set_facecolor('#1e293b')
    ax.tick_params(colors='#94a3b8')
    ax.spines['bottom'].set_color('#475569')
    ax.spines['left'].set_color('#475569')
    ax.spines['top'].set_visible(False)
    ax.spines['right'].set_visible(False)

# Plot 1: Metabolite accumulation
axes[0, 0].plot(weeks, caffeine / (caffeine.max() + 0.001), color='#a3e635', linewidth=2.5, label='Caffeine')
axes[0, 0].plot(weeks, theobromine / (theobromine.max() + 0.001), color='#fb923c', linewidth=2, linestyle='--', label='Theobromine')
axes[0, 0].plot(weeks, intermediate / (intermediate.max() + 0.001), color='#60a5fa', linewidth=2, linestyle=':', label='7-methylxanthine')
axes[0, 0].plot(weeks, xanthosine / xanthosine.max(), color='#94a3b8', linewidth=1.5, linestyle='-.', label='Xanthosine (precursor)')
axes[0, 0].set_xlabel('Weeks after flowering', color='#94a3b8', fontsize=9)
axes[0, 0].set_ylabel('Relative concentration', color='#94a3b8', fontsize=9)
axes[0, 0].set_title('Caffeine Pathway Metabolites During Bean Development', color='#e2e8f0', fontsize=10)
axes[0, 0].legend(fontsize=7.5, facecolor='#1e293b', edgecolor='#475569', labelcolor='#e2e8f0')

# Plot 2: Enzyme activity profiles
axes[0, 1].plot(weeks, xmt_activity, color='#a3e635', linewidth=2, label='XMT (step 1)')
axes[0, 1].plot(weeks, mxmt_activity, color='#f59e0b', linewidth=2, label='MXMT (step 2)')
axes[0, 1].plot(weeks, dxmt_activity, color='#ec4899', linewidth=2, label='DXMT (step 3)')
axes[0, 1].set_xlabel('Weeks after flowering', color='#94a3b8', fontsize=9)
axes[0, 1].set_ylabel('Relative enzyme activity', color='#94a3b8', fontsize=9)
axes[0, 1].set_title('Methyltransferase Activity Profiles', color='#e2e8f0', fontsize=10)
axes[0, 1].legend(fontsize=8, facecolor='#1e293b', edgecolor='#475569', labelcolor='#e2e8f0')

# Plot 3: Michaelis-Menten curves for each enzyme
substrate_conc = np.linspace(0, 1.0, 300)
for (name, Km), color in zip(Km_values.items(), ['#a3e635', '#f59e0b', '#ec4899']):
    Vmax = Vmax_values[name]
    v = Vmax * substrate_conc / (Km + substrate_conc)
    short = name.split('(')[0].strip()
    axes[1, 0].plot(substrate_conc, v, color=color, linewidth=2, label=f'{short} (Km={Km} mM)')
    axes[1, 0].axvline(x=Km, color=color, alpha=0.3, linestyle=':')
axes[1, 0].set_xlabel('Substrate [mM]', color='#94a3b8', fontsize=9)
axes[1, 0].set_ylabel('Velocity (nmol/min/mg)', color='#94a3b8', fontsize=9)
axes[1, 0].set_title('Michaelis-Menten Kinetics of Caffeine Methyltransferases', color='#e2e8f0', fontsize=10)
axes[1, 0].legend(fontsize=7.5, facecolor='#1e293b', edgecolor='#475569', labelcolor='#e2e8f0')

# Plot 4: Alkaloid diversity bar chart
alkaloid_classes = ['Indole alkaloids', 'Benzylisoquinoline', 'Tropane', 'Purine alkaloids', 'Piperidine', 'Quinolizidine']
compound_counts = [3800, 2500, 200, 150, 700, 170]
colors4 = ['#a3e635', '#34d399', '#22d3ee', '#60a5fa', '#a78bfa', '#f472b6']
bars = axes[1, 1].barh(alkaloid_classes, compound_counts, color=colors4, alpha=0.85, edgecolor='#475569')
axes[1, 1].set_xlabel('Estimated compound count', color='#94a3b8', fontsize=9)
axes[1, 1].set_title('Alkaloid Class Diversity', color='#e2e8f0', fontsize=10)
for bar, val in zip(bars, compound_counts):
    axes[1, 1].text(bar.get_width() + 30, bar.get_y() + bar.get_height()/2,
                    str(val), va='center', color='#e2e8f0', fontsize=8)

fig.suptitle('Alkaloid Biosynthesis: Caffeine Pathway Modeling', color='#e2e8f0',
             fontsize=13, fontweight='bold')
plt.tight_layout()
plt.savefig('output.png', dpi=110, bbox_inches='tight',
            facecolor='#0f172a', edgecolor='none')
print("Peak caffeine accumulation at week:", round(weeks[np.argmax(caffeine)], 1))
print("Theobromine:caffeine final ratio:", round(theobromine[-1] / (caffeine[-1] + 0.001), 3))

Click Run to execute the Python code

Code will be executed with Python 3 on the server

← Ch 13: Terpenoid Biosynthesis Ch 15: Flavonoids & Anthocyanins →