Noninvasive inference of the molecular chemotactic response using bacterial trajectories

Kinetic modeling of the chemotactic process in run-and-tumble bacteria

The chemotactic process of run-and-tumble bacteria results from modulating the tumbling rate in response to changes in chemoattractant gradients felt by the bacteria. The response has a characteristic memory time and is subject to important fluctuations. These ingredients are considered in a kinetic description of chemotaxis, allowing the computation of the stationary mobility and the relaxation times needed to reach the steady state. For large memory times, these relaxation times become large, implying that finite-time measurements give rise to nonmonotonic currents as a function of the imposed chemoattractant gradient, contrary to the stationary regime where the response is monotonic. The case of an inhomogeneous signal is analyzed. Contrary to the usual Keller-Segel model, the response is nonlocal, and the bacterial profile is smoothed with a characteristic length that grows with the memory time. Finally, the case of traveling signals is considered, where appreciable differences appear compared to memoryless chemotactic descriptions.

An internal expectation guides Drosophila egg-laying decisions

To better understand how animals make ethologically relevant decisions, we studied egg-laying substrate choice in Drosophila. We found that flies dynamically increase or decrease their egg-laying rates while exploring substrates so as to target eggs to the best, recently visited option. Visiting the best option typically yielded inhibition of egg laying on other substrates for many minutes. Our data support a model in which flies compare the current substrate's value with an internally constructed expectation on the value of available options to regulate the likelihood of laying an egg. We show that dopamine neuron activity is critical for learning and/or expressing this expectation, similar to its role in certain tasks in vertebrates. Integrating sensory experiences over minutes to generate an estimate of the quality of available options allows flies to use a dynamic reference point for judging the current substrate and might be a general way in which decisions are made.

Steady-state running rate sets the speed and accuracy of accumulation of swimming bacteria

We study the chemotaxis of a population of genetically identical swimming bacteria undergoing run and tumble dynamics driven by stochastic switching between clockwise and counterclockwise rotation of the flagellar rotary system, where the steady-state rate of the switching changes in different environments. Understanding chemotaxis quantitatively requires that one links the measured steady-state switching rates of the rotary system, as well as the directional changes of individual swimming bacteria in a gradient of chemoattractant/repellant, to the efficiency of a population of bacteria in moving up/down the gradient. Here we achieve this by using a probabilistic model, parametrized with our experimental data, and show that the response of a population to the gradient is complex. We find the changes to the steady-state switching rate in the absence of gradients affect the average speed of the swimming bacterial population response as well as the width of the distribution. Both must be taken into account when optimizing the overall response of the population in complex environments.

High-throughput automated methods for classical and operant conditioning of Drosophila larvae

Learning which stimuli (classical conditioning) or which actions (operant conditioning) predict rewards or punishments can improve chances of survival. However, the circuit mechanisms that underlie distinct types of associative learning are still not fully understood. Automated, high-throughput paradigms for studying different types of associative learning, combined with manipulation of specific neurons in freely behaving animals, can help advance this field. The Drosophila melanogaster larva is a tractable model system for studying the circuit basis of behaviour, but many forms of associative learning have not yet been demonstrated in this animal. Here, we developed a high-throughput (i.e. multi-larva) training system that combines real-time behaviour detection of freely moving larvae with targeted opto- and thermogenetic stimulation of tracked animals. Both stimuli are controlled in either open- or closed-loop, and delivered with high temporal and spatial precision. Using this tracker, we show for the first time that Drosophila larvae can perform classical conditioning with no overlap between sensory stimuli (i.e. trace conditioning). We also demonstrate that larvae are capable of operant conditioning by inducing a bend direction preference through optogenetic activation of reward-encoding serotonergic neurons. Our results extend the known associative learning capacities of Drosophila larvae. Our automated training rig will facilitate the study of many different forms of associative learning and the identification of the neural circuits that underpin them.

Escherichia coli chemotaxis is information limited

Organisms acquire and use information from their environment to guide their behaviour. However, it is unclear whether this information quantitatively limits their behavioural performance. Here, we relate information to the ability of Escherichia coli to navigate up chemical gradients, the behaviour known as chemotaxis. First, we derive a theoretical limit on the speed with which cells climb gradients, given the rate at which they acquire information. Next, we measure cells' gradient-climbing speeds and the rate of information acquisition by their chemotaxis signaling pathway. We find that E. coli make behavioural decisions with much less than the one bit required to determine whether they are swimming up-gradient. Some of this information is irrelevant to gradient climbing, and some is lost in communication to behaviour. Despite these limitations, E. coli climb gradients at speeds within a factor of two of the theoretical bound. Thus, information can limit the performance of an organism, and sensory-motor pathways may have evolved to efficiently use information acquired from the environment.

Non-Genetic Diversity in Chemosensing and Chemotactic Behavior

Non-genetic phenotypic diversity plays a significant role in the chemotactic behavior of bacteria, influencing how populations sense and respond to chemical stimuli. First, we review the molecular mechanisms that generate phenotypic diversity in bacterial chemotaxis. Next, we discuss the functional consequences of phenotypic diversity for the chemosensing and chemotactic performance of single cells and populations. Finally, we discuss mechanisms that modulate the amount of phenotypic diversity in chemosensory parameters in response to changes in the environment.

Ligand sensing enhances bacterial flagellar motor output via stator recruitment

It is well known that flagellated bacteria, such as Escherichia coli, sense chemicals in their environment by a chemoreceptor and relay the signals via a well-characterized signaling pathway to the flagellar motor. It is widely accepted that the signals change the rotation bias of the motor without influencing the motor speed. Here, we present results to the contrary and show that the bacteria is also capable of modulating motor speed on merely sensing a ligand. Step changes in concentration of non-metabolizable ligand cause temporary recruitment of stator units leading to a momentary increase in motor speeds. For metabolizable ligand, the combined effect of sensing and metabolism leads to higher motor speeds for longer durations. Experiments performed with mutant strains delineate the role of metabolism and sensing in the modulation of motor speed and show how speed changes along with changes in bias can significantly enhance response to changes in its environment.

Single-trajectory characterization of active swimmers in a flow

