Dynamics of cooperativity in chemical sensing among cell-surface receptors

Limits on the accuracy of contact inhibition of locomotion

Cells that collide with each other repolarize away from contact, in a process called contact inhibition of locomotion (CIL), which is necessary for correct development of the embryo. CIL can occur even when cells make a micron-scale contact with a neighbor-much smaller than their size. How precisely can a cell sense cell-cell contact and repolarize in the correct direction? What factors control whether a cell recognizes it has contacted a neighbor? We propose a theoretical model for the limits of CIL where cells recognize the presence of another cell by binding the protein ephrin with the Eph receptor. This recognition is made difficult by the presence of interfering ligands that bind nonspecifically. Both theoretical predictions and simulation results show that it becomes more difficult to sense cell-cell contact when it is difficult to distinguish ephrin from the interfering ligands, or when there are more interfering ligands, or when the contact width decreases. However, the error of estimating contact position remains almost constant when the contact width changes. This happens because the cell gains spatial information largely from the boundaries of cell-cell contact. We study using statistical decision theory the likelihood of a false-positive CIL event in the absence of cell-cell contact, and the likelihood of a false negative where CIL does not occur when another cell is present. Our results suggest that the cell is more likely to make incorrect decisions when the contact width is very small or so large that it nears the cell's perimeter. However, in general, we find that cells have the ability to make reasonably reliable CIL decisions even for very narrow (micron-scale) contacts, even if the concentration of interfering ligands is ten times that of the correct ligands.

How does an organism extract relevant information from transcription factor concentrations?

How does an organism regulate its genes? The involved regulation typically occurs in terms of a signal processing chain: an externally applied stimulus or a maternally supplied transcription factor leads to the expression of some downstream genes, which, in turn, are transcription factors for further genes. Especially during development, these transcription factors are frequently expressed in amounts where noise is still important; yet, the signals that they provide must not be lost in the noise. Thus, the organism needs to extract exactly relevant information in the signal. New experimental approaches involving single-molecule measurements at high temporal precision as well as increased precision in manipulations directly on the genome are allowing us to tackle this question anew. These new experimental advances mean that also from the theoretical side, theoretical advances should be possible. In this review, I will describe, specifically on the example of fly embryo gene regulation, how theoretical approaches, especially from inference and information theory, can help in understanding gene regulation. To do so, I will first review some more traditional theoretical models for gene regulation, followed by a brief discussion of information-theoretical approaches and when they can be applied. I will then introduce early fly development as an exemplary system where such information-theoretical approaches have traditionally been applied and can be applied; I will specifically focus on how one such method, namely the information bottleneck approach, has recently been used to infer structural features of enhancer architecture.

Pleiotropy enables specific and accurate signaling in the presence of ligand cross talk

Living cells sense their environment through the binding of extracellular molecular ligands to cell surface receptors. Puzzlingly, vast numbers of signaling pathways exhibit a high degree of cross talk between different signals whereby different ligands act through the same receptor or shared components downstream. It remains unclear how a cell can accurately process information from the environment in such cross-wired pathways. We show that a feature which commonly accompanies cross talk-signaling pleiotropy (the ability of a receptor to produce multiple outputs)-offers a solution to the cross-talk problem. In a minimal model we show that a single pleiotropic receptor can simultaneously identify and accurately sense the concentrations of arbitrary unknown ligands present individually or in a mixture. We calculate the fundamental limits of the signaling specificity and accuracy of such signaling schemes. The model serves as an elementary "building block" toward understanding more complex cross-wired receptor-ligand signaling networks.

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

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

Multicellular sensing at a feedback-induced critical point

Feedback in sensory biochemical networks can give rise to bifurcations in cells' behavioral response. These bifurcations share many properties with thermodynamic critical points. Evidence suggests that biological systems may operate near these critical points, but the functional benefit of doing so remains poorly understood. Here we investigate a simple biochemical model with nonlinear feedback and multicellular communication to determine if criticality provides a functional benefit in terms of the ability to gain information about a stochastic chemical signal. We find that when signal fluctuations are slow, the mutual information between the signal and the intracellular readout is maximized at criticality, because the benefit of high signal susceptibility outweighs the detriment of high readout noise. When cells communicate, criticality gives rise to long-range correlations in readout molecule number among cells. Consequently, we find that communication increases the mutual information between a given cell's readout and the spatial average of the signal across the population. Finally, we find that both with and without communication, the sensory benefits of criticality compete with critical slowing down, such that the information rate, as opposed to the information itself, is minimized at the critical point. Our results reveal the costs and benefits of feedback-induced criticality for multicellular sensing.

