Module 0: Chemical Foundations of Biodiversity

Life on Earth runs on chemistry. The staggering diversity of species โ€” from chemolithotrophic archaea in deep-sea vents to tropical canopy epiphytes โ€” is underpinned by a shared core of metabolic reactions (glycolysis, TCA cycle, electron transport) overlaid with kingdom-specific secondary metabolism. This module derives the thermodynamic constraints that shape biochemical diversity, introduces the Metabolic Theory of Ecology (MTE) linking body size and temperature to species richness, and explores how chemical niche partitioning generates and maintains biodiversity from first principles.

0.1 Metabolic Diversity Across Kingdoms

All living organisms can be classified by their carbon and energy sources. This fundamental metabolic classification predicts ecological roles, habitat preferences, and biochemical investment strategies.

Trophic Classification

Photoautotrophs (plants, cyanobacteria, algae) fix COโ‚‚ using light energy via photosynthesis. They power nearly all terrestrial and aquatic food webs. The net reaction of oxygenic photosynthesis is:

\[ 6\text{CO}_2 + 6\text{H}_2\text{O} \xrightarrow{h\nu} \text{C}_6\text{H}_{12}\text{O}_6 + 6\text{O}_2 \qquad \Delta G^\circ = +2870 \text{ kJ/mol} \]

Chemoautotrophs (chemolithoautotrophs) derive energy from inorganic redox reactions. These organisms dominate deep-sea hydrothermal vents and subsurface ecosystems. Key examples include sulfur-oxidizing bacteria:

\[ \text{H}_2\text{S} + 2\text{O}_2 \to \text{H}_2\text{SO}_4 \qquad \Delta G^\circ = -798 \text{ kJ/mol} \]

and iron-oxidizing bacteria (Acidithiobacillus ferrooxidans):

\[ 4\text{Fe}^{2+} + \text{O}_2 + 4\text{H}^+ \to 4\text{Fe}^{3+} + 2\text{H}_2\text{O} \qquad \Delta G^\circ = -30 \text{ kJ/mol per Fe}^{2+} \]

Heterotrophs (animals, fungi, most bacteria) obtain both carbon and energy from organic molecules produced by autotrophs. Their metabolic efficiency depends on the electron acceptor used:

\[ \text{C}_6\text{H}_{12}\text{O}_6 + 6\text{O}_2 \to 6\text{CO}_2 + 6\text{H}_2\text{O} \qquad \Delta G^\circ = -2870 \text{ kJ/mol (aerobic)} \]

Mixotrophs combine autotrophic and heterotrophic strategies. Many dinoflagellates and chrysophytes can photosynthesize while also ingesting bacteria or dissolved organic matter. This metabolic flexibility provides advantages in nutrient-poor environments and explains the dominance of mixotrophy in oligotrophic ocean gyres, where up to 50% of bacterivory is performed by photosynthetic flagellates.

Primary vs Secondary Metabolism

Primary metabolism encompasses universally conserved pathways essential for growth and reproduction: glycolysis (Embden-Meyerhof-Parnas pathway), the tricarboxylic acid (TCA/Krebs) cycle, the electron transport chain (ETC), amino acid biosynthesis, and nucleotide metabolism. These pathways are found in virtually all organisms and show remarkable sequence conservation across billions of years of evolution.

Secondary metabolism produces species-specific or lineage-specific compounds with no direct role in growth but critical ecological functions. Major classes include:

  • Alkaloids (plants): nitrogen-containing compounds derived from amino acids โ€” nicotine, morphine, caffeine, strychnine (~27,000 known structures)
  • Terpenoids (plants, fungi): built from C5 isoprene units via MEP or MVA pathways โ€” essential oils, carotenoids, steroids (~55,000 known)
  • Phenolics (plants): aromatic compounds from shikimate/phenylpropanoid pathway โ€” flavonoids, tannins, lignin (~10,000 known)
  • Mycotoxins (fungi): aflatoxins, ochratoxin A, trichothecenes โ€” competitive weapons and food contaminants
  • Antibiotics (bacteria, fungi): streptomycin, penicillin, vancomycin โ€” chemical warfare in soil microbial communities
  • Venoms (animals): complex mixtures of peptides, enzymes, and small molecules โ€” predation and defense

