Confinement-Induced Glassy Dynamics in a Model for Chromosome Organization

Hi-C-guided many-polymer model to decipher 3D genome organization

We propose a high-throughput chromosome conformation capture data-based many-polymer model that allows us to generate an ensemble of multi-scale genome structures. We demonstrate the efficacy of our model by validating the generated structures against experimental measurements and employ them to address key questions regarding genome organization. Our model first confirms a significant correlation between chromosome size and nuclear positioning. Specifically, smaller chromosomes are distributed at the core region, whereas larger chromosomes are at the periphery, interacting with the nuclear envelope. The spatial distribution of A- and B-type compartments, which is nontrivial to infer from the corresponding high-throughput chromosome conformation capture maps alone, can also be elucidated using our model, accounting for an issue such as the effect of chromatin-lamina interaction on the compartmentalization of conventional and inverted nuclei. In accordance with imaging data, the overall shape of the 3D genome structures generated from our model displays significant variation. As a case study, we apply our method to the yellow fever mosquito genome, finding that the predicted morphology displays, on average, a more globular shape than the previously suggested spindle-like organization and that our prediction better aligns with the fluorescence in situ hybridization data. Our model has great potential to be extended to investigate many outstanding issues concerning 3D genome organization.

Nonequilibrium switching of segmental states can influence compaction of chromatin

Knowledge about the dynamic nature of chromatin organization is essential to understand the regulation of processes like DNA transcription and repair. The existing models of chromatin assume that protein organization and chemical states along chromatin are static and the 3D organization is purely a result of protein-mediated intra-chromatin interactions. Here we present a new hypothesis that certain nonequilibrium processes, such as switching of chemical and physical states due to nucleosome assembly/disassembly or gene repression/activation, can also simultaneously influence chromatin configurations. To understand the implications of this inherent nonequilibrium switching, we present a block copolymer model of chromatin, with switching of its segmental states between two states, mimicking active/repressed or protein unbound/bound states. We show that competition between switching timescale Tt, polymer relaxation timescale τp, and segmental relaxation timescale τs can lead to non-trivial changes in chromatin organization, leading to changes in local compaction and contact probabilities. As a function of the switching timescale, the radius of gyration of chromatin shows a non-monotonic behavior with a prominent minimum when Tt ≈ τp and a maximum when Tt ≈ τs. We find that polymers with a small segment length exhibit a more compact structure than those with larger segment lengths. We also find that the switching can lead to higher contact probability and better mixing of far-away segments. Our study also shows that the nature of the distribution of chromatin clusters varies widely as we change the switching rate.

Transcription-induced active forces suppress chromatin motion

The organization of interphase chromosomes in a number of species is starting to emerge thanks to advances in a variety of experimental techniques. However, much less is known about the dynamics, especially in the functional states of chromatin. Some experiments have shown that the motility of individual loci in human interphase chromosome decreases during transcription and increases upon inhibiting transcription. This is a counterintuitive finding because it is thought that the active mechanical force (F) on the order of ten piconewtons, generated by RNA polymerase II (RNAPII) that is presumably transmitted to the gene-rich region of the chromatin, would render it more open, thus enhancing the mobility. We developed a minimal active copolymer model for interphase chromosomes to investigate how F affects the dynamical properties of chromatin. The movements of the loci in the gene-rich region are suppressed in an intermediate range of F and are enhanced at small F values, which has also been observed in experiments. In the intermediate F, the bond length between consecutive loci increases, becoming commensurate with the distance at the minimum of the attractive interaction between nonbonded loci. This results in a transient disorder-to-order transition, leading to a decreased mobility during transcription. Strikingly, the F-dependent change in the locus dynamics preserves the organization of the chromosome at [Formula: see text]. Transient ordering of the loci, which is not found in the polymers with random epigenetic profiles, in the gene-rich region might be a plausible mechanism for nucleating a dynamic network involving transcription factors, RNAPII, and chromatin.

Mesoscale, long-time mixing of chromosomes and its connection to polymer dynamics

Chromosomes are arranged in distinct territories within the nucleus of animal cells. Recent experiments have shown that these territories overlap at their edges, suggesting partial mixing during interphase. Experiments that knock-down of condensin II proteins during interphase indicate increased chromosome mixing, which demonstrates control of the mixing. In this study, we use a generic polymer simulation to quantify the dynamics of chromosome mixing over time. We introduce the chromosome mixing index, which quantifies the mixing of distinct chromosomes in the nucleus. We find that the chromosome mixing index in a small confinement volume (as a model of the nucleus), increases as a power-law of the time, with the scaling exponent varying non-monotonically with self-interaction and volume fraction. By comparing the chromosome mixing index with both monomer subdiffusion due to (non-topological) intermingling of chromosomes as well as even slower reptation, we show that for relatively large volume fractions, the scaling exponent of the chromosome mixing index is related to Rouse dynamics for relatively weak chromosome attractions and to reptation for strong attractions. In addition, we extend our model to more realistically account for the situation of the Drosophila chromosome by including the heterogeneity of the polymers and their lengths to account for microphase separation of euchromatin and heterochromatin and their interactions with the nuclear lamina. We find that the interaction with the lamina further impedes chromosome mixing.

Multivalent binding proteins can drive collapse and reswelling of chromatin in confinement

Collapsed conformations of chromatin have been long suspected of being mediated by interactions with multivalent binding proteins, which can bring together distant sections of the chromatin fiber. In this study, we use Langevin dynamics simulation of a coarse grained chromatin polymer to show that the role of binding proteins can be more nuanced than previously suspected. In particular, for chromatin polymer in confinement, entropic forces can drive reswelling of collapsed chromatin with increasing binder concentrations, and this reswelling transition happens at physiologically relevant binder concentrations. Both the extent of collapse, and also of reswelling depends on the strength of confinement. We also study the kinetics of collapse and reswelling and show that both processes occur in similar timescales. We characterise this reswelling of chromatin in biologically relevant regimes and discuss the non-trivial role of multivalent binding proteins in mediating the spatial organisation of the genome.

Symmetry-based classification of forces driving chromatin dynamics

