Quantitative modeling of Escherichia coli chemotactic motion in environments varying in space and time

Enhanced chemotaxis efficiency of Escherichia coli in viscoelastic solutions

Bacteria inhabit complex environments rich in macromolecular polymers that exhibit viscoelastic properties. While the influence of viscoelasticity on bacterial swimming is recognized, its impact on chemotaxis-a critical behavior for bacterial survival and colonization-remains elusive. In this study, we employed a microfluidic device to establish attractant gradients and observed the chemotactic behavior of Escherichia coli in both viscoelastic solutions containing carboxymethyl cellulose (CMC) and Newtonian buffers. Our results reveal that E. coli demonstrates markedly enhanced chemotactic efficiency in viscoelastic media. Notably, bacteria achieved faster migration velocities and higher steady-state accumulation in areas with higher attractant concentrations compared to those in Newtonian conditions. Through 3D tracking, we determined that changes in bulk motility parameters alone do not account for the observed enhancements. Further investigations through theoretical analysis and stochastic simulations suggested that the main enhancement mechanisms are mitigation of surface hydrodynamic hindrance resulting from solid surfaces commonly present in bacterial habitats, and the induction of a lifting force in viscoelastic solutions. These findings highlight the significant role of the rheological properties of bacterial habitats in shaping their chemotactic strategies, offering deeper insights into bacterial adaptive mechanisms in both natural and clinical settings.

Potassium-mediated bacterial chemotactic response

Bacteria in biofilms secrete potassium ions to attract free swimming cells. However, the basis of chemotaxis to potassium remains poorly understood. Here, using a microfluidic device, we found that Escherichia coli can rapidly accumulate in regions of high potassium concentration on the order of millimoles. Using a bead assay, we measured the dynamic response of individual flagellar motors to stepwise changes in potassium concentration, finding that the response resulted from the chemotaxis signaling pathway. To characterize the chemotactic response to potassium, we measured the dose-response curve and adaptation kinetics via an Förster resonance energy transfer (FRET) assay, finding that the chemotaxis pathway exhibited a sensitive response and fast adaptation to potassium. We further found that the two major chemoreceptors Tar and Tsr respond differently to potassium. Tar receptors exhibit a biphasic response, whereas Tsr receptors respond to potassium as an attractant. These different responses were consistent with the responses of the two receptors to intracellular pH changes. The sensitive response and fast adaptation allow bacteria to sense and localize small changes in potassium concentration. The differential responses of Tar and Tsr receptors to potassium suggest that cells at different growth stages respond differently to potassium and may have different requirements for potassium.

Strong chemotaxis by marine bacteria towards polysaccharides is enhanced by the abundant organosulfur compound DMSP

The ability of marine bacteria to direct their movement in response to chemical gradients influences inter-species interactions, nutrient turnover, and ecosystem productivity. While many bacteria are chemotactic towards small metabolites, marine organic matter is predominantly composed of large molecules and polymers. Yet, the signalling role of these large molecules is largely unknown. Using in situ and laboratory-based chemotaxis assays, we show that marine bacteria are strongly attracted to the abundant algal polysaccharides laminarin and alginate. Unexpectedly, these polysaccharides elicited stronger chemoattraction than their oligo- and monosaccharide constituents. Furthermore, chemotaxis towards laminarin was strongly enhanced by dimethylsulfoniopropionate (DMSP), another ubiquitous algal-derived metabolite. Our results indicate that DMSP acts as a methyl donor for marine bacteria, increasing their gradient detection capacity and facilitating their access to polysaccharide patches. We demonstrate that marine bacteria are capable of strong chemotaxis towards large soluble polysaccharides and uncover a new ecological role for DMSP in enhancing this attraction. These navigation behaviours may contribute to the rapid turnover of polymers in the ocean, with important consequences for marine carbon cycling.

Tethered particle motion of the adaptation enzyme CheR in bacterial chemotaxis

Bacteria perform chemotactic adaptation by sequential modification of multiple modifiable sites on chemoreceptors through stochastic action of tethered adaptation enzymes (CheR and CheB). To study the molecular kinetics of this process, we measured the response to different concentrations of MeAsp for the Tar-only Escherichia coli strain. We found a strong dependence of the methylation rate on the methylation level and established a new mechanism of adaptation kinetics due to tethered particle motion of the methylation enzyme CheR. Experiments with various lengths of the C-terminal flexible chain in the Tar receptor further validated this mechanism. The tethered particle motion resulted in a CheR concentration gradient that ensures encounter-rate matching of the sequential modifiable sites. An analytical model of multisite catalytic reaction showed that this enables robustness of methylation to fluctuations in receptor activity or cell-to-cell variations in the expression of adaptation enzymes and reduces the variation in methylation level among individual receptors.

Two-state swimming: Strategy and survival of a model bacterial predator in response to environmental cues

Bdellovibrio bacteriovorus is a predatory bacterium preying upon Gram-negative bacteria. As such, B. bacteriovorus has the potential to control antibiotic-resistant pathogens and biofilm populations. To survive and reproduce, B. bacteriovorus must locate and infect a host cell. However, in the temporary absence of prey, it is largely unknown how B. bacteriovorus modulate their motility patterns in response to physical or chemical environmental cues to optimize their energy expenditure. To investigate B. bacteriovorus' predation strategy, we track and quantify their motion by measuring speed distributions as a function of starvation time. While an initial unimodal speed distribution relaxing to one for pure diffusion at long times may be expected, instead we observe a bimodal speed distribution with one mode centered around that expected from diffusion and the other centered at higher speeds. What is more, for an increasing amount of time over which B. bacteriovorus is starved, we observe a progressive reweighting from the active swimming state to an apparent diffusive state in the speed distribution. Distributions of trajectory-averaged speeds for B. bacteriovorus are largely unimodal, indicating switching between a faster swim speed and an apparent diffusive state within individual observed trajectories rather than there being distinct active swimming and apparent diffusive populations. We also find that B. bacteriovorus' apparent diffusive state is not merely caused by the diffusion of inviable bacteria as subsequent spiking experiments show that bacteria can be resuscitated and bimodality restored. Indeed, starved B. bacteriovorus may modulate the frequency and duration of active swimming as a means of balancing energy consumption and procurement. Our results thus point to a reweighting of the swimming frequency on a trajectory basis rather than a population level basis.

Agent-Based Modeling of Microbial Communities

Microbial communities are complex living systems that populate the planet with diverse functions and are increasingly harnessed for practical human needs. To deepen the fundamental understanding of their organization and functioning as well as to facilitate their engineering for applications, mathematical modeling has played an increasingly important role. Agent-based models represent a class of powerful quantitative frameworks for investigating microbial communities because of their individualistic nature in describing cells, mechanistic characterization of molecular and cellular processes, and intrinsic ability to produce emergent system properties. This review presents a comprehensive overview of recent advances in agent-based modeling of microbial communities. It surveys the state-of-the-art algorithms employed to simulate intracellular biomolecular events, single-cell behaviors, intercellular interactions, and interactions between cells and their environments that collectively serve as the driving forces of community behaviors. It also highlights three lines of applications of agent-based modeling, namely, the elucidation of microbial range expansion and colony ecology, the design of synthetic gene circuits and microbial populations for desired behaviors, and the characterization of biofilm formation and dispersal. The review concludes with a discussion of existing challenges, including the computational cost of the modeling, and potential mitigation strategies.

Individual-Based Modeling of Spatial Dynamics of Chemotactic Microbial Populations

