🦠 Cell Growth & Morphogenesis

Turing patterns arise when an activator (which promotes both itself and an inhibitor) diffuses more slowly than the inhibitor — a condition called the "diffusion-driven instability". Small random fluctuations in activator concentration are locally amplified but globally suppressed by the faster-spreading inhibitor, spontaneously selecting a characteristic spatial wavelength. This wavelength scales with √(D_activator/reaction_rate) and determines spot or stripe spacing. Real-world confirmation includes melanocyte patterning in zebrafish (manipulated genetically to switch from stripes to spots), digit spacing in vertebrate limbs, and hair follicle arrangement in mammalian skin.

The circular tissue boundary in this simulation implements a biologically realistic constraint: in real embryos, tissues develop within defined spatial boundaries set by epithelial membranes, extracellular matrix scaffolds, and neighbouring cell layers. Confinement restricts the available patterning wavelengths to multiples that "fit" within the tissue geometry, so spot number and arrangement are partly geometrically determined. Gastruloids — synthetic embryos grown from stem cells in confinement — produce remarkably consistent anterior–posterior axis formation, demonstrating that geometry itself instructs morphogenesis.

The packing fraction measures what proportion of the tissue area is occupied by cell bodies (π×r²/tissue_area). Random packing of equal circles reaches a maximum of ~0.906 (hexagonal close packing), but real biological tissues typically achieve packing fractions of 0.7–0.85 because cells have variable sizes and deform under mechanical pressure. In epithelial monolayers, cells pack into Voronoi-like polygons with a characteristic distribution of 5-, 6-, and 7-sided cells described by the topological Euler relation — an order that emerges purely from mechanical energy minimisation.

Morphogens are signalling molecules that form spatial gradients across a tissue and instruct cells to adopt different fates depending on local concentration — a concentration-dependent cell fate decision. Classic morphogens include Bicoid (head-to-tail axis in Drosophila), Sonic Hedgehog (SHH, digit identity), BMP-4 (dorsal–ventral axis), and Wnt (body axis polarity). Growth factors (such as EGF, FGF, PDGF) primarily stimulate cell proliferation and survival rather than fate specification, though the distinction is blurred because many molecules do both. In this simulator, the V species of Gray–Scott acts as the fate-determining morphogen.

In real tissues, cell growth is limited by contact inhibition — cells that are completely surrounded by neighbours stop dividing, a phenomenon discovered by Abercrombie in the 1950s. The molecular basis involves mechanosensitive signalling through the Hippo pathway: when cells are mechanically crowded, the YAP/TAZ transcription factors are phosphorylated and inactivated, turning off growth genes. Cancer cells often lose contact inhibition — a hallmark of malignancy. In this simulator, crowding is modelled simply by counting neighbours within a fixed radius and reducing growth rate proportionally, capturing the phenomenology without the molecular detail.

The characteristic Turing wavelength λ ≈ 2π√(D_activator/degradation_rate) sets the spacing between spots or stripes. Biological systems adjust this scale to match tissue dimensions: as an embryo grows, the same morphogen parameters produce a proportionally scaled pattern (Wolpert's French Flag Problem). Recent work has shown that some morphogen gradients scale with tissue size through feedback loops — the activator drives its own production but also stimulates inhibitor production, and the inhibitor's degradation rate adjusts with tissue volume. This active scaling ensures robustness of pattern formation across a range of embryo sizes.

Yes — there is growing evidence that mechanical forces are as important as chemical signals in morphogenesis. Differential surface tension between cell populations drives cell sorting (Steinberg's differential adhesion hypothesis, 1963); tissue folding in gastrulation and neurulation is driven by apical constriction where cells contract their apical surfaces; and compressive stresses in the developing gut cause it to buckle into characteristic looped shapes. Mechanotransduction — the conversion of mechanical forces into gene expression changes — links the chemical and mechanical aspects of morphogenesis, with the cytoskeleton serving as the key signal transducer.

Organoids are three-dimensional miniature organ models grown from pluripotent or adult stem cells in a gel matrix, pioneered by Hans Clevers' group for the intestine in 2009. They spontaneously self-organise into structures resembling the tissue of origin — complete with crypt–villus architecture in intestinal organoids, or alveoli in lung organoids — driven by the same Turing-type morphogen gradients modelled here. Organoids are revolutionising drug testing, personalised medicine (patient-derived cancer organoids predict chemotherapy response), and disease modelling for conditions including COVID-19 lung infection and cystic fibrosis.