Perspectives for the reconstruction of 3D chromatin conformation using single cell Hi-C data

How fork-length asymmetry affects solvent connectivity and diffusion in grafted polymeric model membranes

The hydrophilic pore morphology and solvent diffusion within model (amphiphilic) polymer membranes are simulated by dissipative particle dynamics (DPD). The polymers are composed of a backbone of 18 covalently bonded A beads to which at regular intervals side chains are attached. The side chains are composed of linear Ap chains (i.e., -A1-A2…Ap) from which two branches, [AsC] and [ArC], split off (s ≤ r). C beads serve as functionalized hydrophilic pendent sites. The branch lengths (s + 1 and r + 1) are varied. Five repeat unit designs (with general formula A3[Ap[AsC][ArC]]) are considered: A2[A3C][A3C] (symmetric branching), A2[A2C][A4C], A2[AC][A5C], A2[C][A6C] (highly asymmetric branching), and A4[AC][A3C]. The distribution of water (W) and W diffusion through nanophase segregated hydrophilic pores is studied. For similar primary length p, an increase in side chain symmetry favors hydrophilic pore connectivity and long-range water transport. C beads located on the longer [ArC] branches reveal the highest C bead mobility and are more strongly associated with water than the C beads on the shorter [AsC] branches. The connectivity of hydrophilic (W and W + C) phases through mapped replica of selected snapshots obtained from Monte Carlo tracer diffusion simulations is in line with trends found from the W bead diffusivities during DPD simulations. The diffusive pathways for protons (H+) in proton exchange membranes and for hydronium (OH-) in anion exchange membranes are the same as for solvents. Therefore, control of the side chain architecture is an interesting design parameter for optimizing membrane conductivities.

Replicating Chromosomes in Whole-Cell Models of Bacteria

Computational models of cells cannot be considered complete unless they include the most fundamental process of life, the replication of genetic material. In a recent study, we presented a computational framework to model systems of replicating bacterial chromosomes as polymers at 10 bp resolution with Brownian dynamics. This approach was used to investigate changes in chromosome organization during replication and extend the applicability of an existing whole-cell model (WCM) for a genetically minimal bacterium, JCVI-syn3A, to the entire cell cycle. To achieve cell-scale chromosome structures that are realistic, we modeled the chromosome as a self-avoiding homopolymer with bending and torsional stiffnesses that capture the essential mechanical properties of dsDNA in Syn3A. Additionally, the polymer interacts with ribosomes distributed according to cryo-electron tomograms of Syn3A. The polymer model was further augmented by computational models of loop extrusion by structural maintenance of chromosomes (SMC) protein complexes and topoisomerase action, and the modeling and analysis of multi-fork replication states.

[Advances in methods and applications of single-cell Hi-C data analysis]

Chromatin three-dimensional genome structure plays a key role in cell function and gene regulation. Single-cell Hi-C techniques can capture genomic structure information at the cellular level, which provides an opportunity to study changes in genomic structure between different cell types. Recently, some excellent computational methods have been developed for single-cell Hi-C data analysis. In this paper, the available methods for single-cell Hi-C data analysis were first reviewed, including preprocessing of single-cell Hi-C data, multi-scale structure recognition based on single-cell Hi-C data, bulk-like Hi-C contact matrix generation based on single-cell Hi-C data sets, pseudo-time series analysis, and cell classification. Then the application of single-cell Hi-C data in cell differentiation and structural variation was described. Finally, the future development direction of single-cell Hi-C data analysis was also prospected.

Dynamics of chromosome organization in a minimal bacterial cell

