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) [CuII(aa)2] 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 CuII(aa)2 could form with l-cysteine (Cys), that is, for bis(l-cysteinato)copper(II) [Cu(Cys)2] 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)2, 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.

Conformational Analyses of Physiological Binary and Ternary Copper(II) Complexes with l-Asparagine and l-Histidine; Study of Tridentate Binding of Copper(II) in Aqueous Solution.

This study explores the structural properties and energy landscapes of the physiologically important bis(l-asparaginato)copper(II) [Cu(l-Asn)2] and (l-histidinato)(l-asparaginato)copper(II) [Cu(l-His)(l-Asn)]. The conformational analyses in the gas phase and implicitly modeled water medium, and magnetic parameters of electron paramagnetic resonance spectra were attained using density functional theory calculations. The apical CuII coordination and hydrogen bonding were analyzed. Predicted lower-energy structures enabled the confirmation and, for apical bonding, also the refinement of structural proposals from literature. Available experimental results were indecisive regarding the amido-group binding in the CuII equatorial plane in solutions, but the examination of the relative stability of Cu(l-Asn)2 conformers in 30 binding modes confirms the glycine-like mode as the most stable one. Previously reported experimental results for Cu(l-His)(l-Asn) were interpreted for l-His to have a tridentate histamine-like mode. However, the aqueous conformers with l-His in the glycinato mode are also predicted to have low energies, which does not contradict the tridentate l-His binding. The predicted magnetic parameters of conformers with an apical oxygen atom (intramolecular or from a water molecule) can reproduce the experimental data. An extent of conformational flexibility and abundance of l-His-containing ternary copper(II) amino acid complexes under physiological conditions may be related.

Calculating the geometry and Raman spectrum of physiological bis(L-histidinato)copper(II): an assessment of DFT functionals for aqueous and isolated systems.

Reliable density functional theory (DFT) calculations can be performed in conjuction with spectroscopic measurements to elucidate the structural properties of physiologically important bis(amino acidato)copper(II) compounds in solutions. They can provide insight into the influence of intermolecular interactions on the molecular geometry in the crystal lattice or solution when compared with a DFT gas-phase minimum. Our previous paper [Markovic et al. (2014) Eur J Inorg Chem 198] reported the DFT-determined geometries and Raman spectra for different conformers of physiological bis(L-histidinato)copper(II) with 20 explicit water molecules, as calculated using the B3LYP functional. The present study examined the reliability of those B3LYP results by applying the M06 functional instead, as it should better account for noncovalent interactions. The water molecules were positioned more compactly around the complex by M06 than by B3LYP. The accuracies of the two functionals when compared to relevant experimental data showed that M06 was better at reproducing in-plane Cu-N bond lengths but B3LYP gave more accurate axial Cu-O distances. Both functionals reproduced the experimental Raman spectrum at pH 8 to similar levels of accuracy and provided precise information on the Cu(II) coordination mode and conformation in aqueous solution. Additionally, we assessed several DFT and DFT-D functionals (BP86, B3LYP, B3LYP-D, M06, M06 L, wB97XD, mPW2PLYPD) by using them to model the geometries of experimental bis(L-histidinato)copper(II) crystalline conformations as isolated systems, and then benchmarking the results against those from high-level second-order pertubation Møller-Plesset (MP2) calculations. Although this assessment resulted in an equivocal conclusion because the MP2 results for the isolated complex were inconsistent with the corresponding DFT outcomes, it does provide new information on future benchmark options.

The Important Role of the Hydroxyl Group on the Conformational Adaptability in Bis(l-threoninato)copper(II) Compared to Bis(l-allo-threoninato)copper(II): Quantum Chemical Study.

