Molecular dynamics simulations of a highly charged peptide from an SH3 domain: possible sequence-function relationship

A seven-residue peptide that is highly conserved in SH3 domains despite being far from the active site has been shown by NMR to be stable in solution. This peptide, biologically important because it is a likely folding nucleus for SH3 domains, provides a challenging subject for molecular dynamics because it is highly charged. We present stable, 10-ns simulations of both the native-like diverging turn structure and a helical model. Free energies of these two conformations, estimated through MM-PBSA analysis using several force fields, suggest a comparable free energy (DeltaDeltaG < or =6 kcal/mol) for native and helix conformations. NOE intensities calculated from the native trajectory reproduce experimental data quite well, suggesting that the conformations sampled by the trajectory reasonably represent those observed in the NMR experiment. The molecular dynamics results, as well as sequence analysis of a diverse 690-member family of SH3 domain proteins, suggest that the presence of two elements is essential for formation of the diverging turn structure: a pair of residues with low helical propensity in the turn region and, as previously recognized, two hydrophobic residues to close the end of the diverging turn. Thus, these two sequence features may form the structural basis for the function of this peptide as a folding nucleus in this family of proteins.