Module 3

Venom Biochemistry

Ion channel toxins, enzyme kinetics, necrotic venoms, and venom as drug leads

Featured Lecture

The Story of Spiders and their Venoms — Dr. Michel Dugon

3.1 Venom Composition

Spider venom is among the most biochemically complex secretions in the animal kingdom. A single species can produce 100 to 1,000+ distinct peptides, each evolved to target specific physiological systems in prey. The venom cocktail is synthesised in paired venom glands located in the chelicerae (or prosoma, in some species) and delivered through hollow or grooved fangs.

The major biochemical classes include:

Neurotoxins

Ion channel modulators (Na$^+$, Ca$^{2+}$, K$^+$ channels). Typically small disulfide-rich peptides (3–8 kDa) with an inhibitor cystine knot (ICK) motif that confers remarkable stability.

Cytotoxins

Membrane-disrupting peptides and phospholipases that lyse cell membranes. Include sphingomyelinase D (brown recluse) and various antimicrobial peptides.

Enzymes

Hyaluronidase (“spreading factor” that degrades connective tissue), phospholipase A$_2$, metalloproteases, and serine proteases that facilitate venom diffusion and tissue destruction.

Dose-Response: The Weibull LD$_{50}$ Model

The toxicity of a venom is characterised by its LD$_{50}$ — the dose lethal to 50% of test animals. The dose-response relationship is well described by a Weibull cumulative distribution function:

$ P(\text{death} | D) = 1 - \exp\!\left[-\left(\frac{D}{\lambda}\right)^k\right] $

• $D$ is the administered dose (mg/kg)
• $\lambda$ is the scale parameter (related to LD$_{50}$)
• $k$ is the shape parameter (steepness of the dose-response curve)

At $P = 0.5$, we solve for the LD$_{50}$:

$ \text{LD}_{50} = \lambda \cdot (\ln 2)^{1/k} $

For spider venoms, LD$_{50}$ values range enormously: from $\sim 0.001$ mg/kg (Sydney funnel-web Atrax robustus, among the most toxic) to $> 100$ mg/kg for many harmless house spiders.

3.2 Ion Channel Toxins

The majority of spider neurotoxins target voltage-gated ion channels in the nervous system of prey. Three exemplary toxins illustrate the diversity:

$\omega$-Agatoxin IVA (Ca$^{2+}$ Channel Blocker)

From the funnel-web spider Agelenopsis aperta. Blocks P/Q-type Ca$^{2+}$ channels at presynaptic nerve terminals, preventing neurotransmitter release. A 48-residue peptide with 4 disulfide bridges. IC$_{50} \approx 1{-}3$ nM — extraordinary potency.

PhTx3 (Na$^+$ Channel from Phoneutria)

The Brazilian wandering spider (Phoneutria nigriventer) produces PhTx3, which delays Na$^+$ channel inactivation, causing persistent depolarisation. In humans, this results in intense pain, salivation, and notably, priapism (persistent erection) — leading to pharmacological interest in the PnTx2-6 subfraction as a potential erectile dysfunction treatment.

$\delta$-Atracotoxin (Na$^+$ Channel from Atrax)

From the Sydney funnel-web spider (Atrax robustus). Slows Na$^+$ channel inactivation, producing massive, sustained neuronal firing. Lethal to primates (including humans) but not to many other mammals — a peculiar selectivity. LD$_{50} \approx 0.16$ mg/kg in mice.

Modified Hodgkin-Huxley Model

The Hodgkin-Huxley model describes action potential generation. For Na$^+$ current:

$ I_{\text{Na}} = g_{\text{Na}} \cdot m^3 \cdot h \cdot (V - E_{\text{Na}}) $

where $m$ is the activation gate, $h$ is the inactivation gate,$g_{\text{Na}}$ is the maximum conductance, and $E_{\text{Na}}$ is the Nernst reversal potential. Toxins like $\delta$-atracotoxin modify the inactivation kinetics by shifting the steady-state inactivation curve:

$ h_\infty(V) = \frac{1}{1 + \exp\!\left(\frac{V - V_{1/2} - \Delta V_{\text{toxin}}}{k_h}\right)} $

The toxin shifts $V_{1/2}$ by $\Delta V_{\text{toxin}}$ (typically +10 to +30 mV), preventing complete inactivation and causing prolonged, repetitive firing.

Figure 3.1: Voltage-gated Na$^+$ channel structure (domains I–IV) with spider toxin binding sites. $\delta$-Atracotoxin binds at Site 3 (DI–DII linker), PhTx3 at Site 4 (DIV S3–S4), both preventing channel inactivation.

3.3 Sphingomyelinase D: The Brown Recluse Mechanism

The brown recluse spider (Loxosceles reclusa) and its relatives produce a unique enzyme — sphingomyelinase D (SMase D) — responsible for the characteristic dermonecrotic lesions. This enzyme is rare in nature, found only in Loxosceles spiders and certain pathogenic bacteria (Corynebacterium).

The enzymatic cascade:

Sphingomyelin $\xrightarrow{\text{SMase D}}$ Ceramide-1-phosphate + Choline

Ceramide-1-phosphate $\rightarrow$ Complement activation (C5a, C3a)

