The development and optimization of classical molecular models is typically challenging and time-consuming. Despite notable progresses in developing efficient methods and optimization approaches, accurately predicting experimental observables and ensuring transferability across varying thermodynamics conditions remains in most cases a significant challenge. A considerable example is the case of water, for which current state-of-art models struggle in matching all the relevant cases of interest at the same time, 12 e.g., bulk properties, free energy of hydration of compounds, stabilization of lipid membranes, interaction with proteins, etc.