1
|
Korsak S, Banecki K, Plewczynski D. Multiscale molecular modeling of chromatin with MultiMM: From nucleosomes to the whole genome. Comput Struct Biotechnol J 2024; 23:3537-3548. [PMID: 39435339 PMCID: PMC11492436 DOI: 10.1016/j.csbj.2024.09.025] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/09/2024] [Revised: 09/28/2024] [Accepted: 09/29/2024] [Indexed: 10/23/2024] Open
Abstract
Motivation: We present a user-friendly 3D chromatin simulation model for the human genome based on OpenMM, addressing the challenges posed by existing models with use-specific implementations. Our approach employs a multi-scale energy minimization strategy, capturing chromatin's hierarchical structure. Initiating with a Hilbert curve-based structure, users can input files specifying nucleosome positioning, loops, compartments, or subcompartments. Results: The model utilizes an energy minimization approach with a large choice of numerical integrators, providing the entire genome's structure within minutes. Output files include the generated structures for each chromosome, offering a versatile and accessible tool for chromatin simulation in bioinformatics studies. Furthermore, MultiMM is capable of producing nucleosome-resolution structures by making simplistic geometric assumptions about the structure and the density of nucleosomes on the DNA. Code availability: Open-source software and the manual are freely available on https://github.com/SFGLab/MultiMM or via pip https://pypi.org/project/MultiMM/.
Collapse
Affiliation(s)
- Sevastianos Korsak
- Faculty of Mathematics and Information Science, Warsaw University of Technology, Warsaw, Poland
- Centre of New Technologies, University of Warsaw, Warsaw, Poland
| | - Krzysztof Banecki
- Faculty of Mathematics and Information Science, Warsaw University of Technology, Warsaw, Poland
- Centre of New Technologies, University of Warsaw, Warsaw, Poland
| | - Dariusz Plewczynski
- Faculty of Mathematics and Information Science, Warsaw University of Technology, Warsaw, Poland
- Centre of New Technologies, University of Warsaw, Warsaw, Poland
| |
Collapse
|
2
|
Dasgupta M, Guha S, Armbruster L, Das D, Mitra MK. Nature of barriers determines first passage times in heterogeneous media. SOFT MATTER 2024. [PMID: 39318347 DOI: 10.1039/d4sm00908h] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 09/26/2024]
Abstract
Intuition suggests that passage times across a region increase with the number of barriers along the path. Can this fail depending on the nature of the barrier? To probe this fundamental question, we exactly solve for the first passage time in general d-dimensions for diffusive transport through a spatially patterned array of obstacles - either entropic or energetic, depending on the nature of the obstacles. For energetic barriers, we show that first passage times vary non-monotonically with the number of barriers, while for entropic barriers it increases monotonically. This non-monotonicity for energetic barriers is further reflected in the behaviour of effective diffusivity as well. We then design a simple experiment where a robotic bug navigates in a heterogeneous environment through a spatially patterned array of obstacles to validate our predictions. Finally, using numerical simulations, we show that this non-monotonic behaviour for energetic barriers is general and extends to even super-diffusive transport.
Collapse
Affiliation(s)
| | - Sougata Guha
- Department of Physics, IIT Bombay, Mumbai 400076, India.
- INFN Napoli, Complesso Universitario di Monte S. Angelo, 80126 Napoli, Italy
| | | | - Dibyendu Das
- Department of Physics, IIT Bombay, Mumbai 400076, India.
| | - Mithun K Mitra
- Department of Physics, IIT Bombay, Mumbai 400076, India.
| |
Collapse
|
3
|
Korsak S, Plewczynski D. LoopSage: An energy-based Monte Carlo approach for the loop extrusion modeling of chromatin. Methods 2024; 223:106-117. [PMID: 38295892 DOI: 10.1016/j.ymeth.2024.01.015] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/01/2023] [Revised: 12/29/2023] [Accepted: 01/10/2024] [Indexed: 02/05/2024] Open
Abstract
The connection between the patterns observed in 3C-type experiments and the modeling of polymers remains unresolved. This paper presents a simulation pipeline that generates thermodynamic ensembles of 3D structures for topologically associated domain (TAD) regions by loop extrusion model (LEM). The simulations consist of two main components: a stochastic simulation phase, employing a Monte Carlo approach to simulate the binding positions of cohesins, and a dynamical simulation phase, utilizing these cohesins' positions to create 3D structures. In this approach, the system's total energy is the combined result of the Monte Carlo energy and the molecular simulation energy, which are iteratively updated. The structural maintenance of chromosomes (SMC) protein complexes are represented as loop extruders, while the CCCTC-binding factor (CTCF) locations on DNA sequence are modeled as energy minima on the Monte Carlo energy landscape. Finally, the spatial distances between DNA segments from ChIA-PET experiments are compared with the computer simulations, and we observe significant Pearson correlations between predictions and the real data. LoopSage model offers a fresh perspective on chromatin loop dynamics, allowing us to observe phase transition between sparse and condensed states in chromatin.
Collapse
Affiliation(s)
- Sevastianos Korsak
- Faculty of Mathematics and Information Science, Warsaw University of Technology, Warsaw, Poland; Center of New Technologies, University of Warsaw, Warsaw, Poland
| | - Dariusz Plewczynski
- Faculty of Mathematics and Information Science, Warsaw University of Technology, Warsaw, Poland; Center of New Technologies, University of Warsaw, Warsaw, Poland.
| |
Collapse
|
4
|
Agarwal A, Korsak S, Choudhury A, Plewczynski D. The dynamic role of cohesin in maintaining human genome architecture. Bioessays 2023; 45:e2200240. [PMID: 37603403 DOI: 10.1002/bies.202200240] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/09/2022] [Revised: 08/03/2023] [Accepted: 08/07/2023] [Indexed: 08/22/2023]
Abstract
Recent advances in genomic and imaging techniques have revealed the complex manner of organizing billions of base pairs of DNA necessary for maintaining their functionality and ensuring the proper expression of genetic information. The SMC proteins and cohesin complex primarily contribute to forming higher-order chromatin structures, such as chromosomal territories, compartments, topologically associating domains (TADs) and chromatin loops anchored by CCCTC-binding factor (CTCF) protein or other genome organizers. Cohesin plays a fundamental role in chromatin organization, gene expression and regulation. This review aims to describe the current understanding of the dynamic nature of the cohesin-DNA complex and its dependence on cohesin for genome maintenance. We discuss the current 3C technique and numerous bioinformatics pipelines used to comprehend structural genomics and epigenetics focusing on the analysis of Cohesin-centred interactions. We also incorporate our present comprehension of Loop Extrusion (LE) and insights from stochastic modelling.
Collapse
Affiliation(s)
- Abhishek Agarwal
- Centre of New Technologies, University of Warsaw, Warsaw, Poland
| | - Sevastianos Korsak
- Faculty of Mathematics and Information Science, Warsaw University of Technology, Warsaw, Poland
| | | | - Dariusz Plewczynski
- Centre of New Technologies, University of Warsaw, Warsaw, Poland
- Faculty of Mathematics and Information Science, Warsaw University of Technology, Warsaw, Poland
| |
Collapse
|
5
|
Mitra D, Pande S, Chatterji A. Polymer architecture orchestrates the segregation and spatial organization of replicating E. coli chromosomes in slow growth. SOFT MATTER 2022; 18:5615-5631. [PMID: 35861071 DOI: 10.1039/d2sm00734g] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/15/2023]
Abstract
The mechanism of chromosome segregation and organization in the bacterial cell cycle of E. coli is one of the least understood aspects in its life cycle. The E. coli chromosome is often modelled as a bead spring ring polymer. We introduce cross-links in the DNA-ring polymer, resulting in the formation of loops within each replicating bacterial chromosome. We use simulations to show that the chosen polymer-topology ensures its self-organization along the cell long-axis, such that various chromosomal loci get spatially localized as seen in vivo. The localization of loci arises due to entropic repulsion between polymer loops within each daughter DNA confined in a cylinder. The cellular addresses of the loci in our model are in fair agreement with those seen in experiments as given in J. A. Cass et al., Biophys. J., 2016, 110, 2597-2609. We also show that the adoption of such modified polymer architectures by the daughter DNAs leads to an enhanced propensity of their spatial segregation. Secondly, we match other experimentally reported results, including observation of the cohesion time and the ter-transition. Additionally, the contact map generated from our simulations reproduces the macro-domain like organization as seen in the experimentally obtained Hi-C map. Lastly, we have also proposed a plausible reconciliation of the 'Train Track' and the 'Replication Factory' models which provide conflicting descriptions of the spatial organization of the replication forks. Thus, we reconcile observations from complementary experimental techniques probing bacterial chromosome organization.
Collapse
|
6
|
Yan H, Surovtsev I, Williams JF, Bailey MLP, King MC, Mochrie SGJ. Extrusion of chromatin loops by a composite loop extrusion factor. Phys Rev E 2021; 104:024414. [PMID: 34525654 DOI: 10.1103/physreve.104.024414] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/02/2021] [Accepted: 06/16/2021] [Indexed: 12/26/2022]
Abstract
Chromatin loop extrusion by structural maintenance of chromosome (SMC) complexes is thought to underlie intermediate-scale chromatin organization inside cells. Motivated by a number of experiments suggesting that nucleosomes may block loop extrusion by SMCs, such as cohesin and condensin complexes, we introduce and characterize theoretically a composite loop extrusion factor (composite LEF) model. In addition to an SMC complex that creates a chromatin loop by encircling two threads of DNA, this model includes a remodeling complex that relocates or removes nucleosomes as it progresses along the chromatin, and nucleosomes that block SMC translocation along the DNA. Loop extrusion is enabled by SMC motion along nucleosome-free DNA, created in the wake of the remodeling complex, while nucleosome rebinding behind the SMC acts as a ratchet, holding the SMC close to the remodeling complex. We show that, for a wide range of parameter values, this collection of factors constitutes a composite LEF that extrudes loops with a velocity, comparable to the velocity of remodeling complex translocation on chromatin in the absence of SMC, and much faster than loop extrusion by an isolated SMC that is blocked by nucleosomes.
Collapse
Affiliation(s)
- Hao Yan
- Integrated Graduate Program in Physical and Engineering Biology, Yale University, New Haven, Connecticut 06511, USA.,Department of Physics, Yale University, New Haven, Connecticut 06511, USA
| | - Ivan Surovtsev
- Department of Physics, Yale University, New Haven, Connecticut 06511, USA.,Department of Cell Biology, Yale School of Medicine, New Haven, Connecticut 06520, USA
| | - Jessica F Williams
- Department of Cell Biology, Yale School of Medicine, New Haven, Connecticut 06520, USA
| | - Mary Lou P Bailey
- Integrated Graduate Program in Physical and Engineering Biology, Yale University, New Haven, Connecticut 06511, USA.,Department of Applied Physics, Yale University, New Haven, Connecticut 06511, USA
| | - Megan C King
- Integrated Graduate Program in Physical and Engineering Biology, Yale University, New Haven, Connecticut 06511, USA.,Department of Cell Biology, Yale School of Medicine, New Haven, Connecticut 06520, USA.,Department of Molecular, Cellular, and Developmental Biology, Yale University, New Haven, Connecticut 06511, USA
| | - Simon G J Mochrie
- Integrated Graduate Program in Physical and Engineering Biology, Yale University, New Haven, Connecticut 06511, USA.,Department of Physics, Yale University, New Haven, Connecticut 06511, USA.,Department of Applied Physics, Yale University, New Haven, Connecticut 06511, USA
| |
Collapse
|