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

Predicted Auxiliary Navigation Mechanism of Peritrichously Flagellated Chemotactic Bacteria

Chemotactic movement of Escherichia coli is one of the most thoroughly studied paradigms of simple behavior. Due to significant competitive advantage conferred by chemotaxis and to high evolution rates in bacteria, the chemotaxis system is expected to be strongly optimized. Bacteria follow gradients by performing temporal comparisons of chemoeffector concentrations along their runs, a strategy which is most efficient given their size and swimming speed. Concentration differences are detected by a sensory system and transmitted to modulate rotation of flagellar motors, decreasing the probability of a tumble and reorientation if the perceived concentration change during a run is positive. Such regulation of tumble probability is of itself sufficient to explain chemotactic drift of a population up the gradient, and is commonly assumed to be the only navigation mechanism of chemotactic E. coli. Here we use computer simulations to predict existence of an additional mechanism of gradient navigation in E. coli. Based on the experimentally observed dependence of cell tumbling angle on the number of switching motors, we suggest that not only the tumbling probability but also the degree of reorientation during a tumble depend on the swimming direction along the gradient. Although the difference in mean tumbling angles up and down the gradient predicted by our model is small, it results in a dramatic enhancement of the cellular drift velocity along the gradient. We thus demonstrate a new level of optimization in E. coli chemotaxis, which arises from the switching of several flagellar motors and a resulting fine tuning of tumbling angle. Similar strategy is likely to be used by other peritrichously flagellated bacteria, and indicates yet another level of evolutionary development of bacterial chemotaxis

Feedback control architecture & the bacterial chemotaxis network

Bacteria move towards favourable and away from toxic environments by changing their swimming pattern. This response is regulated by the chemotaxis signalling pathway, which has an important feature: it uses feedback to ‘reset’ (adapt) the bacterial sensing ability, which allows the bacteria to sense a range of background environmental changes. The role of this feedback has been studied extensively in the simple chemotaxis pathway of Escherichia coli. However it has been recently found that the majority of bacteria have multiple chemotaxis homologues of the E. coli proteins, resulting in more complex pathways. In this paper we investigate the configuration and role of feedback in Rhodobacter sphaeroides, a bacterium containing multiple homologues of the chemotaxis proteins found in E. coli. Multiple proteins could produce different possible feedback configurations, each having different chemotactic performance qualities and levels of robustness to variations and uncertainties in biological parameters and to intracellular noise. We develop four models corresponding to different feedback configurations. Using a series of carefully designed experiments we discriminate between these models and invalidate three of them. When these models are examined in terms of robustness to noise and parametric uncertainties, we find that the non-invalidated model is superior to the others. Moreover, it has a ‘cascade control’ feedback architecture which is used extensively in engineering to improve system performance, including robustness. Given that the majority of bacteria are known to have multiple chemotaxis pathways, in this paper we show that some feedback architectures allow them to have better performance than others. In particular, cascade control may be an important feature in achieving robust functionality in more complex signalling pathways and in improving their performance

Trade-off between Responsiveness and Noise Suppression in Biomolecular System Responses to Environmental Cues

When living systems detect changes in their external environment their response must be measured to balance the need to react appropriately with the need to remain stable, ignoring insignificant signals. Because this is a fundamental challenge of all biological systems that execute programs in response to stimuli, we developed a generalized time-frequency analysis (TFA) framework to systematically explore the dynamical properties of biomolecular networks. Using TFA, we focused on two well-characterized yeast gene regulatory networks responsive to carbon-source shifts and a mammalian innate immune regulatory network responsive to lipopolysaccharides (LPS). The networks are comprised of two different basic architectures. Dual positive and negative feedback loops make up the yeast galactose network; whereas overlapping positive and negative feed-forward loops are common to the yeast fatty-acid response network and the LPS-induced network of macrophages. TFA revealed remarkably distinct network behaviors in terms of trade-offs in responsiveness and noise suppression that are appropriately tuned to each biological response. The wild type galactose network was found to be highly responsive while the oleate network has greater noise suppression ability. The LPS network appeared more balanced, exhibiting less bias toward noise suppression or responsiveness. Exploration of the network parameter space exposed dramatic differences in system behaviors for each network. These studies highlight fundamental structural and dynamical principles that underlie each network, reveal constrained parameters of positive and negative feedback and feed-forward strengths that tune the networks appropriately for their respective biological roles, and demonstrate the general utility of the TFA approach for systems and synthetic biology

A Bayesian perspective on classical control

