A Pili-Driven Bacterial Turbine

Work generated by self-propelled bacteria can be harnessed with the help of microdevices. Such nanofabricated microdevices, immersed in a bacterial bath, may exhibit unidirectional rotational or translational motion. Swimming bacteria that propel with the help of actively rotating flagella are a prototypical example of active agents that can power such microdevices. In this work, we propose a computational model of a micron-sized turbine powered by bacteria that rely on active type IV pili appendages for surface-associated motility. We find that the turbine can rotate persistently over a time scale that significantly exceeds the characteristic times of the single cell motility. The persistent rotation is explained by the collective dynamics of multiple pili of groups of cells attaching to and pulling on turbine. Furthermore, we show that the turbine can rotate permanently in the same direction by altering the pili binding to the turbine surface in an asymmetric fashion. We thus can show that by changing the adhesive properties of the turbine while keeping its symmetric geometry, we can still break the symmetry of its rotation. Altogether, this study widely expands the range of bacteria that can be used to power nanofabricated microdevices, and, due to high pili forces generated by pili retraction, promises to push the harnessed work by several orders of magnitude.


INTRODUCTION
The last years have seen significant advances in nanofabrication, permitting the invention of a wide range of micron-sized artificial devices. A fascinating question is how such devices can be powered by actively moving biological matter, typically consisting of bacteria. Examples are beads that move due to collision or attachment of cells to its surface [1,2], swimming devices due to individual swimming cells and bacterial carpets attached to their surface [3] or confined in the microdevice structure [4,5], and geometrically asymmetric devices immersed in a bath of actively moving cells that rotate due to random collisions with the cells [6][7][8][9][10]. Usually, bacteria that swim with the help of rotating flagella were employed. In this study, we propose to consider twitching bacteria that exhibit surfaceassociated locomotion mediated by type IV pili [11,12]. Pili are microns long active polymeric appendages protruding from the cell membrane. They undergo cycles of protrusion and retraction and can bind to a substrate and the pili of other cells. The combination of these two processes leads to aggregation of cells [13][14][15] and twitching motility on a substrate [16][17][18][19] by a mechanism reminiscent of a grappling hook (see Figure 1A).
Here, we investigate how bacteria that can attach type IV pili to a micron-sized turbine can drive its rotation. This system is especially interesting since multiple pili belonging to an individual bacterium can generate total forces in the range of nano-Newtons [20]. In contrast to that, bacteria that swim with the help of flagella typically have thrust forces in the order of 1-10 pN [21]. The molecular motor involved in the disassembly an hence, in the retraction of an individual pilus, called pilT, can generate forces in the range of 100 -180 pN [22]. This makes pilT the strongest molecular motor known in nature, with forces up to 20 times larger than those generated by kinesins or polymerases [20]. Hence, we are asking if cells with type IV pili that can create so large forces might be more attractive candidates to power microdevices.
In this work, by means of a computational model, we study the dynamics of a rotatable turbine immersed in a bath of twitching bacteria. We investigate how the turbine rotation is affected by the binding and unbinding of the bacterial pili. We find that, due to the adhesion of multiple pili to the turbine over time scales that can strongly exceed the characteristic time scales of the individual pili attachment, the turbine can persistently rotate in one direction over extensive durations. While these persistent rotations have a limited lifetime, one can engineer a system where they become permanently unidirectional by introducing asymmetric binding of the pili to the turbine.

MATERIALS AND METHODS
First, we introduce the computational model of twitching bacteria and their interactions with a turbine. A related version of this computational model was considered previously to describe bacterial surface motility [19] and bacterial aggregates [23][24][25].
While we consider the bacterium Neisseria gonorrhoeae as the primary biological example, the computational model can be easily adapted to account for other bacteria that use type IV pili, e.g., Neisseria meningitidis [14] or Pseudomonas aeruginosa [26]. In the following, we focus on the regime of low cell density. This allows us not to consider cell-cell interactions and the formation of bacterial aggregates mediated by the binding of the pili of different cells. As a result, we do not expect three dimensional aggregates to form [23,24]. This also enables us to only consider a simpler two-dimensional system.

