Invariance properties of bacterial random walks in complex structures

TL;DR

This paper demonstrates that the mean residence time of swimming bacteria in complex microstructures can be predicted by a generalized invariance property of random walks, linking geometry to bacterial behavior.

Contribution

It introduces a generalized invariance property for bacterial random walks that predicts mean residence times based solely on free volume to surface ratio.

Findings

Mean residence time is constrained by free volume to surface ratio.

Structural disorder affects path length distributions but not mean residence time.

Predicts bacterial colonization behavior based on geometric properties.

Abstract

Motile cells often explore natural environments characterized by a high degree of structural complexity. Moreover cell motility is also intrinsically noisy due to spontaneous random reorientation and speed fluctuations. This interplay of internal and external noise sources gives rise to a complex dynamical behavior that can be strongly sensitive to details and hard to model quantitatively. In striking contrast to this general picture we show that the mean residence time of swimming bacteria inside artificial complex microstructures, can be quantitatively predicted by a generalization of a recently discovered invariance property of random walks. We find that variations in geometry and structural disorder have a dramatic effect on the distributions of path length while mean values are strictly constrained by the sole free volume to surface ratio. Biological implications include the possibility of predicting and controlling the colonization of complex natural environments using only geometric informations.

Figures

Equations

$$t=T\frac{S}{\Sigma }$$

$$t=\tau \frac{\rho v}{\pi }PT\frac{1}{N}$$

$$\ell =v\tau =\pi \frac{S}{P}$$

$$\ell =\pi \frac{S}{P}\left(1-\frac{Ns}{S}\right)$$

$$\tau =⟨L/v⟩=⟨L⟩⟨v^{-1}⟩=\ell /\bar{v}$$

Full text

Invariance properties of bacterial random walks in complex structures

Giacomo Frangipane

Dipartimento di Fisica, Sapienza - Università di Roma, Roma, I-00185, Italy

NANOTEC-CNR, Institute of Nanotechnology, Soft and Living Matter Laboratory, Roma, Italy

Gaszton Vizsnyiczai

Dipartimento di Fisica, Sapienza - Università di Roma, Roma, I-00185, Italy

Claudio Maggi

NANOTEC-CNR, Institute of Nanotechnology, Soft and Living Matter Laboratory, Roma, Italy

Romolo Savo

Laboratoire Kastler Brossel, École Normale Supérieure-Paris Sciences et Lettres (PSL) Research University, CNRS, Université Pierre et Marie Curie–Sorbonne Universités, Collège de France, 75005 Paris, France

Optical Nanomaterial Group, Institute for Quantum Electronics, Department of Physics, ETH Zurich, Auguste Piccard Hof 1, 8093, Zurich, Switzerland

Alfredo Sciortino

Dipartimento di Fisica, Sapienza - Università di Roma, Roma, I-00185, Italy

Sylvain Gigan

Laboratoire Kastler Brossel, École Normale Supérieure-Paris Sciences et Lettres (PSL) Research University, CNRS, Université Pierre et Marie Curie–Sorbonne Universités, Collège de France, 75005 Paris, France

Roberto Di Leonardo

Dipartimento di Fisica, Sapienza - Università di Roma, Roma, I-00185, Italy

NANOTEC-CNR, Institute of Nanotechnology, Soft and Living Matter Laboratory, Roma, Italy

roberto.dileonardo@uniroma1.it

**Motile cells often explore natural environments characterized by a high degree of structural complexity [1, 2]. Moreover cell motility is also intrinsically noisy due to spontaneous random reorientation and speed fluctuations [3]. This interplay of internal and external noise sources gives rise to a complex dynamical behavior that can be strongly sensitive to details [4, 5, 6] and hard to model quantitatively. In striking contrast to this general picture we show that the mean residence time of swimming bacteria inside artificial complex microstructures, can be quantitatively predicted by a generalization of a recently discovered invariance property of random walks [7]. We find that variations in geometry and structural disorder have a dramatic effect on the distributions of path length while mean values are strictly constrained by the sole free volume to surface ratio. Biological implications include the possibility of predicting and controlling the colonization of complex natural environments using only geometric informations. **

