Quantitative modeling of sensitivity in bacterial chemotaxis: the role of coupling among different chemoreceptor species

Lattice ultrasensitivity produces large gain in E. coli chemosensing

E. coli use a regular lattice of receptors and attached kinases to detect and amplify faint chemical signals. Kinase output is characterized by precise adaptation to a wide range of background ligand levels and large gain in response to small relative changes in ligand concentration. These characteristics are well described by models which achieve their gain through equilibrium cooperativity. But these models are challenged by two experimental results. First, neither adaptation nor large gain are present in receptor binding assays. Second, in cells lacking adaptation machinery, fluctuations can sometimes be enormous, with essentially all kinases transitioning together. Here we introduce a far-from equilibrium model in which receptors gate the spread of activity between neighboring kinases. This model achieves large gain through a mechanism we term lattice ultrasensitivity (LU). In our LU model, kinase and receptor states are separate degrees of freedom, and kinase kinetics are dominated by chemical rates far-from-equilibrium rather than by equilibrium allostery. The model recapitulates the successes of past models, but also matches the challenging experimental findings. Importantly, unlike past lattice critical models, our LU model does not require parameters to be fine tuned for function.

A nonequilibrium allosteric model for receptor-kinase complexes: The role of energy dissipation in chemotaxis signaling

The Escherichia coli chemotaxis signaling pathway has served as a model system for the adaptive sensing of environmental signals by large protein complexes. The chemoreceptors control the kinase activity of CheA in response to the extracellular ligand concentration and adapt across a wide concentration range by undergoing methylation and demethylation. Methylation shifts the kinase response curve by orders of magnitude in ligand concentration while incurring a much smaller change in the ligand binding curve. Here, we show that the disproportionate shift in binding and kinase response is inconsistent with equilibrium allosteric models. To resolve this inconsistency, we present a nonequilibrium allosteric model that explicitly includes the dissipative reaction cycles driven by adenosine triphosphate (ATP) hydrolysis. The model successfully explains all existing joint measurements of ligand binding, receptor conformation, and kinase activity for both aspartate and serine receptors. Our results suggest that the receptor complex acts as an enzyme: Receptor methylation modulates the ON-state kinetics of the kinase (e.g., phosphorylation rate), while ligand binding controls the equilibrium balance between kinase ON/OFF states. Furthermore, sufficient energy dissipation is responsible for maintaining and enhancing the sensitivity range and amplitude of the kinase response. We demonstrate that the nonequilibrium allosteric model is broadly applicable to other sensor-kinase systems by successfully fitting previously unexplained data from the DosP bacterial oxygen-sensing system. Overall, this work provides a nonequilibrium physics perspective on cooperative sensing by large protein complexes and opens up research directions for understanding their microscopic mechanisms through simultaneous measurements and modeling of ligand binding and downstream responses.

Evolutionary Stability of Small Molecular Regulatory Networks That Exhibit Near-Perfect Adaptation

Large-scale protein regulatory networks, such as signal transduction systems, contain small-scale modules ('motifs') that carry out specific dynamical functions. Systematic characterization of the properties of small network motifs is therefore of great interest to molecular systems biologists. We simulate a generic model of three-node motifs in search of near-perfect adaptation, the property that a system responds transiently to a change in an environmental signal and then returns near-perfectly to its pre-signal state (even in the continued presence of the signal). Using an evolutionary algorithm, we search the parameter space of these generic motifs for network topologies that score well on a pre-defined measure of near-perfect adaptation. We find many high-scoring parameter sets across a variety of three-node topologies. Of all possibilities, the highest scoring topologies contain incoherent feed-forward loops (IFFLs), and these topologies are evolutionarily stable in the sense that, under 'macro-mutations' that alter the topology of a network, the IFFL motif is consistently maintained. Topologies that rely on negative feedback loops with buffering (NFLBs) are also high-scoring; however, they are not evolutionarily stable in the sense that, under macro-mutations, they tend to evolve an IFFL motif and may-or may not-lose the NFLB motif.

Protein-protein complexes can undermine ultrasensitivity-dependent biological adaptation

Robust perfect adaptation (RPA) is a ubiquitously observed signalling response across all scales of biological organization. A major class of network architectures that drive RPA in complex networks is the Opposer module-a feedback-regulated network into which specialized integral-computing 'opposer node(s)' are embedded. Although ultrasensitivity-generating chemical reactions have long been considered a possible mechanism for such adaptation-conferring opposer nodes, this hypothesis has relied on simplified Michaelian models, which neglect the presence of protein-protein complexes. Here we develop complex-complete models of interlinked covalent-modification cycles with embedded ultrasensitivity, explicitly capturing all molecular interactions and protein complexes. Strikingly, we demonstrate that the presence of protein-protein complexes thwarts the network's capacity for RPA in any 'free' active protein form, conferring RPA capacity instead on the concentration of a larger protein pool consisting of two distinct forms of a single protein. We further show that the presence of enzyme-substrate complexes, even at comparatively low concentrations, play a crucial and previously unrecognized role in controlling the RPA response-significantly reducing the range of network inputs for which RPA can obtain, and imposing greater parametric requirements on the RPA response. These surprising results raise fundamental new questions as to the biochemical requirements for adaptation-conferring Opposer modules within complex cellular networks.

Collective behavior and nongenetic inheritance allow bacterial populations to adapt to changing environments

Collective behaviors require coordination among a group of individuals. As a result, individuals that are too phenotypically different from the rest of the group can be left out, reducing heterogeneity, but increasing coordination. If individuals also reproduce, the offspring can have different phenotypes from their parent(s). This raises the question of how these two opposing processes-loss of diversity by collective behaviors and generation of it through growth and inheritance-dynamically shape the phenotypic composition of an isogenic population. We examine this question theoretically using collective migration of chemotactic bacteria as a model system, where cells of different swimming phenotypes are better suited to navigate in different environments. We find that the differential loss of phenotypes caused by collective migration is environment-dependent. With cell growth, this differential loss enables migrating populations to dynamically adapt their phenotype compositions to the environment, enhancing migration through multiple environments. Which phenotypes are produced upon cell division depends on the level of nongenetic inheritance, and higher inheritance leads to larger composition adaptation and faster migration at steady state. However, this comes at the cost of slower responses to new environments. Due to this trade-off, there is an optimal level of inheritance that maximizes migration speed through changing environments, which enables a diverse population to outperform a nondiverse one. Growing populations might generally leverage the selection-like effects provided by collective behaviors to dynamically shape their own phenotype compositions, without mutations.

Ultrasensitivity and bistability in covalent-modification cycles with positive autoregulation

Switch-like behaviours in biochemical networks are of fundamental significance in biological signal processing, and exist as two distinct types: ultra-sensitivity and bistability. Here we propose two new models of a reversible covalent-modification cycle with positive autoregulation (PAR), a motif structure that is thought to be capable of both ultrasensitivity and bistability in different parameter regimes. These new models appeal to a modelling framework that we call complex-complete, which accounts fully for the molecular complexities of the underlying signalling mechanisms. Each of the two new models encodes a specific molecular mechanism for PAR. We demonstrate that the modelling simplifications for PAR models that have been used in previous work, which rely on Michaelian approximations, are unable to accurately recapitulate the qualitative signalling responses supported by our detailed models. Strikingly, we show that complex-complete PAR models are capable of new qualitative responses such as one-way switches and a 'prozone' effect, depending on the specific PAR-encoding mechanism, which are not supported by Michaelian simplifications. Our results highlight the critical importance of accurately representing the molecular details of biochemical signalling mechanisms, and strongly suggest that the Michaelian approximation is inadequate for predictive models of enzyme-mediated chemical reactions with added regulations such as PAR.

