Module 2 · BioGeometry

Symmetry & Chirality

From bilateral animals to pentameric echinoderms, from L-amino acids to dextral snail shells - how symmetry is established and broken in living systems.

1. Symmetry Types in Biology

Biological symmetry is categorized by the group of transformations that preserve the organism's form. Three major classes dominate:

Bilateral symmetry (Bilateria)

A single mirror plane (sagittal plane). Symmetry group \(C_2 = \{e, \sigma\}\). All Bilateria share this property: arthropods, molluscs, annelids, vertebrates, flatworms. Bilateral symmetry is advantageous for directional locomotion - a streamlined body with symmetric paired appendages moves efficiently forward.

Radial symmetry

An n-fold rotation axis, symmetry group \(C_n\) or \(D_n\). Cnidarians (jellyfish, sea anemones) show 4- or 6-fold symmetry. Echinoderms (starfish, sea urchins) are famously pentameric - 5-fold symmetry. Flowers vary: Brassicaceae 4-fold, Rosaceae 5-fold, monocots 3- or 6-fold.

Spiral symmetry

No rotation axis fixes the shape, but a screw motion (rotation + translation) does. Snail shells, foraminifera, climbing plants - all use logarithmic spirals of the form \(r = a e^{b\theta}\).

Tetrameric / pentameric flowers

Flower symmetry is genetically controlled by the ABC model (Coen & Meyerowitz 1991): A-genes (apetala) set perianth identity, B-genes (deficiens, globosa) specify petals, C-genes (agamous) specify reproductive organs. Mutations in these genes can convert pentameric flowers to tetrameric or vice versa.

2. Chirality in Molecules

A molecule is chiral if it is not superimposable on its mirror image. The canonical case is a tetrahedral carbon with four different substituents. The two mirror images are enantiomers, labeled D (dextro) and L (levo) or R and S by CIP rules.

Biological homochirality

Life on Earth is homochiral: all 20 proteinogenic amino acids (except glycine) are L-form; all biological sugars are D-form. This uniformity enables complementary base-pairing, proper enzyme folding and stable secondary structures (alpha helix works for L but not D).

Derivation: homochirality from symmetry breaking

Consider a racemic mixture with autocatalytic production of each enantiomer. Frank's 1953 model:

\[ \frac{d[L]}{dt} = k_1 A [L] - k_2 [L][D], \qquad \frac{d[D]}{dt} = k_1 A [D] - k_2 [L][D] \]

The racemic state [L] = [D] is unstable: any tiny asymmetry amplifies exponentially. Asymmetric catalysis by meteoric chiral minerals (e.g. amino acids extracted from the Murchison meteorite show 5-15% L excess), circularly polarized starlight, or parity violation in weak nuclear interactions could provide the initial bias.

Soai reaction: experimental homochirality

Soai (1995) demonstrated in vitro autocatalytic amplification: starting with 0.00005 %e.e. pyrimidyl alcohol produces 99.5 %e.e. product after a few catalytic cycles - a chemical proof of principle for the Frank mechanism.

3. Snail Shells: Dextral vs Sinistral

About 90% of land snails coil to the right (dextral): viewed from the apex, shell opening is on the right. The other 10% are sinistral. In some species (Partula suturalis), both forms coexist; in others (Lymnaea peregra), sinistral individuals are rare mutants.

Genetic basis

Shell chirality is determined in the fertilized egg by maternal effect of a single gene - the formin-1 locus in Lymnaea. It controls the orientation of the spiral cleavage in 4-8 cell embryos. A single nucleotide change in the maternally inherited mRNA flips chirality 50 generations later.

Raup's shell parameters

Raup (1966) parameterized shell shape by three numbers:

\[ W = \frac{r_2}{r_1}, \qquad D = \frac{d}{r}, \qquad T = \text{translation per turn} \]

where W is the whorl expansion rate (geometric growth between turns), D the distance of the aperture from the coiling axis, T the translation per turn. Nautilus: (W, D, T) = (2, 0, 0); turritella-type snail: (1.5, 0, 0.5); ammonite: (1.3, 1, 0.2). The sign of T determines dextral/sinistral.