$\rightarrow$ Neutrophil infiltration $\rightarrow$ Endothelial damage

$\rightarrow$ Dermonecrosis (tissue death spreading over days)

Michaelis-Menten Kinetics for SMase D

The enzymatic reaction follows standard Michaelis-Menten kinetics:

$ v = \frac{V_{\max} \cdot [S]}{K_m + [S]} = \frac{k_{\text{cat}} \cdot [E]_0 \cdot [S]}{K_m + [S]} $

• $V_{\max} = k_{\text{cat}} \cdot [E]_0$: maximum reaction rate
• $K_m \approx 50{-}200 \;\mu\text{M}$: Michaelis constant for sphingomyelin substrate
• $k_{\text{cat}} \approx 200{-}500 \;\text{s}^{-1}$: catalytic turnover number

The catalytic efficiency $k_{\text{cat}}/K_m$ for Loxosceles SMase D is remarkably high ($\sim 10^6 \;\text{M}^{-1}\text{s}^{-1}$), approaching diffusion-limited rates — meaning even trace amounts of injected venom produce significant tissue damage.

The time-dependent product formation during a venom envenomation:

$ [P](t) = V_{\max} \cdot t - K_m \cdot \ln\!\left(\frac{[S]_0}{[S]_0 - [P](t)}\right) $

3.4 Venom-Derived Pharmaceuticals

Spider venoms represent a vast, largely untapped pharmacological library. With an estimated 10 million+ bioactive peptides across all spider species (fewer than 0.1% characterised), venoms are a rich source of drug leads.

$\omega$-Conotoxin MVIIA (Prialt™)

Although from cone snails (not spiders), this demonstrates the venom-to-drug pipeline. Blocks N-type Ca$^{2+}$ channels. FDA-approved for intractable pain. Spider $\omega$-agatoxins target homologous P/Q-type channels.

Huwentoxin-IV (Analgesic Candidate)

From the Chinese bird spider (Cyriopagopus schmidti). Selectively blocks Na$_v$1.7, a pain-specific sodium channel. Currently in preclinical development as a non-opioid analgesic.

Spider Peptides as Bio-Insecticides

Many spider toxins show remarkable selectivity for insect over mammalian channels. This selectivity arises from subtle structural differences between insect and mammalian channel isoforms, making them ideal insecticide leads.

Selectivity Index

The selectivity index (SI) quantifies how preferentially a toxin targets insect versus mammalian channels:

$ \text{SI} = \frac{\text{IC}_{50}(\text{mammal})}{\text{IC}_{50}(\text{insect})} $

An ideal bio-insecticide has SI $\gg 1$. For example:

• $\omega$-ACTX-Hv1a from Hadronyche versuta: SI $\approx 100{-}1000$
• U$_1$-TRTX-Sp1a: SI $\approx 500$
• Conventional organophosphate insecticides: SI $\approx 1{-}10$ (toxic to mammals too)

3.5 Venom Delivery: Fang Mechanics

Spider fangs (cheliceral teeth) come in two fundamental arrangements:

Araneomorph (Cross-Fang)

Fangs pivot laterally, crossing like pincers. Found in >90% of spider species (the “true spiders”). Allows fine control and smaller prey handling.

Mygalomorph (Parallel Fang)

Fangs strike downward in parallel, like pickaxes. Found in tarantulas, funnel-webs, trapdoor spiders. Requires rearing up to strike. Greater penetration force.

Penetration Force Model

The force required for fang penetration depends on the fang geometry. For a conical fang with tip radius $r_t$, base radius $r_b$, and length $L$:

$ F_{\text{penetration}} = \sigma_y \cdot \pi r_t^2 + 2\pi \mu \int_0^d r(x) \cdot p(x) \, dx $

where $\sigma_y$ is the yield stress of the prey cuticle,$\mu$ is the friction coefficient, $d$ is the penetration depth, and$r(x) = r_t + (r_b - r_t) x/L$ is the local fang radius. The pressure distribution$p(x)$ depends on the elastic properties of the cuticle.

For a sharp fang tip ($r_t \approx 1{-}5\;\mu\text{m}$), the initial puncture force is remarkably small: $\sim 1{-}10$ mN. The fang tip achieves stress concentration factors of $\sigma/\sigma_{\text{far}} \sim 100{-}1000\times$, far exceeding the cuticle yield stress.

Figure 3.2: Venom gland anatomy (left) showing secretory epithelium, compressor muscles, and venom duct leading to the hollow fang. Right: comparison of araneomorph (cross-fang) vs mygalomorph (parallel-fang) strike mechanics.

3.6 Computational Analysis

Dose-Response Curves Across Spider Species

Weibull Dose-Response: LD50 Comparison

Python

script.py70 lines

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

def weibull_response(D, lam, k):
    """Weibull dose-response: P(death | dose D)."""
    return 1.0 - np.exp(-((D / lam) ** k))

def ld50_from_params(lam, k):
    """Calculate LD50 from Weibull parameters."""
    return lam * (np.log(2)) ** (1.0 / k)