One important direction of synthetic biology is to establish desired spatial structures from microbial populations. Underlying this structural development process are different driving factors, among which bacterial motility and chemotaxis serve as a major force. Here, we present an individual-based, biophysical computational framework for mechanistic and multiscale simulation of the spatiotemporal dynamics of motile and chemotactic microbial populations. The framework integrates cellular movement with spatial population growth, mechanical and chemical cellular interactions, and intracellular molecular kinetics. It is validated by a statistical comparison of single-cell chemotaxis simulations with reported experiments. The framework successfully captures colony range expansion of growing isogenic populations and also reveals chemotaxis-modulated, spatial patterns of a two-species amensal community. Partial differential equation-based models subsequently validate these simulation findings. This study provides a versatile computational tool to uncover the fundamentals of microbial spatial ecology as well as to facilitate the design of synthetic consortia for desired spatial patterns.

Numerical Study of the Volcano Effect in Chemotactic Aggregation Based on a Kinetic Transport Equation with Non-instantaneous Tumbling

Aggregation of chemotactic bacteria under a unimodal distribution of chemical cues was investigated by Monte Carlo (MC) simulation based on a kinetic transport equation, which considers an internal adaptation dynamics as well as a finite tumbling duration. It was found that there exist two different regimes of the adaptation time, between which the effect of the adaptation time on the aggregation behavior is reversed; that is, when the adaptation time is as small as the running duration, the aggregation becomes increasingly steeper as the adaptation time increases, while, when the adaptation time is as large as the diffusion time of the population density, the aggregation becomes more diffusive as the adaptation time increases. Moreover, the aggregation profile becomes bimodal (volcano) at the large adaptation-time regime when the tumbling duration is sufficiently large while it is always unimodal at the small adaptation-time regime. A remarkable result of this study is the identification of the parameter regime and scaling for the volcano effect. That is, by comparing the results of MC simulations to the continuum-limit models obtained at each of the small and large adaptation-time scalings, it is clarified that the volcano effect arises due to the coupling of diffusion, adaptation, and finite tumbling duration, which occurs at the large adaptation-time scaling.

Bacterial chemotaxis to saccharides is governed by a trade-off between sensing and uptake

To swim up gradients of nutrients, E. coli senses nutrient concentrations within its periplasm. For small nutrient molecules, periplasmic concentrations typically match extracellular concentrations. However, this is not necessarily the case for saccharides, such as maltose, which are transported into the periplasm via a specific porin. Previous observations have shown that, under various conditions, E. coli limits maltoporin abundance so that, for extracellular micromolar concentrations of maltose, there are predicted to be only nanomolar concentrations of free maltose in the periplasm. Thus, in the micromolar regime, the total uptake of maltose from the external environment into the cytoplasm is limited not by the abundance of cytoplasmic transport proteins but by the abundance of maltoporins. Here we present results from experiments and modeling suggesting that this porin-limited transport enables E. coli to sense micromolar gradients of maltose despite having a high-affinity ABC transport system that is saturated at these micromolar levels. We used microfluidic assays to study chemotaxis of E. coli in various gradients of maltose and methyl-aspartate and leveraged our experimental observations to develop a mechanistic transport-and-sensing chemotaxis model. Incorporating this model into agent-based simulations, we discover a trade-off between uptake and sensing: although high-affinity transport enables higher uptake rates at low nutrient concentrations, it severely limits the range of dynamic sensing. We thus propose that E. coli may limit periplasmic uptake to increase its chemotactic sensitivity, enabling it to use maltose as an environmental cue.

The Effect of the Second Messenger c-di-GMP on Bacterial Chemotaxis in Escherichia coli

c-di-GMP is a ubiquitous bacterial second messenger that plays a central regulatory role in diverse biological processes. c-di-GMP was known to regulate chemotaxis in multiple bacterial species, but its effect on Escherichia coli chemotaxis remained unclear. As an effector of c-di-GMP in E. coli, YcgR when bound with c-di-GMP interacts with the flagellar motor to reduce its speed and its probability of rotating clockwise (CW bias). Here, we found that a significant fraction of the c-di-GMP::YcgR dynamically exchange between the motor and the cytosol. Through fluorescent measurements, we found that there was no competitive binding between the chemotaxis response regulator CheY-P and c-di-GMP::YcgR to the motor. To test the influence of elevated c-di-GMP levels on the chemotaxis pathway, we measured the chemotactic responses of E. coli cells using a FRET assay, finding that elevated c-di-GMP levels had no effect on the upstream part of chemotaxis pathway down to the level of CheY-P concentration. This suggested that the possible effect of elevated c-di-GMP levels on chemotactic motion was through regulation of motor speed and CW bias. Using stochastic simulations of chemotactic swimming, we showed that the effects of reducing motor speed and decreasing CW bias on chemotactic drift velocity are compensating for each other, resulting in minimal effect of elevated c-di-GMP levels on E. coli chemotaxis. Therefore, elevated c-di-GMP levels promote the transition from motile to sedentary forms of bacterial life by reducing the bacterial swimming speed and CW bias, while still maintaining a nearly intact chemotaxis capability in E. coli. IMPORTANCE The ubiquitous bacterial second messenger c-di-GMP was known to regulate chemotaxis in many bacterial species, but its effect on E. coli chemotaxis was unclear. Here we studied the effect of elevated c-di-GMP levels on chemotaxis in E. coli. We found that the binding of c-di-GMP::YcgR (its effector) and the chemotaxis response regulator CheY-P to the flagellar motor are noncompetitive, and elevated c-di-GMP levels do not affect the upstream part of the chemotaxis pathway down to the level of CheY-P concentration. Elevated c-di-GMP levels exert direct effects on the flagellar motor by reducing its speed and CW bias, but the resulting effects on chemotaxis performance are compensating for each other. Our findings here showed that elevated c-di-GMP levels maintain a nearly intact chemotaxis capability when promoting the transition from motile to sedentary forms of bacterial life in E. coli.

A new mode of swimming in singly flagellated Pseudomonas aeruginosa

SignificanceThe monotrichous Pseudomonas aeruginosa was usually thought to swim in a pattern of "run and reverse" (possibly with pauses in between), where straight runs alternated with reverses with angular changes of swimming direction near 180°. Here, by simultaneously tracking the cell swimming and the morphology of its flagellum, we discovered a swimming mode in P. aeruginosa-the wrap mode, during which the flagellar filament wrapped around the cell body and induced large fluctuation of the body orientation. The wrap mode randomized swimming direction, resulting in a broad distribution of angular changes over 0 to 180° with a peak near 90°. This allowed the bacterium to explore the environment more efficiently, which we confirmed by stochastic simulations of P. aeruginosa chemotaxis.

Effect of receptor cooperativity on methylation dynamics in bacterial chemotaxis with weak and strong gradient

We study methylation dynamics of the chemoreceptors as an Escherichia coli cell moves around in a spatially varying chemoattractant environment. We consider attractant concentration with strong and weak spatial gradient. During the uphill and downhill motion of the cell along the gradient, we measure the temporal variation of average methylation level of the receptor clusters. Our numerical simulations we show that the methylation dynamics depends sensitively on the size of the receptor clusters and also on the strength of the gradient. At short times after the beginning of a run, the methylation dynamics is mainly controlled by short runs which are generally associated with high receptor activity. This results in demethylation at short times. But for intermediate or large times, long runs play an important role and depending on receptor cooperativity or gradient strength, the qualitative variation of methylation can be completely different in this time regime. For weak gradient, both for uphill and downhill runs, after the initial demethylation, we find methylation level increases steadily with time for all cluster sizes. Similar qualitative behavior is observed for strong gradient during uphill runs as well. However, the methylation dynamics for downhill runs in strong gradient show highly nontrivial dependence on the receptor cluster size. We explain this behavior as a result of interplay between the sensing and adaptation modules of the signaling network.

A Mathematical Description of Bacterial Chemotaxis in Response to Two Stimuli