We develop a maximum likelihood method to infer relevant physical properties of elongated active particles. Using individual trajectories of advected swimmers as input, we are able to accurately determine their rotational diffusion coefficients and an effective measure of their aspect ratio, also providing reliable estimators for the uncertainties of such quantities. We validate our theoretical construction using numerically generated active trajectories upon no flow, simple shear, and Poiseuille flow, with excellent results. Being designed to rely on single-particle data, our method eases applications in experimental conditions where swimmers exhibit a strong morphological diversity. We briefly discuss some of such ongoing experimental applications, specifically, in the characterization of swimming E. coli in a flow.

Swimming Escherichia coli Cells Explore the Environment by Lévy Walk

Escherichia coli cells swim in aqueous environment in a random walk of alternating runs and tumbles. The diffusion characteristics of this random walk remains unclear. In this study, by tracking the swimming of wild-type cells in a three-dimensional (3D) homogeneous environment, we found that their trajectories are superdiffusive, consistent with Lévy walk behavior. For comparison, we tracked the swimming of mutant cells that lack the chemotaxis signaling noise (the steady-state fluctuation of the concentration of the chemotaxis response regulator CheY-P) and found that their trajectories are normal diffusive. Therefore, wild-type E. coli cells explore the environment by Lévy walk, which originates from the chemotaxis signaling noise. This Lévy walk pattern enhances their efficiency in environmental exploration.IMPORTANCE E. coli cells explore the environment in a random walk of alternating runs and tumbles. By tracking the 3D trajectories of E. coli cells in an aqueous environment, we found that their trajectories are superdiffusive, with a power-law shape for the distribution of run lengths, which is characteristics of Lévy walk. We further show that this Lévy walk behavior is due to the random fluctuation of the output level of the bacterial chemotaxis pathway, and it enhances the efficiency of the bacteria in exploring the environment.

Chemotaxis strategies of bacteria with multiple run modes

Bacterial chemotaxis-a fundamental example of directional navigation in the living world-is key to many biological processes, including the spreading of bacterial infections. Many bacterial species were recently reported to exhibit several distinct swimming modes-the flagella may, for example, push the cell body or wrap around it. How do the different run modes shape the chemotaxis strategy of a multimode swimmer? Here, we investigate chemotactic motion of the soil bacterium Pseudomonas putida as a model organism. By simultaneously tracking the position of the cell body and the configuration of its flagella, we demonstrate that individual run modes show different chemotactic responses in nutrition gradients and, thus, constitute distinct behavioral states. On the basis of an active particle model, we demonstrate that switching between multiple run states that differ in their speed and responsiveness provides the basis for robust and efficient chemotaxis in complex natural habitats.

Statistical Tests for Force Inference in Heterogeneous Environments

We devise a method to detect and estimate forces in a heterogeneous environment based on experimentally recorded stochastic trajectories. In particular, we focus on systems modeled by the heterogeneous overdamped Langevin equation. Here, the observed drift includes a "spurious” force term when the diffusivity varies in space. We show how Bayesian inference can be leveraged to reliably infer forces by taking into account such spurious forces of unknown amplitude as well as experimental sources of error. The method is based on marginalizing the force posterior over all possible spurious force contributions. The approach is combined with a Bayes factor statistical test for the presence of forces. The performance of our method is investigated analytically, numerically and tested on experimental data sets. The main results are obtained in a closed form allowing for direct exploration of their properties and fast computation. The method is incorporated into TRamWAy, an open-source software platform for automated analysis of biomolecule trajectories.

Chemotaxis as a navigation strategy to boost range expansion

Bacterial chemotaxis, the directed movement of cells along gradients of chemoattractants, is among the best-characterized subjects in molecular biology^1-10, but much less is known about its physiological roles¹¹. It is commonly seen as a starvation response when nutrients run out, or as an escape response from harmful situations^12-16. Here we identify an alternative role of chemotaxis by systematically examining the spatiotemporal dynamics of Escherichia coli in soft agar^12,17,18. Chemotaxis in nutrient-replete conditions promotes the expansion of bacterial populations into unoccupied territories well before nutrients run out in the current environment. Low levels of chemoattractants act as aroma-like cues in this process, establishing the direction and enhancing the speed of population movement along the self-generated attractant gradients. This process of navigated range expansion spreads faster and yields larger population gains than unguided expansion following the canonical Fisher-Kolmogorov dynamics^19,20 and is therefore a general strategy to promote population growth in spatially extended, nutrient-replete environments.

Enhanced propagation of motile bacteria on surfaces due to forward scattering

How motile bacteria move near a surface is a problem of fundamental biophysical interest and is key to the emergence of several phenomena of biological, ecological and medical relevance, including biofilm formation. Solid boundaries can strongly influence a cell's propulsion mechanism, thus leading many flagellated bacteria to describe long circular trajectories stably entrapped by the surface. Experimental studies on near-surface bacterial motility have, however, neglected the fact that real environments have typical microstructures varying on the scale of the cells' motion. Here, we show that micro-obstacles influence the propagation of peritrichously flagellated bacteria on a flat surface in a non-monotonic way. Instead of hindering it, an optimal, relatively low obstacle density can significantly enhance cells' propagation on surfaces due to individual forward-scattering events. This finding provides insight on the emerging dynamics of chiral active matter in complex environments and inspires possible routes to control microbial ecology in natural habitats.

Bacterial chemotaxis in a microfluidic T-maze reveals strong phenotypic heterogeneity in chemotactic sensitivity

Many microorganisms have evolved chemotactic strategies to exploit the microscale heterogeneity that frequently characterizes microbial habitats. Chemotaxis has been primarily studied as an average characteristic of a population, with little regard for variability among individuals. Here, we adopt a classic tool from animal ecology - the T-maze - and implement it at the microscale by using microfluidics to expose bacteria to a sequence of decisions, each consisting of migration up or down a chemical gradient. Single-cell observations of clonal Escherichia coli in the maze, coupled with a mathematical model, reveal that strong heterogeneity in the chemotactic sensitivity coefficient exists even within clonal populations of bacteria. A comparison of different potential sources of heterogeneity reveals that heterogeneity in the T-maze originates primarily from the chemotactic sensitivity coefficient, arising from a distribution of pathway gains. This heterogeneity may have a functional role, for example in the context of migratory bet-hedging strategies.

Connecting single-cell properties to collective behavior in multiple wild isolates of the Enterobacter cloacae complex

