The Mechanics and Fluctuation Spectrum of Active Gels

Heat fluctuations in chemically active systems

Chemically active systems such as living cells are maintained out of thermal equilibrium due to chemical events which generate heat and lead to active fluctuations. A key question is to understand on which time and length scales active fluctuations dominate thermal fluctuations. Here, we formulate a stochastic field theory with Poisson white noise to describe the heat fluctuations which are generated by stochastic chemical events and lead to active temperature fluctuations. We find that on large length- and timescales, active fluctuations always dominate thermal fluctuations. However, at intermediate length- and timescales, multiple crossovers exist which highlight the different characteristics of active and thermal fluctuations. Our work provides a framework to characterize fluctuations in active systems and reveals that local equilibrium holds at certain length- and timescales.

Interfacial friction and substrate deformation mediate long-range signal propagation in tissues

Tissue layers can generally slide at the interface, accompanied by the dissipation due to friction. Nevertheless, it remains elusive how force could propagate in a tissue with such interfacial friction. Here, we elaborate the force dynamics in a prototypical multilayer system in which an epithelial monolayer was cultivated upon an elastic substrate in contact with a hard surface, and discover a novel mechanism of pronounced force propagation over a long distance due to interfacial dynamics between substrate layers. We derived an analytical model for the dynamics of the elastic substrate under the shear stress provided by the cell layer at the surface boundary and the friction at bottom. The model reveals that sliding between substrate layers leads to an expanding stretch regime from a shear regime of substrate deformation in time and space. The regime boundary propagating diffusively with a speed depending on the stiffness, thickness, and slipperiness of the substrate, is a robust nature of a deformed elastic sheet with interfacial friction. These results shed new light on force propagation in tissues and our model could serve as a basis for studies of such propagation in a more complex tissue environment.

Active gels, heavy tails, and the cytoskeleton

The eukaryotic cell's cytoskeleton is a prototypical example of an active material: objects embedded within it are driven by molecular motors acting on the cytoskeleton, leading to anomalous diffusive behavior. Experiments tracking the behavior of cell-attached objects have observed anomalous diffusion with a distribution of displacements that is non-Gaussian, with heavy tails. This has been attributed to "cytoquakes" or other spatially extended collective effects. We show, using simulations and analytical theory, that a simple continuum active gel model driven by fluctuating force dipoles naturally creates heavy power-law tails in cytoskeletal displacements. We predict that this power law exponent should depend on the geometry and dimensionality of where force dipoles are distributed through the cell; we find qualitatively different results for force dipoles in a 3D cytoskeleton and a quasi-two-dimensional cortex. We then discuss potential applications of this model both in cells and in synthetic active gels.

Dynamics of undulatory fluctuations of semiflexible filaments in a network

We study the dynamics of a single semiflexible filament coupled to a Hookean spring at its boundary. The spring produces a fluctuating tensile force on the filament, the value of which depends on the filament's instantaneous end-to-end length. The spring thereby introduces a nonlinearity, which mixes the undulatory normal modes of the filament and changes their dynamics. We study these dynamics using the Martin-Siggia-Rose-Janssen-De Dominicis formalism, and compute the time-dependent correlation functions of transverse undulations and of the filament's end-to-end distance. The relaxational dynamics of the modes below a characteristic wavelength sqrt[κ/τ_{R}], set by the filament's bending modulus κ and spring-renormalized tension τ_{R}, are changed by the boundary spring. This occurs near the crossover frequency between tension- and bending-dominated modes of the system. The boundary spring can be used to represent the linear elastic compliance of the rest of the filament network to which the filament is cross linked. As a result, we predict that this nonlinear effect will be observable in the dynamical correlations of constituent filaments of networks and in the networks' collective shear response. The system's dynamic shear modulus is predicted to exhibit the well-known crossover with increasing frequency from ω^{1/2} to ω^{3/4}, but the inclusion of the network's compliance in the analysis of the individual filament dynamics shifts this transition to a higher frequency.

Oriented Active Solids

We present a complete analysis of the linearized dynamics of active solids with uniaxial orientational order, taking into account a hitherto overlooked consequence of rotation invariance. Our predictions include a purely active response of two-dimensional orientationally ordered solids to shear, the possibility of stable active solids with quasi-long-range order in two dimensions and long-range order in three dimensions, generic instability of the solid for one sign of active forcing, and the instability of the uniaxially ordered phase in momentum-conserved systems for large active forcing irrespective of its sign.

Distribution of active forces in the cell cortex

In this work, we study in detail the distribution of stochastic forces generated by the molecular motors activity, in the actin cortex of pre-muscular cells. By combining active and passive rheology experiments, performed on the same micro-bead bound to the actin network through membrane adhesive receptors, we measure the auto-correlation function C_ff(τ) of the average force pulling on the bead. As for any out-of-equilibrium system, the force distribution differs from the thermodynamical equilibrium one, especially at long time scale τ⪆ 1 s where the bead motion becomes partially directed. Thus the fluctuation-dissipation theorem does not apply and one can measure the distance from equilibrium through its violation. We investigate the influence of different parameters on the force distribution, focusing particularly on the role of ligand density: a detailed study shows how the amplitude of active forces increases when the bead is more tightly attached to the cortex. We introduce and study a model, which takes into account the number of bonds between the bead and the cytoskeleton, as well as the viscoelastic properties of the medium. This model faithfully accounts for the experimental observations. Also, it is shown that the amplitude of active forces increases with temperature. Finally, our data confirm that ATP depletion in the cell, or partial inhibition of the actomyosin activity, leads to a decrease of the amplitude of the force distribution. Altogether, we propose a consistent and quantitative description for the motion of a micrometric probe interacting with the actin network, and for the amplitude of the stochastic forces generated by molecular motors in the cortex surrounding this probe.

