Contribution of arginine-glutamate salt bridges to helix stability

Peptide side chain interactions were studied by molecular dynamics simulation using explicit solvent on a peptide with the sequence AAARAAAAEAAEAAAARA. Three different protonation states of the glutamic acid side chains were simulated for four 20 ns runs each, a total simulation time of 240 ns. Two different salt bridge geometries were observed and the preferred geometry was found to depend on Glu - Arg residue spacing. Stable charge clusters were also observed, particularly in the fully charged peptide. Salt bridges were selectively interrupted upon protonation, with concomitant changes in secondary structure. The fully charged peptide was highly helical between residues 9 and 13, although protonation increased helicity near the N-terminus. The contribution of salt bridges to helix stability therefore depends on both position and relative position of charged residues within a sequence.