Engineered chemotaxis core signaling units indicate a constrained kinase-off state

Bacterial chemoreceptors, the histidine kinase CheA, and the coupling protein CheW form transmembrane molecular arrays with remarkable sensing properties. The receptors inhibit or stimulate CheA kinase activity depending on the presence of attractants or repellants, respectively. We engineered chemoreceptor cytoplasmic regions to assume a trimer of receptor dimers configuration that formed well-defined complexes with CheA and CheW and promoted a CheA kinase-off state. These mimics of core signaling units were assembled to homogeneity and investigated by site-directed spin-labeling with pulse-dipolar electron-spin resonance spectroscopy (PDS), small-angle x-ray scattering, targeted protein cross-linking, and cryo-electron microscopy. The kinase-off state was especially stable, had relatively low domain mobility, and associated the histidine substrate and docking domains with the kinase core, thus preventing catalytic activity. Together, these data provide an experimentally restrained model for the inhibited state of the core signaling unit and suggest that chemoreceptors indirectly sequester the kinase and substrate domains to limit histidine autophosphorylation.

Sequential modification of bacterial chemoreceptors is key for achieving both accurate adaptation and high gain

Many regulatory and signaling proteins have multiple modification sites. In bacterial chemotaxis, each chemoreceptor has multiple methylation sites that are responsible for adaptation. However, whether the ordering of the multisite methylation process affects adaptation remains unclear. Furthermore, the benefit of having multiple modification sites is also unclear. Here, we show that sequentially ordered methylation/demethylation is critical for perfect adaptation; adaptation accuracy decreases as randomness in the multisite methylation process increases. A tradeoff between adaptation accuracy and response gain is discovered. We find that this accuracy-gain tradeoff is lifted significantly by having more methylation sites, but only when the multisite modification process is sequential. Our study suggests that having multiple modification sites and a sequential modification process constitute a general strategy to achieve both accurate adaptation and high response gain simultaneously. Our theory agrees with existing data and predictions are made to help identify the molecular mechanism underlying ordered covalent modifications.

Bacterial chemoreceptors have multiple methylation sites, but whether the order of methylation matters is unclear. Here, the authors show that sequentially ordered methylation is critical for perfect adaptation and for attenuating the trade-off between accurate adaptation and high response gain.

Accuracy of position determination in Ca^{2+} signaling

A living cell senses its environment and responds to external signals. In this paper, we study theoretically the precision at which cells can determine the position of a spatially localized transient extracellular signal. To this end, we focus on the case where the stimulus is converted into the release of a small molecule that acts as a second messenger, for example, Ca^{2+}, and activates kinases that change the activity of enzymes by phosphorylating them. We analyze the spatial distribution of phosphorylation events using stochastic simulations as well as a mean-field approach. Kinases that need to bind to the cell membrane for getting activated provide more accurate estimates than cytosolic kinases. Our results could explain why the rate of Ca^{2+} detachment from the membrane-binding conventional protein kinase Cα is larger than its phosphorylation rate.

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.

A dual regulation mechanism of histidine kinase CheA identified by combining network-dynamics modeling and system-level input-output data

It is challenging to decipher molecular mechanisms in biological systems from system-level input-output data, especially for complex processes that involve interactions among multiple components. We addressed this general problem for the bacterial histidine kinase CheA, the activity of which is regulated in chemotaxis signaling complexes by bacterial chemoreceptors. We developed a general network model to describe the dynamics of the system, treating the receptor complex with coupling protein CheW and the P3P4P5 domains of kinase CheA as a regulated enzyme with two substrates, ATP and P1, the phosphoryl-accepting domain of CheA. Our simple network model allowed us to search hypothesis space systematically. For different and progressively more complex regulation schemes, we fit our models to a large set of input-output data with the aim of identifying the simplest possible regulation mechanisms consistent with the data. Our modeling and analysis revealed novel dual regulation mechanisms in which receptor activity regulated ATP binding plus one other process, either P1 binding or phosphoryl transfer between P1 and ATP. Strikingly, in our models receptor control affected the kinetic rate constants of substrate association and dissociation equally and thus did not alter the respective equilibrium constants. We suggest experiments that could distinguish between the two dual-regulation mechanisms. This systems-biology approach of combining modeling and a large input-output dataset should be applicable for studying other complex biological processes.

In complex biological systems, it is often difficult to determine which steps in the underlying biochemical network are regulated by the signal by using direct experimental measurements alone. In this paper, we tackled this general problem in the case of the kinase activity of the multi-domain histidine kinase CheA. We developed a quantitative reaction network model to describe the CheA enzyme kinetics by considering all the key reaction steps explicitly. We used this general model with different regulation schemes of progressively increasing complexities to fit a large input-output dataset. Our modeling revealed novel dual regulation mechanisms in which receptor activity regulated two independent reactions in the network including the ATP binding reaction that was previously unsuspected. Through our quantitative analysis, we found that receptors affected the kinetic rate constants of substrate association and dissociation equally and thus did not alter the respective equilibrium constants. Testable predictions of the kinase activity dynamics are made from our models to further distinguish the different dual regulation mechanisms. Our study shows that combining modeling kinetics of the reaction network and input-output data can help reveal the underlying regulation mechanism in complex networks where probing individual reaction is impossible.

Cooperativity Principles in Self-Assembled Nanomedicine

Nanomedicine is a discipline that applies nanoscience and nanotechnology principles to the prevention, diagnosis, and treatment of human diseases. Self-assembly of molecular components is becoming a common strategy in the design and syntheses of nanomaterials for biomedical applications. In both natural and synthetic self-assembled nanostructures, molecular cooperativity is emerging as an important hallmark. In many cases, interplay of many types of noncovalent interactions leads to dynamic nanosystems with emergent properties where the whole is bigger than the sum of the parts. In this review, we provide a comprehensive analysis of the cooperativity principles in multiple self-assembled nanostructures. We discuss the molecular origin and quantitative modeling of cooperative behaviors. In selected systems, we describe the examples on how to leverage molecular cooperativity to design nanomedicine with improved diagnostic precision and therapeutic efficacy in medicine.

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.

Multiple sources of slow activity fluctuations in a bacterial chemosensory network

Cellular networks are intrinsically subject to stochastic fluctuations, but analysis of the resulting noise remained largely limited to gene expression. The pathway controlling chemotaxis of Escherichia coli provides one example where posttranslational signaling noise has been deduced from cellular behavior. This noise was proposed to result from stochasticity in chemoreceptor methylation, and it is believed to enhance environment exploration by bacteria. Here we combined single-cell FRET measurements with analysis based on the fluctuation-dissipation theorem (FDT) to characterize origins of activity fluctuations within the chemotaxis pathway. We observed surprisingly large methylation-independent thermal fluctuations of receptor activity, which contribute to noise comparably to the energy-consuming methylation dynamics. Interactions between clustered receptors involved in amplification of chemotactic signals are also necessary to produce the observed large activity fluctuations. Our work thus shows that the high response sensitivity of this cellular pathway also increases its susceptibility to noise, from thermal and out-of-equilibrium processes.

Feedback between motion and sensation provides nonlinear boost in run-and-tumble navigation