Active Mechanics Reveal Molecular-Scale Force Kinetics in Living Oocytes

Active diffusion of intracellular components is emerging as an important process in cell biology. This process is mediated by complex assemblies of molecular motors and cytoskeletal filaments that drive force generation in the cytoplasm and facilitate enhanced motion. The kinetics of molecular motors have been precisely characterized in vitro by single molecule approaches, but their in vivo behavior remains elusive. Here, we study the active diffusion of vesicles in mouse oocytes, where this process plays a key role in nuclear positioning during development, and combine an experimental and theoretical framework to extract molecular-scale force kinetics (force, power stroke, and velocity) of the in vivo active process. Assuming a single dominant process, we find that the nonequilibrium activity induces rapid kicks of duration τ ∼ 300 μs resulting in an average force of F ∼ 0.4 pN on vesicles in in vivo oocytes, remarkably similar to the kinetics of in vitro myosin-V. Our results reveal that measuring in vivo active fluctuations allows extraction of the molecular-scale activity in agreement with single-molecule studies and demonstrates a mesoscopic framework to access force kinetics.

Self-assembling outside equilibrium: emergence of structures mediated by dissipation

Self-assembly under non-equilibrium conditions may give rise to the formation of structures not available at equilibrium.

Broken detailed balance and non-equilibrium dynamics in living systems: a review

Living systems operate far from thermodynamic equilibrium. Enzymatic activity can induce broken detailed balance at the molecular scale. This molecular scale breaking of detailed balance is crucial to achieve biological functions such as high-fidelity transcription and translation, sensing, adaptation, biochemical patterning, and force generation. While biological systems such as motor enzymes violate detailed balance at the molecular scale, it remains unclear how non-equilibrium dynamics manifests at the mesoscale in systems that are driven through the collective activity of many motors. Indeed, in several cellular systems the presence of non-equilibrium dynamics is not always evident at large scales. For example, in the cytoskeleton or in chromosomes one can observe stationary stochastic processes that appear at first glance thermally driven. This raises the question how non-equilibrium fluctuations can be discerned from thermal noise. We discuss approaches that have recently been developed to address this question, including methods based on measuring the extent to which the system violates the fluctuation-dissipation theorem. We also review applications of this approach to reconstituted cytoskeletal networks, the cytoplasm of living cells, and cell membranes. Furthermore, we discuss a more recent approach to detect actively driven dynamics, which is based on inferring broken detailed balance. This constitutes a non-invasive method that uses time-lapse microscopy data, and can be applied to a broad range of systems in cells and tissue. We discuss the ideas underlying this method and its application to several examples including flagella, primary cilia, and cytoskeletal networks. Finally, we briefly discuss recent developments in stochastic thermodynamics and non-equilibrium statistical mechanics, which offer new perspectives to understand the physics of living systems.

The Effects of the Interplay between Motor and Brownian Forces on the Rheology of Active Gels

Active gels perform key mechanical roles inside the cell, such as cell division, motion, and force sensing. The unique mechanical properties required to perform such functions arise from the interactions between molecular motors and semiflexible polymeric filaments. Molecular motors can convert the energy released in the hydrolysis of ATP into forces of up to piconewton magnitudes. Moreover, the polymeric filaments that form active gels are flexible enough to respond to Brownian forces but also stiff enough to support the large tensions induced by the motor-generated forces. Brownian forces are expected to have a significant effect especially at motor activities at which stable noncontractile in vitro active gels are prepared for rheological measurements. Here, a microscopic mean-field theory of active gels originally formulated in the limit of motor-dominated dynamics is extended to include Brownian forces. In the model presented here, Brownian forces are included accurately, at real room temperature, even in systems with high motor activity. It is shown that a subtle interplay, or competition, between motor-generated forces and Brownian forces has an important impact on the mass transport and rheological properties of active gels. The model predictions show that at low frequencies the dynamic modulus of active gels is determined mostly by motor protein dynamics. However, Brownian forces significantly increase the breadth of the relaxation spectrum and can affect the shape of the dynamic modulus over a wide frequency range even for ratios of motor to Brownian forces of more than a hundred. Since the ratio between motor and Brownian forces is sensitive to ATP concentration, the results presented here shed some light on how the transient mechanical response of active gels changes with varying ATP concentration.

Integrated Analysis of Intracellular Dynamics of MenaINV Cancer Cells in a 3D Matrix

The intracellular environment is composed of a filamentous network that exhibits dynamic turnover of cytoskeletal components and internal force generation from molecular motors. Particle tracking microrheology enables a means to probe the internal mechanics and dynamics. Here, we develop an analytical model to capture the basic features of the active intracellular mechanical environment, including both thermal and motor-driven effects, and show consistency with a diverse range of experimental microrheology data. We further perform microrheology experiments, integrated with Brownian dynamics simulations of the active cytoskeleton, on metastatic breast cancer cells embedded in a three-dimensional collagen matrix with and without the presence of epidermal growth factor to probe the intracellular mechanical response in a physiologically mimicking scenario. Our results demonstrate that EGF stimulation can alter intracellular stiffness and power output from molecular motor-driven fluctuations in cells overexpressing an invasive isoform of the actin-associated protein Mena.