# Spider venom parameters: (lambda, k, LD50_literature mg/kg)
species = {
    'Atrax robustus\n(Sydney funnel-web)': {'lam': 0.18, 'k': 2.5, 'color': '#ef4444'},
    'Phoneutria nigriventer\n(Brazilian wandering)': {'lam': 0.20, 'k': 2.2, 'color': '#f97316'},
    'Latrodectus mactans\n(Black widow)': {'lam': 0.95, 'k': 2.8, 'color': '#a855f7'},
    'Loxosceles reclusa\n(Brown recluse)': {'lam': 5.5, 'k': 1.8, 'color': '#eab308'},
    'Argiope bruennichii\n(Wasp spider)': {'lam': 80.0, 'k': 2.0, 'color': '#22d3ee'},
}

doses = np.linspace(0.001, 150, 2000)

fig, (ax1, ax2) = plt.subplots(1, 2, figsize=(14, 6))
fig.patch.set_facecolor('#0a0a0a')
for ax in [ax1, ax2]:
    ax.set_facecolor('#111111')
    ax.tick_params(colors='#cccccc')
    for spine in ax.spines.values():
        spine.set_color('#444444')

# Linear scale
for name, params in species.items():
    response = weibull_response(doses, params['lam'], params['k'])
    ld50 = ld50_from_params(params['lam'], params['k'])
    ax1.plot(doses, response * 100, color=params['color'], linewidth=2,
             label=f"{name}\nLD50={ld50:.2f} mg/kg")
    ax1.axhline(y=50, color='#666666', linestyle='--', linewidth=0.8)

ax1.set_xlabel('Dose (mg/kg)', color='#cccccc', fontsize=12)
ax1.set_ylabel('Mortality (%)', color='#cccccc', fontsize=12)
ax1.set_title('Dose-Response: Linear Scale', color='#c4b5fd', fontsize=14, fontweight='bold')
ax1.set_xlim(0, 15)
ax1.set_ylim(0, 105)
ax1.legend(fontsize=8, loc='lower right', facecolor='#1a1a2e', edgecolor='#555555', labelcolor='#cccccc')

# Log scale
doses_log = np.logspace(-3, 2.5, 2000)
for name, params in species.items():
    response = weibull_response(doses_log, params['lam'], params['k'])
    ld50 = ld50_from_params(params['lam'], params['k'])
    ax2.semilogx(doses_log, response * 100, color=params['color'], linewidth=2,
                 label=name.split('\n')[0])
    ax2.axvline(x=ld50, color=params['color'], linestyle=':', linewidth=1, alpha=0.5)

ax2.axhline(y=50, color='#666666', linestyle='--', linewidth=0.8)
ax2.set_xlabel('Dose (mg/kg, log scale)', color='#cccccc', fontsize=12)
ax2.set_ylabel('Mortality (%)', color='#cccccc', fontsize=12)
ax2.set_title('Dose-Response: Log Scale', color='#c4b5fd', fontsize=14, fontweight='bold')
ax2.set_ylim(0, 105)
ax2.legend(fontsize=9, loc='lower right', facecolor='#1a1a2e', edgecolor='#555555', labelcolor='#cccccc')

plt.tight_layout()
plt.savefig('output.png', dpi=150, bbox_inches='tight', facecolor='#0a0a0a')
print("LD50 values (Weibull model):")
for name, params in species.items():
    ld50 = ld50_from_params(params['lam'], params['k'])
    print(f"  {name.replace(chr(10), ' ')}: LD50 = {ld50:.3f} mg/kg")
print("\nNote: Lower LD50 = more toxic. Atrax is ~500x more toxic than Argiope.")

Click Run to execute the Python code

Code will be executed with Python 3 on the server

Hodgkin-Huxley with Toxin Modification

Action Potential: Normal vs Toxin-Modified Na+ Channel

Python

script.py125 lines

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

# Hodgkin-Huxley parameters
C_m = 1.0     # membrane capacitance (uF/cm2)
g_Na = 120.0  # max Na conductance (mS/cm2)
g_K = 36.0    # max K conductance (mS/cm2)
g_L = 0.3     # leak conductance (mS/cm2)
E_Na = 50.0   # Na reversal potential (mV)
E_K = -77.0   # K reversal potential (mV)
E_L = -54.4   # leak reversal potential (mV)

def alpha_m(V): return 0.1 * (V + 40) / (1 - np.exp(-(V + 40) / 10))
def beta_m(V): return 4.0 * np.exp(-(V + 65) / 18)
def alpha_h(V): return 0.07 * np.exp(-(V + 65) / 20)
def beta_h(V): return 1.0 / (1 + np.exp(-(V + 35) / 10))
def alpha_n(V): return 0.01 * (V + 55) / (1 - np.exp(-(V + 55) / 10))
def beta_n(V): return 0.125 * np.exp(-(V + 65) / 80)

# Toxin-modified h-gate: shift V_half of inactivation
def alpha_h_toxin(V, dV):
    return 0.07 * np.exp(-(V + 65 - dV) / 20)
def beta_h_toxin(V, dV):
    return 1.0 / (1 + np.exp(-(V + 35 - dV) / 10))

def simulate_HH(dt, T_max, I_ext, toxin_dV=0.0):
    steps = int(T_max / dt)
    V = np.zeros(steps); m = np.zeros(steps)
    h = np.zeros(steps); n = np.zeros(steps)
    t = np.linspace(0, T_max, steps)
    V[0] = -65.0; m[0] = 0.05; h[0] = 0.6; n[0] = 0.32

