Absolute Affinities of α-Amino Acids for Cu+in the Gas Phase. A Theoretical Study

Accurate and reliable thermochemistry by data analysis of complex thermochemical networks using Active Thermochemical Tables: the case of glycine thermochemistry

Active Thermochemical Tables (ATcT) were successfully used to resolve the existing inconsistencies related to the thermochemistry of glycine, based on statistically analyzing and solving a thermochemical network that includes >3350 chemical species interconnected by nearly 35 000 thermochemically-relevant determinations from experiment and high-level theory. The current ATcT results for the 298.15 K enthalpies of formation are -394.70 ± 0.55 kJ mol-1 for gas phase glycine, -528.37 ± 0.20 kJ mol-1 for solid α-glycine, -528.05 ± 0.22 kJ mol-1 for β-glycine, -528.64 ± 0.23 kJ mol-1 for γ-glycine, -514.22 ± 0.20 kJ mol-1 for aqueous undissociated glycine, and -470.09 ± 0.20 kJ mol-1 for fully dissociated aqueous glycine at infinite dilution. In addition, a new set of thermophysical properties of gas phase glycine was obtained from a fully corrected nonrigid rotor anharmonic oscillator (NRRAO) partition function, which includes all conformers. Corresponding sets of thermophysical properties of α-, β-, and γ-glycine are also presented.

Bis(amino acidato)copper(II) compounds in blood plasma: a review of computed structural properties and amino acid affinities for Cu2+ informing further pharmacological research

Neutral bis(amino acidato)copper(II) [Cu(aa)2] coordination compounds are the physiological species of copper(II) amino acid compounds in blood plasma taking the form of bis(l-histidinato)copper(II) and mixed ternary copper(II)-l-histidine complexes, preferably with l-glutamine, l-threonine, l-asparagine, and l-cysteine. These amino acids have three functional groups that can bind metal ions: the common α-amino and carboxylate groups and a side-chain polar group. In Cu(aa)2, two coordinating groups per amino acid bind to copper(II) in-plane, while the third group can bind apically, which yields many possibilities for axial and planar bonds, that is, for bidentate and tridentate binding. So far, the experimental studies of physiological Cu(aa)2 compounds in solutions have not specified their complete geometries. This paper provides a brief review of my group's research on structural properties of physiological Cu(aa)2 calculated using the density functional theory (DFT) to locate low-energy conformers that can coexist in aqueous solutions. These DFT investigations have revealed high conformational flexibility of ternary Cu(aa)2 compounds for tridentate or bidentate chelation, which may explain copper(II) exchange reactions in the plasma and inform the development of small multifunctional copper(II)-binding drugs with several possible copper(II)-binding groups. Furthermore, our prediction of metal ion affinities for Cu2+ binding with amino-acid ligands in low-energy conformers with different coordination modes of five physiological Cu(aa)2 in aqueous solution supports the findings of their abundance in human plasma obtained with chemical speciation modelling.

Structure prediction of physiological bis(amino acidato)copper(II) species in aqueous solution: The copper(II) compounds with l-glutamine and l-histidine

Neutral (l-histidinato)(l-glutaminato)copper(II) [Cu(His)(Gln)] has been established as the most abundant ternary copper(II) amino acid compound of the exchangeable copper(II) pool in blood plasma. The experimental studies of Cu(His)(Gln) and bis(glutaminato)copper(II) [Cu(Gln)2] in solutions did not specify their complete geometries. To determine the geometries, this paper investigates the conformers, energy landscapes, and a structure-magnetic parameters relation of Cu(Gln)2 and Cu(His)(Gln) by the density functional theory (DFT) calculations. We assume a glycine-like coordination of Gln (other coordination patterns are dismissed because of steric reasons), and three His in-plane copper(II) binding modes. The conformational analyses are performed in the gas phase and implicitly modeled aqueous solution. The reliability of the DFT relative electronic and Gibbs free energies of the Cu(His)(Gln) conformers is confirmed by benchmarking against the corresponding energies obtained by the domain-based local pair natural orbital coupled-cluster method with singles, doubles, and perturbative triples [DLPNO-CCSD(T)]. Several cis- and trans-Cu(His)(Gln) conformers with His in the histaminate-like and glycine-like modes have low Gibbs free energies, and the greatest estimated metal-binding affinities. The DFT-calculated magnetic parameters of the low-energy conformers reproduce best the experimental electron paramagnetic resonance parameters measured in aqueous solutions for trans- and cis-Cu(Gln)2 conformers having two oxygen atoms (either from Gln or water molecules) at the apical positions, and Cu(His)(Gln) conformers having His in the histaminate-like mode with an apically placed carboxylato oxygen atom. The predicted conformational flexibility of His‑copper(II)-amino acid compounds may be connected with their physiological abundance, and the role in copper(II) exchange reactions in blood plasma.

