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

The Binding of Cu2+ to Lipid Membranes Is Not Substantially Influenced by Electrostatic Screening

Gouy-Chapman theory predicts that salt screening and modulating the interfacial charge density should strongly influence the apparent dissociation constant, Kd,app, between Cu2+ and negatively charged phosphatidylserine (PS) lipids in supported lipid bilayers (SLBs). Specifically, Kd,app would be expected to increase (weaken binding) by a factor of 40 when 100 mM NaCl is introduced into the solution because of electrostatic screening between the membrane and Cu2+ cations. Surprisingly, however, fluorescence quenching measurements demonstrate that Kd,app increases by less than a factor of 2 when increasing the salt concentration in the presence of standard buffers, such as tris(hydroxymethyl)aminomethane (Tris). Moreover, increasing the negative surface charge density by a factor of 4 would be predicted to decrease (strengthen binding) Kd,app by 3 orders of magnitude. Instead, Kd,app increases slightly when 15 mol % phosphatidylglycerol (PG), a negatively charged lipid, is introduced into SLBs already containing 5 mol % PS. Such findings indicate that electrostatic double layer theory is not a useful approach for predicting the binding behavior of transition metal cations to negatively charged interfaces. The problem lies with the fact that standard buffers, such as Tris, form a wide variety of coordination complexes in bulk solution with transition metal cations like Cu2+. Typically, dozens of complexes are present simultaneously at any given pH value and the net charge on them ranges from positive to neutral to negative. Such variations in charge on the complexes result in electrostatic screening and interfacial potential effects that are substantially diminished or nonexistent. These results should generally apply to the binding behavior of first row transition metal ions, when the cations predominantly reside in complexes rather than as free ions. This includes in vivo conditions, where the concentration of free transition metal ions is often very low.

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.

Colorimetric Cu2+ Detection of (1E,2E)-1,2-Bis((1H-pyrrol-2-yl)methylene)hydrazine Using a Custom-Built Colorimeter

The compound (1E,2E)-1,2-bis((1H-pyrrol-2-yl)methylene)hydrazine (1) was investigated for its chemosensor application. The colorimetric response of 1 with various ions was investigated, and the selective optical change upon mixing with Cu²⁺ was found. The Cu²⁺ binding stoichiometry of 1 derived from Job's plot and the in silico study give us the tentative structural detail of the binding mode of 1 and Cu²⁺ being 1:1. The binding constant between 1 and Cu²⁺ from the Benesi-Hildebrand plot was 1.49 × 10⁴ M^-1. The limit of detection of 1 in Cu²⁺ detection was 0.64 μM (0.040 ppm), which is much lower than the WHO and US EPA maximum allowable Cu²⁺ level in drinking water (2 and 1.3 ppm, respectively). The custom-built colorimeter demonstrates a good linear relationship between Cu²⁺ concentration and electrical resistance (Ω) upon 1-Cu²⁺ ion binding.

Structural, theoretical and biological activity of mono and binuclear nickel(II) complexes with symmetrical and asymmetrical 4,6-diacetylresorcinol-dithiocarbazate ligands

The present work reports the synthesis and a structural study of two novel dithiocarbazate, the 4,6-diacetylresorcinol-S-benzyldithiocarbazate (H₃L¹) and the 4,6-diacetylresorcinol-bis(S-benzyldithiocarbazate) (H₄L²), and their Ni(II) complexes, [Ni(HL¹)(Py)] (1) and [Ni₂(L²)(PPh₃)₂] (2). Single crystal X-ray analyzes reveal mono and binuclear complexes and the metal centers with distorted square planar geometry. The analyses of the Hirshfeld surface and fingerprints plots revealed intermolecular contacts attributed to the H···H and C···H/H···C bonds. The Density Functional Theory (DFT), with the B3LYP functional and 6-311-G(d,p)/LanL2DZ basis sets, was employed to optimize the geometries of synthesized compounds. From the resulting geometries, the highest occupied and lowest unoccupied molecular orbital maps (HOMO-LUMO), orbital energy gap, electron localization function (ELF), electron density, natural bond orbital (NBO) analysis, and complexation of the ligands with Ni(II) were calculated supporting the experimental data. The ESI (+)-MS/MS data indicated the presence in solution of the characteristic fragmentation with the [H₃L¹]⁺ and [H₄L²]⁺ molecular ions for the ligands. The pharmacological potential of the dithiocarbazate ligands and their Ni(II) complexes were evaluated in vitro against MDA-MB-231 human breast cancer cells. A remarkable cytotoxic activity was observed, more evident for free ligands than complexes at low concentrations; however, this latter showed a better dose-response pattern, being more attractive in terms of pharmacokinetics and therapeutic window.

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.

