Unveiling Nucleosome Dynamics: A Comparative Study Using All-Atom and Coarse-Grained Simulations Enhanced by Principal Component Analysis

Abstract

This study investigates nucleosome dynamics using both all-atom and coarse-grained (CG) molecular dynamics simulations, focusing on the SIRAH force field. Simulations are performed for two nucleosomal DNA sequences--ASP and Widom-601--over six microseconds at physiological salt concentrations. Comparative analysis of structural parameters, such as groove widths and base pair geometries, reveals good agreement between atomistic and CG models, though CG simulations exhibit broader conformational sampling and greater breathing motion of DNA ends. Principal component analysis (PCA) is applied to DNA structural parameters, revealing multiple free energy minima, especially in CG simulations. These findings highlight the potential of the SIRAH CG force field for studying large-scale nucleosome dynamics, offering insights into DNA repositioning and sequence-dependent behavior.