Dissection and integration of bursty transcriptional dynamics for complex systems.

RNA velocity estimation is a potentially powerful tool to reveal the directionality of transcriptional changes in single-cell RNA-sequencing data, but it lacks accuracy, absent advanced metabolic labeling techniques. We developed an approach, TopicVelo, that disentangles simultaneous, yet distinct, dynamics by using a probabilistic topic model, a highly interpretable form of latent space factorization, to infer cells and genes associated with individual processes, thereby capturing cellular pluripotency or multifaceted functionality. Focusing on process-associated cells and genes enables accurate estimation of process-specific velocities via a master equation for a transcriptional burst model accounting for intrinsic stochasticity. The method obtains a global transition matrix by leveraging cell topic weights to integrate process-specific signals. In challenging systems, this method accurately recovers complex transitions and terminal states, while our use of first-passage time analysis provides insights into transient transitions. These results expand the limits of RNA velocity, empowering future studies of cell fate and functional responses.

A unified model for the dynamics of ATP-independent ultrafast contraction.

In nature, several ciliated protists possess the remarkable ability to execute ultrafast motions using protein assemblies called myonemes, which contract in response to Ca2+ ions. Existing theories, such as actomyosin contractility and macroscopic biomechanical latches, do not adequately describe these systems, necessitating development of models to understand their mechanisms. In this study, we image and quantitatively analyze the contractile kinematics observed in two ciliated protists (Vorticella sp. and Spirostomum sp.), and, based on the mechanochemistry of these organisms, we propose a minimal mathematical model that reproduces our observations as well as those published previously. Analyzing the model reveals three distinct dynamic regimes, differentiated by the rate of chemical driving and the importance of inertia. We characterize their unique scaling behaviors and kinematic signatures. Besides providing insights into Ca2+-powered myoneme contraction in protists, our work may also inform the rational design of ultrafast bioengineered systems such as active synthetic cells.

Remembering the Work of Phillip L. Geissler: A Coda to His Scientific Trajectory.

Phillip L. Geissler made important contributions to the statistical mechanics of biological polymers, heterogeneous materials, and chemical dynamics in aqueous environments. He devised analytical and computational methods that revealed the underlying organization of complex systems at the frontiers of biology, chemistry, and materials science. In this retrospective we celebrate his work at these frontiers.

Active oscillatory associative memory.

Traditionally, physical models of associative memory assume conditions of equilibrium. Here, we consider a prototypical oscillator model of associative memory and study how active noise sources that drive the system out of equilibrium, as well as nonlinearities in the interactions between the oscillators, affect the associative memory properties of the system. Our simulations show that pattern retrieval under active noise is more robust to the number of learned patterns and noise intensity than under passive noise. To understand this phenomenon, we analytically derive an effective energy correction due to the temporal correlations of active noise in the limit of short correlation decay time. We find that active noise deepens the energy wells corresponding to the patterns by strengthening the oscillator couplings, where the more nonlinear interactions are preferentially enhanced. Using replica theory, we demonstrate qualitative agreement between this effective picture and the retrieval simulations. Our work suggests that the nonlinearity in the oscillator couplings can improve memory under nonequilibrium conditions.

A strong nonequilibrium bound for sorting of cross-linkers on growing biopolymers.

Understanding the role of nonequilibrium driving in self-organization is crucial for developing a predictive description of biological systems, yet it is impeded by their complexity. The actin cytoskeleton serves as a paradigm for how equilibrium and nonequilibrium forces combine to give rise to self-organization. Motivated by recent experiments that show that actin filament growth rates can tune the morphology of a growing actin bundle cross-linked by two competing types of actin-binding proteins [S. L. Freedman et al., Proc. Natl. Acad. Sci. U.S.A. 116, 16192-16197 (2019)], we construct a minimal model for such a system and show that the dynamics of a growing actin bundle are subject to a set of thermodynamic constraints that relate its nonequilibrium driving, morphology, and molecular fluxes. The thermodynamic constraints reveal the importance of correlations between these molecular fluxes and offer a route to estimating microscopic driving forces from microscopy experiments.

Giant Gateable Osmotic Power Generation from a Goldilocks Two-Dimensional Polymer.

