Hydrogen bonds of glycine complexes were calculated using quantum chemistry calculations at M06L-GD3/def2-TZVPP level and by analyzing the crystal structures from the Cambridge Structural Database (CSD). One hydrogen bond where amino acid plays the role of the H-donor (NH/O), and two where it plays the role of the H-acceptor (O1/HO, O1 is a coordinated oxygen atom, and, O2/HO, O2 is a non-coordinated oxygen atom) were investigated. The calculations were done on octahedral nickel(II), square pyramidal copper(II), square planar copper(II), palladium(II), and platinum(II) glycine complexes with different charges adjusted using water(s) and/or chlorine ion(s) as the remaining ligands. For NH/O hydrogen bond, interaction energies of neutral complexes are the weakest, from −5.2 to −7.2 kcal/mol for neutral, stronger for singly positive, from −8.3 to −12.1 kcal/mol, and the strongest for doubly positive complex, −16.9 kcal/mol. For O1/HO and O2/HO interactions, neutral complexes have weaker interaction energies (from −2.2 to −5.1 kcal/mol for O1/HO, and from −3.7 to −5.0 kcal/mol for O2/HO), than for singly negative complexes (from −6.9 to −8.2 kcal/mol for O1/HO, and from −8.0 to −9.0 kcal/mol for O2/HO). Additionally to the complex charge, metal oxidation number, coordination number, and metal atomic number also influence the hydrogen bond strength, however, the influence is smaller.