Chromatin - the functional form of DNA in the cell - exists in the form of a polymer immersed in a nucleoplasmic fluid inside the cell nucleus. Both chromatin and nucleoplasm are subject to active forces resulting from local biological processes. This activity leads to non-equilibrium phenomena, affecting chromatin organization and dynamics, yet the underlying physics is far from understood. Here, we expand upon a previously developed two-fluid model of chromatin and nucleoplasm by considering three types of activity in the form of force dipoles - two with both forces of the dipole acting on the same fluid (either polymer or nucleoplasm) and a third, with two forces pushing chromatin and solvent in opposite directions. We find that this latter type results in the most significant flows, dominating over most length scales of interest. Due to the friction between the fluids and their viscosity, we observe emergent screening length scales in the active flows of this system. We predict that the presence of different activity types and their relative strengths can be inferred from observing the power spectra of hydrodynamic fluctuations in the chromatin and the nucleoplasm.

Topological gelation of reconnecting polymers

Shuffling of genetic material via reconnection of proximal DNA segments is seen in meiosis and some cancers. In contrast with textbook pictures, reconnections are performed within an entangled environment and can result in knotting or linking, detrimental for cells. By performing Brownian dynamics simulations of reconnecting polymers under confinement, modeling the genome in vivo, we find a topological transition between a gas or liquid of unlinked rings and a gel-like structure with a large number of polydisperse linked rings. This transition can be triggered by increasing polymer stiffness or confinement. Our results suggest ways to design future topological materials, such as DNA-based gels involving recombinase proteins.

DNA recombination is a ubiquitous process that ensures genetic diversity. Contrary to textbook pictures, DNA recombination, as well as generic DNA translocations, occurs in a confined and highly entangled environment. Inspired by this observation, here, we investigate a solution of semiflexible polymer rings undergoing generic cutting and reconnection operations under spherical confinement. Our setup may be realized using engineered DNA in the presence of recombinase proteins or by considering micelle-like components able to form living (or reversibly breakable) polymer rings. We find that in such systems, there is a topological gelation transition, which can be triggered by increasing either the stiffness or the concentration of the rings. Flexible or dilute polymers break into an ensemble of short, unlinked, and segregated rings, whereas sufficiently stiff or dense polymers self-assemble into a network of long, linked, and mixed loops, many of which are knotted. We predict that the two phases should behave qualitatively differently in elution experiments monitoring the escape dynamics from a permeabilized container. Besides shedding some light on the biophysics and topology of genomes undergoing DNA reconnection in vivo, our findings could be leveraged in vitro to design polymeric complex fluids—e.g., DNA-based complex fluids or living polymer networks—with desired topologies.

A simulation model of heterochromatin formation at submolecular detail

Heterochromatin is a physical state of the chromatin fiber that maintains gene repression during cell development. Although evidence exists on molecular mechanisms involved in heterochromatin formation, a detailed structural mechanism of heterochromatin formation needs a better understanding. We made use of a simple Monte Carlo simulation model with explicit representation of key molecular events to observe molecular self-organization leading to heterochromatin formation. Our simulations provide a structural interpretation of several important traits of the heterochromatinization process. In particular, this study provides a depiction of how small amounts of HP1 are able to induce a highly condensed chromatin state through HP1 dimerization and bridging of sequence-remote nucleosomes. It also elucidates structural roots of a yet poorly understood phenomenon of a nondeterministic nature of heterochromatin formation and subsequent gene repression. Experimental chromatin in vivo assay provides an unbiased estimate of time scale of repressive response to a heterochromatin-triggering event.

Formulation of Chromatin Mobility as a Function of Nuclear Size during C. elegans Embryogenesis Using Polymer Physics Theories

During early embryogenesis of the nematode, Caenorhabditis elegans, the chromatin motion markedly decreases. Despite its biological implications, the underlying mechanism for this transition was unclear. By combining theory and experiment, we analyze the mean-square displacement (MSD) of the chromatin loci, and demonstrate that MSD-vs-time relationships in various nuclei collapse into a single master curve by normalizing them with the mesh size and the corresponding time scale. This enables us to identify the onset of the entangled dynamics with the size of tube diameter of chromatin polymer in the C. elegans embryo. Our dynamical scaling analysis predicts the transition between unentangled and entangled dynamics of chromatin polymers, the quantitative formula for MSD as a function of nuclear size and timescale, and provides testable hypotheses on chromatin mobility in other cell types and species.

Dynamic and facilitated binding of topoisomerase accelerates topological relaxation

How type 2 Topoisomerase (TopoII) proteins relax and simplify the topology of DNA molecules is one of the most intriguing open questions in genome and DNA biophysics. Most of the existing models neglect the dynamics of TopoII which is expected of proteins searching their targets via facilitated diffusion. Here, we show that dynamic binding of TopoII speeds up the topological relaxation of knotted substrates by enhancing the search of the knotted arc. Intriguingly, this in turn implies that the timescale of topological relaxation is virtually independent of the substrate length. We then discover that considering binding biases due to facilitated diffusion on looped substrates steers the sampling of the topological space closer to the boundaries between different topoisomers yielding an optimally fast topological relaxation. We discuss our findings in the context of topological simplification in vitro and in vivo.

A Liquid State Perspective on Dynamics of Chromatin Compartments

The interior of the eukaryotic cell nucleus has a crowded and heterogeneous environment packed with chromatin polymers, regulatory proteins, and RNA molecules. Chromatin polymer, assisted by epigenetic modifications, protein and RNA binders, forms multi-scale compartments which help regulate genes in response to cellular signals. Furthermore, chromatin compartments are dynamic and tend to evolve in size and composition in ways that are not fully understood. The latest super-resolution imaging experiments have revealed a much more dynamic and stochastic nature of chromatin compartments than was appreciated before. An emerging mechanism explaining chromatin compartmentalization dynamics is the phase separation of protein and nucleic acids into membraneless liquid condensates. Consequently, concepts and ideas from soft matter and polymer systems have been rapidly entering the lexicon of cell biology. In this respect, the role of computational models is crucial for establishing a rigorous and quantitative foundation for the new concepts and disentangling the complex interplay of forces that contribute to the emergent patterns of chromatin dynamics and organization. Several multi-scale models have emerged to address various aspects of chromatin dynamics, ranging from equilibrium polymer simulations, hybrid non-equilibrium simulations coupling protein binding and chromatin folding, and mesoscopic field-theoretic models. Here, we review these emerging theoretical paradigms and computational models with a particular focus on chromatin’s phase separation and liquid-like properties as a basis for nuclear organization and dynamics.

Extracting multi-way chromatin contacts from Hi-C data