Microfabrication provides an ideal tool to investigate the dynamics of active particles in artificial complex environments having a precise and tunable internal structure [8]. Arrays of scattering obstacles have been used to demonstrate the possibility of rectification and sorting in self-propelled systems [9, 10]. In all these cases small details in dynamical behavior can have significant quantitative consequences. For example, the timescale of spontaneous reorientation or the distribution of scattering angles produced by obstacles have a dramatic effect on the rectification efficiency of microstructures  [4, 5, 6]. In striking contrast, a recently discovered invariance property of random walks implies that the average path length inside a closed domain of arbitrary shape is only proportional to the volume to surface ratio with a numerical prefactor that solely depends on the spatial dimensions. Despite its wide generality, this invariance property has started to find applications only very recently, mainly in the field of light propagation in turbid media [11, 12]. Here we study motile bacteria exploring artificial microstructures with arbitrary geometries and a wide range of obstacle densities. Using the invariance property, we quantitatively predict the mean length of the paths traced by bacteria inside these microstructures. Assuming that the mean swimming speeds are not significantly perturbed by the presence of obstacles, we further predict and experimentally verify the invariance of the mean residence time, a more biologically relevant quantity. As a result, for all obstacle densities, the mean residence time is only determined by the microstructure volume to surface ratio and the mean inverse speed.

We use direct laser writing by two-photon polymerization [13, 14] to build 3D microchambers consisting of a thick “roof” supported by randomly distributed pillars fabricated on a microscope coverslip. A 1.4 $\upmu$m-thick gap between the coverslip and the roof confines bacteria in a quasi-2D geometry. The entire path traced by the cells inside the microchambers is therefore restricted on the focal plane and can be easily tracked by digital video microscopy (see Supplementary Video 1). Moreover, the presence of two close surfaces prevents circular swimming caused by hydrodynamic coupling to a single solid wall [15, 16]. In all experiments we use a smooth-swimming strain of E. coli bacteria expressing a red fluorescent protein. Microstructures are immersed in a low-density cell suspension ($\sim 3\times 10^{6}$ cells $\mathrm{ml}^{-1}$, see Methods) that we track by epifluorescence videomicroscopy. The pillars play the role of obstacles and by varying their number we can tune the degree of internal complexity of the structures. The obstacles are fabricated as randomly distributed and non-overlapping short segments measuring $3.7\leavevmode\nobreak\ \upmu$m in length and with a thickness of $0.7\leavevmode\nobreak\ \upmu$m (see Supplementary Figure 1). Obstacles with such an elongated shape significantly deflect colliding bacteria and can be densely packed inside microstructures.

We fabricate a total number of 11 circular structures with a radius $R=31\leavevmode\nobreak\ \upmu$m and with a varying number of obstacles $N=7,15,25,36,46,55,64,76,86,109,119$. When $N$ is small, bacteria traverse the structures with approximately straight trajectories (see Fig. 1(a)). Since the flux of bacteria through a line element of the boundary is proportional to the cosine of the angle between the cell swimming direction and the local boundary normal, the distribution of path lengths $L$ for straight trajectories is given by the chord-length distribution ${p(L)=L/(2R\sqrt{4R^{2}-L^{2}})}$. As shown in Fig. 1(b), this distribution increases monotonically diverging at the maximum value $L=2R$. The experimental distribution of path lengths, obtained by tracking approximatively 1000 bacteria in a structure with $N=7$ obstacles, closely follows the theoretical distribution with small deviations due to rotational diffusion of bacteria and to scattering by the few isolated obstacles. When the number of obstacles is increased, bacteria are more likely to be back-scattered by obstacles close to the domain boundary giving rise to an excess of short paths that appears as a growing peak at small path lengths. At the same time, dense obstacles produce multiple deflections that result in trajectories with a total length that can exceed the circle diameter producing a tail in the distribution at $L>2R$.

