A fluid salt-bridging cluster and the stabilization of p53

p53 is a homotetrameric tumor suppressor protein that is found to be mutated in most human cancers. Some of these mutations, particularly mutations to R337, fall in the tetramerization domain and cause defects in tetramer formation leading to loss of function. Mutation to His at this site has been found to destabilize the tetramer in a pH-dependent fashion. In structures of the tetramerization domain determined by crystallography, R337 from one monomer makes a salt bridge with D352 from another monomer, apparently helping to stabilize the tetramer. Here we present molecular dynamics simulations of wild-type p53 and the R337His mutant at several different pH and salt conditions. We find that the 337-352 salt bridge is joined by two other charged side chains, R333 and E349. These four residues do not settle into a fixed pattern of salt bridging, but continue to exchange salt-bridging partners on the nanosecond time scale throughout the simulation. This unusual system of fluid salt bridging may explain the previous finding from alanine scanning experiments that R333 contributes significantly to protein stability, even though in the crystal structure it is extended outward into solvent. This extended conformation of R333 appears to be the result of a specific crystal contact and, this contact being absent in the simulation, R333 turns inward to join its interaction partners. When R337 is mutated to His but remains positively charged, it maintains the original interaction with D352, but the newly observed interaction with E349 is weakened, accounting for the reduced stability of R337H even under mildly acidic conditions. When this His is deprotonated, the interaction with D352 is also lost, accounting for the further destabilization observed under mildly alkaline conditions. Simulations were carried out using both explicit and implicit solvent models, and both displayed similar behavior of the fluid salt-bridging cluster, suggesting that implicit solvent models can capture at least the qualitative features of this phenomenon as well as explicit solvent. Simulations under strongly acidic conditions in implicit solvent displayed the beginnings of the unfolding process, a destabilization of the hydrophobic dimer-dimer interface. Computational alanine scanning using the molecular mechanics Poisson-Boltzmann surface area method showed significant correlation to experimental unfolding data for charged and polar residues, but much weaker correlation for hydrophobic residues.