for i in range(1, steps):
        am = alpha_m(V[i-1]); bm = beta_m(V[i-1])
        if toxin_dV == 0:
            ah = alpha_h(V[i-1]); bh = beta_h(V[i-1])
        else:
            ah = alpha_h_toxin(V[i-1], toxin_dV)
            bh = beta_h_toxin(V[i-1], toxin_dV)
        an = alpha_n(V[i-1]); bn = beta_n(V[i-1])

m[i] = m[i-1] + dt * (am * (1 - m[i-1]) - bm * m[i-1])
        h[i] = h[i-1] + dt * (ah * (1 - h[i-1]) - bh * h[i-1])
        n[i] = n[i-1] + dt * (an * (1 - n[i-1]) - bn * n[i-1])

I_Na = g_Na * m[i]**3 * h[i] * (V[i-1] - E_Na)
        I_K = g_K * n[i]**4 * (V[i-1] - E_K)
        I_L = g_L * (V[i-1] - E_L)

V[i] = V[i-1] + dt * (I_ext - I_Na - I_K - I_L) / C_m
    return t, V, m, h

dt = 0.01; T_max = 80.0
I_ext = 10.0

t_norm, V_norm, m_norm, h_norm = simulate_HH(dt, T_max, I_ext, toxin_dV=0.0)
t_mild, V_mild, m_mild, h_mild = simulate_HH(dt, T_max, I_ext, toxin_dV=15.0)
t_severe, V_severe, m_severe, h_severe = simulate_HH(dt, T_max, I_ext, toxin_dV=25.0)

fig, axes = plt.subplots(2, 2, figsize=(14, 10))
fig.patch.set_facecolor('#0a0a0a')
for ax in axes.flat:
    ax.set_facecolor('#111111')
    ax.tick_params(colors='#cccccc')
    for spine in ax.spines.values():
        spine.set_color('#444444')

# Action potentials
axes[0,0].plot(t_norm, V_norm, color='#22d3ee', linewidth=1.5, label='Normal')
axes[0,0].plot(t_mild, V_mild, color='#f97316', linewidth=1.5, label='Mild toxin (+15 mV shift)')
axes[0,0].plot(t_severe, V_severe, color='#ef4444', linewidth=1.5, label='Severe toxin (+25 mV shift)')
axes[0,0].set_xlabel('Time (ms)', color='#cccccc')
axes[0,0].set_ylabel('Membrane potential (mV)', color='#cccccc')
axes[0,0].set_title('Action Potentials', color='#c4b5fd', fontsize=13, fontweight='bold')
axes[0,0].legend(fontsize=9, facecolor='#1a1a2e', edgecolor='#555', labelcolor='#cccccc')

# h-gate dynamics
axes[0,1].plot(t_norm, h_norm, color='#22d3ee', linewidth=1.5, label='Normal h')
axes[0,1].plot(t_mild, h_mild, color='#f97316', linewidth=1.5, label='Toxin h (+15 mV)')
axes[0,1].plot(t_severe, h_severe, color='#ef4444', linewidth=1.5, label='Toxin h (+25 mV)')
axes[0,1].set_xlabel('Time (ms)', color='#cccccc')
axes[0,1].set_ylabel('h (inactivation gate)', color='#cccccc')
axes[0,1].set_title('Inactivation Gate Dynamics', color='#c4b5fd', fontsize=13, fontweight='bold')
axes[0,1].legend(fontsize=9, facecolor='#1a1a2e', edgecolor='#555', labelcolor='#cccccc')

# Steady-state inactivation curves
V_range = np.linspace(-100, 20, 500)
h_inf_normal = alpha_h(V_range) / (alpha_h(V_range) + beta_h(V_range))
h_inf_mild = alpha_h_toxin(V_range, 15) / (alpha_h_toxin(V_range, 15) + beta_h_toxin(V_range, 15))
h_inf_severe = alpha_h_toxin(V_range, 25) / (alpha_h_toxin(V_range, 25) + beta_h_toxin(V_range, 25))

axes[1,0].plot(V_range, h_inf_normal, color='#22d3ee', linewidth=2, label='Normal')
axes[1,0].plot(V_range, h_inf_mild, color='#f97316', linewidth=2, label='+15 mV shift')
axes[1,0].plot(V_range, h_inf_severe, color='#ef4444', linewidth=2, label='+25 mV shift')
axes[1,0].axhline(y=0.5, color='#666', linestyle='--', linewidth=0.8)
axes[1,0].set_xlabel('Membrane potential (mV)', color='#cccccc')
axes[1,0].set_ylabel('h_inf (steady-state)', color='#cccccc')
axes[1,0].set_title('Steady-State Inactivation Curve', color='#c4b5fd', fontsize=13, fontweight='bold')
axes[1,0].legend(fontsize=9, facecolor='#1a1a2e', edgecolor='#555', labelcolor='#cccccc')