Geometry of the Cells, Pili and Turbine
The experimental system we are mimicking are bacterial cells confined in a box with the turbine at its center (see Figure 1B). The cells can move over the substrate via pili binding and unbinding and additionally interact with the turbine by pili attachment and excluded volume effects.
While Neisseria gonorrhoeae cells typically have a diplococcus shape [27], we simplify the in silico cell shape of bacteria as twodimensional circular disks with a radius R. Please note, however, that a diplococcus shape can be also considered [19,[23][24][25]. Each cell possesses exactly N p pili that are homogeneously distributed on the cell outline (see Figure 1A). We approximate pili as straight lines connecting two points: their start point (also called anchor point), located at the circular cell surface, and their end FIGURE 1 | Overview of the computational model of bacteria driving a rotatable turbine by attachment of type IV pili. (A) Sketch of how pili binding and retraction can lead to a force acting on a cell, a process reminiscent of a grappling hook. (B) Overview of the system. Cells are located in a box together with the rotatable turbine. The cells possess pili which they use for the motility on the substrate and for attachment to the turbine. The turbine center is located at r (turb) and has an orientation defined by the angle α. (C) Summary of all forces acting on the cells and turbine. Friction forces with the environment (surrounding liquid and substrate) are not shown. point. The distance between these two points is called the contour length l c of the pilus.
The turbine is described by N turb × 2 disks (with radius R) arranged in a double row array (see Figure 1B). Since the only way how the turbine can move is by rotation, the relative locations of the turbine disks towards each other are fixed. The orientation of the turbine is described by the angle α with respect to the y-axis and the turbine center is located at the position r (turb) (see Figure 1B).

Pili Dynamics and Binding Properties
Initially, a pilus protrudes with the velocity v p in a direction perpendicular to the cell surface. When a pilus reaches a specific length, drawn from an exponential length distribution with mean length l p [17,18], it starts to retract with a velocity v (0) r . This retraction continues until the pilus has a contour length l c = 0. In that case, the pilus is removed and immediately, a new pilus protrudes from the cell from the same position.
A pilus binds stochastically to the substrate or the turbine disks, independently whether it is protruding or retracting. In both cases, a pilus can only bind with its tip and it can bind only to either the substrate or the turbine disks. After binding, a pilus immediately starts to retract [28]. The binding is modeled as a Poisson process with the binding rate κ (sub) att to the substrate and κ (turb) att to any of the turbine disks. Since the model is twodimensional, we ignore that, for pili to bind to a substrate, they first need to be long enough to reach the substrate with their tips. In 3D, this leads to a delay in the initial binding event of a newly protruding pilus, but is less relevant for following attachments of pili. Hence, initially, the binding rate might be smaller to allow a pilus to protrude to reach the substrate. To account for this process, pili that bind to the substrate the first time do so with a rate κ (sub,0) att v p /l p . Again, the binding is modelled as a Poisson process. The rate κ (sub,0) att corresponds to the binding of a pilus that protrude with the velocity v p until it reaches a length determined by an exponential distribution with mean length l p , in which case it binds. On average, this takes the time 1/κ (sub,0) att . An attached pilus will generate a pulling force that acts on the cell and, if attached to a turbine disk, also on the turbine. Each pilus is modeled as a Hookean spring with the spring constant k p . After attachment, a pilus is stretched due to its retraction and hence, mediates a pulling force. This force is proportional to the difference between the contour length l c and the length of the pilus that it would have if it was not attached, here called the free length l f . Here, we consider the case where a pilus can only generate a pulling force and no pushing force, thus the force is zero if l c ≤ l f . Experimentally, it has been shown that the pilus force affects the retraction velocity by v r max 0, v (0) with the stalling force F s [22]. The force of a pilus also affects the unbinding from the substrate. Pilus detachment is modelled by a Poisson process with the rate for detachment from the substrate, as motivated by [18]. Here, we introduce the detachment times τ (sub,1) det and τ (sub,2) det and the detachment forces F (sub,1) det and F (sub,2) det . For detachment from a turbine disk, we describe the rate by a simpler relation with the detachment time τ (turb) det and detachment force F (turb) det . We picked a different expression for the pili-turbine detachment (in comparison to the pili-substrate detachment) because in general, the turbine does not necessarily have to consist of the same material as the substrate, which determines pili binding properties [16,17]. Below, we will explore turbine rotation in dependence on pili-turbine interactions, thus favouring a simpler single exponential form.

Cell Forces and Motility
We model pili as Hookean springs with the spring constant k p . where the pilus tip is attached (see Figure 1C).
Next to the active pili force, passive excluded volume forces are acting on the cells. Cells are located in a two-dimensional box with size L × L and a cell i that overlaps with the boundary wall is exposed to a repulsive force F (wall) i (see Figure 1C). This force is modelled as a harmonic interaction with spring constant k (wall) and acts on the cell center in the normal direction of the boundary wall if the overlap is smaller than the cell radius R. Additionally, intersections of 2 cells i and k lead to a repulsive force F (c−c) ik of the centers of both cells with the spring constant k (c−c) (see Figure 1C). A similar type of repulsive force F (turb m) i is acting between the cells and the turbine disks m with the spring constant k (turb) (see Figure 1C).
The total force of the cell i is given by Additionally, the total torque acting on a cell is In the overdamped limit [29], a force mediates a translational motion of the cell with the velocity and the torque leads to a rotation with the angular velocity vector Frontiers in Physics | www.frontiersin.org May 2022 | Volume 10 | Article 875687 Here, we introduce the translational mobility μ t and the rotational mobility μ r . The same forces and torques cause an equivalent displacement of the pili anchor points r

Turbine Torque and Rotation
The turbine is only able to undergo rotational motion. The total torque acting on the turbine is given by The rotation of the turbine is modelled in the overdamped limit with a mobility μ turb . In that case, the turbine angular velocity vector is

Parameters and Details of Numerical Solution
The simulations were performed on the local computing cluster of the Max Planck Institute for the Physics of Complex Systems (Dresden, Germany), consisting of x86-64 GNU/Linux systems. The code was written in C++. We use an Euler algorithm to solve the equations of motion with a time step δt. While this is one of the most simple numerical schemes to solve the equations of our computational model, it often leads to numerical errors and instabilities when modelling molecular dynamics systems for long times [30]. We do not expect that this is a problem in our system due to the stochastic nature of the pili binding and unbinding, which basically represents our system as a series of many short time events, continuously interrupted by rearrangements in the pili network. Hence, we do not expect any differences in the qualitative outcome of the simulations.
If not stated otherwise, we use the parameters provided in Table 1. Most parameters we use are based on previous studies. The excluded volume spring constants k (c−c) , k (turb) and k (wall) have no effect on the simulation outcome as long as they are chosen large enough to be able to compete with the pili forces. For the remaining parameters, e.g., the turbine mobility, we do not expect a qualitative difference in the results of the simulation.
We initialize the simulation by randomly distributing the cells in the box and only analyse the turbine rotation after an initialization period of 1,000 s. To calculate the angular velocity of the turbine rotation, we compute with Δt = 0.5 s and the turbine orientation α (see Figure 1B).

RESULTS
In the following, we first demonstrate that bacteria binding to the turbine with the help of pili can cause a persistent turbine rotation over time scales that exceed the pili detachment times significantly. Next, we study how the persistent turbine rotation depends on the binding properties of the pili to the turbine and the number of cells in the system. This allows us to unravel the underlying mechanism that causes the persistent rotation of the turbine. Finally, we propose a system where an asymmetric binding of the pili to the turbine causes a permanent unidirectional rotation of the turbine.

Adhesion of Motile Bacteria Drives Turbine Rotation
To investigate how the adhesion of cells to a turbine affects the turbine rotation, we first simulate the computational model for three different cases: 1) pili detach with detachment times of τ (turb) det 50 s (see Figure 2A and Supplementary Movie S1), motivated by binding times in the order of a few minutes inferred from cell trajectories on plastic surfaces [17], 2) pili detach with τ (turb) det 2 s (see Supplementary Movie S2), motivated by considerably smaller binding times experimentally measured for BSA coated beads [18] and 3) pili not bind to the turbine at all (see Supplementary Movie S3) and not the detachment force F (turb) det . We expect that there will be no qualitative difference in the turbine rotation for both cases, an increase (decrease) of τ (turb) det and a increase (decrease) of F (turb) det , pili will bind more (less) strongly to the turbine. For simplicity, we assume that F (turb) det ≔ F s with the pilus stalling force F s , see Eq. 1. In Figure 2B, we show the dynamics of the turbine angle α for the three cases. We find that for the largest detachment time, τ (turb) det 50 s, the rotation is the strongest. If pili are not permitted to attach to the turbine, the rotation appears to be the weakest. Additionally, the distribution of velocities shows that if pili cannot bind to the turbine, the turbine will move with smaller angular velocities ω, while the distributions seem to not depend on the detachment rate τ (turb) det , as long as pili are capable of binding to the turbine. We also find that the turbine does not undergo extended phases in which it is not moving, but instead seems to be continuously rotating (see Supplementary Figure S1).
To quantify how strong the rotation is, we compute the angular mean squared displacement (angular MSD), given by (see Figure 2C). Indeed, the angular MSD is the highest for τ (turb) det 50 s and the lowest for the case where pili do not attach to turbine disks. For τ (turb) det 2 s, the angular MSD is identical to the case τ (turb) det 50 s for small time differences (Δt). In both cases, τ (turb) det 2 s and non-attaching pili, the angular MSD follows a linear scaling, δ α ∝Δt, corresponding to a diffusive regime. For τ (turb) det 50 s, we find a superdiffusive regime for intermediate time differences Δt ≈ 10, −,10 3 s. For particles undergoing translational motion, such superdiffusive regimes usually emerge from ballistic motion. For the rotatable turbine, this corresponds to a regime where the turbine moves persistently in a random direction. Indeed, such a behaviour is observed in Figure 2B, where trajectories of the angle α move persistently in one direction over periods of hundreds to thousands of seconds before the rotation turns towards the opposite direction. For larger To learn more about how persistent the turbine rotation is, we compute the angular velocity autocorrelation function (angular velocity ACF), given by Independently of whether pili can bind to the turbine or not, we find that the correlation function ](Δt) is decaying with time (see Figure 2D). This confirms our previous observation for the angular MSD where we found that for very large time differences the turbine rotation is diffusive and no longer shows signs of persistence. We also find that for τ (turb) det 50 s, the normalised angular velocity ACF ](Δt)/ ](0) has the slowest decay. This behaviour correlates with the persistent rotation for time differences between Δt ≈ 10, −,10 3 . Surprisingly, we find that the decay of the angular velocity ACF is faster for τ (turb) det 2 s than for the case where pili do not bind to the turbine. We will provide an explanation for this behaviour in Section 3.3. Before doing that, we will have a closer look at how the turbine rotation depends on the cell adhesion to the turbine disks and the number of cells in the system.

