Solvent polarization around a polar solute molecule plays an essential role in determining the electronic and thermodynamic properties of solutions. In this study, a solvent-polarizable model in response to solute polarization is proposed, which is coupled with a three-dimensional reference interaction-site model theory. The charge-response kernel is used to describe solvent polarizability, and four different coupling schemes are assessed. The most feasible behavior scheme among them is the one that incorporates responses not only to solute polarization but also to solute-induced solvent polarization. The numerical results indicated that solvent molecules near the polar solute show significant polarization, and therefore, the model proposed here is useful for considering the solvation process and thermodynamics of polar solute molecules.