In this work, we developed and used differentiable molecular dynamics to automate and standardize the determination of force field parameters for HhPF modeling. Our procedure requires only the definition of a loss function, targeting specific molecular properties of interest, which may be based on any microscopic or macroscopic observable, and may target bottom-up data, i.e., from all-atom models, or top-down quantities, i.e., from experimental measurements.