There is a growing realization that multi-way chromatin contacts formed in chromosome structures are fundamental units of gene regulation. However, due to the paucity and complexity of such contacts, it is challenging to detect and identify them using experiments. Based on an assumption that chromosome structures can be mapped onto a network of Gaussian polymer, here we derive analytic expressions for n-body contact probabilities (n > 2) among chromatin loci based on pairwise genomic contact frequencies available in Hi-C, and show that multi-way contact probability maps can in principle be extracted from Hi-C. The three-body (triplet) contact probabilities, calculated from our theory, are in good correlation with those from measurements including Tri-C, MC-4C and SPRITE. Maps of multi-way chromatin contacts calculated from our analytic expressions can not only complement experimental measurements, but also can offer better understanding of the related issues, such as cell-line dependent assemblies of multiple genes and enhancers to chromatin hubs, competition between long-range and short-range multi-way contacts, and condensates of multiple CTCF anchors.

The interplay of chromatin phase separation and lamina interactions in nuclear organization

The genetic material of eukaryotes is segregated into transcriptionally active euchromatin and silent heterochromatin compartments. The spatial arrangement of chromatin compartments evolves over the course of cellular life in a process that remains poorly understood. The latest nuclear imaging experiments reveal a number of dynamical signatures of chromatin that are reminiscent of active multiphase liquids. This includes the observations of viscoelastic response, coherent motions, Ostwald ripening, and coalescence of chromatin compartments. There is also growing evidence that liquid-liquid phase separation of protein and nucleic acid components is the underlying mechanism for the dynamical behavior of chromatin. To dissect the organizational and dynamical implications of chromatin's liquid behavior, we have devised a phenomenological field-theoretic model of the nucleus as a multiphase condensate of liquid chromatin types. Employing the liquid chromatin model of the Drosophila nucleus, we have carried out an extensive set of simulations with an objective to shed light on the dynamics and chromatin patterning observed in the latest nuclear imaging experiments. Our simulations reveal the emergence of experimentally detected mesoscale chromatin channels and spheroidal droplets which arise from the dynamic interplay of chromatin type to type interactions and intermingling of chromosomal territories. We also quantitatively reproduce coherent motions of chromatin domains observed in displacement correlation spectroscopy measurements which are explained within the framework of our model by phase separation of chromatin types operating within constrained intrachromosomal and interchromosomal boundaries. Finally, we illuminate the role of heterochromatin-lamina interactions in the nuclear organization by showing that these interactions enhance the mobility of euchromatin and indirectly introduce correlated motions of heterochromatin droplets.

Multiscale modeling of genome organization with maximum entropy optimization

Three-dimensional (3D) organization of the human genome plays an essential role in all DNA-templated processes, including gene transcription, gene regulation, and DNA replication. Computational modeling can be an effective way of building high-resolution genome structures and improving our understanding of these molecular processes. However, it faces significant challenges as the human genome consists of over 6 × 109 base pairs, a system size that exceeds the capacity of traditional modeling approaches. In this perspective, we review the progress that has been made in modeling the human genome. Coarse-grained models parameterized to reproduce experimental data via the maximum entropy optimization algorithm serve as effective means to study genome organization at various length scales. They have provided insight into the principles of whole-genome organization and enabled de novo predictions of chromosome structures from epigenetic modifications. Applications of these models at a near-atomistic resolution further revealed physicochemical interactions that drive the phase separation of disordered proteins and dictate chromatin stability in situ. We conclude with an outlook on the opportunities and challenges in studying chromosome dynamics.

Dynamic Nuclear Structure Emerges from Chromatin Cross-Links and Motors

The cell nucleus houses the chromosomes, which are linked to a soft shell of lamin protein filaments. Experiments indicate that correlated chromosome dynamics and nuclear shape fluctuations arise from motor activity. To identify the physical mechanisms, we develop a model of an active, cross-linked Rouse chain bound to a polymeric shell. System-sized correlated motions occur but require both motor activity and cross-links. Contractile motors, in particular, enhance chromosome dynamics by driving anomalous density fluctuations. Nuclear shape fluctuations depend on motor strength, cross-linking, and chromosome-lamina binding. Therefore, complex chromosome dynamics and nuclear shape emerge from a minimal, active chromosome-lamina system.

Random First Order Transition Theory for Glassy Dynamics in a Single Condensed Polymer

The number of compact structures of a single condensed polymer (SCP), with similar free energies, grows exponentially with the degree of polymerization. In analogy with structural glasses (SGs), we expect that at low temperatures chain relaxation should occur by activated transitions between the compact metastable states. By evolving the states of the SCP, linearly coupled to a reference state, we show that, below a dynamical transition temperature (T_{d}), the SCP is trapped in a metastable state leading to slow dynamics. At a lower temperature, T_{K}≠0, the configurational entropy vanishes, resulting in a thermodynamic random first order ideal glass transition. The relaxation time obeys the Vogel-Fulcher-Tamman law, diverging at T=T_{0}≈T_{K}. These findings, accord well with the random first order transition theory, establishing that SCP and SG exhibit similar universal characteristics.

A unified framework for inferring the multi-scale organization of chromatin domains from Hi-C

Chromosomes are giant chain molecules organized into an ensemble of three-dimensional structures characterized with its genomic state and the corresponding biological functions. Despite the strong cell-to-cell heterogeneity, the cell-type specific pattern demonstrated in high-throughput chromosome conformation capture (Hi-C) data hints at a valuable link between structure and function, which makes inference of chromatin domains (CDs) from the pattern of Hi-C a central problem in genome research. Here we present a unified method for analyzing Hi-C data to determine spatial organization of CDs over multiple genomic scales. By applying statistical physics-based clustering analysis to a polymer physics model of the chromosome, our method identifies the CDs that best represent the global pattern of correlation manifested in Hi-C. The multi-scale intra-chromosomal structures compared across different cell types uncover the principles underlying the multi-scale organization of chromatin chain: (i) Sub-TADs, TADs, and meta-TADs constitute a robust hierarchical structure. (ii) The assemblies of compartments and TAD-based domains are governed by different organizational principles. (iii) Sub-TADs are the common building blocks of chromosome architecture. Our physically principled interpretation and analysis of Hi-C not only offer an accurate and quantitative view of multi-scale chromatin organization but also help decipher its connections with genome function.