Structure prediction of neutral physiological copper(II) compounds with l-cysteine and l-histidine

Bis(aminoacidato)copper(II) [Cu^II(aa)₂] coordination compounds are the physiological species of copper(II) amino acid compounds in blood plasma. Since there are no experimental data in the literature about the geometries that physiological Cu^II(aa)₂ could form with l-cysteine (Cys), that is, for bis(l-cysteinato)copper(II) [Cu(Cys)₂] and the ternary (l-histidinato)(l-cysteinato)copper(II) [Cu(His)(Cys)], this paper computationally examines the possible conformations that the two compounds could form with the Cys ligand having a protonated sulfur, as in the conventional zwitterion, which was determined to be prevailing in aqueous solution. These two amino acids can bind metals in a tridentate fashion and thus form many possible coordination patterns. Density functional calculations were performed for the conformational analyses in the gas phase and in implicitly modeled aqueous solution using a polarizable continuum model. Additionally, we examine which coordination mode, with thiol or thiolate group, is more stable. The Cys coordination via the amino N and carboxylato O atoms (a glycinato mode) is obtained as the most stable one in aqueous Cu(Cys)₂, and also in Cu(His)(Cys) when the His glycinato or histaminato mode combines with the intact thiol group. Whereas the conformers with N and thiol S as the copper(II) donor atoms are predicted to be the least stable, those with the Cu-N and Cu-S(thiolate) bonding (and protonated carboxylato group) are the most stable. The differences are explained by different covalent and ionic contributions of Cu-S(thiol) vs. Cu-S(thiolate). The study can contribute to the insight into formation and reactivity of the copper(II) cysteinato complexes in solution.

Interaction of l-proline with group IIB (Zn2+, Cd2+, Hg2+) metal cations in the gas and aqueous phases: a quantum computational study

The gas and aqueous phase complexation geometries, electronic interactions, and metal ion affinities of Zn2+, Cd2+, and Hg2+ metal cations with the two most stable conformations of l-proline complexes were studied. The complexes were optimized by density functional theory (B3LYP) using the 6-311++G(d,p) orbital basis set and relativistic pseudopotentials for the metal cations. The interactions of the metal cations at different nucleophilic sites of l-proline were considered as were three modes of interactions including salt bridged, charge solvated 1, and charge solvated 2, which are indicative of binding in a bidentate manner through the carboxylate group, carbonyl and hydroxyl oxygen, and carbonyl oxygen and the nitrogen atom of l-proline. All of the coordination patterns were characterized by both charge transfer and ionic interactions between l-proline and the metal cation. The metal ion affinity (MIA) and interaction energy were also computed for all of the complexes at both the gas and aqueous phases. Results showed that the order of MIA at the gas and aqueous phases are different. MIA order at the gas phase was in the order of Zn2+ > Hg2+ > Cd2+ whereas at the aqueous phase, the order of Zn2+ > Cd2+ > Hg2+ was obtained for MIA. The infrared stretching vibrational modes of the N–H and O–H groups of free l-proline were compared with l-proline–M2+ in both CS1 and CS2 coordination patterns at the gas phase and results showed a considerable shift to lower frequency during complexation process.

Metal-interacted histidine dimer: an ETS-NOCV and XANES study

We have analyzed the metal coordination in a histidine dimer, hydrated with a water molecule, based on the extended transition state scheme with the theory of natural orbitals for chemical valence (ETS-NOCV).