Thermodynamic Efficiency of Photosynthesis vs Chemosynthesis

The maximum thermodynamic efficiency of energy capture can be derived from the ratio of free energy stored in products to energy input. For oxygenic photosynthesis, solar photons at 680 nm carry energy:

\[ E_\text{photon} = \frac{hc}{\lambda} = \frac{(6.626 \times 10^{-34})(3 \times 10^8)}{680 \times 10^{-9}} = 2.92 \times 10^{-19} \text{ J} = 1.82 \text{ eV} \]

Eight photons are needed per COโ‚‚ fixed (4 at PSII + 4 at PSI), so the minimum photon energy input per glucose is:

\[ E_\text{input} = 6 \times 8 \times 1.82 \text{ eV} = 87.4 \text{ eV} = 8430 \text{ kJ/mol glucose} \]

With \(\Delta G_\text{glucose} = 2870\) kJ/mol stored, the maximum thermodynamic efficiency of photosynthesis is:

\[ \eta_\text{photo} = \frac{2870}{8430} \approx 34\% \]

In practice, real-world efficiencies are 1-2% (C3 plants) to 3-4% (C4 plants) due to photorespiration, light saturation, and reflection losses. For chemosynthesis via hydrogen oxidation (\(\text{H}_2 + \frac{1}{2}\text{O}_2 \to \text{H}_2\text{O}\),\(\Delta G = -237\) kJ/mol), the efficiency of carbon fixation depends on coupling to the Calvin cycle, typically achieving 5-10% thermodynamic efficiency.

0.2 Thermodynamic Constraints on Biochemistry

Every biochemical reaction is governed by thermodynamics. The direction of spontaneous change, the equilibrium position, and the temperature sensitivity of reaction rates all follow from fundamental physical chemistry.

Gibbs Free Energy

The Gibbs free energy change determines whether a reaction proceeds spontaneously at constant temperature and pressure:

\[ \Delta G = \Delta H - T\Delta S \]

where \(\Delta H\) is the enthalpy change (bond energy), \(T\) is absolute temperature, and \(\Delta S\) is the entropy change. A reaction is spontaneous when \(\Delta G < 0\). Under non-standard conditions, the actual free energy depends on reactant and product concentrations:

\[ \Delta G = \Delta G^\circ + RT\ln Q \]

where \(Q\) is the reaction quotient. At equilibrium (\(\Delta G = 0\)), this gives the fundamental relationship:

\[ \Delta G^\circ = -RT\ln K_\text{eq} \]

This equation explains why cells maintain metabolite concentrations far from equilibrium โ€” the further from equilibrium, the larger the driving force (\(-\Delta G\)) for the reaction. The ATP/ADP ratio in healthy cells (\(\sim 10:1\)) ensures that ATP hydrolysis provides \(\Delta G \approx -54\) kJ/mol rather than the standard\(\Delta G^\circ = -30.5\) kJ/mol.

Arrhenius Temperature Dependence

The rate constant of an enzyme-catalyzed reaction increases exponentially with temperature (below denaturation), following the Arrhenius equation:

\[ k(T) = A \cdot \exp\!\left(-\frac{E_a}{RT}\right) \]

where \(A\) is the pre-exponential factor (collision frequency), \(E_a\) is the activation energy, \(R = 8.314\) J/(mol K), and \(T\) is absolute temperature. Taking the ratio of rate constants at two temperatures:

\[ \frac{k_2}{k_1} = \exp\!\left[\frac{E_a}{R}\left(\frac{1}{T_1} - \frac{1}{T_2}\right)\right] \]

Derivation of the Q10 Relationship

The Q10 coefficient quantifies how much a reaction rate increases with a 10 C temperature rise. Starting from the Arrhenius ratio, if \(T_2 - T_1 = 10\) K:

\[ Q_{10} \equiv \frac{k(T + 10)}{k(T)} = \exp\!\left[\frac{E_a}{R} \cdot \frac{10}{T(T+10)}\right] \]

For a general temperature interval, the Q10 can be extracted from any two measurements:

\[ Q_{10} = \left(\frac{k_2}{k_1}\right)^{10/(T_2 - T_1)} \]

