Search

Far-from-equilibrium phenomena are critical to all natural and engineered systems, and essential to biological processes responsible for life. For over a century and a half, since Carnot, Clausius, Maxwell, Boltzmann, and Gibbs, among many others, laid the foundation for our understanding of equilibrium processes, scientists and engineers have dreamed of an analogous treatment of non-equilibrium systems. But despite tremendous efforts, a universal theory of non-equilibrium behavior akin to equilibrium statistical mechanics and thermodynamics has evaded description. Several methodologies have proved their ability to accurately describe complex non-equilibrium systems at the macroscopic scale, but their accuracy and predictive capacity is predicated on either phenomenological kinetic equations fit to microscopic data, or on running concurrent simulations at the particle level. Instead, we provide a framework for deriving stand-alone macroscopic thermodynamics models directly from microscopic physics without fitting in overdamped Langevin systems. The only necessary ingredient is a functional form for a parameterized, approximate density of states, in analogy to the assumption of a uniform density of states in the equilibrium microcanonical ensemble. We highlight this framework's effectiveness by deriving analytical approximations for evolving mechanical and thermodynamic quantities in a model of coiled-coil proteins and double stranded DNA, thus producing, to the authors' knowledge, the first derivation of the governing equations for a phase propagating system under general loading conditions without appeal to phenomenology. The generality of our treatment allows for application to any system described by Langevin dynamics with arbitrary interaction energies and external driving, including colloidal macromolecules, hydrogels, and biopolymers.

Constitutive modeling of dielectric elastomers has been of long standing interest in mechanics. Over the last two decades rigorous constitutive models have been developed that couple the electrical response of these polymers with large deformations characteristic of soft solids. A drawback of these models is that unlike classic models of rubber elasticity they do not consider the coupled electromechanical response of single polymer chains which must be treated using statistical mechanics. The objective of this paper is to compute the stretch and polarization of single polymer chains subject to a fixed force and fixed electric field using statistical mechanics. We assume that the dipoles induced by the applied electric field at each link do not interact with each other and compute the partition function using standard techniques. We then calculate the stretch and polarization by taking appropriate derivatives of the partition function and obtain analytical results in various limits. We also perform Markov chain Monte Carlo simulations using the Metropolis and umbrella sampling methods, as well as develop a new sampling method which improves convergence by exploiting a symmetry inherent in dielectric polymer chains. The analytical expressions are shown to agree with the Monte Carlo results over a range of forces and electric fields. Our results complement recent work on the statistical mechanics of electro-responsive chains which obtains analytical expressions in a different ensemble.

This paper aims to quantitatively relate the energy dissipated at a shock wave in a nonlinearly elastic bar to the energy in the oscillations in two related dissipationless, dispersive systems. In contrast to a phase boundary, there is no kinetic relation associated with a shock wave. Three one-dimensional dynamic impact problems are studied: Problem 1 concerns a nonlinearly elastic bar, Problem 2 a discrete chain of particles, and Problem 3 a continuum with a strain gradient term in the constitutive relation. In the impact problem considered, the free boundary of each initially quiescent body is subjected to a sudden velocity that is then held constant for all subsequent time. There is energy dissipation at the shock in Problem 1 but Problems 2 and 3 are conservative. Problem 1 is solved analytically, Problem 2 numerically and an approximate solution to Problem 3 is constructed using modulation theory. The rate of increase of the oscillatory energy in Problems 2 and 3 are calculated and compared with the dissipation rate at the shock in Problem 1. The results indicate that the former is a good measure of the latter.

We systematically study linear and nonlinear wave propagation in a chain composed of piecewise-linear bistable springs. Such bistable systems are ideal testbeds for supporting nonlinear wave dynamical features including transition and (supersonic) solitary waves. We show that bistable chains can support the propagation of subsonic wavepackets which in turn can be trapped by a low-energy phase to induce energy localization. The spatial distribution of these energy foci strongly affects the propagation of linear waves, typically causing scattering, but, in special cases, leading to a reflectionless mode analogous to the Ramsauer-Townsend (RT) effect. Further, we show that the propagation of nonlinear waves can spontaneously generate or remove additional foci, which act as effective "impurities". This behavior serves as a novel mechanism for reversibly programming the dynamic response of bistable chains.

