Multi-population dissolution in confined active fluids.

Autonomous out-of-equilibrium agents or cells in suspension are ubiquitous in biology and engineering. Turning chemical energy into mechanical stress, they generate activity in their environment, which may trigger spontaneous large-scale dynamics. Often, these systems are composed of multiple populations that may reflect the coexistence of multiple species, differing phenotypes, or chemically varying agents in engineered settings. Here, we present a new method for modeling such multi-population active fluids subject to confinement. We use a continuum multi-scale mean-field approach to represent each phase by its first three orientational moments and couple their evolution with those of the suspending fluid. The resulting coupled system is solved using a parallel adaptive level-set-based solver for high computational efficiency and maximal flexibility in the confinement geometry. Motivated by recent experimental work, we employ our method to study the spatiotemporal dynamics of confined bacterial suspensions and swarms dominated by fluid hydrodynamic effects. Our in silico explorations reproduce observed emergent collective patterns, including features of active dissolution in two-population active-passive swarms, with results clearly suggesting that hydrodynamic effects dominate dissolution dynamics. Our work lays the foundation for a systematic characterization and study of collective phenomena in natural or synthetic multi-population systems such as bacteria colonies, bird flocks, fish schools, colloid swimmers, or programmable active matter.

Natural and Engineered Isoforms of the Inflammasome Adaptor ASC Form Noncovalent, pH-Responsive Hydrogels.

The protein ASC polymerizes into intricate filament networks to assemble the inflammasome, a filamentous multiprotein complex that triggers the inflammatory response. ASC carries two Death Domains integrally involved in protein self-association for filament assembly. We have leveraged this behavior to create noncovalent, pH-responsive hydrogels of full-length, folded ASC by carefully controlling the pH as a critical factor in the polymerization process. We show that natural variants of ASC (ASC isoforms) involved in inflammasome regulation also undergo hydrogelation. To further demonstrate this general capability, we engineered proteins inspired by the ASC structure that also form hydrogels. We analyzed the structural network of the natural and engineered protein hydrogels using transmission and scanning electron microscopy and studied their viscoelastic behavior using shear rheology. Our results reveal one of the very few examples of hydrogels created by the self-assembly of globular proteins and domains in their native conformation and show that Death Domains can be used alone or as building blocks to engineer bioinspired hydrogels.

Synchronized oscillations, traveling waves, and jammed clusters induced by steric interactions in active filament arrays.

Autonomous active, elastic filaments that interact with each other to achieve cooperation and synchrony underlie many critical functions in biology. The mechanisms underlying this collective response and the essential ingredients for stable synchronization remain a mystery. Inspired by how these biological entities integrate elasticity with molecular motor activity to generate sustained oscillations, a number of synthetic active filament systems have been developed that mimic oscillations of these biological active filaments. Here, we describe the collective dynamics and stable spatiotemporal patterns that emerge in such biomimetic multi-filament arrays, under conditions where steric interactions may impact or dominate the collective dynamics. To focus on the role of steric interactions, we study the system using Brownian dynamics, without considering long-ranged hydrodynamic interactions. The simulations treat each filament as a connected chain of self-propelling colloids. We demonstrate that short-range steric inter-filament interactions and filament roughness are sufficient - even in the absence of inter-filament hydrodynamic interactions - to generate a rich variety of collective spatiotemporal oscillatory, traveling and static patterns. We first analyze the collective dynamics of two- and three-filament clusters and identify parameter ranges in which steric interactions lead to synchronized oscillations and strongly occluded states. Generalizing these results to large one-dimensional arrays, we find rich emergent behaviors, including traveling metachronal waves, and modulated wavetrains that are controlled by the interplay between the array geometry, filament activity, and filament elasticity. Interestingly, the existence of metachronal waves is non-monotonic with respect to the inter-filament spacing. We also find that the degree of filament roughness significantly affects the dynamics - specifically, filament roughness generates a locking-mechanism that transforms traveling wave patterns into statically stuck and jammed configurations. Taken together, simulations suggest that short-ranged steric inter-filament interactions could combine with complementary hydrodynamic interactions to control the development and regulation of oscillatory collective patterns. Furthermore, roughness and steric interactions may be critical to the development of jammed spatially periodic states; a spatiotemporal feature not observed in purely hydrodynamically interacting systems.

Extensions of the worm-like-chain model to tethered active filaments under tension.

