Yuhai Tu

Weighted Evolving Networks

Many biological, ecological, and economic systems are best described by weighted networks, as the nodes interact with each other with varying strength. However, most evolving network models studied so far are binary, the link strength being either 0 or 1. In this paper we introduce and investigate the scaling properties of a class of models which assign weights to the links as the network evolves. The combined numerical and analytical approach indicates that asymptotically the total weight distribution converges to the scaling behavior of the connectivity distribution, but this convergence is hampered by strong logarithmic corrections.

Modeling the chemotactic response of Escherichia coli to time-varying stimuli

In their natural environment, cells need to extract useful information from complex temporal signals that vary over a wide range of intensities and time scales. Here, we study how such signals are processed by Escherichia coli during chemotaxis by developing a general theoretical model based on receptor adaptation and receptor–receptor cooperativity. Measured responses to various monotonic, oscillatory, and impulsive stimuli are all explained consistently by the underlying adaptation kinetics within this model. For exponential ramp signals, an analytical solution is discovered that reveals a remarkable connection between the dependence of kinase activity on the exponential ramp rate and the receptor methylation rate function. For exponentiated sine-wave signals, spectral analysis shows that the chemotaxis pathway acts as a lowpass filter for the derivative of the signal with the cutoff frequency determined by an intrinsic adaptation time scale. For large step stimuli, we find that the recovery time is determined by the constant maximum methylation rate, which provides a natural explanation for the observed recovery time additivity. Our model provides a quantitative system-level description of the chemotaxis signaling pathway and can be used to predict E. coli chemotaxis responses to arbitrary temporal signals. This model of the receptor system reveals the molecular origin of Webers law in bacterial chemotaxis. We further identify additional constraints required to account for the related observation that the output of this pathway is constant under exponential ramp stimuli, a feature that we call “logarithmic tracking.”

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

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.

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.

The energy-speed-accuracy trade-off in sensory adaptation

Adaptation is the essential process by which an organism becomes better suited to its environment. The benefits of adaptation are well documented, but the cost it incurs remains poorly understood. Here, by analysing a stochastic model of a minimum feedback network underlying many sensory adaptation systems, we show that adaptive processes are necessarily dissipative, and continuous energy consumption is required to stabilize the adapted state. Our study reveals a general relation among energy dissipation rate, adaptation speed and the maximum adaptation accuracy. This energy-speed-accuracy relation is tested in the Escherichia coli chemosensory system, which exhibits near-perfect chemoreceptor adaptation. We identify key requirements for the underlying biochemical network to achieve accurate adaptation with a given energy budget. Moreover, direct measurements confirm the prediction that adaptation slows down as cells gradually de-energize in a nutrient-poor medium without compromising adaptation accuracy. Our work provides a general framework to study cost-performance tradeoffs for cellular regulatory functions and information processing.

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

We propose a general theoretical framework for modeling receptor sensitivity in bacterial chemotaxis, taking into account receptor interactions, including those among different receptor species. We show that our model can quantitatively explain the recent in vivo measurements of receptor sensitivity at different ligand concentrations for both mutant and wild-type strains. For mutant strains, our model can fit the experimental data exactly. For the wild-type cell, our model is capable of achieving high gain while having modest values of Hill coefficient for the response curves. Furthermore, the high sensitivity of the wild-type cell in our model is maintained for a wide range of ambient ligand concentrations, facilitated by near-perfect adaptation and dependence of ligand binding on receptor activity. Our study reveals the importance of coupling among different chemoreceptor species, in particular strong interactions between the aspartate (Tar) and serine (Tsr) receptors, which is crucial in explaining both the mutant and wild-type data. Predictions for the sensitivity of other mutant strains and possible improvements of our model for the wild-type cell are also discussed.

Perfect and Near-Perfect Adaptation in a Model of Bacterial Chemotaxis

The signaling apparatus mediating bacterial chemotaxis can adapt to a wide range of persistent external stimuli. In many cases, the bacterial activity returns to its prestimulus level exactly, and this perfect adaptability is robust against variations in various chemotaxis protein concentrations. We model the bacterial chemotaxis signaling pathway, from ligand binding to CheY phosphorylation. By solving the steady-state equations of the model analytically, we derive a full set of conditions for the system to achieve perfect adaptation. The conditions related to the phosphorylation part of the pathway are discovered for the first time, while other conditions are generalizations of the ones found in previous works. Sensitivity of the perfect adaptation is evaluated by perturbing these conditions. We find that, even in the absence of some of the perfect adaptation conditions, adaptation can be achieved with near-perfect precision as a result of the separation of scales in both chemotaxis protein concentrations and reaction rates, or specific properties of the receptor distribution in different methylation states. Since near-perfect adaptation can be found in much larger regions of the parameter space than that defined by the perfect adaptation conditions, their existence is essential to understand robustness in bacterial chemotaxis.

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

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.

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.