Unveiling the Conformational Dynamics of the Histone Tails Using Markov State Modeling

Abstract

Abstract Biomolecules predominantly exert their function through altering conformational dynamics. The nucleosome core particle (NCP) is the fundamental unit of chromatin. DNA with ∼146 base pairs wrap around the histone octamer to form a nucleosome. The histone octamer is comprised of two copies of each histone protein (H3, H4, H2A, and H2B) with a globular core and disordered N-terminal tails. Epigenetic modifications of the histone N-terminal tails play a critical role in the regulation of chromatin structure and biological processes such as transcription and DNA repair. Here, we report all-atomistic molecular dynamics (MD) simulations of the nucleosome at microsecond timescales to construct Markov state models (MSMs) to elucidate distinct conformations of the histone tails. We employ the time-lagged independent component analysis (tICA) to capture their essential slow dynamics, with k-means clustering used to discretize the conformational space. MSMs unveil distinct states and transition probabilities to characterize the dynamics and kinetics of the tails. Next, we focus on the H2B tail, one of the least studied tails. We show that acetylation increases secondary structure formation, with an increase in transition rates. These findings will aid in understanding the functional implications of tail conformations in nucleosome stability and gene regulation. Significance Statement The nucleosome is the fundamental repeat unit of chromatin, composed of a histone octamer and DNA. Each histone has a globular core and N-terminal tail regions. N-terminal tails are major sites for post-translational modifications, chromatin structure, and gene regulation. Using all-atom molecular dynamics simulations of the nucleosome to examine the dynamics of histone tails and analyze the effects of one of the histone modifications, such as acetylation on histone tails. We explore tail dynamics and kinetics using Markov State Models (MSMs) to gain insight into histone tail structural changes. Acetylation of the H2B tail shows an increased dynamics of the tail based on transition rates. Our study characterizes transition rates of all histone tails that can influence nucleosome plasticity and gene regulation.