# Venom peptide mass distribution
np.random.seed(42)
masses_funnel = np.random.lognormal(mean=np.log(5.0), sigma=0.4, size=300)
masses_widow = np.random.lognormal(mean=np.log(7.0), sigma=0.5, size=200)
masses_recluse = np.random.lognormal(mean=np.log(32.0), sigma=0.3, size=100)

axes[1,1].hist(masses_funnel, bins=40, alpha=0.6, color='#ef4444', label='Funnel-web (neurotoxins)', density=True)
axes[1,1].hist(masses_widow, bins=40, alpha=0.6, color='#a855f7', label='Black widow (latrotoxins)', density=True)
axes[1,1].hist(masses_recluse, bins=30, alpha=0.6, color='#eab308', label='Brown recluse (enzymes)', density=True)
axes[1,1].set_xlabel('Peptide mass (kDa)', color='#cccccc')
axes[1,1].set_ylabel('Density', color='#cccccc')
axes[1,1].set_title('Venom Peptide Mass Distribution', color='#c4b5fd', fontsize=13, fontweight='bold')
axes[1,1].legend(fontsize=9, facecolor='#1a1a2e', edgecolor='#555', labelcolor='#cccccc')
axes[1,1].set_xlim(0, 60)

plt.tight_layout()
plt.savefig('output.png', dpi=150, bbox_inches='tight', facecolor='#0a0a0a')
print("Hodgkin-Huxley toxin simulation results:")
print(f"  Normal: {np.sum(np.diff(np.sign(np.diff(V_norm))) < -1)} action potentials in {T_max} ms")
print(f"  Mild toxin (+15 mV): {np.sum(np.diff(np.sign(np.diff(V_mild))) < -1)} action potentials")
print(f"  Severe toxin (+25 mV): {np.sum(np.diff(np.sign(np.diff(V_severe))) < -1)} action potentials")
print("\nToxin shifts inactivation curve rightward, preventing channel closure.")
print("Result: sustained, repetitive firing -> excitotoxicity -> death.")

Click Run to execute the Python code

Code will be executed with Python 3 on the server

3.7 Advanced Venom Biochemistry

3.7.1 Inhibitor Cystine Knot (ICK) Motif

The inhibitor cystine knot (ICK) is the dominant structural scaffold in spider venom peptides, found in over 1,000 characterized toxins. The ICK motif consists of three disulfide bonds arranged in a distinctive knotted topology:

$\text{Cys}_1{-}\text{Cys}_4, \quad \text{Cys}_2{-}\text{Cys}_5, \quad \text{Cys}_3{-}\text{Cys}_6$

The third disulfide bond (Cys$_3$-Cys$_6$) threads through the macrocyclic ring formed by bonds 1 and 2 and the intervening backbone, creating a true topological knot that cannot be undone without breaking covalent bonds.

This knotted topology confers extraordinary stability. The thermodynamic stability of the ICK fold can be quantified by the unfolding free energy:

$\Delta G_{\text{unfolding}} = \Delta G_{\text{disulfides}} + \Delta G_{\text{backbone}} + \Delta G_{\text{knot}} > 40 \;\text{kJ/mol}$

• $\Delta G_{\text{disulfides}} \approx 3 \times 12{-}15 \;\text{kJ/mol} = 36{-}45 \;\text{kJ/mol}$ from three S-S bonds
• $\Delta G_{\text{backbone}} \approx 5{-}15 \;\text{kJ/mol}$ from hydrogen bonds and hydrophobic packing
• $\Delta G_{\text{knot}} \approx 10{-}20 \;\text{kJ/mol}$ additional kinetic barrier from the threaded topology

The biological consequences of ICK stability are profound:

• Protease resistance: ICK peptides resist degradation by trypsin, chymotrypsin, and pepsin — essential for surviving the proteolytic environment of prey tissues
• Thermal stability: ICK peptides retain function up to 80-100°C ($T_m > 90\;^\circ\text{C}$ for some toxins)
• Target specificity: The loop regions between cysteines form the pharmacophore that contacts ion channels with nanomolar affinity, while the knot core provides a rigid scaffold

3.7.2 Venom Proteomics & Transcriptomics

Modern venomics combines LC-MS/MS mass spectrometry with venom gland transcriptomics to characterize the full molecular complexity of spider venoms. A single species typically produces 200-1,000+ unique peptides and proteins, organized into toxin superfamilies.

The primary evolutionary mechanism driving venom diversity is gene duplication followed by neofunctionalization. After a toxin gene duplicates, one copy is free to accumulate mutations while the other maintains the original function. The rate of adaptive evolution is quantified by the Ka/Ks ratio (also written $dN/dS$ or $\omega$):

$\omega = \frac{K_a}{K_s} = \frac{d_N / d_S}{\text{(nonsynonymous substitutions)} / \text{(synonymous substitutions)}}$

• $\omega < 1$: Purifying selection — mutations are deleterious and removed (housekeeping genes: $\omega \approx 0.1{-}0.2$)
• $\omega = 1$: Neutral evolution — no selective pressure
• $\omega > 1$: Positive selection — mutations are advantageous and favored (toxin genes: $\omega \approx 2{-}5$)

Spider venom genes evolve 3-5× faster than average housekeeping genes, driven by the coevolutionary arms race between predator toxins and prey resistance. This is among the highest rates of adaptive molecular evolution observed in any animal protein family.