Spontaneous Formation of a Globally Connected Contractile Network in a Microtubule-Motor System

Microtubule (MT) networks play key roles in cell division, intracellular transport, and cell motility. These functions of MT networks occur through interactions between MTs and various associated proteins, notably motor proteins that bundle and slide MTs. Our objective in this study was to address the question of how motors determine the nature of MT networks. We conducted in vitro assays using homotetrameric kinesin Eg5, a motor protein involved in the formation and maintenance of the mitotic spindle. The mixing of Eg5 and MTs produced a range of spatiotemporal dynamics depending on the motor/filament ratio. Low motor/filament ratios produced globally connected static MT networks with sparsely distributed contractile active nodes (motor-accumulating points with radially extending MTs). Increasing the motor/filament ratio facilitated the linking of contractile active nodes and led to a global contraction of the network. When the motor/filament ratio was further increased, densely distributed active nodes formed local clusters and segmented the network into pieces with their strong contractile forces. Altering the properties of the motor through the use of chimeric Eg5, which has kinesin-1 heads, resulted in the generation of many isolated asters. These results suggest that the spatial distribution of contractile active nodes determines the dynamics of MT-motor networks. We then developed a coarse-grained model of MT-motor networks and identified two essential features for reproducing the experimentally observed patterns: an accumulation of motors that form the active nodes necessary to generate contractile forces, and a nonlinear dependency of contractile force on motor densities. Our model also enabled us to characterize the mechanical properties of the contractile network. Our study provides insight into how local motor-MT interactions generate the spatiotemporal dynamics of macroscopic network structures.

The effect of active fluctuations on the dynamics of particles, motors and DNA-hairpins

Inspired by recent experiments on the dynamics of particles and polymers in artificial cytoskeletons and in cells, we introduce a modified Langevin equation for a particle in an environment that is a viscoelastic medium and that is brought out of equilibrium by the action of active fluctuations caused by molecular motors. We show that within such a model, the motion of a free particle crosses over from superdiffusive to subdiffusive as observed for tracer particles in an in vitro cytoskeleton or in a cell. We investigate the dynamics of a particle confined by a harmonic potential as a simple model for the motion of the tethered head of kinesin-1. We find that the probability that the head is close to its binding site on the microtubule can be enhanced by a factor of two due to active forces. Finally, we study the dynamics of a particle in a double well potential as a model for the dynamics of DNA-hairpins. We show that the active forces effectively lower the potential barrier between the two minima and study the impact of this phenomenon on the zipping/unzipping rate.

Anomalous diffusion in viscoelastic media with active force dipoles

With the use of the "two-fluid model," we discuss anomalous diffusion induced by active force dipoles in viscoelastic media. Active force dipoles, such as proteins and bacteria, generate nonthermal fluctuating flows that lead to a substantial increment of the diffusion. Using the partial Green's function of the two-fluid model, we first obtain passive (thermal) two-point correlation functions such as the displacement cross-correlation function between the two-point particles separated by a finite distance. We then calculate active (nonthermal) one-point and two-point correlation functions due to active force dipoles. The time correlation of a force dipole is assumed to decay exponentially with a characteristic time scale. We show that the active component of the displacement cross-correlation function exhibits various crossovers from super-diffusive to subdiffusive behaviors depending on the characteristic time scales and the particle separation. Our theoretical results are intimately related to the microrheology technique to detect fluctuations in nonequilibrium environment.

Dynamics of active Rouse chains

We consider how active forces modeled as non-thermal random noise affect the average dynamical properties of a Rouse polymer. As the power spectrum of the noise is not known we keep the analytical treatment as generic as possible and then present results for a few examples of active noise. We discuss the connection between our results and recent experimental studies of dynamics of labeled DNA telomeres in living cells, and propose new chromatin tracking experiments that will allow one to determine the statistical properties of the active forces associated with chromatin remodeling processes.

Broken Detailed Balance of Filament Dynamics in Active Networks

Myosin motor proteins drive vigorous steady-state fluctuations in the actin cytoskeleton of cells. Endogenous embedded semiflexible filaments such as microtubules, or added filaments such as single-walled carbon nanotubes are used as novel tools to noninvasively track equilibrium and nonequilibrium fluctuations in such biopolymer networks. Here, we analytically calculate shape fluctuations of semiflexible probe filaments in a viscoelastic environment, driven out of equilibrium by motor activity. Transverse bending fluctuations of the probe filaments can be decomposed into dynamic normal modes. We find that these modes no longer evolve independently under nonequilibrium driving. This effective mode coupling results in nonzero circulatory currents in a conformational phase space, reflecting a violation of detailed balance. We present predictions for the characteristic frequencies associated with these currents and investigate how the temporal signatures of motor activity determine mode correlations, which we find to be consistent with recent experiments on microtubules embedded in cytoskeletal networks.

MEDYAN: Mechanochemical Simulations of Contraction and Polarity Alignment in Actomyosin Networks

