Module 5: Aquatic Biochemistry & Marine Ecosystems

Ecological Biochemistry & Biodiversity

1. Ocean Acidification

The ocean has absorbed approximately 30% of anthropogenic CO\(_2\) emissions since the Industrial Revolution. This carbon uptake comes at a steep cost: the formation of carbonic acid, driving ocean pH downward in a process known as ocean acidification.

The Carbonate Buffering System

When CO\(_2\) dissolves in seawater, it undergoes a cascade of equilibrium reactions:

\[ \text{CO}_2(\text{g}) \rightleftharpoons \text{CO}_2(\text{aq}) + \text{H}_2\text{O} \rightleftharpoons \text{H}_2\text{CO}_3 \rightleftharpoons \text{H}^+ + \text{HCO}_3^- \rightleftharpoons 2\text{H}^+ + \text{CO}_3^{2-} \]

The equilibrium constants (at 25\(Β°\)C, S=35 psu):

\[ K_1 = \frac{[\text{H}^+][\text{HCO}_3^-]}{[\text{CO}_2]} \approx 1.4 \times 10^{-6} \quad (pK_1 \approx 5.85) \]

\[ K_2 = \frac{[\text{H}^+][\text{CO}_3^{2-}]}{[\text{HCO}_3^-]} \approx 1.2 \times 10^{-9} \quad (pK_2 \approx 8.92) \]

The Henderson-Hasselbalch equation for seawater pH:

\[ \text{pH} = pK_1 + \log\!\left(\frac{[\text{HCO}_3^-]}{[\text{CO}_2]}\right) \]

At current seawater pH ~8.1: \([\text{HCO}_3^-]\) dominates (~90%),\([\text{CO}_3^{2-}]\) ~9%, dissolved \([\text{CO}_2]\) ~1%

pH Change Since Pre-Industrial

Since pre-industrial times (CO\(_2\) ~280 ppm), surface ocean pH has dropped from ~8.2 to ~8.1. This 0.1 unit decrease may seem small, but because pH is logarithmic:

\[ \frac{[\text{H}^+]_{\text{now}}}{[\text{H}^+]_{\text{pre-industrial}}} = 10^{-(8.1 - 8.2)} = 10^{0.1} \approx 1.26 \]

This represents a ~26% increase in hydrogen ion concentration(often rounded to 30% when including spatial variability). Under RCP 8.5, pH could drop to ~7.7 by 2100, representing a 150% increase in H\(^+\) β€” unprecedented in the last 20 million years.

Effect on Calcification: Aragonite Saturation

The critical biological impact is on organisms that build calcium carbonate shells and skeletons (corals, pteropods, coccolithophores, foraminifera):

\[ \Omega_{\text{aragonite}} = \frac{[\text{Ca}^{2+}][\text{CO}_3^{2-}]}{K_{\text{sp}}} \]

  • When \(\Omega > 1\): seawater is supersaturated β€” CaCO\(_3\) precipitation is thermodynamically favorable
  • When \(\Omega < 1\): seawater is undersaturated β€” existing CaCO\(_3\) dissolves
  • Coral reef growth typically requires \(\Omega > 3\); at \(\Omega < 2\), reef erosion exceeds accretion

As CO\(_2\) increases, the reaction \(\text{CO}_2 + \text{H}_2\text{O} + \text{CO}_3^{2-} \rightarrow 2\text{HCO}_3^-\)consumes carbonate ions, directly reducing \(\Omega\). Since pre-industrial times, surface \(\Omega_{\text{aragonite}}\) has decreased from ~4.5 to ~3.5 globally.

Deriving the relationship between CO\(_2\) and \(\Omega\):

From the equilibrium expressions, we can show that the carbonate ion concentration scales inversely with dissolved CO\(_2\):

\[ [\text{CO}_3^{2-}] = \frac{K_1 K_2 \cdot [\text{CO}_2]}{[\text{H}^+]^2} = \frac{K_2 \cdot \text{DIC}}{1 + \frac{[\text{H}^+]}{K_2} + \frac{[\text{H}^+]^2}{K_1 K_2}} \]

Since increasing CO\(_2\) increases \([\text{H}^+]\) (lowers pH), the denominator grows and \([\text{CO}_3^{2-}]\) decreases, thus \(\Omega\) decreases. A doubling of atmospheric CO\(_2\) reduces surface \(\Omega\) by approximately 30%.

2. Coral Bleaching