Some strains of motile bacteria self-organize to form spatial patterns of high and low cell density over length scales that can be observed by eye. One such collective behavior is the formation in semisolid agar media of a high cell density swarm band. We isolated 7 wild strains of the Enterobacter cloacae complex capable of forming this band and found its propagation speed can vary 2.5 fold across strains. To connect such variability in collective motility to strain properties, each strain’s single-cell motility and exponential growth rates were measured. The band speed did not significantly correlate with any individual strain property; however, a multilinear analysis revealed that the band speed was set by a combination of the run speed and tumbling frequency. Comparison of variability in closely-related wild isolates has the potential to reveal how changes in single-cell properties influence the collective behavior of populations.

Multiple CheY Homologs Control Swimming Reversals and Transient Pauses in Azospirillum brasilense

Chemotaxis, together with motility, helps bacteria foraging in their habitat. Motile bacteria exhibit a variety of motility patterns, often controlled by chemotaxis, to promote dispersal. Motility in many bacteria is powered by a bidirectional flagellar motor. The flagellar motor has been known to briefly pause during rotation because of incomplete reversals or stator detachment. Transient pauses were previously observed in bacterial strains lacking CheY, and these events could not be explained by incomplete motor reversals or stator detachment. Here, we systematically analyzed swimming trajectories of various chemotaxis mutants of the monotrichous soil bacterium, Azospirillum brasilense. Like other polar flagellated bacterium, the main swimming pattern in A. brasilense is run and reverse. A. brasilense also uses run-pauses and putative run-reverse-flick-like swimming patterns, although these are rare events. A. brasilense mutant derivatives lacking the chemotaxis master histidine kinase, CheA4, or the central response regulator, CheY7, also showed transient pauses. Strikingly, the frequency of transient pauses increased dramatically in the absence of CheY4. Our findings collectively suggest that reversals and pauses are controlled through signaling by distinct CheY homologs, and thus are likely to be functionally important in the lifestyle of this soil organism.

Applying machine learning to the flagellar motor for biosensing

Escherichia coli detects and follows chemical gradients in its environment in a process known as chemotaxis. The performance of chemotaxis approaches fundamental biosensor speed and sensitivity limits, but there have been relatively few attempts to incorporate the response into a functional biosensor. Toward that end, we have developed software to process digital microscope images of a large number of tethered E. coli responding to different chemical perturbations. Upwards of fifty cells can be recorded in one experiment, allowing for rapid labeling of the chemotactic responses of multiple cells. After we collected hundreds of wild-type chemotactic E. coli motor responses to dilutions of aspartate and leucine, we trained a support vector classifier (SVC) to estimate the order of magnitude of aspartate concentration between 0M, 100nM, and 1μM with a single cell classification subset accuracy of 69%. We trained another SVC to differentiate between aspartate and leucine with a single cell classification subset accuracy of 83%. Using a majority-vote method on a bacterial population of size N, estimates have 95% confidence for N = 27 bacteria for concentration detection and N = 9 bacteria for chemical differentiation. These methods are a step towards adaptable chemotaxis-based biosensing.

Flagellar number governs bacterial spreading and transport efficiency

Peritrichous bacteria synchronize and bundle their flagella to actively swim, while disruption of the bundle leads to a slow motility phase with a weak propulsion. It is still not known whether the number of flagella represents an evolutionary adaptation toward optimizing bacterial navigation. We study the swimming dynamics of differentially flagellated Bacillus subtilis strains in a quasi-two-dimensional system. We find that decreasing the number of flagella N f reduces the average turning angle between two successive run phases and enhances the run time and the directional persistence of the run phase. As a result, having fewer flagella is beneficial for long-distance transport and fast spreading, while having a lot of flagella is advantageous for the processes that require a slower spreading, such as biofilm formation. We develop a two-state random walk model that incorporates spontaneous switchings between the states and yields exact analytical expressions for transport properties, in remarkable agreement with experiments. The results of numerical simulations based on our two-state model suggest that the efficiency of searching and exploring the environment is optimized at intermediate values of N f. The optimal choice of N f, for which the search time is minimized, decreases with increasing the size of the environment in which the bacteria swim.

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.

Stimulus sensing and signal processing in bacterial chemotaxis

Motile bacteria use chemotaxis to migrate towards environments that are favorable for growth and survival. The signaling pathway that mediates this behavior is largely conserved among prokaryotes, with Escherichia coli chemotaxis system being one of the simplest and the best studied. At the core of this pathway are the arrays of clustered chemoreceptors that detect, amplify and integrate various stimuli. Recent work provided deeper understanding of spatial organization and signal processing by these clusters and uncovered the variety of sensory mechanisms used to detect environmental stimuli. Moreover, studies of bacteria with different lifestyles have led to new insights into the diversity and evolutionary conservation of the chemotaxis pathway, as well as the physiological relevance of chemotactic behavior in different environments.

Phenotypic diversity and temporal variability in a bacterial signaling network revealed by single-cell FRET

We present in vivo single-cell FRET measurements in the Escherichia coli chemotaxis system that reveal pervasive signaling variability, both across cells in isogenic populations and within individual cells over time. We quantify cell-to-cell variability of adaptation, ligand response, as well as steady-state output level, and analyze the role of network design in shaping this diversity from gene expression noise. In the absence of changes in gene expression, we find that single cells demonstrate strong temporal fluctuations. We provide evidence that such signaling noise can arise from at least two sources: (i) stochastic activities of adaptation enzymes, and (ii) receptor-kinase dynamics in the absence of adaptation. We demonstrate that under certain conditions, (ii) can generate giant fluctuations that drive signaling activity of the entire cell into a stochastic two-state switching regime. Our findings underscore the importance of molecular noise, arising not only in gene expression but also in protein networks.

Many sophisticated computer programs use random number generators to help solve challenging problems. These problems range from achieving secure communication across the Internet to deciding how best to invest in the stock market. Much research in recent years has found that randomness is also widespread in living cells, where it is often called “noise”. For example, the activity of some genes is so unpredictable to the extent that it appears random. Yet, relatively little is known about how such gene-expression noise propagates up to change how the cell behaves. Many open questions also remain about how cells might exploit these or other fluctuations to achieve complex tasks, like people use random number generators.

Bacteria perform a number of complex tasks. Some bacteria will swim toward chemicals that suggest a potential reward, such as food. Yet they swim away from chemicals that could lead them to harm. This ability is called chemotaxis and it relies on a network of interacting enzymes and other proteins that coordinates a bacterium’s movements with the input from its senses.

