Ionic interactions of cytochrome c play an important role in the electron transfer process. Molecular dynamics simulations of the binding of borate ion, which serves as a model ion, at three different cytochrome C surface sites are performed. This work is motivated by previous NMR studies of cytochrome c in borate solution, which indicate the existence of two types of binding sites, a slow exchange site and a fast exchange site. These two types of binding behavior were observed in the dynamic simulations, offering a molecular interpretation of 'loose' and 'tight' binding. At the 'loose' binding sites (near Lys25/Lys27 and Lys55/Lys73) the ion forms two to three hydrogen bonds to the nearest lysine residue. This binding is transient on the time scale of the simulation, demonstrating the feasibility of fast exchange. At the 'tight' binding site(near Lys13/Lys86), on the other hand, the ion becomes integrated into the protein hydrogen bond network and remains there for the duration of the simulation (exemplifying slow exchange). Binding simulations of the ion at the 'tight' site of H26Q mutant cytochrome c also showed integration of the ion into the protein's hydrogen bond network. However, this integration differs in details from the binding of the ion to the native protein, in agreement with previous NMR observations.