Generating electricity from a salinity gradient, known as osmotic power, provides a sustainable energy source, but it requires precise nanoscale control of membranes for maximum performance. Here, we report an ultrathin membrane, where molecule-specific short-range interactions enable giant gateable osmotic power with a record high power density (2 kW/m2 for 1 M∥1 mM KCl). Our membranes are charge-neutral two-dimensional polymers synthesized from molecular building blocks and operate in a Goldilocks regime that simultaneously maintains high ionic conductivity and permselectivity. Molecular dynamics simulations quantitatively confirm that the functionalized nanopores are small enough for high selectivity through short-range ion-membrane interactions and large enough for fast cross-membrane transport. The short-range mechanism further enables reversible gateable operation, as demonstrated by polarity switching of osmotic power with additional gating ions.

Cluster-Spin-Glass Magnetic Behavior and Morphology in the Coordination Polymer Alloys FeyCo1-yBTT.

We report the synthesis of a series of pseudo-1D coordination polymer (CP) materials with the formula FeyCo1-yBTT (BTT = 1,3,5-benzenetrithiolate). These materials were structurally characterized by PXRD Rietveld, EXAFS, and PDF analyses, revealing that the CP superstructure enables a continuous and isomorphous alloy between the two homometallic compounds. Lower Fe loadings exhibit emergent spin glass magnetic behavior, such as memory effects and composition-dependent spin glass response time constants ranging from 6.9 × 10-9 s to 1.8 × 10-6 s. These data are consistent with the formation of spin clusters within the lattice. The magnetic behavior in these materials was modeled via replica exchange Monte Carlo simulation, which provides a good match for the experimentally measured spin glassing and magnetic phase transitions. These findings underscore how the rigid superstructure of CP and MOF scaffolds can enable the systematic tuning of physical properties, such as the spin glass behavior described here.

Simulating structured fluids with tensorial viscoelasticity.

We consider an immersed elastic body that is actively driven through a structured fluid by a motor or an external force. The behavior of such a system generally cannot be solved analytically, necessitating the use of numerical methods. However, current numerical methods omit important details of the microscopic structure and dynamics of the fluid, which can modulate the magnitudes and directions of viscoelastic restoring forces. To address this issue, we develop a simulation platform for modeling viscoelastic media with tensorial elasticity. We build on the lattice Boltzmann algorithm and incorporate viscoelastic forces, elastic immersed objects, a microscopic orientation field, and coupling between viscoelasticity and the orientation field. We demonstrate our method by characterizing how the viscoelastic restoring force on a driven immersed object depends on various key parameters as well as the tensorial character of the elastic response. We find that the restoring force depends non-monotonically on the rate of diffusion of the stress and the size of the object. We further show how the restoring force depends on the relative orientation of the microscopic structure and the pulling direction. These results imply that accounting for previously neglected physical features, such as stress diffusion and the microscopic orientation field, can improve the realism of viscoelastic simulations. We discuss possible applications and extensions to the method.

Thermodynamic Control of Activity Patterns in Cytoskeletal Networks.

Biological materials, such as the actin cytoskeleton, exhibit remarkable structural adaptability to various external stimuli by consuming different amounts of energy. In this Letter, we use methods from large deviation theory to identify a thermodynamic control principle for structural transitions in a model cytoskeletal network. Specifically, we demonstrate that biasing the dynamics with respect to the work done by nonequilibrium components effectively renormalizes the interaction strength between such components, which can eventually result in a morphological transition. Our work demonstrates how a thermodynamic quantity can be used to renormalize effective interactions, which in turn can tune structure in a predictable manner, suggesting a thermodynamic principle for the control of cytoskeletal structure and dynamics.

Mean-field theory for the structure of strongly interacting active liquids.

Active systems, which are driven out of equilibrium by local non-conservative forces, exhibit unique behaviors and structures with potential utility for the design of novel materials. An important and difficult challenge along the path toward this goal is to precisely predict how the structure of active systems is modified as their driving forces push them out of equilibrium. Here, we use tools from liquid-state theories to approach this challenge for a classic minimal active matter model. First, we construct a nonequilibrium mean-field framework that can predict the structure of systems of weakly interacting particles. Second, motivated by equilibrium solvation theories, we modify this theory to extend it with surprisingly high accuracy to systems of strongly interacting particles, distinguishing it from most existing similarly tractable approaches. Our results provide insight into spatial organization in strongly interacting out-of-equilibrium systems.

From predicting to learning dissipation from pair correlations of active liquids.

