A previous NMR investigation of model decapeptides with identical β-strand sequences and different turn sequences demonstrated that, in these peptide systems, the turn residues played a more predominant role in defining the type of β-hairpin adopted than cross-strand side-chain interactions. This result needed to be tested in longer β-hairpin forming peptides, containing more potentially stabilizing cross-strand hydrogen bonds and side-chain interactions that might counterbalance the influence of the turn sequence. In that direction, we report here on the design and 1H NMR conformational study of three β-hairpin forming pentadecapeptides. The design consists of adding two and three residues at the N- and C-termini, respectively, of the previously studied decapeptides. One of the designed pentadecapeptides includes a potentially stabilizing R-E salt bridge to investigate the influence of this interaction on β-hairpin stability. We suggest that this peptide self-associates by forming intermolecular salt bridges. The other two pentadecapeptides behave as monomers. A conformational analysis of their 1H NMR spectra reveals that they adopt different types of β-hairpin structure despite having identical strand sequences. Hence, the β-turn sequence drives β-hairpin formation in the investigated pentadecapeptides that adopt β-hairpins that are longer than the average protein β-hairpins. These results reinforce our previous suggestion concerning the key role played by the turn sequence in directing the kind of β-hairpin formed by designed peptides.

The GCN4 leucine zipper is a peptide homodimer that has been the subject of a number of experimental and theoretical investigations into the determinants of affinity and specificity. Here, we utilize this model system to investigate electrostatic effects in protein binding using continuum calculations. A particularly novel feature of the computations made here is that they provide an interaction-by-interaction breakdown of the electrostatic contributions to the free energy of docking that includes changes in the interaction of each functional group with solvent and changes in interactions between all pairs of functional groups on binding. The results show that (1) electrostatic effects disfavor binding by roughly 15 kcal/mol due to desolvation effects that are incompletely compensated in the bound state, (2) while no groups strongly stabilize binding, the groups that are most destabilizing are charged and polar side chains at the interface that have been implicated in determining binding specificity, and (3) attractive intramolecular interactions (e.g., backbone hydrogen bonds) that are enhanced on binding due to reduced solvent screening in the bound state contribute significantly to affinity and are likely to be a general effect in other complexes. A comparison is made between the results obtained in an electrostatic analysis carried out calculationally and simulated results corresponding to idealized data from a scanning mutagenesis experiment. It is shown that scanning experiments provide incomplete information on interactions and, if overinterpreted, tend to overestimate the energetic effect of individual side chains that make attractive interactions. Finally, a comparison is made between the results available from a continuum electrostatic model and from a simpler surface-area dependent solvation model. In this case, although the simpler model neglects certain interactions, on average it performs rather well.