Abstract
Intrinsically disordered proteins (IDPs) are key components of cellular signaling and regulatory networks. They frequently remain dynamic even in complexes and thus rely on potentially subtle shifts in the disordered conformational ensemble for function. Understanding the molecular basis of these fascinating mechanisms of IDP function and regulation requires a detailed characterization of dynamic ensembles in various biologically relevant states. Here, we study the phosphorylation dependence of the dynamic interaction between the N-terminal transactivation domain (NTAD) and DNA-binding domain (DBD) of tumor suppressor p53, which plays a key role in the autoinhibition and regulation of p53 activation or termination during various stages of stress response. By extending the hybrid-resolution (HyRes) coarse-grained (CG) protein force field to model phosphorylated side chains, we show that HyRes simulations accurately recapitulate the effects of phosphorylation on the p53 NTAD/DBD interactions. The simulated ensembles show that phosphorylation of Thr55 as well as Ser46 enhances dynamic NTAD/DBD interactions and further induces conformational shifts that promote trans interactions between two p53 dimers to drive dissociation from DNA. These CG simulations thus provide a strong molecular basis in support of previous experimental studies suggesting the central role of dynamic interactions of disordered domains and phosphorylation in the function of p53. The success of this study also suggests that HyRes provides an efficient and viable tool for studying dynamic interactions and post-translational modifications in IDP function and regulation.