A minimal G$\bar{\text{o}}$-model for rebuilding whole genome structures from haploid single-cell Hi-C data

Abstract

We present a minimal computational model, which allows very fast, on-the-fly construction of three dimensional haploid interphase genomes from single cell Hi-C contact maps using the HOOMD-blue molecular dynamics package on graphics processing units. Chromosomes are represented by a string of connected beads, each of which corresponds to 100,000 base pairs, and contacts are mediated via a structure-based harmonic potential. We suggest and test two minimization protocols which consistently fold into conformationally similar low energy states. The latter are similar to previously published structures but are calculated in a fraction of the time. We find evidence that mere fulfillment of contact maps is insufficient to create experimentally relevant structures. Particularly, an excluded volume term is required in our model to induce the formation of chromosome territories. We also observe empirically that contact maps do not capture the chirality of the underlying structures. Depending on starting configurations and protocol details, one of two mirror images emerges. Finally, we analyze the occurrence of knots in a particular chromosome. The same knot appears in (almost) all structures irrespective of minimization protocols or even details of underlying potentials providing further evidence for the existence of knots in interphase chromatin.

arXiv comment18 pages, 7 figures. The updated submission contains additional explanations (for example) to the Go model and why we use this in the context of chromatin, we address the stopping criterion. References have been added as a result to an email received pointing us to publications we had not read before

References in corpus (5)

• Strong scaling of general-purpose molecular dynamics simulations on GPUs
• Equilibration of Long Chain Polymer Melts in Computer Simulations
• Hi-C-constrained physical models of human chromosomes recover functionally-related properties of genome organization
• Shape Transitions and Chiral Symmetry Breaking in the Energy Landscape of the Mitotic Chromosome
• Colocalization of coregulated genes: a steered molecular dynamics study of human chromosome 19

Cited by in corpus (1)

• Single-cell 3D genome reconstruction in the haploid setting using rigidity theory