To demonstrate the utility of BespokeFit in deriving accurate, bespoke torsion parameters at scale, we have chosen to parametrize the entire dataset of ligands compiled by Wang et al. . The dataset comprises 199 drug-like molecules taken from 8 congeneric series, with diverse chemical moieties and a range of net charges. The dataset is often used to validate FF accuracy in the context of free energy calculations , and a possible avenue to accuracy improvements is through the bespoke parametrization of the torsion parameters. Here, we will show that the accuracy of the MM potential energy surface about each torsion angle can be substantially improved via bespoke parametrization, compared to the base FF.