Keegstra et al. set out to find sources of noise that might act as random number generators and help the bacterium E. coli to best perform chemotaxis. An improved version of a technique called in vivo Förster resonance energy transfer (or in vivo FRET for short) was used to give a detectable signal when two proteins involved in the chemotaxis network interacted inside a single bacterium. The experiments showed that this protein network amplifies gene-expression noise for some genes while lessening it for others. In addition, the interactions between proteins encoded by genes acted as an extra source of noise, even when gene-expression noise was eliminated.

Keegstra et al. found that the amount of signaling within the chemotaxis network, as measured by in vivo FRET, varied wildly over time. This revealed two sources of noise at the level of protein signaling. One was due to randomness in the activity of the enzymes involved in tuning the cell’s sensitivity to changes in its environment. The other was due to protein interactions within a large complex that acts as the cell’s sensor. Unexpectedly, this second source of noise under some conditions could be so strong that it flipped the output of the cell’s signaling network back and forth between just two states: “on” and “off”.

Together these findings uncover how signaling networks can not only amplify or lessen gene-expression noise, but can themselves become a source of random events. The new knowledge of how such random events interact with a complex trait in a living cell – namely chemotaxis – could aid future antimicrobial strategies, because many bacteria use chemotaxis to help them establish infections. More generally, the new insights about noise in protein networks could help engineers seeking to build synthetic biochemical networks or produce useful compounds in living cells.

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.

Emergent properties of bacterial chemotaxis pathway

The chemotaxis pathway of Escherichia coli is the most studied sensory system in prokaryotes. The highly conserved general architecture of this pathway consists of two modules which mediate signal transduction and adaptation. The signal transduction module detects and amplifies changes in environmental conditions and rapidly transmits these signals to control bacterial swimming behavior. The adaptation module gradually resets the activity and sensitivity of the first module after initial stimulation and thereby enables the temporal comparisons necessary for bacterial chemotaxis. Recent experimental and theoretical work has unraveled multiple quantitative features emerging from the interplay between these two modules. This has laid the groundwork for rationalization of these emerging properties in the context of the evolutionary optimization of the chemotactic behavior.

Ultrasensitivity and fluctuations in the Barkai-Leibler model of chemotaxis receptors in Escherichia coli

A stochastic version of the Barkai-Leibler model of chemotaxis receptors in Escherichia coli is studied here with the goal of elucidating the effects of intrinsic network noise in their conformational dynamics. The model was originally proposed to explain the robust and near-perfect adaptation of E. coli observed across a wide range of spatially uniform attractant/repellent (ligand) concentrations. In the model, a receptor is either active or inactive and can stochastically switch between the two states. The enzyme CheR methylates inactive receptors while CheB demethylates active receptors and the probability for a receptor to be active depends on its level of methylation and ligand occupation. In a simple version of the model with two methylation sites per receptor (M = 2), we show rigorously, under a quasi-steady state approximation, that the mean active fraction of receptors is an ultrasensitive function of [CheR]/[CheB] in the limit of saturating receptor concentration. Hence the model shows zero-order ultrasensitivity (ZOU), similar to the classical two-state model of covalent modification studied by Goldbeter and Koshland (GK). We also find that in the limits of extremely small and extremely large ligand concentrations, the system reduces to two different two-state GK modules. A quantitative measure of the spontaneous fluctuations in activity is provided by the variance σa2 in the active fraction, which is estimated mathematically under linear noise approximation (LNA). It is found that σa2 peaks near the ZOU transition. The variance is a non-monotonic, but weak function of ligand concentration and a decreasing function of receptor concentration. Gillespie simulations are also performed in models with M = 2, 3 and 4. For M = 2, simulations show excellent agreement with analytical results obtained under LNA. Numerical results for M = 3 and M = 4 are qualitatively similar to our mathematical results in M = 2; while all the models show ZOU in mean activity, the variance is found to be smaller for larger M. The magnitude of receptor noise deduced from available experimental data is consistent with our predictions. A simple analysis of the downstream signaling pathway shows that this noise is large enough to affect the motility of the organism, and may have a beneficial effect on it. The response of mean receptor activity to small time-dependent changes in the external ligand concentration is computed within linear response theory, and found to have a bilobe form, in agreement with earlier experimental observations.

Motility of Escherichia coli in a quasi-two-dimensional porous medium

Bacterial migration through confined spaces is critical for several phenomena, such as biofilm formation, bacterial transport in soils, and bacterial therapy against cancer. In the present work, E. coli (strain K12-MG1655 WT) motility was characterized by recording and analyzing individual bacterium trajectories in a simulated quasi-two-dimensional porous medium. The porous medium was simulated by enclosing, between slide and cover slip, a bacterial-culture sample mixed with uniform 2.98-μm-diameter spherical latex particles. The porosity of the medium was controlled by changing the latex particle concentration. By statistically analyzing several trajectory parameters (instantaneous velocity, turn angle, mean squared displacement, etc.), and contrasting with the results of a random-walk model developed ad hoc, we were able to quantify the effects that different obstacle concentrations have upon bacterial motility.

Inferring the Chemotactic Strategy of P. putida and E. coli Using Modified Kramers-Moyal Coefficients

Many bacteria perform a run-and-tumble random walk to explore their surrounding and to perform chemotaxis. In this article we present a novel method to infer the relevant parameters of bacterial motion from experimental trajectories including the tumbling events. We introduce a stochastic model for the orientation angle, where a shot-noise process initiates tumbles, and analytically calculate conditional moments, reminiscent of Kramers-Moyal coefficients. Matching them with the moments calculated from experimental trajectories of the bacteria E. coli and Pseudomonas putida, we are able to infer their respective tumble rates, the rotational diffusion constants, and the distributions of tumble angles in good agreement with results from conventional tumble recognizers. We also define a novel tumble recognizer, which explicitly quantifies the error in recognizing tumbles. In the presence of a chemical gradient we condition the moments on the bacterial direction of motion and thereby explore the chemotaxis strategy. For both bacteria we recover and quantify the classical chemotactic strategy, where the tumble rate is smallest along the chemical gradient. In addition, for E. coli we detect some cells, which bias their mean tumble angle towards smaller values. Our findings are supported by a scaling analysis of appropriate ratios of conditional moments, which are directly calculated from experimental data.