Coral bleaching is the breakdown of the mutualistic symbiosis between reef-building corals and their endosymbiotic algae (Symbiodiniaceae, formerly zooxanthellae). These dinoflagellate algae provide up to 90% of the coral's energy through photosynthesis, and their pigments give corals their characteristic colors.

The Bleaching Mechanism

Under thermal stress (typically 1-2\(Β°\)C above the local summer maximum for extended periods), a cascade of biochemical events leads to symbiont expulsion:

  1. Photosystem II damage: Excess temperature damages the D1 protein in PSII reaction centers, disrupting electron flow
  2. ROS generation: Blocked electron transport leads to triplet chlorophyll formation, which reacts with O\(_2\) to produce reactive oxygen species (superoxide \(\text{O}_2^{\bullet-}\), singlet oxygen \(^1\text{O}_2\), hydrogen peroxide \(\text{H}_2\text{O}_2\)):\[ {}^3\text{Chl}^* + \text{O}_2 \rightarrow {}^1\text{O}_2 \quad (\text{singlet oxygen}) \]
  3. Antioxidant overwhelm: When ROS production exceeds the detoxification capacity of superoxide dismutase (SOD), catalase, and ascorbate peroxidase, oxidative damage spreads to host tissues
  4. Symbiont expulsion: The coral host initiates exocytosis, apoptosis, or in situ degradation of symbiont cells β€” resulting in the white β€œbleached” appearance

Thermal Threshold Model: Degree Heating Weeks

NOAA's Coral Reef Watch uses Degree Heating Weeks (DHW) to predict bleaching events. DHW accumulates thermal stress above the local bleaching threshold:

\[ \text{DHW} = \frac{1}{7} \sum_{\text{days}} \max(T_{\text{SST}} - T_{\text{threshold}}, 0) \]

where \(T_{\text{threshold}}\) is typically the maximum of the monthly mean (MMM) SST + 1\(Β°\)C, summed over the preceding 12 weeks.

The logistic bleaching probability model:

\[ P(\text{bleaching}) = \frac{1}{1 + \exp\!\left(-k(\text{DHW} - \text{DHW}_{\text{crit}})\right)} \]

\(\text{DHW}_{\text{crit}}\) = species-specific threshold (typically 4-8),\(k\) = sensitivity parameter

NOAA alert thresholds:

  • Bleaching Watch: DHW \(\geq\) 0, expected to exceed threshold
  • Alert Level 1 (DHW \(\geq\) 4): Significant bleaching likely
  • Alert Level 2 (DHW \(\geq\) 8): Severe bleaching and mortality likely

The 2014-2017 global bleaching event (the longest on record) affected 75% of the world's tropical reefs, with DHW values exceeding 16 in many regions.

3. Bioluminescence

Bioluminescence β€” the production of light by living organisms β€” has evolved independently at least 50 times across the tree of life. In the deep ocean (below 200 m), an estimated 76% of organisms are bioluminescent, making it the most common form of communication in the largest habitat on Earth.

The Core Chemistry

\[ \text{Luciferin} + \text{O}_2 \xrightarrow{\text{Luciferase}} \text{Oxyluciferin}^* \rightarrow \text{Oxyluciferin} + h\nu \]

The excited-state oxyluciferin (*) relaxes to ground state by emitting a photon

β€œLuciferin” and β€œluciferase” are generic terms β€” at least 11 chemically distinct luciferin/luciferase systems have been identified. The most well-characterized systems:

Firefly (Photinus pyralis)

Uses D-luciferin + ATP + O\(_2\) \(\rightarrow\) oxyluciferin + AMP + PPi + CO\(_2\) + h\(\nu\). Quantum yield \(\Phi = 0.88\) β€” the highest of any chemiluminescent reaction known. Emission at 562 nm (yellow-green). The ATP requirement makes it unique among bioluminescent systems.

Marine (Coelenterazine-based)

Coelenterazine is the most widespread marine luciferin, used by cnidarians, ctenophores, crustaceans, and fish. In Aequorea victoria, the photoprotein aequorin produces blue light (470 nm) that is then shifted to green (508 nm) by GFP (Green Fluorescent Protein) β€” the basis for the 2008 Nobel Prize in Chemistry.

Chemiluminescent Quantum Yield

The quantum yield is the product of three efficiencies:

\[ \Phi_{CL} = \Phi_C \cdot \Phi_{ES} \cdot \Phi_F \]

  • \(\Phi_C\) = chemical yield (fraction of substrate that reacts to form the excited product)
  • \(\Phi_{ES}\) = excitation yield (fraction of product molecules formed in the excited state)
  • \(\Phi_F\) = fluorescence yield (fraction of excited molecules that emit a photon vs non-radiative decay)

For firefly: \(\Phi_C \approx 1.0\), \(\Phi_{ES} \approx 1.0\),\(\Phi_F \approx 0.88\), giving overall \(\Phi_{CL} = 0.88\). For comparison, the best artificial chemiluminescent systems achieve \(\Phi \approx 0.05\text{-}0.15\).

Dinoflagellate Bioluminescence: Mechanically Triggered

Marine dinoflagellates (e.g., Noctiluca scintillans, Pyrocystis) produce the spectacular blue glow seen in breaking waves and boat wakes. Their bioluminescence is mechanically triggered:

The shear stress threshold model:

\[ L(t) = L_{\max} \cdot \Theta(\tau - \tau_{\text{crit}}) \cdot \exp\!\left(-\frac{t}{\tau_{\text{decay}}}\right) \]

where \(\Theta\) is the Heaviside step function, \(\tau\) is fluid shear stress,\(\tau_{\text{crit}} \approx 0.1\text{-}1\) dyn/cm\(^2\), and\(\tau_{\text{decay}} \approx 100\) ms. The mechanism involves:

  1. Mechanical deformation activates stretch-activated Ca\(^{2+}\) channels in the cell membrane
  2. Ca\(^{2+}\) influx triggers an action potential that propagates along the tonoplast membrane
  3. H\(^+\) is released from vacuole into cytoplasmic β€œscintillons” (specialized organelles containing luciferin and luciferase)
  4. pH drop activates dinoflagellate luciferase (optimal at pH 6, inactive at pH 8)
  5. Flash duration ~100 ms, emitting \(\sim 10^8\) photons at 474 nm (blue)

The ecological function is debated: the β€œburglar alarm” hypothesis suggests that bioluminescence attracts predators of the dinoflagellate's grazers, providing indirect defense. The β€œstartle response” hypothesis proposes that the flash directly startles and deters copepod grazers.

4. Extremophile Biochemistry

Extremophiles thrive in environments that would be lethal to most organisms. Their biochemical adaptations reveal the fundamental limits of molecular stability and have provided biotechnologically invaluable enzymes.

Thermophiles: Proteins That Thrive at 100\(Β°\)C

The most famous thermophilic enzyme is Taq polymerase from Thermus aquaticus (discovered in Yellowstone hot springs, \(T_{\text{opt}} = 72Β°\)C), which revolutionized biology by enabling PCR (polymerase chain reaction). Thermophilic protein stability comes from:

  • Increased ionic interactions (salt bridges): ~30% more than mesophilic homologs
  • Tighter hydrophobic packing: reduced internal cavities
  • Shorter surface loops: less conformational flexibility
  • Higher proline content: restricts backbone flexibility
  • Disulfide bonds in intracellular proteins (rare in mesophiles)

Halophiles: Compatible Solutes

Organisms in hypersaline environments (e.g., Dead Sea, salt evaporation ponds) face osmotic stress β€” water exits the cell down its concentration gradient. Two strategies evolved:

β€œSalt-in” Strategy

Haloarchaea (e.g., Halobacterium) accumulate KCl intracellularly to match external osmolality. All proteins must be adapted to function in ~4 M KCl β€” characterized by excess negative surface charges (acidic proteome, pI ~4.5).

β€œCompatible Solute” Strategy

Most halophilic bacteria synthesize organic osmolytes β€” betaine (glycine betaine), ectoine, trehalose, or proline β€” that balance osmotic pressure without disrupting protein function. These β€œcompatible solutes” are preferentially excluded from protein surfaces, thermodynamically stabilizing the native state.

Membrane Fluidity vs Temperature: Lipid Phase Transition

All organisms must maintain membrane fluidity within a narrow range for proper function of membrane proteins, transport, and signaling. The lipid bilayer undergoes a phase transition from gel (ordered) to liquid-crystalline (disordered) state:

\[ \eta(T) = \eta_0 \cdot \exp\!\left(\frac{E_\eta}{RT}\right) \cdot \frac{1}{1 + \exp\!\left(\frac{\Delta H_m}{R}\left(\frac{1}{T} - \frac{1}{T_m}\right)\right)} \]

