The role of electrostatic interactions in determining the native fold of proteins has been investigated by analyzing the alignment of peptide bond dipole moments with the local electrostatic field generated by the rest of the molecule with and without solvent effects. This alignment was calculated for a set of 112 native proteins by using charges from a gas phase potential. Most of the peptide dipoles in this set of proteins are on average aligned with the electrostatic field. The dipole moments associated with ␣-helical conformations show the best alignment with the electrostatic field, followed by residues in -strand conformations. The dipole moments associated with other secondary structure elements are on average better aligned than in randomly generated conformations. The alignment of a dipole with the local electrostatic field depends on both the topology of the native fold and the charge distribution assumed for all of the residues. The influences of (i) solvent effects, (ii) different sets of charges, and (iii) the charge distribution assumed for the whole molecule were examined with a subset of 22 proteins each of which contains <30 ionizable groups. The results show that alternative charge distribution models lead to significant differences among the associated electrostatic fields, whereas the electrostatic field is less sensitive to the particular set of the adopted charges themselves (empirical conformational energy program for peptides or parameters for solvation energy).charge distribution ͉ dipole moment ͉ electrostatics ͉ protein conformation C harge-charge and ionizable group-peptide dipole interactions have been proposed to play a relevant role in the stabilization of proteins (1) and in determining the native fold (2). Although a quantitative analysis of the relative stability of the native fold requires consideration of a delicate balance among all physical forces, we concentrate here on how electrostatic interactions can be identified as an important determination of the backbone conformation of the native structure.It is well known that the CO and NH dipoles in an ␣-helix are electrostatically aligned nearly parallel to the axis of the helix (3), and that the alignment of peptide dipoles in various types of secondary structure elements is significant in determining the 3D structure (4) and stability of globular proteins (5, 6). The concept of the preferential alignment of the backbone dipole moments ( ) with the local electrostatic field (E) (see Fig. 1) was used by Piela and Scheraga (7) and by Ripoll and Scheraga (8) to develop procedures to surmount the multiple-minima problem in the conformational analysis of peptides and proteins. The charges used in these procedures to define , shown in Fig. 1, are ECEPP͞3 (empirical conformational energy program for peptides) charges (9) that differ from other charges, e.g., PARSE (parameters for solvation energy) charges (10).Approximately 50% of the residues in proteins fit in the most common secondary structure classes (11), such as ␣-helix or...