Bacterial motility depends on a critical flagellum length and energy-optimised assembly.

The flagellum is the most complex macromolecular structure known in bacteria and comprised of around two dozen distinct proteins. The main building block of the long, external flagellar filament, flagellin, is secreted through the flagellar type-III secretion system at a remarkable rate of several tens of thousands amino acids per second, significantly surpassing the rates achieved by other pore-based protein secretion systems. The evolutionary implications and potential benefits of this high secretion rate for flagellum assembly and function, however, have remained elusive. In this study, we provide both experimental and theoretical evidence that the flagellar secretion rate has been evolutionarily optimized to facilitate rapid and efficient construction of a functional flagellum. By synchronizing flagellar assembly, we found that a minimal filament length of 2.5 µm was required for swimming motility. Biophysical modelling revealed that this minimal filament length threshold resulted from an elasto-hydrodynamic instability of the whole swimming cell, dependent on the filament length. Furthermore, we developed a stepwise filament labeling method combined with electron microscopy visualization to validate predicted flagellin secretion rates of up to 10,000 amino acids per second. A biophysical model of flagellum growth demonstrates that the observed high flagellin secretion rate efficiently balances filament elongation and energy consumption, thereby enabling motility in the shortest amount of time. Taken together, these insights underscore the evolutionary pressures that have shaped the development and optimization of the flagellum and type-III secretion system, illuminating the intricate interplay between functionality and efficiency in assembly of large macromolecular structures. Significance statement: Our study demonstrates how protein secretion of the bacterial flagellum is finely tuned to optimize filament assembly rate and flagellum function while minimizing energy consumption. By measuring flagellar filament lengths and bacterial swimming after initiation of flag-ellum assembly, we were able to establish the minimal filament length necessary for swimming motility, which we rationalized physically as resulting from an elasto-hydrodynamic instability of the swimming cell. Our bio-physical model of flagellum growth further illustrates how the physiological flagellin secretion rate is optimized to maximize filament elongation while conserving energy. These findings illuminate the evolutionary pressures that have shaped the function of the bacterial flagellum and type-III secretion system, driving improvements in bacterial motility and overall fitness.

Time-reversal symmetry breaking in the chemosensory array reveals mechanisms for dissipation-enhanced cooperative sensing.

The Escherichia coli chemoreceptors form an extensive array that achieves cooperative and adaptive sensing of extracellular signals. The receptors control the activity of histidine kinase CheA, which drives a nonequilibrium phosphorylation-dephosphorylation reaction cycle for response regulator CheY. Cooperativity and dissipation are both important aspects of chemotaxis signaling, yet their consequences have only been studied separately. Recent single-cell FRET measurements revealed that kinase activity of the array spontaneously switches between active and inactive states, with asymmetric switching times that signify time-reversal symmetry breaking in the underlying dynamics. Here, we present a nonequilibrium lattice model of the chemosensory array, which demonstrates that the observed asymmetric switching dynamics can only be explained by an interplay between the dissipative reactions within individual core units and the cooperative coupling between neighboring units. Microscopically, the switching time asymmetry originates from irreversible transition paths. The model shows that strong dissipation enables sensitive and rapid signaling response by relieving the speed-sensitivity trade-off, which can be tested by future single-cell experiments. Overall, our model provides a general framework for studying biological complexes composed of coupled subunits that are individually driven by dissipative cycles and the rich nonequilibrium physics within.

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.

Bifurcation instructed design of multistate machines.

We propose a design paradigm for multistate machines where transitions from one state to another are organized by bifurcations of multiple equilibria of the energy landscape describing the collective interactions of the machine components. This design paradigm is attractive since, near bifurcations, small variations in a few control parameters can result in large changes to the system's state providing an emergent lever mechanism. Further, the topological configuration of transitions between states near such bifurcations ensures robust operation, making the machine less sensitive to fabrication errors and noise. To design such machines, we develop and implement a new efficient algorithm that searches for interactions between the machine components that give rise to energy landscapes with these bifurcation structures. We demonstrate a proof of concept for this approach by designing magnetoelastic machines whose motions are primarily guided by their magnetic energy landscapes and show that by operating near bifurcations we can achieve multiple transition pathways between states. This proof of concept demonstration illustrates the power of this approach, which could be especially useful for soft robotics and at the microscale where typical macroscale designs are difficult to implement.

Resolving the binding-kinase discrepancy in bacterial chemotaxis: A nonequilibrium allosteric model and the role of energy dissipation.

The Escherichia coli chemotaxis signaling pathway has served as a model system for studying 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 this asymmetric shift in binding and kinase response is inconsistent with equilibrium allosteric models regardless of parameter choices. To resolve this inconsistency, we present a nonequilibrium allosteric model that explicitly includes the dissipative reaction cycles driven by ATP hydrolysis. The model successfully explains all existing measurements for both aspartate and serine receptors. Our results suggest that while ligand binding controls the equilibrium balance between the ON and OFF states of the kinase, receptor methylation modulates the kinetic properties (e.g., the phosphorylation rate) of the ON state. Furthermore, sufficient energy dissipation is necessary 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 new perspective on cooperative sensing by large protein complexes and opens up new research directions for understanding their microscopic mechanisms through simultaneous measurements and modeling of ligand binding and downstream responses.

Asymptotic Absorption-Time Distributions in Extinction-Prone Markov Processes.