Membrane viscosity \(\eta\) as a function of temperature, combining Arrhenius behavior with a phase transition at \(T_m\)

Homeoviscous adaptation strategies:

AdaptationEffect on \(T_m\)Used by
Increase unsaturation (C=C bonds)Lower \(T_m\)Psychrophiles, plants in cold
Shorter acyl chainsLower \(T_m\)Cold-adapted bacteria
Branched chains (iso/anteiso)Lower \(T_m\)Psychrophilic Bacillus
Ether-linked lipidsHigher \(T_m\)Thermophilic archaea
Tetraether monolayer membranesMuch higher \(T_m\)Hyperthermophilic archaea

The desaturase response β€” a molecular thermostat:

When temperature drops, membrane-bound desaturases introduce double bonds into existing fatty acids:

\[ \text{R-CH}_2\text{-CH}_2\text{-R'} + \text{O}_2 + \text{NAD(P)H} \xrightarrow{\Delta\text{-desaturase}} \text{R-CH=CH-R'} + \text{H}_2\text{O} + \text{NAD(P)}^+ \]

The cis double bond introduces a ~30\(Β°\) kink, disrupting packing and lowering \(T_m\) by ~10-15\(Β°\)C per bond

5. Ocean Carbonate Chemistry Diagram

The following diagram illustrates the ocean carbonate buffering system, from atmospheric CO\(_2\) dissolution to aragonite saturation zones at depth:

ATMOSPHERECOβ‚‚(g)420 ppm (↑ from 280 ppm)Sea SurfaceCOβ‚‚ dissolvesCOβ‚‚(aq)~1% of DIC+Hβ‚‚OHβ‚‚CO₃(very unstable)K₁HCO₃⁻~90% of DICKβ‚‚CO₃²⁻~9% of DIC+ H⁺+ H⁺pH ScaleAcidicNeutralBasic8.2Pre-ind.8.1Now7.72100?Aragonite Saturation (Ξ©)Surface: Ξ© β‰ˆ 3.5Supersaturated β†’ shells growCorals, pteropods thriveMid-depth: Ξ© β‰ˆ 1-2Marginal β†’ slow calcificationDeep: Ξ© < 1Undersaturated β†’ shells dissolveSaturation horizon rises with COβ‚‚Increasing depthKey RelationshipspH = pK₁ + log([HCO₃⁻]/[COβ‚‚])Ξ© = [Ca²⁺][CO₃²⁻] / Ksp↑ COβ‚‚ β†’ ↑ H⁺ β†’ ↓ CO₃²⁻ β†’ ↓ Ξ©(The ocean acidification chain reaction)Vulnerable OrganismsCoralsaragonitePteropodsaragoniteCoccolithophorescalciteForaminiferacalciteSaturationHorizonRiseswith COβ‚‚

6. Computational Simulations

The following simulations model: (1) ocean pH projections under three RCP scenarios; (2) coral bleaching probability as a function of degree heating weeks for species with different thermal sensitivities; (3) bioluminescence quantum yield comparison across taxa; and (4) aragonite saturation depth under different CO\(_2\) scenarios.

Python
script.py211 lines

Click Run to execute the Python code

Code will be executed with Python 3 on the server

References

  1. Doney, S.C. et al. (2009). Ocean acidification: the other CO2 problem. Annual Review of Marine Science, 1, 169-192.
  2. Hoegh-Guldberg, O. et al. (2007). Coral reefs under rapid climate change and ocean acidification. Science, 318(5857), 1737-1742.
  3. Hughes, T.P. et al. (2018). Spatial and temporal patterns of mass bleaching of corals in the Anthropocene. Science, 359(6371), 80-83.
  4. Haddock, S.H.D. et al. (2010). Bioluminescence in the sea. Annual Review of Marine Science, 2, 443-493.
  5. Shimomura, O. (2006). Bioluminescence: Chemical Principles and Methods. World Scientific Publishing.
  6. Rothschild, L.J. & Mancinelli, R.L. (2001). Life in extreme environments. Nature, 409(6823), 1092-1101.
  7. Zeebe, R.E. & Wolf-Gladrow, D. (2001). CO2 in Seawater: Equilibrium, Kinetics, Isotopes. Elsevier.
  8. Feely, R.A. et al. (2004). Impact of anthropogenic CO2 on the CaCO3 system in the oceans. Science, 305(5682), 362-366.
←M4: Climate-Dependent BiochemistryM6: Biodiversity-Stabilityβ†’