Energy consumption and cooperation for optimal sensing

The reliable detection of environmental molecules in the presence of noise is an important cellular function, yet the underlying computational mechanisms are not well understood. We introduce a model of two interacting sensors which allows for the principled exploration of signal statistics, cooperation strategies and the role of energy consumption in optimal sensing, quantified through the mutual information between the signal and the sensors. Here we report that in general the optimal sensing strategy depends both on the noise level and the statistics of the signals. For joint, correlated signals, energy consuming (nonequilibrium), asymmetric couplings result in maximum information gain in the low-noise, high-signal-correlation limit. Surprisingly we also find that energy consumption is not always required for optimal sensing. We generalise our model to incorporate time integration of the sensor state by a population of readout molecules, and demonstrate that sensor interaction and energy consumption remain important for optimal sensing.

Cells exhibit exceptional chemical sensitivity, yet we haven’t fully understood how they achieve it. Here the authors consider the mutual information between signals and two coupled sensors as a proxy for sensing performance and show its optimisation depending on noise level and signal statistics.

Maximal information transmission is compatible with ultrasensitive biological pathways

Cells are often considered input-output devices that maximize the transmission of information by converting extracellular stimuli (input) via signaling pathways (communication channel) to cell behavior (output). However, in biological systems outputs might feed back into inputs due to cell motility, and the biological channel can change by mutations during evolution. Here, we show that the conventional channel capacity obtained by optimizing the input distribution for a fixed channel may not reflect the global optimum. In a new approach we analytically identify both input distributions and input-output curves that optimally transmit information, given constraints from noise and the dynamic range of the channel. We find a universal optimal input distribution only depending on the input noise, and we generalize our formalism to multiple outputs (or inputs). Applying our formalism to Escherichia coli chemotaxis, we find that its pathway is compatible with optimal information transmission despite the ultrasensitive rotary motors.

Critical slowing down in biochemical networks with feedback

Near a bifurcation point, the response time of a system is expected to diverge due to the phenomenon of critical slowing down. We investigate critical slowing down in well-mixed stochastic models of biochemical feedback by exploiting a mapping to the mean-field Ising universality class. We analyze the responses to a sudden quench and to continuous driving in the model parameters. In the latter case, we demonstrate that our class of models exhibits the Kibble-Zurek collapse, which predicts the scaling of hysteresis in cellular responses to gradual perturbations. We discuss the implications of our results in terms of the tradeoff between a precise and a fast response. Finally, we use our mapping to quantify critical slowing down in T cells, where the addition of a drug is equivalent to a sudden quench in parameter space.

Receptor crosstalk improves concentration sensing of multiple ligands

Cells need to reliably sense external ligand concentrations to achieve various biological functions such as chemotaxis or signaling. The molecular recognition of ligands by surface receptors is degenerate in many systems, leading to crosstalk between ligand-receptor pairs. Crosstalk is often thought of as a deviation from optimal specific recognition, as the binding of noncognate ligands can interfere with the detection of the receptor's cognate ligand, possibly leading to a false triggering of a downstream signaling pathway. Here we quantify the optimal precision of sensing the concentrations of multiple ligands by a collection of promiscuous receptors. We demonstrate that crosstalk can improve precision in concentration sensing and discrimination tasks. To achieve superior precision, the additional information about ligand concentrations contained in short binding events of the noncognate ligand should be exploited. We present a proofreading scheme to realize an approximate estimation of multiple ligand concentrations that reaches a precision close to the derived optimal bounds. Our results help rationalize the observed ubiquity of receptor crosstalk in molecular sensing.

Precision in a rush: Trade-offs between reproducibility and steepness of the hunchback expression pattern