Bacteria are often exposed to multiple stimuli in complex environments, and their efficient chemotactic decisions are critical to survive and grow in their native environments. Bacterial responses to the environmental stimuli depend on the ratio of their corresponding chemoreceptors. By incorporating the signaling machinery of individual cells, we analyze the collective motion of a population of Escherichia coli bacteria in response to two stimuli, mainly serine and methyl-aspartate (MeAsp), in a one-dimensional and a two-dimensional environment, which is inspired by experimental results in Y. Kalinin et al., J. Bacteriol. 192(7):1796-1800, 2010. Under suitable conditions, we show that if the ratio of the main chemoreceptors of individual cells, namely Tar/Tsr, is less than a specific threshold, the bacteria move to the gradient of serine, and if the ratio is greater than the threshold, the group of bacteria moves toward the gradient of MeAsp. Finally, we examine the theory with Monte Carlo agent-based simulations and verify that our results qualitatively agree well with the experimental results in Y. Kalinin et al. (2010).

Spatial modulation of individual behaviors enables an ordered structure of diverse phenotypes during bacterial group migration

Coordination of diverse individuals often requires sophisticated communications and high-order computational abilities. Microbial populations can exhibit diverse individualistic behaviors, and yet can engage in collective migratory patterns with a spatially sorted arrangement of phenotypes. However, it is unclear how such spatially sorted patterns emerge from diverse individuals without complex computational abilities. Here, by investigating the single-cell trajectories during group migration, we discovered that, despite the constant migrating speed of a group, the drift velocities of individual bacteria decrease from the back to the front. With a Langevin-type modeling framework, we showed that this decreasing profile of drift velocities implies the spatial modulation of individual run-and-tumble random motions, and enables the bacterial population to migrate as a pushed wave front. Theoretical analysis and stochastic simulations further predicted that the pushed wave front can help a diverse population to stay in a tight group, while diverse individuals perform the same type of mean reverting processes around centers orderly aligned by their chemotactic abilities. This mechanism about the emergence of orderly collective migration from diverse individuals is experimentally demonstrated by titration of bacterial chemoreceptor abundance. These results reveal a simple computational principle for emergent ordered behaviors from heterogeneous individuals.

Organisms living in large groups often have to move together in order to navigate, forage for food, and increase their roaming range. Such groups are often made up of distinct individuals that must integrate their different behaviors in order to migrate in the same direction at a similar pace. For instance, for the bacteria Escherichia coli to travel as a condensed group, they must coordinate their response to a set of chemical signals called chemoattractants that tell them where to go.

The chemoattractants surrounding the bacteria are unequally distributed so that there is more of them at the front than the back of the group. During migration, each bacterium moves towards this concentration gradient in a distinct way, spontaneously rotating its direction in a ‘run-and-tumble’ motion that guides it towards areas where there are high levels of these chemical signals. In addition to this variability, how well individual bacteria are able to swim up the gradient also differs within the population. Bacteria that are better at sensing the chemoattractant gradient are placed at the front of the group, while those that are worst are shifted towards the back. This spatial arrangement is thought to help the bacteria migrate together. But how E. coli organize themselves in to this pattern is unclear, especially as they cannot communicate directly with one another and display such diverse, randomized behaviors.

To help answer this question, Bai, He et al. discovered a general principle that describes how single bacterial cells move within a group. The results showed that E. coli alter their run-and-tumble motion depending on where they reside within the population: individuals at the rear drift faster so they can catch up with the group, while those leading the group drift slower to draw themselves back. This ‘reversion behavior’ allows the migrating bacteria to travel at a constant speed around a mean position relative to the group.

A cell’s drifting speed is determined by how well it moves towards the chemoattractant and its response to the concentration gradient. As a result, the mean position around which the bacterium accelerates or deaccelerates will vary depending on how sensitive it is to the chemoattractant gradient. The E. coli therefore spatially arrange themselves so that the more sensitive bacteria are located at the front of the group where the gradient is shallower; and cells that are less sensitive are located towards the back where the gradient is steeper.

These findings suggest a general principle for how bacteria form ordered patterns whilst migrating as a collective group. This behavior could also apply to other populations of distinct individuals, such as ants following a trail or flocks of birds migrating in between seasons.

Effects of internal dynamics on chemotactic aggregation of bacteria

The effects of internal adaptation dynamics on the self-organized aggregation of chemotactic bacteria are investigated by Monte Carlo (MC) simulations based on a two-stream kinetic transport equation coupled with a reaction-diffusion equation of the chemoattractant that bacteria produce. A remarkable finding is a nonmonotonic behavior of the peak aggregation density with respect to the adaptation time; more specifically, aggregation is the most enhanced when the adaptation time is comparable to or moderately larger than the mean run time of bacteria. Another curious observation is the formation of a trapezoidal aggregation profile occurring at a very large adaptation time, where the biased motion of individual cells is rather hindered at the plateau regimes due to the boundedness of the tumbling frequency modulation. Asymptotic analysis of the kinetic transport system is also carried out, and a novel asymptotic equation is obtained at the large adaptation-time regime while the Keller-Segel type equations are obtained when the adaptation time is moderate. Numerical comparison of the asymptotic equations with MC results clarifies that trapezoidal aggregation is well described by the novel asymptotic equation, and the nonmonotonic behavior of the peak aggregation density is interpreted as the transient of the asymptotic solutions between different adaptation time regimes.

Individual based models exhibiting Lévy-flight type movement induced by intracellular noise

We use an individual-based model and its associated kinetic equation to study the generation of long jumps in the motion of E. coli. These models relate the run-and-tumble process to the intracellular reaction where the intrinsic noise plays a central role. Compared with previous work in Perthame et al. (Z Angew Math Phys 69(3):1-15, 2018), in which the parametric assumptions are mainly targeted for mathematical convenience but not well-suited for numerical simulations or comparison with experimental results, our current paper makes use of biologically meaningful pathways and tumbling kernels. The main contribution of this current work is bridging the gap between the theoretical results and experimentally available data. Some particular forms of how the tumbling frequency depends on the internal variable are proposed. Moreover, we propose two individual-based models, one for the tumbling frequency and the other for the receptor activity, and perform numerical simulations. Power-law decay of the run length distribution, which corresponds to Lévy-type motions, is observed in our numerical results. The particular decay rate agrees quantitatively with the analytical result.

Connection between the Bacterial Chemotactic Network and Optimal Filtering

The chemotactic network of Escherichia coli has been studied extensively both biophysically and information theoretically. Nevertheless, connection between these two aspects is still elusive. In this work, we report such a connection. We derive an optimal filtering dynamics under the assumption that E. coli's sensory system optimally infers the binary information whether it is swimming up or down along an exponential ligand gradient from noisy sensory signals. Then we show that a standard biochemical model of the chemotactic network is mathematically equivalent to this information-theoretically optimal dynamics. Moreover, we demonstrate that an experimentally observed nonlinear response relation can be reproduced from the optimal dynamics. These results suggest that the biochemical network of E. coli chemotaxis is designed to optimally extract the binary information along an exponential gradient in a noisy condition.

Collective motion enhances chemotaxis in a two-dimensional bacterial swarm