Detailed structural properties of physiological bis(amino acidato)copper(II) complexes are generally unknown in solutions. This paper examines how stereochemical differences between the essential amino acid l-threonine and its diastereomer l-allo-threonine, which is rarely present in nature, may affect relative stabilities of bis(l-threoninato)copper(II) and bis(l-allo-threoninato)copper(II) in the gas phase and aqueous solution. These amino acids can bind to Cu(II) via the nitrogen and carboxylato oxygen atoms, the nitrogen and hydroxyl oxygen atoms, and the carboxylato and hydroxyl oxygen atoms. We term these coordination modes G, No, and Oo, respectively. The density functional theory (DFT) calculations with the B3LYP functional of the conformational landscapes for all possible coordination modes of both complexes revealed their very similar stability in the gas phase and in aqueous solution. The conformational analyses resulted in 196 and 267 conformers of isolated copper(II) chelates with l-threonine and l-allo-threonine, respectively. The G-G coordination mode is the most stable, both in the gas phase and aqueous solution. Very similar energy values of the lowest-energy solvated cis and trans G-G conformers in implicitly accounted water medium are in accord with the experimental results that these isomers are present in aqueous solution at physiological pH values. The transition-state structures, activation Gibbs free energies, and reaction rates calculated using DFT/B3LYP and MP2 for the transformations from the most stable cis G-G and trans Oo-G conformers to trans G-G ones for the first time reveal several alternate coordination-mode transformation mechanisms in the copper(II) complexes with amino acids other than glycine. The trans Oo-G conformers are kinetically more stable than cis G-G ones in the gas phase. The only significant difference found between the two complexes is a more suitable position of the hydroxyl group in physiological bis(l-threoninato)copper(II) to form intramolecular hydrogen bonds, which may restrain its conformational space.

Solid-State NMR Study of Paramagnetic Bis(alaninato-κ(2)N,O)copper(II) and Bis(1-amino(cyclo)alkane-1-carboxylato-κ(2)N,O)copper(II) Complexes: Reflection of Stereoisomerism and Molecular Mobility in (13)C and (2)H Fast Magic Angle Spinning Spectra.

Solid-state stereochemistry and mobility of paramagnetic copper(II) complexes formed by aliphatic amino acids (l-alanine, d,l-alanine, 1-amino-2-methyl-alanine) and 1-amino(cyclo)alkane-1-carboxylic acids (alkane = propane, butane, pentane, hexane) as bidentate ligands has been studied by (13)C and (2)H solid-state fast magic angle spinning (MAS) NMR spectroscopy. We examined the prospective method to characterize solid-state paramagnetic compounds in a routine way. Both (13)C and (2)H MAS spectra can distinguish d,l and l,l diastereomers of natural and polydeuterated bis([Dn]alaninato)copper(II) (n = 0, 2, 8) complexes with axial and/or equatorial methyl positions (conformations) primarily due to different Fermi-contact (FC) contributions. The three-bond hyperfine couplings clearly show Karplus-like dependence on the torsional angles which turned out to be a useful assignment aid. Density functional theory calculations of the FC term and crystal structures were also used to aid the final assignments. The correlations obtained for bis(alaninato-κ(2)N,O)copper(II) complexes were successfully used to characterize other complexes. The usefulness of the (2)H MAS spectra of the deuterated complexes was underlined. Even the spectra of the easily exchangeable amine protons contained essential stereochemical information. In the case of a dimer structure of bis(1-aminohexane-1-carboxylato-κ(2)N,O)copper(II) both the (13)C and (2)H resolutions were good enough to confirm the presence of the cis and trans forms in the asymmetric unit. With regard to the internal solid-state motions in the crystal lattice, the obtained quadrupolar tensor parameters were similar for the d,l- and l,l-alaninato isomers and also for the cis-trans forms suggesting similar crystal packing effects, static amine deuterons involved in hydrogen bonding, and fast rotating methyl groups.

Combined experimental and computational study of cis-trans isomerism in bis(L-valinato)copper(II).