Fly development amazes us by the precision and reproducibility of gene expression, especially since the initial expression patterns are established during very short nuclear cycles. Recent live imaging of hunchback promoter dynamics shows a stable steep binary expression pattern established within the three minute interphase of nuclear cycle 11. Considering expression models of different complexity, we explore the trade-off between the ability of a regulatory system to produce a steep boundary and minimize expression variability between different nuclei. We show how a limited readout time imposed by short developmental cycles affects the gene’s ability to read positional information along the embryo’s anterior posterior axis and express reliably. Comparing our theoretical results to real-time monitoring of the hunchback transcription dynamics in live flies, we discuss possible regulatory strategies, suggesting an important role for additional binding sites, gradients or non-equilibrium binding and modified transcription factor search strategies.

Despite very limited time, organisms develop in reproducible ways. In the early stages of fly development the information about maternal signals is read out in a few minutes to produce steep and precise gene expression patterns. Motivated by recent live imaging experiments in fly embryos, we explore the consequences of the trade-off between a rushed but reproducible readout and a steep expression pattern on the regulatory modules of gene expression. We show that the current view of one anterior gradient morphogen binding to six binding sites is quantitatively inconsistent with the experimental data given the short readout time, suggesting other regulatory features.

Propagation of regulatory fluctuations induces coordinated switching of flagellar motors in chemotaxis signaling pathway of single bacteria

The random motion of E. coli is driven by multiple flagella motors. When all motors rotate in the counter clockwise direction, the bacteria swims smoothly. A recent experimental report by Terasawa et al. [Biophys J,100,2193,(2011)] demonstrated that a coordination of the motors can occur through signaling pathways, and perturbation of a regulatory molecule disrupted the coordination. Here, we develop a mathematical model to show that a large temporal fluctuation in the regulator concentration can induce a correlated switching of the multiple motors. Such a large fluctuation is generated by a chemotaxis receptor cluster in unilateral cell pole, which then exhibits a spatial propagation through the cytoplasm from the receptor position to the motor around cell periphery. Our numerical simulation successfully reproduces synchronized switching and the lag time in the motions of two distant motors, which has been observed experimentally. We further show that the large fluctuation in the regulator concentration at the motor positions can expand the dynamic range that the motor can respond, which confers robustness to the signaling system.

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.

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.

Predicting Escherichia coli's chemotactic drift under exponential gradient

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

A Boolean approach to bacterial chemotaxis

Bacterium such as Escherichia coli (E. coli) show biased Brownian motion in different chemical concentration gradients. This chemical sensitive motility or chemotaxis has gained considerable interest among scientists for some remarkable features such as chemo-sensory dynamic range, adaptation, diffusion and drift. A Boolean model of the whole chemotaxis process has been developed in this manuscript. The response of the circuit is in accordance with the experimental results available in the literature, providing indirect validation of the model. This simple Boolean network (BN) can be easily integrated into the paradigm of modular whole cell modelling. Another crucial application is in designing bio-inspired micro-robots to detect certain spatio-temporal chemical signatures.

Topologically protected modes in non-equilibrium stochastic systems

Non-equilibrium driving of biophysical processes is believed to enable their robust functioning despite the presence of thermal fluctuations and other sources of disorder. Such robust functions include sensory adaptation, enhanced enzymatic specificity and maintenance of coherent oscillations. Elucidating the relation between energy consumption and organization remains an important and open question in non-equilibrium statistical mechanics. Here we report that steady states of systems with non-equilibrium fluxes can support topologically protected boundary modes that resemble similar modes in electronic and mechanical systems. Akin to their electronic and mechanical counterparts, topological-protected boundary steady states in non-equilibrium systems are robust and are largely insensitive to local perturbations. We argue that our work provides a framework for how biophysical systems can use non-equilibrium driving to achieve robust function.

Energy dissipation characterizes the states far from equilibrium, whilst how it affects the local organization remains elusive. Here, Murugan et al. show that the non-equilibrium systems exhibit topologically protected boundary modes that have been known in electronic and mechanical systems.

Cooperative Clustering Digitizes Biochemical Signaling and Enhances its Fidelity

