Hydration effects on the electrostatic potential around tuftsin.

The electrostatic potential and component dielectric constants from molecular dynamics (MD) trajectories of tuftsin, a tetrapeptide with the amino acid sequence Thr-Lys-Pro-Arg in water and in saline solution are presented. The results obtained from the analysis of the MD trajectories for the total electrostatic potential at points on a grid using the Ewald technique are compared with the solution to the Poisson-Boltzmann (PB) equation. The latter was solved using several sets of dielectric constant parameters. The effects of structural averaging on the PB results were also considered. Solute conformational mobility in simulations gives rise to an electrostatic potential map around the solute dominated by the solute monopole (or lowest order multipole). The detailed spatial variation of the electrostatic potential on the molecular surface brought about by the compounded effects of the distribution of water and ions close to the peptide, solvent mobility, and solute conformational mobility are not qualitatively reproducible from a reparametrization of the input solute and solvent dielectric constants to the PB equation for a single structure or for structurally averaged PB calculations. Nevertheless, by fitting the PB to the MD electrostatic potential surfaces with the dielectric constants as fitting parameters, we found that the values that give the best fit are the values calculated from the MD trajectories. Implications of using such field calculations on the design of tuftsin peptide analogues are discussed.

Salt effects on peptide conformers: a dielectric study of tuftsin.

Four 1-ns molecular dynamics computer simulations of tuftsin, Thr-Lys-Pro-Arg, are analyzed: (1) cis tuftsin in water, (2) trans tuftsin in water, (3) cis tuftsin in 1 M NaCl, and (4) trans tuftsin in 1 M NaCl. Independently of the salt concentration, the trans conformer has a higher dielectric constant than the cis conformer because the former exhibits a more widely distributed charge distribution in space. Independently of the peptide conformation, the presence of salt reduces the dielectric constants of both the peptide and the solvating water molecules because ions, on binding, restrict the motion of other atoms. In contrast to the dielectric constants, neither the peptide conformation nor the salt concentration shows a significant influence on the dielectric relaxation time of water molecules.

Simulations of conformers of tuftsin and a cyclic tuftsin analog.

The conformational properties of the configurational isomers of tuftsin, a linear tetrapeptide with the sequence Thr-Lys-Pro-Arg, were investigated with six 1 ns molecular dynamics simulations in explicit water and in a 1.0 m NaCl solution. The average conformation of the cis isomer is a type VI beta-turn. Our results indicate that water-peptide hydrogen bonding, in addition to intramolecular hydrogen bonds, stabilizes the cis conformer. The trans isomer is neither a beta- nor a gamma-turn. Results are compared with parallel studies on a cyclic analog of tuftsin, cyclo(Thr-Lys-Pro-Arg-Gly). The addition of salt does, not influence the backbone conformation of the peptide. Differences between the structures are confined to the side-chain orientations of the Lys and Arg residues.