The connections between optimal control and Bayesian inference have long been recognised, with the field of stochastic (optimal) control combining these frameworks for the solution of partially observable control problems. In particular, for the linear case with quadratic functions and Gaussian noise, stochastic control has shown remarkable results in different fields, including robotics, reinforcement learning and neuroscience, especially thanks to the established duality of estimation and control processes. Following this idea we recently introduced a formulation of PID control, one of the most popular methods from classical control, based on active inference, a theory with roots in variational Bayesian methods, and applications in the biological and neural sciences. In this work, we highlight the advantages of our previous formulation and introduce new and more general ways to tackle some existing problems in current controller design procedures. In particular, we consider 1) a gradient-based tuning rule for the parameters (or gains) of a PID controller, 2) an implementation of multiple degrees of freedom for independent responses to different types of signals (e.g., two-degree-of-freedom PID), and 3) a novel time-domain formalisation of the performance-robustness trade-off in terms of tunable constraints (i.e., priors in a Bayesian model) of a single cost functional, variational free energy.Comment: 8 pages, Accepted at IJCNN 202

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 affect chemotactic performance, especially when considering that the swimming behavior of the cell determines the input signal. Using analytical models and simulations that incorporate recent experimental evidences about behavioral feedback and flagellar motor adaptation we identify an operational regime of the pathway that maximizes drift velocity for various 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.Comment: Corrected one typo. First two authors contributed equally. Notably, there were various typos in the values of the parameters in the model of motor adaptation. The results remain unchange

Quantitative Modeling of Escherichia coli Chemotactic Motion in Environments Varying in Space and Time

Escherichia coli chemotactic motion in spatiotemporally varying environments is studied by using a computational model based on a coarse-grained description of the intracellular signaling pathway dynamics. We find that the cell's chemotaxis drift velocity vd is a constant in an exponential attractant concentration gradient [L]∝exp(Gx). vd depends linearly on the exponential gradient G before it saturates when G is larger than a critical value GC. We find that GC is determined by the intracellular adaptation rate kR with a simple scaling law: . The linear dependence of vd on G = d(ln[L])/dx directly demonstrates E. coli's ability in sensing the derivative of the logarithmic attractant concentration. The existence of the limiting gradient GC and its scaling with kR are explained by the underlying intracellular adaptation dynamics and the flagellar motor response characteristics. For individual cells, we find that the overall average run length in an exponential gradient is longer than that in a homogeneous environment, which is caused by the constant kinase activity shift (decrease). The forward runs (up the gradient) are longer than the backward runs, as expected; and depending on the exact gradient, the (shorter) backward runs can be comparable to runs in a spatially homogeneous environment, consistent with previous experiments. In (spatial) ligand gradients that also vary in time, the chemotaxis motion is damped as the frequency ω of the time-varying spatial gradient becomes faster than a critical value ωc, which is controlled by the cell's chemotaxis adaptation rate kR. Finally, our model, with no adjustable parameters, agrees quantitatively with the classical capillary assay experiments where the attractant concentration changes both in space and time. Our model can thus be used to study E. coli chemotaxis behavior in arbitrary spatiotemporally varying environments. Further experiments are suggested to test some of the model predictions

Precise Regulation of Gene Expression Dynamics Favors Complex Promoter Architectures

Promoters process signals through recruitment of transcription factors and RNA polymerase, and dynamic changes in promoter activity constitute a major noise source in gene expression. However, it is barely understood how complex promoter architectures determine key features of promoter dynamics. Here, we employ prototypical promoters of yeast ribosomal protein genes as well as simplified versions thereof to analyze the relations among promoter design, complexity, and function. These promoters combine the action of a general regulatory factor with that of specific transcription factors, a common motif of many eukaryotic promoters. By comprehensively analyzing stationary and dynamic promoter properties, this model-based approach enables us to pinpoint the structural characteristics underlying the observed behavior. Functional tradeoffs impose constraints on the promoter architecture of ribosomal protein genes. We find that a stable scaffold in the natural design results in low transcriptional noise and strong co-regulation of target genes in the presence of gene silencing. This configuration also exhibits superior shut-off properties, and it can serve as a tunable switch in living cells. Model validation with independent experimental data suggests that the models are sufficiently realistic. When combined, our results offer a mechanistic explanation for why specific factors are associated with low protein noise in vivo. Many of these findings hold for a broad range of model parameters and likely apply to other eukaryotic promoters of similar structure

Implementing Arithmetic and Other Analytic Operations By Transcriptional Regulation

The transcriptional regulatory machinery of a gene can be viewed as a computational device, with transcription factor concentrations as inputs and expression level as the output. This view begs the question: what kinds of computations are possible? We show that different parameterizations of a simple chemical kinetic model of transcriptional regulation are able to approximate all four standard arithmetic operations: addition, subtraction, multiplication, and division, as well as various equality and inequality operations. This contrasts with other studies that emphasize logical or digital notions of computation in biological networks. We analyze the accuracy and precision of these approximations, showing that they depend on different sets of parameters, and are thus independently tunable. We demonstrate that networks of these “arithmetic” genes can be combined to accomplish yet more complicated computations by designing and simulating a network that detects statistically significant elevations in a time-varying signal. We also consider the much more general problem of approximating analytic functions, showing that this can be achieved by allowing multiple transcription factor binding sites on the promoter. These observations are important for the interpretation of naturally occurring networks and imply new possibilities for the design of synthetic networks