Abstract
Huntington's disease (HD) is a neurological condition caused by an excessive expansion of CAG repeats in the Huntingtin (HTT) gene. Although experiments have shown an altered epigenetic landscape and chromatin architecture upon HD development, the structural consequences on the HTT gene remain elusive. Structural data are only available for model nucleosome systems and yeast systems with human nucleosomes. Here, we use our experimentally validated nucleosome-resolution mesoscale chromatin model to investigate folding changes of the HTT gene associated with HD. We investigate how the histone fold domain of the variant macroH2A1, a biomarker of HD, affects the genome structure by modeling HD-like systems that contain (i) 100% canonical, (ii) 100% macroH2A1, (iii) 50% canonical and 50% macroH2A1, and (iv) 100% hybrid cores (one canonical H2A and one macroH2A1 per nucleosome). Then, we model the mouse HTT gene in healthy and HD conditions by incorporating the CAG expansion and macroH2A1 cores, reducing the linker histone density and tail acetylation levels, and incorporating genomic contacts. Overall, our results show that the histone fold domain of macroH2A1 affects chromatin compaction in a fiber-dependent manner (i.e., nucleosome distribution dependent) and can thus both enhance or repress HTT gene expression. Our modeling of the HTT gene shows that HTT is less compact in the diseased condition, which could accelerate the production of the mutated protein. By suggesting the structural biophysical consequences of the HTT gene under HD conditions, our findings may help in the development of diagnostic and therapeutic treatments for HD.