Gas-phase interactions of organotin compounds with glycine

Gas-phase interactions of organotins with glycine have been studied by combining mass spectrometry experiments and quantum calculations. Positive-ion electrospray spectra show that the interaction of di- and tri-organotins with glycine results in the formation of [(R)2Sn(Gly)-H](+) and [(R)3Sn(Gly)](+) ions, respectively. Di-organotin complexes appear much more reactive than those involving tri-organotins. (MS/MS) spectra of the [(R)3Sn(Gly)](+) ions are indeed simple and only show elimination of intact glycine, generating the [(R)3Sn](+) carbocation. On the other hand, MS/MS spectra of [(R)2Sn(Gly)-H](+) complexes are characterized by numerous fragmentation processes. Six of them, associated with elimination of H2O, CO, H2O + CO and formation of [(R)2SnOH](+) (-57 u),[(R)2SnNH2](+) (-58 u) and [(R)2SnH](+) (-73 u), are systematically observed. Use of labeled glycines notably concludes that the hydrogen atoms eliminated in water and H2O + CO are labile hydrogens. A similar conclusion can be made for hydrogens of [(R2)SnOH](+) and [(R2)SnNH2](+) ions. Interestingly, formation [(R)2SnH](+) ions is characterized by a migration of one the α hydrogen of glycine onto the metallic center. Finally, several dissociation routes are observed and are characteristic of a given organic substituent. Calculations indicated that the interaction between organotins and glycine is mostly electrostatic. For [(R)2Sn(Gly)-H](+) complexes, a preferable bidentate interaction of the type η(2)-O,NH2 is observed, similar to that encountered for other metal ions. [(R)3Sn](+) ions strongly stabilize the zwitterionic form of glycine, which is practically degenerate with respect to neutral glycine. In addition, the interconversion between both forms is almost barrierless. Suitable mechanisms are proposed in order to account for the most relevant fragmentation processes.

Interactions of Zn(II) with single and multiple amino acids. Insights from density functional and ab initio calculations

Calculations were performed to study the interactions of metal ions (M) with (multiple) amino acids (AA) and fill the gap between single AA and proteins. A complete conformational search results in nine and eleven ZnGly isomers at B3P86 and MP2 levels, respectively, and four populated conformers of glycine are responsible for production of these isomers. For all M, the isomers via the OO and NO binding modes are the main constituents, and the OO mode is favored by stronger electrostatic interactions. Binding with more glycines causes larger structural distortions, improves relative stabilities of monodentate binding isomers and generates new binding modes (e.g. ZnB(III) via only the hydroxyl group). The scaling factor of Zn(Gly)(n) structures, the ratio of its binding affinity versus the sum of comprising ZnGly isomers, is linear with glycine number (n), and the linear relationship may not be altered by mutations of glycines and M. It thus allows to estimate M(AA)(n) binding affinities (n ≥ 2) from the comprising MAA structures and analyze their structures with kinetic methods. The DFT and MP2 results become comparable by increasing metal coordination, e.g. the ZnB(III) versus ZnA(I) (zwitterionic) relative energy differs by 41.9 kcal mol(-1) at B3P86 and MP2 levels and is close by addition of three water molecules (4.1 kcal mol(-1)). The presence of water solvent improves the relative stabilities of monodentate binding isomers and results in a broader conformational distribution.

Extensive computational study on coordination of transition metal cations and water molecules to glutamic acid

On the basis of the conformations of glutamic acid (Glu) and analysis of possible metal cation coordination and hydration modes, conformations of Glu metalated with transition metal cations (TMCs), Cu(+/2+), Zn(+/2+), and Fe(+/2+/3+) and hydrations of Glu-Cu(+/2+) and Glu-Zn(+/2+) complexes by up to three water molecules are determined by extensive computational searches. The BHandHLYP functional is chosen as the main computational method as its overall performance for treating the spin multiplicity of TMCs is similar to that of CCSD(T) and better than that of MP2 and B3LYP. All mono- and divalent TMCs prefer tridentate coordination to canonical Glu, while Fe(3+) favors a bidentate coordination to zwitterionic Glu. The ground state of Glu-Fe(+) is found to be a spin sextet. Metal ion affinities of Glu for the TMCs are determined, and an excellent agreement with the experiment for Cu(+) may be obtained if the entropic effect is properly accounted for. Effects of hydration on the stabilities of different Glu-Cu(+/2+)/Zn(+/2+) structures are discussed, and the hydration energies for up to three water molecules are obtained. For the global minimum to take the zwitterionic form, Glu-Zn(+) requires only monohydration, Glu-Zn(2+) needs to be trihydrated, while Glu-Cu(+/2) should be hydrated with four or more water molecules.

