This derivative can be calculated numerically by performing QM/MM calculations with varied MM point charges, and different efficient approximations can be employed. Generally, the most simple approximation (PCVA-E), in which a forward finite difference is used and only the reactant structure is considered, turns out to be sufficient for a qualitative assessment of uncertainties.