These assignment methods allow facile assignment of polarizabilities and compatible partial charges to a wide range of compounds. In addition, we make consistent use of the direct polarization approximation in order to maximize computational efficiency during MD simulations. The direct approximation also avoids the possibility of "polarization catastrophe" inherent to the SCF approach, a numerical instability where two nearby inducible dipoles mutually polarize each other without limit. The current proof-of-concept models are trained for carbon, hydrogen, oxygen, and nitrogen, but our workflows can be easily applied to the full range of elements in drug-like molecules. We demonstrate the utility and accuracy of these models in calculations on isolated molecules and in simulations of organic liquids. We close with a discussion of next steps to further evaluate these methods and to integrate them into a comprehensive polarizable force field.