Inside each one of your cells, there’s enough DNA to stretch more than two metres if laid out end-to-end. Yet somehow, all this DNA fits into a tiny structure called the cell nucleus, so intricately folded and tightly packed that it makes origami seem simple by comparison. But scientists are learning that this incredible folding isn’t just a clever trick to save space; it may also help cells make critical decisions about what kind of cells they become and how they behave.
A study published in the Bulletin of Mathematical Biology by researchers Daria Stepanova (CRM), Meritxell Brunet Guasch (University of Edinburgh), Helen M. Byrne (University of Oxford), and Tomás Alarcón (ICREA–CRM) explores how DNA’s 3D folding influences gene activity. They discovered that the shape of DNA, and how molecules interact with it, may be key to determining whether genes are switched “on” or “off.”
This switching of genes on and off without changing the DNA code is known as epigenetics. Epigenetic changes occur through chemical marks that attach to the DNA molecule or to proteins called histones, around which DNA wraps. Think of these marks as sticky notes placed along the DNA: they don’t change the instructions in our genetic code, but they tell cells which parts of those instructions to read or ignore. These notes are crucial for guiding a cell’s identity, allowing it to specialise into muscle, nerve, skin, or blood cells, even though every cell begins life with the same DNA.
However, a key mystery remains: sometimes these epigenetic marks form smooth patterns across the genome, while other times they become strikingly uneven or “rugged,” marked by abrupt shifts between areas where genes are active and areas where they’re silenced.
Until now, exactly how and why these complex, rugged patterns emerge has not been well understood.
A Virtual Lab Built with Mathematics
To solve this mystery, the research team created a mathematical model, a kind of “virtual lab,” that allows them to test biological rules and see what patterns emerge through computer simulations. Their model specifically explored the interactions between key histone modifications and the enzymes responsible for adding (“writing”) or removing (“erasing”) these chemical marks.

