Part 2: Water, pH & Biological Buffers

Water Properties

Water is arguably the most important molecule in biology. Although deceptively simple — just two hydrogen atoms covalently bonded to one oxygen atom — water's unique physical and chemical properties underpin virtually every biochemical reaction. Life as we know it evolved in an aqueous environment, and no organism can survive without water. In this section, we explore why water behaves the way it does and how its properties create the ideal medium for life.

The Water Molecule: Structure & Polarity

Oxygen is more electronegative (3.44 on the Pauling scale) than hydrogen (2.20), so the shared electrons in each O–H bond are pulled closer to the oxygen nucleus. This unequal electron distribution creates a permanent dipole moment of 1.85 Debye. The molecule adopts a bent geometry with a bond angle of approximately 104.5° (compared to the ideal tetrahedral angle of 109.5°), because oxygen's two lone electron pairs repel each other and the bonding pairs.

Think of water like a tiny magnet: one end (oxygen) carries a partial negative charge (δ−) and the other end (hydrogens) carries a partial positive charge (δ+). This polarity is the master key that unlocks nearly all of water's extraordinary properties.

Hydrogen Bonding

The single most important consequence of water's polarity is its ability to form hydrogen bonds. A hydrogen bond occurs when a hydrogen atom covalently bonded to an electronegative atom (like O or N) is attracted to a lone pair on another electronegative atom.

In liquid water at any instant, each molecule participates in an average of 3.4 hydrogen bonds with neighboring water molecules, creating a dynamic, flickering network that constantly breaks and reforms on a picosecond timescale (~10⁻¹² seconds). Although each individual hydrogen bond is weak (~20 kJ/mol, compared to ~460 kJ/mol for an O–H covalent bond), the collective effect of billions of hydrogen bonds confers extraordinary properties.

Properties Arising from Hydrogen Bonding