Most of morphospace is unoccupied by actual species - a classic finding of Raup that echoes Gould's argument that evolution explores only a narrow corridor of the possible.

4. Left-Right Asymmetry in Vertebrates

Despite bilaterally symmetric external appearance, vertebrates have striking internal asymmetries: the heart is on the left, the liver on the right, gut loops clockwise when viewed from above. The asymmetry is established during early embryogenesis by the node - a transient structure whose nodal cilia break symmetry.

The nodal-flow mechanism

Nodal cilia are motile, tilted posteriorly and rotate clockwise when viewed from the ventral side. This tilted rotation generates a unidirectional leftward fluid flow (Nonaka et al., Cell 1998), carrying signaling vesicles (NVPs) containing Sonic Hedgehog and retinoic acid from right to left. Left-side signaling activates the conserved Nodal-Pitx2 cascade, specifying left-sided organ identity.

Derivation: chiral flow from tilted rotation

Consider a cilium attached at origin, tilted posteriorly by angle θ, rotating about its axis with angular velocity ω. Tip traces an ellipse in x-y plane, but the low-Reynolds-number drag is asymmetric between upstrokes (into fluid) and downstrokes (toward wall). Integrating the Stokes-force field over the rotation:

\[ \langle \vec{v}_{\mathrm{fluid}} \rangle_x = 0, \quad \langle \vec{v}_{\mathrm{fluid}} \rangle_y \neq 0 \]

The y-component (leftward, by convention) is non-zero precisely because of the posterior tilt. Cartright, Piro & Tuval (2004) derived an explicit formula and matched experiments within 15%.

Situs inversus

Mutations in ciliary proteins (Dnah11, Lefty, Nodal) produce situs inversus (heart on right, full mirror reversal) or heterotaxy (random arrangement). Kartagener syndrome combines situs inversus, chronic sinusitis, and bronchiectasis - all caused by immotile cilia.

5. SVG: Symmetry Groups in Biology

6. Simulation: Symmetry Classification

Four panels showing bilateral, pentameric, tetrameric and spiral radial profiles, with automatic FFT-based symmetry detection identifying the dominant harmonic (1-fold, 5-fold, 4-fold, etc.)

Python

script.py94 lines

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

# ======================================================================
# Symmetry classification: detecting bilateral, radial, spiral symmetries
# Algorithm: FFT of radial intensity for radial symmetries
# Reflection symmetry test: correlation with mirrored version
# ======================================================================

np.random.seed(7)
theta = np.linspace(0, 2*np.pi, 600)

# 1) Bilateral: two mirror halves
def bilateral(theta, sigma=0.15):
    return 1 + 0.5*np.cos(theta) + 0.2*np.cos(2*theta) + sigma*np.random.randn(len(theta))*0

# 2) Radial n-fold: cos(n theta)
def nfold(theta, n=5):
    return 1 + 0.3*np.cos(n*theta)

# 3) Spiral: radius grows log-spiral
def spiral(theta_range, b=0.2):
    return np.exp(b*theta_range)

r_bi = bilateral(theta)
r_5 = nfold(theta, 5)  # starfish, flowers
r_4 = nfold(theta, 4)  # butterfly, many insects
r_3 = nfold(theta, 3)
th_sp = np.linspace(0, 8*np.pi, 600)
r_sp = spiral(th_sp, 0.2)

fig, axes = plt.subplots(2, 2, figsize=(14, 12))
fig.patch.set_facecolor('#1a1200')

# Panel 1: Bilateral
ax1 = axes[0, 0]
ax1.plot(r_bi*np.cos(theta), r_bi*np.sin(theta), color='#fbbf24', linewidth=2.3)
ax1.fill(r_bi*np.cos(theta), r_bi*np.sin(theta), color='#fbbf24', alpha=0.2)
ax1.axvline(0, color='#f87171', linewidth=1, linestyle='--', label='Mirror plane')
ax1.set_aspect('equal')
ax1.set_title('Bilateral (most animals)', color='white', fontweight='bold')
ax1.legend(facecolor='#1a1200', edgecolor='#fbbf24', labelcolor='white')
ax1.set_facecolor('#1a1200')
ax1.tick_params(colors='white')
ax1.axis('off')

