The Molecular and Evolutionary Principles of Histone Folding in Eukarya and Archaea

Histones are the dominant proteins to compact and store DNA in both Eukarya and Archaea. For a long time, histones are observed to exist in the unit of dimers but diverge into different formats such as heterodimers in Eukarya or homodimers in Archaea. Here, by studying 11 types of histone proteins, both monomers and their dimeric complexes, using multiscale molecular dynamics (MD) simulations combined with NMR and circular dichroism experiments, we confirm the widely applied "folding upon binding" mechanism of histone structures. A histone dimer appears to form the longest 2 helices followed by other shorter helices and inter-molecular tertiary structures. We report an alternative conformation, namely, the inverted non-native dimer, which has a minimum free energy state. Protein sequence analysis indicates that the inverted conformation can be attributed to a hidden head-tail sequence symmetry underlying all histone proteins. This finding strongly support previously proposed histone evolution hypotheses. Finally, we separately used the MD-based AWSEM and AI-based AlphaFold-Multimer model to predict eukaryotic histone homodimer structures and performed extensive allatom MD simulations to examine their structural stabilities. Our results suggest that eukaryotic histones can also form stable homodimers, whereas their disordered tails-- the structurally asymmetrical region--may tip the balance towards the formation of heterotypic dimers.