For most enzymatic reactions, \(Q_{10} \approx 2\text{-}3\), meaning reaction rates roughly double for each 10 C increase. This has profound ecological consequences:

  • Tropical ecosystems have faster decomposition, nutrient cycling, and primary productivity due to higher ambient temperatures
  • Soil respiration is strongly temperature-dependent (\(Q_{10} \approx 2.4\) globally), creating a positive feedback with climate warming
  • Ectotherms in warm climates have higher metabolic rates, faster growth, shorter generation times, and potentially faster speciation

Water Activity and Enzyme Function

Water activity (\(a_w\)) measures the thermodynamic availability of water for biochemical reactions. Defined as the ratio of vapor pressure of a solution to that of pure water:

\[ a_w = \frac{p}{p_0} \approx \frac{n_\text{water}}{n_\text{water} + n_\text{solute}} \]

Most bacteria require \(a_w > 0.90\) for growth, while xerophilic fungi can grow at \(a_w \approx 0.65\). Halophilic archaea in salt lakes (\(a_w \approx 0.75\)) use compatible solutes (KCl, glycerol, ectoine) to maintain cytoplasmic osmotic balance. The relationship between water activity and the Gibbs free energy of water is:

\[ \mu_w = \mu_w^\circ + RT\ln a_w \]

As \(a_w\) decreases, the chemical potential of water drops, making hydrolysis reactions less favorable and reducing enzyme conformational flexibility. This thermodynamic constraint explains why arid and hypersaline ecosystems have lower microbial diversity despite adequate temperature and nutrient availability.

0.3 The Metabolic Theory of Ecology (MTE)

The Metabolic Theory of Ecology (Brown et al. 2004) unifies body size and temperature as the two master variables governing ecological processes. It builds on the quarter-power allometric scaling of metabolic rate with body mass, combined with the Boltzmann-Arrhenius temperature dependence of biochemical reaction rates.

The Fundamental MTE Equation

Whole-organism metabolic rate \(B\) scales as:

\[ B = B_0 \, M^{3/4} \, \exp\!\left(-\frac{E_a}{kT}\right) \]

where \(B_0\) is a normalization constant, \(M\) is body mass, the exponent \(3/4\) arises from the fractal geometry of resource distribution networks (West, Brown, Enquist 1997), \(E_a \approx 0.6\text{-}0.7\) eV is the average activation energy for rate-limiting metabolic reactions, \(k = 8.617 \times 10^{-5}\) eV/K is Boltzmann's constant, and \(T\) is body temperature in Kelvin.

Deriving Species Richness from MTE

The MTE predicts species richness \(S\) by combining metabolic rate with energy equivalence. Starting from the total energy flux \(J_\text{tot}\) through an ecosystem:

\[ J_\text{tot} = \sum_{i=1}^{S} N_i \, B_i \]

where \(N_i\) is the population size and \(B_i\) is the per-capita metabolic rate of species \(i\). If total energy flux per unit area is approximately constant (determined by resource supply), and if population density scales as:

\[ N \propto M^{-3/4} \cdot \exp\!\left(\frac{E_a}{kT}\right) \]

(the energy equivalence rule: population energy use is independent of body size), then the number of species that can be supported in a community is proportional to the total number of individuals divided by the average population size. The key result is:

\[ \ln S = \ln S_0 - \frac{E_a}{kT_0} + \frac{E_a}{kT} \]

Rearranging:

\[ S \propto \exp\!\left(\frac{E_a}{kT}\right) \]

Wait โ€” this predicts species richness increases as temperature decreases, which contradicts the observed latitudinal gradient. The resolution lies in the speciation rate. Higher temperatures increase mutation rates (via faster metabolism and more reactive oxygen species) and shorten generation times, accelerating speciation. The net prediction, accounting for both energy availability and speciation kinetics, is:

\[ S \propto \exp\!\left(-\frac{E_a}{kT}\right) \cdot t_\text{evol} \]

where \(t_\text{evol}\) reflects evolutionary time (tropics have had longer stable warm periods). Allen et al. (2002) showed that plotting \(\ln S\) vs\(1/kT\) gives a slope of approximately \(-0.65\) eV for diverse taxa, matching the predicted activation energy of metabolism.