# Panel 2: 5-fold radial
ax2 = axes[0, 1]
ax2.plot(r_5*np.cos(theta), r_5*np.sin(theta), color='#2dd4bf', linewidth=2.3)
ax2.fill(r_5*np.cos(theta), r_5*np.sin(theta), color='#2dd4bf', alpha=0.2)
for k in range(5):
    ax2.plot([0, 1.3*np.cos(k*2*np.pi/5)], [0, 1.3*np.sin(k*2*np.pi/5)], color='#f472b6', linewidth=0.6, linestyle='--')
ax2.set_aspect('equal')
ax2.set_title('Pentameric (echinoderms, flowers)', color='white', fontweight='bold')
ax2.set_facecolor('#1a1200')
ax2.tick_params(colors='white')
ax2.axis('off')

# Panel 3: 4-fold
ax3 = axes[1, 0]
ax3.plot(r_4*np.cos(theta), r_4*np.sin(theta), color='#d946ef', linewidth=2.3)
ax3.fill(r_4*np.cos(theta), r_4*np.sin(theta), color='#d946ef', alpha=0.2)
ax3.set_aspect('equal')
ax3.set_title('Tetrameric (flowers, diatoms)', color='white', fontweight='bold')
ax3.set_facecolor('#1a1200')
ax3.tick_params(colors='white')
ax3.axis('off')

# Panel 4: Spiral
ax4 = axes[1, 1]
ax4.plot(r_sp*np.cos(th_sp), r_sp*np.sin(th_sp), color='#fbbf24', linewidth=2.3)
ax4.set_aspect('equal')
ax4.set_title('Spiral (molluscs, foraminifera)', color='white', fontweight='bold')
ax4.set_facecolor('#1a1200')
ax4.tick_params(colors='white')
ax4.axis('off')

plt.tight_layout()
plt.savefig('output.png', facecolor='#1a1200', dpi=110, bbox_inches='tight')
plt.close()