In a dilute liquid environment in which cell-cell interaction is negligible, flagellated bacteria, such as Escherichia coli, perform chemotaxis by biased random walks alternating between run-and-tumble. In a two-dimensional crowded environment, such as a bacterial swarm, the typical behavior of run-and-tumble is absent, and this raises the question whether and how bacteria can perform chemotaxis in a swarm. Here, by examining the chemotactic behavior as a function of the cell density, we showed that chemotaxis is surprisingly enhanced because of cell crowding in a bacterial swarm, and this enhancement is correlated with increase in the degree of cell body alignment. Cells tend to form clusters that move collectively in a swarm with increased effective run length, and we showed analytically that this resulted in increased drift velocity toward attractants. We also explained the enhancement by stochastically simulating bacterial chemotaxis in a swarm. We found that cell crowding in a swarm enhances chemotaxis if the cell-cell interactions used in the simulation induce cell-cell alignment, but it impedes chemotaxis if the interactions are collisions that randomize cell moving direction. Therefore, collective motion in a bacterial swarm enhances chemotaxis.

Theory for the optimal detection of time-varying signals in cellular sensing systems

Living cells often need to measure chemical concentrations that vary in time, yet how accurately they can do so is poorly understood. Here, we present a theory that fully specifies, without any adjustable parameters, the optimal design of a canonical sensing system in terms of two elementary design principles: (1) there exists an optimal integration time, which is determined by the input statistics and the number of receptors; and (2) in the optimally designed system, the number of independent concentration measurements as set by the number of receptors and the optimal integration time equals the number of readout molecules that store these measurements and equals the work to store these measurements reliably; no resource is then in excess and hence wasted. Applying our theory to the Escherichia coli chemotaxis system indicates that its integration time is not only optimal for sensing shallow gradients but also necessary to enable navigation in these gradients.

Direct Mapping from Intracellular Chemotaxis Signaling to Single-Cell Swimming Behavior

Bacterial chemotaxis allows bacteria to sense the chemical environment and modulate their swimming behavior accordingly. Although the intracellular chemotaxis signaling pathway has been studied extensively, experimental studies are still lacking that could provide direct link from the pathway output (the intracellular concentration of the phosphorylated form of the response regulator phosphorylated CheY (CheY-P)) to single-cell swimming behavior. Here, we measured the swimming behavior of individual Escherichia coli cells while simultaneously detecting the intracellular CheY-P concentration, thereby providing a direct relationship between the intracellular CheY-P concentration and the single-cell run-and-tumble behavior. The measured relationship is consistent with the ultrasensitivity of the motor switch and a "veto model" that describes the interaction among individual flagella, although contribution from the voting mechanism could not be ruled out.

A Multi-Scale Approach to Modeling E. coli Chemotaxis

The degree to which we can understand the multi-scale organization of cellular life is tied to how well our models can represent this organization and the processes that drive its evolution. This paper uses Vivarium-an engine for composing heterogeneous computational biology models into integrated, multi-scale simulations. Vivarium's approach is demonstrated by combining several sub-models of biophysical processes into a model of chemotactic E. coli that exchange molecules with their environment, express the genes required for chemotaxis, swim, grow, and divide. This model is developed incrementally, highlighting cross-compartment mechanisms that link E. coli to its environment, with models for: (1) metabolism and transport, with transport moving nutrients across the membrane boundary and metabolism converting them to useful metabolites, (2) transcription, translation, complexation, and degradation, with stochastic mechanisms that read real gene sequence data and consume base pairs and ATP to make proteins and complexes, and (3) the activity of flagella and chemoreceptors, which together support navigation in the environment.

From indication to decision: A hierarchical approach to model the chemotactic behavior of Escherichia coli

Reducing the complex behavior of living entities to its underlying physical and chemical processes is a formidable task in biology. Complex behaviors can be characterized as decision making: the ability to process the incoming information via an intracellular network and act upon this information to choose appropriate strategies. Motility is one such behavior that has been the focus many modeling efforts in the past. Our aim is to reduce the chemotactic behavior in Escherichia coli to its molecular constituents in order to paint a comprehensive and end-to-end picture of this intricate behavior. We utilize a hierarchical approach, consisting of three layers, to achieve this goal: at the first level, chemical reactions involved in chemotaxis are simulated. In the second level, the chemical reactions give rise to the mechanical movement of six independent flagella. At the last layer, the two lower layers are combined to allow a digital bacterium to receive information from its environment and swim through it with verve. Our results are in concert with the experimental studies concerning the motility of E.coli cells. In addition, we show that our detailed model of chemotaxis is reducible to a non-homogeneous Markov process.

Chemotaxis in external fields: Simulations for active magnetic biological matter

The movement of microswimmers is often described by active Brownian particle models. Here we introduce a variant of these models with several internal states of the swimmer to describe stochastic strategies for directional swimming such as run and tumble or run and reverse that are used by microorganisms for chemotaxis. The model includes a mechanism to generate a directional bias for chemotaxis and interactions with external fields (e.g., gravity, magnetic field, fluid flow) that impose forces or torques on the swimmer. We show how this modified model can be applied to various scenarios: First, the run and tumble motion of E. coli is used to establish a paradigm for chemotaxis and investigate how it is affected by external forces. Then, we study magneto-aerotaxis in magnetotactic bacteria, which is biased not only by an oxygen gradient towards a preferred concentration, but also by magnetic fields, which exert a torque on an intracellular chain of magnets. We study the competition of magnetic alignment with active reorientation and show that the magnetic orientation can improve chemotaxis and thereby provide an advantage to the bacteria, even at rather large inclination angles of the magnetic field relative to the oxygen gradient, a case reminiscent of what is expected for the bacteria at or close to the equator. The highest gain in chemotactic velocity is obtained for run and tumble with a magnetic field parallel to the gradient, but in general a mechanism for reverse motion is necessary to swim against the magnetic field and a run and reverse strategy is more advantageous in the presence of a magnetic torque. This finding is consistent with observations that the dominant mode of directional changes in magnetotactic bacteria is reversal rather than tumbles. Moreover, it provides guidance for the design of future magnetic biohybrid swimmers.

In this paper, we propose a modified Active Brownian particle model to describe bacterial swimming behavior under the influence of external forces and torques, in particular of a magnetic torque. This type of interaction is particularly important for magnetic biohybrids (i.e. motile bacteria coupled to a synthetic magnetic component) and for magnetotactic bacteria (i.e. bacteria with a natural intracellular magnetic chain), which perform chemotaxis to swim along chemical gradients, but are also directed by an external magnetic field. The model allows us to investigate the benefits and disadvantages of such coupling between two different directionality mechanisms. In particular we show that the magnetic torque can speed chemotaxis up in some conditions, while it can hinder it in other cases. In addition to an understanding of the swimming strategies of naturally magnetotactic organisms, the results may guide the design of future biomedical devices.

Run-and-tumble motion with steplike responses to a stochastic input

We study a simple run-and-tumble random walk whose switching frequencies between run mode and tumble mode depend on a stochastic signal. We consider a particularly sharp, steplike dependence, where the run-to-tumble switching probability jumps from zero to one as the signal crosses a particular value (say y_{1}) from below. Similarly, tumble-to-run switching probability also shows a jump like this as the signal crosses another value (y_{2}<y_{1}) from above. We are interested in characterizing the effect of signaling noise on the long-time behavior of the random walker. We consider two different time-evolutions of the stochastic signal. In one case, the signal dynamics is an independent stochastic process and does not depend on the run-and-tumble motion. In this case we can analytically calculate the mean value and the complete distribution function of the run duration and tumble duration. In the second case, we assume that the signal dynamics is influenced by the spatial location of the random walker. For this system, we numerically measure the steady state position distribution of the random walker. We discuss some similarities and differences between our system and Escherichia coli chemotaxis, which is another well-known run-and-tumble motion encountered in nature.

Optimal methylation noise for best chemotactic performance of E. coli