Predicting the coordination number within copper chaperones: Atox1 as case study

The concentration of copper ions in biological systems is tightly regulated by metallochaperone proteins which are responsible for Cu(I) delivery to designated locations in the cell. These proteins contain a unique motif (MXCXXC) that binds Cu(I) very tightly and specifically but at the same time allows efficient metal transfer to target proteins that often contain a similar copper binding motif. It was found that binding to Cu(I) is achieved through the two cysteine residues in a low coordination number of 2-3 due to possible binding of a third external ligand. Understanding copper transport requires better understanding of copper coordination. Here we therefore focused on establishing a computational method that can predict the coordination number of copper in copper chaperones. The method is shown to be successful in predicting the coordination of Cu(I) within the human copper chaperone (Atox1). Based on the results, a possible rationale for this unique Cu(I) dicoordination in Atox1 is suggested.

INTERACTION OF THE TAUTOMERIC STATES OF HISTIDINE WITH Cu AND Zn METAL IONS – A THEORETICAL STUDY

A density functional study on the interaction of Cu 2+ and Zn 2+ metal ions at the predominant positions of N τ– H and N π– H histidine tautomers has been performed. The fully optimized energy of the isolated histidine tautomers at B3LYP/6-311++G** level of theory depicts N τ– H tautomer to be much stable compared to that of the N π– H tautomer. The interaction of metal ions forms bidentate and tridentate complexes with N τ– H tautomer, while it is absent in the case of N π– H tautomer emphasizing the role of former in structural determination of liganated proteins. The Zn 2+ ion induces a barrier-free proton transfer when interacted at the carbonyl position of N τ– H histidine tautomer necessitating the tautomer to be in the zwitterionic form for complexation. The thermodynamical analysis predicts a blue shift in the NH and CO stretching vibrational frequencies and suggests N τ– H tautomer to be best suitable for Cu 2+ and Zn 2+ ion interactions. The topological and charge transfer studies in concert with frontier molecular orbital (FMO) analysis confirm the covalent interaction, validating the findings based on the geometrical data.

DFT STUDY OF BOND ENERGIES AND ATTACHMENT SITES OF SAMPLE DIVALENT CATIONS (Mg2+, Ca2+, Zn2+) TO HISTIDINE IN THE GAS PHASE

In view of better understanding interactions of amino acids and peptides with metallic cations in the isolated state, the model system histidine– M 2+ ( M 2+ = Mg 2+, Ca 2+, Zn 2+) has been studied theoretically. The computations have been performed with the help of the density functional theory (DFT) and the B3LYP functional. The extended basis set was the standard 6-311++G**. All the molecular complexes obtained by the interaction between several energetically low-lying tautomers/conformers/zwitterions of histidine and the cations on different binding sites were considered.

Our study shows that complexes of histidine with Mg 2+, Ca 2+, Zn 2+ are rather similar. In the isolated state, the most stable form corresponds to a tridentate complex in which the cation interacts with oxygen and two nitrogen atoms: one from the terminal NH 2 and one from the imidazole ring. All computations indicate that the metal ion affinity (MIA) decreases on going from Zn 2+ to Mg 2+ and Ca 2+, for the considered amino acid. This indicates that histidine prefers to bind to the transition metal cation rather than alkali earth metals. The influence of theses cations on the acidity of histidine were also considered. As expected, upon metal complexation, proton dissociation of histidine becomes much more favorable, that is, its acidity becomes much less endothermic.

Can DFT methods correctly and efficiently predict the coordination number of copper(I) complexes? A case study

