To solve this problem, we and others have advocated making use of QM-to-MM mapping procedures to significantly reduce the number of parameters that require fitting to experimental physical properties, while retaining a small amount of empiricism to ensure accurate condensed phase properties. Though there are isolated examples in the literature, QM-to-MM mapping potentials are not widely available for general use in deriving full MM force fields for organic molecules. In QUBEKit, and its interface with ForceBalance, we provide here a means to train, test and apply these force fields, deriving parameters for all (including bonded) terms in the force field from QM. By using AIM-based charges, Lennard-Jones parameters and virtual sites, we have a consistent set of non-bonded parameters derived from a single set of atomic electron densities, available through widely-used quantum chemistry software packages. By coupling these parameter sets through QUBEKit, we can derive bonded parameters that are specific to the molecule under study, and fully consistent with the non-bonded parameter set. We estimate that for the 51 molecules in the test set, a total of around 20 k parameters are derived from QM, while only seven parameters are fit to experiment (in model 5b). Despite this minimal fitting approach, the model achieves highly competitive accuracy on a challenging test set of liquid property data.