Population Density Predictions

The MTE predicts that population density \(N/A\) (individuals per unit area) scales as:

\[ \frac{N}{A} \propto M^{-3/4} \cdot \exp\!\left(\frac{E_a}{kT}\right) \]

This predicts that small organisms at low temperatures have the highest population densities, consistent with the enormous abundances of soil bacteria (\(\sim 10^9\) cells/g) compared to large mammals (\(\sim 10^{-2}\)/kmยฒ).

Ecosystem Productivity

Total ecosystem productivity (gross primary production, GPP) integrates individual metabolic rates over all autotrophs. The MTE predicts:

\[ \text{GPP} = \sum_i N_i B_i \propto \exp\!\left(-\frac{E_a}{kT}\right) \]

This temperature dependence, with \(E_a \approx 0.32\) eV for photosynthesis (lower than the 0.65 eV for heterotrophic metabolism), creates a fundamental asymmetry: warming increases heterotrophic respiration faster than primary production, potentially converting ecosystems from carbon sinks to carbon sources.

0.4 Chemical Niche Partitioning

Secondary metabolites do not merely serve as defenses โ€” they define ecological niches. Each species' unique biochemical profile determines its interactions with competitors, herbivores, pathogens, and mutualists, effectively creating a "chemical niche" in multidimensional metabolite space.

Tilman's Resource Ratio Theory

David Tilman's (1982) resource ratio theory provides a rigorous framework for understanding coexistence through differential resource use. Consider two species competing for two essential resources (\(R_1, R_2\)). Each species has a zero net growth isocline (ZNGI) โ€” the minimum resource combination needed for population maintenance.

The ZNGI for species \(i\) is defined by:

\[ \frac{dN_i}{dt} = N_i \left[ \min\!\left(\frac{R_1}{R_1 + K_{i1}}, \frac{R_2}{R_2 + K_{i2}}\right) \cdot r_{\max,i} - d_i \right] = 0 \]

where \(K_{ij}\) is the half-saturation constant for species \(i\) on resource \(j\), \(r_{\max,i}\) is maximum growth rate, and\(d_i\) is death rate. Coexistence requires that ZNGI lines cross and that each species consumes proportionally more of the resource that more limits its own growth.

The Competitive Exclusion Principle and Its Biochemical Basis

Gause's competitive exclusion principle states that two species competing for a single limiting resource cannot coexist indefinitely โ€” the superior competitor will drive the other to extinction. For two species competing for one resource with Monod kinetics:

\[ \frac{dN_1}{dt} = N_1\left(\frac{\mu_{\max,1} R}{K_1 + R} - d_1\right), \qquad \frac{dN_2}{dt} = N_2\left(\frac{\mu_{\max,2} R}{K_2 + R} - d_2\right) \]

At steady state, each species draws the resource to its \(R^*\) โ€” the minimum resource concentration supporting a stable population:

\[ R_i^* = \frac{K_i \, d_i}{\mu_{\max,i} - d_i} \]

The species with the lower \(R^*\) wins โ€” it can survive at resource concentrations too low for its competitor. The biochemical basis for different \(R^*\) values lies in enzyme kinetics: species with high-affinity transporters (low \(K_m\)) or efficient metabolic pathways (high yield coefficients) achieve lower \(R^*\).

Secondary Metabolites as Niche Dimensions

In Hutchinson's n-dimensional hypervolume niche concept, secondary metabolites add chemical dimensions to the niche space. A plant producing a unique alkaloid cocktail occupies a distinct chemical niche defined by: (1) the herbivores it deters, (2) the pollinators it attracts, (3) the soil microbes it selects, and (4) the allelopathic interactions with neighboring plants. This chemical niche differentiation explains how hundreds of tree species coexist in tropical forests โ€” each occupying a unique position in chemical defense space, attracting different specialist herbivores and mutualists, and partitioning the "enemy-free space" (Janzen-Connell hypothesis).

The metabolic cost of secondary metabolite production creates trade-offs. Carbon invested in terpenes cannot be used for growth. The optimal defense level \(D^*\) maximizes fitness:

\[ W(D) = W_0 \cdot [1 - h(D)] \cdot [1 - c(D)] \]