An array of square blocks and checkerboard patterns characteristic to Hi-C data reflects the multi-scale organization of the chromatin chain. Deciphering the structures of chromatin domains from Hi-C and associating them with genome function are open problems of great importance in genome research. However, most existing methods are specialized in finding domains at different scales, making it difficult to integrate the solutions. Here we develop a unified framework for modeling and inferring domain structures over multiple scales, based on a physical model of the chromosome that reflects its nature as a three-dimensional object. Our method efficiently explores the space of domain solutions at different genomic scales, and systematically infers the chromatin domains over multiple scales from Hi-C data by employing a single tuning parameter. Our principled interpretation of Hi-C not only offers a quantitative view of multi-scale chromatin organization but also helps understand its connections with genome function.

Genome-in-a-Box: Building a Chromosome from the Bottom Up

Chromosome structure and dynamics are essential for life, as the way that our genomes are spatially organized within cells is crucial for gene expression, differentiation, and genome transfer to daughter cells. There is a wide variety of methods available to study chromosomes, ranging from live-cell studies to single-molecule biophysics, which we briefly review. While these technologies have yielded a wealth of data, such studies still leave a significant gap between top-down experiments on live cells and bottom-up in vitro single-molecule studies of DNA-protein interactions. Here, we introduce "genome-in-a-box" (GenBox) as an alternative in vitro approach to build and study chromosomes, which bridges this gap. The concept is to assemble a chromosome from the bottom up by taking deproteinated genome-sized DNA isolated from live cells and subsequently add purified DNA-organizing elements, followed by encapsulation in cell-sized containers using microfluidics. Grounded in the rationale of synthetic cell research, the approach would enable to experimentally study emergent effects at the global genome level that arise from the collective action of local DNA-structuring elements. We review the various DNA-structuring elements present in nature, from nucleoid-associated proteins and SMC complexes to phase separation and macromolecular crowders. Finally, we discuss how GenBox can contribute to several open questions on chromosome structure and dynamics.

History-dependent nonequilibrium conformations of a highly confined polymer globule in a sphere

Chromatin undergoes condensation-decondensation processes repeatedly during its cell lifetime. The spatial organization of chromatin in nucleus resembles the fractal globule, of which structure significantly differs from an equilibrium polymer globule. There have been efforts to develop a polymer globule model to describe the fractal globulelike structure of tightly packed chromatin in nucleus. However, the transition pathway of a polymer toward a globular state has been often ignored. Because biological systems are intrinsically in nonequilibrium states, the transition pathway that the chromatin would take before reaching the densely packaged globule should be of importance. In this study, by employing a simple polymer model and Langevin dynamics simulations, we investigate the conformational transition of a single polymer from a swollen coil to a compact globule. We aim to elucidate the effect of transition pathways on the final globular structure. We show that a fast collapse induces a nonequilibrium structure even without a specific intramolecular interaction and that its relaxation toward an equilibrium globule is extremely slow. Due to a strong confinement, the fractal globule never relaxes into an equilibrium state during our simulations such that the globular structure becomes dependent on the transition pathway.

Polymer Modeling Predicts Chromosome Reorganization in Senescence

Lamina-associated domains (LADs) cover a large part of the human genome and are thought to play a major role in shaping the nuclear architectural landscape. Here, we perform polymer simulations, microscopy, and mass spectrometry to dissect the roles played by heterochromatin- and lamina-mediated interactions in nuclear organization. Our model explains the conventional organization of heterochromatin and euchromatin in growing cells and the pathological organization found in oncogene-induced senescence and progeria. We show that the experimentally observed changes in the locality of contacts in senescent and progeroid cells can be explained as arising due to phase transitions in the system. Within our simulations, LADs are highly stochastic, as in experiments. Our model suggests that, once established, the senescent phenotype should be metastable even if lamina-mediated interactions were reinstated. Overall, our simulations uncover a generic physical mechanism that can regulate heterochromatin segregation and LAD formation in a wide range of mammalian nuclei.

A Lamin-Associated Chromatin Model for Chromosome Organization

We propose a simple model for chromatin organization based on the interaction of the chromatin fibers with lamin proteins along the nuclear membrane. Lamin proteins are known to be a major factor that influences chromatin organization and hence gene expression in the cells. We provide a quantitative understanding of lamin-associated chromatin organization in a crowded macromolecular environment by systematically varying the heteropolymer segment distribution and the strength of the lamin-chromatin attractive interaction. Our minimal polymer model reproduces the formation of lamin-associated-domains and provides an in silico tool for quantifying domain length distributions for different distributions of heteropolymer segments. We show that a Gaussian distribution of heteropolymer segments, coupled with strong lamin-chromatin interactions, can qualitatively reproduce observed length distributions of lamin-associated-domains. Further, lamin-mediated interaction can enhance the formation of chromosome territories as well as the organization of chromatin into tightly packed heterochromatin and the loosely packed gene-rich euchromatin regions.

Diffusion of DNA-Binding Species in the Nucleus: A Transient Anomalous Subdiffusion Model

Single-particle tracking experiments have measured escape times of DNA-binding species diffusing in living cells: CRISPR-Cas9, TetR, and LacI. The observed distribution is a truncated power law. Working backward from the experimental results, the observed distribution appears inconsistent with a Gaussian distribution of binding energies. Working forward, the observed distribution leads to transient anomalous subdiffusion, in which diffusion is anomalous at short times and normal at long times, here only mildly anomalous. Monte Carlo simulations are used to characterize the time-dependent diffusion coefficient D(t) in terms of the anomalous exponent α, the crossover time tcross, and the limits D(0) and D(∞) and to relate these quantities to the escape time distribution. The simplest interpretations identify the escape time as the actual binding time to DNA or the period of one-dimensional diffusion on DNA in the standard model combining one-dimensional and three-dimensional search, but a more complicated interpretation may be required. The model has several implications for cell biophysics. 1) The initial anomalous regime represents the search of the DNA-binding species for its target DNA sequence. 2) Non-target DNA sites have a significant effect on search kinetics. False positives in bioinformatic searches of the genome are potentially rate-determining in vivo. For simple binding, the search would be speeded if false-positive sequences were eliminated from the genome. 3) Both binding and obstruction affect diffusion. Obstruction ought to be measured directly, using as the primary probe the DNA-binding species with the binding site inactivated and eGFP as a calibration standard among laboratories and cell types. 4) Overexpression of the DNA-binding species reduces anomalous subdiffusion because the deepest binding sites are occupied and unavailable. 5) The model provides a coarse-grained phenomenological description of diffusion of a DNA-binding species, useful in larger-scale modeling of kinetics, FCS, and FRAP.

Confined crowded polymers near attractive surfaces

