Coarse graining of biochemical systems described by discrete stochastic dynamics

Estimator of entropy production for partially accessible Markov networks based on the observation of blurred transitions

A central task in stochastic thermodynamics is the estimation of entropy production for partially accessible Markov networks. We establish an effective transition-based description for such networks with transitions that are not distinguishable and therefore blurred for an external observer. We demonstrate that, in contrast to a description based on fully resolved transitions, this effective description is typically non-Markovian at any point in time. Starting from an information-theoretic bound, we derive an operationally accessible entropy estimator for this observation scenario. We illustrate the operational relevance and the quality of this entropy estimator with a numerical analysis of various representative examples.

Waiting Time Distributions in Hybrid Models of Motor-Bead Assays: A Concept and Tool for Inference

In single-molecule experiments, the dynamics of molecular motors are often observed indirectly by measuring the trajectory of an attached bead in a motor-bead assay. In this work, we propose a method to extract the step size and stalling force for a molecular motor without relying on external control parameters. We discuss this method for a generic hybrid model that describes bead and motor via continuous and discrete degrees of freedom, respectively. Our deductions are solely based on the observation of waiting times and transition statistics of the observable bead trajectory. Thus, the method is non-invasive, operationally accessible in experiments and can, in principle, be applied to any model describing the dynamics of molecular motors. We briefly discuss the relation of our results to recent advances in stochastic thermodynamics on inference from observable transitions. Our results are confirmed by extensive numerical simulations for parameters values of an experimentally realized F1-ATPase assay.

Stochastic paths controlling speed and dissipation

Natural processes occur in a finite amount of time and dissipate energy, entropy, and matter. Near equilibrium, thermodynamic intuition suggests that fast irreversible processes will dissipate more energy and entropy than slow quasistatic processes connecting the same initial and final states. For small systems, recently discovered thermodynamic speed limits suggest that faster processes will dissipate more than slower processes. Here, we test the hypothesis that this relationship between speed and dissipation holds for stochastic paths far from equilibrium. To analyze stochastic paths on finite timescales, we derive an exact expression for the path probabilities of continuous-time Markov chains from the path summation solution to the master equation. We present a minimal model for a driven system in which relative energies of the initial and target states control the speed, and the nonequilibrium currents of a cycle control the dissipation. Although the hypothesis holds near equilibrium, we find that faster processes can dissipate less under far-from-equilibrium conditions because of strong currents. This model serves as a minimal prototype for designing kinetics to sculpt the nonequilibrium path space so that faster paths produce less dissipation.

Inferring entropy production rate from partially observed Langevin dynamics under coarse-graining

The entropy production rate (EPR) measures time-irreversibility in systems operating far from equilibrium. The challenge in estimating the EPR for a continuous variable system is the finite spatiotemporal resolution and the limited accessibility to all of the nonequilibrium degrees of freedom. Here, we estimate the irreversibility in partially observed systems following oscillatory dynamics governed by coupled overdamped Langevin equations. We coarse-grain an observed variable of a nonequilibrium driven system into a few discrete states and estimate a lower bound on the total EPR. As a model system, we use hair-cell bundle oscillations driven by molecular motors, such that the bundle tip position is observed, but the positions of the motors are hidden. In the observed variable space, the underlying driven process exhibits second-order semi-Markov statistics. The waiting time distributions (WTD), associated with transitions among the coarse-grained states, are non-exponential and convey the information on the broken time-reversal symmetry. By invoking the underlying time-irreversibility, we calculate a lower bound on the total EPR from the Kullback-Leibler divergence (KLD) between WTD. We show that the mean dwell-time asymmetry factor - the ratio between the mean dwell-times along the forward direction and the backward direction, can qualitatively measure the degree of broken time reversal symmetry and increases with finer spatial resolution. Finally, we apply our methodology to a continuous-time discrete Markov chain model, coarse-grained into a linear system exhibiting second-order semi-Markovian statistics, and demonstrate the estimation of a lower bound on the total EPR from irreversibility manifested only in the WTD.

Modeling of biomolecular machines in non-equilibrium steady states

Numerical computations have become a pillar of all modern quantitative sciences. Any computation involves modeling-even if often this step is not made explicit-and any model has to neglect details while still being physically accurate. Equilibrium statistical mechanics guides both the development of models and numerical methods for dynamics obeying detailed balance. For systems driven away from thermal equilibrium, such a universal theoretical framework is missing. For a restricted class of driven systems governed by Markov dynamics and local detailed balance, stochastic thermodynamics has evolved to fill this gap and to provide fundamental constraints and guiding principles. The next step is to advance stochastic thermodynamics from simple model systems to complex systems with tens of thousands or even millions of degrees of freedom. Biomolecules operating in the presence of chemical gradients and mechanical forces are a prime example for this challenge. In this Perspective, we give an introduction to isothermal stochastic thermodynamics geared toward the systematic multiscale modeling of the conformational dynamics of biomolecular and synthetic machines, and we outline some of the open challenges.

Excess and loss of entropy production for different levels of coarse graining

We investigate the effect of coarse graining on the thermodynamic properties of a system, focusing on entropy production. As a case of study, we consider a one-dimensional colloidal particle in contact with a thermal bath, moving in a sinusoidal potential and driven out of equilibrium by a small constant force. Different levels of coarse graining are evaluated: At first, we compare the results in the underdamped dynamics with those in the overdamped one (first coarse graining). For large values of the friction coefficient, the two dynamics have the same thermodynamics properties, while, for smaller friction values, the overdamped approximation produces an excess of entropy production with respect to that of the underdamped dynamics. Moreover, for further smaller values of the drag coefficient, the excess of entropy production turns into a loss. These regimes are explained by evaluating the jump statistics, observing that the inertia is able to induce multiple jumps and affect the average jump rate. The periodic shape of the potential allows us to approximate the continuous dynamics via a Markov chain after the introduction of a suitable time and space discretization (second level of coarse graining). This discretization procedure is implemented starting both from the underdamped and the overdamped evolution and is analyzed for different values of the friction coefficient.