Many membrane-bound molecules in cells form small clusters. It has been hypothesized that these clusters convert an analog extracellular signal into a digital intracellular signal and that this conversion increases signaling fidelity. However, the mechanism by which clusters digitize a signal and the subsequent effects on fidelity remain poorly understood. Here we demonstrate using a stochastic model of cooperative cluster formation that sufficient cooperation leads to digital signaling. We show that despite reducing the number of output states, which decreases fidelity, digitization also reduces noise in the system, which increases fidelity. The tradeoff between these effects leads to an optimal cluster size that agrees with experimental measurements.

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.

The principle of conformational signaling

Specific conformations of signaling proteins can serve as “signals” in signal transduction by being recognized by receptors.

Dispersion for two classes of random variables: general theory and application to inference of an external ligand concentration by a cell

We derive expressions for the dispersion for two classes of random variables in Markov processes. Random variables such as current and activity pertain to the first class, which is composed of random variables that change whenever a jump in the stochastic trajectory occurs. The second class corresponds to the time the trajectory spends in a state (or cluster of states). While the expression for the first class follows straightforwardly from known results in the literature, we show that a similar formalism can be used to derive an expression for the second class. As an application, we use this formalism to analyze a cellular two-component network estimating an external ligand concentration. The uncertainty related to this external concentration is calculated by monitoring different random variables related to an internal protein. We show that, inter alia, monitoring the time spent in the phosphorylated state of the protein leads to a finite uncertainty only if there is dissipation, whereas the uncertainty obtained from the activity of the transitions of the internal protein can reach the Berg-Purcell limit even in equilibrium.

Juxtacrine signaling is inherently noisy

Juxtacrine signaling is an important class of signaling systems that plays a crucial role in various developmental processes ranging from coordination of differentiation between neighboring cells to guiding axon growth during neurogenesis. Such signaling systems rely on the interaction between receptors on one cell and trans-membrane ligands on the membrane of a neighboring cell. Like other signaling systems, the ability of signal-receiving cells to accurately determine the concentration of ligands, is affected by stochastic diffusion processes. However, it is not clear how restriction of ligand movement to the two-dimensional (2D) cell membrane in juxtacrine signaling affects the accuracy of ligand sensing. In this study, we use a statistical mechanics approach, to show that long integration times, from around one second to several hours, are required to reach high-sensing accuracy (better than 10%). Surprisingly, the accuracy of sensing cannot be significantly improved, neither by increasing the number of receptors above three to five receptors per contact area, nor by increasing the contact area between cells. We show that these results impose stringent constraints on the dynamics of processes relying on juxtacrine signaling systems, such as axon guidance mediated by Ephrins and developmental patterns mediated by the Notch pathway.

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.

Dynamic membrane patterning, signal localization and polarity in living cells

We review the molecular and physical aspects of the dynamic localization of signaling molecules on the plasma membranes of living cells. At the nanoscale, clusters of receptors and signaling proteins play an essential role in the processing of extracellular signals. At the microscale, "soft" and highly dynamic signaling domains control the interaction of individual cells with their environment. At the multicellular scale, individual polarity patterns control the forces that shape multicellular aggregates and tissues.

Chemodetection in fluctuating environments: receptor coupling, buffering, and antagonism

Variability in the chemical composition of the extracellular environment can significantly degrade the ability of cells to detect rare cognate ligands. Using concepts from statistical detection theory, we formalize the generic problem of detection of small concentrations of ligands in a fluctuating background of biochemically similar ligands binding to the same receptors. We discover that in contrast with expectations arising from considerations of signal amplification, inhibitory interactions between receptors can improve detection performance in the presence of substantial environmental variability, providing an adaptive interpretation to the phenomenon of ligand antagonism. Our results suggest that the structure of signaling pathways responsible for chemodetection in fluctuating and heterogeneous environments might be optimized with respect to the statistics and dynamics of environmental composition. The developed formalism stresses the importance of characterizing nonspecific interactions to understand function in signaling pathways.

Energy dissipation and noise correlations in biochemical sensing