Active systems, which are driven out of equilibrium by local non-conservative forces, can adopt unique behaviors and configurations. An important challenge in the design of novel materials, which utilize such properties, is to precisely connect the static structure of active systems to the dissipation of energy induced by the local driving. Here, we use tools from liquid-state theories and machine learning to take on this challenge. We first analytically demonstrate for an isotropic active matter system that dissipation and pair correlations are closely related when driving forces behave like an active temperature. We then extend a nonequilibrium mean-field framework for predicting these pair correlations, which unlike most existing approaches is applicable even for strongly interacting particles and far from equilibrium, to predicting dissipation in these systems. Based on this theory, we reveal a robust analytic relation between dissipation and structure, which holds even as the system approaches a nonequilibrium phase transition. Finally, we construct a neural network that maps static configurations of particles to their dissipation rate without any prior knowledge of the underlying dynamics. Our results open novel perspectives on the interplay between dissipation and organization out of equilibrium.

Self-organizing motors divide active liquid droplets.

The cytoskeleton is a collection of protein assemblies that dynamically impose spatial structure in cells and coordinate processes such as cell division and mechanical regulation. Biopolymer filaments, cross-linking proteins, and enzymatically active motor proteins collectively self-organize into various precise cytoskeletal assemblies critical for specific biological functions. An outstanding question is how the precise spatial organization arises from the component macromolecules. We develop a system to investigate simple physical mechanisms of self-organization in biological assemblies. Using a minimal set of purified proteins, we create droplets of cross-linked biopolymer filaments. Through the addition of enzymatically active motor proteins, we construct composite assemblies, evocative of cellular structures such as spindles, where the inherent anisotropy drives motor self-organization, droplet deformation, and division into two droplets. These results suggest that simple physical principles underlie self-organization in complex biological assemblies and inform bioinspired materials design.

Energy dissipation and fluctuations in a driven liquid.

Minimal models of active and driven particles have recently been used to elucidate many properties of nonequilibrium systems. However, the relation between energy consumption and changes in the structure and transport properties of these nonequilibrium materials remains to be explored. We explore this relation in a minimal model of a driven liquid that settles into a time periodic steady state. Using concepts from stochastic thermodynamics and liquid state theories, we show how the work performed on the system by various nonconservative, time-dependent forces-this quantifies a violation of time reversal symmetry-modifies the structural, transport, and phase transition properties of the driven liquid.

Topological localization in out-of-equilibrium dissipative systems.

In this paper, we report that notions of topological protection can be applied to stationary configurations that are driven far from equilibrium by active, dissipative processes. We consider two physically disparate systems: stochastic networks governed by microscopic single-particle dynamics, and collections of driven interacting particles described by coarse-grained hydrodynamic theory. We derive our results by mapping to well-known electronic models and exploiting the resulting correspondence between a bulk topological number and the spectrum of dissipative modes localized at the boundary. For the Markov networks, we report a general procedure to uncover the topological properties in terms of the transition rates. For the active fluid on a substrate, we introduce a topological interpretation of fluid dissipative modes at the edge. In both cases, the presence of dissipative couplings to the environment that break time-reversal symmetry are crucial to ensuring topological protection. These examples constitute proof of principle that notions of topological protection do indeed extend to dissipative processes operating out of equilibrium. Such topologically robust boundary modes have implications for both biological and synthetic systems.

Tuning shape and internal structure of protein droplets via biopolymer filaments.

Macromolecules can phase separate to form liquid condensates, which are emerging as critical compartments in fields as diverse as intracellular organization and soft materials design. A myriad of macromolecules, including the protein FUS, form condensates which behave as isotropic liquids. Here, we investigate the influence of filament dopants on the material properties of protein liquids. We find that the short, biopolymer filaments of actin spontaneously partition into FUS droplets to form composite liquid droplets. As the concentration of the filament dopants increases, the coalescence time decreases, indicating that the dopants control viscosity relative to surface tension. The droplet shape is tunable and ranges from spherical to tactoid as the filament length or concentration is increased. We find that the tactoids are well described by a model of a quasi bipolar liquid crystal droplet, where nematic order from the anisotropic actin filaments competes with isotropic interfacial energy from the FUS, controlling droplet shape in a size-dependent manner. Our results demonstrate a versatile approach to construct tunable, anisotropic macromolecular liquids.

Robust oscillations in multi-cyclic Markov state models of biochemical clocks.

