The binding site dependence of binding energy in both metalated and protonated diglycine and triglycine peptides

Abstract

The dependence of binding energy (BE) on the binding sites of both protonated and various metalated (M +/2+ = Li +, Na + or K +, and Be 2+, Mg 2+ or Ca 2+) oligoglycine derivatives, G n H +M x +/2+ (G = glycine and H + = proton; n = 2, 3 and x = 1–3) in the gas phase has been determined at the B3LYP level with different basis sets. Results show that the BEs of these G n H +⋯M x +/2+ complexes will change into negative values when n ⩾ 4 for G n H +Li x +, n ⩾ 5 for G n H +Na x +, n ⩾ 6 for G n H +K x +, n ⩾ 2 for G n H +Mg x 2+, and n ⩾ 3 for G n H +Ca x 2+, respectively. The BEs of both G 2H +Be x 2+ and G 3H +Be x 2+, however, are always negative values. The signs “ n” and “ x” denote the number of glycine residues of linking the two cations (H + and M +/2+) and the serial number of site bound by a metal-ion, respectively. Moreover, the BEs (Δ E) decrease gradually along the increase of oligoglycine size ( n) between two cations for all these different G n H +M x +/2+ systems, i.e., Δ E(GH +M x +/2+) > Δ E(G 2H +M x +/2+) > Δ E(G 3H +M x +/2+). Monohydration at the metal-ion almost keeps these BEs unchangeable. The continuum solvent effect, however, can change them into negative values. Interestingly, the different potential well depths (or activation energy) on the potential energy surface give different information for these monovalent and divalent metal-ion involved systems. In detail, electrostatic effect of monovalent metal-ion on the proton transfer of amino-terminus decreases stepwise along with the number increase of the glycine residues between two cations. Moreover, the well will disappear in the process of proton transfer assisted by a water molecule. Differently, the (N–)H of neighboring the metal-bound site is the most potential proton (hydrogen) to transfer in these divalent metal-ion involved systems.