Heating of polycrystalline cis aquabis(L-valinato)copper(II) at 90 °C resulted in a dehydrated powder. Recrystallization from aqueous solution of the obtained product yielded anhydrous trans bis(L-valinato)copper(II). The X-ray crystal and molecular structures of trans bis(L-valinato)copper(II) and cis aquabis(L-valinato)copper(II) are presented. Molecular modeling calculations were attempted to resolve factors that influenced the isomerization and crystallization of either the aqua cis- or the anhydrous trans-isomer. Conformational analyses of trans- and cis-isomers were completed in vacuo and in crystal by molecular mechanics, and in aqueous solution by molecular dynamics (MD) simulations using the same force field. Although the conformers with trans-configuration are the most stable in vacuo, those with cis-configuration form more favorable intermolecular interactions. Consequently, both cis- and trans-isomers are predicted to be present in aqueous solution. According to the crystal structure simulations and predictions, cis-isomer requires water molecules to form energetically more stable crystal packings than trans-isomer. The MD modeling of the self-assembly of 16 bis(L-valinato)copper(II) complexes in aqueous solution for the first time predicted the crystallization nucleus formation to proceed from monomers to oligomers by Cu-to-O(carboxylato) and/or N-H···O(carboxylato) weak bonds; these oligomers then bind together via water molecules until they acquire the right positions for noncovalent bonding like in the experimental crystal structures. Fifty-nanosecond MD simulations accomplished for a system consisting of equal numbers of complexes and water molecules at 298 and 370 K suggested complete cis-to-trans transformation at the higher temperature. Prevalence of either cis- or trans-conformers in water upon dissolvation may explain the crystallization results.

Structure Prediction of Bis(amino acidato)copper(II) Complexes with a New Force Field for Molecular Modeling.

This article presents a new force field whose parameterization was based on experimental crystal data and quantum chemically obtained vacuum structures of a series of copper(II) complexes with aliphatic α-amino acids and their N-alkyl derivatives, along with the SPC/E water model. The ability of the new force field to reproduce and predict the structural properties of the copper(II) complexes in the gas phase, in simulated crystalline surroundings, and solvated in water is examined. Molecular dynamics (MD) simulations with the new force field yielded time-average structural coordinates of bis(glycinato)copper(II) [the only one of 25 modeled bis(amino acidato)copper(II) systems with published experimental structural data in aqueous solution at room temperature] within the experimental error values. The study of the cis-trans isomerization of bis(glycinato)copper(II) in aqueous medium at 300 K using the quantum chemical polarized continuum model revealed a small energy difference (5 kJ mol(-1)) between the solvated cis and trans minima, in line with the MD energy estimations. The new force field proved promising in predicting the association of the complexes in aqueous solution and formation of a nucleus of crystallization.

Is the enthalpy of fusion of tris(acetylacetonato)metal(III) complexes affected by their potential energy in the crystal state?

