Interaction of neutral and zwitterionic glycine with Zn2+ in gas phase: ab initio and SIBFA molecular mechanics calculations

The interaction of Zn²⁺ with glycine (Gly) in the gas phase is studied by a combination of ab initio and molecular mechanics techniques. The structures and energetics of the various isomers of the Gly–Zn²⁺ complex are first established via high‐level ab initio calculations. Two low‐energy isomers are characterized: one in which the metal ion interacts with the carboxylate end of zwitterionic glycine, and another in which it chelates the amino nitrogen and the carbonyl oygen of neutral glycine. These calculations lead to the first accurate value of the gas‐phase affinity of glycine for Zn²⁺. Ab initio calculations were also used to evaluate the performance of various implementations of the SIBFA force field. To assess the extent of transferability of the distributed multipoles and polarizabilities used in the SIBFA computations, two approaches are followed. In the first, approach (a), these quantities are extracted from the ab initio Hartree–Fock wave functions of glycine or its zwitterion in its entirety, and for each individual Zn²⁺‐binding conformation. In the second, approach (b), they are assembled from the appropriate constitutive fragments, namely methylamine and formic acid for neutral glycine, and protonated methylamine and formate for the zwitterion; they undergo the appropriate vector or matrix rotation to be assembled in the conformation studied. The values of the Zn²⁺–glycine interaction energies are compared to those resulting from ab initio SCF and MP2 computations using both the all‐electron 6‐311+G(2d,2p) basis set and an effective core potential together with the valence CEP 4‐31G(2d) basis set. Approach (a) values closely reproduce the ab initio ones, both in terms of the total interaction energies and of the individual components. Approach (b) can provide a similar match to ab initio interaction energies as does approach (a), provided that the two constitutive Gly building blocks are considered as separate entities having mutual interactions that are computed simultaneously with those occurring with Zn²⁺. Thus, the supermolecule is treated as a three‐body rather than a two‐body system. These results indicate that the current implementation of the SIBFA force field should be adequate to undertake accurate studies on zinc metallopeptides. © 2000 John Wiley & Sons, Inc. J Comput Chem 21: 963–973, 2000