Active matter systems, and in particular the cell cytoskeleton, exhibit complex mechanochemical dynamics that are still not well understood. While prior computational models of cytoskeletal dynamics have lead to many conceptual insights, an important niche still needs to be filled with a high-resolution structural modeling framework, which includes a minimally-complete set of cytoskeletal chemistries, stochastically treats reaction and diffusion processes in three spatial dimensions, accurately and efficiently describes mechanical deformations of the filamentous network under stresses generated by molecular motors, and deeply couples mechanics and chemistry at high spatial resolution. To address this need, we propose a novel reactive coarse-grained force field, as well as a publicly available software package, named the Mechanochemical Dynamics of Active Networks (MEDYAN), for simulating active network evolution and dynamics (available at www.medyan.org). This model can be used to study the non-linear, far from equilibrium processes in active matter systems, in particular, comprised of interacting semi-flexible polymers embedded in a solution with complex reaction-diffusion processes. In this work, we applied MEDYAN to investigate a contractile actomyosin network consisting of actin filaments, alpha-actinin cross-linking proteins, and non-muscle myosin IIA mini-filaments. We found that these systems undergo a switch-like transition in simulations from a random network to ordered, bundled structures when cross-linker concentration is increased above a threshold value, inducing contraction driven by myosin II mini-filaments. Our simulations also show how myosin II mini-filaments, in tandem with cross-linkers, can produce a range of actin filament polarity distributions and alignment, which is crucially dependent on the rate of actin filament turnover and the actin filament’s resulting super-diffusive behavior in the actomyosin-cross-linker system. We discuss the biological implications of these findings for the arc formation in lamellipodium-to-lamellum architectural remodeling. Lastly, our simulations produce force-dependent accumulation of myosin II, which is thought to be responsible for their mechanosensation ability, also spontaneously generating myosin II concentration gradients in the solution phase of the simulation volume.

Active matter systems have the distinct ability to convert energy from their surroundings into mechanical work, which gives rise to them having highly dynamic properties. Modeling active matter systems and capturing their complex behavior has been a great challenge in past years due to the many coupled interactions between their constituent parts, including not only distinct chemical and mechanical properties, but also feedback between them. One of the most intriguing biological active matter systems is the cell cytoskeleton, which can dynamically respond to chemical and mechanical cues to control cell structure and shape, playing a central role in many higher-order cellular processes. To model these systems and reproduce their behavior, we present a new modeling approach which combines the chemical, mechanical, and molecular transport aspects of active matter systems, all represented with equivalent complexity, while also allowing for various forms of mechanochemical feedback. This modeling approach, named MEDYAN, and software implementation is flexible so that a wide range of active matter systems can be simulated with a high level of detail, and ultimately can help to describe active matter phenomena, and in particular, the dynamics of the cell cytoskeleton. In this work, we have used MEDYAN to simulate a cytoskeletal network consisting of actin filaments, cross-linking proteins, and myosin II molecular motors. We found that these systems show rich dynamical behaviors, undergoing alignment and bundling transitions, with an emergent contractility, as the concentrations of myosin II and cross-linking proteins, as well as actin filament turnover rates, are varied.

Force fluctuations in three-dimensional suspended fibroblasts

Cells are sensitive to mechanical cues from their environment and at the same time generate and transmit forces to their surroundings. To test quantitatively forces generated by cells not attached to a substrate, we used a dual optical trap to suspend 3T3 fibroblasts between two fibronectin-coated beads. In this simple geometry, we measured both the cells' elastic properties and the force fluctuations they generate with high bandwidth. Cell stiffness decreased substantially with both myosin inhibition by blebbistatin and serum-starvation, but not with microtubule depolymerization by nocodazole. We show that cortical forces generated by non-muscle myosin II deform the cell from its rounded shape in the frequency regime from 0.1 to 10 Hz. The amplitudes of these forces were strongly reduced by blebbistatin and serum starvation, but were unaffected by depolymerization of microtubules. Force fluctuations show a spectrum that is characteristic for an elastic network activated by random sustained stresses with abrupt transitions.

Dynamics of a polymer in an active and viscoelastic bath

We study the dynamics of an ideal polymer chain in a viscoelastic medium and in the presence of active forces. The motion of the center of mass and of individual monomers is calculated. On time scales that are comparable to the persistence time of the active forces, monomers can move superdiffusively, while on larger time scales subdiffusive behavior occurs. The difference between this subdiffusion and that in the absence of active forces is quantified. We show that the polymer swells in response to active processes and determine how this swelling depends on the viscoelastic properties of the environment. Our results are compared to recent experiments on the motion of chromosomal loci in bacteria.

Molecular control of stress transmission in the microtubule cytoskeleton

In this article, we will summarize recent progress in understanding the mechanical origins of rigidity, strength, resiliency and stress transmission in the MT cytoskeleton using reconstituted networks formed from purified components. We focus on the role of network architecture, crosslinker compliance and dynamics, and molecular determinants of single filament elasticity, while highlighting open questions and future directions for this work.

Bond breaking dynamics in semiflexible networks under load

