Recent years have seen the initial success of the variational implicit-solvent modeling of biomolecules. The dielectric boundary or the solute-solvent interface is the key quantity in such modeling. This work concerns the effective dielectric boundary force that contributes significantly to the conformation and dynamics of an underlying biomolecular system. Such a force is the "normal velocity" in the level-set numerical computation of equilibrium biomolecular structures. Precisely, the dielectric boundary force is defined as the shape derivative of the electrostatic free energy. Both the Poisson-Boltzmann free energy, and its Coulomb-field and Yukawa-field approximations are considered. Analytical formulas of the corresponding dielectric boundary force are derived. This is joint work with Hsiao-Bing Cheng, Li-Tien Cheng, and Xiaoliang Cheng, and Zhengfang Zhang.