To measure chemical concentrations, cells need to extract information from stochastic receptor signals via signaling networks which are also inherently stochastic. Here, we study how the accuracy of sensing depends on the correlations between these extrinsic and intrinsic sources of noise. We find that the sensing precision of signaling networks that are not driven out of equilibrium is fundamentally limited by the fluctuation-dissipation theorem, which generates a tradeoff between the removal of extrinsic and intrinsic noise. As a result, the sensing precision of equilibrium systems is limited by the number of receptors; the downstream network can never improve sensing. To lift the tradeoff, energy dissipation is essential. This allows the receptor to transduce the signal as a catalyst and enables time integration of the receptor state. To beat the sensing limit of equilibrium systems, a canonical nonequilibrium signaling network based on the push-pull motif needs to dissipate at least 1k_{B}T per receptor.

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.

Predicting chemical environments of bacteria from receptor signaling

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

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

How input noise limits biochemical sensing in ultrasensitive systems

Many biological processes are regulated by molecular devices that respond in an ultrasensitive fashion to upstream signals. An important question is whether such ultrasensitivity improves or limits its ability to read out the (noisy) input stimuli. Here, we develop a simple model to study the statistical properties of ultrasensitive signaling systems. We demonstrate that the output sensory noise is always bounded, in contrast to earlier theories using the small noise approximation, which tends to overestimate the impact of noise in ultrasensitive pathways. Our analysis also shows that the apparent sensitivity of the system is ultimately constrained by the input signal-to-noise ratio. Thus, ultrasensitivity can improve the precision of biochemical sensing only to a finite extent. This corresponds to a new limit for ultrasensitive signaling systems, which is strictly tighter than the Berg-Purcell limit.

The biophysical model for accuracy of cellular sensing spatial gradients of multiple chemoattractants

Spatial gradients of surrounding chemoattractants are the key factors in determining the directionality of eukaryotic cell movement. Thus, it is important for cells to accurately measure the spatial gradients of surrounding chemoattractants. Here, we study the precision of sensing the spatial gradients of multiple chemoattractants using cooperative receptor clusters. Cooperative receptors on cells are modeled as an Ising chain of Monod-Wyman-Changeux clusters subject to multiple chemical-gradient fields to study the physical limits of multiple chemoattractants spatial gradients sensing. We found that eukaryotic cells cannot sense each chemoattractant gradient individually. Instead, cells can only sense a weighted sum of surrounding chemical gradients. Moreover, the precision of sensing one chemical gradient is signicantly affected by coexisting chemoattractant concentrations. These findings can provide a further insight into the role of chemoattractants in immune response and help develop novel treatments for inflammatory diseases.

Chemical sensing by nonequilibrium cooperative receptors

Cooperativity arising from local interactions in equilibrium receptor systems provides gain, but does not increase sensory performance, as measured by the signal-to-noise ratio (SNR) due to a fundamental tradeoff between gain and intrinsic noise. Here we allow sensing to be a nonequilibrium process and show that energy dissipation cannot circumvent the fundamental tradeoff, so that the SNR is still optimal for independent receptors. For systems requiring high gain, nonequilibrium 2D-coupled receptors maximize the SNR, revealing a new design principle for biological sensors.

Defining a key receptor-CheA kinase contact and elucidating its function in the membrane-bound bacterial chemosensory array: a disulfide mapping and TAM-IDS Study

The three core components of the ubiquitous bacterial chemosensory array - the transmembrane chemoreceptor, the histidine kinase CheA, and the adaptor protein CheW - assemble to form a membrane-bound, hexagonal lattice in which receptor transmembrane signals regulate kinase activity. Both the regulatory domain of the kinase and the adaptor protein bind to overlapping sites on the cytoplasmic tip of the receptor (termed the protein interaction region). Notably, the kinase regulatory domain and the adaptor protein share the same fold constructed of two SH3-like domains. The present study focuses on the structural interface between the receptor and the kinase regulatory domain. Two models have been proposed for this interface: Model 1 is based on the crystal structure of a homologous Thermotoga complex between a receptor fragment and the CheW adaptor protein. This model has been used in current models of chemosensory array architecture to build the receptor-CheA kinase interface. Model 2 is based on a newly determined crystal structure of a homologous Thermotoga complex between a receptor fragment and the CheA kinase regulatory domain. Both models present unique strengths and weaknesses, and current evidence is unable to resolve which model best describes contacts in the native chemosensory arrays of Escherichia coli, Salmonella typhimurium, and other bacteria. Here we employ disulfide mapping and tryptophan and alanine mutation to identify docking sites (TAM-IDS) to test Models 1 and 2 in well-characterized membrane-bound arrays formed from E. coli and S. typhimurium components. The results reveal that the native array interface between the receptor protein interaction region and the kinase regulatory domain is accurately described by Model 2, but not by Model 1. In addition, the results show that the interface possesses both a structural function that contributes to stable CheA kinase binding in the array and a regulatory function central to transmission of the activation signal from receptor to CheA kinase. On-off switching alters the disulfide formation rates of specific Cys pairs at the interface, but not most Cys pairs, indicating that signaling perturbs localized regions of the interface. The findings suggest a simple model for the rearrangement of the interface triggered by the attractant signal and for longer range transmission of the signal in the chemosensory array.