The movement strategies of bacteria have received increasing attention over the past decade, in particular with respect to the tracking of individual cells and the mathematical description of the resulting trajectories. Bacteria typically move in almost straight runs interrupted by sharp turning events (run-and-tumble). In order to characterize their motion on a single cell level, the tumble events in individual trajectories have to be identified. Traditionally, tumble recognition relies on threshold values that are applied to the swimming speed and the reorientation angle. They are chosen in an ad hoc fashion and introduce a certain degree of arbitrariness to the results of statistical motion analyses. Here, we propose a new stochastic model for the orientation angle of a bacterium and formulate conditonal moments, which we determine both in theory and from experimental trajectories. This provides an alternative way of quantifying the bacterial run-and-tumble strategy and of recognizing tumble events. Our approach no longer relies on arbitrarily chosen segmentation thresholds and rigorously quantifies the uncertainty in tumble recognition. We successfully apply our method not only to the paradigmatic case of E. coli but also to trajectories of the soil bacterium Pseudomonas putida, demonstrating that our approach provides a novel way to reliably characterize the tumbling statistics and chemotaxis strategies of bacterial swimmers across different species.

Rac signal adaptation controls neutrophil mobilization from the bone marrow

Mobilization of neutrophils from the bone marrow determines neutrophil blood counts and thus is medically important. Balanced neutrophil mobilization from the bone marrow depends on the retention-promoting chemokine CXCL12 and its receptor CXCR4 and the egression-promoting chemokine CXCL2 and its receptor CXCR2. Both pathways activate the small guanosine triphosphatase Rac, leaving the role of this signaling event in neutrophil retention and egression ambiguous. On the assumption that active Rac determines persistent directional cell migration, we generated a mathematical model to link chemokine-mediated Rac modulation to neutrophil egression time. Our computer simulation indicated that, in the bone marrow, where the retention signal predominated, egression time strictly depended on the time it took Rac to return to its basal activity (namely, adaptation). This prediction was validated in mice lacking the Rac inhibitor ArhGAP15. Neutrophils in these mice showed prolonged Rac adaptation and cell-autonomous retention in the bone marrow. Our model thus demonstrates that mobilization in the presence of two spatially defined opposing chemotactic cues strictly depends on inhibitors shaping the time course of signal adaptation. Furthermore, our findings might help to find new modes of intervention to treat conditions characterized by excessively low or high circulating neutrophils.

Non-genetic diversity modulates population performance

Biological functions are typically performed by groups of cells that express predominantly the same genes, yet display a continuum of phenotypes. While it is known how one genotype can generate such non-genetic diversity, it remains unclear how different phenotypes contribute to the performance of biological function at the population level. We developed a microfluidic device to simultaneously measure the phenotype and chemotactic performance of tens of thousands of individual, freely swimming Escherichia coli as they climbed a gradient of attractant. We discovered that spatial structure spontaneously emerged from initially well-mixed wild-type populations due to non-genetic diversity. By manipulating the expression of key chemotaxis proteins, we established a causal relationship between protein expression, non-genetic diversity, and performance that was theoretically predicted. This approach generated a complete phenotype-to-performance map, in which we found a nonlinear regime. We used this map to demonstrate how changing the shape of a phenotypic distribution can have as large of an effect on collective performance as changing the mean phenotype, suggesting that selection could act on both during the process of adaptation.

Reconstruction of Spatial Thermal Gradient Encoded in Thermosensory Neuron AFD in Caenorhabditis elegans

During navigation, animals process temporal sequences of sensory inputs to evaluate the surrounding environment. Thermotaxis of Caenorhabditis elegans is a favorable sensory behavior to elucidate how navigating animals process sensory signals from the environment. Sensation and storage of temperature information by a bilaterally symmetric pair of thermosensory neurons, AFD, is essential for the animals to migrate toward the memorized temperature on a thermal gradient. However, the encoding mechanisms of the spatial environment with the temporal AFD activity during navigation remain to be elucidated. Here, we show how the AFD neuron encodes sequences of sensory inputs to perceive spatial thermal environment. We used simultaneous calcium imaging and tracking system for a freely moving animal and characterized the response property of AFD to the thermal stimulus during thermotaxis. We show that AFD neurons respond to shallow temperature increases with intermittent calcium pulses and detect temperature differences with a critical time window of 20 s, which is similar to the timescale of behavioral elements of C. elegans, such as turning. Convolution of a thermal stimulus and the identified response property successfully reconstructs AFD activity. Conversely, deconvolution of the identified response kernel and AFD activity reconstructs the shallow thermal gradient with migration trajectory, indicating that AFD activity and the migration trajectory are sufficient as the encoded signals for thermal environment. Our study demonstrates bidirectional transformation between environmental thermal information and encoded neural activity.

SIGNIFICANCE STATEMENT

Deciphering how information is encoded in the nervous system is an important challenge for understanding the principles of information processing in neural circuits. During navigation behavior, animals transform spatial information to temporal patterns of neural activity. To elucidate how a sensory system achieves this transformation, we focused on a thermosensory neuron in Caenorhabditis elegans called AFD, which plays a major role in a sensory behavior. Using tracking and calcium imaging system for freely moving animals, we identified the response property of the AFD. The identified response property enabled us to reconstruct both neural activity from a temperature stimulus and a spatial thermal environment from neural activity. These results shed light on how a sensory system encodes the environment.

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.

The Role of Adaptation in Bacterial Speed Races

Evolution of biological sensory systems is driven by the need for efficient responses to environmental stimuli. A paradigm among prokaryotes is the chemotaxis system, which allows bacteria to navigate gradients of chemoattractants by biasing their run-and-tumble motion. A notable feature of chemotaxis is adaptation: after the application of a step stimulus, the bacterial running time relaxes to its pre-stimulus level. The response to the amino acid aspartate is precisely adapted whilst the response to serine is not, in spite of the same pathway processing the signals preferentially sensed by the two receptors Tar and Tsr, respectively. While the chemotaxis pathway in E. coli is well characterized, the role of adaptation, its functional significance and the ecological conditions where chemotaxis is selected, are largely unknown. Here, we investigate the role of adaptation in the climbing of gradients by E. coli. We first present theoretical arguments that highlight the mechanisms that control the efficiency of the chemotactic up-gradient motion. We discuss then the limitations of linear response theory, which motivate our subsequent experimental investigation of E. coli speed races in gradients of aspartate, serine and combinations thereof. By using microfluidic techniques, we engineer controlled gradients and demonstrate that bacterial fronts progress faster in equal-magnitude gradients of serine than aspartate. The effect is observed over an extended range of concentrations and is not due to differences in swimming velocities. We then show that adding a constant background of serine to gradients of aspartate breaks the adaptation to aspartate, which results in a sped-up progression of the fronts and directly illustrate the role of adaptation in chemotactic gradient-climbing.