Many organisms navigate gradients by alternating straight motions (runs) with random reorientations (tumbles), transiently suppressing tumbles whenever attractant signal increases. This induces a functional coupling between movement and sensation, since tumbling probability is controlled by the internal state of the organism which, in turn, depends on previous signal levels. Although a negative feedback tends to maintain this internal state close to adapted levels, positive feedback can arise when motion up the gradient reduces tumbling probability, further boosting drift up the gradient. Importantly, such positive feedback can drive large fluctuations in the internal state, complicating analytical approaches. Previous studies focused on what happens when the negative feedback dominates the dynamics. By contrast, we show here that there is a large portion of physiologically-relevant parameter space where the positive feedback can dominate, even when gradients are relatively shallow. We demonstrate how large transients emerge because of non-normal dynamics (non-orthogonal eigenvectors near a stable fixed point) inherent in the positive feedback, and further identify a fundamental nonlinearity that strongly amplifies their effect. Most importantly, this amplification is asymmetric, elongating runs in favorable directions and abbreviating others. The result is a “ratchet-like” gradient climbing behavior with drift speeds that can approach half the maximum run speed of the organism. Our results thus show that the classical drawback of run-and-tumble navigation—wasteful runs in the wrong direction—can be mitigated by exploiting the non-normal dynamics implicit in the run-and-tumble strategy.

Countless bacteria, larvae and even larger organisms (and robots) navigate gradients by alternating periods of straight motion (runs) with random reorientation events (tumbles). Control of the tumble probability is based on previously-encountered signals. A drawback of this run-and-tumble strategy is that occasional runs in the wrong direction are wasteful. Here we show that there is an operating regime within the organism’s internal parameter space where run-and-tumble navigation can be extremely efficient. We characterize how the positive feedback between behavior and sensed signal results in a type of non-equilibrium dynamics, with the organism rapidly tumbling after moving in the wrong direction and extending motion in the right ones. For a distant source, then, the organism can find it fast.

Adaptation through proportion

Adaptation is a ubiquitous feature in biological sensory and signaling networks. It has been suggested that adaptive systems may follow certain simple design principles across diverse organisms, cells and pathways. One class of networks that can achieve adaptation utilizes an incoherent feedforward control, in which two parallel signaling branches exert opposite but proportional effects on the output at steady state. In this paper, we generalize this adaptation mechanism by establishing a steady-state proportionality relationship among a subset of nodes in a network. Adaptation can be achieved by using any two nodes in the sub-network to respectively regulate the output node positively and negatively. We focus on enzyme networks and first identify basic regulation motifs consisting of two and three nodes that can be used to build small networks with proportional relationships. Larger proportional networks can then be constructed modularly similar to LEGOs. Our method provides a general framework to construct and analyze a class of proportional and/or adaptation networks with arbitrary size, flexibility and versatile functional features.

Information processing in bacteria: memory, computation, and statistical physics: a key issues review

Living systems have to constantly sense their external environment and adjust their internal state in order to survive and reproduce. Biological systems, from as complex as the brain to a single E. coli cell, have to process these data in order to make appropriate decisions. How do biological systems sense external signals? How do they process the information? How do they respond to signals? Through years of intense study by biologists, many key molecular players and their interactions have been identified in different biological machineries that carry out these signaling functions. However, an integrated, quantitative understanding of the whole system is still lacking for most cellular signaling pathways, not to say the more complicated neural circuits. To study signaling processes in biology, the key thing to measure is the input-output relationship. The input is the signal itself, such as chemical concentration, external temperature, light (intensity and frequency), and more complex signals such as the face of a cat. The output can be protein conformational changes and covalent modifications (phosphorylation, methylation, etc), gene expression, cell growth and motility, as well as more complex output such as neuron firing patterns and behaviors of higher animals. Due to the inherent noise in biological systems, the measured input-output dependence is often noisy. These noisy data can be analysed by using powerful tools and concepts from information theory such as mutual information, channel capacity, and the maximum entropy hypothesis. This information theory approach has been successfully used to reveal the underlying correlations between key components of biological networks, to set bounds for network performance, and to understand possible network architecture in generating observed correlations. Although the information theory approach provides a general tool in analysing noisy biological data and may be used to suggest possible network architectures in preserving information, it does not reveal the underlying mechanism that leads to the observed input-output relationship, nor does it tell us much about which information is important for the organism and how biological systems use information to carry out specific functions. To do that, we need to develop models of the biological machineries, e.g. biochemical networks and neural networks, to understand the dynamics of biological information processes. This is a much more difficult task. It requires deep knowledge of the underlying biological network-the main players (nodes) and their interactions (links)-in sufficient detail to build a model with predictive power, as well as quantitative input-output measurements of the system under different perturbations (both genetic variations and different external conditions) to test the model predictions to guide further development of the model. Due to the recent growth of biological knowledge thanks in part to high throughput methods (sequencing, gene expression microarray, etc) and development of quantitative in vivo techniques such as various florescence technology, these requirements are starting to be realized in different biological systems. The possible close interaction between quantitative experimentation and theoretical modeling has made systems biology an attractive field for physicists interested in quantitative biology. In this review, we describe some of the recent work in developing a quantitative predictive model of bacterial chemotaxis, which can be considered as the hydrogen atom of systems biology. Using statistical physics approaches, such as the Ising model and Langevin equation, we study how bacteria, such as E. coli, sense and amplify external signals, how they keep a working memory of the stimuli, and how they use these data to compute the chemical gradient. In particular, we will describe how E. coli cells avoid cross-talk in a heterogeneous receptor cluster to keep a ligand-specific memory. We will also study the thermodynamic costs of adaptation for cells to maintain an accurate memory. The statistical physics based approach described here should be useful in understanding design principles for cellular biochemical circuits in general.

Theory in Biology: Figure 1 or Figure 7?

The pace of modern science is staggering. The quantities of data now flowing from DNA sequencers, fluorescence and electron microscopes, mass spectrometers, and other mind-blowing instruments leave us faced with information overload. This explosion in data has brought on its heels a concomitant need for efforts at the kinds of synthesis and unification we see in theoretical physics. Often in cell biology, when theoretical modeling takes place, it is as a figure 7 reflection on experiments that have already been done, with data fitting providing a metric of success. Figure 1 theory, by way of contrast, is about living dangerously by turning our thinking into formal mathematical predictions and confronting that math with experiments that have not yet been done.

Notes on stochastic (bio)-logic gates: computing with allosteric cooperativity

Recent experimental breakthroughs have finally allowed to implement in-vitro reaction kinetics (the so called enzyme based logic) which code for two-inputs logic gates and mimic the stochastic AND (and NAND) as well as the stochastic OR (and NOR). This accomplishment, together with the already-known single-input gates (performing as YES and NOT), provides a logic base and paves the way to the development of powerful biotechnological devices. However, as biochemical systems are always affected by the presence of noise (e.g. thermal), standard logic is not the correct theoretical reference framework, rather we show that statistical mechanics can work for this scope: here we formulate a complete statistical mechanical description of the Monod-Wyman-Changeaux allosteric model for both single and double ligand systems, with the purpose of exploring their practical capabilities to express noisy logical operators and/or perform stochastic logical operations. Mixing statistical mechanics with logics, and testing quantitatively the resulting findings on the available biochemical data, we successfully revise the concept of cooperativity (and anti-cooperativity) for allosteric systems, with particular emphasis on its computational capabilities, the related ranges and scaling of the involved parameters and its differences with classical cooperativity (and anti-cooperativity).

