Electrostatic properties of protein molecules have essential roles in their specific functions and their structural stabilization. However, there are not any effective methods to detect the local field directly. Moreover, since the protein-solvent system is quite heterogeneous from a dielectric point of view, it has been difficult to analyze quantitatively by any theoretical approaches. In the present article, these problems are reviewed and the recent developments are introduced.