Despite the strong dependence on $N$ of the path length distribution, the invariance theorem by Blanco and Fournier [7] (BFT) can be used to obtain a strikingly simple formula for the mean path length that quantitatively matches observed values for all $N$. To this aim we now illustrate an alternative derivation of the BFT result by means of a very simple argument. Let’s imagine to observe for a total time $T$ the motion of $N$ non-interacting random walkers moving with constant speed $v$ over a region of space of total area $\Sigma$. If reorientation dynamics is isotropic, the walkers will fill the region with a uniform concentration $\rho=N/\Sigma$ and an everywhere isotropic velocity distribution. The time $t$ spent by each walker inside a generic subregion of area $S$ will be given by:

[TABLE]

Alternatively we can obtain $t$ as the mean residence time $\tau$ multiplied by the number of times each walker is expected to enter the subregion during the observation time $T$. This number can be obtained by counting the number of bacteria that enter the subregion in a total time $T$ and then divide by $N$. Since, for an isotropic velocity distribution, the inward flux of walkers is $\rho v/\pi$ over the entire perimeter $P$ of the sub-region, we have:

[TABLE]

Equating (1) and (2) we obtain:

[TABLE]

When the subregion contains rigid obstacles, the above derivation remains valid as long as their presence does not significantly perturb the homogeneity of the density field in the accessible space and the isotropy of random walks along the domain boundary [17]. This condition is fulfilled for small and convex obstacles producing scattering events lasting for a short interaction time. In this regime, however, the total area $S$ in Eq. (3) must be replaced with the total accessible area $S^{\prime}=S-Ns$ where $s$ is the excluded area due to a single obstacle and $Ns$ the total excluded area. Therefore the mean path length $\ell$ is expected to decrease linearly with the number of obstacles according to the law:

[TABLE]

In Fig. 1(c) we report as solid circles the experimental values of the mean path lengths as a function of the number of obstacles. Experimental data closely match the predicted linear behaviour in (Eq. 4) with an intercept fixed at $\pi S/P=\pi R/2=49.3\;\upmu$m and a value for ${s=6.3\leavevmode\nobreak\ \upmu\mathrm{m}^{2}}$ that is independently estimated by dilating the physical boundary of the obstacle by half the cell width ($0.4\leavevmode\nobreak\ \upmu$m [18]) (see Supplementary Figure 1).

There are many situations in which the residence time is more relevant then the path length. For example, if adhesion to surfaces is an activated process [19] the probability that a cell will remain stuck inside a structure will be proportional to the total time spent inside it. The population average residence time can be related to the mean path length under the reasonable assumption that cell speeds $v$ are uncorrelated to path lengths $L$ so that:

[TABLE]

with $\bar{v}=1/\langle v^{-1}\rangle$. Although speeds can be broadly distributed in a cell population we can again use the invariance property to quantitatively predict mean residence times once the characteristic speed $\bar{v}$ is known. The measured mean residence time as a function of $N$ is plotted in Fig. 1(c) showing again a close agreement with the theoretical prediction in Eq. (5) where the parameter $\bar{v}$ has been independently measured from experimental trajectories.

