Prediction of secondary ionization of the phosphate group in phosphotyrosine peptides.

A computational approach, based on a continuum molecular electrostatics model, for the calculation of the pK(a) values of secondary ionization of the phosphate group in phenyl phosphate derivatives is described. The method uses the ESP atomic charges of the mono-anionic and di-anionic forms of the ionizable phosphate group, computed with the use of the density functional method, and applies a new concept of the model group, being the reference state for the pK(a) calculations. Both conformational flexibility and tautomeric degrees of freedom are taken into account in the calculations. The method was parameterized using experimentally available pK(a) values of four derivatives of phenyl phosphates, and phosphotyrosine. Subsequently this parameterization was used to predict pK(a) of the phosphate group in a short peptide Gly-Gly-Tyr(P)-Ala, and in a longer peptide consisting of 12 residues, the latter in water, and in a complex with a protein-phospholipase. The agreement between the computed and the experimental pK(a) values is better than +/-0.3 pH units for the optimized solute dielectric constant of 11-13. This approach is promising and its extension to other phospho-amino acids is in progress.

Binding of SPXK- and APXK-peptide motifs to AT-rich DNA. Experimental and theoretical studies.

The binding properties of the SPXK- and APXK-type peptides to the AT-rich DNA fragments of different length were studied by measuring the competition of peptides with Hoechst 33258 dye for DNA binding and by the gel shift assay analysis. In parallel to the experimental studies, molecular modeling techniques were used to analyze possible binding modes of the SPXZ and APXK motifs to the AT-rich DNA. The results of the competition measurements and gel shift assays suggest that serine at the i-1 position (i is proline) can be replaced by alanine without affecting the binding properties of the motif. Thus, the presence of the conserved serine in this motif in many DNA-binding proteins is probably not dictated by structural requirements. Based on the results of molecular modeling studies we propose that the binding mode of the SPXK-type motifs to the AT-rich DNA resembles closely that between the N-terminal arm of the homeodomain and DNA. This model confirms that serine in the SPXK motifs is not essential for the DNA binding. The model also indicates that if X in the motif is glutamic acid, this residue is probably protonated in the complex with DNA.