Integral equation model for aqueous solvatlon of polyatomic solutes: Application to the determination of the free energy surface for the internal motion of biomolecules