A version of configuration bias Monte Carlo suitable for the insertion and deletion of ions in grand canonical Monte Carlo simulations of ionic solutions has been developed. The method has been optimized and applied to a system where the straightforward insertion and deletion of ions is prohibitively impractical. The present approach is likely to be useful in other systems as well.