Sedimentation in active fluids has come into focus due to the ubiquity of swimming micro-organisms in natural and industrial processes. Here, we investigate sedimentation dynamics of passive particles in a fluid as a function of bacteria E. coli concentration. Results show that the presence of swimming bacteria significantly reduces the speed of the sedimentation front even in the dilute regime, in which the sedimentation speed is expected to be independent of particle concentration. Furthermore, bacteria increase the dispersion of the passive particles, which determines the width of the sedimentation front. For short times, particle sedimentation speed has a linear dependence on bacterial concentration. Mean square displacement data shows, however, that bacterial activity decays over long experimental (sedimentation) times. An advection-diffusion equation coupled to bacteria population dynamics seems to capture concentration profiles relatively well. A single parameter, the ratio of single particle speed to the bacteria flow speed can be used to predict front sedimentation speed., Comment: Soft Matter, (2021). arXiv admin note: substantial text overlap with arXiv:1710.04068

Polymerization of dendritic actin networks underlies important mechanical processes in cell biology such as the protrusion of lamellipodia, propulsion of growth cones in dendrites of neurons, intracellular transport of organelles and pathogens, among others. The forces required for these mechanical functions have been deduced from mechano-chemical models of actin polymerization; most models are focused on single growing filaments, and only a few address polymerization of filament networks through simulations. Here we propose a continuum model of surface growth and filament nucleation to describe polymerization of dendritic actin networks. The model describes growth and elasticity in terms of macroscopic stresses, strains and filament density rather than focusing on individual filaments. The microscopic processes underlying polymerization are subsumed into kinetic laws characterizing the change of filament density and the propagation of growing surfaces. This continuum model can predict the evolution of actin networks in disparate experiments. A key conclusion of the analysis is that existing laws relating force to polymerization speed of single filaments cannot predict the response of growing networks. Therefore a new kinetic law, consistent with the dissipation inequality, is proposed to capture the evolution of dendritic actin networks under different loading conditions. This model may be extended to other settings involving a more complex interplay between mechanical stresses and polymerization kinetics, such as the growth of networks of microtubules, collagen filaments, intermediate filaments and carbon nanotubes.

Allosteric interactions in DNA are crucial for various biological processes. These interactions are quantified by measuring the change in free energy as a function of the distance between the binding sites for two ligands. Here we show that trends in the interaction energy of ligands binding to DNA can be explained within an elastic birod model. The birod model accounts for the deformation of each strand as well as the change in stacking energy due to perturbations in position and orientation of the bases caused by the binding of ligands. The strain fields produced by the ligands decay with distance from the binding site. The interaction energy of two ligands decays exponentially with the distance between them and oscillates with the periodicity of the double helix in quantitative agreement with experimental measurements. The trend in the computed interaction energy is similar to that in the perturbation of groove width produced by the binding of a single ligand which is consistent with molecular simulations. Our analysis provides a new framework to understand allosteric interactions in DNA and can be extended to other rod-like macromolecules whose elasticity plays a role in biological functions.

Allosteric interactions between molecules bound to DNA at distant locations have been known for a long time. The phenomenon has been studied via experiments and numerical simulations, but a comprehensive understanding grounded in a theory of DNA elasticity remains a challenge. Here we quantify allosteric interactions between two entities bound to DNA by using the theory of birods. We recognize that molecules bound to DNA cause local deformations that can be captured in a birod model which consists of two elastic strands interacting via an elastic web representing the base-pairs. We show that the displacement field caused by bound entities decays exponentially with distance from the binding site. We compute the interaction energy between two proteins on DNA as a function of distance between them and find that it decays exponentially while oscillating with the periodicity of the double-helix, in excellent agreement with experiments. The decay length of the interaction energy can be determined in terms of the mechanical properties of the strands and the webbing in our birod model, and it varies with the GC content of the DNA. Our model provides a framework for viewing allosteric interactions in DNA within the ambit of configurational forces of continuum elasticity.