We characterize absorption-time distributions for birth-death Markov chains with an absorbing boundary. For "extinction-prone" chains (which drift on average toward the absorbing state) the asymptotic distribution is Gaussian, Gumbel, or belongs to a family of skewed distributions. The latter two cases arise when the dynamics slow down dramatically near the boundary. Several models of evolution, epidemics, and chemical reactions fall into these classes; in each case we establish new results for the absorption-time distribution. Applications to African sleeping sickness are discussed.

Cellular Signaling beyond the Wiener-Kolmogorov Limit.

Accurate propagation of signals through stochastic biochemical networks involves significant expenditure of cellular resources. The same is true for regulatory mechanisms that suppress fluctuations in biomolecular populations. Wiener-Kolmogorov (WK) optimal noise filter theory, originally developed for engineering problems, has recently emerged as a valuable tool to estimate the maximum performance achievable in such biological systems for a given metabolic cost. However, WK theory has one assumption that potentially limits its applicability: it relies on a linear, continuum description of the reaction dynamics. Despite this, up to now no explicit test of the theory in nonlinear signaling systems with discrete molecular populations has ever seen performance beyond the WK bound. Here we report the first direct evidence of the bound being broken. To accomplish this, we develop a theoretical framework for multilevel signaling cascades, including the possibility of feedback interactions between input and output. In the absence of feedback, we introduce an analytical approach that allows us to calculate exact moments of the stationary distribution for a nonlinear system. With feedback, we rely on numerical solutions of the system's master equation. The results show WK violations in two common network motifs: a two-level signaling cascade and a negative feedback loop. However, the magnitude of the violation is biologically negligible, particularly in the parameter regime where signaling is most effective. The results demonstrate that while WK theory does not provide strict bounds, its predictions for performance limits are excellent approximations, even for nonlinear systems.

Myosin V executes steps of variable length via structurally constrained diffusion.

The molecular motor myosin V transports cargo by stepping on actin filaments, executing a random diffusive search for actin binding sites at each step. A recent experiment suggests that the joint between the myosin lever arms may not rotate freely, as assumed in earlier studies, but instead has a preferred angle giving rise to structurally constrained diffusion. We address this controversy through comprehensive analytical and numerical modeling of myosin V diffusion and stepping. When the joint is constrained, our model reproduces the experimentally observed diffusion, allowing us to estimate bounds on the constraint energy. We also test the consistency between the constrained diffusion model and previous measurements of step size distributions and the load dependence of various observable quantities. The theory lets us address the biological significance of the constrained joint and provides testable predictions of new myosin behaviors, including the stomp distribution and the run length under off-axis force.

Fitness dependence of the fixation-time distribution for evolutionary dynamics on graphs.

Evolutionary graph theory models the effects of natural selection and random drift on structured populations of competing mutant and nonmutant individuals. Recent studies have found that fixation times in such systems often have right-skewed distributions. Little is known, however, about how these distributions and their skew depend on mutant fitness. Here we calculate the fitness dependence of the fixation-time distribution for the Moran Birth-death process in populations modeled by two extreme networks: the complete graph and the one-dimensional ring lattice, obtaining exact solutions in the limit of large network size. We find that with non-neutral fitness, the Moran process on the ring has normally distributed fixation times, independent of the relative fitness of mutants and nonmutants. In contrast, on the complete graph, the fixation-time distribution is a fitness-weighted convolution of two Gumbel distributions. When fitness is neutral, the fixation-time distribution jumps discontinuously and becomes highly skewed on both the complete graph and the ring. Even on these simple networks, the fixation-time distribution exhibits a rich fitness dependence, with discontinuities and regions of universality. Extensions of our results to two-fitness Moran models, times to partial fixation, and evolution on random networks are discussed.

Modeling the network dynamics of pulse-coupled neurons.

We derive a mean-field approximation for the macroscopic dynamics of large networks of pulse-coupled theta neurons in order to study the effects of different network degree distributions and degree correlations (assortativity). Using the ansatz of Ott and Antonsen [Chaos 18, 037113 (2008)], we obtain a reduced system of ordinary differential equations describing the mean-field dynamics, with significantly lower dimensionality compared with the complete set of dynamical equations for the system. We find that, for sufficiently large networks and degrees, the dynamical behavior of the reduced system agrees well with that of the full network. This dimensional reduction allows for an efficient characterization of system phase transitions and attractors. For networks with tightly peaked degree distributions, the macroscopic behavior closely resembles that of fully connected networks previously studied by others. In contrast, networks with highly skewed degree distributions exhibit different macroscopic dynamics due to the emergence of degree dependent behavior of different oscillators. For nonassortative networks (i.e., networks without degree correlations), we observe the presence of a synchronously firing phase that can be suppressed by the presence of either assortativity or disassortativity in the network. We show that the results derived here can be used to analyze the effects of network topology on macroscopic behavior in neuronal networks in a computationally efficient fashion.

Experimental evidence for the photoisomerization of higher fullerenes.

The initial relaxation dynamics of the photoexcited fullerenes C60, C70, C76, C84, C86, and C90 were investigated by dispersion-free femtosecond pump-probe spectroscopy. Under identical excitation conditions, the formation of the lowest excited state slows down for the larger fullerenes. This trend in dynamics, monitored throughout the visible and NIR range, is found to correlate with the number of isomers in accordance with the isomerization mechanism suggested by Stone and Wales. The Stone-Wales isomerization was calculated as thermally inaccessible but photoinduced barrierless. The energy difference of the isomers is in the 1 meV range, and back-isomerization is observed on the picosecond time scale. The characteristic spectrally broad transient absorption of the investigated fullerenes is promising for fast optical gating applications.