Interaction Modes and Absolute Affinities of α-Amino Acids for Mn2+: A Comprehensive Picture

Interaction Structure and Affinity of Zwitterionic Amino Acids with Important Metal Cations (Cd2+, Cu2+, Fe3+, Hg2+, Mn2+, Ni2+ and Zn2+) in Aqueous Solution: A Theoretical Study

Heavy metals are non-biodegradable and carcinogenic pollutants with great bio-accumulation potential. Their ubiquitous occurrence in water and soils has caused serious environmental concerns. Effective strategies that can eliminate the heavy metal pollution are urgently needed. Here the adsorption potential of seven heavy metal cations (Cd2+, Cu2+, Fe3+, Hg2+, Mn2+, Ni2+ and Zn2+) with 20 amino acids was systematically investigated with Density Functional Theory method. The binding energies calculated at B3LYP-D3/def2TZVP level showed that the contribution order of amino acid side chains to the binding affinity was carboxyl > benzene ring > hydroxyl > sulfhydryl > amino group. The affinity order was inversely proportional to the radius and charge transfer of heavy metal cations, approximately following the order of: Ni2+ > Fe3+ > Cu2+ > Hg2+ > Zn2+ > Cd2+ > Mn2+. Compared to the gas-phase in other researches, the water environment has a significant influence on structures and binding energies of the heavy metal and amino acid binary complexes. Collectively, the present results will provide a basis for the design of a chelating agent (e.g., adding carboxyl or a benzene ring) to effectively remove heavy metals from the environment.

The mechanisms of α-H and proton transfers of glycine induced by Mg2+

A MP2/6-31++G(d,p)//B3LYP/6-31++G(d,p) method was used to investigate the mechanisms of α- H and proton transfers of glycine induced by Mg 2+. Eight complexes were obtained, six of which were neutral and the other two were zwitterionic. Among them, the zwitterion with a binding energy of 159.4 kcal/mol was the most stable structure. Conformation transformations of the complexes caused by the rotation of single bond and the transfers of α- H and proton were completed via seven transition states. The inductive effect of Mg 2+ made the electron cloud of glycine deviate to Mg 2+, which activated the covalent bond involving the transferred proton. The neutral complex can be turned into the zwitterionic one by the transfers of both carboxyl hydrogen and α- H , and the energy barrier of each reaction was less than 9.2 kcal/mol. After the transfer of α- H , a delocalized π bond was formed in glycine skeleton and the α- C atom took 0.19 positive charges. So the chemical activity of the glycine enhanced, and glycine was readily available for addition and nucleophilic substitution reactions. The path from the most stable glycine conformer G1 to the zwitterionic conformation I is G1 → G1–G3 → G3 → G3–G4 → G4 → G2–G4 → G2 → VI → I–VI → I, and the highest energy barrier of this path is 9.2 kcal/mol.

Structural and computational insights into the versatility of cadmium binding to proteins

Cadmium is a highly toxic group XII metal, similar to zinc and mercury. Unlike zinc, which is one of the most common metal cofactors in biology, cadmium is highly toxic. Many Zn(2+)-binding proteins can bind Cd(2+)-ions without significantly affecting their structures. Here, the protein data bank is analysed with regard to protein-cadmium interactions, which shows that cadmium can bind to a variety of ion binding sites in proteins. Statistical analysis of Cd(2+)-side chain interactions is compared with a similar analysis of other ions. This analysis reveals that with regard to amino acid side-chain preference, Cd(2+) is more similar to Mn(2+) than to Zn(2+) or Hg(2+). Finally, the interaction energies of three native metal binding proteins are calculated where Cd(2+) binds instead of Zn(2+), Ca(2+) or Cu(2+). The interaction energies are decomposed into individual components whose contributions are discussed.

Performance of M06, M06-2X, and M06-HF density functionals for conformationally flexible anionic clusters: M06 functionals perform better than B3LYP for a model system with dispersion and ionic hydrogen-bonding interactions

We present a comparative assessment of the performance of the M06 suite of density functionals (M06, M06-2X, and M06-HF) against an MP2 benchmark for calculating the relative energies and geometric structures of the Cl(-)·arginine and Br(-)·arginine halide ion-amino acid clusters. Additional results are presented for the popular B3LYP density functional. The Cl(-)·arginine and Br(-)·arginine complexes are important prototypes for the phenomenon of anion-induced zwitterion formation. Results are presented for the canonical (noncharge separated) and zwitterionic (charge separated) tautomers of the clusters, as well as the numerous conformational isomers of the clusters. We find that all of the M06 functions perform well in terms of predicting the general trends in the conformer relative energies and identifying the global minimum conformer. This is in contrast to the B3LYP functional, which performed significantly less well for the canonical tautomers of the clusters where dispersion interactions contribute more significantly to the conformer energetics. We find that the M06 functional gave the lowest mean unsigned error for the relative energies of the canonical conformers (2.10 and 2.36 kJ/mol for Br(-)·arginine and Cl(-)·arginine), while M06-2X gave the lowest mean unsigned error for the zwitterionic conformers (0.85 and 1.23 kJ/mol for Br(-)·arginine and Cl(-)·arginine), thus providing insight into the types of physical systems where each of these functionals should perform best.

Gas-phase salt bridge interactions between glutamic acid and arginine

The gas-phase side chain-side chain (SC-SC) interaction and possible proton transfer between glutamic acid (Glu) and arginine (Arg) residues are studied under low-temperature conditions in an overall neutral peptide. Conformation-specific IR spectra, obtained with the free electron laser FELIX, in combination with density functional theory (DFT) calculations, provide insight into the occurrence of intramolecular proton transfer and detailed information on the conformational preferences of the peptides Z-Glu-Alan-Arg-NHMe (n = 0,1,3). Low-energy structures are obtained using molecular dynamics simulations via the simulated annealing approach, resulting in three types of SC-SC interactions, in particular two types of pair-wise interactions and one bifurcated interaction. These low-energy structures are optimized and frequency calculations are performed using the B3LYP functional, for structural analysis, and the M05-2x functional, for relative energies, employing the 6-311+G(d,p) basis set. Comparison of experimental and computed spectra suggests that only a single conformation was present for each of the three peptides. Despite the increasing spacing between the Glu and Arg residues, the peptides have several types of interactions in common, in particular specific SC-SC and dispersion interactions between the Arg side chain and the phenyl ring of the Z-cap. Comparison with previous experiments on Ac-Glu-Ala-Phe-Ala-Arg-NHMe as well as molecular dynamics simulations further suggest that the pairwise interaction observed here is indeed energetically most favorable for short peptide sequences.