Turbine Rotation Is Controlled by Bacterial Adhesion Strength and Cell Number
We begin with a systematic analysis of how the angular MSD δ α (Δt) depends on the detachment time of pili from the turbine disks τ (turb) det , shown in Figure 3A. For small time differences Δt < 10 s, we find that δ α is independent of τ (turb) det and exhibits a linear scaling, δ α ∝Δt, corresponding to a diffusive regime. For high enough values of the detachment time τ (turb) det and Δt > 10 s, the angular MSD becomes superdiffusive, δ α ∝Δt 2 . The duration of this superdiffusive regime increases with increasing τ (turb) det and is barely observable for detachment times τ (turb) det < 10 s. For even longer time differences Δt, the MSD becomes diffusive again, δ α ∝Δt.
Next, we investigate the normalised angular velocity ACF ](Δt)/](0) (see Figure 3B) and find that for Δt < 2 s, it rapidly decreases. For Δt ≥ 2 s, it exponentially decreases with the Δt. By fitting a function of the form ](Δt) ∝ exp(−Δt/τ) to this later regime, we identify the characteristic time τ of the exponential decay of ](Δt) and can investigate its dependence on τ (turb) det . We find that it is increasing with increasing detachment times τ (turb) det . This is in accordance with the increasing duration of the superdiffusive regime of the angular MSD. The stronger pili of a cell bind to the turbine, the more persistently the turbine rotates in one direction before it inverts its direction.
Additionally, we investigated how the rotation of the turbine is affected by the density ρ of cells in the system. We find that with increasing density, the angular MSD δ α becomes smaller (see Figure 3C). While we observe the diffusive regime of δ α (Δt) for Δt < 10 s, a pronounced superdiffusive regime for larger Δt vanishes if the cell density becomes too large. For the normalised angular velocity ACF ](Δt)/](0) (see Figure 3D), the characteristic decay time τ is initially increasing with cell density ρ and then decreases when ρ increases further.
In the following section, we provide qualitative arguments that explain the observed behaviours.