The model considers epigenetic modifications (H3K27me3, H3K4me3, H3K27ac) on histone H3 residues. It divides chromatin into fixed-size genomic regions characterised by these modifications, linked through reinforcing and inhibitory enzyme-driven feedback.
What set their approach apart was the inclusion of two critical elements usually overlooked in other studies. First, they included the actual 3D folding structure of chromatin (DNA wrapped around histones) because how chromatin folds changes which parts of the genome interact. Second, they factored in the limited availability and competition among enzymes. These enzymes constantly “fight” over which sites on the chromatin they modify, creating a microscopic molecular tug-of-war.
“There are not so many models of epigenetic regulation, especially those which also account for chromatin conformation; therefore, previous models have been formulated in a more simplified setting.” As Stepanova explains, “enzyme competition can be important in shaping epigenetic landscapes. It is a more realistic assumption than a constant or uniform concentration of enzyme along the chromatin. Chromatin is folded, and enzyme molecules are ‘floating’ around, ready to bind to any genomic site in their vicinity.” Moreover, she notes, this enzyme competition can “break the symmetry and eventually lead to the emergence of rugged epigenetic patterns.”
By combining chromatin folding and enzyme competition, the team found that remarkably uneven patterns emerged. When distant segments of DNA come into close contact due to folding, they influence each other’s epigenetic states. At the same time, because the enzymes needed to modify these regions are limited, only some areas get marked, leaving others untouched. The result is a patchwork, or rugged landscape, of active and inactive zones that mirrors what scientists see in real cells.
Stepanova further highlights why identifying these rugged and uniform regimes is crucial: “Epigenetic marks are known to have self-reinforcing feedback mechanisms, promoting the addition of the same or similar marks in neighbouring chromatin regions. If the region is extended, these loops would make the entire region uniform. This does not coincide with reality, since we know that, especially in differentiated cells, there are chromatin regions which are silenced and chromatin domains which are active.”
She notes that this raises an important question: “What is the ‘stopping factor’ that maintains the boundaries between these activating and inhibitory regions?” Their research provides “a possible explanation for one of the factors contributing to the emergence of alternating silenced and active genomic regions.”
Interestingly, the team’s model also captured another aspect observed in real cells called bivalent chromatin. These are regions carrying both activating and silencing epigenetic regulators at once, like genes being kept on standby, ready to quickly turn “on” or “off” in response to specific signals. These bivalent states help cells transition smoothly from uniform epigenetic landscapes to more rugged patterns, or vice versa.
Bridging Scales
The researchers designed their model to function effectively across different scales, from individual gene regions to larger chromatin structures known as Topologically Associating Domains (TADs). This mesoscale approach is particularly important. “Higher-order structures, such as loops and small TADs, can be visualised more precisely, and we can investigate how these structures influence the formation of epigenetic patterning,” Stepanova explains. “If we go to large scales of the entire chromosome, these details are also lost. In my view, the mesoscopic approach is ideal for exploring epigenetic landscapes.”
Stepanova also highlights the unique advantages mathematics offers in studying epigenetics. “Using mathematical modelling, we can investigate various aspects of biological processes which cannot be manipulated, visualised, isolated, measured, or controlled in an experimental setting,” she notes. While experimental studies are powerful, they can be costly and time-consuming. Theoretical modelling, she points out, is essentially cost-free, though closer collaboration with biologists is essential.
As an example of this type of collaboration, Stepanova mentions the current work of her master’s student Roger Bosch, who is investigating how various methods of processing chromatin conformation data affect our understanding of epigenetic patterns. These techniques aim to remove experimental artefacts and clarify the true chromatin folding patterns. However, Stepanova remarks, “It is still unclear how much these methods give us a prettier picture and how much they genuinely bring us closer to the accurate structure of cellular DNA.”
Questions like these underline both the importance and potential of integrating theoretical modelling closely with experimental research to improve our understanding of biological data.
Citation:
Stepanova, D., Brunet Guasch, M., Byrne, H.M. et al. Understanding How Chromatin Folding and Enzyme Competition Affect Rugged Epigenetic Landscapes. Bull Math Biol 87, 59 (2025). https://doi.org/10.1007/s11538-025-01434-0
crm researchers
Tomás Alarcón studied Physics, obtaining his PhD in 2000. He then pursued postdoctoral research at Oxford’s Centre for Mathematical Biology, collaborating with Helen Byrne and Philip Maini on mathematical models of tumour growth, a topic he continues researching today. After further postdoctoral positions at University College London and Imperial College, focusing on tumour dormancy, receptor dynamics, and epigenetic evolution, he led the Mathematical Biology group at the Basque Centre for Applied Mathematics (BCAM). Since 2010, he has headed the Computational & Mathematical Biology group at CRM.
Daria Stepanova is a researcher specialising in mathematical modelling applied to theoretical biology. She completed her PhD at the CRM (UAB) under the supervision of Tomás Alarcón (ICREA–CRM), Helen Byrne, and Philip Maini (University of Oxford), focusing on blood vessel growth (angiogenesis). After a research stay at Oxford studying cellular interactions linked to cancer, she joined the Laboratorio Subterráneo de Canfranc (LSC), providing mathematical support for experiments such as Hyper-Kamiokande and investigating how low-background radiation affects biological systems. Currently back at CRM as a postdoc, she continues exploring the complex interactions between epigenetic regulation and the 3D chromatin structure.
Subscribe for more CRM News
|
|
CRM CommPau Varela
|
CRM July Newsletter
CAVARET 2 brings together 83 researchers on curves and abelian varieties at the Universitat de Barcelona
The second edition of CAVARET, "Curves, Abelian VArieties and RElated Topics", took place at the Faculty of Mathematics and Computer Science of the Universitat de Barcelona from 13 to 17 July 2026. Sixteen invited lectures and ten posters brought together 83...
Mathematics on board: CRM’s demand forecasting model enters service on Barcelona’s H12 corridor
The H12 has run along Gran Via for years. What TMB switched on in July is the set of upgrades designed for it inside eBRT2030, a Horizon Europe project testing electric Bus Rapid Transit in six European cities. Behind it, forecasting passenger demand at every stop...
Dídac Gil Rams defends his doctoral thesis on the splitting of separatrices in generalized standard maps
Supervised by Inmaculada Baldomá and Pau Martín, the thesis proves the existence of chaotic dynamics in a broad family of area-preserving maps and combines classical analysis with computer-assisted proof. The defence took place on 3 July 2026 at the FME (UPC).Dídac...
How to Protect Bacteriophages Without Trapping Them? Mathematics to Improve an Alternative to Antibiotics
Bacteriophages — viruses capable of infecting and destroying bacteria — are considered one of the most promising alternatives to tackle the growing problem of antibiotic resistance. However, for these treatments to work, phages must overcome several obstacles before...
The CRM hosts nine international students for the 10th Barcelona International Youth Science Challenge
Nine secondary school students from six countries spent eight days at the CRM in July, working through four minicourses on fluid dynamics, perturbation methods, iteration and chaos, and randomness from Markov chains to machine learning, plus a session on how...
The Mathematics Behind the Three CRM Lectures at ICM 2026
Three CRM-affiliated researchers, Xavier Tolsa (plenary speaker), Xavier Cabré and Joaquim Ortega-Cerdà, lectured at the ICM 2026 in Philadelphia. Cabré surveyed the regularity of stable solutions to reaction-diffusion equations, including his optimal dimension-9...
Crossing Dimensional Frontiers: a breakthrough in the shape of stable solutions to semilinear elliptic equations
More than twenty years of efforts were required to solve a central problem in the theory of partial differential equations. A body of work developed over decades —and recognized with the Frontiers of Science Award (2023)— culminates in the resolution of an open...
Anna Florio, third CRM Chair of Excellence: billiards, friction, and the value of staying
The dynamicist Anna Florio (CEREMADE, Université Paris Dauphine-PSL) spends the summer at the CRM as the centre's third María de Maeztu Chair of Excellence. She is working with Immaculada Baldomà and Tere Martínez-Seara on billiards in the analytic setting, and on 22...
Hong Wang wins the Fields Medal, a year after her Kakeya lecture at the CRM
Hong Wang (NYU Courant and IHES) has received the Fields Medal at the International Congress of Mathematicians in Philadelphia, alongside Yu Deng (University of Chicago), John Pardon (Stony Brook University) and Jacob Tsimerman (University of Toronto). In June 2025,...
A Conference on Nonlinear PDEs on the Occasion of Xavier Cabré’s Anniversary
One anniversary, eight speakers, July and PDES (People, Discussions, Emotions and Summer) in BarcelonaThe conference Nonlinear PDEs was held at the Institut d'Estudis Catalans (IEC) on 6 and 7 July 2026. While PDEs usually stand for Partial Differential Equations, the...
Jezabel Curbelo receives the 2026 ECCOMAS Jacques-Louis Lions Award
Jezabel Curbelo, full professor at the Universitat Politècnica de Catalunya and affiliated researcher at the Centre de Recerca Matemàtica, has been awarded the ECCOMAS Jacques-Louis Lions Young Investigator Award 2026, in the field of computational mathematics. The...