Transition Metal Ion Interactions with Disordered Amyloid-β Peptides in the Pathogenesis of Alzheimer's Disease: Insights from Computational Chemistry Studies

Monomers and oligomers of the amyloid-β peptide aggregate to form the fibrils found in the brains of Alzheimer's disease patients. These monomers and oligomers are largely disordered and can interact with transition metal ions, affecting the mechanism and kinetics of amyloid-β aggregation. Due to the disordered nature of amyloid-β, its rapid aggregation, as well as solvent and paramagnetic effects, experimental studies face challenges in the characterization of transition metal ions bound to amyloid-β monomers and oligomers. The details of the coordination chemistry between transition metals and amyloid-β obtained from experiments remain debated. Furthermore, the impact of transition metal ion binding on the monomeric or oligomeric amyloid-β structures and dynamics are still poorly understood. Computational chemistry studies can serve as an important complement to experimental studies and can provide additional knowledge on the binding between amyloid-β and transition metal ions. Many research groups conducted first-principles calculations, ab initio molecular dynamics simulations, quantum mechanics/classical mechanics simulations, and classical molecular dynamics simulations for studying the interplay between transition metal ions and amyloid-β monomers and oligomers. This review summarizes the current understanding of transition metal interactions with amyloid-β obtained from computational chemistry studies. We also emphasize the current view of the coordination chemistry between transition metal ions and amyloid-β. This information represents an important foundation for future metal ion chelator and drug design studies aiming to combat Alzheimer's disease.

Experimental and theoretical investigations of cyclometalated ruthenium(ii) complex containing CCC-pincer and anti-inflammatory drugs as ligands: synthesis, characterization, inhibition of cyclooxygenase and in vitro cytotoxicity activities in various cancer cell lines

The cyclometalated ruthenium(ii) complex was synthesized and studied for cytotoxicity. The interaction of Ru(ii) complex with COX-2 was studied by experimental and molecular docking.

Complexes of arzanol with a Cu2+ ion: a DFT study

Arzanol (C₂₂H₂₆O₇) is a naturally occurring acylphloroglucinol largely responsible for the anti-inflammatory, anti-oxidant, antibiotic and antiviral activities of Helichrysum italicum. Like all acylphloroglucinols, the molecule contains a carboxylic substituent (-COR group); for arzanol, this is a -COCH₃ group. The molecule is further characterized by the presence of an α-pyrone ring bonded in meta to -COR through a methylene bridge, and of a prenyl chain bonded to the other meta position. The molecule can form up to three intramolecular hydrogen bonds (IHB) simultaneously, and their presence and patterns are the major stabilizing factors. This work considers complexes of representative conformers of arzanol with a Cu²⁺ ion, taking into account the different possibilities for the binding of the Cu²⁺ ion to the electron-density rich sites of the molecule and including simultaneous coordination to two geometrically suitable sites. Calculations were performed at the DFT/B3LYP level, using the 6-31+G(d,p) basis set for the C, O and H atoms and the LANL2DZ pseudopotential for the Cu²⁺ ion. Interaction energies show preference for simultaneous binding of Cu²⁺ to two sites. Simultaneous binding to the O of a phenol OH neighboring the prenyl chain and to the π bond of the prenyl chain appears to be the most favorable option, followed by simultaneous binding to the sp² O of the α-pyrone ring and the O of the phenol OH ortho to -COR on the side of the α-pyrone ring. The charge of the Cu²⁺ ion is reduced to +1 or slightly less in the complexes, which is consistent with the molecules' antioxidant (reducing) ability. Graphical abstract The copper ion prefers to attach to two sites of the arzanol molecule simultaneously. The arzanol molecule reduces the charge of the copper ion from +2 to +1 by transferring an electron to it; it becomes a radical molecular cation. The distribution of the unpaired electron in the molecule (as highlighted by the spin density maps) depends on the site/s to which the Cu²⁺ ion binds and on the molecule's conformer.

A DFT Study of Structural and Bonding Properties of Complexes Obtained from First-Row Transition Metal Chelation by 3-Alkyl-4-phenylacetylamino-4,5-dihydro-1H-1,2,4-triazol-5-one and Its Derivatives

