🔬 physics

Rate-Independent Epigenetics: a thermodynamically consistent framework for modelling epigenetic response

This paper proposes a thermodynamically consistent mathematical framework called Rate-Independent Epigenetics, which models heritable chromatin modifications as rate-independent dissipative systems to rigorously capture phenomena like memory, hysteresis, and catastrophic switches while establishing the existence, uniqueness, and numerical convergence of solutions.

Original authors: Jacobo Ayensa-Jiménez, Ignacio ROmero

Published 2026-07-17

📖 4 min read☕ Coffee break read

CC BY 4.0

Original authors: Jacobo Ayensa-Jiménez, Ignacio ROmero

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine your body as a bustling city where every cell is a unique neighborhood. Some neighborhoods are libraries (storing information), others are factories (making proteins), and some are parks (relaxing). But here's the twist: the DNA blueprint for every neighborhood is identical. So, how does a skin cell know to stay a skin cell and not suddenly decide to become a brain cell? The answer lies in epigenetics. Think of epigenetics as a set of sticky notes and highlighters placed on the DNA blueprint. These notes don't change the words written in the book, but they tell the cell which chapters to read and which to ignore. Once a cell puts a sticky note on a chapter to say "This is a factory," it remembers that instruction even when the original signal is gone. This "cellular memory" is why a scar stays a scar and why a cell doesn't randomly flip-flop between identities.

However, these sticky notes are tricky. They have a habit of being stubborn. If you push a cell to change its mind, it might resist until you push hard enough, and then—snap—it flips to a new state. Even if you stop pushing, it often stays in that new state. It's like a heavy door with a sticky latch: you have to shove it past a certain point to get it open, and once it's open, it won't swing back shut on its own. Scientists have long wanted a mathematical way to describe this "sticky door" behavior without getting lost in a maze of complex chemical equations. They needed a rulebook that explains how cells remember, how they resist change, and how they suddenly switch, all while obeying the fundamental laws of physics (like energy conservation).

This is where the paper "Rate-Independent Epigenetics" comes in. The authors, Jacobo Ayensa-Jiménez and Ignacio Romero, propose a new mathematical framework they call Rate-Independent Epigenetics (RIE). They treat the cell's epigenetic state not as a fluid that flows smoothly, but as a system with "stickiness" and "memory" that behaves the same way regardless of how fast you push it. Whether you push the cell's environment slowly over a year or quickly over a second, the result is the same: the cell resists until a specific threshold is crossed, then it jumps to a new state and stays there.

The paper builds a "thermodynamically consistent" skeleton for this behavior. In simple terms, this means they created a model that automatically follows the laws of energy. They didn't just guess how the cell moves; they started with three basic ingredients:

The Landscape: A map of energy hills and valleys where the cell wants to sit (like a ball rolling to the bottom of a valley).
The Push: The external environment (like a hand pushing the ball).
The Friction: A resistance that stops the ball from moving unless the push is strong enough.

By combining these, the authors derived a set of rules that predict exactly when a cell will stay put, when it will move, and when it will get "stuck" in a new position. They proved that their model works mathematically and that it never violates the laws of physics. They also built a computer program to simulate these rules. When they tested it, the program perfectly recreated known behaviors, like a cell remembering a past stressor or getting locked into a new state after a strong push.

The paper explicitly rules out the idea that these switches are just slow, smooth flows that depend on speed. Instead, it argues that for these specific types of epigenetic changes, the speed of the push doesn't matter; only the strength of the push matters. The model shows that if you push hard enough to cross a threshold, the cell jumps. If you don't, it stays put. There is no "in-between" sliding.

The authors are very confident in their mathematical proofs, showing that solutions to their equations exist and are unique under certain conditions. They also ran computer simulations to show that their method works for both simple "one-way" switches and complex "two-way" switches (like a light switch that can be on or off). While they don't claim to have solved every mystery of biology, they have provided a robust, physics-based tool that can be used to understand how cells make permanent decisions. It's like giving scientists a new, reliable ruler to measure the "stickiness" of cellular memory, ensuring that any future theories about how cells change their minds are built on a foundation that respects the laws of energy.

Technical Summary: Rate-Independent Epigenetics (RIE)

Problem Statement
Epigenetic changes are heritable, long-lived modifications of chromatin states (e.g., DNA methylation, histone modification) that exhibit distinct phenomenological features: bistability (discrete phenotypic states separated by barriers), memory (persistence of marks after stimulus removal), and hysteresis (path dependence where erasure requires crossing a distinct threshold). Current quantitative models often rely on ad-hoc evolution laws or kinetic bistable switches where hysteresis is rate-dependent. The authors identify that these features—threshold activation, path dependence, and residual states—are the hallmarks of rate-independent dissipative evolution. The paper proposes a mathematical framework, Rate-Independent Epigenetics (RIE), to model these phenomena within the rigorous theory of rate-independent systems, ensuring thermodynamic consistency by construction rather than post-hoc imposition.