Sedimentation in active fluids has come into focus due to the ubiquity of swimming micro-organisms in natural and artificial environments. Here, we experimentally investigate sedimentation of passive particles in water containing various concentrations of the bacterium E. coli. Results show that the presence of living bacteria reduces the velocity of the sedimentation front even in the dilute regime, where the sedimentation velocity is expected to be independent of particle concentration. Bacteria increase the effective diffusion coefficient of the passive particles, which determines the width of the sedimentation front. For higher bacteria concentration, we find the development of two sedimentation fronts due to bacterial death. A model in which an advection-diffusion equation describing the settling of particles under gravity is coupled to the population dynamics of the bacteria seems to capture the experimental trends relatively well.

We carried out dynamic force manipulations $in$ $silico$ on a variety of superhelical protein fragments from myosin, chemotaxis receptor, vimentin, fibrin, and phenylalanine zippers that vary in size and topology of their $\alpha$-helical packing. When stretched along the superhelical axis, all superhelices show elastic, plastic, and inelastic elongation regimes, and undergo a dynamic transition from the $\alpha$-helices to the $\beta$-sheets, which marks the onset of plastic deformation. Using Abeyaratne-Knowles formulation of phase transitions, we developed a theory to model mechanical and kinetic properties of protein superhelices under mechanical non-equilibrium conditions and to map their energy landscapes. The theory was validated by comparing the simulated and theoretical force-strain spectra. Scaling laws for the elastic force and the force for $\alpha$-to-$\beta$ transition to plastic deformation can be used to rationally design new materials of required mechanical strength with desired balance between stiffness and plasticity.

The rheological properties of fibrin networks have been of long-standing interest. As such there is a wealth of studies of their shear and tensile responses, but their compressive behavior remains unexplored. Here, by characterization of the network structure with synchronous measurement of the fibrin storage and loss moduli at increasing degrees of compression, we show that the compressive behavior of fibrin networks is similar to that of cellular solids. A non-linear stress-strain response of fibrin consists of three regimes: 1) an initial linear regime, in which most fibers are straight, 2) a plateau regime, in which more and more fibers buckle and collapse, and 3) a markedly non-linear regime, in which network densification occurs {{by bending of buckled fibers}} and inter-fiber contacts. Importantly, the spatially non-uniform network deformation included formation of a moving "compression front" along the axis of strain, which segregated the fibrin network into compartments with different fiber densities and structure. The Young's modulus of the linear phase depends quadratically on the fibrin volume fraction while that in the densified phase depends cubically on it. The viscoelastic plateau regime corresponds to a mixture of these two phases in which the fractions of the two phases change during compression. We model this regime using a continuum theory of phase transitions and analytically predict the storage and loss moduli which are in good agreement with the experimental data. Our work shows that fibrin networks are a member of a broad class of natural cellular materials which includes cancellous bone, wood and cork.

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 Peclet 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., Comment: 5 Figures

The effects of fluid viscosity on the kinematics of a small swimmer at low Reynolds number are investigated in both experiments and in a simple model. The swimmer is the nematode Caenorhabditis elegans, which is an undulating roundworm approximately 1 mm long. Experiments show that the nematode maintains a highly periodic swimming behavior as the fluid viscosity is varied from 1.0 mPa-s to 12 mPa-s. Surprisingly, the nematode's swimming speed (~0.35 mm/s) is nearly insensitive to the range of fluid viscosities investigated here. However, the nematode's beating frequency decreases to an asymptotic value (~1.7 Hz) with increasing fluid viscosity. A simple model is used to estimate the nematode's Young's modulus and tissue viscosity. Both material properties increase with increasing fluid viscosity. It is proposed that the increase in Young's modulus may be associated with muscle contraction in response to larger mechanical loading while the increase in effective tissue viscosity may be associated with the energy necessary to overcome increased fluid drag forces., Comment: submitted to Experimental Mechanics

