The autocatalytic reaction system with a small number of molecules is studied numerically by stochastic particle simulations. A novel state due to fluctuation and discreteness in molecular numbers is found, characterized as an extinction of molecule species alternately in the autocatalytic reaction loop. Phase transition to this state with changes of the system size and flow is studied, while a single-molecule switch of the molecule distributions is reported. The relevance of the results to intracellular processes is briefly discussed.

The mammalian genome is organized into submegabase-sized chromatin domains (CDs) including topologically associating domains, which have been identified using chromosome conformation capture-based methods. Single-nucleosome imaging in living mammalian cells has revealed subdiffusively dynamic nucleosome movement. It is unclear how single nucleosomes within CDs fluctuate and how the CD structure reflects the nucleosome movement. Here, we present a polymer model wherein CDs are characterized by fractal dimensions and the nucleosome fibers fluctuate in a viscoelastic medium with memory. We analytically show that the mean-squared displacement (MSD) of nucleosome fluctuations within CDs is subdiffusive. The diffusion coefficient and the subdiffusive exponent depend on the structural information of CDs. This analytical result enabled us to extract information from the single-nucleosome imaging data for HeLa cells. Our observation that the MSD is lower at the nuclear periphery region than the interior region indicates that CDs in the heterochromatin-rich nuclear periphery region are more compact than those in the euchromatin-rich interior region with respect to the fractal dimensions as well as the size. Finally, we evaluated that the average size of CDs is in the range of 100–500 nm and that the relaxation time of nucleosome movement within CDs is a few seconds. Our results provide physical and dynamic insights into the genome architecture in living cells. The mammalian genome is partitioned into topological chromatin domains (CDs) in the living cell nuclei. Gene expression is highly regulated within CDs according to their structure, whereas chromatin itself is highly dynamic. This raises the following question: how is the CD structure in such dynamic chromatin? We developed a conceptual framework that unifies chromatin dynamics and structure. Using a polymer model with a fractal domain structure in a viscoelastic medium, we analytically show that nucleosome movement is subdiffusive and depends on CD structure. Hence, structural information can be extracted based on nucleosome movement in living cells with single-particle tracking experiments. This framework provides physical insights into the relationship between dynamic genome organization and gene expression.

Single-molecule tracking of GPCR diffusion within the plasma membrane reveals the type of ligand with which it interacts.

Analyzing nonlinear conformational relaxation dynamics in elastic networks corresponding to two classical motor proteins, we find that they respond by well defined internal mechanical motions to various initial deformations and that these motions are robust against external perturbations. We show that this behavior is not characteristic for random elastic networks. However, special network architectures with such properties can be designed by evolutionary optimization methods. Using them, an example of an artificial elastic network, operating as a cyclic machine powered by ligand binding, is constructed.

We investigate the effects of the spatial discreteness of molecules in reaction-diffusion systems. It is found that discreteness within the so-called Kuramoto length can lead to a localization of molecules, resulting in novel steady states that do not exist in the continuous case. These states are analyzed theoretically as the fixed points of accelerated localized reactions, an approach that was verified to be in good agreement with stochastic particle simulations. The relevance of this discreteness-induced state to biological intracellular processes is discussed.

Genomes are spatiotemporally organized within the cell nucleus. Genome-wide chromosome conformation capture (Hi-C) technologies have uncovered the 3D genome organization. Furthermore, live-cell imaging experiments have revealed that genomes are functional in 4D. Although computational modeling methods can convert 2D Hi-C data into population-averaged static 3D genome models, exploring 4D genome nature based on 2D Hi-C data remains lacking. Here, we describe a 4D simulation method, PHi-C (polymer dynamics deciphered from Hi-C data), that depicts 4D genome features from 2D Hi-C data by polymer modeling. PHi-C allows users to interpret 2D Hi-C data as physical interaction parameters within single chromosomes. The physical interaction parameters can then be used in the simulations and analyses to demonstrate dynamic characteristics of genomic loci and chromosomes as observed in live-cell imaging experiments. PHi-C is available at https://github.com/soyashinkai/PHi-C.