Organisms often use cyclic changes in the concentrations of chemical species to precisely time biological functions. Underlying these biochemical clocks are chemical reactions and transport processes, which are inherently stochastic. Understanding the physical basis for robust biochemical oscillations in the presence of fluctuations has thus emerged as an important problem. In a previous paper [C. del Junco and S. Vaikuntanathan, Phys. Rev. E 101, 012410 (2020)], we explored this question using the non-equilibrium statistical mechanics of single-ring Markov state models of biochemical networks that support oscillations. Our finding was that they can exploit non-equilibrium driving to robustly maintain the period and coherence of oscillations in the presence of randomness in the rates. Here, we extend our work to Markov state models consisting of a large cycle decorated with multiple small cycles. These additional cycles are intended to represent alternate pathways that the oscillator may take as it fluctuates about its average path. Combining a mapping to single-cycle networks based on first passage time distributions with our previously developed theory, we are able to make analytical predictions for the period and coherence of oscillations in these networks. One implication of our predictions is that a high energy budget can make different network topologies and arrangements of rates degenerate as far as the period and coherence of oscillations are concerned. Excellent agreement between analytical and numerical results confirms that this is the case. Our results suggest that biochemical oscillators can be more robust to fluctuations in the path of the oscillator when they have a high energy budget.

Nucleation and shape dynamics of model nematic tactoids around adhesive colloids.

Recent experiments have shown how nematically ordered tactoid shaped actin droplets can be reorganized and divided by the action of myosin molecular motors. In this paper, we consider how similar morphological changes can potentially be achieved under equilibrium conditions. Using simulations, both atomistic and continuum, and a simple macroscopic model, we explore how the nucleation dynamics, shape changes, and the final steady state of a nematic tactoid droplet can be modified by interactions with model adhesive colloids that mimic a myosin motor cluster. We show how tactoid reorganization may occur in an equilibrium colloidal-nematic setting. We then suggest based on the simple macroscopic model how the simulation models may be extended to potentially stabilize divided tactoids.

Liquid behavior of cross-linked actin bundles.

The actin cytoskeleton is a critical regulator of cytoplasmic architecture and mechanics, essential in a myriad of physiological processes. Here we demonstrate a liquid phase of actin filaments in the presence of the physiological cross-linker, filamin. Filamin condenses short actin filaments into spindle-shaped droplets, or tactoids, with shape dynamics consistent with a continuum model of anisotropic liquids. We find that cross-linker density controls the droplet shape and deformation timescales, consistent with a variable interfacial tension and viscosity. Near the liquid-solid transition, cross-linked actin bundles show behaviors reminiscent of fluid threads, including capillary instabilities and contraction. These data reveal a liquid droplet phase of actin, demixed from the surrounding solution and dominated by interfacial tension. These results suggest a mechanism to control organization, morphology, and dynamics of the actin cytoskeleton.

Topological Waves in Fluids with Odd Viscosity.

Fluids in which both time reversal and parity are broken can display a dissipationless viscosity that is odd under each of these symmetries. Here, we show how this odd viscosity has a dramatic effect on topological sound waves in fluids, including the number and spatial profile of topological edge modes. Odd viscosity provides a short-distance cutoff that allows us to define a bulk topological invariant on a compact momentum space. As the sign of odd viscosity changes, a topological phase transition occurs without closing the bulk gap. Instead, at the transition point, the topological invariant becomes ill defined because momentum space cannot be compactified. This mechanism is unique to continuum models and can describe fluids ranging from electronic to chiral active systems.

Interface height fluctuations and surface tension of driven liquids with time-dependent dynamics.

Interfaces in phase-separated driven liquids are one example of how energy input at the single-particle level changes the long-length-scale material properties of nonequilibrium systems. Here, we measure interfacial fluctuations in simulations of two liquids driven by time-dependent forces, one with repulsive interactions and one with attractive interactions. The time-dependent forces lead to currents along the interface, which can modify the scaling of interface height fluctuations with respect to predictions from capillary wave theory (CWT). We therefore characterize the whole spectrum of fluctuations to determine whether CWT applies. In the system with repulsive interactions, we find that the interface fluctuations are well-described by CWT at one amplitude of the driving forces but not at others. In the system with attractive interactions, they obey CWT for all amplitudes of driving, allowing us to extract an effective surface tension. The surface tension increases linearly over two orders of magnitude of the driving forces, more than doubling its equilibrium value. Our results show how the interfaces of nonequilibrium liquids with time-dependent forces are modified by energy input.