Motility of Escherichia coli in a quasi-two-dimensional porous medium

Dispersion of run-and-tumble microswimmers through disordered media

Understanding the transport properties of microorganisms and self-propelled particles in porous media has important implications for human health as well as microbial ecology. In free space, most microswimmers perform diffusive random walks as a result of the interplay of self-propulsion and orientation decorrelation mechanisms such as run-and-tumble dynamics or rotational diffusion. In an unstructured porous medium, collisions with the microstructure result in a decrease in the effective spatial diffusivity of the particles from its free-space value. Here, we analyze this problem for a simple model system consisting of noninteracting point particles performing run-and-tumble dynamics through a two-dimensional disordered medium composed of a random distribution of circular obstacles, in the absence of Brownian diffusion or hydrodynamic interactions. The particles are assumed to collide with the obstacles as hard spheres and subsequently slide on the obstacle surface with no frictional resistance while maintaining their orientation, until they either escape or tumble. We show that the variations in the long-time diffusivity can be described by a universal dimensionless hindrance function f(ϕ,Pe) of the obstacle area fraction ϕ and Péclet number Pe, or ratio of the swimmer run length to the obstacle size. We analytically derive an asymptotic expression for the hindrance function valid for dilute media (Peϕ≪1), and its extension to denser media is obtained using stochastic simulations. As we explain, the model is also easily generalized to describe dispersion in three dimensions.

Dichotomic ratchet in a two-dimensional corrugated channel

We consider a particle diffusing in a two-dimensional (2D) channel of varying width h(x). It is driven by a force of constant magnitude f but random orientation there or back along the channel. We derive the effective generalized Fick-Jacobs equation for this system, which describes the dynamics of such a particle in the longitudinal coordinate x. Aside from the effective diffusion coefficient D(x), our mapping also generates an additional effective potential -γ(x) added to the entropic potential -log[h(x)]. It acquires an increasing or decreasing component in asymmetric periodic channels, and thus it explains appearance of the ratchet current. We study this effect on a trial example and compare the results of our true 2D theory with a commonly used effective one-dimensional description; the data are verified by the numerical solution of the full 2D problem.

Polymer looping kinetics in active heterogeneous environments

A typical biological environment is usually featured by crowding and heterogeneity, leading to complex reaction kinetics of the immersed macromolecules. In the present work, we adopt Langevin dynamics simulations to systematically investigate polymer looping kinetics in active heterogeneous media crowded with a mixture of mobile active particles and immobile obstacles. For comparison, a parallel study is also performed in the passive heterogeneous media. We explicitly analyze the change of looping time and looping probability with the variation of obstacle ratio, volume fraction and crowder size. We reveal the novel phenomena of inhibition-facilitation transition of the looping rate induced by heterogeneity, crowdedness and activity. In addition, our results demonstrate a very non-trivial crowder size effect on the looping kinetics. The underlying mechanism is rationalized by the interplay of polymer diffusion, conformational change and looping free-energy barrier. The competing effect arising from active particles and obstacles on structural and dynamical properties of the polymer yields a consistent scenario for our observations. Lastly, the non-exponential kinetics of the looping process is also analyzed. We find that both activity and crowding can strengthen the heterogeneity degree of the looping kinetics.

Impact of confined geometries on hopping and trapping of motile bacteria in porous media

We use a random walk particle-tracking (RWPT) approach to elucidate the impact of porous media confinement and cell-cell interactions on bacterial transport. The model employs stochastic alternating motility states consisting of hopping movement and trapping reorientation. The stochastic motility patterns are defined based on direct visualization of individual trajectory data. We validate our model against experimental data, at single-cell resolution, of bacterial E. coli motion in three-dimensional confined porous media. Results show that the model is able to efficiently simulate the spreading dynamics of motile bacteria as it captures the impact of cell-cell interaction and pore confinement, which marks the transition to a late-time subdiffusive regime. Furthermore, the model is able to qualitatively reproduce the observed directional persistence. Our RWPT model constitutes a meshless simple method which is easy to implement and does not invoke ad hoc assumptions but represents the basis for a multiscale approach to the study of bacterial dispersal in porous systems.

Critical patch size reduction by heterogeneous diffusion