Unbinding Dynamics of Attached Bacteria Explains Characteristic Time of Turbine Rotation
To understand how the persistent rotation of the turbine in an otherwise symmetric system can emerge, we first consider how strongly cells bind to the turbine disks with the help of their pili and how pili binding affects substrate attachment of the remaining pili. Our hypothesis is the following: if cells stay attached to the same position on the turbine surface for an extended time, which is considerably larger than average attachment times of individual pili, the cells will continuously pull the turbine in the same direction. This mechanism is dramatically different to the previously reported rotary microdevices driven by swimming bacteria [6,7], where cells collide with the microdevice and hence, push it. Instead, in our simulations, cells pull on the turbine. In Figure 4A, we show sketch of a cell which is attached to the turbine disk with some of its pili, while the other pili bind or unbind from the substrate. Only pili that protrude towards the turbine disks can attach, while pili protruding away from the turbine can only attach to the substrate. Thus, cells preferentially bind to the substrate in the direction away from the turbine. Since the cell is also attached with some of its pili to the turbine disks, it is pulling the turbine along the direction pointing from the turbine towards the cell. Thus, the cell pulls on the turbine, leading to a rotation of the turbine in the direction of the attached cell. This process is enhanced by the cooperation of multiple pili of a cell. To clarify this, we investigate the durations of cells attaching to the turbine disk before they detach again. To this aim, we solve a stochastic model of pili binding and unbinding, where a cell, possessing in total N 0 pili that are all growing in the direction of the turbine, has n pili attached to the turbine. With a rate k att (N 0 − n) a non-attached pilus attaches to the turbine, while with the rate k det n a pilus detaches from the turbine (see Figure 4B). By describing the attachment and detachment as Poisson processes, we can numerically solve this system with the help of a Gillespie algorithm and investigate the mean time for a complete unbinding of all pili of a cell from the turbine (see Figure 4C). We find that the mean unbinding time is increasing rapidly with the pili detachment time τ off exp(F s /F (turb) det )/k det , reaching around 10 s for τ off = 2 s and 10 7 s for τ off = 50 s. For simplicity, we assume that all pili pull with their stalling force and that the detachment force is identical to the stalling force, F s F (turb) det . The duration of how long cells bind to the turbine corresponds to the duration of superdiffusive behaviour in Figure 3A and also exhibits the same qualitative behaviour as the decay time of the angular velocity autocorrelation (see Figure 3B). We do not expect a perfect quantitative agreement between the two time scales ( Figure 3B and Figure 4C) as the simplified model of pili binding and unbinding of a cell to the turbine provided here ignores that cells can indeed move even if its pili are bound to the turbine, e.g., parallel to the turbine surface. Additionally, when multiple cells bind to the same region of the turbine, they will interact via excluded volume effects. This will lead to additional forces acting on the involved pili and thus, might enhance their detachment.
We can now provide an explanation of the three time regimes observed in the angular MSD in Figure 2C and angular velocity ACF in Figure 2D. For short times, Δt < 10 s, individual pili stochastically bind and unbind to the substrate, pulling on the turbine, and additionally, cells randomly collide with the turbine. This causes small displacements of the turbine angle α in one direction, also explaining the very sudden drop of the angular velocity ACF for short times since the life time of such displacements is very small. This leads to random fluctuations of the turbine, and as a result, to a diffusive scaling. For larger times, 10 s < Δt < 10 3 s, asymmetries in the distribution of cells on the surface of the turbine lead to a unidirectional rotation of the turbine and as a result, to a superdiffusive regime. The persistence in the turbine rotation also explains the positive angular velocity ACF. The duration of this regime depends on how strong pili bind to the turbine disk (see Figure 4C). Finally, for Δt > 10 3 s, the distribution of cells on the turbine surface is re-arranged up to such a degree that the turbine can change its direction of motion. This will then lead to a diffusive regime again. This time also corresponds to the characteristic time of the angular velocity ACF decay, confirming that for large enough times, the turbine forgets its initial direction of motion.
For cases where pili bind stronger to the turbine, we expect to not only have more pili to bind to the turbine, but also more cells to be bound to the turbine with its pili. To this aim, we check the number of cells and pili bound to the turbine as a function of time (see Supplementary Figure S2). We find that after an initial phase of pili attachment, for larger pili-turbine detachment times τ (turb) det more cells and more pili are bound to the turbine on average. This is in agreement with our previously introduced hypothesis that the turbine exhibits persistent rotational motion due to cells that are bound to the turbine surface.
Next, we consider the origin of the cell density dependent rotation of the turbine, see Figure 3C and Figure 3D. If the cell density is very small, often no cell will bind to the turbine most of the time and thus, the angular MSD is initially increasing with cell density. For moderate cell densities the difference between cells attaching to different sides of the turbine (see Figure 4D) will be significant and due to the random asymmetry, the turbine will rotate persistently in one direction until cells randomly unbind from the turbine. If the cell density gets too large, more and more cells are pulling in the opposite direction of the turbine rotation, reducing its persistence.
There are additional processes that can have a significant impact on the turbine rotation: 1) Due to the rotation of the turbine, it will constantly "pick up" cells it collides with. For the traditional rotary microdevices that are driven by swimming bacteria that collide with the device and push it [6][7][8], this would lead to a torque that acts in the opposite direction of the original rotation. When cells pull on the turbine instead, the persistent rotation is being enhanced since asymmetric distribution of cells bound to the turbine disks is getting even stronger (see Figure 4E). 2) In Section 3.1, we found that the angular velocity ACF seems to decay faster for the case with pili attachment and τ (turb) det 2 s, than for the case without pili attachment (see Figure 2D). A possible explanation for this is that even without pili-turbine attachment, cells can still rotate the turbine due to collisions and resulting excluded volume forces. That way, we reproduce a system where cells push the turbine, instead of pulling on it. This process possesses its own characteristic time scales that are linked to the substrate motility of the cells.
To summarise, the persistent rotation of the turbine originates from persistent attachment of cells to the turbine over times scales much longer than the characteristic detachment times of individual pili. While the resulting rotation is indeed persistent and, depending on the binding properties, can go on over multiple revolutions, it is not permanently unidirectional. There is no asymmetry in the initial rotational direction of the turbine and even though the turbine can move in one direction over an extended time, this rotation will at some point reverse direction. Next, we will provide an example of how cells using type IV pili can drive a permanent unidirectional rotation of a turbine.

