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

Overview

Journal bioRxiv

Date 2025 Feb 3

PMID 39896498

Authors

Rutika Patel

Sharon M Loverde

Affiliations

Soon will be listed here.

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.

References

Biswas M, Voltz K, Smith J, Langowski J . Role of histone tails in structural stability of the nucleosome. PLoS Comput Biol. 2011; 7(12):e1002279. PMC: 3240580. DOI: 10.1371/journal.pcbi.1002279. View

Beauchamp K, McGibbon R, Lin Y, Pande V . Simple few-state models reveal hidden complexity in protein folding. Proc Natl Acad Sci U S A. 2012; 109(44):17807-13. PMC: 3497769. DOI: 10.1073/pnas.1201810109. View

Potoyan D, Papoian G . Regulation of the H4 tail binding and folding landscapes via Lys-16 acetylation. Proc Natl Acad Sci U S A. 2012; 109(44):17857-62. PMC: 3497739. DOI: 10.1073/pnas.1201805109. View

Huertas J, MacCarthy C, Scholer H, Cojocaru V . Nucleosomal DNA Dynamics Mediate Oct4 Pioneer Factor Binding. Biophys J. 2020; 118(9):2280-2296. PMC: 7202942. DOI: 10.1016/j.bpj.2019.12.038. View

Lee J, Wei S, Lee T . Effects of histone acetylation by Piccolo NuA4 on the structure of a nucleosome and the interactions between two nucleosomes. J Biol Chem. 2011; 286(13):11099-109. PMC: 3064164. DOI: 10.1074/jbc.M110.192047. View

Song F, Chen P, Sun D, Wang M, Dong L, Liang D . Cryo-EM study of the chromatin fiber reveals a double helix twisted by tetranucleosomal units. Science. 2014; 344(6182):376-80. DOI: 10.1126/science.1251413. View

Chang L, Takada S . Histone acetylation dependent energy landscapes in tri-nucleosome revealed by residue-resolved molecular simulations. Sci Rep. 2016; 6:34441. PMC: 5048180. DOI: 10.1038/srep34441. View

Bascom G, Schlick T . Chromatin Fiber Folding Directed by Cooperative Histone Tail Acetylation and Linker Histone Binding. Biophys J. 2018; 114(10):2376-2385. PMC: 6129558. DOI: 10.1016/j.bpj.2018.03.008. View

Malmstrom R, Kornev A, Taylor S, Amaro R . Allostery through the computational microscope: cAMP activation of a canonical signalling domain. Nat Commun. 2015; 6:7588. PMC: 4504738. DOI: 10.1038/ncomms8588. View

10.

Voelz V, Bowman G, Beauchamp K, Pande V . Molecular simulation of ab initio protein folding for a millisecond folder NTL9(1-39). J Am Chem Soc. 2010; 132(5):1526-8. PMC: 2835335. DOI: 10.1021/ja9090353. View

11.

Bowman G, Pande V . Protein folded states are kinetic hubs. Proc Natl Acad Sci U S A. 2010; 107(24):10890-5. PMC: 2890711. DOI: 10.1073/pnas.1003962107. View

12.

Patel R, Onyema A, Tang P, Loverde S . Conformational Dynamics of the Nucleosomal Histone H2B Tails Revealed by Molecular Dynamics Simulations. J Chem Inf Model. 2024; 64(12):4709-4726. PMC: 11200259. DOI: 10.1021/acs.jcim.4c00059. View

13.

Reynoird N, Schwartz B, Delvecchio M, Sadoul K, Meyers D, Mukherjee C . Oncogenesis by sequestration of CBP/p300 in transcriptionally inactive hyperacetylated chromatin domains. EMBO J. 2010; 29(17):2943-52. PMC: 2944051. DOI: 10.1038/emboj.2010.176. View

14.

Prinz J, Wu H, Sarich M, Keller B, Senne M, Held M . Markov models of molecular kinetics: generation and validation. J Chem Phys. 2011; 134(17):174105. DOI: 10.1063/1.3565032. View

15.

Tian C, Kasavajhala K, Belfon K, Raguette L, Huang H, Migues A . ff19SB: Amino-Acid-Specific Protein Backbone Parameters Trained against Quantum Mechanics Energy Surfaces in Solution. J Chem Theory Comput. 2019; 16(1):528-552. DOI: 10.1021/acs.jctc.9b00591. View

16.

Galindo-Murillo R, Robertson J, Zgarbova M, Sponer J, Otyepka M, Jurecka P . Assessing the Current State of Amber Force Field Modifications for DNA. J Chem Theory Comput. 2016; 12(8):4114-27. PMC: 4980684. DOI: 10.1021/acs.jctc.6b00186. View

17.

Singhal N, Snow C, Pande V . Using path sampling to build better Markovian state models: predicting the folding rate and mechanism of a tryptophan zipper beta hairpin. J Chem Phys. 2004; 121(1):415-25. DOI: 10.1063/1.1738647. View

18.

Kim T, Nosella M, Bolik-Coulon N, Harkness R, Huang S, Kay L . Correlating histone acetylation with nucleosome core particle dynamics and function. Proc Natl Acad Sci U S A. 2023; 120(15):e2301063120. PMC: 10104578. DOI: 10.1073/pnas.2301063120. View

19.

Kulic I, Schiessel H . Chromatin dynamics: nucleosomes go mobile through twist defects. Phys Rev Lett. 2003; 91(14):148103. DOI: 10.1103/PhysRevLett.91.148103. View

20.

Lee K, HAYES J . The N-terminal tail of histone H2A binds to two distinct sites within the nucleosome core. Proc Natl Acad Sci U S A. 1997; 94(17):8959-64. PMC: 22982. DOI: 10.1073/pnas.94.17.8959. View