Undulatory locomotion, as seen in the nematode \emph{Caenorhabditis elegans}, is a common swimming gait of organisms in the low Reynolds number regime, where viscous forces are dominant. While the nematode's motility is expected to be a strong function of its material properties, measurements remain scarce. Here, the swimming behavior of \emph{C.} \emph{elegans} are investigated in experiments and in a simple model. Experiments reveal that nematodes swim in a periodic fashion and generate traveling waves which decay from head to tail. The model is able to capture the experiments' main features and is used to estimate the nematode's Young's modulus $E$ and tissue viscosity $\eta$. For wild-type \emph{C. elegans}, we find $E\approx 3.77$ kPa and $\eta \approx-860$ Pa$\cdot$s; values of $\eta$ for live \emph{C. elegans} are negative because the tissue is generating rather than dissipating energy. Results show that material properties are sensitive to changes in muscle functional properties, and are useful quantitative tools with which to more accurately describe new and existing muscle mutants., Comment: To appear in Biophysical Journal

DNA loop formation is one of several mechanisms used by organisms to regulate genes. The free energy of forming a loop is an important factor in determining whether the associated gene is switched on or off. In this paper we use an elastic rod model of DNA to determine the free energy of forming short (50--100 basepair), protein mediated DNA loops. Superhelical stress in the DNA of living cells is a critical factor determining the energetics of loop formation, and we explicitly account for it in our calculations. The repressor protein itself is regarded as a rigid coupler; its geometry enters the problem through the boundary conditions it applies on the DNA. We show that a theory with these ingredients is sufficient to explain certain features observed in modulation of in vivo gene activity as a function of the distance between operator sites for the lac repressor. We also use our theory to make quantitative predictions for the dependence of looping on superhelical stress, which may be testable both in vivo and in single-molecule experiments such as the tethered particle assay and the magnetic bead assay., Comment: Accepted for publication in Phys Rev E. Contains the forbidden appendix they refused to publish

A variety of viruses tightly pack their genetic material into protein capsids that are barely large enough to enclose the genome. In particular, in bacteriophages, forces as high as 60 pN are encountered during packaging and ejection, produced by DNA bending elasticity and self-interactions. The high forces are believed to be important for the ejection process, though the extent of their involvement is not yet clear. As a result, there is a need for quantitative models and experiments that reveal the nature of the forces relevant to DNA ejection. Here we report measurements of the ejection forces for two different mutants of bacteriophage lambda, lambda b221cI26 and lambda cI60, which differ in genome length by ~30%. As expected for a force-driven ejection mechanism, the osmotic pressure at which DNA release is completely inhibited varies with the genome length: we find inhibition pressures of 15 atm and 25 atm, respectively, values that are in agreement with our theoretical calculations.

The conjunction of insights from structural biology, solution biochemistry, genetics and single molecule biophysics has provided a renewed impetus for the construction of quantitative models of biological processes. One area that has been a beneficiary of these experimental techniques is the study of viruses. In this paper we describe how the insights obtained from such experiments can be utilized to construct physical models of processes in the viral life cycle. We focus on dsDNA bacteriophages and show that the bending elasticity of DNA and its electrostatics in solution can be combined to determine the forces experienced during packaging and ejection of the viral genome. Furthermore, we quantitatively analyze the effect of fluid viscosity and capsid expansion on the forces experienced during packaging. Finally, we present a model for DNA ejection from bacteriophages based on the hypothesis that the energy stored in the tightly packed genome within the capsid leads to its forceful ejection. The predictions of our model can be tested through experiments in vitro where DNA ejection is inhibited by the application of external osmotic pressure.

Biophysicists and structural biologists increasingly acknowledge the role played by the mechanical properties of macromolecules as a critical element in many biological processes. This change has been brought about, in part, by the advent of single molecule biophysics techniques that have made it possible to exert piconewton forces on key macromolecules and observe their deformations at nanometer length scales, as well as to observe the mechanical action of macromolecules such as molecular motors. This has opened up immense possibilities for a new generation of mechanical investigations that will respond to such measurements in an attempt to develop a coherent theory for the mechanical behavior of macromolecules under conditions where thermal and chemical effects are on an equal footing with deterministic forces. This paper presents an application of the principles of mechanics to the problem of DNA packaging, one of the key events in the life cycle of bacterial viruses with special reference to the nature of the internal forces that are built up during the DNA packaging process., Comment: 18 pages, 7 figures, To appear in the Journal of Mechanics and Physics of Solids