We have demonstrated that, for a fixed shape, the density of obstacles has a strong effect on the distribution of path lengths while leaving the mean value unchanged as long as the total internal accessible area remains the same. We will now investigate the effects of shape and size on the path length distributions and again conclude that, although path distributions are very sensitive to the detailed geometry of the domain, mean path lengths and residence times can be accurately predicted using only the surface to perimeter ratio. We fabricate microstructures of circular, squared and triangular shapes all in four different sizes as shown in Fig. 2(a). In all structures the obstacle density is fixed to ${N/S=0.016\leavevmode\nobreak\ \upmu\mathrm{m}^{-2}}$. All recorded trajectories are also plotted in Fig. 2(a) with a color map encoding the scaled path length $L/\ell$. In Fig. 2(b) we plot for each structure the histogram of path lengths normalized by the expected value $\ell$ in Eq. (4). As we move from circles to squares and then to triangles we observe that the corresponding histograms are characterized by larger frequencies for both short and long path lengths. The probability of short paths increases as we introduce convex corners and reduce their angle, as evidenced by the blue-violet color of the corners in squared and triangular structures. The presence of longer paths in triangles can be also understood noting that the normalized maximum straight paths for circles, square and triangles are respectively $4/\pi,4\sqrt{2}/\pi$ and $4\sqrt{3}/\pi$ (see Supplementary Figure 2). Despite these qualitative and quantitative differences between histograms, mean path lengths are always close to 1 when normalized by the predicted value $\ell$ as shown in Fig 2(b). The observed path length distributions can be very well reproduced by a Lorentz gas model where collisions with obstacles are described as instantaneous reorientation events [20]. These events are assumed to occur with a constant probability per unit path length given by $\rho\,\sigma$, where $\sigma=2b/\pi$ is the rotational average cross section of line obstacles with length $b$ (see Methods). For perfectly straight trajectories, circles and squares have path length distributions with a vertical asymptote. A strong peak is found at the location of this divergence and is progressively smeared as the linear size of the domain becomes larger than the mean free path ${\lambda=(\rho\sigma)^{-1}}$. Upon increasing the size of the region a peak around $L=0$ appears due to bacteria that scatter close to the boundary, reorient and exit. The mix of these two types of trajectories gives rise to the observed two-peak structure. Differently for the triangle no asymptote is present in the $P(L)$ for large $\lambda$, therefore only a peak at $L=0$ progressively appears. These results show that $P(L)$ depends strongly on the random walk properties (as the mean free path), and on the domain features (shape and size). Differently the mean values $\ell$ and $\tau$ depend only on the surface to perimeter ratio of the domain as expressed by Eq. (4). Experimental values of $\ell$ are shown Fig. 3(a) for all shapes and sizes. We find an excellent agreement with theory using the same value ${s=6.3\leavevmode\nobreak\ \upmu\mathrm{m}^{2}}$ that was used previously to predict the experimental values of $\ell$ as a function of the number of obstacles (see Fig. 1(c)). This value of $s$ has also been used to calculate $\tau$ according to Eq. (5) which captures quantitatively well the experimental data as shown in Fig. 3(b).

In conclusion, we have shown that bacterial random walks inside complex and crowded microstructures do obey a general invariance property that constraints the mean value of the path length to a simple and predictable value uniquely determined by the volume to surface ratio of the structure. More interestingly the mean residence time is also invariant as long as the volume to surface ratio remains constant. In a counterintuitive way, a greater number of obstacles does not increase the average residence time, but shortens it by decreasing the total accessible volume. The robustness of our results opens the way to different application perspectives as for example using bacteria to gather information about the internal structure of natural microenvironments. Alternatively, any significant deviation of measured residence times from the predicted value may signal the presence of physical (e.g. adhesion) and biological (e.g. taxis) mechanisms breaking the uniformity and isotropy of random walks.

Methods

Microfabrication

The microchamber structures were fabricated with a custom built two-photon polymerization setup [21] from SU-8 photoresist (MicroChem Corp). A single fabrication focus was used with 8 mW laser power and 100 $\upmu\mathrm{m\;s}^{-1}$ scanning speed. After the laser fabrication scan the SU-8 photoresist sample was baked at 100 *∘*C for 7 minutes, then developed by its standard developer solvent and finally rinsed in a 1:1 mixture of water and ethanol. Strong adhesion of the microchamber structures to the carrier coverglass was ensured by a layer of OmniCoat adhesion promoter (MicroChem Corp).

