各位老师好，最近我在做一个蛋白质-小分子相互作用相关的模拟，文章投出去后审稿人有一条回复意见说：For instance, it has already been shown by Trout and co-workers that using standard forcefields for arginine does not reproduce the correct protein-arginine preferential interactions.意在指出我是用的精氨酸力场参数需要单独被优化。我查阅了审稿人提到的这篇文章，文章里是通过实验测量的优先相互作用系数和计算得到的优先相互作用系数进行比较，并以实验为基础对小分子中各原子LJ势的σ和ε进行同比例调整，具体操作可见下图。
但是审稿人使用的是CHARMM36m下未经修改的精氨酸拓扑文件，我在文章里使用的方法是通过M062x/def2TZVP计算RESP2电荷，力常数则是先计算的Hessian矩阵，用sobtop根据DRIH方法生成。请问这种方法还需要像审稿人说的那样再去单独针对每个蛋白质来拟合小分子的力常数吗？