Intracellular elastic filaments such as microtubules are subject to thermal Brownian noise and active noise generated by molecular motors that convert chemical energy into mechanical work. Similarly, polymers in living fluids such as bacterial suspensions and swarms suffer bending deformations as they interact with single bacteria or with cell clusters. Often, these filaments perform mechanical functions and interact with their networked environment through cross-links or have other similar constraints placed on them. Here, we examine the mechanical properties-under tension-of such constrained active filaments under canonical boundary conditions motivated by experiments. Fluctuations in the filament shape are a consequence of two types of random forces-thermal Brownian forces and activity derived forces with specified time and space correlation functions. We derive force-extension relationships and expressions for the mean square deflections for tethered filaments under various boundary conditions including hinged and clamped constraints. The expressions for hinged-hinged boundary conditions are reminiscent of the worm-like-chain model and feature effective bending moduli and mode-dependent non-thermodynamic effective temperatures controlled by the imposed force and by the activity. Our results provide methods to estimate the activity by measurements of the force-extension relation of the filaments or their mean square deflections, which can be routinely performed using optical traps, tethered particle experiments, or other single molecule techniques.

Particle diffusion in active fluids is non-monotonic in size.

We experimentally investigate the effect of particle size on the motion of passive polystyrene spheres in suspensions of Escherichia coli. Using particles covering a range of sizes from 0.6 to 39 microns, we probe particle dynamics at both short and long time scales. In all cases, the particles exhibit super-diffusive ballistic behavior at short times before eventually transitioning to diffusive behavior. Surprisingly, we find a regime in which larger particles can diffuse faster than smaller particles: the particle long-time effective diffusivity exhibits a peak in particle size, which is a deviation from classical thermal diffusion. We also find that the active contribution to particle diffusion is controlled by a dimensionless parameter, the Péclet number. A minimal model qualitatively explains the existence of the effective diffusivity peak and its dependence on bacterial concentration. Our results have broad implications on characterizing active fluids using concepts drawn from classical thermodynamics.

Facile Determination of the Poisson's Ratio and Young's Modulus of Polyacrylamide Gels and Polydimethylsiloxane.

Polyacrylamide hydrogels (PAH gel) and polydimethylsiloxane (PDMS, an elastomer) are two soft materials often used in cell mechanics and mechanobiology, in manufacturing lab-on-a-chip applications, among others. This is partly due to the ability to tune their elasticity with ease in addition to various chemical modifications. For affine polymeric networks, two (of three) elastic constants, Young's modulus (E), the shear modulus (G), and Poisson's ratio (ν), describe the purely elastic response to external forces. However, the literature addressing the experimental determination of ν for PAH (sometimes called PAA gels in the literature) and the PDMS elastomer is surprisingly limited when compared to the literature that reports values of the elastic moduli, E and G. Here, we present a facile method to obtain the Poisson's ratio and Young's modulus for PAH gel and PDMS elastomer based on static tensile tests. The value of ν obtained from the deformation of the sample is compared to the value determined by comparing E and G via a second independent method that utilizes small amplitude shear rheology. We show that the Poisson's ratio may vary significantly from the value for incompressible materials (ν = 0.5), often assumed in the literature even for soft compressible hydrogels. Surprisingly, we find a high degree of agreement between elastic constants obtained by shear rheology and macroscopic static tension test data for polyacrylamide hydrogels but not for elastomeric PDMS.

Clots reveal anomalous elastic behavior of fiber networks.

The adaptive mechanical properties of soft and fibrous biological materials are relevant to their functionality. The emergence of the macroscopic response of these materials to external stress and intrinsic cell traction from local deformations of their structural components is not well understood. Here, we investigate the nonlinear elastic behavior of blood clots by combining microscopy, rheology, and an elastic network model that incorporates the stretching, bending, and buckling of constituent fibrin fibers. By inhibiting fibrin cross-linking in blood clots, we observe an anomalous softening regime in the macroscopic shear response as well as a reduction in platelet-induced clot contractility. Our model explains these observations from two independent macroscopic measurements in a unified manner, through a single mechanical parameter, the bending stiffness of individual fibers. Supported by experimental evidence, our mechanics-based model provides a framework for predicting and comprehending the nonlinear elastic behavior of blood clots and other active biopolymer networks in general.

Size-Dependent Diffusion and Dispersion of Particles in Mucin.

