In contrast to our previous work, we have supplemented the atoms-in-molecule nonbonded parameters with molecule-specific bonded parameters derived from the QM Hessian matrix and torsional scans. In agreement with our previous study, we showed that the so-called modified Seminario method is able to reproduce QM normal mode vibrational frequencies to high accuracy (6.5% here). Closer examination of bond and angle force field parameters for widely used atom types reveals that these parameters are reasonably transferable between closely-related molecules. Such analyses of more complex molecules could be used to identify problems with standard force fields where bonded parameters may require re-fitting or the inclusion of more atom types. In addition, we have shown that for three molecules, QM relative energies of an ensemble of structures are modelled reasonably well with the QUBE FF when combined with torsional fitting. It should be noted that torsional fitting is the major computational expense in QUBE (since it requires a constrained QM optimization at each torsion angle), and methods to reduce this expense are under investigation.