Biological sensory pathways are presumed to evolve for the processing of environmental information, yet quantitative evidence is scant. Chemotaxis allows bacteria to sense chemical gradients but their ecological distribution, e.g. whether natural gradients sensed by E. coli change slowly or rapidly in space and time, is unknown. That distribution matters, as it controls constraints and selective pressure acting on the pathway. We used microfluidic devices to generate controlled chemoattractant gradients and measure the speed of bacterial climbing of those gradients. We could thereby assay the impact of adaptation properties of the chemotaxis pathway onto the progression of gradient climbing. We specifically show that loss of adaptation, induced by adding a background of serine to gradients of aspartate, leads to a faster progression of the bacteria along the chemoattractant gradient. We finally discuss why our experiments suggest that ecological conditions are likely to involve chemoattractant profiles more complex than constant gradients usually considered in the laboratory.

Bacterial chemotaxis: information processing, thermodynamics, and behavior

Escherichia coli has long been used as a model organism due to the extensive experimental characterization of its pathways and molecular components. Take chemotaxis as an example, which allows bacteria to sense and swim in response to chemicals, such as nutrients and toxins. Many of the pathway's remarkable sensing and signaling properties are now concisely summarized in terms of design (or engineering) principles. More recently, new approaches from information theory and stochastic thermodynamics have begun to address how pathways process environmental stimuli and what the limiting factors are. However, to fully capitalize on these theoretical advances, a closer connection with single-cell experiments will be required.

Distinct predictive performance of Rac1 and Cdc42 in cell migration

We propose a new computation-based approach for elucidating how signaling molecules are decoded in cell migration. In this approach, we performed FRET time-lapse imaging of Rac1 and Cdc42, members of Rho GTPases which are responsible for cell motility, and quantitatively identified the response functions that describe the conversion from the molecular activities to the morphological changes. Based on the identified response functions, we clarified the profiles of how the morphology spatiotemporally changes in response to local and transient activation of Rac1 and Cdc42, and found that Rac1 and Cdc42 activation triggers laterally propagating membrane protrusion. The response functions were also endowed with property of differentiator, which is beneficial for maintaining sensitivity under adaptation to the mean level of input. Using the response function, we could predict the morphological change from molecular activity, and its predictive performance provides a new quantitative measure of how much the Rho GTPases participate in the cell migration. Interestingly, we discovered distinct predictive performance of Rac1 and Cdc42 depending on the migration modes, indicating that Rac1 and Cdc42 contribute to persistent and random migration, respectively. Thus, our proposed predictive approach enabled us to uncover the hidden information processing rules of Rho GTPases in the cell migration.

High-throughput 3D tracking of bacteria on a standard phase contrast microscope

Bacteria employ diverse motility patterns in traversing complex three-dimensional (3D) natural habitats. 2D microscopy misses crucial features of 3D behaviour, but the applicability of existing 3D tracking techniques is constrained by their performance or ease of use. Here we present a simple, broadly applicable, high-throughput 3D bacterial tracking method for use in standard phase contrast microscopy. Bacteria are localized at micron-scale resolution over a range of 350 × 300 × 200 μm by maximizing image cross-correlations between their observed diffraction patterns and a reference library. We demonstrate the applicability of our technique to a range of bacterial species and exploit its high throughput to expose hidden contributions of bacterial individuality to population-level variability in motile behaviour. The simplicity of this powerful new tool for bacterial motility research renders 3D tracking accessible to a wider community and paves the way for investigations of bacterial motility in complex 3D environments.

Microscopy techniques used to study the movement of swimming microbes are limited to two dimensions or require sophisticated devices. Here, Taute et al. present a simple method for high-throughput 3D tracking of bacteria using standard phase contrast microscopy.

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.

Variation of swimming speed enhances the chemotactic migration of Escherichia coli

Studies on chemotaxis of Escherichia coli have shown that modulation of tumble frequency causes a net drift up the gradient of attractants. Recently, it has been demonstrated that the bacteria is also capable of varying its runs speed in uniform concentration of attractant. In this study, we investigate the role of swimming speed on the chemotactic migration of bacteria. To this end, cells are exposed to gradients of a non-metabolizable analogue of glucose which are sensed via the Trg sensor. When exposed to a gradient, the cells modulate their tumble duration, which is accompanied with variation in swimming speed leading to drift velocities that are much higher than those achieved through the modulation of the tumble duration alone. We use an existing intra-cellular model developed for the Tar receptor and incorporate the variation of the swimming speed along with modulation of tumble frequency to predict drift velocities close to the measured values. The main implication of our study is that E. coli not only modulates the tumble frequency, but may also vary the swimming speed to affect chemotaxis and thereby efficiently sample its nutritionally rich environment.

Estimating cellular parameters through optimization procedures: elementary principles and applications

Construction of quantitative models is a primary goal of quantitative biology, which aims to understand cellular and organismal phenomena in a quantitative manner. In this article, we introduce optimization procedures to search for parameters in a quantitative model that can reproduce experimental data. The aim of optimization is to minimize the sum of squared errors (SSE) in a prediction or to maximize likelihood. A (local) maximum of likelihood or (local) minimum of the SSE can efficiently be identified using gradient approaches. Addition of a stochastic process enables us to identify the global maximum/minimum without becoming trapped in local maxima/minima. Sampling approaches take advantage of increasing computational power to test numerous sets of parameters in order to determine the optimum set. By combining Bayesian inference with gradient or sampling approaches, we can estimate both the optimum parameters and the form of the likelihood function related to the parameters. Finally, we introduce four examples of research that utilize parameter optimization to obtain biological insights from quantified data: transcriptional regulation, bacterial chemotaxis, morphogenesis, and cell cycle regulation. With practical knowledge of parameter optimization, cell and developmental biologists can develop realistic models that reproduce their observations and thus, obtain mechanistic insights into phenomena of interest.