Mucus, composed significantly of glycosylated mucins, is a soft and rheologically complex material that lines respiratory, reproductive, and gastrointestinal tracts in mammals. Mucus may present as a gel, as a highly viscous fluid, or as a viscoelastic fluid. Mucus acts as a barrier to the transport of harmful microbes and inhaled atmospheric pollutants to underlying cellular tissue. Studies on mucin gels have provided critical insights into the chemistry of the gels, their swelling kinetics, and the diffusion and permeability of molecular constituents such as water. The transport and dispersion of micron and sub-micron particles in mucin gels and solutions, however, differs from the motion of small molecules since the much larger tracers may interact with microstructure of the mucin network. Here, using brightfield and fluorescence microscopy, high-speed particle tracking, and passive microrheology, we study the thermally driven stochastic movement of 0.5-5.0 µm tracer particles in 10% mucin solutions at neutral pH, and in 10% mucin mixed with industrially relevant dust; specifically, unmodified limestone rock dust, modified limestone, and crystalline silica. Particle trajectories are used to calculate mean square displacements and the displacement probability distributions; these are then used to assess tracer diffusion and transport. Complex moduli are concomitantly extracted using established microrheology techniques. We find that under the conditions analyzed, the reconstituted mucin behaves as a weak viscoelastic fluid rather than as a viscoelastic gel. For small- to moderately sized tracers with a diameter of lessthan 2 µm, we find that effective diffusion coefficients follow the classical Stokes-Einstein relationship. Tracer diffusivity in dust-laden mucin is surprisingly larger than in bare mucin. Probability distributions of mean squared displacements suggest that heterogeneity, transient trapping, and electrostatic interactions impact dispersion and overall transport, especially for larger tracers. Our results motivate further exploration of physiochemical and rheological mechanisms mediating particle transport in mucin solutions and gels.

Natural and engineered isoforms of the inflammasome adaptor ASC form non-covalent, pH-responsive hydrogels.

The protein ASC polymerizes into intricate filament networks to assemble the inflammasome, a filamentous multiprotein complex that triggers the inflammatory response. ASC carries two Death Domains integrally involved in protein self-association for filament assembly. We have leveraged this behavior to create non-covalent, pH-responsive hydrogels of full-length, folded ASC by carefully controlling the pH as a critical factor in the polymerization process. We show that natural variants of ASC (ASC isoforms) involved in inflammasome regulation also undergo hydrogelation. To further demonstrate this general capability, we engineered proteins inspired in the ASC structure that successfully form hydrogels. We analyzed the structural network of the natural and engineered protein hydrogels using transmission and scanning electron microscopy, and studied their viscoelastic behavior by shear rheology. Our results reveal one of the very few examples of hydrogels created by the self-assembly of globular proteins and domains in their native conformation and show that Death Domains can be used alone or as building blocks to engineer bioinspired hydrogels.

How things get stuck: kinetics, elastohydrodynamics, and soft adhesion.

We consider the sticking of a fluid-immersed colloidal particle with a substrate coated by polymeric tethers, a model for soft, wet adhesion in many natural and artificial systems. Our theory accounts for the kinetics of binding, the elasticity of the tethers, and the hydrodynamics of fluid drainage between the colloid and the substrate, characterized by three dimensionless parameters: the ratio of the viscous drainage time to the kinetics of binding, the ratio of elastic to thermal energies, and the size of the particle relative to the height of the polymer brush. For typical experimental parameters and discrete families of tethers, we find that adhesion proceeds via punctuated steps, where rapid transitions to increasingly bound states are separated by slow aging transients, consistent with recent observations. Our results also suggest that the bound particle is susceptible to fluctuation-driven instabilities parallel to the substrate.

Glycosaminoglycans and glycoproteins influence the elastic response of synovial fluid nanofilms on model oxide surfaces.

Synovial fluid (SF) is the natural lubricant found in articulated joints, providing unique cartilage surface protecting films under confinement and relative motion. While it is known that the synergistic interactions of the macromolecular constituents provide its unique load-bearing and tribological performance, it is not fully understood how two of the main constituents, glycosaminoglycans (GAGs) and glycoproteins, regulate the formation and mechanics of robust load-bearing films. Here, we present evidence that the load-bearing capabilities, rather than the tribological performance, of the formed SF films depend strongly on its components' integrity. For this purpose, we used a combination of enzymatic treatments, quartz crystal microbalance with dissipation (QCM-D), and the surface forces apparatus (SFA) to characterize the formation and load-bearing capabilities of SF films on model oxide (i.e., silicates) surfaces. We find that, upon cleavage of proteins, the elasticity of the films is reduced and that cleaving GAGs results in irreversible (plastic) molecular re-arrangements of the film constituents when subjected to confinement. Understanding thin film mechanics of SF can provide insight into the progression of diseases, such as arthritis, but may also be applicable to the development of new implant surface treatments or new biomimetic lubricants.

