Adaptation dynamics in densely clustered chemoreceptors

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.

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

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

Non-Genetic Diversity in Chemosensing and Chemotactic Behavior

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

Multi-omics response of Pannonibacter phragmitetus BB to hexavalent chromium

The release of hexavalent chromium [Cr(VI)] into water bodies poses a major threat to the environment and human health. However, studies of the biological response to Cr(VI) are limited. In this study, a toxic bacterial mechanism of Cr(VI) was investigated using Pannonibacter phragmitetus BB (hereafter BB), which was isolated from chromate slag. The maximum Cr(VI) concentrations with respect to the resistance and reduction by BB are 4000 mg L^-1 and 2500 mg L^-1, respectively. In the BB genome, more genes responsible for Cr(VI) resistance and reduction are observed compared with other P. phragmitetus strains. A total of 361 proteins were upregulated to respond to Cr(VI) exposure, including enzymes for Cr(VI) uptake, intracellular reduction, ROS detoxification, DNA repair, and Cr(VI) efflux and proteins associated with novel mechanisms involving extracellular reduction mediated by electron transfer, quorum sensing, and chemotaxis. Based on metabolomic analysis, 174 metabolites were identified. Most of the upregulated metabolites are involved in amino acid, glucose, lipid, and energy metabolisms. The results show that Cr(VI) induces metabolite production, while metabolites promote Cr(VI) reduction. Overall, multi-enzyme expression and metabolite production by BB contribute to its high ability to resist/reduce Cr(VI). This study provides details supporting the theory of Cr(VI) reduction and a theoretical basis for the efficient bioremoval of Cr(VI) from the environment.

Optimal methylation noise for best chemotactic performance of E. coli

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

Behavioral Variability and Phenotypic Diversity in Bacterial Chemotaxis

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

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.

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

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

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

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

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

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

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

Emergent properties of bacterial chemotaxis pathway

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

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

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

Networked Chemoreceptors Benefit Bacterial Chemotaxis Performance

Motile bacteria use large receptor arrays to detect and follow chemical gradients in their environment. Extended receptor arrays, composed of networked signaling complexes, promote cooperative stimulus control of their associated signaling kinases. Here, we used structural lesions at the communication interface between core complexes to create an Escherichia coli strain with functional but dispersed signaling complexes. This strain allowed us to directly study how networking of signaling complexes affects chemotactic signaling and gradient-tracking performance. We demonstrate that networking of receptor complexes provides bacterial cells with about 10-fold-heightened detection sensitivity to attractants while maintaining a wide dynamic range over which receptor adaptational modifications can tune response sensitivity. These advantages proved especially critical for chemotaxis toward an attractant source under conditions in which bacteria are unable to alter the attractant gradient.

Chemoreceptor arrays are found in many motile bacteria. However, although our understanding of bacterial chemotaxis is quite detailed, the signaling and behavioral advantages of networked receptor arrays had not been directly studied in cells. We have recently shown that lesions in a key interface of the E. coli receptor array diminish physical connections and functional coupling between core signaling complexes while maintaining their basic signaling capacity. In this study, we exploited an interface 2 mutant to show, for the first time, that coupling between core complexes substantially enhances stimulus detection and chemotaxis performance.

Direct Correlation between Motile Behavior and Protein Abundance in Single Cells

Understanding how stochastic molecular fluctuations affect cell behavior requires the quantification of both behavior and protein numbers in the same cells. Here, we combine automated microscopy with in situ hydrogel polymerization to measure single-cell protein expression after tracking swimming behavior. We characterized the distribution of non-genetic phenotypic diversity in Escherichia coli motility, which affects single-cell exploration. By expressing fluorescently tagged chemotaxis proteins (CheR and CheB) at different levels, we quantitatively mapped motile phenotype (tumble bias) to protein numbers using thousands of single-cell measurements. Our results disagreed with established models until we incorporated the role of CheB in receptor deamidation and the slow fluctuations in receptor methylation. Beyond refining models, our central finding is that changes in numbers of CheR and CheB affect the population mean tumble bias and its variance independently. Therefore, it is possible to adjust the degree of phenotypic diversity of a population by adjusting the global level of expression of CheR and CheB while keeping their ratio constant, which, as shown in previous studies, confers functional robustness to the system. Since genetic control of protein expression is heritable, our results suggest that non-genetic diversity in motile behavior is selectable, supporting earlier hypotheses that such diversity confers a selective advantage.

Cell-to-cell variations in protein numbers due to random fluctuations at the molecular level lead to cell-to-cell variations in behavior. To maintain predictable responses, signaling networks have evolved robustness against noise, but in some situations phenotypic diversity in a clonal population can be beneficial as a bet hedging or division of labor strategy. Investigating of how random molecular fluctuations affect cell behavior requires to measure biological parameters at different scales. Here, we report a new experiment that allows the measure of both protein numbers and behavior in cells that are free to move in their environment. Using Escherichia coli, a model system for the study of cellular behavior, we investigated the effects variations in the numbers of the chemo-receptor modification enzymes on single-cell swimming behavior. We found that the mean and variance of the behavior can be adjusted independently in the population by adjusting protein expression. This mechanism allows for the genetic selection of phenotypic diversity without disrupting correlations in protein expression that are important for the overall robustness of the chemotaxis system.

Allosteric proteins as logarithmic sensors

Many sensory systems, from vision and hearing in animals to signal transduction in cells, respond to fold changes in signal relative to background. Responding to fold change requires that the system senses signal on a logarithmic scale, responding identically to a change in signal level from 1 to 3, or from 10 to 30. It is an ongoing search in the field to understand the ways in which a logarithmic sensor can be implemented at the molecular level. In this work, we present evidence that logarithmic sensing can be implemented with a single protein, by means of allosteric regulation. Specifically, we find that mathematical models show that allosteric proteins can respond to stimuli on a logarithmic scale. Next, we present evidence from measurements in the literature that some allosteric proteins do operate in a parameter regime that permits logarithmic sensing. Finally, we present examples suggesting that allosteric proteins are indeed used in this capacity: allosteric proteins play a prominent role in systems where fold-change detection has been proposed. This finding suggests a role as logarithmic sensors for the many allosteric proteins across diverse biological processes.

Fundamental constraints on the abundances of chemotaxis proteins

Flagellated bacteria, such as Escherichia coli, perform directed motion in gradients of concentration of attractants and repellents in a process called chemotaxis. The E. coli chemotaxis signaling pathway is a model for signal transduction, but it has unique features. We demonstrate that the need for fast signaling necessitates high abundances of the proteins involved in this pathway. We show that further constraints on the abundances of chemotaxis proteins arise from the requirements of self-assembly both of flagellar motors and of chemoreceptor arrays. All these constraints are specific to chemotaxis, and published data confirm that chemotaxis proteins tend to be more highly expressed than their homologs in other pathways. Employing a chemotaxis pathway model, we show that the gain of the pathway at the level of the response regulator CheY increases with overall chemotaxis protein abundances. This may explain why, at least in one E. coli strain, the abundance of all chemotaxis proteins is higher in media with lower nutrient content. We also demonstrate that the E. coli chemotaxis pathway is particularly robust to abundance variations of the motor protein FliM.

Adaptability of non-genetic diversity in bacterial chemotaxis

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

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

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

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

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

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

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

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

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

Limits of feedback control in bacterial chemotaxis

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

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