Restrains vs RMSD Benchmark
Here we compare different torsion fitting parameters for a series of fragments from the base series of ligands. Fragments with multiple rotatable bonds were fit at the same time to replicate a common bespoke fit use case.
Setup
All bespoke fits were run with fully expanded torsion k values and unique torsion smirks. The torsion k prior was 6.
For the Restraint case, after optimisation, the restraints were removed and a forcebalance torsionprofile single point was calculated, these values are reported in the table.
The RMSD is calculated at the QM minimum on the torsion drive.
RMSE = kcal/mol
RMSD = Angstrom
Molecule | Restaint RMSE | No Restraint RMSE | No Restraint + RMSD RMSE | Restraint RMSD | No Restraint RMSD | No Restraint + RMSD RMSD |
---|---|---|---|---|---|---|
| 0.0969 | 0.0638 | 0.0539 | 0.024 | 0.024 | 0.024 |
| 0.0818 | 0.0820 | 0.0834 | 0.008 | 0.008 | 0.008 |
| 0.0639 | 0.0598 | 0.0629 | 0.078 | 0.078 | 0.078 |
| 0.0850 | 0.0823 | 0.0844 | 0.019 | 0.019 | 0.019 |
| 0.0997 | 0.0996 | 0.1002 | 0.012 | 0.012 | 0.012 |
| 0.2052 | 0.1809 | 0.1102 | 0.163 | 0.165 | 0.163 |
| 0.1635 | 0.1911 | 0.2206 | 0.372 | 0.366 | 0.410 |
| 0.2181 | 0.1878 | 0.1622 | 0.147 | 0.144 | 0.134 |
| 0.2365 | 0.2495 | 0.1858 | 0.128 | 0.129 | 0.125 |
| 0.1770 | 0.1682 | 0.1719 | 0.167 | 0.165 | 0.166 |
| 0.2028 | 0.1889 | 0.1766 | 0.163 | 0.159 | 0.153 |
| 0.1262 | 0.1381 | 0.0949 | 0.138 | 0.138 | 0.138 |
| 0.0611 | 0.0616 | 0.0612 | 0.028 | 0.028 | 0.028 |
| 0.1655 | 0.1704 | 0.1328 | 0.209 | 0.205 | 0.160 |
| 0.1352 | 0.1013 | 0.0810 | 0.346 | 0.347 | 0.316 |
| 0.0860 | 0.0861 | 0.0859 | 0.018 | 0.018 | 0.018 |
| 0.1820 | 0.1371 | 0.1257 | 0.195 | 0.192 | 0.188 |
| 0.2939 | 0.2726 | 0.2636 | 0.163 | 0.155 | 0.158 |
| 0.2842 | 0.1816 | 0.1338 | 0.138 | 0.137 | 0.135 |
| 0.4331 | 0.3848 | 0.3835 | 0.315 | 0.300 | 0.264 |
| 0.1042 | 0.0916 | 0.0642 | 0.142 | 0.142 | 0.142 |
| 0.1537 | 0.1555 | 0.1310 | 0.171 | 0.181 | 0.171 |
| 0.1294 | 0.1328 | 0.1211 | 0.194 | 0.188 | 0.187 |
| 0.0842 | 0.0843 | 0.0868 | 0.013 | 0.013 | 0.013 |
| 0.0757 | 0.0790 | 0.0490 | 0.166 | 0.166 | 0.166 |
| 0.1718 | 0.1808 | 0.1027 | 0.195 | 0.195 | 0.183 |
| 0.4739 | 0.4807 | 0.4561 | 0.139 | 0.138 | 0.138 |
| 0.1725 | 0.1752 | 0.1096 | 0.185 | 0.186 | 0.180 |
| 0.0797 | 0.0795 | 0.0793 | 0.020 | 0.020 | 0.020 |
| 0.1820 | 0.1782 | 0.1711 | 0.170 | 0.168 | 0.159 |
| 0.1877 | 0.1883 | 0.1764 | 0.139 | 0.139 | 0.139 |
| 0.1064 | 0.1050 | 0.1029 | 0.032 | 0.032 | 0.032 |
| 0.2217 | 0.2128 | 0.1937 | 0.151 | 0.153 | 0.150 |
| 0.2433 | 0.2597 | 0.2089 | 0.138 | 0.137 | 0.141 |
| 0.0790 | 0.0789 | 0.0785 | 0.017 | 0.017 | 0.017 |
| 0.0825 | 0.0824 | 0.0820 | 0.019 | 0.019 | 0.019 |
| 0.2388 | 0.2136 | 0.2228 | 0.221 | 0.223 | 0.216 |
| 0.2702 | 0.2531 | 0.2130 | 0.260 | 0.260 | 0.262 |
| 0.2211 | 0.1167 | 0.0932 | 0.139 | 0.148 | 0.154 |
| 0.3787 | 0.3378 | 0.3377 | 0.251 | 0.244 | 0.245 |
| 0.2448 | 0.2328 | 0.2028 | 0.144 | 0.144 | 0.147 |
| 0.2740 | 0.2775 | 0.2539 | 0.136 | 0.136 | 0.136 |
| 0.0771 | 0.0767 | 0.0765 | 0.022 | 0.022 | 0.022 |
| 0.2053 | 0.1823 | 0.1863 | 0.024 | 0.024 | 0.024 |
Average | 0.1785 | 0.1669 | 0.1494 | 0.1345 | 0.1337 | 0.1309 |
STD | 0.0966 | 0.09192 | 0.0886 | 0.0931 | 0.09178 | 0.0905 |