where \(h(D)\) is the fraction of tissue lost to herbivory (decreasing with defense\(D\)) and \(c(D)\) is the metabolic cost of defense (increasing with\(D\)). The optimum satisfies \(dW/dD = 0\), which yields:

\[ -h'(D^*)[1 - c(D^*)] = c'(D^*)[1 - h(D^*)] \]

At the optimum, the marginal reduction in herbivory equals the marginal increase in metabolic cost, weighted by current fitness.

Metabolic Diversity Tree

This diagram illustrates how primary metabolism (universally conserved) branches into kingdom-specific secondary metabolism, generating the chemical diversity that underlies ecological niche partitioning.

Metabolic Diversity: Primary โ†’ Secondary MetabolismPrimary MetabolismGlycolysis | TCA | ETCAmino acid | Nucleotide synth.Universally conservedPlantsSecondary MetabolismAlkaloids~27,000 knownNicotine (nAChR block)Morphine (opioid agonist)Caffeine (A2A antagonist)Terpenoids~55,000 knownIsoprene (thermotolerance)Menthol (monoterpene)Artemisinin (sesqui.)Phenolics~10,000 knownTannins (protein precip.)Flavonoids (UV screen)Lignin (structural)FungiSecondary MetabolismMycotoxins~400 knownAflatoxin B1 (carcinogen)Ochratoxin A (nephrotox.)Ergotamine (vasospasm)ฮฒ-LactamsAntibioticsPenicillin (cell wall)CephalosporinLovastatin (HMG-CoA)BacteriaSecondary MetabolismAntibiotics~23,000 knownStreptomycin (ribosome)Vancomycin (cell wall)Tetracycline (30S ribo.)SiderophoresIron scavengingEnterobactin (catechol)Pyoverdine (fluoresc.)DesferrioxamineAnimalsSecondary MetabolismVenomsPeptide toxinsConotoxins (ion ch.)ฮฑ-Bungarotoxin (nAChR)Melittin (bee, membr.)PheromonesSignaling moleculesBombykol (moth sex)Alarm pheromonesTrail pheromones (ants)KEY: Primary metabolism is universally conserved across all domains of lifeSecondary metabolites are lineage-specific and create unique chemical nichesCombined ~100,000+ known secondary metabolites drive ecological interactions and biodiversity

Python: MTE Predictions & Q10 Effects

This simulation explores four key predictions of the Metabolic Theory of Ecology: (1) metabolic rate vs body mass at different temperatures, showing the \(M^{3/4}\) scaling and Boltzmann temperature factor; (2) the latitudinal diversity gradient predicted by the temperature dependence of speciation rates; (3) Q10 effects on decomposition rates showing exponential sensitivity; and (4) Tilman's resource ratio theory showing how ZNGI lines determine competitive outcomes.

Python
script.py131 lines

Click Run to execute the Python code

Code will be executed with Python 3 on the server

References

  1. Brown, J.H., Gillooly, J.F., Allen, A.P., Savage, V.M. & West, G.B. (2004). Toward a metabolic theory of ecology. Ecology, 85(7), 1771-1789.
  2. West, G.B., Brown, J.H. & Enquist, B.J. (1997). A general model for the origin of allometric scaling laws in biology. Science, 276(5309), 122-126.
  3. Allen, A.P., Brown, J.H. & Gillooly, J.F. (2002). Global biodiversity, biochemical kinetics, and the energetic-equivalence rule. Science, 297(5586), 1545-1548.
  4. Tilman, D. (1982). Resource Competition and Community Structure. Princeton University Press.
  5. Davidson, E.A. & Janssens, I.A. (2006). Temperature sensitivity of soil carbon decomposition and feedbacks to climate change. Nature, 440(7081), 165-173.
  6. Hartmann, T. (2007). From waste products to ecochemicals: fifty years research of plant secondary metabolism. Phytochemistry, 68(22-24), 2831-2846.
  7. Gause, G.F. (1934). The Struggle for Existence. Williams & Wilkins, Baltimore.
  8. Hutchinson, G.E. (1957). Concluding remarks. Cold Spring Harbor Symposia on Quantitative Biology, 22, 415-427.