Our optimization methodology achieved remarkable success in reproducing potential energy surfaces that closely match high-level DFT calculations. The effectiveness of our approach is particularly evident in the accurate reproduction of reaction barriers for both systems. Starting from parameters optimized for CCR, we were able to develop accurate parameters for DHFR by incorporating a modest set of additional training data, demonstrating an efficient strategy for extending parameters to mechanistically similar reactions.