Microscopy

Epifluorescence imaging were performed on an inverted optical microscope (Nikon TE-2000U) equipped with a 60$\times$ (NA=1.27) water immersion objective and a high-sensitivity CMOS camera (Hamamatsu Orca Flash 4.0).

Cell cultures

The E. coli strain used is the HCB437[22], a smoooth-swimmer, transformed by a plasmid expressing the red fluorescent protein mRFP1 under the control of the lacI inducible promoter (BioBricks, BBa_J04450 coding device inserted in pSB1C3 plasmid backbone, iGEM Catalogue). Cells were grown overnight in 10 mL of LB supplemented with kanamycin (30 $\upmu\mathrm{g\;mL^{-1}}$) and chloramphenicol (20 $\upmu\,\mathrm{g\;mL^{-1}}$) in a shaking incubator at 33*∘*C and 200 rpm. In the morning 50 $\upmu$L of the saturated culture was diluted (1:100) into 5 mL of tryptone broth supplemented with antibiotics and incubated in the same condition for 4 hours (until an OD${}_{590}\sim$ 0.2). The production of mRFP1 was induced by adding 1 mM IPTG, and the cultures were grown for 2 further hours until exponential phase (OD${}_{590}\sim$ 0.6-0.8). Cells were collected and washed three times by centrifugation (1500 rcf, 5’) in motility buffer (MB: 0.1 mM EDTA, 10 mM phosphate buffer, 1 mM glucose and 0.02%Tween20).

Simulations

We simulate bacterial trajectories as straight runs intercalated by random reorientations (run and tumble dynamics) that mimic collisions with obstacles. Run lengths are exponentially distributed with a mean free path ${\lambda=(\rho\sigma)^{-1}}$. The obstacle density $\rho$ is fixed to the experimental value ${\rho=0.016\leavevmode\nobreak\ \upmu\mathrm{m}^{-2}}$. The obstacle cross section $\sigma=2.9\leavevmode\nobreak\ \upmu\mathrm{m}$ is evaluated as the average cross section of a randomly oriented line with length $b=4.5\leavevmode\nobreak\ \upmu$m, i.e. ${\sigma=\pi^{-1}\int d\theta\,b\cos(\theta)=2b/\pi}$. Assuming that an elongated particle fully aligns with a line-shaped obstacle, collisions will produce angular deflections $\delta$ that are only determined by the obstacle orientation and restricted to the interval $[-\pi/2,\pi/2]$. Weighting each outgoing direction with the cross section of the correspondingly tilted obstacle we obtain the distribution of deflection angles $p(\delta)=\sin|\delta|$ for $|\delta|\leq\pi/2$ and zero elsewhere. We also account for a small reduction in $\delta$ that could be caused by hydrodynamic interactions, causing particles to curve around the obstacle edges [23], or by the fact that reorientation during collision takes a finite time [16], while the obstacle has finite length. The total deflection angle is then set to $\delta^{\prime}=\delta-\mathrm{sign}(\delta)\,\alpha$, where $\alpha$ is set to $\alpha=0.39\,(\approx\pi/8)$ to fit the histograms in Fig. 2(b).

Code availability

The codes that support the findings of this study are available from the corresponding authors upon reasonable request.

Acknowledgements

The research leading to these results has received funding from the European Research Council under the European Union’s Seventh Framework Programme (FP7/ 2007–2013)/ERC Grant Agreement No. 307940. S.G. acknowledges support from Institut Universitaire de France. R.S. acknowledges funding from the European Union’s Horizon 2020 research and innovation programme under the Marie Sklodowska-Curie grant agreement No. 800487.

Data Availability

The data that support the plots within this paper and other findings of this study are available from the corresponding authors upon reasonable request.

Additional information

The authors declare no competing financial interests.

The corresponding author is responsible for submitting a competing financial interests statement on behalf of all authors of the paper. This statement must be included in the submitted article file.