Properties of cooperatively induced phases in sensing models

A large number of eukaryotic cells are able to directly detect external chemical gradients with great accuracy and the ultimate limit to their sensitivity has been a topic of debate for many years. Previous work has been done to understand many aspects of this process but little attention has been paid to the possibility of emergent sensing states. Here we examine how cooperation between sensors existing in a two-dimensional network, as they do on the cell's surface, can both enhance and fundamentally alter the response of the cell to a spatially varying signal. We show that weakly interacting sensors linearly amplify the cell's response to an external gradient while a network of strongly interacting sensors form a collective nonlinear response with two separate domains of active and inactive sensors forming what have called a "1/2-state." In our analysis we examine the cell's ability to sense the direction of a signal and pay special attention to the substantially different behavior realized in the strongly interacting regime.

Relation between chemotaxis and consumption of amino acids in bacteria

Chemotaxis enables bacteria to navigate chemical gradients in their environment, accumulating toward high concentrations of attractants and avoiding high concentrations of repellents. Although finding nutrients is likely to be an important function of bacterial chemotaxis, not all characterized attractants are nutrients. Moreover, even for potential nutrients, the exact relation between the metabolic value of chemicals and their efficiency as chemoattractants has not been systematically explored. Here we compare the chemotactic response of amino acids with their use by bacteria for two well‐established models of chemotactic behavior, Escherichia coli and Bacillus subtilis. We demonstrate that in E. coli chemotaxis toward amino acids indeed strongly correlates with their utilization. However, no such correlation is observed for B. subtilis, suggesting that in this case, the amino acids are not followed because of their nutritional value but rather as environmental cues.

Napoleon Is in Equilibrium

It has been said that the cell is the test tube of the twenty-first century. If so, the theoretical tools needed to quantitatively and predictively describe what goes on in such test tubes lag sorely behind the stunning experimental advances in biology seen in the decades since the molecular biology revolution began. Perhaps surprisingly, one of the theoretical tools that has been used with great success on problems ranging from how cells communicate with their environment and each other to the nature of the organization of proteins and lipids within the cell membrane is statistical mechanics. A knee-jerk reaction to the use of statistical mechanics in the description of cellular processes is that living organisms are so far from equilibrium that one has no business even thinking about it. But such reactions are probably too hasty given that there are many regimes in which, because of a separation of timescales, for example, such an approach can be a useful first step. In this article, we explore the power of statistical mechanical thinking in the biological setting, with special emphasis on cell signaling and regulation. We show how such models are used to make predictions and describe some recent experiments designed to test them. We also consider the limits of such models based on the relative timescales of the processes of interest.

The role of membrane-mediated interactions in the assembly and architecture of chemoreceptor lattices

In vivo fluorescence microscopy and electron cryo-tomography have revealed that chemoreceptors self-assemble into extended honeycomb lattices of chemoreceptor trimers with a well-defined relative orientation of trimers. The signaling response of the observed chemoreceptor lattices is remarkable for its extreme sensitivity, which relies crucially on cooperative interactions among chemoreceptor trimers. In common with other membrane proteins, chemoreceptor trimers are expected to deform the surrounding lipid bilayer, inducing membrane-mediated anisotropic interactions between neighboring trimers. Here we introduce a biophysical model of bilayer-chemoreceptor interactions, which allows us to quantify the role of membrane-mediated interactions in the assembly and architecture of chemoreceptor lattices. We find that, even in the absence of direct protein-protein interactions, membrane-mediated interactions can yield assembly of chemoreceptor lattices at very dilute trimer concentrations. The model correctly predicts the observed honeycomb architecture of chemoreceptor lattices as well as the observed relative orientation of chemoreceptor trimers, suggests a series of “gateway” states for chemoreceptor lattice assembly, and provides a simple mechanism for the localization of large chemoreceptor lattices to the cell poles. Our model of bilayer-chemoreceptor interactions also helps to explain the observed dependence of chemotactic signaling on lipid bilayer properties. Finally, we consider the possibility that membrane-mediated interactions might contribute to cooperativity among neighboring chemoreceptor trimers.

The chemotaxis system allows bacteria to respond to minute changes in chemical concentration, and serves as a paradigm for biological signal processing and the self-assembly of large protein lattices in living cells. The sensitivity of the chemotaxis system relies crucially on cooperative interactions among chemoreceptor trimers, which are organized into intricate honeycomb lattices. Chemoreceptors are membrane proteins and, hence, are expected to deform the surrounding lipid bilayer, leading to membrane-mediated interactions between chemoreceptor trimers. Using a biophysical model of bilayer-chemoreceptor interactions we show that the membrane-mediated interactions induced by chemoreceptor trimers provide a mechanism for the observed self-assembly of chemoreceptor lattices. We find that the directionality of membrane-mediated interactions between trimers complements protein-protein interactions in the stabilization of the observed honeycomb architecture of chemoreceptor lattices. Our results suggest that the symmetry of membrane protein complexes such as chemoreceptor trimers is reflected in the anisotropy of membrane-mediated interactions, yielding a general mechanism for the self-assembly of ordered protein lattices in cell membranes.

Precision sensing by two opposing gradient sensors: how does Escherichia coli find its preferred pH level?

It is essential for bacteria to find optimal conditions for their growth and survival. The optimal levels of certain environmental factors (such as pH and temperature) often correspond to some intermediate points of the respective gradients. This requires the ability of bacteria to navigate from both directions toward the optimum location and is distinct from the conventional unidirectional chemotactic strategy. Remarkably, Escherichia coli cells can perform such a precision sensing task in pH taxis by using the same chemotaxis machinery, but with opposite pH responses from two different chemoreceptors (Tar and Tsr). To understand bacterial pH sensing, we developed an Ising-type model for a mixed cluster of opposing receptors based on the push-pull mechanism. Our model can quantitatively explain experimental observations in pH taxis for various mutants and wild-type cells. We show how the preferred pH level depends on the relative abundance of the competing sensors and how the sensory activity regulates the behavioral response. Our model allows us to make quantitative predictions on signal integration of pH and chemoattractant stimuli. Our study reveals two general conditions and a robust push-pull scheme for precision sensing, which should be applicable in other adaptive sensory systems with opposing gradient sensors.

Excitation and adaptation in bacteria-a model signal transduction system that controls taxis and spatial pattern formation

The machinery for transduction of chemotactic stimuli in the bacterium E. coli is one of the most completely characterized signal transduction systems, and because of its relative simplicity, quantitative analysis of this system is possible. Here we discuss models which reproduce many of the important behaviors of the system. The important characteristics of the signal transduction system are excitation and adaptation, and the latter implies that the transduction system can function as a "derivative sensor" with respect to the ligand concentration in that the DC component of a signal is ultimately ignored if it is not too large. This temporal sensing mechanism provides the bacterium with a memory of its passage through spatially- or temporally-varying signal fields, and adaptation is essential for successful chemotaxis. We also discuss some of the spatial patterns observed in populations and indicate how cell-level behavior can be embedded in population-level descriptions.

Prolonged stimuli alter the bacterial chemosensory clusters