3.7.3 Synergistic Venom Action

Individual venom peptides are typically less toxic than the whole venom. This synergism is a key feature of spider venom pharmacology — the complex mixture is evolved to be more effective than the sum of its parts. Synergism is quantified using the Combination Index (CI) of Chou and Talalay:

$\text{CI} = \frac{D_1}{D_{x,1}} + \frac{D_2}{D_{x,2}}$

• $D_1, D_2$: doses of components 1 and 2 in the combination that achieve effect x
• $D_{x,1}, D_{x,2}$: doses of each component alone that achieve the same effect x
• CI < 1: Synergism — the combination requires less total drug
• CI = 1: Additive effect
• CI > 1: Antagonism

A canonical example of venom synergism involves three components acting in concert:

Neurotoxin

Blocks ion channels to paralyze prey. Requires access to nerve terminals, which is facilitated by the spreading factor.

Cytotoxin

Disrupts cell membranes, causing tissue damage that releases intracellular contents and amplifies the neurotoxin's access to targets.

Hyaluronidase

The “spreading factor” — degrades hyaluronic acid in the extracellular matrix, dramatically accelerating diffusion of neurotoxins and cytotoxins through tissue.

3.7.4 Antivenoms & Recombinant Toxins

Traditional antivenom production involves immunizing horses (or sheep) with whole venom and harvesting the resulting polyclonal IgG antibodies. While effective, this approach has significant limitations: batch variability, risk of serum sickness, and high cost ($50-$500+ per vial).

The next generation of antivenoms uses recombinant toxin fragments — specific toxin epitopes produced in E. coli or yeast expression systems — to generate targeted antibodies. The binding kinetics of antibody-toxin neutralization follow:

$K_d = \frac{k_{\text{off}}}{k_{\text{on}}} = \frac{[\text{Ab}][\text{Tx}]}{[\text{Ab:Tx}]}$

• $K_d$: dissociation constant (lower = tighter binding). Therapeutic antibodies need $K_d < 1\;\text{nM}$
• $k_{\text{on}}$: association rate constant, typically $\sim 10^5{-}10^6 \;\text{M}^{-1}\text{s}^{-1}$
• $k_{\text{off}}$: dissociation rate constant. For effective neutralization, $k_{\text{off}} < 10^{-3} \;\text{s}^{-1}$ (half-life > 10 min)

The neutralization efficiency depends on the antibody reaching sufficient concentration to sequester the toxin before it binds its target:

$\text{Neutralization}(\%) = \frac{[\text{Ab}]_0 / K_d}{1 + [\text{Ab}]_0 / K_d} \times 100$

Recombinant approaches enable cocktails of monoclonal antibodies targeting the most medically significant toxins, potentially replacing broad-spectrum horse-derived antivenoms with more precise, safer therapeutics.

3.8 Computational Analysis: Advanced Venom Biochemistry

Left: Free energy landscape for ICK fold unfolding compared to simpler disulfide architectures — the knotted topology creates an additional kinetic barrier of ~10-20 kJ/mol. Center: Ka/Ks ratio analysis showing positive selection ($\omega > 1$) in venom toxin genes vs purifying selection in housekeeping genes. Right: Combination Index analysis demonstrating synergistic interactions between venom components (CI < 1).

ICK Stability, Venom Proteomics & Synergism Analysis

Python

script.py138 lines

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

fig, axes = plt.subplots(1, 3, figsize=(17, 6))
fig.patch.set_facecolor('#030712')
fig.suptitle("Advanced Venom Biochemistry: ICK Stability, Proteomics & Synergism", fontsize=14, color='white', fontweight='bold', y=1.02)

# Panel 1: ICK fold stability - disulfide bond energetics
ax1 = axes[0]
ax1.set_facecolor('#1a0a2e')

# Free energy landscape for ICK unfolding
reaction_coord = np.linspace(0, 1, 500)

# Native ICK (3 disulfide bonds, knotted topology)
dG_ick = 45 * (4 * reaction_coord**2 - 4 * reaction_coord**3 + reaction_coord**4)
# Add kinetic barrier
dG_ick += 15 * np.exp(-((reaction_coord - 0.3) / 0.08)**2)

# Linear 3-disulfide (no knot)
dG_linear = 25 * (4 * reaction_coord**2 - 4 * reaction_coord**3 + reaction_coord**4)
dG_linear += 8 * np.exp(-((reaction_coord - 0.3) / 0.1)**2)

# Single disulfide peptide
dG_single = 12 * (4 * reaction_coord**2 - 4 * reaction_coord**3 + reaction_coord**4)

