1
|
Lao Z, Kamat K, Jiang Z, Zhang B. OpenNucleome for high resolution nuclear structural and dynamical modeling. bioRxiv 2024:2023.10.16.562451. [PMID: 37905090 PMCID: PMC10614770 DOI: 10.1101/2023.10.16.562451] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/02/2023]
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, MA, USA
| | - Kartik Kamat
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, MA, USA
| | - Zhongling Jiang
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, MA, USA
| | - Bin Zhang
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, MA, USA
| |
Collapse
|
2
|
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: 8] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [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
|
3
|
Kamat K, Naullage PM, Molinero V, Peters B. Oriented attachment kinetics for rod-like particles at a flat surface: Buffon's needle at the nanoscale. J Chem Phys 2022; 157:214113. [PMID: 36511557 DOI: 10.1063/5.0124531] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/23/2022] Open
Abstract
The adsorption of large rod-like molecules or crystallites on a flat crystal face, similar to Buffon's needle, requires the rods to "land," with their binding sites in precise orientational alignment with matching sites on the surface. An example is provided by long, helical antifreeze proteins (AFPs), which bind at specific facets and orientations on the ice surface. The alignment constraint for adsorption, in combination with the loss in orientational freedom as the molecule diffuses toward the surface, results in an entropic barrier that hinders the adsorption. Prior kinetic models do not factor in the complete geometry of the molecule, nor explicitly enforce orientational constraints for adsorption. Here, we develop a diffusion-controlled adsorption theory for AFP molecules binding at specific orientations to flat ice surfaces. We formulate the diffusion equation with relevant boundary conditions and present analytical solutions to the attachment rate constant. The resulting rate constant is a function of the length and aspect ratio of the AFP, the distance threshold associated with binding, and solvent conditions such as temperature and viscosity. These results and methods of calculation may also be useful for predicting the kinetics of crystal growth through oriented attachment.
Collapse
Affiliation(s)
- Kartik Kamat
- Department of Chemical Engineering, University of California Santa Barbara, Santa Barbara, California 93106, USA
| | - Pavithra M Naullage
- Department of Chemistry, The University of Utah, Salt Lake City, Utah 84112, USA
| | - Valeria Molinero
- Department of Chemistry, The University of Utah, Salt Lake City, Utah 84112, USA
| | - Baron Peters
- Chemical and Biomolecular Engineering, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, USA
| |
Collapse
|
4
|
Abstract
The human genome is arranged in the cell nucleus nonrandomly, and phase separation has been proposed as an important driving force for genome organization. However, the cell nucleus is an active system, and the contribution of nonequilibrium activities to phase separation and genome structure and dynamics remains to be explored. We simulated the genome using an energy function parametrized with chromosome conformation capture (Hi-C) data with the presence of active, nondirectional forces that break the detailed balance. We found that active forces that may arise from transcription and chromatin remodeling can dramatically impact the spatial localization of heterochromatin. When applied to euchromatin, active forces can drive heterochromatin to the nuclear envelope and compete with passive interactions among heterochromatin that tend to pull them in opposite directions. Furthermore, active forces induce long-range spatial correlations among genomic loci beyond single chromosome territories. We further showed that the impact of active forces could be understood from the effective temperature defined as the fluctuation-dissipation ratio. Our study suggests that nonequilibrium activities can significantly impact genome structure and dynamics, producing unexpected collective phenomena.
Collapse
Affiliation(s)
- Zhongling Jiang
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139-4307, United States
| | - Yifeng Qi
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139-4307, United States
| | - Kartik Kamat
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139-4307, United States
| | - Bin Zhang
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139-4307, United States
| |
Collapse
|
5
|
Abstract
Some of the most potent antifreeze proteins (AFPs) are approximately rigid helical structures that bind with one side in contact with the ice surface at specific orientations. These AFPs take random orientations in solution; however, most orientations become sterically inaccessible as the AFP approaches the ice surface. The effect of these inaccessible orientations on the rate of adsorption of AFP to ice has never been explored. Here, we present a diffusion-controlled theory of adsorption kinetics that accounts for these orientational restrictions to predict a rate constant for adsorption (kon, in m/s) as a function of the length and width of the AFP molecules. We find that kon decreases with length and diameter of the AFP and is almost proportional to the inverse of the area of the binding surface. We demonstrate that the restricted orientations create an entropic barrier to AFP adsorption, which we compute to be approximately 7 kBT for most AFPs and up to 9 kBT for Maxi, the largest known AFP. We compare the entropic resistance 1/kon to resistances for diffusion through boundary layers and across typical distances in the extracellular matrix and find that these entropic and diffusion resistances could become comparable in the small confined spaces of biological environments.
Collapse
Affiliation(s)
- Kartik Kamat
- Department of Chemical Engineering, University of California, Santa Barbara, Santa Barbara, California 93106, United States
| | - Pavithra M Naullage
- Department of Chemistry, The University of Utah, Salt Lake City, Utah 84112, United States
| | - Valeria Molinero
- Department of Chemistry, The University of Utah, Salt Lake City, Utah 84112, United States
| | - Baron Peters
- Chemical and Biomolecular Engineering, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, United States.,Department of Chemistry and Biochemistry, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, United States
| |
Collapse
|
6
|
Kamat K, Guo R, Reutzel-Edens SM, Price SL, Peters B. Diabat method for polymorph free energies: Extension to molecular crystals. J Chem Phys 2020; 153:244105. [PMID: 33380078 DOI: 10.1063/5.0024727] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Lattice-switch Monte Carlo and the related diabat methods have emerged as efficient and accurate ways to compute free energy differences between polymorphs. In this work, we introduce a one-to-one mapping from the reference positions and displacements in one molecular crystal to the positions and displacements in another. Two features of the mapping facilitate lattice-switch Monte Carlo and related diabat methods for computing polymorph free energy differences. First, the mapping is unitary so that its Jacobian does not complicate the free energy calculations. Second, the mapping is easily implemented for molecular crystals of arbitrary complexity. We demonstrate the mapping by computing free energy differences between polymorphs of benzene and carbamazepine. Free energy calculations for thermodynamic cycles, each involving three independently computed polymorph free energy differences, all return to the starting free energy with a high degree of precision. The calculations thus provide a force field independent validation of the method and allow us to estimate the precision of the individual free energy differences.
Collapse
Affiliation(s)
- Kartik Kamat
- Department of Chemical Engineering, University of California-Santa Barbara, Santa Barbara, California 93106, USA
| | - Rui Guo
- Department of Chemistry, University College London, London WC1H 0AJ, United Kingdom
| | - Susan M Reutzel-Edens
- Small Molecule Design and Development, Eli Lilly and Company, Indianapolis, Indiana 46285, USA
| | - Sarah L Price
- Department of Chemistry, University College London, London WC1H 0AJ, United Kingdom
| | - Baron Peters
- Department of Chemical and Biomolecular Engineering, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, USA
| |
Collapse
|
7
|
Yappert R, Kamat K, Peters B. The overdamped transmission coefficient: Recovering the true mean first passage time from an inaccurate reaction coordinate. J Chem Phys 2019; 151:184108. [DOI: 10.1063/1.5117237] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022] Open
Affiliation(s)
- Ryan Yappert
- Department of Chemical and Biomolecular Engineering, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, USA
| | - Kartik Kamat
- Department of Chemical Engineering, University of California, Santa Barbara, California 93106, USA
| | - Baron Peters
- Department of Chemical and Biomolecular Engineering, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, USA
| |
Collapse
|
8
|
Abstract
Polymorph free-energy differences are critical to several applications. A recently proposed diabat interpolation framework estimated free-energy differences between polymorphs by quadratic interpolation of diabats. This work extends the Zwanzig-Bennett relation to the NPT ensemble so that the diabats directly give Gibbs free-energy differences. We also demonstrate how the approach can be used in cases where the diabats are not parabolic. We illustrate the diabat method for Gibbs free-energy difference of zirconium (BCC and HCP phases) and compare it with the conventional lattice switch Monte Carlo approach.
Collapse
Affiliation(s)
- Kartik Kamat
- Department of Chemical Engineering, University of California Santa Barbara, Santa Barbara, California 93106, USA
| | - Baron Peters
- Department of Chemical Engineering, University of California Santa Barbara, Santa Barbara, California 93106, USA
| |
Collapse
|
9
|
Abstract
Existing methods to compute free-energy differences between polymorphs use harmonic approximations, advanced non-Boltzmann bias sampling techniques, and/or multistage free-energy perturbations. This work demonstrates how Bennett's diabat interpolation method ( J. Comput. Phys. 1976, 22, 245 ) can be combined with energy gaps from lattice-switch Monte Carlo techniques ( Phys. Rev. E 2000, 61, 906 ) to swiftly estimate polymorph free-energy differences. The new method requires only two unbiased molecular dynamics simulations, one for each polymorph. To illustrate the new method, we compute the free-energy difference between face-centered cubic and body-centered cubic polymorphs for a Gaussian core solid. We discuss the justification for parabolic models of the free-energy diabats and similarities to methods that have been used in studies of electron transfer.
Collapse
Affiliation(s)
- Kartik Kamat
- Department of Chemical Engineering and ‡Department of Chemistry and Biochemistry, University of California , Santa Barbara, California 93106, United States
| | - Baron Peters
- Department of Chemical Engineering and ‡Department of Chemistry and Biochemistry, University of California , Santa Barbara, California 93106, United States
| |
Collapse
|