Density functional calculations were used to explore the complexation of 3-alkyl-4-phenylacetylamino-4,5-dihydro-1h-1,2,4-triazol-5-one (ADPHT) derivatives by first-row transition metal cations. Neutral ADPHT ligand and mono deprotonated ligands have been used. Geometry optimizations have been performed in gas-phase and solution-phase (water, benzene, and N,N-dimethylformamide (DMF)) with B3LYP/Mixed I (LanL2DZ for metal atom and 6-31+G(d,p) for C, N, O, and H atoms) and with B3LYP/Mixed II (6-31G(d) for metal atom and 6-31+G(d,p) for C, N, O, and H atoms) especially in the gas-phase. Single points have also been carried out at CCSD(T) level. The B3LYP/Mixed I method was used to calculate thermodynamic energies (energies, enthalpies, and Gibb energies) of the formation of the complexes analyzed. The B3LYP/Mixed I complexation energies in the gas phase are therefore compared to those obtained using B3LYP/Mixed II and CCSD(T) calculations. Our results pointed out that the deprotonation of the ligand increases the binding affinity independently of the metal cation used. The topological parameters yielded from Quantum Theory of Atom in Molecules (QTAIM) indicate that metal-ligand bonds are partly covalent. The significant reduction of the proton affinity (PA) observed when passing from ligands to complexes in gas-phase confirms the notable enhancement of antioxidant activities of neutral ligands.

A new and recyclable system based on tropin ionic liquids for resolution of several racemic amino acids

A new, recyclable solid-liquid resolution system was developed based on tropin ionic liquids [C_nDTr][L-Pro]₂ for the enantiomeric resolution of racemic phenylalanine and other α-substituted carboxylic acids including tryptophan, tyrosine, benzene glycine and mandelic acid. With racemic phenylalanine as resolution model, effect factors were investigated for better resolution conditions. On the conditions, some efficient resolution were achieved, for instance the e.e. (98%) and product yield (76%) in solid phase for phenylalanine, and the e.e. 99.71% in solid phase for tryptophan. Chiral product was verified with fourier transform infrared spectroscopy (FT-IR), Raman spectrum, thermal gravity analysis (TG), elemental analysis (EA) and chiral HPLC. Further, the resolution mechanism was studied with computer molecular dynamic simulations and UV-vis. The resolution was closely related to the formation of complexes (L-Phe-Cu(Ⅱ)-L-pro^-) and the spatial configuration of D/L-Phe. The system is characteristic of high resolution, no organic solvent, easy isolation of solid-liquid and recycle of all chemical materials as much as possible.

High resolution of racemic phenylalanine with dication imidazolium-based chiral ionic liquids in a solid-liquid two-phase system

A novel solid-liquid two-phase system was developed for the chiral separation of racemic phenylalanine with new dication imidazolium-based chiral ionic liquids. Preliminary experiments showed distinct enantioselectivity in amino acid extraction with the novel solid-liquid two-phase system, more L-enantiomer of amino acid cooperatively interacted with ionic liquids and copper ions to be the solid phase. Various factors, including the alkyl chain length of cations of ionic liquids, the amount of copper acetate, the ratio of n_(ILs)/n_(Cu²⁺), the amount of water and racemic phenylalanine, the resolution time together with the resolution temperature, were systematically investigated for their influence on resolution efficiency. The results showed that, under a certain condition, the enantiomeric excess value and the yield of phenylalanine in liquid phase (mainly containing D-enantiomer) were 67.8% and 96.5%, the enantiomeric excess value and the yield of phenylalanine in solid phase (mainly containing L-enantiomer) were 99.2% and 85.2%. Finally, 2D NMR technology, infrared spectroscopy and molecular simulation method were used to study the interaction mechanism. The results indicated that L-enantiomer of phenylalanine interacts more strongly with chiral ILs and Cu²⁺. The novel system has characteristics of free-organic solvent, simple operation, fast separation process and very high resolution efficiency for racemic phenylalanine. This work could provide a new and alternative resolution approach for other chiral separations.

Synthesis, Biological, and Quantum Chemical Studies of Zn(II) and Ni(II) Mixed-Ligand Complexes Derived from N,N-Disubstituted Dithiocarbamate and Benzoic Acid

Some mixed-ligand complexes of Zn(II) and Ni(II) derived from the sodium salt ofN-alkyl-N-phenyl dithiocarbamate and benzoic acid have been prepared. The complexes are represented as ZnMDBz, ZnEDBz, NiMDBz, and NiEDBz (MD:N-methyl-N-phenyl dithiocarbamate, ED:N-ethyl-N-phenyl dithiocarbamate, and Bz: benzoate); and their coordination behavior was characterized on the basis of elemental analyses, IR, electronic spectra, magnetic and conductivity measurements, and quantum chemical calculations. The magnetic moment measurement and electronic spectra were in agreement with the four proposed coordinate geometries for nickel and zinc complexes and were corroborated by the theoretical quantum chemical calculations. The quantum chemically derived thermodynamics parameters revealed that the formation ofN-methyl-N-phenyl dithiocarbamate complexes is more thermodynamically favourable than that of theN-ethyl-N-phenyl dithiocarbamate complexes. The bioefficacy of the mixed-ligand complexes examined against different microbes showed moderate to high activity against the test microbes. The anti-inflammatory and antioxidant studies of the metal complexes showed that the ethyl substituted dithiocarbamate complexes exhibited better anti-inflammatory and antioxidant properties than the methyl substituted dithiocarbamate complexes.

