Electrostatic properties at the interface of p53 Transactivation domain binding

Intrinsically disordered proteins (IDPs) are an abundant class of proteins and protein regions which rapidly change between multiple structures without an equilibrium position. IDPs exist as a series of conformational ensembles of semi-stable conformations that can be adopted based on a hilly landscape of shallow free energy minima. Disordered sequences share characteristic features differentiating them from globular proteins, including low sequence complexity, low occurrence of hydrophobic residues, high polar and charged residue content, and high flexibility. IDPs are commonly involved in regulation in the cell, and frequently function as, or interact with, hub proteins in protein-protein interaction networks, making them an important class of macromolecules for understanding regulatory and other processes. Given their functional importance, these proteins are widely studied. Many analytical techniques are used, though rapid conformational sampling by IDPs makes it difficult to detect many states with NMR or other techniques. Computational approaches such as molecular dynamics are increasingly used to probe the binding and conformational sampling of these proteins, allowing for observation of factors that cannot be observed with traditional analytical methods such as NMR, such as differing conformational ensembles and the dipoles of individual residues. Here, we studied the role of electrostatic interactions in IDP protein-protein interaction using molecular dynamics simulations performed with the Drude-2019 force field (FF), a polarizable model that allows for more accurate representation of electrostatics, an important factor for highly charged systems like IDPs. For this project, a prototypical protein with disordered regions, p53, was simulated with two protein partners, the nuclear coactivator domain of the CREB binding protein (CBP), and the E3 ubiquitin-protein ligase mouse double minute 2 (MDM2). p53 is widely studied, and the p53 transactivation domain (TAD) is disordered and binds to many structurally diverse partners, making this protein domain a useful model for probing the role of electrostatic interactions formed by IDPs at protein-protein binding interfaces. We found that the Drude-2019 FF allows for simulation of the p53 TAD with Cα chemical shifts comparable to those observed with NMR, supporting that the Drude-2019 FF performs well in simulating IDPs. We observed large relative change in sidechain dipole moments when comparing the p53 TAD alone and when bound to either CBP or MDM2. We observed that aliphatic and aromatic amino acids experienced the largest relative change in sidechain dipole moments, and that there is sensitivity to binding shown in this dipole response. The largest percent changes in sidechain dipole moment were found to localize at and around the binding interface. Understanding the binding interactions of IDPs at a fundamental level, including the role of electrostatic interactions, may help with targeting IDPs or their partners for drug design.