The Becke3LYP functional of DFT theory was used to investigate molecular structure and sodium affinity of the systems CH 3 CO 2 Na (1), CH 3 -O-SO 3 Na (2), CH 3 -NH-SO 3 Na (3), saccharide_1Na 2 (4), saccharide_2Na (5), saccharide_3Na 3 (6), saccharide_4Na 2 (7), and saccharide_5Na 2 (8), respectively, which are models of N-and Osulfate glycosaminoglycans. Interaction enthalpies, entropies and Gibbs energies of the sodium-coordinated systems in the gas phase were determined with the B3LYP/6-311+G(d,p) and B3LYP/6-31+G(d) methods. The computed Gibbs energies, DG o , of model systems 1-3 are negative and span a rather broad energy interval (from -500 to -1,500 kJ mol -1 ). Gibbs interaction energies for sodium acetate, sodium sulfate and sodium sulfamate functions of the five saccharides, systems 4-8 are always lower than those values found for the model compounds 1-3. The ionization of sodium salts of saccharides studied in gas phase is in most cases connected with considerable conformational rearrangement of the ionic species. This rearrangement causes an additional energetic stabilization of anionic species and is connected with the substantial release of entropy.