We present results from molecular dynamics simulations of a spherically confined neutral polymer in the presence of crowding particles, studying polymer shapes and conformations as a function of the strength of the attraction to the confining wall, solvent quality, and the density of crowders. The conformations of the polymer under good solvent conditions are weakly dependent on crowder particle density, even when the polymer is strongly confined. In contrast, under poor solvent conditions, when the polymer assumes a collapsed conformation when unconfined, it can exhibit transitions to two different adsorbed phases, when either the interaction with the wall or the density of crowder particles is changed. One such transition involves a desorbed collapsed phase change to an adsorbed extended phase as the attraction of the polymer towards the confining wall is increased. Such an adsorbed extended phase can exhibit a second transition to an ordered adsorbed collapsed phase as the crowder particle density is increased. The ordered adsorbed collapsed phase of the polymer differs significantly in its structure from the desorbed collapsed phase. We revisit the earlier understanding of the adsorption of confined polymers on attractive surfaces in light of our results.

Nonequilibrium Biophysical Processes Influence the Large-Scale Architecture of the Cell Nucleus

Model approaches to nuclear architecture have traditionally ignored the biophysical consequences of ATP-fueled active processes acting on chromatin. However, transcription-coupled activity is a source of stochastic forces that are substantially larger than the Brownian forces present at physiological temperatures. Here, we describe an approach to large-scale nuclear architecture in metazoans that incorporates cell-type-specific active processes. The model predicts the statistics of positional distributions, shapes, and overlaps of each chromosome. Simulations of the model reproduce common organizing principles underlying large-scale nuclear architecture across human cell nuclei in interphase. These include the differential positioning of euchromatin and heterochromatin, the territorial organization of chromosomes (including both gene-density-based and size-based chromosome radial positioning schemes), the nonrandom locations of chromosome territories, and the shape statistics of individual chromosomes. We propose that the biophysical consequences of the distribution of transcriptional activity across chromosomes should be central to any chromosome positioning code.

Dynamic Crowding Regulates Transcription

The nuclear environment is highly crowded by biological macromolecules, including chromatin and mobile proteins, which alter the kinetics and efficiency of transcriptional machinery. These alterations have been described, both theoretically and experimentally, for steady-state crowding densities; however, temporal changes in crowding density ("dynamic crowding") have yet to be integrated with gene expression. Dynamic crowding is pertinent to nuclear biology because processes such as chromatin translocation and protein diffusion lend to highly mobile biological crowders. Therefore, to capture such dynamic crowding and investigate its influence on transcription, we employ a three-pronged, systems-molecular approach. A system of chemical reactions represents the transcription pathway, the rates of which are determined by molecular-scale simulations; Brownian dynamics and Monte Carlo simulations quantify protein diffusion and DNA-protein binding affinity, dependent on macromolecular density. Altogether, this approach shows that transcription depends critically on dynamic crowding as the gene expression resultant from dynamic crowding can be profoundly different than that of steady-state crowding. In fact, expression levels can display both amplification and suppression and are notably different for genes or gene populations with different chemical and structural properties. These properties can be exploited to impose circadian expression, which is asymmetric and varies in strength, or to explain expression in cells under biomechanical stress. Therefore, this work demonstrates that dynamic crowding nontrivially alters transcription kinetics and presents dynamic crowding within the bulk nuclear nanoenvironment as a novel regulatory framework for gene expression.

Brownian Dynamics Simulations of Biological Molecules

Brownian dynamics (BD) is a technique for carrying out computer simulations of physical systems that are driven by thermal fluctuations. Biological systems at the macromolecular and cellular level, while falling in the gap between well-established atomic-level models and continuum models, are especially suitable for such simulations. We present a brief history, examples of important biological processes that are driven by thermal motion, and those that have been profitably studied by BD. We also present some of the challenges facing developers of algorithms and software, especially in the attempt to simulate larger systems more accurately and for longer times.

Compressing Θ-Chain in Slit Geometry

When compressed in a slit of width D, a Θ-chain that displays the scaling of size R₀ (diameter) with respect to the number of monomers N, R₀ ∼ aN^1/2, expands in the lateral direction as R_|| ∼ aN^ν(a/D)^2ν-1. Provided that the Θ condition is strictly maintained throughout the compression, the well-known scaling exponent of Θ-chain in two dimensions, ν = 4/7, is anticipated in a perfect confinement. However, numerics shows that upon increasing compression from R₀/D < 1 to R₀/D ≫ 1, ν gradually deviates from ν = 1/2 and plateaus at ν = 3/4, the exponent associated with the self-avoiding walk in two dimensions. Using both theoretical considerations and numerics, we argue that it is highly nontrivial to maintain the Θ condition under confinement because of two major effects. First, as the dimension is reduced from three to two dimensions, the contributions of higher order virial terms, which can be ignored in three dimensions at large N, become significant, making the perturbative expansion used in Flory-type approach inherently problematic. Second and more importantly, the geometrical confinement, which is regarded as an applied external field, alters the second virial coefficient (B₂) changes from B₂ = 0 (Θ condition) in free space to B₂ > 0 (good-solvent condition) in confinement. Our study provides practical insight into how confinement affects the conformation of a single polymer chain.

Heterogeneous Loop Model to Infer 3D Chromosome Structures from Hi-C

Adapting a well-established formalism in polymer physics, we develop a minimalist approach to infer three-dimensional folding of chromatin from Hi-C data. The three-dimensional chromosome structures generated from our heterogeneous loop model (HLM) are used to visualize chromosome organizations that can substantiate the measurements from fluorescence in situ hybridization, chromatin interaction analysis by paired-end tag sequencing, and RNA-seq signals. We demonstrate the utility of the HLM with several case studies. Specifically, the HLM-generated chromosome structures, which reproduce the spatial distribution of topologically associated domains from fluorescence in situ hybridization measurement, show the phase segregation between two types of topologically associated domains explicitly. We discuss the origin of cell-type-dependent gene-expression level by modeling the chromatin globules of α-globin and SOX2 gene loci for two different cell lines. We also use the HLM to discuss how the chromatin folding and gene-expression level of Pax6 loci, associated with mouse neural development, are modulated by interactions with two enhancers. Finally, HLM-generated structures of chromosome 19 of mouse embryonic stem cells, based on single-cell Hi-C data collected over each cell-cycle phase, visualize changes in chromosome conformation along the cell-cycle. Given a contact frequency map between chromatic loci supplied from Hi-C, HLM is a computationally efficient and versatile modeling tool to generate chromosome structures that can complement interpreting other experimental data.