In response to a concentration gradient of chemoattractant, E. coli bacterium modulates the rotational bias of flagellar motors which control its run-and-tumble motion, to migrate towards regions of high chemoattractant concentration. Presence of stochastic noise in the biochemical pathway of the cell has important consequences on the switching mechanism of motor bias, which in turn affects the runs and tumbles of the cell in a significant way. We model the intracellular reaction network in terms of coupled time evolution of three stochastic variables-kinase activity, methylation level, and CheY-P protein level-and study the effect of methylation noise on the chemotactic performance of the cell. In presence of a spatially varying nutrient concentration profile, a good chemotactic performance allows the cell to climb up the concentration gradient quickly and localize in the nutrient-rich regions in the long time limit. Our simulations show that the best performance is obtained at an optimal noise strength. While it is expected that chemotaxis will be weaker for very large noise, it is counterintuitive that the performance worsens even when noise level falls below a certain value. We explain this striking result by detailed analysis of CheY-P protein level statistics for different noise strengths. We show that when the CheY-P level falls below a certain (noise-dependent) threshold the cell tends to move down the concentration gradient of the nutrient, which has a detrimental effect on its chemotactic response. This threshold value decreases as noise is increased, and this effect is responsible for noise-induced enhancement of chemotactic performance. In a harsh chemical environment, when the nutrient degrades with time, the amount of nutrient intercepted by the cell trajectory is an effective performance criterion. In this case also, depending on the nutrient lifetime, we find an optimum noise strength when the performance is at its best.

Behavioral Variability and Phenotypic Diversity in Bacterial Chemotaxis

Living cells detect and process external signals using signaling pathways that are affected by random fluctuations. These variations cause the behavior of individual cells to fluctuate over time (behavioral variability) and generate phenotypic differences between genetically identical individuals (phenotypic diversity). These two noise sources reduce our ability to predict biological behavior because they diversify cellular responses to identical signals. Here, we review recent experimental and theoretical advances in understanding the mechanistic origin and functional consequences of such variation in Escherichia coli chemotaxis-a well-understood model of signal transduction and behavior. After briefly summarizing the architecture and logic of the chemotaxis system, we discuss determinants of behavior and chemotactic performance of individual cells. Then, we review how cell-to-cell differences in protein abundance map onto differences in individual chemotactic abilities and how phenotypic variability affects the performance of the population. We conclude with open questions to be addressed by future research.

A computational model for how cells choose temporal or spatial sensing during chemotaxis

Cell size is thought to play an important role in choosing between temporal and spatial sensing in chemotaxis. Large cells are thought to use spatial sensing due to large chemical difference at its ends whereas small cells are incapable of spatial sensing due to rapid homogenization of proteins within the cell. However, small cells have been found to polarize and large cells like sperm cells undergo temporal sensing. Thus, it remains an open question what exactly governs spatial versus temporal sensing. Here, we identify the factors that determines sensing choices through mathematical modeling of chemotactic circuits. Comprehensive computational search of three-node signaling circuits has identified the negative integral feedback (NFB) and incoherent feedforward (IFF) circuits as capable of adaptation, an important property for chemotaxis. Cells are modeled as one-dimensional circular system consisting of diffusible activator, inactivator and output proteins, traveling across a chemical gradient. From our simulations, we find that sensing outcomes are similar for NFB or IFF circuits. Rather than cell size, the relevant parameters are the 1) ratio of cell speed to the product of cell diameter and rate of signaling, 2) diffusivity of the output protein and 3) ratio of the diffusivities of the activator to inactivator protein. Spatial sensing is favored when all three parameters are low. This corresponds to a cell moving slower than the time it takes for signaling to propagate across the cell diameter, has an output protein that is polarizable and has a local-excitation global-inhibition system to amplify the chemical gradient. Temporal sensing is favored otherwise. We also find that temporal sensing is more robust to noise. By performing extensive literature search, we find that our prediction agrees with observation in a wide range of species and cell types ranging from E. coli to human Fibroblast cells and propose that our result is universally applicable.

Unicellular organisms and other single cells often have to migrate towards food sources or away from predators by sensing chemicals present in the environment. There are two ways for a cell to sense these external chemicals: temporal sensing, where the cell senses the external chemical at two different time points after it has moved through a certain distance, or spatial sensing, where the cell senses the external chemical at two different locations on its cellular surface (e.g., the front and rear of the cell) simultaneously. It has been thought that small unicellular organisms employ temporal sensing as their small size prohibits sensing at two different locations on the cellular surface. Using computational modeling, we find that the choice between temporal and spatial sensing is determined by the ratio of cell velocity to the product of cell diameter and rate of signaling, as well as the diffusivities of the signaling proteins. Predictions from our model agree with experimental observations over a wide range of cells, where a fast-moving, small cell performs better comparing the chemoattractant at different times in its trajectory; whereas, a slow-moving, big cell performs better by comparing the chemoattractant concentration at its two ends.

Adaptation of Living Systems

Adaptation refers to the biological phenomenon where living systems change their internal states in response to changes in their environments in order to maintain certain key functions critical for their survival and fitness. Adaptation is one of the most ubiquitous and arguably one of the most fundamental properties of living systems. It occurs throughout all biological scales, from adaptation of populations of species over evolutionary time to adaptation of a single cell to different environmental stresses during its life span. In this article, we review some of the recent progress made in understanding molecular mechanisms of cellular level adaptation. We take the minimalist (or the physicist) approach and study the simplest systems that exhibit generic adaptive behaviors. We focus on understanding the basic biochemical interaction networks in living matter that are responsible for adaptation dynamics. By combining theoretical modeling with quantitative experimentation, we demonstrate universal features in adaptation as well as important differences in different cellular systems, including chemotaxis in bacterium cells (Escherichia coli) and eukaryotic cells (Dictyostelium). Future work in extending the modeling framework to study adaptation in more complex systems such as sensory neurons are discussed.

Statistical signatures of a targeted search by bacteria

Chemoattractant gradients are rarely well-controlled in nature and recent attention has turned to bacterial chemotaxis toward typical bacterial food sources such as food patches or even bacterial prey. In environments with localized food sources reminiscent of a bacterium's natural habitat, striking phenomena-such as the volcano effect or banding-have been predicted or expected to emerge from chemotactic models. However, in practice, from limited bacterial trajectory data it is difficult to distinguish targeted searches from an untargeted search strategy for food sources. Here we use a theoretical model to identify statistical signatures of a targeted search toward point food sources, such as prey. Our model is constructed on the basis that bacteria use temporal comparisons to bias their random walk, exhibit finite memory and are subject to random (Brownian) motion as well as signaling noise. The advantage with using a stochastic model-based approach is that a stochastic model may be parametrized from individual stochastic bacterial trajectories but may then be used to generate a very large number of simulated trajectories to explore average behaviors obtained from stochastic search strategies. For example, our model predicts that a bacterium's diffusion coefficient increases as it approaches the point source and that, in the presence of multiple sources, bacteria may take substantially longer to locate their first source giving the impression of an untargeted search strategy.

Applications of Microfluidics in Quantitative Biology

Quantitative biology is dedicated to taking advantage of quantitative reasoning and advanced engineering technologies to make biology more predictable. Microfluidics, as an emerging technique, provides new approaches to precisely control fluidic conditions on small scales and collect data in high-throughput and quantitative manners. In this review, the authors present the relevant applications of microfluidics to quantitative biology based on two major categories (channel-based microfluidics and droplet-based microfluidics), and their typical features. We also envision some other microfluidic techniques that may not be employed in quantitative biology right now, but have great potential in the near future.

Drift and Behavior of E. coli Cells