Synthesis, DFT Calculation, and Antimicrobial Studies of Novel Zn(II), Co(II), Cu(II), and Mn(II) Heteroleptic Complexes Containing Benzoylacetone and Dithiocarbamate

Heteroleptic complexes of zinc(II), copper(II), manganese(II), and cobalt(II) of the types [MLL'(H2O)2]·nH2O and [MLL']·nH2O have been synthesized using sodium N-methyl-N-phenyldithiocarbamate (L) and benzoylacetone (L'). The metal complexes were characterized by elemental analysis, electrical conductance, magnetic susceptibility, infrared (IR), and UV-visible spectroscopic studies. The electrical conductance measurements revealed the nonelectrolytic nature of the synthesized complexes. The results of the elemental analyses, magnetic susceptibility measurements, and electronic spectra inferred that the Zn(II) complex adopted a four-coordinate geometry while the Co(II), Cu(II), and Mn(II) complexes assumed octahedral geometries. The IR spectra showed that the metal ions coordinated with the ligands via the S- and O-donor atoms. The geometry, electronic, and thermodynamic parameters of the complexes were obtained from density functional theory (DFT) calculations. The spin density distributions, relative strength of H-bonds, and thermodynamic parameters revealed that the order of stability of the metal complexes is Mn < Co < Cu > Zn. The agar diffusion methods were used to study the antimicrobial activity of the complexes against two Gram positive bacteria (S. aureus and S. pneumoniae), one Gram negative bacterium (E. coli), and two fungi organisms (A. niger and A. candida) and the complexes showed a broad spectrum of activities against the microbes.

Evaluation of the antiradical activity of hyperjovinol-A utilizing donor-acceptor maps

Hyperjovinol-A ((2-methyl-1-(2,4,6-trihydroxy-3-(3-hydroxy-3,7-dimethyloct-6-enyl)phen yl)propan-1-one) is an acylated phloroglucinol isolated from Hypericum Jovis and exhibiting good antioxidant activity. The study investigates the compound’s antiradical ability on the basis of the electron-donor and electron-acceptor abilities of its conformers, deriving donor and acceptor indexes and mapping them in terms of donor-acceptor maps (DAM). The DAMs of vitamins C and E and of carotene astaxantine are used as comparison references. Calculations were performed at the DFT/BPW91/6-311+G(d,p) level, with optimizations on fully relaxed geometries to obtain the conformers of the neutral molecule in vacuo, and single point calculations to obtain the energies of the cationic and anionic species in vacuo and of the neutral, cationic, and anionic species in water, ethanol, and pentylethanoate. The calculations in solution utilized the polarizable continuum model (PCM). The results indicate that hyperjovinol-A may have better antiradical activity than vitamin C. This is in agreement with experimental results, showing that the antioxidant activity of hyperjovinol-A is comparable with that of the best drugs currently in clinical use. The activity is maintained in solution. The Fukui function f(·) was also calculated for all the conformers of hyperjovinol-A, to identify the regions of highest reactivity.

Sodium and potassium ion directed self-assembled multinuclear assembly of divalent nickel or copper and L-leucine derived ligand

L-leucine derived ligand (H(2)L(L-leu)), KOH, and Ni(II) salt in 2:2:1 ratio self-assembled into a rather large (approximately 13 A) supramolecular assembly with the formula [K{Ni(HL(L-leu))(2)}(3)](+) (1). Structural characterization showed three [Ni(HL(L-leu))(2)] units encapsulated K(+) similar to organic crown ethers/cryptand. Substituting Ni(II) with Cu(II) and K(+) with Na(+) in the above reaction resulted in a set of structurally identical assemblies with the general formula [M'{M(HL(L-leu))(2)}(3)](+), where M' is either K(+) or Na(+) and M is either Cu(II) or Ni(II); [Na{Ni(HL(L-leu))(2)}(3)]ClO(4) (2), Na{Ni(HL(L-leu))(2)}(3)]OTf (3), [K{Cu(HL(L-leu))(2)}(3)]ClO(4) (4), [Na{Cu(HL(L-leu))(2)}(3)]ClO(4) (5), [K{Cu(HL(L-leu))(2)}(3)]NO(3) (6). Electrospray Ionization (ESI)-mass spectra of the assemblies in MeOH showed the retention of assemblies in solution. Visible spectroscopic studies showed retention of assembly 1 in N,N-dimethylformamide (DMF) which is stable even after the addition of 5 equiv of [18]-crown-6. The assemblies in 2-6 show various degrees of dissociation to [M(HL(L-leu))(2)] and M', in stronger H-bonding methanol. The dissociation can be reversed upon addition of excess KNO(3)/NaNO(3) salt. Structural characterization of [Cu(HL(L-leu))(2)(MeCN)] (7) along with its transformation to [K{Cu(HL(L-leu))(2)}(3)](+) in the presence of K(+) salt demonstrated that the assembly formation proceeds through an alkali metal ion induced ligand reorientation within the [Cu(HL(L-leu))(2)] units which is further stabilized by six strong H-bonds holding the assembly. Interestingly, visible spectra of 1 and 2 shows that minor structural changes caused by replacing K(+) with Na(+) is sufficient to shift the d-d transition of Ni(II) by approximately 70 nm, thereby providing an indirect way of distinguishing K(+) and Na(+), none of which have spectroscopic signature in the visible range.

Copper Cation Interactions with Biologically Essential Types of Ligands: A Computational DFT Study

This work presents a systematic theoretical study on Cu(I) and Cu(II) cations in variable hydrogen sulfide-aqua-ammine ligand fields. These ligands model the biologically most common environment for Cu ions. Molecular structures of the complexes were optimized at the density functional theory (DFT) level. Subsequent thorough energy analyses revealed the following trends: (i) The ammine complexes are the most stable, followed by those containing the aqua and hydrogen sulfide ligands, which are characterized by similar stabilization energies. (ii) The most preferred Cu(I) coordination number is 2 in ammine or aqua ligand fields. A qualitatively different binding picture was obtained for complexes with H(2)S ligands where the 4-coordination is favored. (iii) The 4- and 5-coordinated structures belong to the most stable complexes for Cu(II), regardless of the ligand types. Vertical and adiabatic ionization potentials of Cu(I) complexes were calculated. Charge distribution (using the natural population analysis (NPA) method) and molecular orbital analyses were performed to elucidate the nature of bonding in the examined systems. The results provide in-depth insight into the Cu-binding properties and can be, among others, used for the calibration of bioinorganic force fields.

Chemical Exchange Reaction of Glycinatocopper(II) Complex in Water: A Theoretical Study

The axial water exchange on glycinatocopper(II) complexes was theoretically investigated by using density functional theory (DFT) calculations. Glycinatocopper(II) complexes are well-known by the diffusion controlled exchange of axial ligands. Calculations using explicitly coordinating water molecules and solvent models showed that bis-glycinatocopper(II) complexes have a four-coordinate planar structure, in which waters are excluded from the axial positions of Cu(II) due to the Jahn-Teller effect. This may be because coordinating axial waters induce the discrepancy in the most stable ligand field splittings of inner 3d and outer 4d orbitals of the Cu(II) cation. To estimate the reactivity of the axial water exchange, we calculated the rate constant by calculating Gibbs free energies for the activation. As a result, we obtained the rate constant as k = 3.61 x 10(10) s(-1) in aqueous solution at T = 298.15 K. This rate constant is slightly larger than that of the diffusion controlled exchange of axial waters, which is experimentally observed in the order of 10(9) s(-1). Finally, we determined the structures of tris-glycinatocopper(II) complexes. It was consequently found that the third glycine is coordinated to Cu with the amino groups as experimentally observed.

Theoretical model of the aqua-copper [Cu(H2O)5]+cation interactions with guanine

Pentaaqua complexes of Cu(I) with guanine were optimized at the DFT B3PW91/6-31G(d) level. For the most stable structures, vibration frequencies and NBO charges were computed followed by energy analyses. The order of individual conformers was very sensitive to the method and basis sets used for the calculation. Several conformers are practically degenerated in energy. The inclusion of an entropy term changes the order of the conformers' stability. Water molecules associated at the N9 position of guanine are favored by the inclusion of the entropy correction. Bonding energies of Cu-O(aqua) interactions were estimated to be about 60 kcal mol(-1) and for Cu-N7 bonding in the range of 75-83 kcal mol(-1). The broad range in Cu-N interaction energies demonstrates the role of induction effects caused by water molecules associated at the various sites of guanine. The charge distribution of the guanine molecule is changed remarkably by the coordination of a Cu(I) cation, which can also change the base-pairing pattern of the guanine.