Analogy — The Velcro Effect: Imagine each hydrogen bond as a single hook-and-loop fastener (Velcro). One hook is easy to pull apart. But if you cover an entire surface with thousands of tiny hooks (like water's extensive H-bond network), the collective grip becomes formidable. That's why water has such unusually high boiling points and heat capacities compared to molecules of similar molecular weight.

Solvent Properties

Water is often called the "universal solvent" — a slight exaggeration, but one rooted in truth. Water dissolves more substances than any other common liquid. Its ability to dissolve compounds depends on the type of solute:

Ionic Compounds (Electrolytes)

When NaCl dissolves, water molecules surround each Na⁺ and Cl⁻ ion in an organized arrangement called a hydration shell (or solvation shell). The δ− oxygen atoms point toward cations, while the δ+ hydrogen atoms point toward anions. The energy released from forming these ion-dipole interactions (the hydration enthalpy) must be comparable to or greater than the lattice energy to dissolve the crystal.

Polar Molecules

Sugars, amino acids, and small alcohols dissolve via hydrogen bonds and dipole-dipole interactions with water. Glucose, for instance, has five hydroxyl (–OH) groups that readily hydrogen-bond with water, making it highly soluble (~900 g/L at 25°C).

Nonpolar Molecules — The Hydrophobic Effect

Nonpolar molecules like fats and oils cannot participate in hydrogen bonding. When forced into water, they disrupt the H-bond network. Water molecules organize into ordered "cages" (clathrate-like structures) around the nonpolar solute, which is entropically unfavorable. The system minimizes this penalty by driving nonpolar molecules together — the hydrophobic effect.

Amphipathic Molecules

Molecules with both polar and nonpolar regions (amphipathic) orient themselves at water interfaces. Phospholipids, for example, arrange their polar head groups toward water and their fatty acid tails away from water, spontaneously forming bilayer membranes. Detergents and bile salts exploit the same principle, forming micelles that emulsify fats. This self-assembly behavior is entirely driven by water's properties — no biological machinery required.

import numpy as np

# Demonstrate hydration energy calculation
# Hydration enthalpy vs lattice energy determines solubility

compounds = {
    'NaCl': {'lattice_energy': 786, 'hydration_energy': 783, 'soluble': True},
    'KBr': {'lattice_energy': 672, 'hydration_energy': 679, 'soluble': True},
    'AgCl': {'lattice_energy': 905, 'hydration_energy': 851, 'soluble': False},
    'BaSO4': {'lattice_energy': 2423, 'hydration_energy': 2295, 'soluble': False},
    'CaCl2': {'lattice_energy': 2195, 'hydration_energy': 2247, 'soluble': True},
}

print("Compound | Lattice (kJ/mol) | Hydration (kJ/mol) | ΔE (kJ/mol) | Soluble?")
print("-" * 78)
for name, data in compounds.items():
    delta = data['hydration_energy'] - data['lattice_energy']
    prediction = "Yes" if delta >= 0 else "No"
    actual = "Yes" if data['soluble'] else "No"
    match = "✓" if prediction == actual else "≈"
    print(f"{name:8s} | {data['lattice_energy']:17d} | {data['hydration_energy']:18d} | {delta:+11d} | {actual:3s} {match}")

print("\nNote: Positive ΔE (hydration > lattice) favors dissolution")
print("Entropy contributions also play a role in real solubility")

The pH Scale

The concentration of hydrogen ions (H⁺ — or more precisely, hydronium ions H₃O⁺) in a solution has a profound effect on virtually every biochemical process. Because these concentrations span many orders of magnitude, the Danish chemist Søren Sørensen introduced the pH scale in 1909 while working at the Carlsberg Laboratory in Copenhagen (yes, the beer company funded foundational biochemistry research!).

The Autoionization of Water

Pure water undergoes a very slight autoionization (self-dissociation):

H₂O ⇌ H⁺ + OH⁻

At 25°C, only about 1 in every 555 million water molecules is ionized at any given moment. The ion product of water (K_w) equals 1.0 × 10⁻¹⁴, giving pure water a [H⁺] of 1.0 × 10⁻⁷ M, corresponding to pH 7.00 — neutral.

Analogy: Imagine a stadium with 555 million spectators (water molecules). At any instant, only ONE person stands up (ionizes). Yet that one person's behavior affects the entire stadium — just as tiny changes in [H⁺] dramatically alter enzyme activity and protein structure.

The Logarithmic Nature of pH

Because pH is logarithmic, each one-unit change in pH represents a 10-fold change in [H⁺]. A solution at pH 3 has 10× more H⁺ than one at pH 4, and 10,000× more H⁺ than one at pH 7. This has a critical clinical implication:

Biological pH Ranges

Different cellular compartments maintain distinct pH values tailored to their function. This pH compartmentalization is a hallmark of eukaryotic cell biology:

Why pH Matters for Enzymes

Enzymes have amino acid residues in their active sites whose ionization states depend on pH. A histidine residue (pK_a ≈ 6.0) may need to be protonated to catalyze a reaction; at pH 7.4 most histidines are deprotonated. This explains why each enzyme has a characteristic pH optimum — the pH at which it achieves maximum catalytic rate:

import numpy as np

# Calculate pH from [H+] and vice versa
# Demonstrate the logarithmic relationship

h_concentrations = [1e-1, 1e-3, 4e-8, 1e-7, 1e-10, 1e-14]
labels = ['Stomach acid', 'Vinegar', 'Blood (7.4)',
          'Pure water', 'Baking soda', 'Drain cleaner']

print("Solution           | [H⁺] (M)      | pH    | [OH⁻] (M)     | pOH")
print("-" * 80)
for conc, label in zip(h_concentrations, labels):
    pH = -np.log10(conc)
    pOH = 14 - pH
    oh_conc = 10**(-pOH)
    print(f"{label:20s}| {conc:14.2e} | {pH:5.2f} | {oh_conc:14.2e} | {pOH:5.2f}")

print("\n--- pH Change Impact ---")
# Show that small pH changes mean large [H+] changes
ph_normal = 7.40
ph_acidosis = 7.20
ph_severe = 7.00

h_normal = 10**(-ph_normal)
h_acidosis = 10**(-ph_acidosis)
h_severe = 10**(-ph_severe)

print(f"Normal blood pH {ph_normal}: [H⁺] = {h_normal*1e9:.1f} nM")
print(f"Mild acidosis pH {ph_acidosis}: [H⁺] = {h_acidosis*1e9:.1f} nM ({h_acidosis/h_normal:.1f}x increase)")
print(f"Severe acidosis pH {ph_severe}: [H⁺] = {h_severe*1e9:.1f} nM ({h_severe/h_normal:.1f}x increase)")

Henderson-Hasselbalch Equation

The Henderson-Hasselbalch equation is perhaps the most practically important equation in acid-base biochemistry. Derived independently by Lawrence Joseph Henderson (1908) and Karl Albert Hasselbalch (1917), it relates the pH of a solution to the pK_a of an acid and the ratio of its conjugate base to acid forms.

Derivation from the K_a Expression

For a weak acid HA dissociating in water:

HA ⇌ H⁺ + A⁻

The acid dissociation constant is:

K_a = [H⁺][A⁻] / [HA]

Taking the negative logarithm of both sides:

−log K_a = −log [H⁺] − log([A⁻] / [HA])

Substituting pK_a = −log K_a and pH = −log[H⁺]:

Key Predictions of the Equation

Applications

Application 1: Laboratory Buffer Preparation

When preparing a buffer solution for a biochemistry experiment, you need to choose an acid/base pair whose pK_a is close to your desired pH. For example, to make a pH 7.4 buffer, you might choose:

• HEPES (pK_a = 7.55) — widely used in cell culture
• Tris (pK_a = 8.06) — common in molecular biology, but note its pK_a is temperature-dependent
• Phosphate (pK_a2 = 6.86) — inexpensive, but precipitates Ca²⁺ and Mg²⁺

Application 2: Predicting Drug Absorption (Ion Trapping)

Application 3: Amino Acid Charge Prediction

Amino acids have multiple ionizable groups. Using H-H, you can predict the charge state at any pH. For glycine (pK_a1 = 2.34 for −COOH; pK_a2 = 9.60 for −NH₃⁺):

• pH 1: Both groups protonated → net charge +1 (cation)
• pH 5.97 (pI): −COOH deprotonated, −NH₃⁺ protonated → net charge 0 (zwitterion)
• pH 11: Both groups deprotonated → net charge −1 (anion)

import numpy as np

# Henderson-Hasselbalch Calculator
# Calculate pH from pKa and ratio, or ratio from pH and pKa

def henderson_hasselbalch_pH(pKa, ratio_base_acid):
    """Calculate pH given pKa and [A-]/[HA] ratio"""
    return pKa + np.log10(ratio_base_acid)

def henderson_hasselbalch_ratio(pH, pKa):
    """Calculate [A-]/[HA] ratio given pH and pKa"""
    return 10**(pH - pKa)

# Example 1: Bicarbonate buffer in blood
pKa_carbonic = 6.1
ratio_blood = 20  # [HCO3-]/[H2CO3] in normal blood
pH_blood = henderson_hasselbalch_pH(pKa_carbonic, ratio_blood)
print(f"Bicarbonate buffer: pKa={pKa_carbonic}, [HCO3-]/[H2CO3]={ratio_blood}")
print(f"  → pH = {pH_blood:.2f}")

# Example 2: Aspirin absorption
pKa_aspirin = 3.5
pH_stomach = 1.5
pH_blood_plasma = 7.4
ratio_stomach = henderson_hasselbalch_ratio(pH_stomach, pKa_aspirin)
ratio_plasma = henderson_hasselbalch_ratio(pH_blood_plasma, pKa_aspirin)
pct_HA_stomach = 100 / (1 + ratio_stomach)
pct_HA_plasma = 100 / (1 + ratio_plasma)
print(f"\nAspirin (pKa={pKa_aspirin}):")
print(f"  Stomach (pH {pH_stomach}): {pct_HA_stomach:.1f}% uncharged (absorbable)")
print(f"  Plasma  (pH {pH_blood_plasma}): {pct_HA_plasma:.4f}% uncharged (trapped as ion)")

# Example 3: Buffer capacity across pH range
pKa_phosphate = 6.86
print(f"\nPhosphate buffer (pKa={pKa_phosphate}) — percent base form:")
for pH in np.arange(5.0, 9.0, 0.5):
    ratio = henderson_hasselbalch_ratio(pH, pKa_phosphate)
    pct_base = 100 * ratio / (1 + ratio)
    in_range = "✓ buffer zone" if abs(pH - pKa_phosphate) <= 1 else "  outside range"
    print(f"  pH {pH:.1f}: {pct_base:6.1f}% HPO4²⁻   {in_range}")

Biological Buffer Systems

Living organisms produce enormous quantities of acid continuously. Cellular metabolism generates roughly 15,000 mmol of CO₂ per day (volatile acid) and approximately 50–100 mEq of non-volatile acids (sulfuric, phosphoric, and organic acids) daily. Without robust buffering systems, blood pH would plummet within minutes. The body employs three major buffer systems, each operating in complementary compartments and timescales.

Bicarbonate Buffer

The bicarbonate buffer system is the most important extracellular buffer in the body, accounting for about 75% of blood's buffering capacity. Its apparent paradox — a pK_a of 6.1 buffering at pH 7.4, well outside the ±1 rule — is resolved by the fact that CO₂ is an open system.

The Equilibrium

CO₂ + H₂O ⇌ H₂CO₃ ⇌ H⁺ + HCO₃⁻

The enzyme carbonic anhydrase (one of the fastest enzymes known, with a k_cat of ~10⁶ s⁻¹) rapidly catalyzes the first step. In the Henderson-Hasselbalch form:

pH = 6.1 + log([HCO₃⁻] / [0.03 × P_CO₂])

where 0.03 is the solubility coefficient of CO₂ in mmol/L per mmHg. Normal values: [HCO₃⁻] = 24 mEq/L, P_CO₂ = 40 mmHg, giving [H₂CO₃] = 0.03 × 40 = 1.2 mM:

pH = 6.1 + log(24/1.2) = 6.1 + log(20) = 6.1 + 1.30 = 7.40

Phosphate Buffer

The phosphate buffer system (H₂PO₄⁻/HPO₄²⁻, pK_a2 = 6.86) is the primary intracellular buffer and plays a major role in renal acid excretion. Its pK_a of 6.86 is much closer to intracellular pH (~7.2) than the bicarbonate pK_a, making it a more "traditional" buffer within cells.

Why Phosphate Dominates Inside Cells

• High intracellular concentration: Total intracellular phosphate is ~75 mM (vs ~1 mM in plasma), providing excellent buffering capacity
• Ideal pK_a: At pH 7.2, approximately 61% is HPO₄²⁻ (base) and 39% is H₂PO₄⁻ (acid) — close to the 1:1 ratio for optimal buffering
• Important in urine: As urine pH drops to 4.5–6.0 in the renal tubules, the phosphate buffer is in its prime buffering range, helping the kidneys excrete excess H⁺ as titratable acid

Protein Buffers

Proteins are the most abundant buffers in the body, accounting for approximately 60% of total body buffering capacity. Their buffering power comes from ionizable amino acid side chains, particularly histidine (pK_a ≈ 6.0 of the imidazole group), which operates near physiological pH.

Hemoglobin: The Master Blood Buffer

Hemoglobin (Hb) is the most important protein buffer in blood, not just because it's abundant (150 g/L) but because its buffering capacity changes depending on oxygenation state:

import numpy as np

# Compare buffer systems at different pH values
# Calculate buffering capacity (β) = dB/dpH

def buffer_capacity(C_total, pH, pKa):
    """Calculate buffer capacity at given pH
    β = 2.303 * C * Ka * [H+] / (Ka + [H+])^2
    """
    H = 10**(-pH)
    Ka = 10**(-pKa)
    return 2.303 * C_total * Ka * H / (Ka + H)**2

# Buffer systems with physiological concentrations
buffers = {
    'Bicarbonate': {'C': 24e-3, 'pKa': 6.1},    # 24 mM in plasma
    'Phosphate':   {'C': 1e-3, 'pKa': 6.86},     # ~1 mM in plasma
    'Phosphate (intracellular)': {'C': 75e-3, 'pKa': 6.86},  # ~75 mM
    'Histidine (protein)': {'C': 10e-3, 'pKa': 6.0},  # approximate
}

pH_range = np.arange(5.5, 8.5, 0.1)

print("Buffer Capacity (β, mol/L per pH unit) at Key pH Values")
print("-" * 75)
print(f"{'Buffer System':30s} | {'pH 6.1':>8s} | {'pH 6.9':>8s} | {'pH 7.2':>8s} | {'pH 7.4':>8s}")
print("-" * 75)
for name, params in buffers.items():
    vals = [buffer_capacity(params['C'], ph, params['pKa'])
            for ph in [6.1, 6.9, 7.2, 7.4]]
    print(f"{name:30s} | {vals[0]:8.4f} | {vals[1]:8.4f} | {vals[2]:8.4f} | {vals[3]:8.4f}")

print("\nNote: Bicarbonate's actual capacity in vivo is much higher")
print("because it operates as an OPEN system (CO₂ exhaled by lungs)")

Clinical Acid-Base Disorders

Acid-base disorders occur when the body's pH regulation systems are overwhelmed or malfunction. These disorders are among the most common and critical problems in emergency medicine, intensive care, and nephrology. Understanding them requires integrating everything we've learned about pH, buffers, and Henderson-Hasselbalch.

The Four Primary Disorders

Arterial Blood Gas (ABG) Interpretation

The ABG is the gold standard test for assessing acid-base status. A systematic approach:

1. Check pH: <7.35 = acidemia; >7.45 = alkalemia
2. Check P_CO₂: If moving in the opposite direction of pH → respiratory disorder; if same direction → metabolic disorder
3. Check HCO₃⁻: If moving in the same direction as pH → metabolic component
4. Assess compensation: Is the opposing system responding? (Respiratory responds in minutes; renal takes hours-days)
5. Calculate anion gap: AG = Na⁺ − (Cl⁻ + HCO₃⁻); normal 8–12 mEq/L

The Anion Gap: Diagnostic Power

The anion gap is a calculated value that helps narrow the differential diagnosis of metabolic acidosis. When unmeasured anions (like lactate, ketoacids, or toxins) accumulate, they "replace" bicarbonate, widening the gap:

Compensation Rules (Winter's Formula & Others)

The body never overcompensates — compensation brings the pH closer to normal but doesn't fully correct it. Expected compensation formulas help us detect mixed disorders (when two primary disorders coexist):

import numpy as np

# ABG Interpreter - Systematic acid-base analysis
def interpret_abg(pH, pCO2, HCO3, Na=140, Cl=104):
    """Interpret arterial blood gas results"""
    print(f"=== ABG Analysis ===")
    print(f"pH: {pH:.2f}  |  PCO₂: {pCO2} mmHg  |  HCO₃⁻: {HCO3} mEq/L")

    # Step 1: Acidemia or alkalemia?
    if pH < 7.35:
        status = "ACIDEMIA"
    elif pH > 7.45:
        status = "ALKALEMIA"
    else:
        status = "Normal pH"
    print(f"\n1. pH status: {status}")

    # Step 2: Primary disorder
    if pH < 7.35:
        if pCO2 > 45:
            primary = "Respiratory Acidosis"
        elif HCO3 < 22:
            primary = "Metabolic Acidosis"
        else:
            primary = "Mixed disorder"
    elif pH > 7.45:
        if pCO2 < 35:
            primary = "Respiratory Alkalosis"
        elif HCO3 > 26:
            primary = "Metabolic Alkalosis"
        else:
            primary = "Mixed disorder"
    else:
        primary = "No primary disorder (or fully compensated)"
    print(f"2. Primary disorder: {primary}")

    # Step 3: Anion gap
    AG = Na - (Cl + HCO3)
    print(f"3. Anion Gap: {Na} - ({Cl} + {HCO3}) = {AG} mEq/L", end="")
    if AG > 12:
        print(" (ELEVATED - suggests organic acid accumulation)")
    else:
        print(" (Normal)")

    # Step 4: Compensation check (Winter's formula for metabolic acidosis)
    if "Metabolic Acidosis" in primary:
        expected_pCO2 = 1.5 * HCO3 + 8
        print(f"4. Winter's formula: Expected PCO₂ = {expected_pCO2:.0f} ± 2 mmHg")
        if abs(pCO2 - expected_pCO2) <= 2:
            print("   → Appropriate respiratory compensation")
        elif pCO2 > expected_pCO2 + 2:
            print("   → Concurrent respiratory acidosis!")
        else:
            print("   → Concurrent respiratory alkalosis!")

    return primary

# Case 1: DKA
print("--- Case 1: Diabetic Ketoacidosis ---")
interpret_abg(pH=7.15, pCO2=20, HCO3=7, Na=140, Cl=105)

print("\n--- Case 2: COPD with CO₂ retention ---")
interpret_abg(pH=7.32, pCO2=60, HCO3=30, Na=138, Cl=98)

print("\n--- Case 3: Anxiety hyperventilation ---")
interpret_abg(pH=7.52, pCO2=25, HCO3=20, Na=140, Cl=104)

Exercises & Practice Problems

Interactive Worksheet Tool

Conclusion & Next Steps

Water is far more than a passive solvent — it is an active participant in biochemistry. Its polarity and hydrogen bonding network endow it with remarkable properties: high heat capacity for thermal buffering, high heat of vaporization for evaporative cooling, surface tension for capillary action, and the ability to dissolve a vast range of biological molecules. The hydrophobic effect, driven by water's H-bond network, is the dominant force behind protein folding and membrane assembly.

The pH scale provides a compact way to express hydrogen ion concentrations that span 14 orders of magnitude. The Henderson-Hasselbalch equation connects pH to pK_a and concentration ratios, enabling buffer preparation, drug absorption predictions, and amino acid charge calculations. The body's three major buffer systems — bicarbonate (extracellular), phosphate (intracellular), and protein/hemoglobin — work in concert with respiratory and renal compensatory mechanisms to maintain blood pH within the narrow 7.35–7.45 range essential for life.