The rational modification of protein stability is an important goal of protein design. Protein surface electrostatic interactions are not evolutionarily optimized for stability and are an attractive target for the rational redesign of proteins. We show that surface charge mutants can exert stabilizing effects in distinct and unanticipated ways, including ones that are not predicted by existing methods, even when only solvent-exposed sites are targeted. Individual mutation of three solvent-exposed lysines in the villin headpiece subdomain significantly stabilizes the protein, but the mechanism of stabilization is very different in each case. One mutation destabilizes native-state electrostatic interactions but has a larger destabilizing effect on the denatured state, a second removes the desolvation penalty paid by the charged residue, whereas the third introduces unanticipated native-state interactions but does not alter electrostatics. Our results show that even seemingly intuitive mutations can exert their effects through unforeseen and complex interactions.atomistic simulations | pH-dependent stability | protein pKa measurements | nuclear magnetic resonance | protein biophysics T he mutation of charged surface residues to enhance electrostatics is a popular approach in protein engineering (1). Target residues are often selected using estimates of the protein electrostatic potential (2). This approach relies on identifying residues that are involved in unfavorable electrostatic interactions, typically with residues of the same charge, or on identifying sites where new favorable electrostatic interactions can be introduced (3). Solvent-exposed charged residues are thought to not be involved in critical packing interactions in the native state, to not suffer large desolvation penalties upon protein folding, nor to make significant interactions in the denatured-state ensemble (DSE). Thus, residues to be targeted are typically chosen on the basis of calculations of protein native-state ensemble (NSE) electrostatics. Methods ranging from simple inspection of the protein surface, modified Tanford-Kirkwood approaches (3, 4) and Poisson-Boltzmann (PB) calculations have been used (5-7). Irrespective of the details, the general strategy is based on the assumption that surface electrostatic interactions can be modified without altering other native-state interactions and without altering DSE energetics. An attractive feature of these approaches is that they are computationally inexpensive and, unlike selectionbased methods or directed evolution, involve the generation of a limited number of mutants. Any increase in stability is generally assumed to result from modification of NSE electrostatics. We show that this approach leads to complicated and unanticipated results, even for very simple proteins.We use the villin headpiece subdomain HP36 as our model system. HP36 is a small three-helix protein that has become an extremely popular model system for experimental and computational studies of protein folding, owing to its ...