Permanent Unidirectional Rotation due to Asymmetric Cell-Turbine Attachment
In order to produce a permanent unidirectional rotation of the turbine, an asymmetry of the turbine is required. Typically, such assymetries are created by altering the turbine geometry [7,8,12], but here we have a chance to exploit asymmetries in the adhesion of cells to the turbine instead. In Figure 5A, Figure 5B and Supplementary Movie S4, we provide an example of how altered pili binding properties on one side of the turbine wings can lead to a unidirectional rotation. Here, we consider different cases: 1) pili bind stronger to the manipulated turbine disks, corresponding to a larger detachment time τ (turb) det or 2) pili bind weaker to the manipulated turbine disks, with a lower value of τ (turb) det . Depending on how pili bind to these regions, the direction of the turbine is affected as, on average, cells are more or less strongly adhered to the manipulated region and thus, more or less cells are bound to these sides of the turbine and mediate the turbine rotation. In case (i), the turbine rotates in the direction of manipulated region because more cells are attached there, for 2) less cells are attached in the manipulated region and the turbine rotates in the opposite direction.
We also estimate the average angular velocity ω of the turbine and find values around 2-4°/s. For the used rotational mobility of the turbine μ turb (see Table 1), this corresponds to a turbine torque around 10 pN μm, comparable to previously published values for rotary devices in bacterial baths of hundreds of swimming cells [7]. Here, however, the system consists of 40 cells only. Note that this is only the lower limit of the turbine torque and the same turbine might be able to generate larger torques if it experiences an opposing torque. In the absence of a counteracting torque, the rotation speed of the turbine is limited by the retraction velocity of pili, around 2 μm/s [22].
We see evidence that the rotation is permanently unidirectional in the angular MSD (see Figure 5C) which is superdiffusive (ballistic) for arbitrary time differences Δt, equivalent to a rotation in one direction. Additionally, the normalised angular velocity ACF (see Figure 5D) is no longer decaying with time Δt, suggesting a constant rotation in the same direction.
In experiments, such manipulated regions on the surface of the turbine could be generated by coating it with a chemical such as BSA, which has been shown to alter the detachment time of the pili [17,18,32].