Spreading rates of bacterial colonies depend on substrate stiffness and permeability.

The ability of bacteria to colonize and grow on different surfaces is an essential process for biofilm development. Here, we report the use of synthetic hydrogels with tunable stiffness and porosity to assess physical effects of the substrate on biofilm development. Using time-lapse microscopy to track the growth of expanding Serratia marcescens colonies, we find that biofilm colony growth can increase with increasing substrate stiffness, unlike what is found on traditional agar substrates. Using traction force microscopy-based techniques, we find that biofilms exert transient stresses correlated over length scales much larger than a single bacterium, and that the magnitude of these forces also increases with increasing substrate stiffness. Our results are consistent with a model of biofilm development in which the interplay between osmotic pressure arising from the biofilm and the poroelastic response of the underlying substrate controls biofilm growth and morphology.

Matrix Stiffness Modulates Mechanical Interactions and Promotes Contact between Motile Cells.

The mechanical micro-environment of cells and tissues influences key aspects of cell structure and function, including cell motility. For proper tissue development, cells need to migrate, interact, and form contacts. Cells are known to exert contractile forces on underlying soft substrates and sense deformations in them. Here, we propose and analyze a minimal biophysical model for cell migration and long-range cell-cell interactions through mutual mechanical deformations of the substrate. We compute key metrics of cell motile behavior, such as the number of cell-cell contacts over a given time, the dispersion of cell trajectories, and the probability of permanent cell contact, and analyze how these depend on a cell motility parameter and substrate stiffness. Our results elucidate how cells may sense each other mechanically and generate coordinated movements and provide an extensible framework to further address both mechanical and short-range biophysical interactions.

Three-dimensional nonlinear dynamics of prestressed active filaments: Flapping, swirling, and flipping.

Initially straight slender elastic filaments or rods with constrained ends buckle and form stable two-dimensional shapes when prestressed by bringing the ends together. Beyond a critical value of this prestress, rods can also deform off plane and form twisted three-dimensional equilibrium shapes. Here, we analyze the three-dimensional instabilities and dynamics of such deformed filaments subject to nonconservative active follower forces and fluid drag. We find that softly constrained filaments that are clamped at one end and pinned at the other exhibit stable two-dimensional planar flapping oscillations when active forces are directed toward the clamped end. Reversing the directionality of the forces quenches the instability. For strongly constrained filaments with both ends clamped, computations reveal an instability arising from the twist-bend-activity coupling. Planar oscillations are destabilized by off-planar perturbations resulting in twisted three-dimensional swirling patterns interspersed with periodic flipping or reversal of the swirling direction. These striking swirl-flip transitions are characterized by two distinct timescales: the time period for a swirl (rotation) and the time between flipping events. We interpret these reversals as relaxation oscillation events driven by accumulation of torsional energy. Each cycle is initiated by a fast jump in torsional deformation with a subsequent slow decrease in net torsion until the next cycle. Our work reveals the rich tapestry of spatiotemporal patterns when weakly inertial strongly damped rods are deformed by nonconservative active forces. Taken together, our results suggest avenues by which prestress, elasticity, and activity may be used to design synthetic macroscale pumps or mixers.

Buckling instabilities and spatio-temporal dynamics of active elastic filaments.

Biological filaments driven by molecular motors tend to experience tangential propulsive forces also known as active follower forces. When such a filament encounters an obstacle, it deforms, which reorients its follower forces and alters its entire motion. If the filament pushes a cargo, the friction on the cargo can be enough to deform the filament, thus affecting the transport properties of the cargo. Motivated by cytoskeletal filament motility assays, we study the dynamic buckling instabilities of a two-dimensional slender elastic filament driven through a dissipative medium by tangential propulsive forces in the presence of obstacles or cargo. We observe two distinct instabilities. When the filament's head is pinned or experiences significant translational but little rotational drag from its cargo, it buckles into a steadily rotating coiled state. When it is clamped or experiences both significant translational and rotational drag from its cargo, it buckles into a periodically beating, overall translating state. Using minimal analytically tractable models, linear stability theory and fully nonlinear computations, we study the onset of each buckling instability, characterize each buckled state, and map out the phase diagram of the system. Finally, we use particle-based Brownian dynamics simulations to show our main results are robust to moderate noise and steric repulsion. Overall, our results provide a unified framework to understand the dynamics of tangentially propelled filaments and filament-cargo assemblies.

