We used a novel charge optimization technique to study the small ribonuclease barnase and to analyze its interaction with a natural tight binding inhibitor, the protein barstar. The approach uses a continuum model to explicitly determine the charge distributions that lead to the most favorable electrostatic contribution to binding when competing desolvation and interaction effects are included. Given its backbone fold, barstar is electrostatically optimized for tight binding to barnase when compared with mutants where residues have been substituted with one of the 20 common amino acids. Natural proteins thus appear to use optimization of electrostatic interactions as one strategy for achieving tight binding.