The clustering of membrane-bound receptors plays an essential role in various biological systems. A notable model system for studying this phenomenon is the bacterial chemosensory cluster that allows motile bacteria to navigate along chemical gradients in their environment. While the basic structure of these chemosensory clusters is becoming clear, their dynamic nature and operation are not yet understood. By measuring the fluorescence polarization of tagged receptor clusters in live Escherichia coli cells, we provide evidence for stimulus-induced dynamics in these sensory clusters. We find that when a stimulus is applied, the packing of the receptors slowly decreases and that the process reverses when the stimulus is removed. Consistent with these physical changes we find that the effective cooperativity of the kinase response slowly evolves in the presence of a stimulus. Time-lapse fluorescence imaging indicates that, despite these changes, the receptor clusters do not generally dissociate upon ligand binding. These data reveal stimulus-dependent plasticity in chemoreceptor clusters.

Coordinated switching of bacterial flagellar motors: evidence for direct motor-motor coupling?

The swimming of Escherichia coli is powered by its multiple flagellar motors. Each motor spins either clockwise or counterclockwise, under the control of an intracellular regulator, CheY-P. There can be two mechanisms (extrinsic and intrinsic) to coordinate the switching of bacterial motors. The extrinsic one arises from the fact that different motors in the same cell sense a common input (CheY-P) which fluctuates near the motors' response threshold. An alternative, intrinsic mechanism is direct motor-motor coupling which makes synchronized switching energetically favorable. Here, we develop simple models for both mechanisms and uncover their different hallmarks. A quantitative comparison to the recent experiments suggests that the direct coupling mechanism may be accountable for the observed sharp correlation between motors in a single Escherichia coli. Possible origins of this coupling (e.g., hydrodynamic interaction) are discussed.

Statistical mechanics of Monod-Wyman-Changeux (MWC) models

The 50th anniversary of the classic Monod-Wyman-Changeux (MWC) model provides an opportunity to survey the broader conceptual and quantitative implications of this quintessential biophysical model. With the use of statistical mechanics, the mathematical implementation of the MWC concept links problems that seem otherwise to have no ostensible biological connection including ligand-receptor binding, ligand-gated ion channels, chemotaxis, chromatin structure and gene regulation. Hence, a thorough mathematical analysis of the MWC model can illuminate the performance limits of a number of unrelated biological systems in one stroke. The goal of our review is twofold. First, we describe in detail the general physical principles that are used to derive the activity of MWC molecules as a function of their regulatory ligands. Second, we illustrate the power of ideas from information theory and dynamical systems for quantifying how well the output of MWC molecules tracks their sensory input, giving a sense of the "design" constraints faced by these receptors.

Quantitative modeling of bacterial chemotaxis: signal amplification and accurate adaptation

We review the recent developments in understanding the bacterial chemotaxis signaling pathway by using quantitative modeling methods. The models developed are based on structural information of the signaling complex and the dynamics of the underlying biochemical network. We focus on two important functions of the bacterial chemotaxis signaling pathway: signal amplification and adaptation. We describe in detail the structure and the dynamics of the mathematical models and how they compare with existing experiments, emphasizing the predictability of the models. Finally, we outline future directions for developing the modeling approach to better understand the bacterial chemosensory system.

A "trimer of dimers"-based model for the chemotactic signal transduction network in bacterial chemotaxis

The network that controls chemotaxis in Escherichia coli is one of the most completely characterized signal transduction systems to date. Receptor clustering accounts for characteristics such as high sensitivity, precise adaptation over a wide dynamic range of ligand concentrations, and robustness to variations in the amounts of intracellular proteins. To gain insights into the structure-function relationship of receptor clusters and understand the mechanism behind the high-performance signaling, we develop and analyze a model for a single trimer of dimers. This new model extends an earlier model (Spiro et al. in Proc. Natl. Acad. Sci. 94:7263-7268, 1997) to incorporate the recent experimental findings that the core structure of receptor clusters is the trimer of receptor dimers. We show that the model can reproduce most of the experimentally-observed behaviors, including excitation, adaptation, high sensitivity, and robustness to parameter variations. In addition, the model makes a number of new predictions as to how the adaptation time varies with the expression level of various proteins involved in signal transduction. Our results provide a more mechanistically-based description of the structure-function relationship for the signaling system, and show the key role of the interaction among dimer members of the trimer in the chemotactic response of cells.

Pathway-based mean-field model for Escherichia coli chemotaxis

We develop a mean-field theory for Escherichia coli chemotaxis based on the coupled spatiotemporal dynamics of the cell population and the mean receptor methylation level field. This multiscale model connects the cells' population level motility behavior with the molecular level pathway dynamics. It reveals a simple scaling dependence of the chemotaxis velocity on the adaptation rate in exponential gradients. It explains the molecular origin of a maximum chemotaxis velocity. Simulations of our model in various spatiotemporal stimuli profiles show quantitative agreements with experiments. Moreover, it predicts a surprising reversal of chemotaxis group velocity in traveling wave environments. Our approach may be used to bridge molecular level pathway dynamics with cellular behavior in other biological systems.

The physics of protein-DNA interaction networks in the control of gene expression

Protein-DNA interaction networks play a central role in many fundamental cellular processes. In gene regulation, physical interactions and reactions among the molecular components together with the physical properties of DNA control how genes are turned on and off. A key player in all these processes is the inherent flexibility of DNA, which provides an avenue for long-range interactions between distal DNA elements through DNA looping. Such versatility enables multiple interactions and results in additional complexity that is remarkably difficult to address with traditional approaches. This topical review considers recent advances in statistical physics methods to study the assembly of protein-DNA complexes with loops, their effects in the control of gene expression, and their explicit application to the prototypical lac operon genetic system of the E. coli bacterium. In the last decade, it has been shown that the underlying physical properties of DNA looping can actively control transcriptional noise, cell-to-cell variability, and other properties of gene regulation, including the balance between robustness and sensitivity of the induction process. These physical properties are largely dependent on the free energy of DNA looping, which accounts for DNA bending and twisting effects. These new physical methods have also been used in reverse to uncover the actual in vivo free energy of looping double-stranded DNA in living cells, which was not possible with existing experimental techniques. The results obtained for DNA looping by the lac repressor inside the E. coli bacterium showed a more malleable DNA than expected as a result of the interplay of the simultaneous presence of two distinct conformations of looped DNA.

Responding to chemical gradients: bacterial chemotaxis

Chemotaxis allows bacteria to follow gradients of nutrients and other environmental stimuli. The bacterium Escherichia coli performs chemotaxis via a run-and-tumble strategy in which sensitive temporal comparisons lead to a biased random walk, with longer runs in the preferred gradient direction. The chemotaxis network of E. coli has developed over the years into one of the most thoroughly studied model systems for signal transduction and behavior, yielding general insights into such properties of cellular networks as signal amplification, signal integration, and robustness. Despite its relative simplicity, the operation of the E. coli chemotaxis network is highly refined and evolutionarily optimized at many levels. For example, recent studies revealed that the network adjusts its signaling properties dependent on the extracellular environment, apparently to optimize chemotaxis under particular conditions. The network can even utilize potentially detrimental stochastic fluctuations in protein levels and reaction rates to maximize the chemotactic performance of the population.

Adapt locally and act globally: strategy to maintain high chemoreceptor sensitivity in complex environments