The breakdown of the local thermal equilibrium approximation for a polymer chain during packaging

The conformational relaxation of a polymer chain often slows down in various biological and engineering processes. The polymer, then, may stay in nonequilibrium states throughout the process such that one may not invoke the local thermal equilibrium (LTE) approximation, which has been usually employed to describe the kinetics of various processes. In this work, motivated by recent single-molecule experiments on DNA packaging into a viral capsid, we investigate how the nonequilibrium conformations and the LTE approximation would affect the packaging of a polymer chain into small confinement. We employ a simple but generic coarse-grained model and Langevin dynamics simulations to investigate the packaging kinetics. The polymer segments (both inside and outside the confinement) stay away from equilibrium under strong external force. We devise a simulation scheme to invoke the LTE approximation during packaging and find that the relaxation of nonequilibrium conformations plays a critical role in regulating the packaging rate.

Bacterial chromosome organization. II. Few special cross-links, cell confinement, and molecular crowders play the pivotal roles

Using a coarse-grained bead-spring model of bacterial chromosomes of Caulobacter crescentus and Escherichia coli, we show that just 33 and 38 effective cross-links in 4017 and 4642 monomer chains at special positions along the chain contour can lead to the large-scale organization of the DNA polymer, where confinement effects of the cell walls play a key role in the organization. The positions of the 33/38 cross-links along the chain contour are chosen from the Hi-C contact map of bacteria C. crescentus and E. coli. We represent 1000 base pairs as a coarse-grained monomer in our bead-spring flexible ring polymer model of the DNA polymer. Thus, 4017/4642 beads on a flexible ring polymer represent the C. crescentus/E. coli DNA polymer with 4017/4642 kilo-base pairs. Choosing suitable parameters from Paper I, we also incorporate the role of compaction of the polymer coil due to the presence of molecular crowders and the ability of the chain to release topological constraints. We validate our prediction of the organization of the bacterial chromosomes with available experimental data and also give a prediction of the approximate positions of different segments within the cell. In the absence of confinement, the minimal number of effective cross-links required to organize the DNA chains of 4017/4642 monomers was 60/82 [Agarwal et al., Europhys. Lett. 121, 18004 (2018) and Agarwal et al., J. Phys.: Condens. Matter 30, 034003 (2018)].

Investigating epithelial-to-mesenchymal transition with integrated computational and experimental approaches

The transition between epithelial and mesenchymal (EMT) is a fundamental cellular process that plays critical roles in development, cancer metastasis, and tissue wound healing. EMT is not a binary process but involves multiple partial EMT states that give rise to a high degree of cell state plasticity. Here, we first reviewed several studies on theoretical predictions and experimental verification of these intermediate states, the role of partial EMT on kidney fibrosis development, and how quantitative signaling information controls cell commitment to partial or full EMT upon transient signals. Next, we summarized existing knowledge and open questions on the coupling between EMT and other biological processes, such as the cell cycle, epigenetic regulation, stemness, and apoptosis. Taken together, EMT is a model system that has attracted increasing interests for quantitative experimental and theoretical studies.

Molecular Mechanism of Spontaneous Nucleosome Unraveling

Meters of DNA wrap around histone proteins to form nucleosomes and fit inside the micron-diameter nucleus. For the genetic information encoded in the DNA to become available for transcription, replication, and repair, the DNA-histone assembly must be disrupted. Experiment has indicated that the outer stretches of nucleosomal DNA "breathe" by spontaneously detaching from and reattaching to the histone core. Here, we report direct observation of spontaneous DNA breathing in atomistic molecular dynamics simulations, detailing a microscopic mechanism of the DNA breathing process. According to our simulations, the outer stretches of nucleosomal DNA detach in discrete steps involving 5 or 10 base pairs, with the detachment process being orchestrated by the motion of several conserved histone residues. The inner stretches of nucleosomal DNA are found to be more stably associated with the histone core by more abundant nonspecific DNA-protein contacts, providing a microscopic interpretation of nucleosome unraveling experiments. The CG content of nucleosomal DNA is found to anticorrelate with the extent of unwrapping, supporting the possibility that AT-rich segments may signal the start of transcription by forming less stable nucleosomes.

Chain organization of human interphase chromosome determines the spatiotemporal dynamics of chromatin loci

We investigate spatiotemporal dynamics of human interphase chromosomes by employing a heteropolymer model that incorporates the information of human chromosomes inferred from Hi-C data. Despite considerable heterogeneities in the chromosome structures generated from our model, chromatins are organized into crumpled globules with space-filling (SF) statistics characterized by a single universal scaling exponent (ν = 1/3), and this exponent alone can offer a quantitative account of experimentally observed, many different features of chromosome dynamics. The local chromosome structures, whose scale corresponds to that of topologically associated domains (∼ 0.1 − 1 Mb), display dynamics with a fast relaxation time (≲ 1 − 10 sec); in contrast, the long-range spatial reorganization of the entire chromatin (≳O(102) Mb) occurs on a much slower time scale (≳ hour), providing the dynamic basis of cell-to-cell variability and glass-like behavior of chromosomes. Biological activities, modeled using stronger isotropic white noises added to active loci, accelerate the relaxation dynamics of chromatin domains associated with the low frequency modes and induce phase segregation between the active and inactive loci. Surprisingly, however, they do not significantly change the dynamics at local scales from those obtained under passive conditions. Our study underscores the role of chain organization of chromosome in determining the spatiotemporal dynamics of chromatin loci.

Chromosomes are giant chain molecules made of hundreds of megabase-long DNA intercalated with proteins. Structure and dynamics of interphase chromatin in space and time hold the key to understanding the cell type-dependent gene regulation. In this study, we establish that the crumpled and space-filling (SF) organization of chromatin fiber in the chromosome territory, characterized by a single scaling exponent, is sufficient to explain the complex spatiotemporal hierarchy in chromatin dynamics as well as the subdiffusive motion of the chromatin loci. While seemingly a daunting problem at a first glance, our study shows that relatively simple principles, rooted in polymer physics, can be used to grasp the essence of dynamical properties of the interphase chromatin.

Cell Volume Controls Protein Stability and Compactness of the Unfolded State