Chemotaxis of the bacterium Escherichia coli is well understood in shallow chemical gradients, but its swimming behavior remains difficult to interpret in steep gradients. By focusing on single-cell trajectories from simulations, we investigated the dependence of the chemotactic drift velocity on attractant concentration in an exponential gradient. Whereas maxima of the average drift velocity can be interpreted within analytical linear-response theory of chemotaxis in shallow gradients, limits in drift due to steep gradients and finite number of receptor-methylation sites for adaptation go beyond perturbation theory. For instance, we found a surprising pinning of the cells to the concentration in the gradient at which cells run out of methylation sites. To validate the positions of maximal drift, we recorded single-cell trajectories in carefully designed chemical gradients using microfluidics.

Predicting Escherichia coli's chemotactic drift under exponential gradient

Bacterial species are known to show chemotaxis, i.e., the directed motions in the presence of certain chemicals, whereas the motion is random in the absence of those chemicals. The bacteria modulate their run time to induce chemotactic drift towards the attractant chemicals and away from the repellent chemicals. However, the existing theoretical knowledge does not exhibit a proper match with experimental validation, and hence there is a need for developing alternate models and validating experimentally. In this paper a more robust theoretical model is proposed to investigate chemotactic drift of peritrichous Escherichia coli under an exponential nutrient gradient. An exponential gradient is used to understand the steady state behavior of drift because of the logarithmic functionality of the chemosensory receptors. Our theoretical estimations are validated through the experimentation and simulation results. Thus, the developed model successfully delineates the run time, run trajectory, and drift velocity as measured from the experiments.

Propulsion and Chemotaxis in Bacteria-Driven Microswimmers

Despite the large body of experimental work recently on biohybrid microsystems, few studies have focused on theoretical modeling of such systems, which is essential to understand their underlying functioning mechanisms and hence design them optimally for a given application task. Therefore, this study focuses on developing a mathematical model to describe the 3D motion and chemotaxis of a type of widely studied biohybrid microswimmer, where spherical microbeads are driven by multiple attached bacteria. The model is developed based on the biophysical observations of the experimental system and is validated by comparing the model simulation with experimental 3D swimming trajectories and other motility characteristics, including mean squared displacement, speed, diffusivity, and turn angle. The chemotaxis modeling results of the microswimmers also agree well with the experiments, where a collective chemotactic behavior among multiple bacteria is observed. The simulation result implies that such collective chemotaxis behavior is due to a synchronized signaling pathway across the bacteria attached to the same microswimmer. Furthermore, the dependencies of the motility and chemotaxis of the microswimmers on certain system parameters, such as the chemoattractant concentration gradient, swimmer body size, and number of attached bacteria, toward an optimized design of such biohybrid system are studied. The optimized microswimmers would be used in targeted cargo, e.g., drug, imaging agent, gene, and RNA, transport and delivery inside the stagnant or low-velocity fluids of the human body as one of their potential biomedical applications.

Multi-fractal characterization of bacterial swimming dynamics: a case study on real and simulated Serratia marcescens

To add to the current state of knowledge about bacterial swimming dynamics, in this paper, we study the fractal swimming dynamics of populations of Serratia marcescens bacteria both in vitro and in silico, while accounting for realistic conditions like volume exclusion, chemical interactions, obstacles and distribution of chemoattractant in the environment. While previous research has shown that bacterial motion is non-ergodic, we demonstrate that, besides the non-ergodicity, the bacterial swimming dynamics is multi-fractal in nature. Finally, we demonstrate that the multi-fractal characteristic of bacterial dynamics is strongly affected by bacterial density and chemoattractant concentration.

Noise-Induced Increase of Sensitivity in Bacterial Chemotaxis

Flagellated bacteria, like Escherichia coli, can swim toward beneficial environments by modulating the rotational direction of their flagellar motors through a chemotaxis signal transduction network. The noise of this network, the random fluctuation of the intracellular concentration of the signal protein CheY-P with time, has been identified in studies of single cell behavioral variability, and found to be important in coordination of multiple motors in a bacterium and in enhancement of bacterial drift velocity in chemical gradients. Here, by comparing the behavioral difference between motors of wild-type E. coli and mutants without signal noise, we measured the magnitude of this noise in wild-type cells, and found that the noise increases the sensitivity of the bacterial chemotaxis network downstream at the level of the flagellar motor. This provided a simple mechanism for the noise-induced enhancement of chemotactic drift, which we confirmed by simulating the E. coli chemotactic motion in various spatial profiles of chemo-attractant concentration.

Growth-dependent behavioral difference in bacterial chemotaxis

Cells can adjust to their growth environments and regulate their behavior accordingly. To study how cells accomplish this growth-dependent adjustment from the molecular to the behavioral level, we used bacterial chemotaxis as a model system to explore the behavioral difference for bacteria grown in nutrient-rich and nutrient-poor media. We found that bacteria grown in a nutrient-poor medium exhibit faster chemotaxis adaptation, and this enables them to respond more rapidly to a changing environment and increases their ability to localize to a nutrient concentration peak. We identified the molecular mechanisms behind this behavioral difference through coarse-grained modeling, and demonstrated its physiological consequences by simulating bacterial chemotactic motion in spatiotemporally varying environments and in a static environment with a nutrient concentration peak.

Barrier Crossing in Escherichia coli Chemotaxis

We study cell navigation in spatiotemporally complex environments by developing a microfluidic racetrack device that creates a traveling wave with multiple peaks and a tunable wave speed. We find that while the population-averaged chemotaxis drift speed increases with wave speed for low wave speed, it decreases sharply for high wave speed. This reversed dependence of population-averaged chemotaxis drift speed on wave speed is caused by a "barrier-crossing" phenomenon, where a cell hops backwards from one peak attractant location to the peak behind by crossing an unfavorable (barrier) region with low attractant concentrations. By using a coarse-grained model of chemotaxis, we map bacterial motility in an attractant field to the random motion of an overdamped particle in an effective potential. The observed barrier-crossing phenomenon of living cells and its dependence on the spatiotemporal profile of attractant concentration are explained quantitatively by Kramers reaction rate theory.

Physical limits on bacterial navigation in dynamic environments

Many chemotactic bacteria inhabit environments in which chemicals appear as localized pulses and evolve by processes such as diffusion and mixing. We show that, in such environments, physical limits on the accuracy of temporal gradient sensing govern when and where bacteria can accurately measure the cues they use to navigate. Chemical pulses are surrounded by a predictable dynamic region, outside which bacterial cells cannot resolve gradients above noise. The outer boundary of this region initially expands in proportion to the square root of time before rapidly contracting. Our analysis also reveals how chemokinesis-the increase in swimming speed many bacteria exhibit when absolute chemical concentration exceeds a threshold-may serve to enhance chemotactic accuracy and sensitivity when the chemical landscape is dynamic. More generally, our framework provides a rigorous method for partitioning bacteria into populations that are 'near' and 'far' from chemical hotspots in complex, rapidly evolving environments such as those that dominate aquatic ecosystems.

Regional Homogeneity: A Multimodal, Multiscale Neuroimaging Marker of the Human Connectome

Much effort has been made to understand the organizational principles of human brain function using functional magnetic resonance imaging (fMRI) methods, among which resting-state fMRI (rfMRI) is an increasingly recognized technique for measuring the intrinsic dynamics of the human brain. Functional connectivity (FC) with rfMRI is the most widely used method to describe remote or long-distance relationships in studies of cerebral cortex parcellation, interindividual variability, and brain disorders. In contrast, local or short-distance functional interactions, especially at a scale of millimeters, have rarely been investigated or systematically reviewed like remote FC, although some local FC algorithms have been developed and applied to the discovery of brain-based changes under neuropsychiatric conditions. To fill this gap between remote and local FC studies, this review will (1) briefly survey the history of studies on organizational principles of human brain function; (2) propose local functional homogeneity as a network centrality to characterize multimodal local features of the brain connectome; (3) render a neurobiological perspective on local functional homogeneity by linking its temporal, spatial, and individual variability to information processing, anatomical morphology, and brain development; and (4) discuss its role in performing connectome-wide association studies and identify relevant challenges, and recommend its use in future brain connectomics studies.

