We propose an extended treatment of the charge response kernel (CRK), (∂Qa/∂Vb), which describes the response of partial charges on atomic sites to external electrostatic potential, on the basis of the density functional theory (DFT) via the coupled perturbed Kohn-Sham equations. The present CRK theory incorporates regulation procedures in the definition of partial charges to avoid unphysical large fluctuation of the CRK on "buried" sites. The CRKs of some alcohol and organic molecules, methanol, ethanol, propanol, butanol, dimethylsulfoxide (DMSO), and tetrahydrofuran (THF) were calculated, demonstrating that the new CRK model at the DFT level has greatly improved the performance of accuracy in comparison with that at the Hartree-Fock level previously proposed. The CRK model was also applied to investigate spatial nonlocality of the charge response through alkyl chain sequences. The CRK model at the DFT level enables us to construct a nonempirical strategy for polarizable molecular modeling, with practical reliability and robustness.