Fast, high-throughput measurement of collective behaviour in a bacterial population

Swimming bacteria explore their environment by performing a random walk, which is biased in response to, for example, chemical stimuli, resulting in a collective drift of bacterial populations towards 'a better life'. This phenomenon, called chemotaxis, is one of the best known forms of collective behaviour in bacteria, crucial for bacterial survival and virulence. Both single-cell and macroscopic assays have investigated bacterial behaviours. However, theories that relate the two scales have previously been difficult to test directly. We present an image analysis method, inspired by light scattering, which measures the average collective motion of thousands of bacteria simultaneously. Using this method, a time-varying collective drift as small as 50 nm s(-1) can be measured. The method, validated using simulations, was applied to chemotactic Escherichia coli bacteria in linear gradients of the attractant α-methylaspartate. This enabled us to test a coarse-grained minimal model of chemotaxis. Our results clearly map the onset of receptor methylation, and the transition from linear to logarithmic sensing in the bacterial response to an external chemoeffector. Our method is broadly applicable to problems involving the measurement of collective drift with high time resolution, such as cell migration and fluid flows measurements, and enables fast screening of tactic behaviours.

Universal response-adaptation relation in bacterial chemotaxis

The bacterial strategy of chemotaxis relies on temporal comparisons of chemical concentrations, where the probability of maintaining the current direction of swimming is modulated by changes in stimulation experienced during the recent past. A short-term memory required for such comparisons is provided by the adaptation system, which operates through the activity-dependent methylation of chemotaxis receptors. Previous theoretical studies have suggested that efficient navigation in gradients requires a well-defined adaptation rate, because the memory time scale needs to match the duration of straight runs made by bacteria. Here we demonstrate that the chemotaxis pathway of Escherichia coli does indeed exhibit a universal relation between the response magnitude and adaptation time which does not depend on the type of chemical ligand. Our results suggest that this alignment of adaptation rates for different ligands is achieved through cooperative interactions among chemoreceptors rather than through fine-tuning of methylation rates for individual receptors. This observation illustrates a yet-unrecognized function of receptor clustering in bacterial chemotaxis.

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.

Variation in swimming speed of Escherichia coli in response to attractant

It is well known that Escherichia coli executes chemotactic motion in response to chemical cues by modulating the flagellar motor bias alone. However, previous studies have reported the possibility of variation in run speed in the presence of attractants although it is unclear whether bacteria can deliberately modulate their swimming speeds in response to environmental cues or if the motor speeds are hardwired. By studying the detailed motion of cells in a uniform concentration of glucose and its non-metabolizable analogue, we show that changing concentrations may be accompanied by variation in the swimming speed. For a fixed run duration, cells exposed to the attractants achieved a higher peak-swimming speed after a tumble compared with that in plain motility buffer. Our experiments using the mutant strain lacking the Trg sensor show no change in swimming speed with varying concentrations of the non-metabolizable analogue, suggesting that sensing may play a role in the observed variation of swimming speed.

Microfluidics expanding the frontiers of microbial ecology

Microfluidics has significantly contributed to the expansion of the frontiers of microbial ecology over the past decade by allowing researchers to observe the behaviors of microbes in highly controlled microenvironments, across scales from a single cell to mixed communities. Spatially and temporally varying distributions of organisms and chemical cues that mimic natural microbial habitats can now be established by exploiting physics at the micrometer scale and by incorporating structures with specific geometries and materials. In this article, we review applications of microfluidics that have resulted in insightful discoveries on fundamental aspects of microbial life, ranging from growth and sensing to cell-cell interactions and population dynamics. We anticipate that this flexible multidisciplinary technology will continue to facilitate discoveries regarding the ecology of microorganisms and help uncover strategies to control microbial processes such as biofilm formation and antibiotic resistance.

Logarithmic and power law input-output relations in sensory systems with fold-change detection

Two central biophysical laws describe sensory responses to input signals. One is a logarithmic relationship between input and output, and the other is a power law relationship. These laws are sometimes called the Weber-Fechner law and the Stevens power law, respectively. The two laws are found in a wide variety of human sensory systems including hearing, vision, taste, and weight perception; they also occur in the responses of cells to stimuli. However the mechanistic origin of these laws is not fully understood. To address this, we consider a class of biological circuits exhibiting a property called fold-change detection (FCD). In these circuits the response dynamics depend only on the relative change in input signal and not its absolute level, a property which applies to many physiological and cellular sensory systems. We show analytically that by changing a single parameter in the FCD circuits, both logarithmic and power-law relationships emerge; these laws are modified versions of the Weber-Fechner and Stevens laws. The parameter that determines which law is found is the steepness (effective Hill coefficient) of the effect of the internal variable on the output. This finding applies to major circuit architectures found in biological systems, including the incoherent feed-forward loop and nonlinear integral feedback loops. Therefore, if one measures the response to different fold changes in input signal and observes a logarithmic or power law, the present theory can be used to rule out certain FCD mechanisms, and to predict their cooperativity parameter. We demonstrate this approach using data from eukaryotic chemotaxis signaling.

One of the first measurements an experimentalist makes to understand a sensory system is to explore the relation between input signal and the systems response amplitude. Here, we show using mathematical models that this measurement can give important clues about the possible mechanism of sensing. We use models that incorporate the nearly-universal features of sensory systems, including hearing and vision, and the sensing pathways of individual cells. These nearly-universal features include exact adaptation-the ability to ignore prolonged input stimuli and return to basal activity, and fold-change detection- response to relative changes in input, not absolute changes. Together with information on the input-output relationship-e.g. is it a logarithmic or a power law relationship-we show that these conditions provide enough constraints to allow the researcher to reject certain circuit designs; it also predicts, if one assumes a given design, one of its key parameters. This study can thus help unify our understanding of sensory systems, and help pinpoint the possible biological circuits based on physiological measurements.

Limits of feedback control in bacterial chemotaxis