Natural search algorithms as a bridge between organisms, evolution, and ecology

The ability to navigate is a hallmark of living systems, from single cells to higher animals. Searching for targets, such as food or mates in particular, is one of the fundamental navigational tasks many organisms must execute to survive and reproduce. Here, we argue that a recent surge of studies of the proximate mechanisms that underlie search behavior offers a new opportunity to integrate the biophysics and neuroscience of sensory systems with ecological and evolutionary processes, closing a feedback loop that promises exciting new avenues of scientific exploration at the frontier of systems biology.

Speed-dependent chemotactic precision in marine bacteria

Chemotaxis underpins important ecological processes in marine bacteria, from the association with primary producers to the colonization of particles and hosts. Marine bacteria often swim with a single flagellum at high speeds, alternating "runs" with either 180° reversals or ∼90° "flicks," the latter resulting from a buckling instability of the flagellum. These adaptations diverge from Escherichia coli's classic run-and-tumble motility, yet how they relate to the strong and rapid chemotaxis characteristic of marine bacteria has remained unknown. We investigated the relationship between swimming speed, run-reverse-flick motility, and high-performance chemotaxis by tracking thousands of Vibrio alginolyticus cells in microfluidic gradients. At odds with current chemotaxis models, we found that chemotactic precision-the strength of accumulation of cells at the peak of a gradient-is swimming-speed dependent in V. alginolyticus Faster cells accumulate twofold more tightly by chemotaxis compared with slower cells, attaining an advantage in the exploitation of a resource additional to that of faster gradient climbing. Trajectory analysis and an agent-based mathematical model revealed that this unexpected advantage originates from a speed dependence of reorientation frequency and flicking, which were higher for faster cells, and was compounded by chemokinesis, an increase in speed with resource concentration. The absence of any one of these adaptations led to a 65-70% reduction in the population-level resource exposure. These findings indicate that, contrary to what occurs in E. coli, swimming speed can be a fundamental determinant of the gradient-seeking capabilities of marine bacteria, and suggest a new model of bacterial chemotaxis.

Singly Flagellated Pseudomonas aeruginosa Chemotaxes Efficiently by Unbiased Motor Regulation

Pseudomonas aeruginosa is an opportunistic human pathogen that has long been known to chemotax. More recently, it has been established that chemotaxis is an important factor in the ability of P. aeruginosa to make biofilms. Genes that allow P. aeruginosa to chemotax are homologous with genes in the paradigmatic model organism for chemotaxis, Escherichia coli. However, P. aeruginosa is singly flagellated and E. coli has multiple flagella. Therefore, the regulation of counterclockwise/clockwise flagellar motor bias that allows E. coli to efficiently chemotax by runs and tumbles would lead to inefficient chemotaxis by P. aeruginosa, as half of a randomly oriented population would respond to a chemoattractant gradient in the wrong sense. How P. aeruginosa regulates flagellar rotation to achieve chemotaxis is not known. Here, we analyze the swimming trajectories of single cells in microfluidic channels and the rotations of cells tethered by their flagella to the surface of a variable-environment flow cell. We show that P. aeruginosa chemotaxes by symmetrically increasing the durations of both counterclockwise and clockwise flagellar rotations when swimming up the chemoattractant gradient and symmetrically decreasing rotation durations when swimming down the chemoattractant gradient. Unlike the case for E. coli, the counterclockwise/clockwise bias stays constant for P. aeruginosa. We describe P. aeruginosa’s chemotaxis using an analytical model for symmetric motor regulation. We use this model to do simulations that show that, given P. aeruginosa’s physiological constraints on motility, its distinct, symmetric regulation of motor switching optimizes chemotaxis.

Chemotaxis has long been known to strongly affect biofilm formation by the opportunistic human pathogen P. aeruginosa, whose essential chemotaxis genes have homologues in E. coli, which achieves chemotaxis by biasing the relative probability of counterclockwise and clockwise flagellar rotation. However, the physiological difference between multiflagellated E. coli and singly flagellated P. aeruginosa implies that biased motor regulation should prevent P. aeruginosa populations from chemotaxing efficiently. Here, we used experiments, analytical modeling, and simulations to demonstrate that P. aeruginosa uses unbiased, symmetric regulation of the flagellar motor to maximize its chemotaxis efficiency. This mode of chemotaxis was not previously known and demonstrates a new variant of a paradigmatic signaling system in an important human pathogen.

Derivation of the bacterial run-and-tumble kinetic equation from a model with biochemical pathway

Kinetic-transport equations are, by now, standard models to describe the dynamics of populations of bacteria moving by run-and-tumble. Experimental observations show that bacteria increase their run duration when encountering an increasing gradient of chemotactic molecules. This led to a first class of models which heuristically include tumbling frequencies depending on the path-wise gradient of chemotactic signal. More recently, the biochemical pathways regulating the flagellar motors were uncovered. This knowledge gave rise to a second class of kinetic-transport equations, that takes into account an intra-cellular molecular content and which relates the tumbling frequency to this information. It turns out that the tumbling frequency depends on the chemotactic signal, and not on its gradient. For these two classes of models, macroscopic equations of Keller-Segel type, have been derived using diffusion or hyperbolic rescaling. We complete this program by showing how the first class of equations can be derived from the second class with molecular content after appropriate rescaling. The main difficulty is to explain why the path-wise gradient of chemotactic signal can arise in this asymptotic process. Randomness of receptor methylation events can be included, and our approach can be used to compute the tumbling frequency in presence of such a noise.

Inferring Models of Bacterial Dynamics toward Point Sources

Experiments have shown that bacteria can be sensitive to small variations in chemoattractant (CA) concentrations. Motivated by these findings, our focus here is on a regime rarely studied in experiments: bacteria tracking point CA sources (such as food patches or even prey). In tracking point sources, the CA detected by bacteria may show very large spatiotemporal fluctuations which vary with distance from the source. We present a general statistical model to describe how bacteria locate point sources of food on the basis of stochastic event detection, rather than CA gradient information. We show how all model parameters can be directly inferred from single cell tracking data even in the limit of high detection noise. Once parameterized, our model recapitulates bacterial behavior around point sources such as the “volcano effect”. In addition, while the search by bacteria for point sources such as prey may appear random, our model identifies key statistical signatures of a targeted search for a point source given any arbitrary source configuration.

Predicting chemical environments of bacteria from receptor signaling

Sensory systems have evolved to respond to input stimuli of certain statistical properties, and to reliably transmit this information through biochemical pathways. Hence, for an experimentally well-characterized sensory system, one ought to be able to extract valuable information about the statistics of the stimuli. Based on dose-response curves from in vivo fluorescence resonance energy transfer (FRET) experiments of the bacterial chemotaxis sensory system, we predict the chemical gradients chemotactic Escherichia coli cells typically encounter in their natural environment. To predict average gradients cells experience, we revaluate the phenomenological Weber's law and its generalizations to the Weber-Fechner law and fold-change detection. To obtain full distributions of gradients we use information theory and simulations, considering limitations of information transmission from both cell-external and internal noise. We identify broad distributions of exponential gradients, which lead to log-normal stimuli and maximal drift velocity. Our results thus provide a first step towards deciphering the chemical nature of complex, experimentally inaccessible cellular microenvironments, such as the human intestine.

