Isothermal titration calorimetry (ITC) is a technique to measure the stoichiometry and thermodynamics from binding experiments. Identifying an appropriate mathematical model to evaluate titration curves of receptors with multiple sites is challenging, particularly when the stoichiometry or binding mechanism is not available. In a recent theoretical study, we presented a differential binding model (DBM) to study calorimetry titrations independently of the interaction among the binding sites (Herrera, I.; Winnik, M. A. J. Phys. Chem. B 2013, 117, 8659-8672): Here, we build upon our DBM and show its practical application to evaluate calorimetry titrations of receptors with multiple sites independently of the titration direction. Specifically, we present a set of ordinary differential equations (ODEs) with the general form d[S]/dV that can be integrated numerically to calculate the equilibrium concentrations of free and bound species S at every injection step and, subsequently, to evaluate the volume-normalized heat signal (delta Q(V) = delta q/dV) of direct and reverse calorimetry titrations. Additionally, we identify factors that influence the shape of the titration curve and can be used to optimize the initial concentrations of titrant and analyte. We demonstrate the flexibility of our updated DBM by applying these differentials and a global regression analysis to direct and reverse calorimetric titrations of gadolinium ions with multidentate ligands of increasing denticity) namely, diglycolic acid (DGA), citric acid (CIT), and nitrilotriacetic acid (NTA), and use statistical tests to validate the stoichiometries for the metal-ligand pairs studied.