We examine the bond-breaking dynamics of transiently cross-linked semiflexible networks using a single filament model in which that filament is peeled from an array of cross-linkers. We examine the effect of quenched disorder in the placement of the linkers along the filament and the effect of stochastic bond-breaking (assuming Bell model unbinding kinetics) on the dynamics of filament cross-linker dissociation and the statistics of ripping events. We find that bond forces decay exponentially away from the point of loading and that bond breaking proceeds sequentially down the linker array from the point of loading in a series of stochastic ripping events. We compare these theoretical predictions to the observed trajectories of large beads in a cross-linked microtubule network and identify the observed jumps of the bead with the linker rupture events predicted by the single filament model.

Nonthermal fluctuations of the mitotic spindle

We present direct measurements of fluctuations in the nucleus of yeast cells. While prior work has shown these fluctuations to be active and non-thermal in character, their origin and time dependence are not understood. We show that the nuclear fluctuations we observe are quantitatively consistent with uncorrelated, active force fluctuations driving a nuclear medium that is dominated by an uncondensed DNA solution, for which we perform rheological measurements on an in vitro model system under similar conditions to what is expected in the nucleus.

Active cell mechanics: Measurement and theory

Living cells are active mechanical systems that are able to generate forces. Their structure and shape are primarily determined by biopolymer filaments and molecular motors that form the cytoskeleton. Active force generation requires constant consumption of energy to maintain the nonequilibrium activity to drive organization and transport processes necessary for their function. To understand this activity it is necessary to develop new approaches to probe the underlying physical processes. Active cell mechanics incorporates active molecular-scale force generation into the traditional framework of mechanics of materials. This review highlights recent experimental and theoretical developments towards understanding active cell mechanics. We focus primarily on intracellular mechanical measurements and theoretical advances utilizing the Langevin framework. These developing approaches allow a quantitative understanding of nonequilibrium mechanical activity in living cells. This article is part of a Special Issue entitled: Mechanobiology.

The mechanical properties of early Drosophila embryos measured by high-speed video microrheology

In early development, Drosophila melanogaster embryos form a syncytium, i.e., multiplying nuclei are not yet separated by cell membranes, but are interconnected by cytoskeletal polymer networks consisting of actin and microtubules. Between division cycles 9 and 13, nuclei and cytoskeleton form a two-dimensional cortical layer. To probe the mechanical properties and dynamics of this self-organizing pre-tissue, we measured shear moduli in the embryo by high-speed video microrheology. We recorded position fluctuations of injected micron-sized fluorescent beads with kHz sampling frequencies and characterized the viscoelasticity of the embryo in different locations. Thermal fluctuations dominated over nonequilibrium activity for frequencies between 0.3 and 1000 Hz. Between the nuclear layer and the yolk, the cytoplasm was homogeneous and viscously dominated, with a viscosity three orders of magnitude higher than that of water. Within the nuclear layer we found an increase of the elastic and viscous moduli consistent with an increased microtubule density. Drug-interference experiments showed that microtubules contribute to the measured viscoelasticity inside the embryo whereas actin only plays a minor role in the regions outside of the actin caps that are closely associated with the nuclei. Measurements at different stages of the nuclear division cycle showed little variation.

Substrate stiffness-modulated registry phase correlations in cardiomyocytes map structural order to coherent beating

Recent experiments show that both striation, an indication of the structural registry in muscle fibres, as well as the contractile strains produced by beating cardiac muscle cells can be optimized by substrate stiffness. Here we show theoretically how the substrate rigidity dependence of the registry data can be mapped onto that of the strain measurements. We express the elasticity-mediated structural registry as a phase-order parameter using a statistical physics approach that takes the noise and disorder inherent in biological systems into account. By assuming that structurally registered myofibrils also tend to beat in phase, we explain the observed dependence of both striation and strain measurements of cardiomyocytes on substrate stiffness in a unified manner. The agreement of our ideas with experiment suggests that the correlated beating of heart cells may be limited by the structural order of the myofibrils, which in turn is regulated by their elastic environment.

The role of filament length, finite-extensibility and motor force dispersity in stress relaxation and buckling mechanisms in non-sarcomeric active gels

After relaxing some assumptions we apply a single-chain mean-field mathematical model recently introduced [RSC Adv. (2014)] to describe the role of molecular motors in the mechanical properties of active gels. The model allows physics that are not available in models postulated on coarser levels of description. Moreover it proposes a level of description that allows the prediction of observables at time scales too difficult to achieve in multi-chain simulations for realistic filament lengths and densities. We model the semiflexible filaments that compose the active gel as bead-spring chains; molecular motors are accounted for by using a mean-field approach, in which filaments undergo transitions of one motor attachment state depending on the state of the probe filament. The level of description includes the end-to-end distance and attachment state of the filaments, and the motor-generated forces, as stochastic state variables which evolve according to a proposed differential Chapman-Kolmogorov equation. The motor-generated forces are drawn from a stationary distribution of motor stall forces. We consider bead-spring chains with multiple beads, explore the effect of finite-extensibility of the strands and incorporate into the model motor force distributions that have been measured experimentally. The model can no longer be solved analytically but is amenable to numerical simulation. This version of the model allows a more quantitative description of buckling dynamics [Lenz et. al. PRL, 2012, 108, 238107] and the dynamic modulus of active gels. The effect of finite extensibility of the filament strands on the dynamic modulus was also found to be in agreement with the microrheology experiments of Mizuno et. al., [Science, 2007, 315, 370-373].

Impact of dimensionality and network disruption on microrheology of cancer cells in 3D environments