Macromolecular crowding is widely accepted as one of the factors that can alter protein stability, structure, and function inside cells. Less often considered is that crowding can be dynamic: as cell volume changes, either as a result of external duress or in the course of the cell cycle, water moves in or out through membrane channels, and crowding changes in tune. Both theory and in vitro experiments predict that protein stability will be altered as a result of crowding changes. However, it is unclear how much the structural ensemble is altered as crowding changes in the cell. To test this, we look at the response of a FRET-labeled kinase to osmotically induced volume changes in live cells. We examine both the folded and unfolded states of the kinase by changing the temperature of the media surrounding the cell. Our data reveals that crowding compacts the structure of its unfolded ensemble but stabilizes the folded protein. We propose that the structure of proteins lacking a rigid, well-defined tertiary structure could be highly sensitive to both increases and decreases in cell volume. Our findings present a possible mechanism for disordered proteins to act as sensors and actuators of cell cycle or external stress events that coincide with a change in macromolecular crowding.

Interphase human chromosome exhibits out of equilibrium glassy dynamics

Fingerprints of the three-dimensional organization of genomes have emerged using advances in Hi-C and imaging techniques. However, genome dynamics is poorly understood. Here, we create the chromosome copolymer model (CCM) by representing chromosomes as a copolymer with two epigenetic loci types corresponding to euchromatin and heterochromatin. Using novel clustering techniques, we establish quantitatively that the simulated contact maps and topologically associating domains (TADs) for chromosomes 5 and 10 and those inferred from Hi-C experiments are in good agreement. Chromatin exhibits glassy dynamics with coherent motion on micron scale. The broad distribution of the diffusion exponents of the individual loci, which quantitatively agrees with experiments, is suggestive of highly heterogeneous dynamics. This is reflected in the cell-to-cell variations in the contact maps. Chromosome organization is hierarchical, involving the formation of chromosome droplets (CDs) on genomic scale, coinciding with the TAD size, followed by coalescence of the CDs, reminiscent of Ostwald ripening.

Computational construction of 3D chromatin ensembles and prediction of functional interactions of alpha-globin locus from 5C data

Conformation capture technologies measure frequencies of interactions between chromatin regions. However, understanding gene-regulation require knowledge of detailed spatial structures of heterogeneous chromatin in cells. Here we describe the nC-SAC (n-Constrained-Self Avoiding Chromatin) method that transforms experimental interaction frequencies into 3D ensembles of chromatin chains. nC-SAC first distinguishes specific from non-specific interaction frequencies, then generates 3D chromatin ensembles using identified specific interactions as spatial constraints. Application to α-globin locus shows that these constraints (∼20%) drive the formation of ∼99% all experimentally captured interactions, in which ∼30% additional to the imposed constraints is found to be specific. Many novel specific spatial contacts not captured by experiments are also predicted. A subset, of which independent ChIA-PET data are available, is validated to be RNAPII-, CTCF-, and RAD21-mediated. Their positioning in the architectural context of imposed specific interactions from nC-SAC is highly important. Our results also suggest the presence of a many-body structural unit involving α-globin gene, its enhancers, and POL3RK gene for regulating the expression of α-globin in silent cells.

DNA Internal Motion Likely Accelerates Protein Target Search in a Packed Nucleoid

Transcription factors must diffuse through densely packed and coiled DNA to find their binding sites. Using a coarse-grained model of DNA and lac repressor (LacI) in the Escherichia coli nucleoid, simulations were performed to examine how LacI diffuses in such a space. Despite the canonical picture of LacI diffusing rather freely, in reality the DNA is densely packed, is not rigid but highly mobile, and the dynamics of DNA dictates to a great extent the LacI motion. A possibly better picture of unbound LacI motion is that of gated diffusion, where DNA confines LacI in a cage, but LacI can move between cages when hindering DNA strands move out of the way. Three-dimensional diffusion constants for unbound LacI computed from simulations closely match those for unbound LacI in vivo reported in the literature. The internal motions of DNA appear to be governed by strong internal forces arising from being crowded into the small space of the nucleoid. A consequence of the DNA internal motion is that protein target search may be accelerated.

Spatial organization of the budding yeast genome in the cell nucleus and identification of specific chromatin interactions from multi-chromosome constrained chromatin model

Nuclear landmarks and biochemical factors play important roles in the organization of the yeast genome. The interaction pattern of budding yeast as measured from genome-wide 3C studies are largely recapitulated by model polymer genomes subject to landmark constraints. However, the origin of inter-chromosomal interactions, specific roles of individual landmarks, and the roles of biochemical factors in yeast genome organization remain unclear. Here we describe a multi-chromosome constrained self-avoiding chromatin model (mC-SAC) to gain understanding of the budding yeast genome organization. With significantly improved sampling of genome structures, both intra- and inter-chromosomal interaction patterns from genome-wide 3C studies are accurately captured in our model at higher resolution than previous studies. We show that nuclear confinement is a key determinant of the intra-chromosomal interactions, and centromere tethering is responsible for the inter-chromosomal interactions. In addition, important genomic elements such as fragile sites and tRNA genes are found to be clustered spatially, largely due to centromere tethering. We uncovered previously unknown interactions that were not captured by genome-wide 3C studies, which are found to be enriched with tRNA genes, RNAPIII and TFIIS binding. Moreover, we identified specific high-frequency genome-wide 3C interactions that are unaccounted for by polymer effects under landmark constraints. These interactions are enriched with important genes and likely play biological roles.

The architecture of the cell nucleus and the spatial organization of the genome are important in determining nuclear functions. Single-cell imaging techniques and chromosome conformation capture (3C) based methods have provided a wealth of information on the spatial organization of chromosomes. Here we describe a multi-chromosome ensemble model of chromatin chains for understanding the folding principles of budding yeast genome. By overcoming severe challenges in sampling self-avoiding chromatin chains in nuclear confinement, we succeed in generating a large number of model genomes of budding yeast. Our model predicts chromatin interactions that have good correlation with experimental measurements. Our results showed that the spatial confinement of cell nucleus and excluded-volume effect are key determinants of the folding behavior of yeast chromosomes, and largely account for the observed intra-chromosomal interactions. Furthermore, we determined the specific roles of individual nuclear landmarks and biochemical factors, and our analysis showed that centromere tethering largely determines inter-chromosomal interactions. In addition, we were able to infer biological properties from the organization of modeled genomes. We found that the spatial locations of important elements such as fragile sites and tRNA genes are largely determined by the tethering of centromeres to the Spindle Pole Body. We further showed that many of these spatial locations can be predicted by using the genomic distances to the centromeres. Overall, our results revealed important insight into the organizational principles of the budding yeast genome and predicted a number of important biological findings that are fully experimentally testable.

