Semiempirical quantum calculations of biomolecular systems in solution were performed using recently developed linear-scaling methods to examine the role of solute polarization in the process of solvation. The solvation free energy of several protein and DNA molecules and complexes were computed and decomposed to asses the relative magnitude of electrostatic and polarization contributions. The effect of solvation and complex formation on the electronic density of states was also studied.