In bacterial chemotaxis, several types of ligand-specific receptors form mixed clusters, wherein receptor-receptor interactions lead to signal amplification and integration. However, it remains unclear how a mixed receptor cluster adapts to individual stimuli and whether it can differentiate between different types of ligands. Here, we combine theoretical modeling with experiments to reveal the adaptation dynamics of the mixed chemoreceptor cluster in Escherichia coli. We show that adaptation occurs locally and is ligand-specific: only the receptor that binds the external ligand changes its methylation level when the system adapts, whereas other types of receptors change methylation levels transiently. Permanent methylation crosstalk occurs when the system fails to adapt accurately. This local adaptation mechanism enables cells to differentiate individual stimuli by encoding them into the methylation levels of corresponding types of chemoreceptors. It tunes each receptor to its most responsive state to maintain high sensitivity in complex environments and prevents saturation of the cluster by one signal.

Lateral density of receptor arrays in the membrane plane influences sensitivity of the E. coli chemotaxis response

In chemotactic bacteria, transmembrane chemoreceptors, CheA and CheW form the core signalling complex of the chemotaxis sensory apparatus. These complexes are organized in extended arrays in the cytoplasmic membrane that allow bacteria to respond to changes in concentration of extracellular ligands via a cooperative, allosteric response that leads to substantial amplification of the signal induced by ligand binding. Here, we have combined cryo-electron tomographic studies of the 3D spatial architecture of chemoreceptor arrays in intact E. coli cells with computational modelling to develop a predictive model for the cooperativity and sensitivity of the chemotaxis response. The predictions were tested experimentally using fluorescence resonance energy transfer (FRET) microscopy. Our results demonstrate that changes in lateral packing densities of the partially ordered, spatially extended chemoreceptor arrays can modulate the bacterial chemotaxis response, and that information about the molecular organization of the arrays derived by cryo-electron tomography of intact cells can be translated into testable, predictive computational models of the chemotaxis response.

Optimal receptor-cluster size determined by intrinsic and extrinsic noise

Biological cells sense external chemical stimuli in their environment using cell-surface receptors. To increase the sensitivity of sensing, receptors often cluster. This process occurs most noticeably in bacterial chemotaxis, a paradigm for sensing and signaling in general. While amplification of weak stimuli is useful in the absence of noise, its usefulness is less clear in the presence of extrinsic input noise and intrinsic signaling noise. Here, exemplified in a bacterial chemotaxis system, we combine the allosteric Monod-Wyman-Changeux model for signal amplification by receptor complexes with calculations of noise to study their interconnectedness. Importantly, we calculate the signal-to-noise ratio, describing the balance of beneficial and detrimental effects of clustering for the cell. Interestingly, we find that there is no advantage for the cell to build receptor complexes for noisy input stimuli in the absence of intrinsic signaling noise. However, with intrinsic noise, an optimal complex size arises in line with estimates of the size of chemoreceptor complexes in bacteria and protein aggregates in lipid rafts of eukaryotic cells.

A modular gradient-sensing network for chemotaxis in Escherichia coli revealed by responses to time-varying stimuli

Combining in vivo FRET with time-varying stimuli, such as steps, ramps, and sinusoids allowed deduction of the molecular mechanisms underlying cellular signal processing.

The bacterial chemotaxis pathway can be described as a two-module feedback circuit, the transfer functions of which we have characterized quantitatively by experiment. Model-driven experimental design allowed the use of a single FRET pair for measurements of both transfer functions of the pathway.

The adaptation module's transfer function revealed that feedback near steady state is weak, consistent with high sensitivity to shallow gradients, but also strong steady-state fluctuations in pathway output.

The measured response to oscillatory stimuli defines the frequency band over which the chemotaxis system can compute time derivatives.

In searching for better environments, bacteria sample their surroundings by random motility, and make temporal comparisons of experienced sensory cues to bias their movement toward favorable directions (Berg and Brown, 1972). Thus, the problem of sensing spatial gradients is reduced to time-derivative computations, carried out by a signaling pathway that is well characterized at the molecular level in Escherichia coli. Here, we study the physiology of this signal processing system in vivo by fluorescence resonance energy transfer (FRET) experiments in which live cells are stimulated by time-varying chemoeffector signals. By measuring FRET between the active response regulator of the pathway CheY-P and its phosphatase CheZ, each labeled with GFP variants, we obtain a readout that is directly proportional to pathway activity (Sourjik et al, 2007). We analyze the measured response functions in terms of mechanistic models of signaling, and discuss functional consequences of the observed quantitative characteristics.

Experiments are guided by a coarse-grained modular model (Tu et al, 2008) of the sensory network (Figure 1), in which we identify two important ‘transfer functions': one corresponding to the receptor–kinase complex, which responds to changes in input ligand concentration on a fast time scale, and another corresponding to the adaptation system, which provides negative feedback, opposing the effect of ligand on a slower time scale. For the receptor module, we calibrate an allosteric MWC-type model of the receptor–kinase complex by FRET measurements of the ‘open-loop' transfer function G([L],m) using step stimuli. This calibration provides a basis for using the same FRET readout (between CheY-P and CheZ) to further study properties of the adaptation module.

It is well known that adaptation in E. coli's chemotaxis system uses integral feedback, which guarantees exact restoration of the baseline activity after transient responses to step stimuli (Barkai and Leibler, 1997; Yi et al, 2000). However, the output of time-derivative computations during smoothly varying stimuli depends not only on the presence of integral feedback, but also on what is being integrated. As this integrand can in general be any function of the output, we represent it by a black-box function F(a) in our model, and set out to determine its shape by experiments with time-varying stimuli.

We first apply exponential ramp stimuli—waveforms in which the logarithm of the stimulus level varies linearly with time, at a fixed rate r. It was shown many years ago that during such a stimulus, the kinase output of the pathway changes to a new constant value, a_c that is dependent on the applied ramp rate, r (Block et al, 1983). A plot of a_c versus r (Figure 5A) can thus be considered as an output of time-derivative computations by the network, and could also be used to study the ‘gradient sensitivity' of bacteria traveling at constant speeds.

To obtain the feedback transfer function, F(a), we apply a simple coordinate transformation, identified using our model, to the same ramp-response data (Figure 5B). This function reveals how the temporal rate of change of the feedback signal m depends on the current output signal a. The shape of this function is analyzed using a biochemical reaction scheme, from which in vivo kinetic parameters of the feedback enzymes, CheR and CheB, are extracted. The fitted Michaelis constants for these enzymatic reactions are small compared with the steady-state abundance of their substrates, thus indicating that these enzymes operate close to saturation in vivo. The slope of the function near steady state can be used to assess the strength of feedback, and to compute the relaxation time of the system, τ_m. Relaxation is found to be slow (i.e. large τ_m), consistent with large fluctuations about the steady-state activity caused by the near-saturation kinetics of the feedback enzymes (Emonet and Cluzel, 2008).

Finally, exponential sine-wave stimuli are used to map out the system's frequency response (Figure 5C). The measured data points for both the amplitude and phase of the response are found to be in excellent agreement with model predictions based on parameters from the independently measured step and ramp responses. No curve fitting was required to obtain this agreement. Although the amplitude response as a function of frequency resembles a first-order high-pass filter with a well-defined cutoff frequency, ν_m, we point out that the chemotaxis pathway is actually a low-pass filter if the time derivative of the input is viewed as the input signal. In this latter perspective, ν_m defines an upper bound for the frequency band over which time-derivative computations can be carried out.