Contact Statistics Highlight Distinct Organizing Principles of Proteins and RNA

Although both RNA and proteins have densely packed native structures, chain organizations of these two biopolymers are fundamentally different. Motivated by the recent discoveries in chromatin folding that interphase chromosomes have territorial organization with signatures pointing to metastability, we analyzed the biomolecular structures deposited in the Protein Data Bank and found that the intrachain contact probabilities, P(s) as a function of the arc length s, decay in power-law ∼s(-γ) over the intermediate range of s, 10 ≲ s ≲ 110. We found that the contact probability scaling exponent is γ ≈ 1.11 for large RNA (N > 110), γ ≈ 1.41 for small-sized RNA (N < 110), and γ ≈ 1.65 for proteins. Given that Gaussian statistics is expected for a fully equilibrated chain in polymer melts, the deviation of γ-value from γ = 1.5 for the subchains of large RNA in the native state suggests that the chain configuration of RNA is not fully equilibrated. It is visually clear that folded structures of large-sized RNA (N ≳ 110) adopt crumpled structures, partitioned into modular multidomains assembled by proximal sequences along the chain, whereas the polypeptide chain of folded proteins looks better mixed with the rest of the structure. Our finding of γ ≈ 1 for large RNA might be an ineluctable consequence of the hierarchical ordering of the secondary to tertiary elements in the folding process.

Physical properties of the chromosomes and implications for development

Remarkable progress has been made in understanding chromosome structures inside the cell nucleus. Recent advances in Hi-C technologies enable the detection of genome-wide chromatin interactions, providing insight into three-dimensional (3D) genome organization. Advancements in the spatial and temporal resolutions of imaging as well as in molecular biological techniques allow the tracking of specific chromosomal loci, improving our understanding of chromosome movements. From these data, we are beginning to understand how the intra-nuclear locations of chromatin loci and the 3D genome structure change during development and differentiation. This emerging field of genome structure and dynamics research requires an interdisciplinary approach including efficient collaborations between experimental biologists and physicists, informaticians, or engineers. Quantitative and mathematical analyses based on polymer physics are becoming increasingly important for processing and interpreting experimental data on 3D chromosome structures and dynamics. In this review, we aim to provide an overview of recent research on the physical aspects of chromosome structure and dynamics oriented for biologists. These studies have mainly focused on chromosomes at the cellular level, using unicellular organisms and cultured cells. However, physical parameters that change during development, such as nuclear size, may impact genome structure and dynamics. Here, we discuss how chromatin dynamics and genome structures in early embryos change during development, which we expect will be a hot topic in the field of chromatin dynamics in the near future. We hope this review helps developmental biologists to quantitatively investigate the physical natures of chromosomes in developmental biology research.

As a Nucleus Enters a Small Pore, Chromatin Stretches and Maintains Integrity, Even with DNA Breaks

As a cell pushes or pulls its nucleus through a small constriction, the chromatin must distort and somehow maintain genomic stability despite ever-present double-strand breaks in the DNA. Here we visualize within a living cell the pore-size dependent deformation of a specific locus engineered into chromosome-1 and cleaved. An mCherry-tagged nuclease targets the submicron locus, causing DNA cleavage and recruiting repair factors such as GFP-53BP1 to a large region around the locus. Aspiration of a cell and its nucleus into a micropipette shows that chromatin aligns and stretches parallel to the pore. Extension is largest in small pores, increasing >10-fold but remaining 30-fold shorter than the DNA contour length in the locus. Brochard and de Gennes' blob model for tube geometry fits the data, with a simple modification for chromatin crowding. Continuity of the highly extended, cleaved chromatin is also maintained, consistent with folding and cross bridging of the DNA. Surprisingly, extensional integrity is unaffected by an inhibitor of the DNA repair scaffold.

Shape Transitions and Chiral Symmetry Breaking in the Energy Landscape of the Mitotic Chromosome

We derive an unbiased information theoretic energy landscape for chromosomes at metaphase using a maximum entropy approach that accurately reproduces the details of the experimentally measured pairwise contact probabilities between genomic loci. Dynamical simulations using this landscape lead to cylindrical, helically twisted structures reflecting liquid crystalline order. These structures are similar to those arising from a generic ideal homogenized chromosome energy landscape. The helical twist can be either right or left handed so chiral symmetry is broken spontaneously. The ideal chromosome landscape when augmented by interactions like those leading to topologically associating domain formation in the interphase chromosome reproduces these behaviors. The phase diagram of this landscape shows that the helical fiber order and the cylindrical shape persist at temperatures above the onset of chiral symmetry breaking, which is limited by the topologically associating domain interaction strength.

A topologically driven glass in ring polymers

The static and dynamic properties of ring polymers in concentrated solutions remains one of the last deep unsolved questions in polymer physics. At the same time, the nature of the glass transition in polymeric systems is also not well understood. In this work, we study a novel glass transition in systems made of circular polymers by exploiting the topological constraints that are conjectured to populate concentrated solutions of rings. We show that such rings strongly interpenetrate through one another, generating an extensive network of topological interactions that dramatically affects their dynamics. We show that a kinetically arrested state can be induced by randomly pinning a small fraction of the rings. This occurs well above the classical glass transition temperature at which microscopic mobility is lost. Our work both demonstrates the existence of long-lived inter-ring penetrations and realizes a novel, topologically induced, glass transition.

Multiscale Modeling of Cellular Epigenetic States: Stochasticity in Molecular Networks, Chromatin Folding in Cell Nuclei, and Tissue Pattern Formation of Cells

Genome sequences provide the overall genetic blueprint of cells, but cells possessing the same genome can exhibit diverse phenotypes. There is a multitude of mechanisms controlling cellular epigenetic states and that dictate the behavior of cells. Among these, networks of interacting molecules, often under stochastic control, depending on the specific wirings of molecular components and the physiological conditions, can have a different landscape of cellular states. In addition, chromosome folding in three-dimensional space provides another important control mechanism for selective activation and repression of gene expression. Fully differentiated cells with different properties grow, divide, and interact through mechanical forces and communicate through signal transduction, resulting in the formation of complex tissue patterns. Developing quantitative models to study these multi-scale phenomena and to identify opportunities for improving human health requires development of theoretical models, algorithms, and computational tools. Here we review recent progress made in these important directions.