# Symmetry detection via FFT
def detect_n_fold(r):
    F = np.abs(np.fft.fft(r - np.mean(r)))[:len(r)//2]
    peaks = np.argsort(F)[-5:][::-1]
    return peaks

print('Symmetry detection via FFT of radial profile:')
print(f'  Bilateral (1-fold): dominant harmonic {detect_n_fold(r_bi)[0]}')
print(f'  Pentameric (5-fold): dominant harmonic {detect_n_fold(r_5)[0]}')
print(f'  Tetrameric (4-fold): dominant harmonic {detect_n_fold(r_4)[0]}')

Click Run to execute the Python code

Code will be executed with Python 3 on the server

7. Simulation: Chirality Detection via Signed Volume

L- and D-alanine rendered in 3D. The signed tetrahedral volume of the four substituents distinguishes the enantiomers: positive for L, negative for D - an algorithmic approach to handedness.

Python

script.py71 lines

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

# ======================================================================
# Chirality detection: determinant of 3D configuration
# L-alanine vs D-alanine
# For a tetrahedral carbon with 4 distinct substituents, the oriented
# volume (triple product) sign distinguishes enantiomers.
# ======================================================================

# Alanine: central C with NH2, COOH, H, CH3
# L-form (natural): counter-clockwise priority NH2 > COOH > CH3 > H
def alanine(chirality=1):
    # chirality = +1 for L, -1 for D
    C = np.array([0, 0, 0])
    NH2 = np.array([1, 1, 0])
    COOH = chirality * np.array([-1, 1, 0])
    CH3 = np.array([0, -1, 1])
    H = np.array([0, -1, -1])
    return C, NH2, COOH, CH3, H

L = alanine(+1)
D = alanine(-1)

# Compute signed volume
def signed_volume(v1, v2, v3, v4):
    a = v2 - v1
    b = v3 - v1
    c = v4 - v1
    return np.linalg.det(np.array([a, b, c]))

vol_L = signed_volume(*L[:4])
vol_D = signed_volume(*D[:4])
print(f'Signed volume L-alanine: {vol_L:+.3f}')
print(f'Signed volume D-alanine: {vol_D:+.3f}')

fig = plt.figure(figsize=(15, 6))
fig.patch.set_facecolor('#1a1200')

ax1 = fig.add_subplot(121, projection='3d')
labels = ['C', 'NH2', 'COOH', 'CH3', 'H']
colors = ['#fbbf24', '#60a5fa', '#f87171', '#2dd4bf', '#f472b6']
for pt, lab, c in zip(L, labels, colors):
    ax1.scatter(*pt, s=200, c=c, edgecolors='white', linewidth=1.5)
    ax1.text(pt[0]+0.15, pt[1]+0.15, pt[2]+0.15, lab, color=c, fontsize=11)
for atom in L[1:]:
    ax1.plot([L[0][0], atom[0]], [L[0][1], atom[1]], [L[0][2], atom[2]], color='#fbbf24', linewidth=1)
ax1.set_title(f'L-alanine (signed vol = {vol_L:+.2f})', color='white', fontweight='bold')
ax1.set_facecolor('#1a1200')
ax1.tick_params(colors='white')
ax1.set_xlim(-2, 2); ax1.set_ylim(-2, 2); ax1.set_zlim(-2, 2)

ax2 = fig.add_subplot(122, projection='3d')
for pt, lab, c in zip(D, labels, colors):
    ax2.scatter(*pt, s=200, c=c, edgecolors='white', linewidth=1.5)
    ax2.text(pt[0]+0.15, pt[1]+0.15, pt[2]+0.15, lab, color=c, fontsize=11)
for atom in D[1:]:
    ax2.plot([D[0][0], atom[0]], [D[0][1], atom[1]], [D[0][2], atom[2]], color='#fbbf24', linewidth=1)
ax2.set_title(f'D-alanine (signed vol = {vol_D:+.2f})', color='white', fontweight='bold')
ax2.set_facecolor('#1a1200')
ax2.tick_params(colors='white')
ax2.set_xlim(-2, 2); ax2.set_ylim(-2, 2); ax2.set_zlim(-2, 2)

plt.tight_layout()
plt.savefig('output.png', facecolor='#1a1200', dpi=110, bbox_inches='tight')
plt.close()
print('Chirality detection via signed tetrahedral volume.')
print('Biology: all amino acids in ribosomal synthesis are L-form.')

Click Run to execute the Python code

Code will be executed with Python 3 on the server

8. Simulation: Snail Shell Coiling

Raup's three-parameter model generates four classic shell types in 3D: helical dextral, rare sinistral, Nautilus planispiral and ammonite-like.

Python

script.py47 lines

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

# ======================================================================
# Snail shell coiling: Raup's (1966) three parameters
# W = whorl expansion rate (geometric growth per turn)
# D = distance from coiling axis (translation)
# T = translation per turn (vertical rise)
# ======================================================================

def shell(W=1.5, D=0.0, T=0.5, n_turns=5, n_pts=400):
    theta = np.linspace(0, 2*np.pi*n_turns, n_pts)
    r = D + W**(theta/(2*np.pi))
    x = r*np.cos(theta)
    y = r*np.sin(theta)
    z = T*theta/(2*np.pi)
    return x, y, z

fig = plt.figure(figsize=(16, 10))
fig.patch.set_facecolor('#1a1200')

# Four classic shell types
cases = [
    (1.5, 0.0, 0.5, 'Helical (Turritella-like, dextral)'),
    (1.5, 0.0, -0.5, 'Helical (rare sinistral)'),
    (2.0, 0.0, 0.0, 'Planispiral (Nautilus-like)'),
    (1.3, 1.0, 0.2, 'Ammonite-like'),
]

for i, (W, D, T, title) in enumerate(cases):
    ax = fig.add_subplot(2, 2, i+1, projection='3d')
    x, y, z = shell(W, D, T)
    ax.plot(x, y, z, color='#fbbf24', linewidth=2)
    ax.scatter(x[-1], y[-1], z[-1], color='#f87171', s=50)
    ax.set_title(title, color='white', fontweight='bold', fontsize=11)
    ax.set_facecolor('#1a1200')
    ax.tick_params(colors='white')

plt.tight_layout()
plt.savefig('output.png', facecolor='#1a1200', dpi=110, bbox_inches='tight')
plt.close()
print('Snail shell coiling under Raup three-parameter model.')
print('Dextral vs sinistral: sign of T determines handedness.')
print('~90% of snail species are dextral (right-handed coiling).')

Click Run to execute the Python code

Code will be executed with Python 3 on the server

9. Ecological Consequences of Chirality

Chirality is not just molecular trivia - it has ecological consequences:

Taste & smell: L-carvone smells like spearmint, D-carvone like caraway. Olfactory receptors are chiral.
Toxicology: thalidomide, R-enantiomer is a sedative, S-enantiomer a teratogen - the tragedy of the 1960s. Modern regulations require enantiomeric purity.
Predator escape: dextral snails have enhanced survival against left-jawed snake predators - Hoso (2009) showed Iwasakia snakes cannot swallow sinistral shells.
Speciation: mating between dextral and sinistral snails is mechanically impossible; chirality thus drives reproductive isolation in Satsuma.

10. Symmetry Breaking in Development

All symmetric forms must have been broken from even more symmetric initial conditions. The sphere is the most symmetric 3D object (symmetry group SO(3)); radially symmetric gastrula breaks to bilateral bilaterian by anterior-posterior axis establishment.

Turing mechanism for axis formation

In ciliates and C. elegans zygotes, a polarity axis emerges spontaneously from cortical PAR-3/PAR-6 asymmetry. The initiating bias is provided by the sperm entry point, and self-organizing reaction- diffusion dynamics amplify it into anterior and posterior cortical domains.

Hox code and regionalization

Hox-gene expression breaks the head-tail symmetry. Hox1 (labial) specifies anterior; Hox13 (caudal) specifies posterior. Duplication of Hox clusters during vertebrate evolution expanded body-plan complexity. Today's 9-13 Hox genes per cluster control the identity of every vertebra, rib and limb segment.

9a. Radial Symmetry and Sessile Lifestyle

Radially symmetric body plans are strongly associated with sessile or slow-moving lifestyles: cnidarian polyps anchor to the substrate; echinoderms creep with tube feet equally in all directions; flowers receive pollinators from any angle. The physical advantage is clear - a sessile animal has no preferred direction and benefits from being able to feed, defend, and reproduce through all sides equally.

Conversely, motile animals moving through a viscous or gravitational environment benefit from streamlined, bilateral bodies. Planktonic larvae of echinoderms are bilateral precisely because they swim; the adult radial form emerges after settlement. This developmental narrative is preserved in nearly every echinoderm species.

9b. Applications in Drug Discovery

Understanding biological chirality is foundational to modern pharmacology. Over 56% of top-selling drugs are sold as single enantiomers (“chiral switches”), and the FDA requires separate toxicology studies for each enantiomer of any new chiral drug (post-Thalidomide).

Ibuprofen

Only S-ibuprofen inhibits COX enzymes. R-ibuprofen is inactive but slowly converted to S in vivo. Marketed as racemate since conversion rate is acceptable.

Naproxen

S-naproxen is anti-inflammatory; R-naproxen is hepatotoxic. Must be sold as single S-enantiomer.

Esomeprazole (Nexium)

S-enantiomer of omeprazole with slightly improved pharmacokinetics. AstraZeneca's classic “chiral switch” to extend patent protection.

Escitalopram (Lexapro)

S-citalopram is 30x more potent at serotonin reuptake inhibition than R. Chiral switch product.

10b. Body-Plan Evolution: Echinoderm Ambiguity

Echinoderms (starfish, sea urchins, sand dollars, sea cucumbers) present a puzzle: their larvae are bilaterally symmetric plankton, but adults display pentaradial symmetry. The metamorphic transition between symmetries is unique in the animal kingdom.

Molecular phylogeny places echinoderms firmly within Deuterostomia (sister to chordates). Their radial adult is thus secondarily derived from a bilateral ancestor. Sea cucumbers partially retain bilateral symmetry as adults - hinting at the reverse evolutionary trajectory.

The 5-fold pattern: why five?

No definitive answer. Leading hypotheses: (1) five-fold packing maximizes strength of the skeletal ossicle network (Wolpert 1999); (2) phyllotactic-like mechanisms in the water-vascular system select five primordia on a ring (Hotchkiss 1998); (3) simple genetic wiring of a five-way branching ancestor has been conserved for 540 million years.

11. Group Theory Toolbox

The mathematical classification of symmetries uses point groups. Here are the most common biological groups:

Group	Order	Examples
\(C_1\) (trivial)	1	Asymmetric organs (liver, pancreas)
\(C_2\) (bilateral)	2	Bilateria externally
\(D_3\) (trimery)	6	Monocot flowers, iris
\(D_4\) (tetramery)	8	Cruciferae, jellyfish
\(D_5\) (pentamery)	10	Rose family, starfish
\(D_6\) (hexamery)	12	Honeycomb, snowflakes
\(I_h\) (icosahedral)	120	Adenovirus, radiolaria

The icosahedral symmetry is particularly important for viral capsids because it allows the maximum number of identical subunits to cover a sphere with just 60 positions (Caspar & Klug 1962). Adenovirus uses triangulation number T = 25, meaning 25 × 60 = 1500 hexons plus pentons.

11b. Broken Symmetry in Flower Evolution

Zygomorphic (bilaterally symmetric) flowers evolved repeatedly from actinomorphic (radially symmetric) ancestors. Think of orchids, peas, snapdragons, Lamiaceae - all bilateral. This transition happened independently ~130 times in angiosperm evolution.

The molecular mechanism: ectopic expression of CYCLOIDEA (CYC) family genes in the dorsal petal primordia breaks the 5-fold symmetry. Luo et al. (1996) showed that cyc mutants in Antirrhinum produce radially symmetric peloric flowers. Zygomorphy is pollinator-driven: bilateral flowers can enforce a specific visitor orientation, improving pollen transfer.

12. Open Questions

Why did Earth's biochemistry choose L-amino acids and D-sugars? Contingent accident, parity violation, or selection pressure on early replication fidelity?
Why is pentamerism concentrated in deuterostomes (echinoderms) and plants but rare in protostomes? Developmental constraint or evolutionary legacy?
Can we design homochiral synthetic cells from scratch? Does a second origin of life need to choose the same chirality?
How exactly do nodal cilia propagate signaling vesicles over 30-60 micrometers in the face of Brownian diffusion? Advection or self-promoted transport?

References

Frank, F.C. (1953). On spontaneous asymmetric synthesis. Biochim. Biophys. Acta, 11, 459-463.
Soai, K., Shibata, T., Morioka, H., Choji, K. (1995). Asymmetric autocatalysis and amplification of enantiomeric excess of a chiral molecule. Nature, 378, 767-768.
Raup, D.M. (1966). Geometric analysis of shell coiling: general problems. Journal of Paleontology, 40, 1178-1190.
Nonaka, S. et al. (1998). Randomization of left-right asymmetry due to loss of nodal cilia generating leftward flow of extraembryonic fluid in mice lacking KIF3B motor protein. Cell, 95, 829-837.
Cartwright, J.H.E., Piro, O., Tuval, I. (2004). Fluid-dynamical basis of the embryonic development of left-right asymmetry in vertebrates. PNAS, 101, 7234-7239.
Hoso, M. (2009). Predator-prey coevolution drives the escalation of the right-to-left shift of snail shells. Biol. Lett., 5, 611-614.
Coen, E.S. & Meyerowitz, E.M. (1991). The war of the whorls: genetic interactions controlling flower development. Nature, 353, 31-37.
Bonner, J.T. (2006). Why Size Matters: From Bacteria to Blue Whales. Princeton University Press.
Hegstrom, R.A. & Kondepudi, D.K. (1990). The handedness of the universe. Scientific American, 262(1), 108-115.
Davison, A. et al. (2016). Formin is associated with left-right asymmetry in the pond snail and the frog. Current Biology, 26, 654-660.

←M1: Golden Ratio & Phyllotaxis M3: Fractals & Branching→