The two types of measurements yield complementary information regarding time-derivative computations by E. coli. The ramp-responses characterize the asymptotically constant output when a temporal gradient is held fixed over extended periods. Interestingly, the ramp responses do not depend on receptor cooperativity, but only on properties of the adaptation system, and thus can be used to reveal the in vivo adaptation kinetics, even outside the linear regime of the kinase response. The frequency response is highly relevant in considering spatial searches in the real world, in which experienced gradients are not held fixed in time. The characteristic cutoff frequency ν_m is found by working within the linear regime of the kinase response, and depends on parameters from both modules (it increases with both cooperativity in the receptor module, and the strength of feedback in the adaptation module).

Both ramp responses and sine-wave responses were measured at two different temperatures (22 and 32°C), and found to differ significantly. Both the slope of F(a) near steady state, from ramp experiments, and the characteristic cutoff frequency, from sine-wave experiments, were higher by a factor of ∼3 at 32°C. Fits of the enzymatic model to F(a) suggest that temperature affects the maximal velocity (V_max) more strongly than the Michaelis constants (K_m) for CheR and CheB.

Successful application of inter-molecular FRET in live cells using GFP variants always requires some degree of serendipity. Genetic fusions to these bulky fluorophores can impair the function of the original proteins, and even when fusions are functional, efficient FRET still requires the fused fluorophores to come within the small (<10 nm) Förster radius on interactions between the labeled proteins. Thus, when a successful FRET pair is identified, it is desirable to make the most of it. We have shown here that combined with careful temporal control of input stimuli, and appropriately calibrated models, a single FRET pair can be used to study the structure of multiple transfer functions within a signaling network.

The Escherichia coli chemotaxis-signaling pathway computes time derivatives of chemoeffector concentrations. This network features modules for signal reception/amplification and robust adaptation, with sensing of chemoeffector gradients determined by the way in which these modules are coupled in vivo. We characterized these modules and their coupling by using fluorescence resonance energy transfer to measure intracellular responses to time-varying stimuli. Receptor sensitivity was characterized by step stimuli, the gradient sensitivity by exponential ramp stimuli, and the frequency response by exponential sine-wave stimuli. Analysis of these data revealed the structure of the feedback transfer function linking the amplification and adaptation modules. Feedback near steady state was found to be weak, consistent with strong fluctuations and slow recovery from small perturbations. Gradient sensitivity and frequency response both depended strongly on temperature. We found that time derivatives can be computed by the chemotaxis system for input frequencies below 0.006 Hz at 22°C and below 0.018 Hz at 32°C. Our results show how dynamic input–output measurements, time honored in physiology, can serve as powerful tools in deciphering cell-signaling mechanisms.

Physical limits on cellular sensing of spatial gradients

Many eukaryotic cells are able to detect chemical gradients by directly measuring spatial concentration differences. The precision of such gradient sensing is limited by fluctuations in the binding of diffusing particles to specific receptors on the cell surface. Here, we explore the physical limits of the spatial sensing mechanism by modeling the chemotactic cell as an Ising spin chain subject to a spatially varying field. Our results demonstrate that the accuracy to sense the gradient direction not only increases dramatically with the cell size but also can be improved significantly by introducing receptor cooperativity. Thus, receptor coupling may open the possibility for small bacteria to perform spatial measurements of gradients, as supported by a recent experimental finding.

Chemotactic response and adaptation dynamics in Escherichia coli

Adaptation of the chemotaxis sensory pathway of the bacterium Escherichia coli is integral for detecting chemicals over a wide range of background concentrations, ultimately allowing cells to swim towards sources of attractant and away from repellents. Its biochemical mechanism based on methylation and demethylation of chemoreceptors has long been known. Despite the importance of adaptation for cell memory and behavior, the dynamics of adaptation are difficult to reconcile with current models of precise adaptation. Here, we follow time courses of signaling in response to concentration step changes of attractant using in vivo fluorescence resonance energy transfer measurements. Specifically, we use a condensed representation of adaptation time courses for efficient evaluation of different adaptation models. To quantitatively explain the data, we finally develop a dynamic model for signaling and adaptation based on the attractant flow in the experiment, signaling by cooperative receptor complexes, and multiple layers of feedback regulation for adaptation. We experimentally confirm the predicted effects of changing the enzyme-expression level and bypassing the negative feedback for demethylation. Our data analysis suggests significant imprecision in adaptation for large additions. Furthermore, our model predicts highly regulated, ultrafast adaptation in response to removal of attractant, which may be useful for fast reorientation of the cell and noise reduction in adaptation.

Bacterial chemotaxis is a paradigm for sensory systems, and thus has attracted immense interest from biologists and modelers alike. Using this pathway, cells can sense chemical molecules in their environment, and bias their movement towards nutrients and away from toxins. To avoid over- or understimulation of the signaling pathway, receptors adapt to current external conditions by covalent receptor modification, ultimately allowing cells to chemotax over a wide range of background concentrations. While the robustness and precision in adaptation was previously explained, we quantify the dynamics of adaptation, important for cell memory and behavior, as well as noise filtering in the pathway. Specifically, we study the intracellular signaling response and subsequent adaptation to concentration step changes in attractant chemicals. We combine measurements of signaling in living cells with a dynamic model for strongly coupled receptors, even including the effects of concentration flow in the experiment. Using a novel way of summarizing time-dependent data, we derive a new adaptation model, predicting additional layers of feedback regulation. As a consequence, adaptation to sudden exposure of unfavorable conditions is very fast, which may be useful for a quick reorientation and escape of the cell.

An agent-based model of signal transduction in bacterial chemotaxis

We report the application of agent-based modeling to examine the signal transduction network and receptor arrays for chemotaxis in Escherichia coli, which are responsible for regulating swimming behavior in response to environmental stimuli. Agent-based modeling is a stochastic and bottom-up approach, where individual components of the modeled system are explicitly represented, and bulk properties emerge from their movement and interactions. We present the Chemoscape model: a collection of agents representing both fixed membrane-embedded and mobile cytoplasmic proteins, each governed by a set of rules representing knowledge or hypotheses about their function. When the agents were placed in a simulated cellular space and then allowed to move and interact stochastically, the model exhibited many properties similar to the biological system including adaptation, high signal gain, and wide dynamic range. We found the agent based modeling approach to be both powerful and intuitive for testing hypotheses about biological properties such as self-assembly, the non-linear dynamics that occur through cooperative protein interactions, and non-uniform distributions of proteins in the cell. We applied the model to explore the role of receptor type, geometry and cooperativity in the signal gain and dynamic range of the chemotactic response to environmental stimuli. The model provided substantial qualitative evidence that the dynamic range of chemotactic response can be traced to both the heterogeneity of receptor types present, and the modulation of their cooperativity by their methylation state.

Phototactic and chemotactic signal transduction by transmembrane receptors and transducers in microorganisms

Microorganisms show attractant and repellent responses to survive in the various environments in which they live. Those phototaxic (to light) and chemotaxic (to chemicals) responses are regulated by membrane-embedded receptors and transducers. This article reviews the following: (1) the signal relay mechanisms by two photoreceptors, Sensory Rhodopsin I (SRI) and Sensory Rhodopsin II (SRII) and their transducers (HtrI and HtrII) responsible for phototaxis in microorganisms; and (2) the signal relay mechanism of a chemoreceptor/transducer protein, Tar, responsible for chemotaxis in E. coli. Based on results mainly obtained by our group together with other findings, the possible molecular mechanisms for phototaxis and chemotaxis are discussed.

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 v(d) is a constant in an exponential attractant concentration gradient [L] proportional, variantexp(Gx). v(d) depends linearly on the exponential gradient G before it saturates when G is larger than a critical value G(C). We find that G(C) is determined by the intracellular adaptation rate k(R) with a simple scaling law: G(C) infinity k(1/2)(R). The linear dependence of v(d) 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 G(C) and its scaling with k(R) 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 omega of the time-varying spatial gradient becomes faster than a critical value omega(c), which is controlled by the cell's chemotaxis adaptation rate k(R). 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.