DISCUSSION
In this study, we investigated how bacteria that use type IV pili for surface motility may drive the rotation of a micron-sized turbine. We found that due to spatial asymmetries in the amount of cells that are bound to the turbine, a persistent rotation is observed over time scales that can strongly exceed the characteristic time scales of the pili. The persistent rotation is enhanced when pili bind more strongly to the turbine and weaker when the density of cells becomes larger. The observed persistent rotation of a symmetric turbine has a characteristic time scale. For larger times, the rotation direction is reversed stochastically and there is no preferred direction of rotation. A persistent and unidirectional rotation can be generated by altering the binding properties of the type IV pili on parts of the turbine in an asymmetric manner.
To build a rotatable micron-sized device, typically microfabrication technologies are used [5,6,8]. Established technologies make it easy to construct the turbine body, but normally the microdevices are not attached to the substrate [9]. In that case, a geometrically symmetric turbine becomes comparable to a passive rod or polymer in a bath of active particles [33,34], with two major differences: forces are exerted by pulling instead of pushing, and the attachment of a pulling cell to the microdevice can last longer compared to the single pilus retraction time. In this case, we may expect to see rotational and translational diffusion but also some more complex trajectories, e.g., where the turbine rotates itself and its center of mass is on a circular track. Interestingly, if the rod is not straight, but has a concave shape, it can lead to trapping of particles and persistent motion of the turbine [33,35,36] Harnessing the power of the turbine is considerably harder if the turbine exhibits translational motion. Hence, it is required to immobilise the center of rotation of the turbine to only allow for rotation. One way to do that would be to plug the turbine onto a immobile column. The column must have been microfabricated such that it is tightly connected to the substrate and can withstand the forces generated by the cells. An alternative solution could be to trap the turbine in a circular depression [6] that has the same radius as the turbine, so that the only degree of freedom is the turbine rotation.
Here, we have shown that microdevices immersed in a bath of bacteria that use pili, such as Neisseria gonorrhoeae or Pseudomonas aeruginosa, can efficiently harness the power of the type IV pilus machinery. This system is particularly superior to previously reported microdevices driven by swimming bacteria [6][7][8], since the involved molecular motor pilT is the strongest known molecular motor [22] and a single cell is capable of generating forces in the nano-Newton range.
In the future, it will be interesting to investigate the interplay of asymmetries in the geometrical and adhesive properties of the turbine. Additionally, we expect that the binding of pili with the substrate will also have an effect on the turbine rotation. Higher values of the pili-substrate detachment times τ (sub,1) det and τ (sub,2) det will lead to an increase in the torque exerted by attached cells on the turbine and will likely exceed the observed lower limit of the torque of around 10 pN μm. This will enable the turbine to more efficiently harness the work generated by type IV pili. Importantly, one needs to consider that if pili bind too strongly to the substrate, their substrate mobility will be weakened [19] and hence, it will be less likely that cells will be in the vicinity of the turbine. Furthermore, cells pulling in the direction opposite to the turbine rotation will exert larger forces opposing the rotation. In the future, it will be interesting to investigate the optimal substrate binding properties to maximize the turbine torque. Additionally, bacteria use type IV pili to form aggregates consisting of hundreds to thousands of cells [15,23,24,37,38]. The impact of this aggregation process on the turbine rotation remains unclear.

DATA AVAILABILITY STATEMENT
The source code of the computational model can be accessed via https://github.com/wolframponisch/Pili-Driven-Bacterial-Turbine. The Matlab scripts used to analyse the data and generate the figures are available upon request.