Abstract
Human small Heat Shock Protein 1 (HSPB1) belongs to the Small Heat Shock Protein (sHSP) superfamily, a group of ATP-independent molecular chaperones essential for cellular stress responses and protein quality control. These proteins share a conserved domain organization, with a structured Alpha-Crystallin domain (ACD) flanked by disordered N-terminal and C-terminal regions (NTR and CTR). While the prevailing evolutionary hypothesis for the sHSP family suggests that the disordered regions evolved independently and at a faster rate than the ACD, this study provides, for the first time, evidence of coevolution between these regions in human HSPB1, introducing new insights into the evolutionary mechanisms that sustain critical regulatory interactions. By integrating evolutionary and structural approaches, we estimated evolutionary rates per region and position, analyzed the composition of key interacting motifs, and employed structural modeling with AlphaFold 2 to assess the prevalence of these interactions. Our findings reveal that while the disordered regions globally evolve faster than the ACD, specific motifs involved in regulatory interactions exhibit lower-than-average evolutionary rates, reflecting evolutionary constraints imposed by their functional importance. This coevolutionary mechanism may also extend to other small Heat Shock Proteins featuring interacting motifs in the NTR, CTR, or both, offering a new perspective for studying their molecular evolution. Furthermore, the analysis presented in this work could be applied to assess coevolution in other proteins with intrinsically disordered regions.