Outside the laboratory, bacteria live in complex microenvironments characterized by competition for space and available nutrients. Although often inaccessible by experiments, understanding the spatio-temporal dynamics of bacterial microenvironments is biomedically important. For instance, the chemical environment that symbiotic Escherichia coli encounter in the human gut relates to health of the gastrointestinal tract, gut metabolism, immune response, and tissue homeostasis. Other complex microenvironments include soil and biofilms. Assuming that bacterial sensory systems have evolved to optimally sense typical gradients, we treat signaling data, the signaling pathway with its architecture and reaction rates, and computer simulations of swimming bacteria in different gradients as “prior knowledge” to “reverse engineer” E. coli's habitat. Our identified gradients are exponentially shaped with wide-ranging rate values. These microenvironments most likely stem from local fluctuating nutrient sources and degradation by competing species, in which bacteria have evolved to swim with optimal performance.

Adaptability of non-genetic diversity in bacterial chemotaxis

Bacterial chemotaxis systems are as diverse as the environments that bacteria inhabit, but how much environmental variation can cells tolerate with a single system? Diversification of a single chemotaxis system could serve as an alternative, or even evolutionary stepping-stone, to switching between multiple systems. We hypothesized that mutations in gene regulation could lead to heritable control of chemotactic diversity. By simulating foraging and colonization of E. coli using a single-cell chemotaxis model, we found that different environments selected for different behaviors. The resulting trade-offs show that populations facing diverse environments would ideally diversify behaviors when time for navigation is limited. We show that advantageous diversity can arise from changes in the distribution of protein levels among individuals, which could occur through mutations in gene regulation. We propose experiments to test our prediction that chemotactic diversity in a clonal population could be a selectable trait that enables adaptation to environmental variability.

DOI:http://dx.doi.org/10.7554/eLife.03526.001

Bacterial colonies are generally made up of genetically identical cells. Despite this, a closer look at the members of a bacterial colony shows that these cells can have very different behaviors. For example, some cells may grow more quickly than others, or be more resistant to antibiotics. The mechanisms driving this diversity are only beginning to be identified and understood.

Escherichia coli bacteria can move towards, or away from, certain chemicals in their surrounding environment to help them navigate toward favorable conditions. This behavior is known as chemotaxis. The signals from all of these chemicals are processed in E. coli by just one set of proteins, which control the different behaviors that are needed for the bacteria to follow them. Different numbers of these proteins are found in different—but genetically identical—bacteria, and the number of proteins is linked to how the bacteria perform these behaviors.

It has been suggested that diversity can be beneficial to the overall bacterial population, as it helps the population survive environmental changes. This suggests that the level of diversity in the population should adapt to the level of diversity in the environment. However, it remains unknown how this adaptation occurs.

Frankel et al. developed and combined several models and simulations to investigate whether differences in chemotaxis protein production help an E. coli colony to survive. The models show that in different environments, it can be beneficial for the population as a whole if different cells have different responses to the chemicals present. For example, if a lot of a useful chemical is present, bacteria are more likely to survive by heading straight to the source. If not much chemical is detected, the bacteria may need to move in a more exploratory manner.

Frankel et al. find that different amounts of chemotaxis proteins produce these different behaviors. To survive in a changing environment, it is therefore best for the E. coli colony to contain cells that have different amounts of these proteins. Frankel et al. propose that the variability of chemotaxis protein levels between genetically identical cells can change through mutations in the genes that control how many of the proteins are produced, and predict that such mutations allow populations to adapt to environmental changes.

The environments simulated in the model were much simpler than would be found in the real world, and Frankel et al. describe experiments that are now being performed to confirm and expand on their results. The model could be used in the future to shed light on the behavior of other cells that are genetically identical but exhibit diverse behaviors, from other bacterial species to more complex cancer cells.

DOI:http://dx.doi.org/10.7554/eLife.03526.002

Behaviors and strategies of bacterial navigation in chemical and nonchemical gradients

Navigation of cells to the optimal environmental condition is critical for their survival and growth. Escherichia coli cells, for example, can detect various chemicals and move up or down those chemical gradients (i.e., chemotaxis). Using the same signaling machinery, they can also sense other external factors such as pH and temperature and navigate from both sides toward some intermediate levels of those stimuli. This mode of precision sensing is more sophisticated than the (unidirectional) chemotaxis strategy and requires distinctive molecular mechanisms to encode and track the preferred external conditions. To systematically study these different bacterial taxis behaviors, we develop a continuum model that incorporates microscopic signaling events in single cells into macroscopic population dynamics. A simple theoretical result is obtained for the steady state cell distribution in general. In particular, we find the cell distribution is controlled by the intracellular sensory dynamics as well as the dependence of the cells' speed on external factors. The model is verified by available experimental data in various taxis behaviors (including bacterial chemotaxis, pH taxis, and thermotaxis), and it also leads to predictions that can be tested by future experiments. Our analysis help reveal the key conditions/mechanisms for bacterial precision-sensing behaviors and directly connects the cellular taxis performances with the underlying molecular parameters. It provides a unified framework to study bacterial navigation in complex environments with chemical and non-chemical stimuli.

Bacteria, such as E. coli, live in a complex environment with varying chemical and/or non-chemical stimuli. They constantly seek for and migrate to optimal environmental conditions. A well-known example is E. coli chemotaxis which direct cell movements up or down chemical gradients. Using the same machinery, E. coli can also respond to non-chemical factors (e.g., pH and temperature) and navigate toward certain intermediate, optimal levels of those stimuli. Such taxis behaviors are more sophisticated and require distinctive sensing mechanisms. In this paper, we develop a unified model for different bacterial taxis strategies. This multiscale model incorporates intracellular signaling pathways into population dynamics and leads to a simple theoretical result regarding the steady-state population distribution. Our model can be applied to reveal the key mechanisms for different taxis behaviors and quantitatively account for various experimental data. New predictions can be made within this new model framework to direct future experiments.

Fold-change detection in a whole-pathway model of Escherichia coli chemotaxis

There has been recent interest in sensory systems that are able to display a response which is proportional to a fold change in stimulus concentration, a feature referred to as fold-change detection (FCD). Here, we demonstrate FCD in a recent whole-pathway mathematical model of Escherichia coli chemotaxis. FCD is shown to hold for each protein in the signalling cascade and to be robust to kinetic rate and protein concentration variation. Using a sensitivity analysis, we find that only variations in the number of receptors within a signalling team lead to the model not exhibiting FCD. We also discuss the ability of a cell with multiple receptor types to display FCD and explain how a particular receptor configuration may be used to elucidate the two experimentally determined regimes of FCD behaviour. All findings are discussed in respect of the experimental literature.

Analytical modeling and experimental characterization of chemotaxis in Serratia marcescens

This paper presents a modeling and experimental framework to characterize the chemotaxis of Serratia marcescens (S. marcescens) relying on two-dimensional and three-dimensional tracking of individual bacteria. Previous studies mainly characterized bacterial chemotaxis based on population density analysis. Instead, this study focuses on single-cell tracking and measuring the chemotactic drift velocity V(C) from the biased tumble rate of individual bacteria on exposure to a concentration gradient of l-aspartate. The chemotactic response of S. marcescens is quantified over a range of concentration gradients (10^{-3} to 5 mM/mm) and average concentrations (0.5 × 10(-3) to 2.5 mM). Through the analysis of a large number of bacterial swimming trajectories, the tumble rate is found to have a significant bias with respect to the swimming direction. We also verify the relative gradient sensing mechanism in the chemotaxis of S. marcescens by measuring the change of V(C) with the average concentration and the gradient. The applied full pathway model with fitted parameters matches the experimental data. Finally, we show that our measurements based on individual bacteria lead to the determination of the motility coefficient μ (7.25 × 10(-6) cm(2)/s) of a population. The experimental characterization and simulation results for the chemotaxis of this bacterial species contribute towards using S. marcescens in chemically controlled biohybrid systems.