Charge substitutions generated by site-directed mutagenesis at the termini of adjacent anti-parallel a-helices in a four-helix bundle protein were used to determine a precise value for the contribution of indirect charge-charge interactions to overall protein stability, and to simulate the electrostatic effects of a-helix macrodipoles. Thermodynamic double mutant cycles were constructed to measure the interaction energy between such charges on adjacent anti-parallel helices in the four-helix bundle cytochrome b562 from Escherichia coli. Previously, theoretical calculations of helix macrodipole interactions using modeled four-helix bundle proteins have predicted values ranging over an order of magnitude from 0.2 to 2.5 kcal/mol. Our system represents the first experimental evidence for electrostatic interactions such as those between partial charges due to helix macrodipole charges. At the positions mutated, we have measured a favorable interaction energy of 0.6 kcal/mol between opposite charges simulating an anti-parallel helix pair. Pairs of negative or positive charges simulating a parallel orientation of helices produce an unfavorable interaction of similar magnitude. The interaction energies show a strong dependence upon ionic strength, consistent with an electrostatic effect. Indirect electrostatic contacts do appear to confer a limited stabilization upon the association of anti-parallel packing of helices, favoring this orientation by as much as 1 kcal/mol at 20 mM K phosphate.Keywords: a-helix macrodipoles; cytochrome b562; electrostatic interactions; four-helix bundle; protein stabilityThe simplicity of the four-helix bundle structure motif makes it ideal for experimental and theoretical studies on protein folding and stability. In addition to its widespread use in model calculations (Sheridan et al., 1982;Chou et al., 1988;Gilson & Honig, 1989; Carlacci & Chou, 1990a,b), it is the template for at least three attempts at de novo protein design (Regan & DeGrado, 1988;Hahn et al., 1990; Hecht et al., 1990). Examples of this structure family include cytochrome c', hemerythrin, myohemerythrin, apoferritin, bovine somatotropin, tobacco mosaic virus coat protein, interleukin 2, the ColE1 Rop protein, and cytochrome b562 (Chou et al., 1988). We have chosen the four-helix bundle to study the role of electrostatic interactions such as those between a-helix macrodipoles in protein folding and stability. Abbreviations: cyt b562, cytochrome bS62; Guan-HCI, guanidine hydrochloride; AG, change in Gibbs free energy for unfolding; AAG, difference in AG values; AGin,, interaction free energy; T,, temperature of half denaturation; AS,, change in entropy at the T,,,; AH,, change in enthalpy at the T,; AT,, change in T,; wt, wild type; Keq, equilibrium constant for unfolding; PCR, polymerase chain reaction.