ax1.plot(reaction_coord, dG_ick, color='#c084fc', linewidth=2.5, label='ICK motif (knotted, 3 S-S)')
ax1.plot(reaction_coord, dG_linear, color='#22d3ee', linewidth=2.5, label='Linear 3-disulfide')
ax1.plot(reaction_coord, dG_single, color='#f87171', linewidth=2.5, linestyle='--', label='Single disulfide')
ax1.fill_between(reaction_coord, dG_ick, alpha=0.1, color='#c084fc')
ax1.axhline(y=40, color='#a78bfa', linestyle=':', linewidth=1, alpha=0.5)
ax1.annotate('ICK threshold\n(>40 kJ/mol)', xy=(0.7, 41), fontsize=9, color='#a78bfa', ha='center')
ax1.set_xlabel('Unfolding Reaction Coordinate', color='white', fontsize=11)
ax1.set_ylabel('Free Energy (kJ/mol)', color='white', fontsize=11)
ax1.set_title('ICK Fold Stability\nDisulfide Bond Energetics', color='#d8b4fe', fontsize=11)
ax1.tick_params(colors='white'); ax1.grid(True, alpha=0.2, color='#a78bfa')
for sp in ax1.spines.values(): sp.set_color('#7c3aed')
ax1.legend(fontsize=8, facecolor='#1a0a2e', edgecolor='#7c3aed', labelcolor='white', loc='upper left')

# Panel 2: Venom proteome diversity (Ka/Ks analysis)
ax2 = axes[1]
ax2.set_facecolor('#1a0a2e')

np.random.seed(42)
# Toxin genes: high Ka/Ks (positive selection)
n_toxin = 80
ka_ks_toxin = np.random.lognormal(mean=np.log(2.5), sigma=0.6, size=n_toxin)
ka_ks_toxin = np.clip(ka_ks_toxin, 0.3, 12)
diversity_toxin = np.random.uniform(0.1, 1.0, n_toxin)

# Housekeeping genes: low Ka/Ks (purifying selection)
n_house = 60
ka_ks_house = np.random.lognormal(mean=np.log(0.15), sigma=0.5, size=n_house)
ka_ks_house = np.clip(ka_ks_house, 0.01, 0.8)
diversity_house = np.random.uniform(0.1, 1.0, n_house)

ax2.scatter(diversity_toxin, ka_ks_toxin, c='#c084fc', s=40, alpha=0.7, edgecolors='white',
            linewidth=0.3, label='Venom toxin genes', zorder=5)
ax2.scatter(diversity_house, ka_ks_house, c='#34d399', s=40, alpha=0.7, edgecolors='white',
            linewidth=0.3, label='Housekeeping genes', zorder=5)
ax2.axhline(y=1.0, color='#fbbf24', linestyle='--', linewidth=1.5, alpha=0.8)
ax2.annotate('Ka/Ks = 1 (neutral)', xy=(0.5, 1.1), fontsize=9, color='#fbbf24', ha='center')
ax2.fill_between([0, 1.1], 1.0, 15, alpha=0.08, color='#ef4444')
ax2.fill_between([0, 1.1], 0, 1.0, alpha=0.08, color='#34d399')
ax2.text(0.95, 8, 'Positive\nselection', fontsize=9, color='#f87171', ha='right', va='center')
ax2.text(0.95, 0.3, 'Purifying\nselection', fontsize=9, color='#34d399', ha='right', va='center')
ax2.set_yscale('log')
ax2.set_xlabel('Gene Sequence Diversity', color='white', fontsize=11)
ax2.set_ylabel('Ka/Ks Ratio', color='white', fontsize=11)
ax2.set_title('Venom Gene Evolution\nKa/Ks Positive Selection', color='#d8b4fe', fontsize=11)
ax2.set_xlim(0, 1.05); ax2.set_ylim(0.01, 15)
ax2.tick_params(colors='white'); ax2.grid(True, alpha=0.2, color='#a78bfa')
for sp in ax2.spines.values(): sp.set_color('#7c3aed')
ax2.legend(fontsize=9, facecolor='#1a0a2e', edgecolor='#7c3aed', labelcolor='white', loc='upper left')

# Panel 3: Synergism - Combination Index (CI)
ax3 = axes[2]
ax3.set_facecolor('#1a0a2e')

# Dose fractions for two-component mixture
f1 = np.linspace(0.01, 0.99, 200)
f2 = 1 - f1

# Individual IC50 values
Dx1 = 10.0  # neurotoxin IC50 (ug/mL)
Dx2 = 25.0  # cytotoxin IC50 (ug/mL)

# Synergistic combination (CI < 1)
D1_syn = f1 * 5.0   # actual dose needed in combo
D2_syn = f2 * 8.0
CI_syn = D1_syn / Dx1 + D2_syn / Dx2

# Additive combination (CI = 1)
D1_add = f1 * Dx1 * 0.5
D2_add = f2 * Dx2 * 0.5
CI_add = D1_add / Dx1 + D2_add / Dx2

# Antagonistic (CI > 1)
D1_ant = f1 * 15.0
D2_ant = f2 * 35.0
CI_ant = D1_ant / Dx1 + D2_ant / Dx2