Dimensionality is a fundamental component that can have profound implications on the characteristics of physical systems. In cell biology, however, the majority of studies on cell physical properties, from rheology to force generation to migration, have been performed on 2D substrates, and it is not clear how a more realistic 3D environment influences cell properties. Here, we develop an integrated approach and demonstrate the combination of mitochondria-tracking microrheology, microfluidics, and Brownian dynamics simulations to explore the impact of dimensionality on intracellular mechanics and on the effects of intracellular disruption. Additionally, we consider both passive thermal and active motor-driven processes within the cell and demonstrate through modeling how active internal fluctuations are modulated via dimensionality. Our results demonstrate that metastatic breast cancer cells (MDA-MB-231) exhibit more solid-like internal motions in 3D compared to 2D, and actin network disruption via Cytochalasin D has a more pronounced effect on internal cell fluctuations in 2D. Our computational results and modeling show that motor-induced active stress fluctuations are enhanced in 2D, leading to increased local intracellular particle fluctuations and apparent fluid-like behavior.

Biomechanical properties at the cellular and subcellular levels are important in providing proper biological functions, from cell migratory capabilities to intracellular transport. Deregulation in these properties can lead to disease states such as cancer metastasis. We develop and demonstrate an integrated experimental and computational approach to study intracellular mechanics. We demonstrate that a key environmental factor, dimensionality, plays a significant role in modulating intracellular mechanical behavior. This is important as typical cell biology and mechanics experiments are performed on 2D substrates, which do not capture the physiological features of 3D matrices and may not induce physiologically accurate cell properties. We further develop an effective temperature model to describe how dimensionality changes intracellular particle motion by altering the activity of molecular motors.

High-resolution mapping of intracellular fluctuations using carbon nanotubes

Cells are active systems with molecular force generation that drives complex dynamics at the supramolecular scale. We present a quantitative study of molecular motions in cells over times from milliseconds to hours. Noninvasive tracking was accomplished by imaging highly stable near-infrared luminescence of single-walled carbon nanotubes targeted to kinesin-1 motor proteins in COS-7 cells. We observed a regime of active random "stirring" that constitutes an intermediate mode of transport, different from both thermal diffusion and directed motor activity. High-frequency motion was found to be thermally driven. At times greater than 100 milliseconds, nonequilibrium dynamics dominated. In addition to directed transport along microtubules, we observed strong random dynamics driven by myosins that result in enhanced nonspecific transport. We present a quantitative model connecting molecular mechanisms to mesoscopic fluctuations.

Dynamics of elastic interactions in soft and biological matter

Cells probe their mechanical environment and can change the organization of their cytoskeletons when the elastic and viscous properties of their environment are modified. We use a model in which the forces exerted by small, contractile acto-myosin filaments (e.g., nascent stress fibers in stem cells) on the extracellular matrix are modeled as local force dipoles. In some cases, the strain field caused by these force dipoles propagates quickly enough so that only static elastic interactions need be considered. On the other hand, in the case of significant energy dissipation, strain propagation is slower and may be eliminated completely by the relaxation of the cellular cytoskeleton (e.g., by cross-link dissociation). Here, we consider several dissipative mechanisms that affect the propagation of the strain field in adhered cells and consider these effects on the interaction between force dipoles and their resulting mutual orientations. This is a first step in understanding the development of orientational (nematic) or layering (smectic) order in the cytoskeleton. We use the theory to estimate the propagation time of the strain fields over a cellular distance for different mechanisms and find that in some cases it can be of the order of seconds, thus competing with the cytoskeletal relaxation time. Furthermore, for a simple system of two force dipoles, we predict that in some cases the orientation of force dipoles might change significantly with time, e.g., for short times the dipoles exhibit parallel alignment while for later times they align perpendicularly.

Microrheology of highly crosslinked microtubule networks is dominated by force-induced crosslinker unbinding

We determine the time- and force-dependent viscoelastic responses of reconstituted networks of microtubules that have been strongly crosslinked by biotin-streptavidin bonds. To measure the microscale viscoelasticity of such networks, we use a magnetic tweezers device to apply localized forces. At short time scales, the networks respond nonlinearly to applied force, with stiffening at small forces, followed by a reduction in the stiffening response at high forces, which we attribute to the force-induced unbinding of crosslinks. At long time scales, force-induced bond unbinding leads to local network rearrangement and significant bead creep. Interestingly, the network retains its elastic modulus even under conditions of significant plastic flow, suggesting that crosslinker breakage is balanced by the formation of new bonds. To better understand this effect, we developed a finite element model of such a stiff filament network with labile crosslinkers obeying force-dependent Bell model unbinding dynamics. The coexistence of dissipation, due to bond breakage, and the elastic recovery of the network is possible because each filament has many crosslinkers. Recovery can occur as long as a sufficient number of the original crosslinkers are preserved under the loading period. When these remaining original crosslinkers are broken, plastic flow results.

Stochastic thermodynamics, fluctuation theorems and molecular machines