Population survival depends on a large set of factors and on how they are distributed in space. Due to landscape heterogeneity, species can occupy particular regions that provide the ideal scenario for development, working as a refuge from harmful environmental conditions. Survival occurs if population growth overcomes the losses caused by adventurous individuals that cross the patch edge. In this work, we consider a single species dynamics in a patch with a space-dependent diffusion coefficient. We show analytically, within the Stratonovich framework, that heterogeneous diffusion reduces the minimal patch size for population survival when contrasted with the homogeneous case with the same average diffusivity. Furthermore, this result is robust regardless of the particular choice of the diffusion coefficient profile. We also discuss how this picture changes beyond the Stratonovich framework. Particularly, the Itô case, which is nonanticipative, can promote the opposite effect, while Hänggi-Klimontovich interpretation reinforces the reduction effect.

First passage of an active particle in the presence of passive crowders

We experimentally study the stochastic transport of a self-propelled camphor boat, driven by Marangoni forces, through a crowd of passive paper discs floating on water. We analyze the statistics of the first passage times of the active particle to travel from the center of a circular container to its boundary. While the mean times rise monotonically as a function of the covered area fraction φ of the passive paper discs, their fluctuations show a non-monotonic behavior - being higher at low and high value of φ compared to intermediate values. The reason is traced to an interplay of two distinct sources of fluctuations - one intrinsic to the dynamics, while the other due to the crowding.

Chemotactic behaviour of Escherichia coli at high cell density

At high cell density, swimming bacteria exhibit collective motility patterns, self-organized through physical interactions of a however still debated nature. Although high-density behaviours are frequent in natural situations, it remained unknown how collective motion affects chemotaxis, the main physiological function of motility, which enables bacteria to follow environmental gradients in their habitats. Here, we systematically investigate this question in the model organism Escherichia coli, varying cell density, cell length, and suspension confinement. The characteristics of the collective motion indicate that hydrodynamic interactions between swimmers made the primary contribution to its emergence. We observe that the chemotactic drift is moderately enhanced at intermediate cell densities, peaks, and is then strongly suppressed at higher densities. Numerical simulations reveal that this suppression occurs because the collective motion disturbs the choreography necessary for chemotactic sensing. We suggest that this physical hindrance imposes a fundamental constraint on high-density behaviours of motile bacteria, including swarming and the formation of multicellular aggregates and biofilms.

Enhanced propagation of motile bacteria on surfaces due to forward scattering

How motile bacteria move near a surface is a problem of fundamental biophysical interest and is key to the emergence of several phenomena of biological, ecological and medical relevance, including biofilm formation. Solid boundaries can strongly influence a cell's propulsion mechanism, thus leading many flagellated bacteria to describe long circular trajectories stably entrapped by the surface. Experimental studies on near-surface bacterial motility have, however, neglected the fact that real environments have typical microstructures varying on the scale of the cells' motion. Here, we show that micro-obstacles influence the propagation of peritrichously flagellated bacteria on a flat surface in a non-monotonic way. Instead of hindering it, an optimal, relatively low obstacle density can significantly enhance cells' propagation on surfaces due to individual forward-scattering events. This finding provides insight on the emerging dynamics of chiral active matter in complex environments and inspires possible routes to control microbial ecology in natural habitats.

Single-species fragmentation: The role of density-dependent feedback

Internal feedback is commonly present in biological populations and can play a crucial role in the emergence of collective behavior. To describe the temporal evolution of the distribution of a single-species population, we consider a generalization of the Fisher-KPP equation. This equation includes the elementary processes of random motion, reproduction, and, importantly, nonlocal interspecific competition, which introduces a spatial scale of interaction. In addition, we take into account feedback mechanisms in diffusion and growth processes, mimicked by power-law density dependencies. This feedback includes, for instance, anomalous diffusion, reaction to overcrowding or to the rarefaction of the population, as well as Allee-like effects. We show that, depending on the kind of feedback that takes place, the population can self-organize splitting into disconnected subpopulations, in the absence of external constraints. Through extensive numerical simulations, we investigate the temporal evolution and the characteristics of the stationary population distribution in the one-dimensional case. We discuss the crucial role that density-dependence has on pattern formation, particularly on fragmentation, which can bring important consequences to processes such as epidemic spread and speciation.

Dynamics of Active Filaments in Porous Media

