A linear-scaling self-consistent field method for calculation of the electronic structure of biological macromolecules in solution is presented. The method is applied at the semiempirical Hartree-Fock level to the determination of heats of formation, solvation free energies, and density of electronic states for several protein and DNA systems.