ax3.plot(f1, CI_syn, color='#c084fc', linewidth=2.5, label='Synergistic (spider venom)')
ax3.plot(f1, CI_add, color='#fbbf24', linewidth=2, linestyle='--', label='Additive')
ax3.plot(f1, CI_ant, color='#f87171', linewidth=2, linestyle=':', label='Antagonistic')
ax3.axhline(y=1.0, color='#666', linestyle='-', linewidth=1, alpha=0.5)
ax3.fill_between(f1, 0, 1.0, alpha=0.08, color='#c084fc')
ax3.fill_between(f1, 1.0, 2.5, alpha=0.08, color='#f87171')
ax3.annotate('CI < 1: Synergism', xy=(0.5, 0.3), fontsize=10, color='#c084fc', ha='center', fontweight='bold')
ax3.annotate('CI > 1: Antagonism', xy=(0.5, 1.8), fontsize=10, color='#f87171', ha='center')
ax3.set_xlabel('Fraction of Component 1 (Neurotoxin)', color='white', fontsize=11)
ax3.set_ylabel('Combination Index (CI)', color='white', fontsize=11)
ax3.set_title('Venom Synergism\nCombination Index Analysis', color='#d8b4fe', fontsize=11)
ax3.set_ylim(0, 2.5); ax3.set_xlim(0, 1)
ax3.tick_params(colors='white'); ax3.grid(True, alpha=0.2, color='#a78bfa')
for sp in ax3.spines.values(): sp.set_color('#7c3aed')
ax3.legend(fontsize=9, facecolor='#1a0a2e', edgecolor='#7c3aed', labelcolor='white')

plt.tight_layout()
plt.savefig('output.png', dpi=150, bbox_inches='tight', facecolor='#030712')

print("ICK Fold Stability:")
print("  ICK motif unfolding barrier: >40 kJ/mol (extraordinary protease resistance)")
print("  Linear 3-disulfide: ~25 kJ/mol")
print("  Single disulfide: ~12 kJ/mol")
print("\nVenom Gene Evolution (Ka/Ks):")
toxin_median = np.median(ka_ks_toxin)
house_median = np.median(ka_ks_house)
print(f"  Toxin genes median Ka/Ks: {toxin_median:.2f} (positive selection)")
print(f"  Housekeeping genes median Ka/Ks: {house_median:.2f} (purifying selection)")
print(f"  Toxin genes evolve {toxin_median/house_median:.1f}x faster")
print("\nSynergism Analysis:")
print(f"  Minimum CI (synergistic): {CI_syn.min():.3f}")
print("  Spider venom components act synergistically: neurotoxin + cytotoxin + spreading factor")

Click Run to execute the Python code

Code will be executed with Python 3 on the server

References

• King, G.F. & Hardy, M.C. (2013). Spider-venom peptides: structure, pharmacology, and potential for control of insect pests. Annual Review of Entomology, 58, 475–496.
• Nicholson, G.M. (2007). Insect-selective spider toxins targeting voltage-gated sodium channels. Toxicon, 49(4), 490–512.
• Binford, G.J. et al. (2009). Sphingomyelinase D from venoms of Loxosceles spiders: evolutionary insights from cDNA sequences and gene structure. Toxicon, 53(5), 547–557.
• Hodgkin, A.L. & Huxley, A.F. (1952). A quantitative description of membrane current and its application to conduction and excitation in nerve. Journal of Physiology, 117(4), 500–544.
• Escoubas, P. et al. (2000). Structure and pharmacology of spider venom neurotoxins. Biochimie, 82(9–10), 893–907.
• Pineda, S.S. et al. (2020). ArachnoServer 3.0: an online resource for automated discovery, analysis and annotation of spider toxins. Bioinformatics, 34(6), 1074–1076.
• Lewis, R.J. & Garcia, M.L. (2003). Therapeutic potential of venom peptides. Nature Reviews Drug Discovery, 2(10), 790–802.
• Foelix, R.F. (2011). Biology of Spiders, 3rd edition. Oxford University Press.

Venom Biochemistry

Featured Lecture

3.1 Venom Composition

Neurotoxins

Cytotoxins

Enzymes

Dose-Response: The Weibull LD\(_{50}\) Model

3.2 Ion Channel Toxins

\(\omega\)-Agatoxin IVA (Ca\(^{2+}\) Channel Blocker)

PhTx3 (Na\(^+\) Channel from Phoneutria)

\(\delta\)-Atracotoxin (Na\(^+\) Channel from Atrax)

Modified Hodgkin-Huxley Model

3.3 Sphingomyelinase D: The Brown Recluse Mechanism

Michaelis-Menten Kinetics for SMase D

3.4 Venom-Derived Pharmaceuticals

\(\omega\)-Conotoxin MVIIA (Prialt™)

Huwentoxin-IV (Analgesic Candidate)

Spider Peptides as Bio-Insecticides

Selectivity Index

3.5 Venom Delivery: Fang Mechanics

Araneomorph (Cross-Fang)

Mygalomorph (Parallel Fang)

Penetration Force Model

3.6 Computational Analysis

Dose-Response Curves Across Spider Species

Weibull Dose-Response: LD50 Comparison

Hodgkin-Huxley with Toxin Modification

Action Potential: Normal vs Toxin-Modified Na+ Channel

3.7 Advanced Venom Biochemistry

3.7.1 Inhibitor Cystine Knot (ICK) Motif

3.7.2 Venom Proteomics & Transcriptomics

3.7.3 Synergistic Venom Action

Neurotoxin

Cytotoxin

Hyaluronidase

3.7.4 Antivenoms & Recombinant Toxins

3.8 Computational Analysis: Advanced Venom Biochemistry

ICK Stability, Venom Proteomics & Synergism Analysis

References