Stochastic thermodynamics as reviewed here systematically provides a framework for extending the notions of classical thermodynamics such as work, heat and entropy production to the level of individual trajectories of well-defined non-equilibrium ensembles. It applies whenever a non-equilibrium process is still coupled to one (or several) heat bath(s) of constant temperature. Paradigmatic systems are single colloidal particles in time-dependent laser traps, polymers in external flow, enzymes and molecular motors in single molecule assays, small biochemical networks and thermoelectric devices involving single electron transport. For such systems, a first-law like energy balance can be identified along fluctuating trajectories. For a basic Markovian dynamics implemented either on the continuum level with Langevin equations or on a discrete set of states as a master equation, thermodynamic consistency imposes a local-detailed balance constraint on noise and rates, respectively. Various integral and detailed fluctuation theorems, which are derived here in a unifying approach from one master theorem, constrain the probability distributions for work, heat and entropy production depending on the nature of the system and the choice of non-equilibrium conditions. For non-equilibrium steady states, particularly strong results hold like a generalized fluctuation-dissipation theorem involving entropy production. Ramifications and applications of these concepts include optimal driving between specified states in finite time, the role of measurement-based feedback processes and the relation between dissipation and irreversibility. Efficiency and, in particular, efficiency at maximum power can be discussed systematically beyond the linear response regime for two classes of molecular machines, isothermal ones such as molecular motors, and heat engines such as thermoelectric devices, using a common framework based on a cycle decomposition of entropy production.

Active, motor-driven mechanics in a DNA gel

Cells are capable of a variety of dramatic stimuli-responsive mechanical behaviors. These capabilities are enabled by the pervading cytoskeletal network, an active gel composed of structural filaments (e.g., actin) that are acted upon by motor proteins (e.g., myosin). Here, we describe the synthesis and characterization of an active gel using noncytoskeletal components. We use methods of base-pair-templated DNA self assembly to create a hybrid DNA gel containing stiff tubes and flexible linkers. We then activate the gel by adding the motor FtsK50C, a construct derived from the bacterial protein FtsK that, in vitro, has a strong and processive DNA contraction activity. The motors stiffen the gel and create stochastic contractile events that affect the positions of attached beads. We quantify the fluctuations of the beads and show that they are comparable both to measurements of cytoskeletal systems and to theoretical predictions for active gels. Thus, we present a DNA-based active gel whose behavior highlights the universal aspects of nonequilibrium, motor-driven networks.

Tensegrity and motor-driven effective interactions in a model cytoskeleton

Actomyosin networks are major structural components of the cell. They provide mechanical integrity and allow dynamic remodeling of eukaryotic cells, self-organizing into the diverse patterns essential for development. We provide a theoretical framework to investigate the intricate interplay between local force generation, network connectivity, and collective action of molecular motors. This framework is capable of accommodating both regular and heterogeneous pattern formation, arrested coarsening and macroscopic contraction in a unified manner. We model the actomyosin system as a motorized cat's cradle consisting of a crosslinked network of nonlinear elastic filaments subjected to spatially anti-correlated motor kicks acting on motorized (fibril) crosslinks. The phase diagram suggests there can be arrested phase separation which provides a natural explanation for the aggregation and coalescence of actomyosin condensates. Simulation studies confirm the theoretical picture that a nonequilibrium many-body system driven by correlated motor kicks can behave as if it were at an effective equilibrium, but with modified interactions that account for the correlation of the motor driven motions of the actively bonded nodes. Regular aster patterns are observed both in Brownian dynamics simulations at effective equilibrium and in the complete stochastic simulations. The results show that large-scale contraction requires correlated kicking.

Active contractility in actomyosin networks

Contractile forces are essential for many developmental processes involving cell shape change and tissue deformation. Recent experiments on reconstituted actomyosin networks, the major component of the contractile machinery, have shown that active contractility occurs above a threshold motor concentration and within a window of cross-link concentration. We present a microscopic dynamic model that incorporates two essential aspects of actomyosin self-organization: the asymmetric load response of individual actin filaments and the correlated motor-driven events mimicking myosin-induced filament sliding. Using computer simulations, we examine how the concentration and susceptibility of motors contribute to their collective behavior and interplay with the network connectivity to regulate macroscopic contractility. Our model is shown to capture the formation and dynamics of contractile structures and agree with the observed dependence of active contractility on microscopic parameters, including the contractility onset. Cooperative action of load-resisting motors in a force-percolating structure integrates local contraction/buckling events into a global contractile state via an active coarsening process, in contrast to the flow transition driven by uncorrelated kicks of susceptible motors.

Effective temperature of red-blood-cell membrane fluctuations

Biologically driven nonequilibrium fluctuations are often characterized by their non-Gaussianity or by an "effective temperature", which is frequency dependent and higher than the ambient temperature. We address these two measures theoretically by examining a randomly kicked particle, with a variable number of kicking motors, and show how these two indicators of nonequilibrium behavior can contradict. Our results are compared with new experiments on shape fluctuations of red-blood cell membranes, and demonstrate how the physical nature of the motors in this system can be revealed using these global measures of nonequilibrium.

Physics met biology, and the consequence was

We summarise the contributions to the discussion and the links between them. The complex relationship between the physical and biological sciences demonstrates three "axes of tension": the role of simulation, the interplay between levels of explanation, and the generality of "laws". We identify examples of true synergy between approaches that genuinely explore new research territory, and underscore the contemporary value of the type of discussions contained in this volume.

Cytoskeleton fluidization versus resolidification: prestress effect