Computational models of cells cannot be considered complete unless they include the most fundamental process of life, the replication and inheritance of genetic material. By creating a computational framework to model systems of replicating bacterial chromosomes as polymers at 10 bp resolution with Brownian dynamics, we investigate changes in chromosome organization during replication and extend the applicability of an existing whole-cell model (WCM) for a genetically minimal bacterium, JCVI-syn3A, to the entire cell-cycle. To achieve cell-scale chromosome structures that are realistic, we model the chromosome as a self-avoiding homopolymer with bending and torsional stiffnesses that capture the essential mechanical properties of dsDNA in Syn3A. In addition, the conformations of the circular DNA must avoid overlapping with ribosomes identitied in cryo-electron tomograms. While Syn3A lacks the complex regulatory systems known to orchestrate chromosome segregation in other bacteria, its minimized genome retains essential loop-extruding structural maintenance of chromosomes (SMC) protein complexes (SMC-scpAB) and topoisomerases. Through implementing the effects of these proteins in our simulations of replicating chromosomes, we find that they alone are sufficient for simultaneous chromosome segregation across all generations within nested theta structures. This supports previous studies suggesting loop-extrusion serves as a near-universal mechanism for chromosome organization within bacterial and eukaryotic cells. Furthermore, we analyze ribosome diffusion under the influence of the chromosome and calculate in silico chromosome contact maps that capture inter-daughter interactions. Finally, we present a methodology to map the polymer model of the chromosome to a Martini coarse-grained representation to prepare molecular dynamics models of entire Syn3A cells, which serves as an ultimate means of validation for cell states predicted by the WCM.

Biology and Model Predictions of the Dynamics and Heterogeneity of Chromatin-Nuclear Lamina Interactions

Associations of chromatin with the nuclear lamina, at the nuclear periphery, help shape the genome in 3 dimensions. The genomic landscape of lamina-associated domains (LADs) is well characterized, but much remains unknown on the physical and mechanistic properties of chromatin conformation at the nuclear lamina. Computational models of chromatin folding at, and interactions with, a surface representing the nuclear lamina are emerging in attempts to characterize these properties and predict chromatin behavior at the lamina in health and disease. Here, we highlight the heterogeneous nature of the nuclear lamina and LADs, outline the main 3-dimensional chromatin structural modeling methods, review applications of modeling chromatin-lamina interactions and discuss biological insights inferred from these models in normal and disease states. Lastly, we address perspectives on future developments in modeling chromatin interactions with the nuclear lamina.

Super-resolution visualization of chromatin loop folding in human lymphoblastoid cells using interferometric photoactivated localization microscopy

The three-dimensional (3D) genome structure plays a fundamental role in gene regulation and cellular functions. Recent studies in 3D genomics inferred the very basic functional chromatin folding structures known as chromatin loops, the long-range chromatin interactions that are mediated by protein factors and dynamically extruded by cohesin. We combined the use of FISH staining of a very short (33 kb) chromatin fragment, interferometric photoactivated localization microscopy (iPALM), and traveling salesman problem-based heuristic loop reconstruction algorithm from an image of the one of the strongest CTCF-mediated chromatin loops in human lymphoblastoid cells. In total, we have generated thirteen good quality images of the target chromatin region with 2–22 nm oligo probe localization precision. We visualized the shape of the single chromatin loops with unprecedented genomic resolution which allowed us to study the structural heterogeneity of chromatin looping. We were able to compare the physical distance maps from all reconstructed image-driven computational models with contact frequencies observed by ChIA-PET and Hi-C genomic-driven methods to examine the concordance between single cell imaging and population based genomic data.

Persistent spectral simplicial complex-based machine learning for chromosomal structural analysis in cellular differentiation

The three-dimensional (3D) chromosomal structure plays an essential role in all DNA-templated processes, including gene transcription, DNA replication and other cellular processes. Although developing chromosome conformation capture (3C) methods, such as Hi-C, which can generate chromosomal contact data characterized genome-wide chromosomal structural properties, understanding 3D genomic nature-based on Hi-C data remains lacking. Here, we propose a persistent spectral simplicial complex (PerSpectSC) model to describe Hi-C data for the first time. Specifically, a filtration process is introduced to generate a series of nested simplicial complexes at different scales. For each of these simplicial complexes, its spectral information can be calculated from the corresponding Hodge Laplacian matrix. PerSpectSC model describes the persistence and variation of the spectral information of the nested simplicial complexes during the filtration process. Different from all previous models, our PerSpectSC-based features provide a quantitative global-scale characterization of chromosome structures and topology. Our descriptors can successfully classify cell types and also cellular differentiation stages for all the 24 types of chromosomes simultaneously. In particular, persistent minimum best characterizes cell types and Dim (1) persistent multiplicity best characterizes cellular differentiation. These results demonstrate the great potential of our PerSpectSC-based models in polymeric data analysis.