Quenching active swarms: effects of light exposure on collective motility in swarming Serratia marcescens.

Swarming colonies of the light-responsive bacteria Serratia marcescens grown on agar exhibit robust fluctuating large-scale flows that include arrayed vortices, jets and sinuous streamers. We study the immobilization and quenching of these collective flows when the moving swarm is exposed to intense wide-spectrum light with a substantial ultraviolet component. We map the emergent response of the swarm to light in terms of two parameters-light intensity and duration of exposure-and identify the conditions under which collective motility is impacted. For small exposure times and/or low intensities, we find collective motility to be negligibly affected. Increasing exposure times and/or intensity to higher values suppresses collective motility but only temporarily. Terminating exposure allows bacteria to recover and eventually reestablish collective flows similar to that seen in unexposed swarms. For long exposure times or at high intensities, exposed bacteria become paralysed and form aligned, jammed regions where macroscopic speeds reduce to zero. The effective size of the quenched region increases with time and saturates to approximately the extent of the illuminated region. Post-exposure, active bacteria dislodge immotile bacteria; initial dissolution rates are strongly dependent on duration of exposure. Based on our experimental observations, we propose a minimal Brownian dynamics model to examine the escape of exposed bacteria from the region of exposure. Our results complement studies on planktonic bacteria, inform models of patterning in gradated illumination and provide a starting point for the study of specific wavelengths on swarming bacteria.

The propagation of active-passive interfaces in bacterial swarms.

Propagating interfaces are ubiquitous in nature, underlying instabilities and pattern formation in biology and material science. Physical principles governing interface growth are well understood in passive settings; however, our understanding of interfaces in active systems is still in its infancy. Here, we study the evolution of an active-passive interface using a model active matter system, bacterial swarms. We use ultra-violet light exposure to create compact domains of passive bacteria within Serratia marcescens swarms, thereby creating interfaces separating motile and immotile cells. Post-exposure, the boundary re-shapes and erodes due to self-emergent collective flows. We demonstrate that the active-passive boundary acts as a diffuse interface with mechanical properties set by the flow. Intriguingly, interfacial velocity couples to local swarm speed and interface curvature, raising the possibility that an active analogue to classic Gibbs-Thomson-Stefan conditions may control this boundary propagation.

Flagellar swimming in viscoelastic fluids: role of fluid elastic stress revealed by simulations based on experimental data.

Many important biological functions depend on microorganisms' ability to move in viscoelastic fluids such as mucus and wet soil. The effects of fluid elasticity on motility remain poorly understood, partly because the swimmer strokes depend on the properties of the fluid medium, which obfuscates the mechanisms responsible for observed behavioural changes. In this study, we use experimental data on the gaits of Chlamydomonas reinhardtii swimming in Newtonian and viscoelastic fluids as inputs to numerical simulations that decouple the swimmer gait and fluid type in order to isolate the effect of fluid elasticity on swimming. In viscoelastic fluids, cells employing the Newtonian gait swim faster but generate larger stresses and use more power, and as a result the viscoelastic gait is more efficient. Furthermore, we show that fundamental principles of swimming based on viscous fluid theory miss important flow dynamics: fluid elasticity provides an elastic memory effect that increases both the forward and backward speeds, and (unlike purely viscous fluids) larger fluid stress accumulates around flagella moving tangent to the swimming direction, compared with the normal direction.

Flagellar dynamics of a connected chain of active, polar, Brownian particles.

We show that active, self-propelled particles that are connected together to form a single chain that is anchored at one end can produce the graceful beating motions of flagella. Changing the boundary condition from a clamp to a pivot at the anchor leads to steadily rotating tight coils. Strong noise in the system disrupts the regularity of the oscillations. We use a combination of detailed numerical simulations, mean-field scaling analysis and first passage time theory to characterize the phase diagram as a function of the filament length, passive elasticity, propulsion force and noise. Our study suggests minimal experimental tests for the onset of oscillations in an active polar chain.

Dynamical self-regulation in self-propelled particle flows.

We study a continuum model of overdamped self-propelled particles with aligning interactions in two dimensions. Combining analytical theory and computations, we map out the phase diagram for the parameter space covered by the model. We find that the system self-organizes into two robust structures in different regions of parameter space: solitary waves composed of ordered swarms moving through a low density disordered background, and stationary radially symmetric asters. The self-regulating nature of the flow yields phase separation, ubiquitous in this class of systems, and controls the formation of solitary waves. Self-propulsion and the associated active convection play a crucial role in aster formation.