The electrostatic allostery could be the trigger for the changes in dynamics for the PDZ domain of PICK1

The PDZ domain is a highly abundant protein-protein interaction domain that exists in many signaling proteins, such as PICK1. Despite the highly conserved structure of the PDZ family, the PDZ family has an extremely low sequence identity, making each PDZ domain unique. PICK1 is the only protein in the human genome that is comprised of a PDZ domain and a BAR domain. PICK1 regulates surface membrane proteins and has been identified as an integral player in drug addiction. Like many PDZ-containing proteins, PICK1 is positively regulated by its PDZ domain and has thus drawn attention to be a potential drug target to curb the effects of substance abuse. The goal of this study is to use all-atom molecular dynamics simulations and the electrostatic analysis program, DelPhi, to better understand the unique interactions and dynamic changes in the PICK1 PDZ domain upon complex formation. Our results demonstrated that the PICK1 PDZ domain shares similar canonical PDZ-ligand hydrogen bonding networks and fluctuations of the carboxylate-binding loop to other PDZ domains. Furthermore, our results are unique to the PICK1 PDZ domain as we reveal that the binding of ligand opens up the binding pocket and, at the same time, reduces the fluctuations of both the central part of the binding pocket and the short loop region between the A-helix and {beta}C-strand. More importantly, the binding of ligand resulted in charge redistribution at the binding pocket region as well as the N- and C-termini of the PDZ domain that are not a part of the binding pocket. These results suggest that the electrostatic allostery resulted from ligand binding could be the key factor leading to the changes in dynamics which may be associated with the activation of PICK1. Based on these results, an effective drug to target PDZ domain must not only stably bind to the PICK1 PDZ domain but also prevent the electrostatic allostery of the PDZ domain.