Multiscale modeling reveals the ion-mediated phase separation of nucleosome core particles.

Due to the vast length scale inside the cell nucleus, multiscale models are required to understand chromatin folding, structure, and dynamics and how they regulate genomic activities such as DNA transcription, replication, and repair. We study the interactions and structure of condensed phases formed by the universal building block of chromatin, the nucleosome core particle (NCP), using bottom-up multiscale coarse-grained (CG) simulations with a model extracted from all-atom MD simulations. In the presence of the multivalent cations Mg(H2O)62+ or CoHex3+, we analyze the internal structures of the NCP aggregates and the contributions of histone tails and ions to the aggregation patterns. We then derive a "super" coarse-grained (SCG) NCP model to study the macroscopic scale phase separation of NCPs. The SCG simulations show the formation of NCP aggregates with Mg(H2O)62+ concentration-dependent densities and sizes. Variation of the CoHex3+ concentrations results in highly ordered lamellocolumnar and hexagonal columnar phases in agreement with experimental data. The results give detailed insights into nucleosome interactions and for understanding chromatin folding in the cell nucleus.

A Bottom-Up Coarse-Grained Model for Nucleosome-Nucleosome Interactions with Explicit Ions.

The nucleosome core particle (NCP) is a large complex of 145-147 base pairs of DNA and eight histone proteins and is the basic building block of chromatin that forms the chromosomes. Here, we develop a coarse-grained (CG) model of the NCP derived through a systematic bottom-up approach based on underlying all-atom MD simulations to compute the necessary CG interactions. The model produces excellent agreement with known structural features of the NCP and gives a realistic description of the nucleosome-nucleosome attraction in the presence of multivalent cations (Mg(H2O)62+ or Co(NH3)63+) for systems comprising 20 NCPs. The results of the simulations reveal structural details of the NCP-NCP interactions unavailable from experimental approaches, and this model opens the prospect for the rigorous modeling of chromatin fibers.

Bottom-Up Coarse-Grained Modeling of DNA.

Recent advances in methodology enable effective coarse-grained modeling of deoxyribonucleic acid (DNA) based on underlying atomistic force field simulations. The so-called bottom-up coarse-graining practice separates fast and slow dynamic processes in molecular systems by averaging out fast degrees of freedom represented by the underlying fine-grained model. The resulting effective potential of interaction includes the contribution from fast degrees of freedom effectively in the form of potential of mean force. The pair-wise additive potential is usually adopted to construct the coarse-grained Hamiltonian for its efficiency in a computer simulation. In this review, we present a few well-developed bottom-up coarse-graining methods, discussing their application in modeling DNA properties such as DNA flexibility (persistence length), conformation, "melting," and DNA condensation.

Modeling DNA Flexibility: Comparison of Force Fields from Atomistic to Multiscale Levels.

Accurate parametrization of force fields (FFs) is of ultimate importance for computer simulations to be reliable and to possess a predictive power. In this work, we analyzed, in multi-microsecond simulations of a 40-base-pair DNA fragment, the performance of four force fields, namely, the two recent major updates of CHARMM and two from the AMBER family. We focused on a description of double-helix DNA flexibility and dynamics both at atomistic and at mesoscale level in coarse-grained (CG) simulations. In addition to the traditional analysis of different base-pair and base-step parameters, we extended our analysis to investigate the ability of the force field to parametrize a CG DNA model by structure-based bottom-up coarse-graining, computing DNA persistence length as a function of ionic strength. Our simulations unambiguously showed that the CHARMM36 force field is unable to preserve DNA's structural stability at over-microsecond time scale. Both versions of the AMBER FF, parmbsc0 and parmbsc1, showed good agreement with experiment, with some bias of parmbsc0 parameters for intermediate A/B form DNA structures. The CHARMM27 force field provides stable atomistic trajectories and overall (among the considered force fields) the best fit to experimentally determined DNA flexibility parameters both at atomistic and at mesoscale level.

A multiscale analysis of DNA phase separation: from atomistic to mesoscale level.

DNA condensation and phase separation is of utmost importance for DNA packing in vivo with important applications in medicine, biotechnology and polymer physics. The presence of hexagonally ordered DNA is observed in virus capsids, sperm heads and in dinoflagellates. Rigorous modelling of this process in all-atom MD simulations is presently difficult to achieve due to size and time scale limitations. We used a hierarchical approach for systematic multiscale coarse-grained (CG) simulations of DNA phase separation induced by the three-valent cobalt(III)-hexammine (CoHex3+). Solvent-mediated effective potentials for a CG model of DNA were extracted from all-atom MD simulations. Simulations of several hundred 100-bp-long CG DNA oligonucleotides in the presence of explicit CoHex3+ ions demonstrated aggregation to a liquid crystalline hexagonally ordered phase. Following further coarse-graining and extraction of effective potentials, we conducted modelling at mesoscale level. In agreement with electron microscopy observations, simulations of an 10.2-kb-long DNA molecule showed phase separation to either a toroid or a fibre with distinct hexagonal DNA packing. The mechanism of toroid formation is analysed in detail. The approach used here is based only on the underlying all-atom force field and uses no adjustable parameters and may be generalised to modelling chromatin up to chromosome size.