In higher organisms, all cells share the same genome, but every cell expresses only a limited and specific set of genes that defines the cell type. During cell division, not only the genome, but also the cell type is inherited by the daughter cells. This intriguing phenomenon is achieved by a variety of processes that have been collectively termed epigenetics: the stable and inheritable changes in gene expression patterns. This article reviews the extremely rich and exquisitely multiscale physical mechanisms that govern the biological processes behind the initiation, spreading, and inheritance of epigenetic states. These include not only the changes in the molecular properties associated with the chemical modifications of DNA and histone proteins, such as methylation and acetylation, but also less conventional changes, typically in the physics that governs the three-dimensional organization of the genome in cell nuclei. Strikingly, to achieve stability and heritability of epigenetic states, cells take advantage of many different physical principles, such as the universal behavior of polymers and copolymers, the general features of dynamical systems, and the electrostatic and mechanical properties related to chemical modifications of DNA and histones. By putting the complex biological literature in this new light, the emerging picture is that a limited set of general physical rules play a key role in initiating, shaping, and transmitting this crucial “epigenetic landscape.” This new perspective not only allows one to rationalize the normal cellular functions, but also helps to understand the emergence of pathological states, in which the epigenetic landscape becomes dysfunctional.

Localized functional domains within chromosomes, known as topologically associating domains (TADs), have been recently highlighted. In Drosophila, TADs are biochemically defined by epigenetic marks, this suggesting that the 3D arrangement may be the “missing link” between epigenetics and gene activity. Recent observations (Boettiger et al. in Nature 529(7586):418–422, 2016) provide access to structural features of these domains with unprecedented resolution thanks to super-resolution experiments. In particular, they give access to the distribution of the radii of gyration for domains of different linear length and associated with different transcriptional activity states: active, inactive or repressed. Intriguingly, the observed scaling laws lack consistent interpretation in polymer physics. We develop a new methodology conceived to extract the best information from such super-resolution data by exploiting the whole distribution of gyration radii, and to place these experimental results on a theoretical framework. We show that the experimental data are compatible with the finite-size behavior of a self-attracting polymer. The same generic polymer model leads to quantitative differences between active, inactive and repressed domains. Active domains behave as pure polymer coils, while inactive and repressed domains both lie at the coil–globule crossover. For the first time, the “color-specificity” of both the persistence length and the mean interaction energy are estimated, leading to important differences between epigenetic states. These results point toward a crucial role of criticality to enhance the system responsivity, resulting in both energy transitions and structural rearrangements. We get strong indications that epigenetically induced changes in nucleosome–nucleosome interaction can cause chromatin to shift between different activity states.

We develop a new powerful method to reproduce in silico single-molecule manipulation experiments. We demonstrate that flexible polymers such as DNA can be simulated using rigid body dynamics thanks to an original implementation of Langevin dynamics in an open source library called Open Dynamics Engine. We moreover implement a global thermostat which accelerates the simulation sampling by two orders of magnitude. We reproduce force-extension as well as rotation-extension curves of reference experimental studies. Finally, we extend the model to simulations where the control parameter is no longer the torsional strain but instead the torque, and predict the expected behavior for this case which is particularly challenging theoretically and experimentally. Video game techniques are designed to simulate rigid body dynamics of macroscopic bodies, e.g. characters or vehicles, in a realistic manner. However they are not able to deal with temperature effects, hence they are not able to deal with molecules. In order to extend these powerful techniques to molecular modeling, we implement here Langevin Dynamics in an open source library called Open Dynamics Engine. Moreover we add a “global thermostat” to this Langevin Dynamics, which accelerates the simulation sampling by two orders of magnitude. With these radically new simulation techniques, we prove that we can accurately reproduce single-molecule manipulation experiments in silico, in particular force-extension as well as rotation-extension curves of reference experimental studies. The method developed here represents an unparalleled tool for the study of more complex single molecule manipulation experiments, notably when DNA interacts with proteins. Furthermore the simulation technique that we propose here has all the functionalities required to tackle the nuclear organization of chromosomes at every length scale, from DNA to whole nuclei.