In this Article we present enthalpies of fusion and melting points obtained from new thermochemical measurements of tris(acetylacetonato)metal(III), M(acac)(3), complexes (M = Fe, Al, Cr, Mn, Co) using differential scanning calorimetry (DSC) and evaluate them in relation to their different values found in the literature. An enthalpy of fusion of 27.67 kJ mol(-)(1) was derived for Mn(acac)(3) from a symmetrical DSC thermogram captured for the first time. The enthalpy value was indirectly confirmed with the solubility measurements of Mn(acac)(3) in acetylacetone. A hypothesis has been stated that the enthalpy of fusion and the potential energy of M(acac)(3) in the crystal state may be related. To calculate molecular in-crystal potential energy, in this Article we proposed a molecular mechanics model for the M(acac)(3) class of compounds. Nine X-ray crystal structures of M(acac)(3) complexes (M = Fe, Al, V, Mn, Co, Cr, Sc) were included in the modeling. The conformational potential energy was minimized for a molecule surrounded by other molecules in the crystal lattice. The partial charges from two schemes, the electrostatic potential (ESP) fit and the natural population analysis (NPA), were used to construct two types of force fields to examine which force field type would yield a better fit with the experimental thermal properties. The final force fields were named FF-ESP and FF-NPA. Both force field sets reproduced well the experimental crystal data of nine M(acac)(3) complexes as well as of tris(3-methyl-2,4-pentanedionato-O,O')cobalt(III). Only in-crystal potential energies derived by FF-NPA yielded a significant correlation (correlation coefficient R = -0.71) with the measured enthalpies of fusion. The enthalpy of fusion for Co(acac)(3) could not be determined experimentally because of simultaneous decomposition and fusion, and it is predicted to be 33.2 kJ mol(-)(1) from the correlation regression line.

Modeling anhydrous and aqua copper(II) amino acid complexes: a new molecular mechanics force field parametrization based on quantum chemical studies and experimental crystal data.

This paper presents the vacuum structures of aquacopper(II) bis(amino acid) complexes with glycine, sarcosine, N,N-dimethylglycine, and N-tert-butyl-N-methylglycine estimated using the B3LYP method. The differences between the B3LYP vacuum structures and experimental crystal structures suggested considerable influence of crystal lattice packing effects on the changes in the complexes' geometries. A previously developed molecular mechanics force field for modeling anhydrous copper(II) amino acidates was reoptimized to simulate these changes and predict the properties of both trans and cis anhydrous and aqua copper(II) amino acid complexes. The modeling included experimental molecular and crystal structures of 13 anhydrous and 10 aqua copper(II) amino acidates with the same atom types (Cu(II), C, H, N, and O) but various copper(II) coordination polyhedron geometries, crystal symmetries, and intermolecular interactions. The empirical parameters of the selected potential energy functions were optimized on the B3LYP vacuum copper(II) coordination geometries of three anhydrous copper(II) amino acidates and on experimental crystalline internal coordinates and unit cell dimensions of six anhydrous and six aqua copper(II) amino acid complexes. The respective equilibrium structures were calculated in vacuo and in simulated crystalline environment. The efficacy of the final force field, FFW, was examined. The total root-mean-square deviations between the experimental and theoretical crystal values were 0.018 A in the bond lengths, 2.2 degrees in the valence angles, 5.5 degrees in the torsion angles, and 0.395 A in the unit cell lengths. FFW reproduced the unit cell volumes in the range from -8.1 to 9.6%. The means of Cu to axial water oxygen distances were 2.4 +/- 0.1 A (experiment) and 2.6 +/- 0.1 A (FFW). This paper describes the ability of the molecular mechanics model and FFW force field to simulate the flexibility of the metal coordination polyhedron. The new force field proved effective in predicting the most stable molecular conformation of copper(II) amino acidato systems in vacuo.

The influence of amino acid side chains on water binding to the copper(II) in bis(N,N-dimethyl-L-alpha-isoleucinato)-copper(II): an EPR and molecular mechanics study.

Simulations were done of the electron paramagnetic resonance (EPR) spectra for bis(N,N-dimethyl-L-alpha-isoleucinato)copper(II) dissolved in deuterated methanol as a function of temperature. They indicated different behaviour of the complex below and above 300 degrees K. The effect was examined by the conformational analysis of the copper(II) complex with a new molecular mechanics force field.

Why Are Copper(II) Amino Acid Complexes Not Planar in Their Crystal Structures? An ab Initio and Molecular Mechanics Study.

This paper presents geometries of copper(II) chelates with L-alanine, L-leucine, and L-N,N-dimethylvaline optimized by the hybrid density functional method B3LYP. According to the molecular quantum mechanics results, a square-planar copper(II) coordination geometry is electronically favored in vacuo. Deviations from the planar configuration observed in the crystal state should be attributed to sterical intramolecular and/or intermolecular effects. This paper proposes a new molecular mechanics model for tetracoordinated copper(II) amino acidates to investigate these effects in detail. The empirical parameter set for the selected potential energy functions was optimized both with respect to the X-ray crystal structures (internal coordinates and unit cell constants) and with respect to the quantum mechanically derived valence angles around copper. To test this newly developed force field (FF), the equilibrium geometries of 10 molecules are predicted in vacuo and in approximate crystalline surrounding. The results were compared with their ab initio and experimental crystal structures, respectively. The unit cell volumes were reproduced in a range from -7.0% to 2.1%. The total root-mean-square deviations between the experimental and FF in crystal internal coordinates were 0.017 Å in the bond lengths, 2.2 degrees in the valence angles, and 3.6 degrees in the torsion angles. The force field is capable of reproducing the changes in the chelate rings' torsion angles caused by the crystal packing forces and successfully explains the nonplanarity of Cu(II) amino acid complexes in their crystal structures.