New Computational Models for Electrostatics of Macromolecules in Solvents