Responses of Escherichia coli bacteria to two opposing chemoattractant gradients depend on the chemoreceptor ratio

Escherichia coli chemotaxis has long served as a simple model of environmental signal processing, and bacterial responses to single chemical gradients are relatively well understood. Less is known about the chemotactic behavior of E. coli in multiple chemical gradients. In their native environment, cells are often exposed to multiple chemical stimuli. Using a recently developed microfluidic chemotaxis device, we exposed E. coli cells to two opposing but equally potent gradients of major attractants, methyl-aspartate and serine. The responses of E. coli cells demonstrated that chemotactic decisions depended on the ratio of the respective receptor number of Tar/Tsr. In addition, the ratio of Tar to Tsr was found to vary with cells' growth conditions, whereby it depended on the culture density but not on the growth duration. These results provide biological insights into the decision-making processes of chemotactic bacteria that are subjected to multiple chemical stimuli and demonstrate the importance of the cellular microenvironment in determining phenotypic behavior.

Chemotaxis: how bacteria use memory

Bacterial chemotaxis represents one of the simplest and best studied examples of unicellular behavior. Chemotaxis allows swimming bacterial cells to follow chemical gradients in the environment by performing temporal comparisons of ligand concentrations. The process of chemotaxis in the model bacteriumEscherichia colihas been studied in great molecular detail over the past 40 years, using a large range of experimental tools to investigate physiology, genetics and biochemistry of the system. The abundance of quantitative experimental data enabled detailed computational modeling of the pathway and theoretical analyses of such properties as robustness and signal amplification. Because of the temporal mode of gradient sensing in bacterial chemotaxis, molecular memory is an essential component of the chemotaxis pathway. Recent studies suggest that the memory time scale has been evolutionary optimized to perform optimal comparisons of stimuli while swimming in the gradient. Moreover, noise in the adaptation system, which results from variations of the adaptation rate both over time and among cells, might be beneficial for the overall chemotactic performance of the population.

Defining network topologies that can achieve biochemical adaptation

Many signaling systems show adaptation-the ability to reset themselves after responding to a stimulus. We computationally searched all possible three-node enzyme network topologies to identify those that could perform adaptation. Only two major core topologies emerge as robust solutions: a negative feedback loop with a buffering node and an incoherent feedforward loop with a proportioner node. Minimal circuits containing these topologies are, within proper regions of parameter space, sufficient to achieve adaptation. More complex circuits that robustly perform adaptation all contain at least one of these topologies at their core. This analysis yields a design table highlighting a finite set of adaptive circuits. Despite the diversity of possible biochemical networks, it may be common to find that only a finite set of core topologies can execute a particular function. These design rules provide a framework for functionally classifying complex natural networks and a manual for engineering networks. For a video summary of this article, see the PaperFlick file with the Supplemental Data available online.

Signal amplification in a lattice of coupled protein kinases

The bacterium Escherichia coli detects chemical attractants and repellents by means of a cluster of transmembrane receptors and associated molecules. Experiments have shown that this cluster amplifies the signal about 35-fold and current models attribute this amplification to cooperative interactions between neighbouring receptors. However, when applied to the mixed population of receptors of wild-type E. coli, these models lead to indiscriminate methylation of all receptor types rather than the selective methylation observed experimentally. In this paper, we propose that cooperative interactions occur not between receptors but in the underlying lattice of CheA molecules. In our model, each CheA molecule is stimulated by its neighbours via their flexible P1 domains and modulated by the ligand binding and methylation states of associated receptors. We test this idea with detailed, molecular-based stochastic simulations and show that it gives an accurate reproduction of signalling in this system, including ligand-specific adaptation.

Logarithmic sensing in Escherichia coli bacterial chemotaxis

We studied the response of swimming Escherichia coli (E. coli) bacteria in a comprehensive set of well-controlled chemical concentration gradients using a newly developed microfluidic device and cell tracking imaging technique. In parallel, we carried out a multi-scale theoretical modeling of bacterial chemotaxis taking into account the relevant internal signaling pathway dynamics, and predicted bacterial chemotactic responses at the cellular level. By measuring the E. coli cell density profiles across the microfluidic channel at various spatial gradients of ligand concentration grad[L] and the average ligand concentration [L] near the peak chemotactic response region, we demonstrated unambiguously in both experiments and model simulation that the mean chemotactic drift velocity of E. coli cells increased monotonically with grad [L]/[L] or approximately grad(log[L])--that is E. coli cells sense the spatial gradient of the logarithmic ligand concentration. The exact range of the log-sensing regime was determined. The agreements between the experiments and the multi-scale model simulation verify the validity of the theoretical model, and revealed that the key microscopic mechanism for logarithmic sensing in bacterial chemotaxis is the adaptation kinetics, in contrast to explanations based directly on ligand occupancy.

Differential affinity and catalytic activity of CheZ in E. coli chemotaxis

Push–pull networks, in which two antagonistic enzymes control the activity of a messenger protein, are ubiquitous in signal transduction pathways. A classical example is the chemotaxis system of the bacterium Escherichia coli, in which the kinase CheA and the phosphatase CheZ regulate the phosphorylation level of the messenger protein CheY. Recent experiments suggest that both the kinase and the phosphatase are localized at the receptor cluster, and Vaknin and Berg recently demonstrated that the spatial distribution of the phosphatase can markedly affect the dose–response curves. We argue, using mathematical modeling, that the canonical model of the chemotaxis network cannot explain the experimental observations of Vaknin and Berg. We present a new model, in which a small fraction of the phosphatase is localized at the receptor cluster, while the remainder freely diffuses in the cytoplasm; moreover, the phosphatase at the cluster has a higher binding affinity for the messenger protein and a higher catalytic activity than the phosphatase in the cytoplasm. This model is consistent with a large body of experimental data and can explain many of the experimental observations of Vaknin and Berg. More generally, the combination of differential affinity and catalytic activity provides a generic mechanism for amplifying signals that could be exploited in other two-component signaling systems. If this model is correct, then a number of recent modeling studies, which aim to explain the chemotactic gain in terms of the activity of the receptor cluster, should be reconsidered.

In both prokaryotes and eukaryotes, extra- and intracellular signals are often processed by biochemical networks in which two enzymes together control the activity of a messenger protein via opposite modification reactions. A well-known example is the chemotaxis network of Escherichia coli that controls the swimming behavior of the bacterium in response to chemical stimuli. Recent experiments suggest that the two counteracting enzymes in this network are colocalized at the receptor cluster, while experiments by Vaknin and Berg indicate that the spatial distribution of the enzymes by itself can markedly affect the response of the network. We argue using mathematical modeling that the most widely used model of the chemotaxis network is inconsistent with these experimental observations. We then present an alternative model in which part of one enzyme is colocalized with the other enzyme at the receptor cluster, while the remainder freely diffuses in the cytoplasm; moreover, the fraction at the cluster both binds more strongly to the messenger protein and modifies it faster. This model is consistent with a large number of experimental observations and provides a generic mechanism for amplifying signals.