Methodology
The RIE framework models epigenetic evolution as a rate-independent system defined by a triple $(Q, \mathcal{E}, \Psi)$ :

State Space ( $Q$ ): A closed, convex subset of $\mathbb{R}^n$ representing epigenetic configurations.
Stored Energy ( $\mathcal{E}$ ): A functional $\mathcal{E}(q, S) = F(q) - q \cdot \ell(S)$ , where $F$ is the configurational free energy (the "landscape," potentially non-convex/bistable) and $S(t)$ is an external micro-environmental loading path.
Dissipation Potential ( $\Psi$ ): A convex, lower semicontinuous, 1-homogeneous map $\Psi: \mathbb{R}^n \to \mathbb{R}$ encoding the resistance of the epigenetic machinery to change. The 1-homogeneity ( $\Psi(\lambda v) = \lambda \Psi(v)$ ) is the defining feature distinguishing rate-independent from viscous (rate-dependent) evolution.

The evolution is governed by an energetic principle (global stability and energy balance), leading to two equivalent formulations:

Energetic Formulation: A global stability condition coupled with an exact energy balance equation.
Subdifferential Formulation: A generalized Biot inclusion of the form $0 \in \partial\Psi(\dot{q}) + \partial_q \mathcal{E}(q, S)$ , which reduces to a sweeping process or play operator in specific cases.

Key Contributions

Thermodynamic Consistency: The framework derives governing equations that inherently satisfy the first and second laws of thermodynamics. The energy balance is an exact conservation law (work equals change in free energy plus dissipation), and the 1-homogeneity of $\Psi$ enforces non-negative, minimal (economical) dissipation, satisfying the second law without additional assumptions.
Well-Posedness Theory:
- Existence: Proven for energetic solutions under standard coercivity and continuity assumptions.
- Uniqueness: Established under convexity of the energy landscape. For the biologically relevant non-convex (bistable) case, the paper establishes the existence and scalar uniqueness of Balanced-Viscosity (BV) solutions. These solutions resolve the ambiguity of energetic jumps at switches by selecting the physically correct transition path via a vanishing-viscosity limit.
Numerical Method: The authors construct a variational time integrator based on incremental minimization. In the scalar case, this reduces to a two-sided return map (analogous to computational plasticity).
- The method is proven to converge to energetic solutions.
- It satisfies a global energy-consistency estimate of order $O(h)$ , ensuring the discrete energy balance approximates the continuous one.
- For scalar problems with monotone loading, the integrator exhibits nodal exactness, reproducing the exact solution at time nodes.

Results and Validation
The framework and numerical method were validated against three benchmarks with closed-form solutions:

Linear Play Operator (Convex): A minimal model of epigenetic memory. The results reproduce both reversible (partial erasure) and irreversible (locked-in) responses depending on the ratio of the dissipation threshold to the loading amplitude. The integrator achieved machine-precision accuracy at nodes and $O(h)$ energy consistency.
Bistable Double-Well (Non-Convex): A piecewise-quadratic potential representing competing phenotypes. The model successfully captured discontinuous jumps (catastrophic switches) between wells. The numerical scheme resolved the jump threshold ( $\ell = \rho$ ) and demonstrated that the energetic solution jumps earlier than a locally stable (balanced-viscosity) solution would, highlighting the necessity of the BV selection for physical accuracy in non-convex landscapes.
Nonlinear Convex Example: A problem with a transcendental force ( $g(q) = q + \sinh q$ ) requiring iterative Newton solvers. The method successfully recovered the exact solution at nodes and maintained first-order convergence, confirming the robustness of the approach even when closed-form inversion is impossible.

Significance and Claims
The paper claims to be the first transfer of the theory of energetic and balanced-viscosity solutions to the domain of epigenetics. Its primary significance lies in:

Conceptual Shift: Recasting epigenetic memory, thresholds, and hysteresis as the rate-independent evolution of an internal variable, moving away from ad-hoc kinetic laws.
Structural Robustness: Providing a "skeleton" for a thermodynamically consistent theory where the first and second laws are built into the variational structure.
Predictive Capability: Offering analytical thresholds and residual state predictions (epigenetic scars) that can be calibrated against experimental data (e.g., hypoxia-driven remodelling).
Extensibility: The framework serves as a foundation for future extensions, including viscous dissipation (to model rate-dependent relaxation and ratcheting) and vectorial state variables (for coupled epigenetic marks).

The authors maintain a modest stance, acknowledging that rate independence is an idealization (the quasi-static limit) and that the thermodynamic statements are "effective" within the model's internal variable framework rather than literal molecular energetics. The work is presented as a rigorous mathematical foundation for quantitative epigenetic modeling.

Technical Summary: Rate-Independent Epigenetics (RIE)

More like this