1
|
Xu Y, Das P, McCord RP, Shen T. Node features of chromosome structure networks and their connections to genome annotation. Comput Struct Biotechnol J 2024; 23:2240-2250. [PMID: 38827231 PMCID: PMC11140560 DOI: 10.1016/j.csbj.2024.05.026] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/13/2023] [Revised: 05/14/2024] [Accepted: 05/15/2024] [Indexed: 06/04/2024] Open
Abstract
The 3D conformations of chromosomes can encode biological significance, and the implications of such structures have been increasingly appreciated recently. Certain chromosome structural features, such as A/B compartmentalization, are frequently extracted from Hi-C pairwise genome contact information (physical association between different regions of the genome) and compared with linear annotations of the genome, such as histone modifications and lamina association. We investigate how additional properties of chromosome structure can be deduced using an abstract graph representation of the contact heatmap, and describe specific network properties that can have a strong connection with some of these biological annotations. We constructed chromosome structure networks (CSNs) from bulk Hi-C data and calculated a set of site-resolved (node-based) network properties. These properties are useful for characterizing certain aspects of chromosomal structure. We examined the ability of network properties to differentiate several scenarios, such as haploid vs diploid cells, partially inverted nuclei vs conventional architecture, depletion of chromosome architectural proteins, and structural changes during cell development. We also examined the connection between network properties and a series of other linear annotations, such as histone modifications and chromatin states including poised promoter and enhancer labels. We found that semi-local network properties exhibit greater capability in characterizing genome annotations compared to diffusive or ultra-local node features. For example, the local square clustering coefficient can be a strong classifier of lamina-associated domains. We demonstrated that network properties can be useful for highlighting large-scale chromosome structure differences that emerge in different biological situations.
Collapse
Affiliation(s)
- Yingjie Xu
- Graduate School of Genome Science & Technology, University of Tennessee, Knoxville, TN 37996, USA
| | - Priyojit Das
- Department of Genetics, Harvard Medical School, Boston, MA 02115, USA
- Department of Molecular Biology, Massachusetts General Hospital, Boston, MA 02114, USA
| | - Rachel Patton McCord
- Department of Biochemistry & Cellular and Molecular Biology, University of Tennessee, Knoxville, TN 37996, USA
| | - Tongye Shen
- Department of Biochemistry & Cellular and Molecular Biology, University of Tennessee, Knoxville, TN 37996, USA
| |
Collapse
|
2
|
Conte M, Abraham A, Esposito A, Yang L, Gibcus JH, Parsi KM, Vercellone F, Fontana A, Di Pierno F, Dekker J, Nicodemi M. Polymer Physics Models Reveal Structural Folding Features of Single-Molecule Gene Chromatin Conformations. Int J Mol Sci 2024; 25:10215. [PMID: 39337699 PMCID: PMC11432541 DOI: 10.3390/ijms251810215] [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: 07/14/2024] [Revised: 09/17/2024] [Accepted: 09/22/2024] [Indexed: 09/30/2024] Open
Abstract
Here, we employ polymer physics models of chromatin to investigate the 3D folding of a 2 Mb wide genomic region encompassing the human LTN1 gene, a crucial DNA locus involved in key cellular functions. Through extensive Molecular Dynamics simulations, we reconstruct in silico the ensemble of single-molecule LTN1 3D structures, which we benchmark against recent in situ Hi-C 2.0 data. The model-derived single molecules are then used to predict structural folding features at the single-cell level, providing testable predictions for super-resolution microscopy experiments.
Collapse
Affiliation(s)
- Mattia Conte
- Dipartimento di Fisica, Università di Napoli Federico II, and INFN Napoli, Complesso Universitario di Monte Sant'Angelo, 80126 Naples, Italy
| | - Alex Abraham
- Dipartimento di Fisica, Università di Napoli Federico II, and INFN Napoli, Complesso Universitario di Monte Sant'Angelo, 80126 Naples, Italy
| | - Andrea Esposito
- Dipartimento di Fisica, Università di Napoli Federico II, and INFN Napoli, Complesso Universitario di Monte Sant'Angelo, 80126 Naples, Italy
| | - Liyan Yang
- Department of Systems Biology, University of Massachusetts Chan Medical School, Worcester, MA 01605, USA
| | - Johan H Gibcus
- Department of Systems Biology, University of Massachusetts Chan Medical School, Worcester, MA 01605, USA
| | - Krishna M Parsi
- Diabetes Center of Excellence and Program in Molecular Medicine, University of Massachusetts Chan Medical School, Worcester, MA 01655, USA
| | - Francesca Vercellone
- DIETI, Università di Napoli Federico II, Via Claudio 21, 80125 Naples, Italy
- INFN Napoli, Complesso Universitario di Monte Sant'Angelo, 80126 Naples, Italy
| | - Andrea Fontana
- Dipartimento di Fisica, Università di Napoli Federico II, and INFN Napoli, Complesso Universitario di Monte Sant'Angelo, 80126 Naples, Italy
| | - Florinda Di Pierno
- DIETI, Università di Napoli Federico II, Via Claudio 21, 80125 Naples, Italy
- INFN Napoli, Complesso Universitario di Monte Sant'Angelo, 80126 Naples, Italy
| | - Job Dekker
- Department of Systems Biology, University of Massachusetts Chan Medical School, Worcester, MA 01605, USA
- Howard Hughes Medical Institute, Chevy Chase, MD 20815, USA
| | - Mario Nicodemi
- Dipartimento di Fisica, Università di Napoli Federico II, and INFN Napoli, Complesso Universitario di Monte Sant'Angelo, 80126 Naples, Italy
| |
Collapse
|
3
|
Lao Z, Kamat KD, Jiang Z, Zhang B. OpenNucleome for high-resolution nuclear structural and dynamical modeling. eLife 2024; 13:RP93223. [PMID: 39146200 PMCID: PMC11326778 DOI: 10.7554/elife.93223] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/17/2024] Open
Abstract
The intricate structural organization of the human nucleus is fundamental to cellular function and gene regulation. Recent advancements in experimental techniques, including high-throughput sequencing and microscopy, have provided valuable insights into nuclear organization. Computational modeling has played significant roles in interpreting experimental observations by reconstructing high-resolution structural ensembles and uncovering organization principles. However, the absence of standardized modeling tools poses challenges for furthering nuclear investigations. We present OpenNucleome-an open-source software designed for conducting GPU-accelerated molecular dynamics simulations of the human nucleus. OpenNucleome offers particle-based representations of chromosomes at a resolution of 100 KB, encompassing nuclear lamina, nucleoli, and speckles. This software furnishes highly accurate structural models of nuclear architecture, affording the means for dynamic simulations of condensate formation, fusion, and exploration of non-equilibrium effects. We applied OpenNucleome to uncover the mechanisms driving the emergence of 'fixed points' within the nucleus-signifying genomic loci robustly anchored in proximity to specific nuclear bodies for functional purposes. This anchoring remains resilient even amidst significant fluctuations in chromosome radial positions and nuclear shapes within individual cells. Our findings lend support to a nuclear zoning model that elucidates genome functionality. We anticipate OpenNucleome to serve as a valuable tool for nuclear investigations, streamlining mechanistic explorations and enhancing the interpretation of experimental observations.
Collapse
Affiliation(s)
- Zhuohan Lao
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, United States
| | - Kartik D Kamat
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, United States
| | - Zhongling Jiang
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, United States
| | - Bin Zhang
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, United States
| |
Collapse
|
4
|
Conte M, Abraham A, Esposito A, Yang L, Gibcus JH, Parsi KM, Vercellone F, Fontana A, Pierno FD, Dekker J, Nicodemi M. Polymer physics models reveal structural folding features of single-molecule gene chromatin conformations. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2024:2024.07.16.603769. [PMID: 39071404 PMCID: PMC11275793 DOI: 10.1101/2024.07.16.603769] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 07/30/2024]
Abstract
Here, we employ polymer physics models of chromatin to investigate the 3D folding of a 2Mb wide genomic region encompassing the human LTN1 gene, a crucial DNA locus involved in key cellular functions. Through extensive Molecular Dynamics simulations, we reconstruct in-silico the ensemble of single-molecule LTN1 3D structures, which we benchmark against recent in-situ Hi-C 2.0 data. The model-derived single molecules are then used to predict structural folding features at the single-cell level, providing testable predictions for super-resolution microscopy experiments.
Collapse
Affiliation(s)
- Mattia Conte
- Dipartimento di Fisica, Università di Napoli Federico II, and INFN Napoli, Complesso Universitario di Monte Sant’Angelo, 80126 Naples, Italy
| | - Alex Abraham
- Dipartimento di Fisica, Università di Napoli Federico II, and INFN Napoli, Complesso Universitario di Monte Sant’Angelo, 80126 Naples, Italy
| | - Andrea Esposito
- Dipartimento di Fisica, Università di Napoli Federico II, and INFN Napoli, Complesso Universitario di Monte Sant’Angelo, 80126 Naples, Italy
| | - Liyan Yang
- Department of Systems Biology, University of Massachusetts Chan Medical School, Worcester, MA 01605, USA
| | - Johan H. Gibcus
- Department of Systems Biology, University of Massachusetts Chan Medical School, Worcester, MA 01605, USA
| | - Krishna M. Parsi
- Diabetes Center of Excellence and Program in Molecular Medicine, University of Massachusetts Chan Medical School, Worcester, MA 01655
| | - Francesca Vercellone
- DIETI, Università di Napoli Federico II, Via Claudio 21, 80125 Naples, and INFN Napoli, Complesso Universitario di Monte Sant’Angelo, 80126 Naples, Italy
| | - Andrea Fontana
- Dipartimento di Fisica, Università di Napoli Federico II, and INFN Napoli, Complesso Universitario di Monte Sant’Angelo, 80126 Naples, Italy
| | - Florinda Di Pierno
- DIETI, Università di Napoli Federico II, Via Claudio 21, 80125 Naples, and INFN Napoli, Complesso Universitario di Monte Sant’Angelo, 80126 Naples, Italy
| | - Job Dekker
- Department of Systems Biology, University of Massachusetts Chan Medical School, Worcester, MA 01605, USA
- Howard Hughes Medical Institute, Chevy Chase, MD 20815, USA
| | - Mario Nicodemi
- Dipartimento di Fisica, Università di Napoli Federico II, and INFN Napoli, Complesso Universitario di Monte Sant’Angelo, 80126 Naples, Italy
| |
Collapse
|
5
|
Sandoval-Velasco M, Dudchenko O, Rodríguez JA, Pérez Estrada C, Dehasque M, Fontsere C, Mak SST, Khan R, Contessoto VG, Oliveira Junior AB, Kalluchi A, Zubillaga Herrera BJ, Jeong J, Roy RP, Christopher I, Weisz D, Omer AD, Batra SS, Shamim MS, Durand NC, O'Connell B, Roca AL, Plikus MV, Kusliy MA, Romanenko SA, Lemskaya NA, Serdyukova NA, Modina SA, Perelman PL, Kizilova EA, Baiborodin SI, Rubtsov NB, Machol G, Rath K, Mahajan R, Kaur P, Gnirke A, Garcia-Treviño I, Coke R, Flanagan JP, Pletch K, Ruiz-Herrera A, Plotnikov V, Pavlov IS, Pavlova NI, Protopopov AV, Di Pierro M, Graphodatsky AS, Lander ES, Rowley MJ, Wolynes PG, Onuchic JN, Dalén L, Marti-Renom MA, Gilbert MTP, Aiden EL. Three-dimensional genome architecture persists in a 52,000-year-old woolly mammoth skin sample. Cell 2024; 187:3541-3562.e51. [PMID: 38996487 DOI: 10.1016/j.cell.2024.06.002] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/30/2023] [Revised: 03/07/2024] [Accepted: 06/03/2024] [Indexed: 07/14/2024]
Abstract
Analyses of ancient DNA typically involve sequencing the surviving short oligonucleotides and aligning to genome assemblies from related, modern species. Here, we report that skin from a female woolly mammoth (†Mammuthus primigenius) that died 52,000 years ago retained its ancient genome architecture. We use PaleoHi-C to map chromatin contacts and assemble its genome, yielding 28 chromosome-length scaffolds. Chromosome territories, compartments, loops, Barr bodies, and inactive X chromosome (Xi) superdomains persist. The active and inactive genome compartments in mammoth skin more closely resemble Asian elephant skin than other elephant tissues. Our analyses uncover new biology. Differences in compartmentalization reveal genes whose transcription was potentially altered in mammoths vs. elephants. Mammoth Xi has a tetradic architecture, not bipartite like human and mouse. We hypothesize that, shortly after this mammoth's death, the sample spontaneously freeze-dried in the Siberian cold, leading to a glass transition that preserved subfossils of ancient chromosomes at nanometer scale.
Collapse
Affiliation(s)
| | - Olga Dudchenko
- The Center for Genome Architecture and Department of Molecular and Human Genetics, Baylor College of Medicine, Houston, TX 77030, USA; Center for Theoretical Biological Physics, Rice University, Houston, TX 77030, USA.
| | - Juan Antonio Rodríguez
- Center for Evolutionary Hologenomics, University of Copenhagen, DK-1353 Copenhagen, Denmark; Centre Nacional d'Anàlisi Genòmica, CNAG, 08028 Barcelona, Spain
| | - Cynthia Pérez Estrada
- The Center for Genome Architecture and Department of Molecular and Human Genetics, Baylor College of Medicine, Houston, TX 77030, USA; Center for Theoretical Biological Physics, Rice University, Houston, TX 77030, USA
| | - Marianne Dehasque
- Centre for Palaeogenetics, SE-106 91 Stockholm, Sweden; Department of Bioinformatics and Genetics, Swedish Museum of Natural History, 10405 Stockholm, Sweden; Department of Zoology, Stockholm University, SE-106 91 Stockholm, Sweden
| | - Claudia Fontsere
- Center for Evolutionary Hologenomics, University of Copenhagen, DK-1353 Copenhagen, Denmark
| | - Sarah S T Mak
- Center for Evolutionary Hologenomics, University of Copenhagen, DK-1353 Copenhagen, Denmark
| | - Ruqayya Khan
- The Center for Genome Architecture and Department of Molecular and Human Genetics, Baylor College of Medicine, Houston, TX 77030, USA
| | | | | | - Achyuth Kalluchi
- Department of Genetics, Cell Biology and Anatomy, University of Nebraska Medical Center, Omaha, NE 68198, USA
| | - Bernardo J Zubillaga Herrera
- Department of Physics, Northeastern University, Boston, MA 02115, USA; Center for Theoretical Biological Physics, Northeastern University, Boston, MA 02215, USA
| | - Jiyun Jeong
- The Center for Genome Architecture and Department of Molecular and Human Genetics, Baylor College of Medicine, Houston, TX 77030, USA
| | - Renata P Roy
- The Center for Genome Architecture and Department of Molecular and Human Genetics, Baylor College of Medicine, Houston, TX 77030, USA; Center for Theoretical Biological Physics, Rice University, Houston, TX 77030, USA; Departments of Biology and Physics, Texas Southern University, Houston, TX 77004, USA
| | - Ishawnia Christopher
- The Center for Genome Architecture and Department of Molecular and Human Genetics, Baylor College of Medicine, Houston, TX 77030, USA
| | - David Weisz
- The Center for Genome Architecture and Department of Molecular and Human Genetics, Baylor College of Medicine, Houston, TX 77030, USA
| | - Arina D Omer
- The Center for Genome Architecture and Department of Molecular and Human Genetics, Baylor College of Medicine, Houston, TX 77030, USA
| | - Sanjit S Batra
- The Center for Genome Architecture and Department of Molecular and Human Genetics, Baylor College of Medicine, Houston, TX 77030, USA
| | - Muhammad S Shamim
- The Center for Genome Architecture and Department of Molecular and Human Genetics, Baylor College of Medicine, Houston, TX 77030, USA
| | - Neva C Durand
- The Center for Genome Architecture and Department of Molecular and Human Genetics, Baylor College of Medicine, Houston, TX 77030, USA; Broad Institute of MIT and Harvard, Cambridge, MA 02142, USA
| | - Brendan O'Connell
- Department of Biomolecular Engineering, University of California, Santa Cruz, Santa Cruz, CA 95064, USA; Department of Molecular and Medical Genetics, Oregon Health & Science University, Portland, OR 97239, USA
| | - Alfred L Roca
- Department of Animal Sciences and Carl R. Woese Institute for Genomic Biology, University of Illinois at Urbana-Champaign, Urbana, IL 61801, USA
| | - Maksim V Plikus
- Department of Developmental and Cell Biology, University of California, Irvine, Irvine, CA 92697, USA
| | - Mariya A Kusliy
- Institute of Molecular and Cellular Biology SB RAS, Novosibirsk 630090, Russia
| | | | - Natalya A Lemskaya
- Institute of Molecular and Cellular Biology SB RAS, Novosibirsk 630090, Russia
| | | | - Svetlana A Modina
- Institute of Molecular and Cellular Biology SB RAS, Novosibirsk 630090, Russia
| | - Polina L Perelman
- Institute of Molecular and Cellular Biology SB RAS, Novosibirsk 630090, Russia
| | - Elena A Kizilova
- Institute of Cytology and Genetics SB RAS, Novosibirsk 630090, Russia
| | | | - Nikolai B Rubtsov
- Institute of Cytology and Genetics SB RAS, Novosibirsk 630090, Russia
| | - Gur Machol
- The Center for Genome Architecture and Department of Molecular and Human Genetics, Baylor College of Medicine, Houston, TX 77030, USA
| | - Krisha Rath
- The Center for Genome Architecture and Department of Molecular and Human Genetics, Baylor College of Medicine, Houston, TX 77030, USA
| | - Ragini Mahajan
- The Center for Genome Architecture and Department of Molecular and Human Genetics, Baylor College of Medicine, Houston, TX 77030, USA; Center for Theoretical Biological Physics, Rice University, Houston, TX 77030, USA; Department of Biosciences, Rice University, Houston, TX 77005, USA
| | - Parwinder Kaur
- UWA School of Agriculture and Environment, University of Western Australia, Perth, WA 6009, Australia
| | - Andreas Gnirke
- Broad Institute of MIT and Harvard, Cambridge, MA 02142, USA
| | | | - Rob Coke
- San Antonio Zoo, San Antonio, TX 78212, USA
| | | | | | - Aurora Ruiz-Herrera
- Departament de Biologia Cel·lular, Fisiologia i Immunologia and Genome Integrity and Instability Group, Institut de Biotecnologia i Biomedicina, Universitat Autònoma de Barcelona, 08193 Cerdanyola del Vallès, Spain
| | | | | | - Naryya I Pavlova
- Institute of Biological Problems of Cryolitezone SB RAS, Yakutsk 677000, Russia
| | - Albert V Protopopov
- Academy of Sciences of Sakha Republic, Yakutsk 677000, Russia; North-Eastern Federal University, Yakutsk 677027, Russia
| | - Michele Di Pierro
- Department of Physics, Northeastern University, Boston, MA 02115, USA; Center for Theoretical Biological Physics, Northeastern University, Boston, MA 02215, USA
| | | | - Eric S Lander
- Broad Institute of MIT and Harvard, Cambridge, MA 02142, USA; Department of Systems Biology, Harvard Medical School, Boston, MA 02115, USA; Department of Biology, Massachusetts Institute of Technology, Cambridge, MA 02139, USA
| | - M Jordan Rowley
- Department of Genetics, Cell Biology and Anatomy, University of Nebraska Medical Center, Omaha, NE 68198, USA
| | - Peter G Wolynes
- Center for Theoretical Biological Physics, Rice University, Houston, TX 77030, USA; Department of Biosciences, Rice University, Houston, TX 77005, USA; Departments of Physics, Astronomy, & Chemistry, Rice University, Houston, TX 77005, USA
| | - José N Onuchic
- Center for Theoretical Biological Physics, Rice University, Houston, TX 77030, USA; Department of Biosciences, Rice University, Houston, TX 77005, USA; Departments of Physics, Astronomy, & Chemistry, Rice University, Houston, TX 77005, USA
| | - Love Dalén
- Centre for Palaeogenetics, SE-106 91 Stockholm, Sweden; Department of Bioinformatics and Genetics, Swedish Museum of Natural History, 10405 Stockholm, Sweden; Department of Zoology, Stockholm University, SE-106 91 Stockholm, Sweden
| | - Marc A Marti-Renom
- Centre Nacional d'Anàlisi Genòmica, CNAG, 08028 Barcelona, Spain; Centre for Genomic Regulation, The Barcelona Institute for Science and Technology, 08003 Barcelona, Spain; ICREA, 08010 Barcelona, Spain; Universitat Pompeu Fabra, 08002 Barcelona, Spain.
| | - M Thomas P Gilbert
- Center for Evolutionary Hologenomics, University of Copenhagen, DK-1353 Copenhagen, Denmark; University Museum NTNU, 7012 Trondheim, Norway.
| | - Erez Lieberman Aiden
- The Center for Genome Architecture and Department of Molecular and Human Genetics, Baylor College of Medicine, Houston, TX 77030, USA; Center for Theoretical Biological Physics, Rice University, Houston, TX 77030, USA; Broad Institute of MIT and Harvard, Cambridge, MA 02142, USA.
| |
Collapse
|
6
|
Cao Z, Wolynes PG. Motorized chain models of the ideal chromosome. Proc Natl Acad Sci U S A 2024; 121:e2407077121. [PMID: 38954553 PMCID: PMC11252987 DOI: 10.1073/pnas.2407077121] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/08/2024] [Accepted: 06/06/2024] [Indexed: 07/04/2024] Open
Abstract
An array of motor proteins consumes chemical energy in setting up the architectures of chromosomes. Here, we explore how the structure of ideal polymer chains is influenced by two classes of motors. The first class which we call "swimming motors" acts to propel the chromatin fiber through three-dimensional space. They represent a caricature of motors such as RNA polymerases. Previously, they have often been described by adding a persistent flow onto Brownian diffusion of the chain. The second class of motors, which we call "grappling motors" caricatures the loop extrusion processes in which segments of chromatin fibers some distance apart are brought together. We analyze these models using a self-consistent variational phonon approximation to a many-body Master equation incorporating motor activities. We show that whether the swimming motors lead to contraction or expansion depends on the susceptibility of the motors, that is, how their activity depends on the forces they must exert. Grappling motors in contrast to swimming motors lead to long-ranged correlations that resemble those first suggested for fractal globules and that are consistent with the effective interactions inferred by energy landscape analyses of Hi-C data on the interphase chromosome.
Collapse
Affiliation(s)
- Zhiyu Cao
- Center for Theoretical Biological Physics, Rice University, Houston, TX77005
- Department of Chemical Physics, University of Science and Technology of China, Hefei, Anhui230026, China
| | - Peter G. Wolynes
- Center for Theoretical Biological Physics, Rice University, Houston, TX77005
- Department of Chemistry, Rice University, Houston, TX77005
- Department of Physics, Rice University, Houston, TX77005
| |
Collapse
|
7
|
Liu S, Athreya A, Lao Z, Zhang B. From Nucleosomes to Compartments: Physicochemical Interactions Underlying Chromatin Organization. Annu Rev Biophys 2024; 53:221-245. [PMID: 38346246 PMCID: PMC11369498 DOI: 10.1146/annurev-biophys-030822-032650] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/18/2024]
Abstract
Chromatin organization plays a critical role in cellular function by regulating access to genetic information. However, understanding chromatin folding is challenging due to its complex, multiscale nature. Significant progress has been made in studying in vitro systems, uncovering the structure of individual nucleosomes and their arrays, and elucidating the role of physicochemical forces in stabilizing these structures. Additionally, remarkable advancements have been achieved in characterizing chromatin organization in vivo, particularly at the whole-chromosome level, revealing important features such as chromatin loops, topologically associating domains, and nuclear compartments. However, bridging the gap between in vitro and in vivo studies remains challenging. The resemblance between in vitro and in vivo chromatin conformations and the relevance of internucleosomal interactions for chromatin folding in vivo are subjects of debate. This article reviews experimental and computational studies conducted at various length scales, highlighting the significance of intrinsic interactions between nucleosomes and their roles in chromatin folding in vivo.
Collapse
Affiliation(s)
- Shuming Liu
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts, USA;
| | - Advait Athreya
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts, USA;
| | - Zhuohan Lao
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts, USA;
| | - Bin Zhang
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts, USA;
| |
Collapse
|
8
|
Kadlof M, Banecki K, Chiliński M, Plewczynski D. Chromatin image-driven modelling. Methods 2024; 226:54-60. [PMID: 38636797 DOI: 10.1016/j.ymeth.2024.04.006] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/07/2023] [Revised: 03/13/2024] [Accepted: 04/05/2024] [Indexed: 04/20/2024] Open
Abstract
The challenge of modelling the spatial conformation of chromatin remains an open problem. While multiple data-driven approaches have been proposed, each has limitations. This work introduces two image-driven modelling methods based on the Molecular Dynamics Flexible Fitting (MDFF) approach: the force method and the correlational method. Both methods have already been used successfully in protein modelling. We propose a novel way to employ them for building chromatin models directly from 3D images. This approach is termed image-driven modelling. Additionally, we introduce the initial structure generator, a tool designed to generate optimal starting structures for the proposed algorithms. The methods are versatile and can be applied to various data types, with minor modifications to accommodate new generation imaging techniques.
Collapse
Affiliation(s)
- Michał Kadlof
- Faculty of Mathematics and Information Science, Warsaw University of Technology, 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
| | - Mateusz Chiliński
- Faculty of Mathematics and Information Science, Warsaw University of Technology, Warsaw, Poland; Centre of New Technologies, University of Warsaw, Warsaw, Poland; Section for Computational and RNA Biology, Department of Biology, University of Copenhagen, Copenhagen, Denmark
| | - Dariusz Plewczynski
- Faculty of Mathematics and Information Science, Warsaw University of Technology, Warsaw, Poland; Centre of New Technologies, University of Warsaw, Warsaw, Poland
| |
Collapse
|
9
|
Wei J, Xue Y, Liu Y, Tian H, Shao Y, Gao YQ. Steric repulsion introduced by loop constraints modulates the microphase separation of chromatins. J Chem Phys 2024; 160:054904. [PMID: 38341710 DOI: 10.1063/5.0189692] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/30/2023] [Accepted: 01/15/2024] [Indexed: 02/13/2024] Open
Abstract
Within the confines of a densely populated cell nucleus, chromatin undergoes intricate folding, forming loops, domains, and compartments under the governance of topological constraints and phase separation. This coordinated process inevitably introduces interference between different folding strategies. In this study, we model interphase chromatins as block copolymers with hetero-hierarchical loops within a confined system. Employing dissipative particle dynamics simulations and scaling analysis, we aim to explain how the structure and distribution of loop domains modulate the microphase separation of chromatins. Our results highlight the correlation between the microphase separation of the copolymer and the length, heterogeneity, and hierarchically nested levels of the loop domains. This correlation arises from steric repulsion intrinsic to loop domains. The steric repulsion induces variations in chain stiffness (including local orientation correlations and the persistence length), thereby influencing the degree of phase separation. Through simulations of block copolymers with distinct groups of hetero-hierarchical loop anchors, we successfully reproduce changes in phase separation across diverse cell lines, under fixed interaction parameters. These findings, in qualitative alignment with Hi-C data, suggest that the variations of loop constraints alone possess the capacity to regulate higher-order structures and the gene expressions of interphase chromatins.
Collapse
Affiliation(s)
| | - Yue Xue
- Beijing National Laboratory for Molecular Sciences, College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, China
- Biomedical Pioneering Innovation Center (BIOPIC), Peking University, Beijing 100871, China
| | - Yawei Liu
- Beijing Key Laboratory of Ionic Liquids Clean Process, CAS Key Laboratory of Green Process and Engineering, State Key Laboratory of Multiphase Complex Systems, Institute of Process Engineering, Chinese Academy of Sciences, Beijing 100190, China
| | - Hao Tian
- Beijing National Laboratory for Molecular Sciences, College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, China
- Biomedical Pioneering Innovation Center (BIOPIC), Peking University, Beijing 100871, China
| | - Yingfeng Shao
- State Key Laboratory of Nonlinear Mechanics, Institute of Mechanics, Chinese Academy of Sciences, Beijing 100190, China
| | - Yi Qin Gao
- Changping Laboratory, Beijing 102206, China
- Beijing National Laboratory for Molecular Sciences, College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, China
- Biomedical Pioneering Innovation Center (BIOPIC), Peking University, Beijing 100871, China
- Shenzhen Bay Laboratory, 5F, No. 9 Duxue Rd., Nanshan District, Shenzhen 518055, Guangdong, China
| |
Collapse
|
10
|
Zhang Y, Boninsegna L, Yang M, Misteli T, Alber F, Ma J. Computational methods for analysing multiscale 3D genome organization. Nat Rev Genet 2024; 25:123-141. [PMID: 37673975 PMCID: PMC11127719 DOI: 10.1038/s41576-023-00638-1] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 07/12/2023] [Indexed: 09/08/2023]
Abstract
Recent progress in whole-genome mapping and imaging technologies has enabled the characterization of the spatial organization and folding of the genome in the nucleus. In parallel, advanced computational methods have been developed to leverage these mapping data to reveal multiscale three-dimensional (3D) genome features and to provide a more complete view of genome structure and its connections to genome functions such as transcription. Here, we discuss how recently developed computational tools, including machine-learning-based methods and integrative structure-modelling frameworks, have led to a systematic, multiscale delineation of the connections among different scales of 3D genome organization, genomic and epigenomic features, functional nuclear components and genome function. However, approaches that more comprehensively integrate a wide variety of genomic and imaging datasets are still needed to uncover the functional role of 3D genome structure in defining cellular phenotypes in health and disease.
Collapse
Affiliation(s)
- Yang Zhang
- Computational Biology Department, School of Computer Science, Carnegie Mellon University, Pittsburgh, PA, USA
| | - Lorenzo Boninsegna
- Department of Microbiology, Immunology and Molecular Genetics and Institute for Quantitative and Computational Biosciences, University of California Los Angeles, Los Angeles, CA, USA
| | - Muyu Yang
- Computational Biology Department, School of Computer Science, Carnegie Mellon University, Pittsburgh, PA, USA
| | - Tom Misteli
- Center for Cancer Research, National Cancer Institute, Bethesda, MD, USA.
| | - Frank Alber
- Department of Microbiology, Immunology and Molecular Genetics and Institute for Quantitative and Computational Biosciences, University of California Los Angeles, Los Angeles, CA, USA.
| | - Jian Ma
- Computational Biology Department, School of Computer Science, Carnegie Mellon University, Pittsburgh, PA, USA.
| |
Collapse
|
11
|
Feng C, Wang J, Chu X. Large-scale data-driven and physics-based models offer insights into the relationships among the structures, dynamics, and functions of chromosomes. J Mol Cell Biol 2023; 15:mjad042. [PMID: 37365687 PMCID: PMC10782906 DOI: 10.1093/jmcb/mjad042] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/06/2023] [Revised: 03/22/2023] [Accepted: 06/25/2023] [Indexed: 06/28/2023] Open
Abstract
The organized three-dimensional chromosome architecture in the cell nucleus provides scaffolding for precise regulation of gene expression. When the cell changes its identity in the cell-fate decision-making process, extensive rearrangements of chromosome structures occur accompanied by large-scale adaptations of gene expression, underscoring the importance of chromosome dynamics in shaping genome function. Over the last two decades, rapid development of experimental methods has provided unprecedented data to characterize the hierarchical structures and dynamic properties of chromosomes. In parallel, these enormous data offer valuable opportunities for developing quantitative computational models. Here, we review a variety of large-scale polymer models developed to investigate the structures and dynamics of chromosomes. Different from the underlying modeling strategies, these approaches can be classified into data-driven ('top-down') and physics-based ('bottom-up') categories. We discuss their contributions to offering valuable insights into the relationships among the structures, dynamics, and functions of chromosomes and propose the perspective of developing data integration approaches from different experimental technologies and multidisciplinary theoretical/simulation methods combined with different modeling strategies.
Collapse
Affiliation(s)
- Cibo Feng
- Advanced Materials Thrust, Function Hub, The Hong Kong University of Science and Technology (Guangzhou), Guangzhou 511400, China
- Green e Materials Laboratory, The Hong Kong University of Science and Technology (Guangzhou), Guangzhou 511400, China
- College of Physics, Jilin University, Changchun 130012, China
| | - Jin Wang
- Department of Chemistry and Physics, The State University of New York at Stony Brook, Stony Brook, NY 11794, USA
| | - Xiakun Chu
- Advanced Materials Thrust, Function Hub, The Hong Kong University of Science and Technology (Guangzhou), Guangzhou 511400, China
- Green e Materials Laboratory, The Hong Kong University of Science and Technology (Guangzhou), Guangzhou 511400, China
- Division of Life Science, The Hong Kong University of Science and Technology, Hong Kong SAR 999077, China
- Guangzhou Municipal Key Laboratory of Materials Informatics, The Hong Kong University of Science and Technology (Guangzhou), Guangzhou 511400, China
| |
Collapse
|
12
|
Chu X, Wang J. Quantifying the large-scale chromosome structural dynamics during the mitosis-to-G1 phase transition of cell cycle. Open Biol 2023; 13:230175. [PMID: 37907089 PMCID: PMC10618054 DOI: 10.1098/rsob.230175] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/09/2023] [Accepted: 09/27/2023] [Indexed: 11/02/2023] Open
Abstract
Cell cycle is known to be regulated by the underlying gene network. Chromosomes, which serve as the scaffold for gene expressions, undergo significant structural reorganizations during mitosis. Understanding the mechanism of the cell cycle from the chromosome structural perspective remains a grand challenge. In this study, we applied an integrated theoretical approach to investigate large-scale chromosome structural dynamics during the mitosis-to-G1 phase transition. We observed that the chromosome structural expansion and adaptation of the structural asphericity do not occur synchronously and attributed this behaviour to the unique unloading sequence of the two types of condensins. Furthermore, we observed that the coherent motions between the chromosomal loci are primarily enhanced within the topologically associating domains (TADs) as cells progress to the G1 phase, suggesting that TADs can be considered as both structural and dynamical units for organizing the three-dimensional chromosome. Our analysis also reveals that the quantified pathways of chromosome structural reorganization during the mitosis-to-G1 phase transition exhibit high stochasticity at the single-cell level and show nonlinear behaviours in changing TADs and contacts formed at the long-range regions. Our findings offer valuable insights into large-scale chromosome structural dynamics after mitosis.
Collapse
Affiliation(s)
- Xiakun Chu
- Department of Chemistry, State University of New York at Stony Brook, Stony Brook, NY 11794, USA
| | - Jin Wang
- Department of Chemistry, State University of New York at Stony Brook, Stony Brook, NY 11794, USA
- Department of Physics and Astronomy, State University of New York at Stony Brook, Stony Brook, NY 11794, USA
| |
Collapse
|
13
|
Junier I, Ghobadpour E, Espeli O, Everaers R. DNA supercoiling in bacteria: state of play and challenges from a viewpoint of physics based modeling. Front Microbiol 2023; 14:1192831. [PMID: 37965550 PMCID: PMC10642903 DOI: 10.3389/fmicb.2023.1192831] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/23/2023] [Accepted: 09/25/2023] [Indexed: 11/16/2023] Open
Abstract
DNA supercoiling is central to many fundamental processes of living organisms. Its average level along the chromosome and over time reflects the dynamic equilibrium of opposite activities of topoisomerases, which are required to relax mechanical stresses that are inevitably produced during DNA replication and gene transcription. Supercoiling affects all scales of the spatio-temporal organization of bacterial DNA, from the base pair to the large scale chromosome conformation. Highlighted in vitro and in vivo in the 1960s and 1970s, respectively, the first physical models were proposed concomitantly in order to predict the deformation properties of the double helix. About fifteen years later, polymer physics models demonstrated on larger scales the plectonemic nature and the tree-like organization of supercoiled DNA. Since then, many works have tried to establish a better understanding of the multiple structuring and physiological properties of bacterial DNA in thermodynamic equilibrium and far from equilibrium. The purpose of this essay is to address upcoming challenges by thoroughly exploring the relevance, predictive capacity, and limitations of current physical models, with a specific focus on structural properties beyond the scale of the double helix. We discuss more particularly the problem of DNA conformations, the interplay between DNA supercoiling with gene transcription and DNA replication, its role on nucleoid formation and, finally, the problem of scaling up models. Our primary objective is to foster increased collaboration between physicists and biologists. To achieve this, we have reduced the respective jargon to a minimum and we provide some explanatory background material for the two communities.
Collapse
Affiliation(s)
- Ivan Junier
- CNRS, UMR 5525, VetAgro Sup, Grenoble INP, TIMC, Université Grenoble Alpes, Grenoble, France
| | - Elham Ghobadpour
- CNRS, UMR 5525, VetAgro Sup, Grenoble INP, TIMC, Université Grenoble Alpes, Grenoble, France
- École Normale Supérieure (ENS) de Lyon, CNRS, Laboratoire de Physique and Centre Blaise Pascal de l'ENS de Lyon, Lyon, France
| | - Olivier Espeli
- Center for Interdisciplinary Research in Biology (CIRB), Collège de France, CNRS, INSERM, Université PSL, Paris, France
| | - Ralf Everaers
- École Normale Supérieure (ENS) de Lyon, CNRS, Laboratoire de Physique and Centre Blaise Pascal de l'ENS de Lyon, Lyon, France
| |
Collapse
|
14
|
Jain N, Thakur S. Structure and dynamics of chemically active ring polymers: swelling to collapse. SOFT MATTER 2023; 19:7358-7369. [PMID: 37740385 DOI: 10.1039/d3sm00839h] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 09/24/2023]
Abstract
The ring structures are common in many synthetic or natural systems and experience both local and long-range forces by chemical sensing. This work is an effort to investigate the structural and dynamical properties of a chemically active ring in an explicit solvent bath utilizing hybrid molecular dynamics (MD) and multiparticle collision dynamics (MPCD) simulation techniques. We show that by tuning the chemical properties of the ring, it can be converted from a chemo-attractant to a chemo-repellent, thereby changing the steady state to be either collapsed or swelled as compared to its passive limit. We quantify these observations by comparing the scaling laws, local structures and the dynamics of active and passive rings. Furthermore, we show the impact of varying numbers of active sites by calculating the contact probability of the collapse state that highlights diverse structures. We also analyze the dynamics of the ring by finding the relaxation time and the mean square displacement of the centre of mass. A faster relaxation with enhanced diffusion is observed for the active rings.
Collapse
Affiliation(s)
- Namita Jain
- Department of Physics, Indian Institute of Science Education and Research Bhopal, Bhopal 462066, India.
| | - Snigdha Thakur
- Department of Physics, Indian Institute of Science Education and Research Bhopal, Bhopal 462066, India.
| |
Collapse
|
15
|
Schuette G, Ding X, Zhang B. Efficient Hi-C inversion facilitates chromatin folding mechanism discovery and structure prediction. Biophys J 2023; 122:3425-3438. [PMID: 37496267 PMCID: PMC10502442 DOI: 10.1016/j.bpj.2023.07.017] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/17/2023] [Revised: 07/10/2023] [Accepted: 07/24/2023] [Indexed: 07/28/2023] Open
Abstract
Genome-wide chromosome conformation capture (Hi-C) experiments have revealed many structural features of chromatin across multiple length scales. Further understanding genome organization requires relating these discoveries to the mechanisms that establish chromatin structures and reconstructing these structures in three dimensions, but both objectives are difficult to achieve with existing algorithms that are often computationally expensive. To alleviate this challenge, we present an algorithm that efficiently converts Hi-C data into contact energies, which measure the interaction strength between genomic loci brought into proximity. Contact energies are local quantities unaffected by the topological constraints that correlate Hi-C contact probabilities. Thus, extracting contact energies from Hi-C contact probabilities distills the biologically unique information contained in the data. We show that contact energies reveal the location of chromatin loop anchors, support a phase separation mechanism for genome compartmentalization, and parameterize polymer simulations that predict three-dimensional chromatin structures. Therefore, we anticipate that contact energy extraction will unleash the full potential of Hi-C data and that our inversion algorithm will facilitate the widespread adoption of contact energy analysis.
Collapse
Affiliation(s)
- Greg Schuette
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts
| | - Xinqiang Ding
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts
| | - Bin Zhang
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts.
| |
Collapse
|
16
|
Yildirim A, Hua N, Boninsegna L, Zhan Y, Polles G, Gong K, Hao S, Li W, Zhou XJ, Alber F. Evaluating the role of the nuclear microenvironment in gene function by population-based modeling. Nat Struct Mol Biol 2023; 30:1193-1206. [PMID: 37580627 PMCID: PMC10442234 DOI: 10.1038/s41594-023-01036-1] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/03/2021] [Accepted: 06/16/2023] [Indexed: 08/16/2023]
Abstract
The nuclear folding of chromosomes relative to nuclear bodies is an integral part of gene function. Here, we demonstrate that population-based modeling-from ensemble Hi-C data-provides a detailed description of the nuclear microenvironment of genes and its role in gene function. We define the microenvironment by the subnuclear positions of genomic regions with respect to nuclear bodies, local chromatin compaction, and preferences in chromatin compartmentalization. These structural descriptors are determined in single-cell models, thereby revealing the structural variability between cells. We demonstrate that the microenvironment of a genomic region is linked to its functional potential in gene transcription, replication, and chromatin compartmentalization. Some chromatin regions feature a strong preference for a single microenvironment, due to association with specific nuclear bodies in most cells. Other chromatin shows high structural variability, which is a strong indicator of functional heterogeneity. Moreover, we identify specialized nuclear microenvironments, which distinguish chromatin in different functional states and reveal a key role of nuclear speckles in chromosome organization. We demonstrate that our method produces highly predictive three-dimensional genome structures, which accurately reproduce data from a variety of orthogonal experiments, thus considerably expanding the range of Hi-C data analysis.
Collapse
Affiliation(s)
- Asli Yildirim
- Institute for Quantitative and Computational Biosciences, University of California Los Angeles, Los Angeles, CA, USA
- Department of Microbiology, Immunology, and Molecular Genetics, University of California Los Angeles, Los Angeles, CA, USA
| | - Nan Hua
- Institute for Quantitative and Computational Biosciences, University of California Los Angeles, Los Angeles, CA, USA
- Department of Microbiology, Immunology, and Molecular Genetics, University of California Los Angeles, Los Angeles, CA, USA
| | - Lorenzo Boninsegna
- Institute for Quantitative and Computational Biosciences, University of California Los Angeles, Los Angeles, CA, USA
- Department of Microbiology, Immunology, and Molecular Genetics, University of California Los Angeles, Los Angeles, CA, USA
| | - Yuxiang Zhan
- Institute for Quantitative and Computational Biosciences, University of California Los Angeles, Los Angeles, CA, USA
- Department of Microbiology, Immunology, and Molecular Genetics, University of California Los Angeles, Los Angeles, CA, USA
- Department of Quantitative and Computational Biology, University of Southern California, Los Angeles, CA, USA
| | - Guido Polles
- Institute for Quantitative and Computational Biosciences, University of California Los Angeles, Los Angeles, CA, USA
- Department of Microbiology, Immunology, and Molecular Genetics, University of California Los Angeles, Los Angeles, CA, USA
| | - Ke Gong
- Institute for Quantitative and Computational Biosciences, University of California Los Angeles, Los Angeles, CA, USA
- Department of Microbiology, Immunology, and Molecular Genetics, University of California Los Angeles, Los Angeles, CA, USA
| | - Shengli Hao
- Institute for Quantitative and Computational Biosciences, University of California Los Angeles, Los Angeles, CA, USA
- Department of Microbiology, Immunology, and Molecular Genetics, University of California Los Angeles, Los Angeles, CA, USA
| | - Wenyuan Li
- Department of Pathology, David Geffen School of Medicine, University of California Los Angeles, Los Angeles, CA, USA
| | - Xianghong Jasmine Zhou
- Department of Pathology, David Geffen School of Medicine, University of California Los Angeles, Los Angeles, CA, USA
| | - Frank Alber
- Institute for Quantitative and Computational Biosciences, University of California Los Angeles, Los Angeles, CA, USA.
- Department of Microbiology, Immunology, and Molecular Genetics, University of California Los Angeles, Los Angeles, CA, USA.
- Department of Quantitative and Computational Biology, University of Southern California, Los Angeles, CA, USA.
| |
Collapse
|
17
|
Schuette G, Ding X, Zhang B. Efficient Hi-C inversion facilitates chromatin folding mechanism discovery and structure prediction. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2023:2023.03.17.533194. [PMID: 36993500 PMCID: PMC10055272 DOI: 10.1101/2023.03.17.533194] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 03/31/2023]
Abstract
Genome-wide chromosome conformation capture (Hi-C) experiments have revealed many structural features of chromatin across multiple length scales. Further understanding genome organization requires relating these discoveries to the mechanisms that establish chromatin structures and reconstructing these structures in three dimensions, but both objectives are difficult to achieve with existing algorithms that are often computationally expensive. To alleviate this challenge, we present an algorithm that efficiently converts Hi-C data into contact energies, which measure the interaction strength between genomic loci brought into proximity. Contact energies are local quantities unaffected by the topological constraints that correlate Hi-C contact probabilities. Thus, extracting contact energies from Hi-C contact probabilities distills the biologically unique information contained in the data. We show that contact energies reveal the location of chromatin loop anchors, support a phase separation mechanism for genome compartmentalization, and parameterize polymer simulations that predict three-dimensional chromatin structures. Therefore, we anticipate that contact energy extraction will unleash the full potential of Hi-C data and that our inversion algorithm will facilitate the widespread adoption of contact energy analysis.
Collapse
Affiliation(s)
- Greg Schuette
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, MA 02139, USA
| | - Xinqiang Ding
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, MA 02139, USA
| | - Bin Zhang
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, MA 02139, USA
| |
Collapse
|
18
|
Chan B, Rubinstein M. Theory of chromatin organization maintained by active loop extrusion. Proc Natl Acad Sci U S A 2023; 120:e2222078120. [PMID: 37253009 PMCID: PMC10266055 DOI: 10.1073/pnas.2222078120] [Citation(s) in RCA: 6] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/02/2023] [Accepted: 04/13/2023] [Indexed: 06/01/2023] Open
Abstract
The active loop extrusion hypothesis proposes that chromatin threads through the cohesin protein complex into progressively larger loops until reaching specific boundary elements. We build upon this hypothesis and develop an analytical theory for active loop extrusion which predicts that loop formation probability is a nonmonotonic function of loop length and describes chromatin contact probabilities. We validate our model with Monte Carlo and hybrid Molecular Dynamics-Monte Carlo simulations and demonstrate that our theory recapitulates experimental chromatin conformation capture data. Our results support active loop extrusion as a mechanism for chromatin organization and provide an analytical description of chromatin organization that may be used to specifically modify chromatin contact probabilities.
Collapse
Affiliation(s)
- Brian Chan
- Department of Biomedical Engineering, Duke University, Durham, NC27708
| | - Michael Rubinstein
- Department of Biomedical Engineering, Duke University, Durham, NC27708
- Department of Mechanical Engineering and Materials Science, Duke University, Durham, NC27708
- Department of Chemistry, Duke University, Durham, NC27708
- Department of Physics, Duke University, Durham, NC27708
- Institute for Chemical Reaction Design and Discovery (World Premier International Research Center Initiative-ICReDD), Hokkaido University, Sapporo001-0021, Japan
| |
Collapse
|
19
|
Tolokh IS, Kinney NA, Sharakhov IV, Onufriev AV. Strong interactions between highly dynamic lamina-associated domains and the nuclear envelope stabilize the 3D architecture of Drosophila interphase chromatin. Epigenetics Chromatin 2023; 16:21. [PMID: 37254161 PMCID: PMC10228000 DOI: 10.1186/s13072-023-00492-9] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/19/2023] [Accepted: 05/04/2023] [Indexed: 06/01/2023] Open
Abstract
BACKGROUND Interactions among topologically associating domains (TADs), and between the nuclear envelope (NE) and lamina-associated domains (LADs) are expected to shape various aspects of three-dimensional (3D) chromatin structure and dynamics; however, relevant genome-wide experiments that may provide statistically significant conclusions remain difficult. RESULTS We have developed a coarse-grained dynamical model of D. melanogaster nuclei at TAD resolution that explicitly accounts for four distinct epigenetic classes of TADs and LAD-NE interactions. The model is parameterized to reproduce the experimental Hi-C map of the wild type (WT) nuclei; it describes time evolution of the chromatin over the G1 phase of the interphase. The simulations include an ensemble of nuclei, corresponding to the experimentally observed set of several possible mutual arrangements of chromosomal arms. The model is validated against multiple structural features of chromatin from several different experiments not used in model development. Predicted positioning of all LADs at the NE is highly dynamic-the same LAD can attach, detach and move far away from the NE multiple times during interphase. The probabilities of LADs to be in contact with the NE vary by an order of magnitude, despite all having the same affinity to the NE in the model. These probabilities are mostly determined by a highly variable local linear density of LADs along the genome, which also has the same strong effect on the predicted positioning of individual TADs -- higher probability of a TAD to be near NE is largely determined by a higher linear density of LADs surrounding this TAD. The distribution of LADs along the chromosome chains plays a notable role in maintaining a non-random average global structure of chromatin. Relatively high affinity of LADs to the NE in the WT nuclei substantially reduces sensitivity of the global radial chromatin distribution to variations in the strength of TAD-TAD interactions compared to the lamin depleted nuclei, where a small (0.5 kT) increase of cross-type TAD-TAD interactions doubles the chromatin density in the central nucleus region. CONCLUSIONS A dynamical model of the entire fruit fly genome makes multiple genome-wide predictions of biological interest. The distribution of LADs along the chromatin chains affects their probabilities to be in contact with the NE and radial positioning of highly mobile TADs, playing a notable role in creating a non-random average global structure of the chromatin. We conjecture that an important role of attractive LAD-NE interactions is to stabilize global chromatin structure against inevitable cell-to-cell variations in TAD-TAD interactions.
Collapse
Affiliation(s)
- Igor S. Tolokh
- Department of Computer Science, Virginia Tech, Blacksburg, VA 24061 USA
| | - Nicholas Allen Kinney
- Department of Computer Science, Virginia Tech, Blacksburg, VA 24061 USA
- Department of Entomology, Virginia Tech, Blacksburg, VA 24061 USA
- Edward Via College of Osteopathic Medicine, 2265 Kraft Drive, Blacksburg, VA 24060 USA
| | | | - Alexey V. Onufriev
- Department of Computer Science, Virginia Tech, Blacksburg, VA 24061 USA
- Department of Physics, Virginia Tech, Blacksburg, VA 24061 USA
- Center for Soft Matter and Biological Physics, Virginia Tech, Blacksburg, VA 24061 USA
| |
Collapse
|
20
|
Bajpai G, Safran S. Mesoscale, long-time mixing of chromosomes and its connection to polymer dynamics. PLoS Comput Biol 2023; 19:e1011142. [PMID: 37228178 DOI: 10.1371/journal.pcbi.1011142] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/06/2022] [Accepted: 05/01/2023] [Indexed: 05/27/2023] Open
Abstract
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.
Collapse
Affiliation(s)
- Gaurav Bajpai
- Department of Chemical and Biological Physics, Weizmann Institute of Science, Rehovot, Israel
| | - Samuel Safran
- Department of Chemical and Biological Physics, Weizmann Institute of Science, Rehovot, Israel
| |
Collapse
|
21
|
Goychuk A, Kannan D, Chakraborty AK, Kardar M. Polymer folding through active processes recreates features of genome organization. Proc Natl Acad Sci U S A 2023; 120:e2221726120. [PMID: 37155885 PMCID: PMC10194017 DOI: 10.1073/pnas.2221726120] [Citation(s) in RCA: 17] [Impact Index Per Article: 17.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/25/2022] [Accepted: 04/02/2023] [Indexed: 05/10/2023] Open
Abstract
From proteins to chromosomes, polymers fold into specific conformations that control their biological function. Polymer folding has long been studied with equilibrium thermodynamics, yet intracellular organization and regulation involve energy-consuming, active processes. Signatures of activity have been measured in the context of chromatin motion, which shows spatial correlations and enhanced subdiffusion only in the presence of adenosine triphosphate. Moreover, chromatin motion varies with genomic coordinate, pointing toward a heterogeneous pattern of active processes along the sequence. How do such patterns of activity affect the conformation of a polymer such as chromatin? We address this question by combining analytical theory and simulations to study a polymer subjected to sequence-dependent correlated active forces. Our analysis shows that a local increase in activity (larger active forces) can cause the polymer backbone to bend and expand, while less active segments straighten out and condense. Our simulations further predict that modest activity differences can drive compartmentalization of the polymer consistent with the patterns observed in chromosome conformation capture experiments. Moreover, segments of the polymer that show correlated active (sub)diffusion attract each other through effective long-ranged harmonic interactions, whereas anticorrelations lead to effective repulsions. Thus, our theory offers nonequilibrium mechanisms for forming genomic compartments, which cannot be distinguished from affinity-based folding using structural data alone. As a first step toward exploring whether active mechanisms contribute to shaping genome conformations, we discuss a data-driven approach.
Collapse
Affiliation(s)
- Andriy Goychuk
- Institute for Medical Engineering and Science, Massachusetts Institute of Technology, Cambridge, MA02139
| | - Deepti Kannan
- Department of Physics, Massachusetts Institute of Technology, Cambridge, MA02139
| | - Arup K. Chakraborty
- Institute for Medical Engineering and Science, Massachusetts Institute of Technology, Cambridge, MA02139
- Department of Physics, Massachusetts Institute of Technology, Cambridge, MA02139
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, MA02139
- Ragon Institute of Massachusetts General Hospital, Massachusetts Institute of Technology and Harvard University, Cambridge, MA02139
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, MA02139
| | - Mehran Kardar
- Department of Physics, Massachusetts Institute of Technology, Cambridge, MA02139
| |
Collapse
|
22
|
Ruben BS, Brahmachari S, Contessoto VG, Cheng RR, Oliveira Junior AB, Di Pierro M, Onuchic JN. Structural reorganization and relaxation dynamics of axially stressed chromosomes. Biophys J 2023; 122:1633-1645. [PMID: 36960531 PMCID: PMC10183323 DOI: 10.1016/j.bpj.2023.03.029] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/13/2022] [Revised: 02/06/2023] [Accepted: 03/17/2023] [Indexed: 03/25/2023] Open
Abstract
Chromosomes endure mechanical stresses throughout the cell cycle; for example, resulting from the pulling of chromosomes by spindle fibers during mitosis or deformation of the nucleus during cell migration. The response to physical stress is closely related to chromosome structure and function. Micromechanical studies of mitotic chromosomes have revealed them to be remarkably extensible objects and informed early models of mitotic chromosome organization. We use a data-driven, coarse-grained polymer modeling approach to explore the relationship between the spatial organization of individual chromosomes and their emergent mechanical properties. In particular, we investigate the mechanical properties of our model chromosomes by axially stretching them. Simulated stretching led to a linear force-extension curve for small strain, with mitotic chromosomes behaving about 10-fold stiffer than interphase chromosomes. Studying their relaxation dynamics, we found that chromosomes are viscoelastic solids with a highly liquid-like, viscous behavior in interphase that becomes solid-like in mitosis. This emergent mechanical stiffness originates from lengthwise compaction, an effective potential capturing the activity of loop-extruding SMC complexes. Chromosomes denature under large strains via unraveling, which is characterized by opening of large-scale folding patterns. By quantifying the effect of mechanical perturbations on the chromosome's structural features, our model provides a nuanced understanding of in vivo mechanics of chromosomes.
Collapse
Affiliation(s)
- Benjamin S Ruben
- Center for Theoretical Biological Physics, Rice University, Houston, Texas; Biophysics PhD Program, Harvard University, Cambridge, Massachusetts.
| | | | | | - Ryan R Cheng
- Center for Theoretical Biological Physics, Rice University, Houston, Texas; Department of Chemistry, University of Kentucky, Lexington, Kentucky
| | | | - Michele Di Pierro
- Department of Physics, Northeastern University, Boston, Massachusetts; Center for Theoretical Biological Physics, Northeastern University, Boston, Massachusetts
| | - José N Onuchic
- Center for Theoretical Biological Physics, Rice University, Houston, Texas; Department of Physics and Astronomy, Department of Chemistry, Department of BioSciences, Rice University, Houston, Texas
| |
Collapse
|
23
|
Kamat K, Lao Z, Qi Y, Wang Y, Ma J, Zhang B. Compartmentalization with nuclear landmarks yields random, yet precise, genome organization. Biophys J 2023; 122:1376-1389. [PMID: 36871158 PMCID: PMC10111368 DOI: 10.1016/j.bpj.2023.03.003] [Citation(s) in RCA: 9] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/05/2022] [Revised: 02/19/2023] [Accepted: 03/01/2023] [Indexed: 03/06/2023] Open
Abstract
The 3D organization of eukaryotic genomes plays an important role in genome function. While significant progress has been made in deciphering the folding mechanisms of individual chromosomes, the principles of the dynamic large-scale spatial arrangement of all chromosomes inside the nucleus are poorly understood. We use polymer simulations to model the diploid human genome compartmentalization relative to nuclear bodies such as nuclear lamina, nucleoli, and speckles. We show that a self-organization process based on a cophase separation between chromosomes and nuclear bodies can capture various features of genome organization, including the formation of chromosome territories, phase separation of A/B compartments, and the liquid property of nuclear bodies. The simulated 3D structures quantitatively reproduce both sequencing-based genomic mapping and imaging assays that probe chromatin interaction with nuclear bodies. Importantly, our model captures the heterogeneous distribution of chromosome positioning across cells while simultaneously producing well-defined distances between active chromatin and nuclear speckles. Such heterogeneity and preciseness of genome organization can coexist due to the nonspecificity of phase separation and the slow chromosome dynamics. Together, our work reveals that the cophase separation provides a robust mechanism for us to produce functionally important 3D contacts without requiring thermodynamic equilibration that can be difficult to achieve.
Collapse
Affiliation(s)
- Kartik Kamat
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts
| | - Zhuohan Lao
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts
| | - Yifeng Qi
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts
| | - Yuchuan Wang
- Computational Biology Department, School of Computer Science, Carnegie Mellon University, Pittsburgh, Pennsylvania
| | - Jian Ma
- Computational Biology Department, School of Computer Science, Carnegie Mellon University, Pittsburgh, Pennsylvania
| | - Bin Zhang
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts.
| |
Collapse
|
24
|
Shi G, Thirumalai D. A maximum-entropy model to predict 3D structural ensembles of chromatin from pairwise distances with applications to interphase chromosomes and structural variants. Nat Commun 2023; 14:1150. [PMID: 36854665 PMCID: PMC9974990 DOI: 10.1038/s41467-023-36412-4] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/26/2022] [Accepted: 01/31/2023] [Indexed: 03/02/2023] Open
Abstract
The principles that govern the organization of genomes, which are needed for an understanding of how chromosomes are packaged and function in eukaryotic cells, could be deciphered if the three-dimensional (3D) structures are known. Recently, single-cell imaging techniques have been developed to determine the 3D coordinates of genomic loci in vivo. Here, we introduce a computational method (Distance Matrix to Ensemble of Structures, DIMES), based on the maximum entropy principle, with experimental pairwise distances between loci as constraints, to generate a unique ensemble of 3D chromatin structures. Using the ensemble of structures, we quantitatively account for the distribution of pairwise distances, three-body co-localization, and higher-order interactions. The DIMES method can be applied to both small and chromosome-scale imaging data to quantify the extent of heterogeneity and fluctuations in the shapes across various length scales. We develop a perturbation method in conjunction with DIMES to predict the changes in 3D structures from structural variations. Our method also reveals quantitative differences between the 3D structures inferred from Hi-C and those measured in imaging experiments. Finally, the physical interpretation of the parameters extracted from DIMES provides insights into the origin of phase separation between euchromatin and heterochromatin domains.
Collapse
Affiliation(s)
- Guang Shi
- Department of Chemistry, University of Texas at Austin, Austin, Texas, 78712, USA. .,Department of Materials Science, University of Illinois, Urbana, Illinois, 61801, USA.
| | - D Thirumalai
- Department of Chemistry, University of Texas at Austin, Austin, Texas, 78712, USA. .,Department of Physics, University of Texas at Austin, Austin, Texas, 78712, USA.
| |
Collapse
|
25
|
Unveiling the Machinery behind Chromosome Folding by Polymer Physics Modeling. Int J Mol Sci 2023; 24:ijms24043660. [PMID: 36835064 PMCID: PMC9967178 DOI: 10.3390/ijms24043660] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/14/2023] [Revised: 02/06/2023] [Accepted: 02/09/2023] [Indexed: 02/16/2023] Open
Abstract
Understanding the mechanisms underlying the complex 3D architecture of mammalian genomes poses, at a more fundamental level, the problem of how two or multiple genomic sites can establish physical contacts in the nucleus of the cells. Beyond stochastic and fleeting encounters related to the polymeric nature of chromatin, experiments have revealed specific, privileged patterns of interactions that suggest the existence of basic organizing principles of folding. In this review, we focus on two major and recently proposed physical processes of chromatin organization: loop-extrusion and polymer phase-separation, both supported by increasing experimental evidence. We discuss their implementation into polymer physics models, which we test against available single-cell super-resolution imaging data, showing that both mechanisms can cooperate to shape chromatin structure at the single-molecule level. Next, by exploiting the comprehension of the underlying molecular mechanisms, we illustrate how such polymer models can be used as powerful tools to make predictions in silico that can complement experiments in understanding genome folding. To this aim, we focus on recent key applications, such as the prediction of chromatin structure rearrangements upon disease-associated mutations and the identification of the putative chromatin organizing factors that orchestrate the specificity of DNA regulatory contacts genome-wide.
Collapse
|
26
|
Contessoto VG, Dudchenko O, Aiden EL, Wolynes PG, Onuchic JN, Di Pierro M. Interphase chromosomes of the Aedes aegypti mosquito are liquid crystalline and can sense mechanical cues. Nat Commun 2023; 14:326. [PMID: 36658127 PMCID: PMC9852290 DOI: 10.1038/s41467-023-35909-2] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/16/2022] [Accepted: 01/06/2023] [Indexed: 01/20/2023] Open
Abstract
We use data-driven physical simulations to study the three-dimensional architecture of the Aedes aegypti genome. Hi-C maps exhibit both a broad diagonal and compartmentalization with telomeres and centromeres clustering together. Physical modeling reveals that these observations correspond to an ensemble of 3D chromosomal structures that are folded over and partially condensed. Clustering of the centromeres and telomeres near the nuclear lamina appears to be a necessary condition for the formation of the observed structures. Further analysis of the mechanical properties of the genome reveals that the chromosomes of Aedes aegypti, by virtue of their atypical structural organization, are highly sensitive to the deformation of the nuclei. This last finding provides a possible physical mechanism linking mechanical cues to gene regulation.
Collapse
Affiliation(s)
- Vinícius G Contessoto
- Center for Theoretical Biological Physics, Rice University, Houston, TX, USA.
- Instituto de Biociências, Letras e Ciências Exatas, UNESP - Univ. Estadual Paulista, Departamento de Física, São José do Rio Preto, SP, Brazil.
| | - Olga Dudchenko
- Center for Theoretical Biological Physics, Rice University, Houston, TX, USA
- The Center for Genome Architecture, Department of Molecular and Human Genetics, Baylor College of Medicine, Houston, TX, USA
| | - Erez Lieberman Aiden
- Center for Theoretical Biological Physics, Rice University, Houston, TX, USA
- The Center for Genome Architecture, Department of Molecular and Human Genetics, Baylor College of Medicine, Houston, TX, USA
| | - Peter G Wolynes
- Center for Theoretical Biological Physics, Rice University, Houston, TX, USA
- Department of Physics & Astronomy, Rice University, Houston, TX, USA
- Department of Chemistry, Rice University, Houston, TX, USA
- Department of Biosciences, Rice University, Houston, TX, USA
| | - José N Onuchic
- Center for Theoretical Biological Physics, Rice University, Houston, TX, USA.
- Department of Physics & Astronomy, Rice University, Houston, TX, USA.
- Department of Chemistry, Rice University, Houston, TX, USA.
- Department of Biosciences, Rice University, Houston, TX, USA.
| | - Michele Di Pierro
- Department of Physics, Northeastern University, Boston, MA, USA.
- Center for Theoretical Biological Physics, Northeastern University, Boston, MA, USA.
| |
Collapse
|
27
|
Chu X, Wang J. Insights into the cell fate decision-making processes from chromosome structural reorganizations. BIOPHYSICS REVIEWS 2022; 3:041402. [PMID: 38505520 PMCID: PMC10914134 DOI: 10.1063/5.0107663] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/06/2022] [Accepted: 11/25/2022] [Indexed: 03/21/2024]
Abstract
The cell fate decision-making process, which provides the capability of a cell transition to a new cell type, involves the reorganizations of 3D genome structures. Currently, the high temporal resolution picture of how the chromosome structural rearrangements occur and further influence the gene activities during the cell-state transition is still challenging to acquire. Here, we study the chromosome structural reorganizations during the cell-state transitions among the pluripotent embryonic stem cell, the terminally differentiated normal cell, and the cancer cell using a nonequilibrium landscape-switching model implemented in the molecular dynamics simulation. We quantify the chromosome (de)compaction pathways during the cell-state transitions and find that the two pathways having the same destinations can merge prior to reaching the final states. The chromosomes at the merging states have similar structural geometries but can differ in long-range compartment segregation and spatial distribution of the chromosomal loci and genes, leading to cell-type-specific transition mechanisms. We identify the irreversible pathways of chromosome structural rearrangements during the forward and reverse transitions connecting the same pair of cell states, underscoring the critical roles of nonequilibrium dynamics in the cell-state transitions. Our results contribute to the understanding of the cell fate decision-making processes from the chromosome structural perspective.
Collapse
Affiliation(s)
- Xiakun Chu
- Advanced Materials Thrust, Function Hub, The Hong Kong University of Science and Technology (Guangzhou), Nansha, Guangzhou, Guangdong 511400, China
| | - Jin Wang
- Department of Chemistry and Physics, State University of New York at Stony Brook, Stony Brook, New York 11794, USA
| |
Collapse
|
28
|
Chu X, Wang J. Quantifying Chromosome Structural Reorganizations during Differentiation, Reprogramming, and Transdifferentiation. PHYSICAL REVIEW LETTERS 2022; 129:068102. [PMID: 36018639 DOI: 10.1103/physrevlett.129.068102] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/21/2021] [Accepted: 06/08/2022] [Indexed: 06/15/2023]
Abstract
We developed a nonequilibrium model to study chromosome structural reorganizations within a simplified cell developmental system. From the chromosome structural perspective, we predicted that the neural progenitor cell is on the neural developmental path and very close to the transdifferentiation path from the fibroblast to the neuron cell. We identified an early bifurcation of stem cell differentiation processes and the cell-of-origin-specific reprogramming pathways. Our theoretical results are in good agreement with available experimental evidence, promoting future applications of our approach.
Collapse
Affiliation(s)
- Xiakun Chu
- Advanced Materials Thrust, The Hong Kong University of Science and Technology (Guangzhou), Nansha, Guangzhou, Guangdong 511400, China
- Department of Chemistry and Physics, State University of New York at Stony Brook, Stony Brook, New York 11794, USA
| | - Jin Wang
- Center for Theoretical Interdisciplinary Sciences, Wenzhou Institute, University of Chinese Academy of Sciences, Wenzhou, Zhejiang 325001, China
- Department of Chemistry and Physics, State University of New York at Stony Brook, Stony Brook, New York 11794, USA
| |
Collapse
|
29
|
Kumari K, Ravi Prakash J, Padinhateeri R. Heterogeneous interactions and polymer entropy decide organization and dynamics of chromatin domains. Biophys J 2022; 121:2794-2812. [PMID: 35672951 PMCID: PMC9382282 DOI: 10.1016/j.bpj.2022.06.008] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/08/2022] [Revised: 04/28/2022] [Accepted: 06/01/2022] [Indexed: 11/02/2022] Open
Abstract
Chromatin is known to be organized into multiple domains of varying sizes and compaction. While these domains are often imagined as static structures, they are highly dynamic and show cell-to-cell variability. Since processes such as gene regulation and DNA replication occur in the context of these domains, it is important to understand their organization, fluctuation, and dynamics. To simulate chromatin domains, one requires knowledge of interaction strengths among chromatin segments. Here, we derive interaction-strength parameters from experimentally known contact maps and use them to predict chromatin organization and dynamics. Taking two domains on the human chromosome as examples, we investigate its three-dimensional organization, size/shape fluctuations, and dynamics of different segments within a domain, accounting for hydrodynamic effects. Considering different cell types, we quantify changes in interaction strengths and chromatin shape fluctuations in different epigenetic states. Perturbing the interaction strengths systematically, we further investigate how epigenetic-like changes can alter the spatio-temporal nature of the domains. Our results show that heterogeneous weak interactions are crucial in determining the organization of the domains. Computing effective stiffness and relaxation times, we investigate how perturbations in interactions affect the solid- and liquid-like nature of chromatin domains. Quantifying dynamics of chromatin segments within a domain, we show how the competition between polymer entropy and interaction energy influence the timescales of loop formation and maintenance of stable loops.
Collapse
Affiliation(s)
- Kiran Kumari
- IITB-Monash Research Academy, Indian Institute of Technology Bombay, Mumbai, Maharashtra, 400076, India; Department of Biosciences and Bioengineering, Indian Institute of Technology Bombay, Mumbai, Maharashtra 400076, India; Department of Chemical Engineering, Monash University, Melbourne, VIC 3800, Australia.
| | - J Ravi Prakash
- Department of Chemical Engineering, Monash University, Melbourne, VIC 3800, Australia.
| | - Ranjith Padinhateeri
- Department of Biosciences and Bioengineering, Indian Institute of Technology Bombay, Mumbai, Maharashtra 400076, India.
| |
Collapse
|
30
|
Conte M, Irani E, Chiariello AM, Abraham A, Bianco S, Esposito A, Nicodemi M. Loop-extrusion and polymer phase-separation can co-exist at the single-molecule level to shape chromatin folding. Nat Commun 2022; 13:4070. [PMID: 35831310 PMCID: PMC9279381 DOI: 10.1038/s41467-022-31856-6] [Citation(s) in RCA: 31] [Impact Index Per Article: 15.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/30/2021] [Accepted: 07/06/2022] [Indexed: 11/09/2022] Open
Abstract
Loop-extrusion and phase-separation have been proposed as mechanisms that shape chromosome spatial organization. It is unclear, however, how they perform relative to each other in explaining chromatin architecture data and whether they compete or co-exist at the single-molecule level. Here, we compare models of polymer physics based on loop-extrusion and phase-separation, as well as models where both mechanisms act simultaneously in a single molecule, against multiplexed FISH data available in human loci in IMR90 and HCT116 cells. We find that the different models recapitulate bulk Hi-C and average multiplexed microscopy data. Single-molecule chromatin conformations are also well captured, especially by phase-separation based models that better reflect the experimentally reported segregation in globules of the considered genomic loci and their cell-to-cell structural variability. Such a variability is consistent with two main concurrent causes: single-cell epigenetic heterogeneity and an intrinsic thermodynamic conformational degeneracy of folding. Overall, the model combining loop-extrusion and polymer phase-separation provides a very good description of the data, particularly higher-order contacts, showing that the two mechanisms can co-exist in shaping chromatin architecture in single cells.
Collapse
Affiliation(s)
- Mattia Conte
- Dipartimento di Fisica, Università di Napoli Federico II, and INFN Napoli, Complesso Universitario di Monte Sant'Angelo, 80126, Naples, Italy
| | - Ehsan Irani
- Berlin Institute for Medical Systems Biology, Max-Delbrück Centre (MDC) for Molecular Medicine, Berlin, Germany
- Berlin Institute of Health (BIH), MDC-Berlin, Berlin, Germany
| | - Andrea M Chiariello
- Dipartimento di Fisica, Università di Napoli Federico II, and INFN Napoli, Complesso Universitario di Monte Sant'Angelo, 80126, Naples, Italy
| | - Alex Abraham
- Dipartimento di Fisica, Università di Napoli Federico II, and INFN Napoli, Complesso Universitario di Monte Sant'Angelo, 80126, Naples, Italy
| | - Simona Bianco
- Berlin Institute for Medical Systems Biology, Max-Delbrück Centre (MDC) for Molecular Medicine, Berlin, Germany
| | - Andrea Esposito
- Dipartimento di Fisica, Università di Napoli Federico II, and INFN Napoli, Complesso Universitario di Monte Sant'Angelo, 80126, Naples, Italy
| | - Mario Nicodemi
- Dipartimento di Fisica, Università di Napoli Federico II, and INFN Napoli, Complesso Universitario di Monte Sant'Angelo, 80126, Naples, Italy.
- Berlin Institute for Medical Systems Biology, Max-Delbrück Centre (MDC) for Molecular Medicine, Berlin, Germany.
- Berlin Institute of Health (BIH), MDC-Berlin, Berlin, Germany.
| |
Collapse
|
31
|
Contessoto VG, Cheng RR, Onuchic JN. Uncovering the statistical physics of 3D chromosomal organization using data-driven modeling. Curr Opin Struct Biol 2022; 75:102418. [PMID: 35839701 DOI: 10.1016/j.sbi.2022.102418] [Citation(s) in RCA: 8] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/28/2022] [Revised: 05/04/2022] [Accepted: 06/03/2022] [Indexed: 11/03/2022]
Abstract
In recent years, much effort has been devoted to understanding the three-dimensional (3D) organization of the genome and how genomic structure mediates nuclear function. The development of experimental techniques that combine DNA proximity ligation with high-throughput sequencing, such as Hi-C, have substantially improved our knowledge about chromatin organization. Numerous experimental advancements, not only utilizing DNA proximity ligation but also high-resolution genome imaging (DNA tracing), have required theoretical modeling to determine the structural ensembles consistent with such data. These 3D polymer models of the genome provide an understanding of the physical mechanisms governing genome architecture. Here, we present an overview of the recent advances in modeling the ensemble of 3D chromosomal structures by employing the maximum entropy approach combined with polymer physics. Particularly, we discuss the minimal chromatin model (MiChroM) along with the "maximum entropy genomic annotations from biomarkers associated with structural ensembles" (MEGABASE) model, which have been remarkably successful in the accurate modeling of chromosomes consistent with both Hi-C and DNA-tracing data.
Collapse
Affiliation(s)
- Vinícius G Contessoto
- Center for Theoretical Biological Physics, Rice University, Houston, TX, USA. https://twitter.com/Vini_Contessoto
| | - Ryan R Cheng
- Center for Theoretical Biological Physics, Rice University, Houston, TX, USA. https://twitter.com/ryanrcheng
| | - José N Onuchic
- Center for Theoretical Biological Physics, Rice University, Houston, TX, USA; Department of Physics and Astronomy, Rice University, Houston, TX, USA; Department of Chemistry, Rice University, Houston, TX, USA; Department of Biosciences, Rice University, Houston, TX, USA.
| |
Collapse
|
32
|
Lamberti WF, Zang C. Extracting physical characteristics of higher-order chromatin structures from 3D image data. Comput Struct Biotechnol J 2022; 20:3387-3398. [PMID: 35832633 PMCID: PMC9260447 DOI: 10.1016/j.csbj.2022.06.018] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/15/2022] [Revised: 06/09/2022] [Accepted: 06/09/2022] [Indexed: 11/25/2022] Open
Abstract
Higher-order chromatin structures have functional impacts on gene regulation and cell identity determination. Using high-throughput sequencing (HTS)-based methods like Hi-C, active or inactive compartments and open or closed topologically associating domain (TAD) structures can be identified on a cell population level. Recently developed high-resolution three-dimensional (3D) molecular imaging techniques such as 3D electron microscopy with in situ hybridization (3D-EMSIH) and 3D structured illumination microscopy (3D-SIM) enable direct detection of physical representations of chromatin structures in a single cell. However, computational analysis of 3D image data with explainability and interpretability on functional characteristics of chromatin structures is still challenging. We developed Extracting Physical-Characteristics from Images of Chromatin Structures (EPICS), a machine-learning based computational method for processing high-resolution chromatin 3D image data. Using EPICS on images produced by 3D-EMISH or 3D-SIM techniques, we generated more direct 3D representations of higher-order chromatin structures, identified major chromatin domains, and determined the open or closed status of each domain. We identified several high-contributing features from the model as the major physical characteristics that define the open or closed chromatin domains, demonstrating the explainability and interpretability of EPICS. EPICS can be applied to the analysis of other high-resolution 3D molecular imaging data for spatial genomics studies. The R and Python codes of EPICS are available at https://github.com/zang-lab/epics.
Collapse
Affiliation(s)
- William Franz Lamberti
- Center for Public Health Genomics, University of Virginia, Charlottesville, VA 22908, USA
| | - Chongzhi Zang
- Center for Public Health Genomics, University of Virginia, Charlottesville, VA 22908, USA
- Department of Public Health Sciences, University of Virginia, Charlottesville, VA 22908, USA
| |
Collapse
|
33
|
Yildirim A, Boninsegna L, Zhan Y, Alber F. Uncovering the Principles of Genome Folding by 3D Chromatin Modeling. Cold Spring Harb Perspect Biol 2022; 14:a039693. [PMID: 34400556 PMCID: PMC9248826 DOI: 10.1101/cshperspect.a039693] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/25/2022]
Abstract
Our understanding of how genomic DNA is tightly packed inside the nucleus, yet is still accessible for vital cellular processes, has grown dramatically over recent years with advances in microscopy and genomics technologies. Computational methods have played a pivotal role in the structural interpretation of experimental data, which helped unravel some organizational principles of genome folding. Here, we give an overview of current computational efforts in mechanistic and data-driven 3D chromatin structure modeling. We discuss strengths and limitations of different methods and evaluate the added value and benefits of computational approaches to infer the 3D structural and dynamic properties of the genome and its underlying mechanisms at different scales and resolution, ranging from the dynamic formation of chromatin loops and topological associated domains to nuclear compartmentalization of chromatin and nuclear bodies.
Collapse
Affiliation(s)
- Asli Yildirim
- Institute for Quantitative and Computational Biosciences, Department of Microbiology, Immunology and Molecular Genetics, University of California Los Angeles, Los Angeles, California 90095, USA
| | - Lorenzo Boninsegna
- Institute for Quantitative and Computational Biosciences, Department of Microbiology, Immunology and Molecular Genetics, University of California Los Angeles, Los Angeles, California 90095, USA
| | - Yuxiang Zhan
- Institute for Quantitative and Computational Biosciences, Department of Microbiology, Immunology and Molecular Genetics, University of California Los Angeles, Los Angeles, California 90095, USA
- Quantitative and Computational Biology, Department of Biological Sciences, University of Southern California, Los Angeles, California 90089, USA
| | - Frank Alber
- Institute for Quantitative and Computational Biosciences, Department of Microbiology, Immunology and Molecular Genetics, University of California Los Angeles, Los Angeles, California 90095, USA
- Quantitative and Computational Biology, Department of Biological Sciences, University of Southern California, Los Angeles, California 90089, USA
| |
Collapse
|
34
|
Bera P, Wasim A, Mondal J. Hi-C embedded polymer model of Escherichia coli reveals the origin of heterogeneous subdiffusion in chromosomal loci. Phys Rev E 2022; 105:064402. [PMID: 35854496 DOI: 10.1103/physreve.105.064402] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/29/2021] [Accepted: 05/10/2022] [Indexed: 06/15/2023]
Abstract
Underneath its apparently simple architecture, the circular chromosome of Escherichia coli is known for displaying complex dynamics in its cytoplasm, with past investigations hinting at inherently diverse mobilities of chromosomal loci across the genome. To decipher its origin, we simulate the dynamics of genome-wide spectrum of E. coli chromosomal loci, via integrating its experimentally derived Hi-C interaction matrix within a polymer-based model. Our analysis demonstrates that, while the dynamics of the chromosome is subdiffusive in a viscoelastic media, the diffusion constants are strongly dependent of chromosomal loci coordinates and diffusive exponents (α) are widely heterogenous with α ≈ 0.36-0.60. The loci-dependent heterogeneous dynamics and mean first-passage times of interloci encounter were found to be modulated via genetically distant interloci communications and is robust even in the presence of active, ATP-dependent noises. Control investigations reveal that the absence of Hi-C-derived interactions in the model would have abolished the traits of heterogeneous loci diffusion, underscoring the key role of loci-specific genetically distant interaction in modulating the underlying heterogeneity of the loci diffusion.
Collapse
Affiliation(s)
- Palash Bera
- Tata Institute of Fundamental Research, Hyderabad 500046, India
| | - Abdul Wasim
- Tata Institute of Fundamental Research, Hyderabad 500046, India
| | | |
Collapse
|
35
|
Generation of dynamic three-dimensional genome structure through phase separation of chromatin. Proc Natl Acad Sci U S A 2022; 119:e2109838119. [PMID: 35617433 PMCID: PMC9295772 DOI: 10.1073/pnas.2109838119] [Citation(s) in RCA: 17] [Impact Index Per Article: 8.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
Significance DNA functions in living cells are crucially affected by the three-dimensional genome structure and dynamics. We analyze the whole genome of human cells by developing a polymer model of interphase nuclei. The model reveals the essential importance of the unfolding process of chromosomes from the condensed mitotic state for describing the interphase nuclei; through the unfolding process, heterogeneous repulsive interactions among chromatin chains induce phase separation of chromatin, which quantitatively explains the experimentally observed various genomic data. We can use this model structure as a platform to analyze the relationship among genome structure, dynamics, and functions.
Collapse
|
36
|
Amini Farsani Z, Schmid VJ. Maximum Entropy Technique and Regularization Functional for Determining the Pharmacokinetic Parameters in DCE-MRI. J Digit Imaging 2022; 35:1176-1188. [PMID: 35618849 PMCID: PMC9582183 DOI: 10.1007/s10278-022-00646-3] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/04/2021] [Revised: 04/15/2022] [Accepted: 04/21/2022] [Indexed: 10/31/2022] Open
Abstract
This paper aims to solve the arterial input function (AIF) determination in dynamic contrast-enhanced MRI (DCE-MRI), an important linear ill-posed inverse problem, using the maximum entropy technique (MET) and regularization functionals. In addition, estimating the pharmacokinetic parameters from a DCE-MR image investigations is an urgent need to obtain the precise information about the AIF-the concentration of the contrast agent on the left ventricular blood pool measured over time. For this reason, the main idea is to show how to find a unique solution of linear system of equations generally in the form of [Formula: see text] named an ill-conditioned linear system of equations after discretization of the integral equations, which appear in different tomographic image restoration and reconstruction issues. Here, a new algorithm is described to estimate an appropriate probability distribution function for AIF according to the MET and regularization functionals for the contrast agent concentration when applying Bayesian estimation approach to estimate two different pharmacokinetic parameters. Moreover, by using the proposed approach when analyzing simulated and real datasets of the breast tumors according to pharmacokinetic factors, it indicates that using Bayesian inference-that infer the uncertainties of the computed solutions, and specific knowledge of the noise and errors-combined with the regularization functional of the maximum entropy problem, improved the convergence behavior and led to more consistent morphological and functional statistics and results. Finally, in comparison to the proposed exponential distribution based on MET and Newton's method, or Weibull distribution via the MET and teaching-learning-based optimization (MET/TLBO) in the previous studies, the family of Gamma and Erlang distributions estimated by the new algorithm are more appropriate and robust AIFs.
Collapse
Affiliation(s)
- Zahra Amini Farsani
- Bayesian Imaging and Spatial Statistics Group, Institute of Statistics, Ludwig-Maximilian-Universität München, Ludwigstraße 33, 80539, Munich, Germany. .,Statistics Department, School of Science, Lorestan University, 68151-44316, Khorramabad, Iran.
| | - Volker J Schmid
- Bayesian Imaging and Spatial Statistics Group, Institute of Statistics, Ludwig-Maximilian-Universität München, Ludwigstraße 33, 80539, Munich, Germany
| |
Collapse
|
37
|
Gong W, Wee J, Wu MC, Sun X, Li C, Xia K. Persistent spectral simplicial complex-based machine learning for chromosomal structural analysis in cellular differentiation. Brief Bioinform 2022; 23:6583209. [PMID: 35536545 DOI: 10.1093/bib/bbac168] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/21/2022] [Revised: 04/12/2022] [Accepted: 03/13/2022] [Indexed: 11/13/2022] Open
Abstract
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.
Collapse
Affiliation(s)
- Weikang Gong
- Faculty of Environmental and Life Sciences, Beijing University of Technology, Beijing, China 100124.,Division of Mathematical Sciences, School of Physical and Mathematical Sciences, Nanyang Technological University, Singapore 637371
| | - JunJie Wee
- Division of Mathematical Sciences, School of Physical and Mathematical Sciences, Nanyang Technological University, Singapore 637371
| | - Min-Chun Wu
- Division of Mathematical Sciences, School of Physical and Mathematical Sciences, Nanyang Technological University, Singapore 637371
| | - Xiaohan Sun
- Faculty of Environmental and Life Sciences, Beijing University of Technology, Beijing, China 100124
| | - Chunhua Li
- Faculty of Environmental and Life Sciences, Beijing University of Technology, Beijing, China 100124
| | - Kelin Xia
- Division of Mathematical Sciences, School of Physical and Mathematical Sciences, Nanyang Technological University, Singapore 637371
| |
Collapse
|
38
|
Chu W, Chu X, Wang J. Uncovering the Quantitative Relationships Among Chromosome Fluctuations, Epigenetics, and Gene Expressions of Transdifferentiation on Waddington Landscape. ADVANCED SCIENCE (WEINHEIM, BADEN-WURTTEMBERG, GERMANY) 2022; 9:e2103617. [PMID: 35104056 PMCID: PMC8981899 DOI: 10.1002/advs.202103617] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 08/27/2021] [Revised: 12/21/2021] [Indexed: 06/14/2023]
Abstract
The 3D spatial organization of the chromosomes appears to be linked to the gene function, which is cell type-specific. The chromosome structural ensemble switching model (CSESM) is developed by employing a heteropolymer model on different cell types and the important quantitative relationships among the chromosome ensemble, the epigenetic marks, and the gene expressions are uncovered, that both chromosome fluctuation and epigenetic marks have strong linear correlations with the gene expressions. The results support that the two compartments have different behaviors, corresponding to the relatively sparse and fluctuating phase (compartment A) and the relatively dense and stable phase (compartment B). Importantly, through the investigation of the transdifferentiation processes between the peripheral blood mononuclear cell (PBMC) and the bipolar neuron (BN), a quantitative description for the transdifferentiation is provided, which can be linked to the Waddington landscape. In addition, compared to the direct transdifferentiation between PBMC and BN, the transdifferentiation via the intermediate state neural progenitor cell (NPC) follows a different path (an "uphill" followed by a "downhill"). These theoretical studies bridge the gap among the chromosome fluctuations/ensembles, the epigenetics, and gene expressions in determining the cell fate.
Collapse
Affiliation(s)
- Wen‐Ting Chu
- State Key Laboratory of Electroanalytical ChemistryChangchun Institute of Applied ChemistryChinese Academy of SciencesChangchunJilin130022China
| | - Xiakun Chu
- Department of Chemistry & PhysicsState University of New York at Stony BrookStony BrookNY11794USA
| | - Jin Wang
- Department of Chemistry & PhysicsState University of New York at Stony BrookStony BrookNY11794USA
| |
Collapse
|
39
|
Perspectives on the landscape and flux theory for describing emergent behaviors of the biological systems. J Biol Phys 2022; 48:1-36. [PMID: 34822073 PMCID: PMC8866630 DOI: 10.1007/s10867-021-09586-5] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/20/2021] [Accepted: 09/07/2021] [Indexed: 10/19/2022] Open
Abstract
We give a review on the landscape theory of the equilibrium biological systems and landscape-flux theory of the nonequilibrium biological systems as the global driving force. The emergences of the behaviors, the associated thermodynamics in terms of the entropy and free energy and dynamics in terms of the rate and paths have been quantitatively demonstrated. The hierarchical organization structures have been discussed. The biological applications ranging from protein folding, biomolecular recognition, specificity, biomolecular evolution and design for equilibrium systems as well as cell cycle, differentiation and development, cancer, neural networks and brain function, and evolution for nonequilibrium systems, cross-scale studies of genome structural dynamics and experimental quantifications/verifications of the landscape and flux are illustrated. Together, this gives an overall global physical and quantitative picture in terms of the landscape and flux for the behaviors, dynamics and functions of biological systems.
Collapse
|
40
|
Chu X, Wang J. Dynamics and Pathways of Chromosome Structural Organizations during Cell Transdifferentiation. JACS AU 2022; 2:116-127. [PMID: 35098228 PMCID: PMC8791059 DOI: 10.1021/jacsau.1c00416] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 09/21/2021] [Indexed: 06/14/2023]
Abstract
Direct conversion of one differentiated cell type into another is defined as cell transdifferentiation. In avoidance of forming pluripotency, cell transdifferentiation can reduce the potential risk of tumorigenicity, thus offering significant advantages over cell reprogramming in clinical applications. Until now, the mechanism of cell transdifferentiation is still largely unknown. It has been well recognized that cell transdifferentiation is determined by the underlying gene expression regulation, which relies on the accurate adaptation of the chromosome structure. To dissect the transdifferentiation at the molecular level, we develop a nonequilibrium landscape-switching model to investigate the chromosome structural dynamics during the state transitions between the human fibroblast and neuron cells. We uncover the high irreversibility of the transdifferentiation at the local chromosome structural ranges, where the topologically associating domains form. In contrast, the pathways in the two opposite directions of the transdifferentiation projected onto the chromosome compartment profiles are highly overlapped, indicating that the reversibility vanishes at the long-range chromosome structures. By calculating the contact strengths in the chromosome at the states along the paths, we observe strengthening contacts in compartment A concomitant with weakening contacts in compartment B at the early stages of the transdifferentiation. This further leads to adapting contacts toward the ones at the embryonic stem cell. In light of the intimate structure-function relationship at the chromosomal level, we suggest an increase of "stemness" during the transdifferentiation. In addition, we find that the neuron progenitor cell (NPC), a cell developmental state, is located on the transdifferentiation pathways projected onto the long-range chromosome contacts. The findings are consistent with the previous single-cell RNA sequencing experiment, where the NPC-like cell states were observed during the direct conversion of the fibroblast to neuron cells. Thus, we offer a promising microscopic and physical approach to study the cell transdifferentiation mechanism from the chromosome structural perspective.
Collapse
Affiliation(s)
- Xiakun Chu
- Department
of Chemistry, State University of New York
at Stony Brook, Stony
Brook, New York 11794, United States
| | - Jin Wang
- Department
of Chemistry, State University of New York
at Stony Brook, Stony
Brook, New York 11794, United States
- Department
of Physics and Astronomy, State University
of New York at Stony Brook, Stony
Brook, New York 11794, United States
| |
Collapse
|
41
|
Modified Maximum Entropy Method and Estimating the AIF via DCE-MRI Data Analysis. ENTROPY 2022; 24:e24020155. [PMID: 35205451 PMCID: PMC8871336 DOI: 10.3390/e24020155] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 12/06/2021] [Revised: 01/16/2022] [Accepted: 01/17/2022] [Indexed: 02/06/2023]
Abstract
Background: For the kinetic models used in contrast-based medical imaging, the assignment of the arterial input function named AIF is essential for the estimation of the physiological parameters of the tissue via solving an optimization problem. Objective: In the current study, we estimate the AIF relayed on the modified maximum entropy method. The effectiveness of several numerical methods to determine kinetic parameters and the AIF is evaluated—in situations where enough information about the AIF is not available. The purpose of this study is to identify an appropriate method for estimating this function. Materials and Methods: The modified algorithm is a mixture of the maximum entropy approach with an optimization method, named the teaching-learning method. In here, we applied this algorithm in a Bayesian framework to estimate the kinetic parameters when specifying the unique form of the AIF by the maximum entropy method. We assessed the proficiency of the proposed method for assigning the kinetic parameters in the dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI), when determining AIF with some other parameter-estimation methods and a standard fixed AIF method. A previously analyzed dataset consisting of contrast agent concentrations in tissue and plasma was used. Results and Conclusions: We compared the accuracy of the results for the estimated parameters obtained from the MMEM with those of the empirical method, maximum likelihood method, moment matching (“method of moments”), the least-square method, the modified maximum likelihood approach, and our previous work. Since the current algorithm does not have the problem of starting point in the parameter estimation phase, it could find the best and nearest model to the empirical model of data, and therefore, the results indicated the Weibull distribution as an appropriate and robust AIF and also illustrated the power and effectiveness of the proposed method to estimate the kinetic parameters.
Collapse
|
42
|
Boninsegna L, Yildirim A, Zhan Y, Alber F. Integrative approaches in genome structure analysis. Structure 2022; 30:24-36. [PMID: 34963059 PMCID: PMC8959402 DOI: 10.1016/j.str.2021.12.003] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/17/2021] [Revised: 11/13/2021] [Accepted: 12/01/2021] [Indexed: 12/17/2022]
Abstract
New technological advances in integrated imaging, sequencing-based assays, and computational analysis have revolutionized our view of genomes in terms of their structure and dynamics in space and time. These advances promise a deeper understanding of genome functions and mechanistic insights into how the nucleus is spatially organized and functions. These wide arrays of complementary data provide an opportunity to produce quantitative integrative models of nuclear organization. In this article, we highlight recent key developments and discuss the outlook for these fields.
Collapse
Affiliation(s)
- Lorenzo Boninsegna
- Institute for Quantitative and Computational Biosciences, University of California Los Angeles, Los Angeles, CA 90095, USA; Department of Microbiology, Immunology and Molecular Genetics, University of California Los Angeles, Los Angeles, CA 90095, USA
| | - Asli Yildirim
- Institute for Quantitative and Computational Biosciences, University of California Los Angeles, Los Angeles, CA 90095, USA; Department of Microbiology, Immunology and Molecular Genetics, University of California Los Angeles, Los Angeles, CA 90095, USA
| | - Yuxiang Zhan
- Institute for Quantitative and Computational Biosciences, University of California Los Angeles, Los Angeles, CA 90095, USA; Department of Microbiology, Immunology and Molecular Genetics, University of California Los Angeles, Los Angeles, CA 90095, USA; Department of Quantitative and Computational Biology, University of Southern California, Los Angeles, CA 90089, USA
| | - Frank Alber
- Institute for Quantitative and Computational Biosciences, University of California Los Angeles, Los Angeles, CA 90095, USA; Department of Microbiology, Immunology and Molecular Genetics, University of California Los Angeles, Los Angeles, CA 90095, USA; Department of Quantitative and Computational Biology, University of Southern California, Los Angeles, CA 90089, USA.
| |
Collapse
|
43
|
Natesan R, Gowrishankar K, Kuttippurathu L, Kumar PBS, Rao M. Active Remodeling of Chromatin and Implications for In Vivo Folding. J Phys Chem B 2021; 126:100-109. [DOI: 10.1021/acs.jpcb.1c08655] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Ramakrishnan Natesan
- Department of Cancer Biology, University of Pennsylvania, Philadelphia, Pennsylvania 19104, United States
- Department of Physics, Indian Institute of Technology Madras, Chennai 600036, India
| | | | - Lakshmi Kuttippurathu
- Department of Physics, Indian Institute of Technology Madras, Chennai 600036, India
- Daniel Baugh Institute for Functional Genomics and Computational Biology, Department of Pathology, Anatomy and Cell Biology, Thomas Jefferson University, Philadelphia, Pennsylvania 19107, United States
| | - P. B. Sunil Kumar
- Department of Physics, Indian Institute of Technology Madras, Chennai 600036, India
- Department of Physics, Indian Institute of Technology Palakkad, Palakkad 668557, Kerala, India
| | - Madan Rao
- Simons Centre for the Study of Living Machines, National Centre for Biological Sciences (TIFR), Bengaluru 560065, India
| |
Collapse
|
44
|
Chu WT, Yan Z, Chu X, Zheng X, Liu Z, Xu L, Zhang K, Wang J. Physics of biomolecular recognition and conformational dynamics. REPORTS ON PROGRESS IN PHYSICS. PHYSICAL SOCIETY (GREAT BRITAIN) 2021; 84:126601. [PMID: 34753115 DOI: 10.1088/1361-6633/ac3800] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/26/2021] [Accepted: 11/09/2021] [Indexed: 06/13/2023]
Abstract
Biomolecular recognition usually leads to the formation of binding complexes, often accompanied by large-scale conformational changes. This process is fundamental to biological functions at the molecular and cellular levels. Uncovering the physical mechanisms of biomolecular recognition and quantifying the key biomolecular interactions are vital to understand these functions. The recently developed energy landscape theory has been successful in quantifying recognition processes and revealing the underlying mechanisms. Recent studies have shown that in addition to affinity, specificity is also crucial for biomolecular recognition. The proposed physical concept of intrinsic specificity based on the underlying energy landscape theory provides a practical way to quantify the specificity. Optimization of affinity and specificity can be adopted as a principle to guide the evolution and design of molecular recognition. This approach can also be used in practice for drug discovery using multidimensional screening to identify lead compounds. The energy landscape topography of molecular recognition is important for revealing the underlying flexible binding or binding-folding mechanisms. In this review, we first introduce the energy landscape theory for molecular recognition and then address four critical issues related to biomolecular recognition and conformational dynamics: (1) specificity quantification of molecular recognition; (2) evolution and design in molecular recognition; (3) flexible molecular recognition; (4) chromosome structural dynamics. The results described here and the discussions of the insights gained from the energy landscape topography can provide valuable guidance for further computational and experimental investigations of biomolecular recognition and conformational dynamics.
Collapse
Affiliation(s)
- Wen-Ting Chu
- State Key Laboratory of Electroanalytical Chemistry, Changchun Institute of Applied Chemistry, Chinese Academy of Sciences, Changchun 130022, People's Republic of China
| | - Zhiqiang Yan
- State Key Laboratory of Electroanalytical Chemistry, Changchun Institute of Applied Chemistry, Chinese Academy of Sciences, Changchun 130022, People's Republic of China
| | - Xiakun Chu
- Department of Chemistry & Physics, State University of New York at Stony Brook, Stony Brook, NY 11794, United States of America
| | - Xiliang Zheng
- State Key Laboratory of Electroanalytical Chemistry, Changchun Institute of Applied Chemistry, Chinese Academy of Sciences, Changchun 130022, People's Republic of China
| | - Zuojia Liu
- State Key Laboratory of Electroanalytical Chemistry, Changchun Institute of Applied Chemistry, Chinese Academy of Sciences, Changchun 130022, People's Republic of China
| | - Li Xu
- State Key Laboratory of Electroanalytical Chemistry, Changchun Institute of Applied Chemistry, Chinese Academy of Sciences, Changchun 130022, People's Republic of China
| | - Kun Zhang
- State Key Laboratory of Electroanalytical Chemistry, Changchun Institute of Applied Chemistry, Chinese Academy of Sciences, Changchun 130022, People's Republic of China
| | - Jin Wang
- Department of Chemistry & Physics, State University of New York at Stony Brook, Stony Brook, NY 11794, United States of America
| |
Collapse
|
45
|
Abstract
Nuclear bodies are membraneless condensates that may form via liquid-liquid phase separation. The viscoelastic chromatin network could impact their stability and may hold the key for understanding experimental observations that defy predictions of classical theories. However, quantitative studies on the role of the chromatin network in phase separation have remained challenging. Using a diploid human genome model parameterized with chromosome conformation capture (Hi-C) data, we study the thermodynamics and kinetics of nucleoli formation. Dynamical simulations predict the formation of multiple droplets for nucleolar particles that experience specific interactions with nucleolus-associated domains (NADs). Coarsening dynamics, surface tension, and coalescence kinetics of the simulated droplets are all in quantitative agreement with experimental measurements for nucleoli. Free energy calculations further support that a two-droplet state, often observed for nucleoli in somatic cells, is metastable and separated from the single-droplet state with an entropic barrier. Our study suggests that nucleoli-chromatin interactions facilitate droplets' nucleation but hinder their coarsening due to the coupled motion between droplets and the chromatin network: as droplets coalesce, the chromatin network becomes increasingly constrained. Therefore, the chromatin network supports a nucleation and arrest mechanism to stabilize the multi-droplet state for nucleoli and possibly for other nuclear bodies.
Collapse
|
46
|
Stochastic chromatin packing of 3D mitotic chromosomes revealed by coherent X-rays. Proc Natl Acad Sci U S A 2021; 118:2109921118. [PMID: 34750262 DOI: 10.1073/pnas.2109921118] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 10/01/2021] [Indexed: 11/18/2022] Open
Abstract
DNA molecules are atomic-scale information storage molecules that promote reliable information transfer via fault-free repetitions of replications and transcriptions. Remarkable accuracy of compacting a few-meters-long DNA into a micrometer-scale object, and the reverse, makes the chromosome one of the most intriguing structures from both physical and biological viewpoints. However, its three-dimensional (3D) structure remains elusive with challenges in observing native structures of specimens at tens-of-nanometers resolution. Here, using cryogenic coherent X-ray diffraction imaging, we succeeded in obtaining nanoscale 3D structures of metaphase chromosomes that exhibited a random distribution of electron density without characteristics of high-order folding structures. Scaling analysis of the chromosomes, compared with a model structure having the same density profile as the experimental results, has discovered the fractal nature of density distributions. Quantitative 3D density maps, corroborated by molecular dynamics simulations, reveal that internal structures of chromosomes conform to diffusion-limited aggregation behavior, which indicates that 3D chromatin packing occurs via stochastic processes.
Collapse
|
47
|
Chu X, Wang J. Deciphering the molecular mechanism of the cancer formation by chromosome structural dynamics. PLoS Comput Biol 2021; 17:e1009596. [PMID: 34752443 PMCID: PMC8631624 DOI: 10.1371/journal.pcbi.1009596] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/02/2021] [Revised: 11/30/2021] [Accepted: 10/28/2021] [Indexed: 12/15/2022] Open
Abstract
Cancer reflects the dysregulation of the underlying gene network, which is strongly related to the 3D genome organization. Numerous efforts have been spent on experimental characterizations of the structural alterations in cancer genomes. However, there is still a lack of genomic structural-level understanding of the temporal dynamics for cancer initiation and progression. Here, we use a landscape-switching model to investigate the chromosome structural transition during the cancerization and reversion processes. We find that the chromosome undergoes a non-monotonic structural shape-changing pathway with initial expansion followed by compaction during both of these processes. Furthermore, our analysis reveals that the chromosome with a more expanding structure than those at both the normal and cancer cell during cancerization exhibits a sparse contact pattern, which shows significant structural similarity to the one at the embryonic stem cell in many aspects, including the trend of contact probability declining with the genomic distance, the global structural shape geometry and the spatial distribution of loci on the chromosome. In light of the intimate structure-function relationship at the chromosomal level, we further describe the cell state transition processes by the chromosome structural changes, suggesting an elevated cell stemness during the formation of the cancer cells. We show that cell cancerization and reversion are highly irreversible processes in terms of the chromosome structural transition pathways, spatial repositioning of chromosomal loci and hysteresis loop of contact evolution analysis. Our model draws a molecular-scale picture of cell cancerization from the chromosome structural perspective. The process contains initial reprogramming towards the stem cell followed by the differentiation towards the cancer cell, accompanied by an initial increase and subsequent decrease of the cell stemness.
Collapse
Affiliation(s)
- Xiakun Chu
- Department of Chemistry, State University of New York at Stony Brook, Stony Brook, New York, United States of America
| | - Jin Wang
- Department of Chemistry, State University of New York at Stony Brook, Stony Brook, New York, United States of America
- Department of Physics and Astronomy, State University of New York at Stony Brook, Stony Brook, New York, United States of America
| |
Collapse
|
48
|
Liang J, Perez-Rathke A. Minimalistic 3D chromatin models: Sparse interactions in single cells drive the chromatin fold and form many-body units. Curr Opin Struct Biol 2021; 71:200-214. [PMID: 34399301 DOI: 10.1016/j.sbi.2021.06.017] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/05/2021] [Revised: 06/27/2021] [Accepted: 06/29/2021] [Indexed: 11/26/2022]
Abstract
Computational three-dimensional chromatin modeling has helped uncover principles of genome organization. Here, we discuss methods for modeling three-dimensional chromatin structures, with focus on a minimalistic polymer model which inverts population Hi-C into single-cell conformations. Utilizing only basic physical properties, this model reveals that a few specific Hi-C interactions can fold chromatin into conformations consistent with single-cell imaging, Dip-C, and FISH measurements. Aggregated single-cell chromatin conformations also reproduce Hi-C frequencies. This approach allows quantification of structural heterogeneity and discovery of many-body interaction units and has revealed additional insights, including (1) topologically associating domains as a byproduct of folding driven by specific interactions, (2) cell subpopulations with different structural scaffolds are developmental stage dependent, and (3) the functional landscape of many-body units within enhancer-rich regions. We also discuss these findings in relation to the genome structure-function relationship.
Collapse
Affiliation(s)
- Jie Liang
- Center for Bioinformatics and Quantitative Biology & Richard and Loan Hill Department of Bioengineering, University of Illinois at Chicago, Chicago, IL, 60612, USA.
| | - Alan Perez-Rathke
- Center for Bioinformatics and Quantitative Biology & Richard and Loan Hill Department of Bioengineering, University of Illinois at Chicago, Chicago, IL, 60612, USA
| |
Collapse
|
49
|
Wei J, Tian H, Zhou R, Shao Y, Song F, Gao YQ. Topological Constraints with Optimal Length Promote the Formation of Chromosomal Territories at Weakened Degree of Phase Separation. J Phys Chem B 2021; 125:9092-9101. [PMID: 34351763 DOI: 10.1021/acs.jpcb.1c03523] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
It is generally agreed that the nuclei of eukaryotic cells at interphase are partitioned into disjointed territories, with distinct regions occupied by certain chromosomes. However, the underlying mechanism for such territorialization is still under debate. Here we model chromosomes as coarse-grained block copolymers and to investigate the effect of loop domains (LDs) on the formation of compartments and territories based on dissipative particle dynamics. A critical length of LDs, which depends sensitively on the length of polymeric blocks, is obtained to minimize the degree of phase separation. This also applies to the two-polymer system: The critical length not only maximizes the degree of territorialization but also minimizes the degree of phase separation. Interestingly, by comparing with experimental data, we find the critical length for LDs and the corresponding length of blocks to be respectively very close to the mean length of topologically associating domains (TADs) and chromosomal segments with different densities of CpG islands for human chromosomes. The results indicate that topological constraints with optimal length can contribute to the formation of territories by weakening the degree of phase separation, which likely promotes the chromosomal flexibility in response to genetic regulations.
Collapse
Affiliation(s)
- Jiachen Wei
- State Key Laboratory of Nonlinear Mechanics and Beijing Key Laboratory of Engineered Construction and Mechanobiology, Institute of Mechanics, Chinese Academy of Sciences, Beijing 100190, China.,Shenzhen Bay Laboratory, 5F, No. 9 Duxue Road, Nanshan District, 518055 Shenzhen, Guangdong, China
| | - Hao Tian
- Beijing National Laboratory for Molecular Sciences, College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, China.,Biomedical Pioneering Innovation Center (BIOPIC), Peking University, Beijing 100871, China.,Beijing Advanced Innovation Center for Genomics (ICG), Peking University, Beijing 100871, China
| | - Rui Zhou
- Beijing National Laboratory for Molecular Sciences, College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, China.,Biomedical Pioneering Innovation Center (BIOPIC), Peking University, Beijing 100871, China.,Beijing Advanced Innovation Center for Genomics (ICG), Peking University, Beijing 100871, China
| | - Yingfeng Shao
- State Key Laboratory of Nonlinear Mechanics and Beijing Key Laboratory of Engineered Construction and Mechanobiology, Institute of Mechanics, Chinese Academy of Sciences, Beijing 100190, China.,School of Engineering Science, University of Chinese Academy of Sciences, Beijing 100049, China
| | - Fan Song
- State Key Laboratory of Nonlinear Mechanics and Beijing Key Laboratory of Engineered Construction and Mechanobiology, Institute of Mechanics, Chinese Academy of Sciences, Beijing 100190, China.,School of Engineering Science, University of Chinese Academy of Sciences, Beijing 100049, China
| | - Yi Qin Gao
- Shenzhen Bay Laboratory, 5F, No. 9 Duxue Road, Nanshan District, 518055 Shenzhen, Guangdong, China.,Beijing National Laboratory for Molecular Sciences, College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, China.,Biomedical Pioneering Innovation Center (BIOPIC), Peking University, Beijing 100871, China.,Beijing Advanced Innovation Center for Genomics (ICG), Peking University, Beijing 100871, China
| |
Collapse
|
50
|
Oliveira Junior AB, Estrada CP, Aiden EL, Contessoto VG, Onuchic JN. Chromosome Modeling on Downsampled Hi-C Maps Enhances the Compartmentalization Signal. J Phys Chem B 2021; 125:8757-8767. [PMID: 34319725 DOI: 10.1021/acs.jpcb.1c04174] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The human genome is organized within a nucleus where chromosomes fold into an ensemble of different conformations. Chromosome conformation capture techniques such as Hi-C provide information about the genome architecture by creating a 2D heat map. Initially, Hi-C map experiments were performed in human interphase cell lines. Recently, efforts were expanded to several different organisms, cell lines, tissues, and cell cycle phases where obtaining high-quality maps is challenging. Poor sampled Hi-C maps present high sparse matrices where compartments located far from the main diagonal are difficult to observe. Aided by recently developed models for chromatin folding and dynamics investigation, we introduce a framework to enhance the compartments' information far from the diagonal observed in experimental sparse matrices. The simulations were performed using the Open-MiChroM platform aided by new trained parameters in the minimal chromatin model (MiChroM) energy function. The simulations optimized on a downsampled experimental map (10% of the original data) allow the prediction of a contact frequency similar to that of the complete (100%) experimental Hi-C. The modeling results open a discussion on how simulations and modeling can increase the statistics and help fill in some Hi-C regions not captured by poor sampling experiments. Open-MiChroM simulations allow us to explore the 3D genome organization of different organisms, cell lines, and cell phases that often do not produce high-quality Hi-C maps.
Collapse
Affiliation(s)
| | - Cynthia Perez Estrada
- The Center for Genome Architecture, Department of Molecular and Human Genetics, Baylor College of Medicine, Houston, Texas, United States
| | - Erez Lieberman Aiden
- The Center for Genome Architecture, Department of Molecular and Human Genetics, Baylor College of Medicine, Houston, Texas, United States
| | - Vinícius G Contessoto
- Instituto de Biociências, Letras e Ciências Exatas, UNESP - Univ. Estadual Paulista, Departamento de Física, São José do Rio Preto, SP, Brazil
| | | |
Collapse
|