The motion of active polymers in a two-dimensional porous medium is shown to depend critically on flexibility, activity, and degree of polymerization. For a given Péclet number, we observe a transition from localization to diffusion as the stiffness of the chains is increased. Whereas stiff chains move almost unhindered through the porous medium, flexible ones spiral and get stuck. Their motion can be accounted for by the model of a continuous time random walk with a renewal process corresponding to unspiraling. The waiting time distribution is shown to develop heavy tails for decreasing stiffness, resulting in subdiffusive and ultimately caged behavior.

Efficiency of the flagellar propulsion of Escherichia coli in confined microfluidic geometries

Bacterial movement in confined spaces is routinely encountered either in a natural environment or in artificial structures. Consequently, the ability to understand and predict the behavior of motile bacterial cells in confined geometries is essential to many applications, spanning from the more classical, such as the management complex microbial networks involved in diseases, biomanufacturing, mining, and environment, to the more recent, such as single cell DNA sequencing and computation with biological agents. Fortunately, the development of this understanding can be helped by the decades-long advances in semiconductor microfabrication, which allow the design and the construction of complex confining structures used as test beds for the study of bacterial motility. To this end, here we use microfabricated channels with varying sizes to study the interaction of Escherichia coli with solid confining spaces. It is shown that an optimal channel size exists for which the hydrostatic potential allows the most efficient movement of the cells. The improved understanding of how bacteria move will result in the ability to design better microfluidic structures based on their interaction with bacterial movement.

Ideal circle microswimmers in crowded media

Microswimmers are exposed in nature to crowded environments and their transport properties depend in a subtle way on the interaction with obstacles. Here, we investigate a model for a single ideal circle swimmer exploring a two-dimensional disordered array of impenetrable obstacles. The microswimmer moves on circular orbits in the freely accessible space and follows the surface of an obstacle for a certain time upon collision. Depending on the obstacle density and the radius of the circular orbits, the microswimmer displays either long-range transport or is localized in a finite region. We show that there are transitions from two localized states to a diffusive state each driven by an underlying static percolation transition. We determine the non-equilibrium state diagram and calculate the mean-square displacements and diffusivities by computer simulations. Close to the transition lines transport becomes subdiffusive which is rationalized as a dynamic critical phenomenon.

Confined Flow: Consequences and Implications for Bacteria and Biofilms

Bacteria overwhelmingly live in geometrically confined habitats that feature small pores or cavities, narrow channels, or nearby interfaces. Fluid flows through these confined habitats are ubiquitous in both natural and artificial environments colonized by bacteria. Moreover, these flows occur on time and length scales comparable to those associated with motility of bacteria and with the formation and growth of biofilms, which are surface-associated communities that house the vast majority of bacteria to protect them from host and environmental stresses. This review describes the emerging understanding of how flow near surfaces and within channels and pores alters physical processes that control how bacteria disperse, attach to surfaces, and form biofilms. This understanding will inform the development and deployment of technologies for drug delivery, water treatment, and antifouling coatings and guide the structuring of bacterial consortia for production of chemicals and pharmaceuticals.

Nonlinear population dynamics in a bounded habitat

A key issue in ecology is whether a population will survive long term or go extinct. This is the question we address in this paper for a population in a bounded habitat. We will restrict our study to the case of a single species in a one-dimensional habitat of length L. The evolution of the population density distribution ρ(x, t), where x is the position and t the time, is governed by elementary processes such as growth and dispersal, which, in standard models, are typically described by a constant per capita growth rate and normal diffusion, respectively. However, feedbacks in the regulatory mechanisms and external factors can produce density-dependent rates. Therefore, we consider a generalization of the standard evolution equation, which, after dimensional scaling and assuming large carrying capacity, becomes ∂_tρ=∂_x(ρ^ν-1∂_xρ)+ρ^μ, where μ,ν∈R. This equation is complemented by absorbing boundaries, mimicking adverse conditions outside the habitat. For this nonlinear problem, we obtain, analytically, exact expressions of the critical habitat size L_c for population survival, as a function of the exponents and initial conditions. We find that depending on the values of the exponents (ν, μ), population survival can occur for either L ≥ L_c, L ≤ L_c or for any L. This generalizes the usual statement that L_c represents the minimum habitat size. In addition, nonlinearities introduce dependence on the initial conditions, affecting L_c.