Inputs to signaling pathways can have complex statistics that depend on the environment and on the behavioral response to previous stimuli. Such behavioral feedback is particularly important in navigation. Successful navigation relies on proper coupling between sensors, which gather information during motion, and actuators, which control behavior. Because reorientation conditions future inputs, behavioral feedback can place sensors and actuators in an operational regime different from the resting state. How then can organisms maintain proper information transfer through the pathway while navigating diverse environments? In bacterial chemotaxis, robust performance is often attributed to the zero integral feedback control of the sensor, which guarantees that activity returns to resting state when the input remains constant. While this property provides sensitivity over a wide range of signal intensities, it remains unclear how other parameters such as adaptation rate and adapted activity affect chemotactic performance, especially when considering that the swimming behavior of the cell determines the input signal. We examine this issue using analytical models and simulations that incorporate recent experimental evidences about behavioral feedback and flagellar motor adaptation. By focusing on how sensory information carried by the response regulator is best utilized by the motor, we identify an operational regime that maximizes drift velocity along chemical concentration gradients for a wide range of environments and sensor adaptation rates. This optimal regime is outside the dynamic range of the motor response, but maximizes the contrast between run duration up and down gradients. In steep gradients, the feedback from chemotactic drift can push the system through a bifurcation. This creates a non-chemotactic state that traps cells unless the motor is allowed to adapt. Although motor adaptation helps, we find that as the strength of the feedback increases individual phenotypes cannot maintain the optimal operational regime in all environments, suggesting that diversity could be beneficial.

The biased random walk is a fundamental strategy used by many organisms to navigate their environment. Drift along the desired direction is achieved by reducing the probability to reorient whenever conditions improve. In the chemotaxis system of Escherichia coli, this is accomplished with a sensory module that implements negative integral feedback control, the output of which is relayed to the flagellar motors (the actuators) by a response regulator to control the probability to change direction. The proper dynamical coupling between sensor and actuator is critical for the performance of the random walker. Here, we identify an optimal regime for this coupling that maximizes drift velocity in the direction of the gradient in multiple environments. Our analysis reveals that feedback of the behavior onto the system in steep gradients can constrain individual cell performance, by causing bi-stable behavior that can trap cells in non-chemotactic states. These limitations are inherent in the biased random walk strategy with integral feedback control, but can be alleviated if the output of the pathway adapts, as recently characterized for the flagellar motors in Escherichia coli.

A bacterial swimmer with two alternating speeds of propagation

We recorded large data sets of swimming trajectories of the soil bacterium Pseudomonas putida. Like other prokaryotic swimmers, P. putida exhibits a motion pattern dominated by persistent runs that are interrupted by turning events. An in-depth analysis of their swimming trajectories revealed that the majority of the turning events is characterized by an angle of ϕ1 = 180° (reversals). To a lesser extent, turning angles of ϕ2 = 0° are also found. Remarkably, we observed that, upon a reversal, the swimming speed changes by a factor of two on average-a prominent feature of the motion pattern that, to our knowledge, has not been reported before. A theoretical model, based on the experimental values for the average run time and the rotational diffusion, recovers the mean-square displacement of P. putida if the two distinct swimming speeds are taken into account. Compared to a swimmer that moves with a constant intermediate speed, the mean-square displacement is strongly enhanced. We furthermore observed a negative dip in the directional autocorrelation at intermediate times, a feature that is only recovered in an extended model, where the nonexponential shape of the run-time distribution is taken into account.

Synthetic circuit for exact adaptation and fold-change detection

Biological organisms use their sensory systems to detect changes in their environment. The ability of sensory systems to adapt to static inputs allows wide dynamic range as well as sensitivity to input changes including fold-change detection, a response that depends only on fold changes in input, and not on absolute changes. This input scale invariance underlies an important strategy for search that depends solely on the spatial profile of the input. Synthetic efforts to reproduce the architecture and response of cellular circuits provide an important step to foster understanding at the molecular level. We report the bottom-up assembly of biochemical systems that show exact adaptation and fold-change detection. Using a malachite green aptamer as the output, a synthetic transcriptional circuit with the connectivity of an incoherent feed-forward loop motif exhibits pulse generation and exact adaptation. A simple mathematical model was used to assess the amplitude and duration of pulse response as well as the parameter regimes required for fold-change detection. Upon parameter tuning, this synthetic circuit exhibits fold-change detection for four successive rounds of two-fold input changes. The experimental realization of fold-change detection circuit highlights the programmability of transcriptional switches and the ability to obtain predictive dynamical systems in a cell-free environment for technological applications.

Efficient fold-change detection based on protein-protein interactions

Various biological sensory systems exhibit a response to a relative change of the stimulus, often referred to as fold-change detection. In the past few years, fold-change detecting mechanisms, based on transcriptional networks, have been proposed. Here we present a fold-change detecting mechanism, based on protein-protein interactions, consisting of two interacting proteins. This mechanism does not consume chemical energy and is not subject to transcriptional and translational noise, in contrast to previously proposed mechanisms. We show by analytical and numerical calculations that the mechanism is robust and can have a fast, precise, and efficient response for parameters that are relevant to eukaryotic cells.

How the motility pattern of bacteria affects their dispersal and chemotaxis

Most bacteria at certain stages of their life cycle are able to move actively; they can swim in a liquid or crawl on various surfaces. A typical path of the moving cell often resembles the trajectory of a random walk. However, bacteria are capable of modifying their apparently random motion in response to changing environmental conditions. As a result, bacteria can migrate towards the source of nutrients or away from harmful chemicals. Surprisingly, many bacterial species that were studied have several distinct motility patterns, which can be theoretically modeled by a unifying random walk approach. We use this approach to quantify the process of cell dispersal in a homogeneous environment and show how the bacterial drift velocity towards the source of attracting chemicals is affected by the motility pattern of the bacteria. Our results open up the possibility of accessing additional information about the intrinsic response of the cells using macroscopic observations of bacteria moving in inhomogeneous environments.

Evaluation of bacterial motility from non-Gaussianity of finite-sample trajectories using the large deviation principle

Motility of bacteria is usually recognized in the trajectory data and compared with Brownian motion, but the diffusion coefficient is insufficient to evaluate it. In this paper, we propose a method based on the large deviation principle. We show that it can be used to evaluate the non-Gaussian characteristics of model Escherichia coli motions and to distinguish combinations of the mean running duration and running speed that lead to the same diffusion coefficient. Our proposed method does not require chemical stimuli to induce the chemotaxis in a specific direction, and it is applicable to various types of self-propelling motions for which no a priori information of, for example, threshold parameters for run and tumble or head/tail direction is available. We also address the issue of the finite-sample effect on the large deviation quantities, but we propose to make use of it to characterize the nature of motility.