Spatial partitioning improves the reliability of biochemical signaling

Spatial heterogeneity is a hallmark of living systems, even at the molecular scale in individual cells. A key example is the partitioning of membrane-bound proteins via lipid domain formation or cytoskeleton-induced corralling. However, the impact of this spatial heterogeneity on biochemical signaling processes is poorly understood. Here, we demonstrate that partitioning improves the reliability of biochemical signaling. We exactly solve a stochastic model describing a ubiquitous motif in membrane signaling. The solution reveals that partitioning improves signaling reliability via two effects: it moderates the nonlinearity of the switching response, and it reduces noise in the response by suppressing correlations between molecules. An optimal partition size arises from a trade-off between minimizing the number of proteins per partition to improve signaling reliability and ensuring sufficient proteins per partition to maintain signal propagation. The predicted optimal partition size agrees quantitatively with experimentally observed systems. These results persist in spatial simulations with explicit diffusion barriers. Our findings suggest that molecular partitioning is not merely a consequence of the complexity of cellular substructures, but also plays an important functional role in cell signaling.

Trapping of diffusing particles by clusters of absorbing disks on a reflecting wall with disk centers on sites of a square lattice

A simple approximate formula is derived for the rate constant that describes steady-state flux of diffusing particles through a cluster of perfectly absorbing disks on the otherwise reflecting flat wall, assuming that the disk centers occupy neighboring sites of a square lattice. A distinctive feature of trapping by a disk cluster is that disks located at the cluster periphery shield the disks in the center of the cluster. This competition of the disks for diffusing particles makes it impossible to find an exact analytical solution for the rate constant in the general case. To derive the approximate formula, we use a recently suggested approach [A. M. Berezhkovskii, L. Dagdug, V. A. Lizunov, J. Zimmerberg, and S. M. Bezrukov, J. Chem. Phys. 136, 211102 (2012)], which is based on the replacement of the disk cluster by an effective uniform partially absorbing spot. The formula shows how the rate constant depends on the size and shape of the cluster. To check the accuracy of the formula, we compare its predictions with the values of the rate constant obtained from Brownian dynamics simulations. The comparison made for 18 clusters of various shapes and sizes shows good agreement between the theoretical predictions and numerical results.

Fundamental limits on sensing chemical concentrations with linear biochemical networks

Living cells often need to extract information from biochemical signals that are noisy. We study how accurately cells can measure chemical concentrations with signaling networks that are linear. For stationary signals of long duration, they can reach, but not beat, the Berg-Purcell limit, which relies on uniformly averaging in time the fluctuations in the input signal. For short times or nonstationary signals, however, they can beat the Berg-Purcell limit, by nonuniformly time averaging the input. We derive the optimal weighting function for time averaging and use it to provide the fundamental limit of measuring chemical concentrations with linear signaling networks.

Communication: Clusters of absorbing disks on a reflecting wall: competition for diffusing particles

Trapping of diffusing particles by a cluster of absorbing disks on the otherwise reflecting wall is a manifestly many-body problem because of the disk competition for the particles. By replacing the cluster with an effective uniformly absorbing spot, we derive a simple formula for the rate constant that characterizes the trapping. The formula shows how the rate constant depends on the size and shape of the cluster.