The Keller–Segel system is a class of mathematical models describing the collective movement of cells or organisms in response to chemical signals, a process known as chemotaxis. It was first introduced in the 1970s by Evelyn Fox Keller and Lee Segel to explain the aggregation behavior of Dictyostelium discoideum, a slime mold that migrates and forms clusters in response to chemoattractants.
Mathematical formulation The most common form of the Keller–Segel model is a system of coupled nonlinear partial differential equations (PDEs). In its simplest parabolic–parabolic version, it is written as:
∂ u ∂ t = D u Δ u − χ ∇ ⋅ ( u ∇ v ) , ∂ v ∂ t = D v Δ v − α v + β u , {\displaystyle {\begin{aligned}{\frac {\partial u}{\partial t}}&=D_{u}\Delta u-\chi \nabla \cdot (u\nabla v),\\{\frac {\partial v}{\partial t}}&=D_{v}\Delta v-\alpha v+\beta u,\end{aligned}}}
where
u ( x , t ) {\displaystyle u(x,t)} represents the density of cells,
v ( x , t ) {\displaystyle v(x,t)} is the concentration of the chemoattractant,
D u , D v > 0 {\displaystyle D_{u},D_{v}>0} are diffusion coefficients,
χ > 0 {\displaystyle \chi >0} denotes chemotactic sensitivity,
α , β {\displaystyle \alpha ,\beta } are parameters for signal decay and production.
Applications The Keller–Segel system has been widely used in mathematical biology to study:
bacterial chemotaxis, slime mold aggregation, tumor angiogenesis, and ecological population dynamics. It has also served as a prototype model in applied mathematics, illustrating how nonlinear PDEs can capture pattern formation, blow-up phenomena, and self-organization.
Mathematical properties A major research focus concerns the global existence and blow-up of solutions. In two dimensions, the system exhibits a critical mass phenomenon: below a certain threshold of initial cell density, solutions remain globally bounded, while above it, solutions may blow up in finite time, modeling cell aggregation into singular clusters.
Experimental applications The Keller-Segel equations have been used to describe chemotaxis of bacterial populations. To accurately describe experimental data, generalizations or modifications are often required, such as explicitly accounting for the chemotaxis log-sensing (Weber's law) capabilities within a finite dynamic range, non-linear signal production and decay, and hydrodynamic effects at high cell density.
Range expansion When cells are inoculated at the center of a soft-agar Petri dish, they grow and consume nutrients, collectively generating nutrient gradients that point toward uncolonized regions. Following these self-generated gradients leads to rapid range expansion that settles at a constant velocity. Models based on the Keller–Segel equations accurately describe such front propagation, including the shape of the outgoing bacterial front and its expansion rate.
Condensation and pattern formation Bacteria may modify their chemical environment to attract one another and form macroscopic condensates. By secreting attractant molecules, bacteria chemotactically attract neighboring cells, creating positive feedback between signal production and cell accumulation. In a distinct mechanism, bacteria embedded in a uniform repellent environment spontaneously condense into millimeter-sized focal points through a chemotaxis-driven instability: cells locally remove repellent molecules, thereby attracting others. The Keller-Segel equations provide a framework to explain the instability onset, the shape of condensates, and their coalescence.
See also Reaction–diffusion system Chemotaxis
References