The coordination number of various experimentally known Cu(I) compounds is studied using density functional theory. Various basis sets are tested, aiming to establish a reliable level for prediction of the coordination number of these and other Cu(I) complexes. It is found that most levels exhibit correct trends, namely, the bulkier ligands demonstrate larger preference for coordination of two ligands. Proper absolute values are obtained when dispersion corrections are also included in the calculations. It is concluded that the fairly small modified 6-31+G* basis set due to Pulay represents a good compromise between accuracy and efficiency, followed by Balabanov and Peterson's all-electron aug-cc-pVDZ basis set. The overall energy is decomposed into various components whose relative contribution to the overall tendency of forming a complex with a particular coordination is examined. It is shown that two opposing contributions play a major role: the interaction energy of the ligand being added and the deformation energy of the copper's coordination sphere prior to the ligand addition. The former being a stabilizing contribution, leads to higher coordination numbers while the later, a destabilizing contribution, is shown to favor lower coordination numbers.

First-principle conformational analysis of glycine residues in the alphabeta-tubulin dimer

Amino acids, especially glycine, have been extensively studied whereas their conformational behaviors in proteins are rather limited. In this work, all the polypeptides containing glycine residues are truncated from the alphabeta-tubulin dimer and refined with the partial optimization technique, where the backbone atoms of the previous and posterior residues to the glycine residues are fixed at the experimental Cartesian positions whereas the others fully relaxed. The combination of the polypeptide models and partial optimization technique is validated by twolayer ONIOM calculations, being effective to retain the local structures of proteins and meanwhile optimizing the concerned glycine residues towards energy minima. Owing to the lack of side chains, various types of hydrogen bonding interactions are detected in the glycine residues. The conformational analyses of proteins are mainly based on the dihedrals. Theta(1)( angleH(5)C(2)C(3)O(4)), Theta(2)( angleH(6)C(2)C(3)O(4)) and Theta(3)( angleN(1)C(2)C(3)O(4)) are the three key dihedrals to determine the glycine conformations, which were found to change synchronously and correlate with each other by four linear equations. It thus provided the first evidence of the correlation of intra-dihedrals for the amino acid residues in proteins. In addition, the three dihedrals of the gas-phase glycine conformations were also found to suit well with these linear equations, elucidating the rationality of using amino acids as the computational models of proteins.

Ca2+ selectivity of the sarcoplasmic reticulum Ca2+-ATPase at the enzyme-water interface and in the Ca2+ entrance channel

The sarcoplasmic reticulum (SR) Ca(2+)-ATPase, a P-type transmembrane protein, can transport Ca(2+) from the cytoplasmic to the luminal side over other cations specifically. The proposed Ca(2+) entrance channel, composed of the main-chain carbonyl oxygen and side-chain carboxyl oxygen atoms of the amino acids, opens on the enzyme surface, just above the biphospholipid layer membrane-water interface, where Trp residues are frequently found. In this work, the physicochemical nature of Ca(2+) selectivity over Mg(2+) on the surface of the SR Ca(2+)-ATPase has been investigated using the density functional theory (DFT) method. The selection process can be regarded as the first step of the specificity of the enzyme to transport Ca(2+). Subsequently, the specificity of the entrance channel to conduct Ca(2+) over other cations has also been explored. As revealed by thermodynamic analyses, either the aromatic or the aliphatic amino acid residues distributed on the surface of Ca(2+)-ATPase have a bigger affinity to Mg(2+) than to Ca(2+), resulting in a concentration decrease of free Mg(2+) in the local region. Thus, Ca(2+) can transport into the Ca(2+)-entrance channel more easily. Whereafter, for a small quantity of Mg(2+) entering this channel accompanying the Ca(2+) current, the strong electrostatic interactions between Mg(2+) and the ligands will limit the activity of this metal ion, which facilitates the weakly bonded Ca(2+) passing through the channel at a relatively high rate, as suggested by the "sticky-pore" hypothesis. Furthermore, the corresponding theoretical investigations have demonstrated that the increase of the ligand electronegativity can enhance their discrimination between these two cations effectively.