Dual Role of Protein Phosphorylation in DNA Activator/Coactivator Binding

Binding free energies are calculated for the phosphorylated and unphosphorylated complexes between the kinase inducible domain (KID) of the DNA transcriptional activator cAMP response element binding (CREB) protein and the KIX domain of its coactivator, CREB-binding protein (CBP). To our knowledge, this is the first application of a method based on a potential of mean force (PMF) with restraining potentials to compute the binding free energy of protein-protein complexes. The KID:KIX complexes are chosen here because of their biological relevance to the DNA transcription process and their relatively small size (81 residues for the KIX domain of CBP, and 28 residues for KID). The results for pKID:KIX and KID:KIX are −9.55 and −4.96 kcal/mol, respectively, in good agreement with experimental estimates (−8.8 and −5.8 kcal/mol, respectively). A comparison between specific contributions to protein-protein binding for the phosphorylated and unphosphorylated complexes reveals a dual role for the phosphorylation of KID at Ser-133 in effecting a more favorable free energy of the bound system: 1), stabilization of the unbound conformation of phosphorylated KID due to favorable intramolecular interactions of the phosphate group of Ser-133 with the charged groups of an arginine-rich region spanning both α-helices, which lowers the configurational entropy; and 2), more favorable intermolecular electrostatic interactions between pSer-133 and Arg-131 of KID, and Lys-662, Tyr-658, and Glu-666 of KIX. Charge reduction through ligand phosphorylation emerges as a possible mechanism for controlling the unbound state conformation of KID and, ultimately, gene expression. This work also demonstrates that the PMF-based method with restraining potentials provides an added benefit in that important elements of the binding pathway are evidenced. Furthermore, the practicality of the PMF-based method for larger systems is validated by agreement with experiment. In addition, we provide a somewhat differently structured exposition of the PMF-based method with restraining potentials and outline its generalization to systems in which both protein and ligand may adopt unbound conformations that are different from those of the bound state

Decomposition of Protein Experimental Compressibility into Intrinsic and Hydration Shell Contributions

The experimental determination of protein compressibility reflects both the protein intrinsic compressibility and the difference between the compressibility of water in the protein hydration shell and bulk water. We use molecular dynamics simulations to explore the dependence of the isothermal compressibility of the hydration shell surrounding globular proteins on differential contributions from charged, polar, and apolar protein-water interfaces. The compressibility of water in the protein hydration shell is accounted for by a linear combination of contributions from charged, polar, and apolar solvent-accessible surfaces. The results provide a formula for the deconvolution of experimental data into intrinsic and hydration contributions when a protein of known structure is investigated. The physical basis for the model is the variation in water density shown by the surface-specific radial distribution functions of water molecules around globular proteins. The compressibility of water hydrating charged atoms is lower than bulk water compressibility, the compressibility of water hydrating apolar atoms is somewhat larger than bulk water compressibility, and the compressibility of water around polar atoms is about the same as the compressibility of bulk water. We also assess whether hydration water compressibility determined from small compound data can be used to estimate the compressibility of hydration water surrounding proteins. The results, based on an analysis from four dipeptide solutions, indicate that small compound data cannot be used directly to estimate the compressibility of hydration water surrounding proteins