The differential elastic modulus of an active actomyosin network is computed as a function of applied stress, taking into account both thermal and motor contributions to filament compliance in the low-frequency domain. It is shown that, due to a dual nature of motor activity, increasing motor concentration may either stiffen the network due to stronger prestress or soften it due to motor agitation, in accordance with experimental data. Prestress anisotropy, which may be induced by redistribution of motors triggered by external force, causes anisotropy of the elastic moduli. This helps to explain the contradictory phenomena of cell fluidization and resolidification in response to transient stretch observed in recent experiments.

Motor-driven effective temperature and viscoelastic response of active matter

We consider dynamic response of a cytoskeletal network to both thermal and motor-induced fluctuations. The latter are viewed in two independent ways, as either additive or multiplicative colored noise. Due to a natural upper frequency limit of the motor agitation, the response of a living cell is similar to that of an equilibrium system in the high-frequency domain. At lower frequencies, the role of motor agitation manifests itself in intensified network fluctuations, which is equivalent to effective growth of the environment temperature. The effective temperature becomes frequency dependent, which signifies violation of the conventional fluctuation-dissipation theorem. The motor action affects the dynamic shear modulus in two opposite ways: by stiffening the network through filament prestress and softening it through increased agitation. The latter tendency is isolated when only single-headed motors are present. The theory is in good agreement with experimental measurements of the amplitude of the shear modulus under these conditions.

In vivo determination of fluctuating forces during endosome trafficking using a combination of active and passive microrheology

Background

Regulation of intracellular trafficking is a central issue in cell biology. The forces acting on intracellular vesicles (endosomes) can be assessed in living cells by using a combination of active and passive microrheology.

Methodology/Principal Findings

This dual approach is based on endosome labeling with magnetic nanoparticles. The resulting magnetic endosomes act both as probes that can be manipulated with external magnetic fields to infer the viscoelastic modulus of their surrounding microenvironment, and as biological vehicles that are trafficked along the microtubule network by means of forces generated by molecular motors. The intracellular viscoelastic modulus exhibits power law dependence with frequency, which is microtubule and actin-dependent. The mean square displacements of endosomes do not follow the predictions of the fluctuation-dissipation theorem, which offers evidence for active force generation. Microtubule disruption brings the intracellular medium closer to thermal equilibrium: active forces acting on the endosomes depend on microtubule-associated motors. The power spectra of these active forces, deduced through the use of a generalized Langevin equation, show a power law decrease with frequency and reveal an actin-dependent persistence of the force with time. Experimental spectra have been reproduced by a simple model consisting in a series of force steps power-law distributed in time. This model enlightens the role of the cytoskeleton dependent force exerted on endosomes to perform intracellular trafficking.

Conclusions/Significance

In this work, the influence of cytoskeleton components and molecular motors on intracellular viscoelasticity and transport is addressed. The use of an original probe, the magnetic endosome, allows retrieving the power spectrum of active forces on these organelles thanks to interrelated active and passive measures.

Finally a computational model gives estimates of the force itself and hence of the number of the motors pulling on endosomes.

Nonlocal fluctuation correlations in active gels

Many active materials and biological systems are driven far from equilibrium by embedded agents that spontaneously generate forces and distort the surrounding material. Probing and characterizing these athermal fluctuations are essential to understand the properties and behaviors of such systems. Here we present a mathematical procedure to estimate the local action of force-generating agents from the observed fluctuating displacement fields. The active agents are modeled as oriented force dipoles or isotropic compression foci, and the matrix on which they act is assumed to be either a compressible elastic continuum or a coupled network-solvent system. Correlations at a single point and between points separated by an arbitrary distance are obtained, giving a total of three independent fluctuation modes that can be tested with microrheology experiments. Since oriented dipoles and isotropic compression foci give different contributions to these fluctuation modes, ratiometric analysis allows us characterize the force generators. We also predict and experimentally find a high-frequency ballistic regime, arising from individual force-generating events in the form of the slow buildup of stress followed by rapid but finite decay. Finally, we provide a quantitative statistical model to estimate the mean filament tension from these athermal fluctuations, which leads to stiffening of active networks.

Multiphase flow models of biogels from crawling cells to bacterial biofilms

This article reviews multiphase descriptions of the fluid mechanics of cytoplasm in crawling cells and growing bacterial biofilms. These two systems involve gels, which are mixtures composed of a polymer network permeated by water. The fluid mechanics of these systems is essential to their biological function and structure. Their mathematical descriptions must account for the mechanics of the polymer, the water, and the interaction between these two phases. This review focuses on multiphase flow models because this framework is natural for including the relative motion between the phases, the exchange of material between phases, and the additional stresses within the network that arise from nonspecific chemical interactions and the action of molecular motors. These models have been successful in accounting for how different forces are generated and transmitted to achieve cell motion and biofilm growth and they have demonstrated how emergent structures develop though the interactions of the two phases. A short description of multiphase flow models of tumor growth is included to highlight the flexibility of the model in describing diverse biological applications.

DNA condensation by field-induced non-equilibrium noise

DNA molecules under good solvent conditions condense when subjected to strong AC electrical fields. It is shown that AC electrical fields couple charge fluctuations to molecular shape fluctuations, which produces a non-equilibrium noise source that can dominate over equilibrium thermal noise in the long-wavelength limit. The field-induced excess charge fluctuations amplify the Asakura-Oosawa fluctuation attraction force between neighboring chains, providing a mechanism for field-induced DNA condensation.