Isothermal titration calorimetry (ITC) was used to study the binding of Cd(2+) and Zn(2+) by glutathione (GSH) and phytochelatins (PC(n)), the metal sequestering compounds in plants and algae. The results are compared with those obtained by differential pulse polarography (DPP) assisted by multivariate curve resolution with alternating least squares (MCR-ALS) and by electrospray ionization mass spectrometry (ESI-MS). ITC allows one to determine (i) the stoichiometries of the different complexes (and confirms those found by DPP/MCR-ALS and ESI-MS) and (ii) their binding and thermodynamic parameters. For Cd-PC(4), the sequential binding sites model with two identical sites yields the best fitting of ITC curves and confirms the presence of CdPC(4) and Cd(2)PC(4) complexes. For Zn-PC(4), exothermic formation of ZnPC(4) is reported. Conditional stability and formation constants for Cd-GSH and Zn-GSH are determined from the fitting of the proper model to experimental ITC curves. The effect of different buffers in the complexation processes shows the key role of the choice of the buffer in calorimetric study.