The formation of a tetrameric assembly is essential for the ability of the tumor suppressor protein p53 to act as a transcription factor. Such a quaternary conformation is driven by a specific tetramerization domain, separated from the central DNA-binding domain by a flexible linker. Despite the distance, functional crosstalk between the two domains has been reported. This phenomenon can explain the pathogenicity of some inherited or somatically acquired mutations in the tetramerization domain, including the widespread R337H missense mutation present in the population in south Brazil. In this work, we combined computational predictions through extended all-atom molecular dynamics simulations with functional assays in a genetically defined yeast-based model system to reveal structural features of p53 tetramerization domains and their transactivation capacity and specificity. In addition to the germline and cancer-associated R337H and R337C, other rationally designed missense mutations targeting a significant salt-bridge interaction that stabilizes the p53 tetramerization domain were studied (i.e., R337D, D352R, and the double-mutation R337D plus D352R). The simulations revealed a destabilizing effect of the pathogenic mutations within the p53 tetramerization domain and highlighted the importance of electrostatic interactions between residues 337 and 352. The transactivation assay, performed in yeast by tuning the expression of wild-type and mutant p53 proteins, revealed that p53 tetramerization mutations could decrease the transactivation potential and alter transactivation specificity, in particular by better tolerating negative features in weak DNA-binding sites. These results establish the effect of naturally occurring variations at positions 337 and 352 on p53's conformational stability and function.
Keywords: R337H; TP53; molecular dynamics; tetramerization domain; transactivation assays.