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

Microsolvation of cobalt, nickel, and copper atoms with ammonia: a theoretical study of the solvated electron precursors

CONTEXT

The s-block metals dissolved in ammonia form metal-ammonia complexes with diffuse electrons which could be used for redox catalysis. In this theoretical paper, we investigated the possibility of the d-bloc transition metals (Mn, Fe, Co, Ni, and Cu) solvated by ammonia. It has been demonstrated that both Mn and Fe atoms undergo into an oxidative reaction with NH3 forming an inserted species, HMNH2. On the contrary, the Co, Ni, and Cu atoms can accommodate four NH3, via the coordination bond, to form the first solvation sphere within C2v, D2d, and Td point groups, respectively. Addition of a fifth NH3 constitute the second solvation shell by forming hydrogen bond with the other NH3s. Interestingly, M(NH3)4 (M = Co, Ni, and Cu) is a so-called solvated electron precursor and should be considered as a monocation M(NH3)4+ kernel in tight contact with one electron distributed over its periphery. This nearly free electron could be used to capture a CO2 molecule and engages in a reduction reaction.

METHODS

Geometry optimization of the stationary points on the potential energy surface was performed using density functional theory - CAM-B3LYP functional including the GD3BJ dispersion contribution - in combination with the 6-311 +  + G(2d, 2p) basis set for all the atoms. All first-principles calculations were performed using the Gaussian 09 quantum chemical packages. The natural electron configuration of transition atom engaged in the compounds has been found using the natural bond orbital (NBO) method. We used the EDR (electron delocalization range) approach to analyze the structure of solvated electrons in real space. We also used the electron localization function (ELF) to measure the degree of electronic localization within a chemical compound. The EDR and ELF analyses are done using the TopMod and Multiwfn packages, respectively.

Structures of the solvated copper(ii) ion in ammonia at various temperatures

We investigated theoretically the structures and relative stabilities of the solvated copper(ii) ion in ammonia, Cu²⁺(NH₃)_n, n = 1–10.

Molecular Facts on the Structure and Dynamics of Electrolyte Species in Cu-Cl Cycle for Hydrogen Generation: An Insight from Molecular Dynamic Simulations

The Cu complex, which is the key chemical species in well-known Cu-Cl hybrid thermochemical cycles and also in numerous metal hydrometallurgical and sedimentary deposit processes, displays a wide variety of structural and dynamical characteristics that are further complicated by the presence of multiple oxidation states of Cu ions with different coordination chemistries, therefore they are difficult to explore from experiments alone. In this article, an attempt has been made to understand the coordination behavior of the Cu complex using MD simulations. The study provides compelling evidence of the experimentally observed multiple stoichiometries of Cu ions, i.e., 1:6:0, 1:5:1, and 1:4:2 for Cu⁺:H₂O:Cl^- and 1:6:0 for Cu²⁺:H₂O:Cl^-. The presence of the anionic Cu complex, [Cu⁺Cl₂]^-·2H₂O, [Cu⁺Cl₂]^-·3H₂O, [Cu²⁺Cl₃]^-·H₂O, and [Cu²⁺Cl₃]^-·2H₂O, was captured in the presence of excess chloride ions. Furthermore, the probability distribution profiles have been estimated to determine the most possible complex in the considered systems. The results establish structural and dynamical reformation of the Cu complex with change in the salt concentration or variation in the solvent medium in which they are dissolved. Moreover, the structure and kinetics of the Cu ions in the Cu-Cl electrolyzer have been explored over a large range of the electric field by extending the simulated systems for varied strengths of the electric fields. It has been observed that with an increase in the strength of the electric field, the water molecules lose their coordination strength with central Cu ions, which, on the other hand, results in a significant change in the structure of the captured complex. The diffusion dynamics of the ions is altered while applying the electric field, which is furthermore modified while increasing the strength of electric field beyond a critical limit. In fact, the diffusion mechanism of the ions was seen to be transformed from Brownian-like to linear motion and then to hopping diffusion with the increasing strength of the electric field. To the best of our knowledge, this is the first time when the multiple oxidation states of the Cu ion are explored using MD simulations, and the coexisting pictures of the multiple coordinations and the solvent effects have been clearly revealed. Also to date, the present article is the first one to report the insights of the structure and the dynamics of the ions in the Cu-Cl electrolyzer over a wide range of the electric field. The present studies will be very helpful in understanding the mechanism involved in numerous metal hydrometallurgical and sedimentary deposit processes and to comprehend the analogies involved in the electrode reactions of the Cu-Cl cycle for hydrogen generation.

Side Reactions with an Equilibrium Constraint: Detailed Mechanism of the Substitution Reaction of Tetraplatin with dGMP as a Starting Step of the Platinum(IV) Reduction Process

Two possible pathways of the substitution reaction within the reduction process of the Pt^IV(DACH)Cl₄ by dGMP are compared: associative reaction course and autocatalytic Basolo-Pearson mechanisms. Since two forms: single-protonated and fully deprotonated phosphate group of dGMP are present in equilibrium at neutral and mildly acidic solutions, consideration of a side reactions scheme with acido-basic equilibrium-constraint is a very important model for obtaining reliable results. The examined complexes are optimized at the B3LYP-GD3BJ/6-31G(d) level with the COSMO implicit solvation model and Klamt's radii used for cavity construction. Energy characteristics and thermodynamics for all reaction branches are determined using the B3LYP-GD3BJ/6-311++G(2df,2pd)/IEF-PCM/scaled-UAKS level with Wertz's entropy corrections. Rate constants are estimated for each individual branch according to Eyring's transition state theory (TST), averaged according to equilibrium constraint and compared with available experimental data. The determined reaction barriers of the autocatalytic pathway fairly correspond with experimental values. Furthermore, autocatalytic reaction of tetraplatin and its two analogues complexes [Pt^IV(en)Cl₄ and Pt^IV(NH₃)₂Cl₄] are explored and compared with measured data in order to examined general reaction descriptors.

Interactions of the "piano-stool" [ruthenium(II)(η(6) -arene)(quinolone)Cl](+) complexes with water; DFT computational study

Full optimizations of stationary points along the reaction coordinate for the hydration of several quinolone Ru(II) half-sandwich complexes were performed in water environment using the B3PW91/6-31+G(d)/PCM/UAKS method. The role of diffuse functions (especially on oxygen) was found crucial for correct geometries along the reaction coordinate. Single-point (SP) calculations were performed at the B3LYP/6-311++G(2df,2pd)/DPCM/saled-UAKS level. In the first part, two possible reaction mechanisms-associative and dissociative were compared. It was found that the dissociative mechanism of the hydration process is kinetically slightly preferred. Another important conclusion concerns the reaction channels. It was found that substitution of chloride ligand (abbreviated in the text as dechlorination reaction) represents energetically and kinetically the most feasible pathway. In the second part the same hydration reaction was explored for reactivity comparison of the Ru(II)-complexes with several derivatives of nalidixic acid: cinoxacin, ofloxacin, and (thio)nalidixic acid. The hydration process is about four orders of magnitude faster in a basic solution compared to neutral/acidic environment with cinoxacin and nalidixic acid as the most reactive complexes in the former and latter environments, respectively. The explored hydration reaction is in all cases endergonic; nevertheless the endergonicity is substantially lower (by ∼6 kcal/mol) in basic environment. © 2016 Wiley Periodicals, Inc.

Intensified effects of multi-Cu modification on the electronic properties of the modified base pairs containing hetero-ring-expanded pyrimidine bases

Novel DNA base pair derivatives (A2CunU, A2CunC, G3CunU, and G3CunC) are designed by aromatic expansion of pyrimidine bases with four kinds of hetero-rings (denoted by nC and nU, n = 1, 2, 3, and 4) and metal-decoration through Cu replacement of hydrogens in the Watson-Crick hydrogen bond region. Their structures and properties are calculated for examining the cooperative effects of the two modification ways. The calculated results reveal that multiple Cu decoration makes up the deficiencies of size-expansion, and exhibits not only increase of structural stability and reduction of ionization potentials, but also ideal shrink of the HOMO-LUMO gaps, notable enhancement of interbase coupling as well as remarkable redshifts of π → π* transitions for all M-x modified base pairs. The decrease extents of the gaps and ionization potentials follow the same order G3CunU > G3CunC > A2CunU > A2CunC, and in each series (denoted by different n), the gaps, ionization potentials and first π → π* transition energies have an order of 4 < 1 < 2 < 3. The Cu d orbitals function as bridges for π electron delocalization on the conjugated aromatic rings of two bases, leading to an enhancement of transverse electronic communication, as verified by spin density delocalization, orbital composition changes, redshift of the π → π* transition and also advocated by the electron-sharing indexes such as delocalization index, Mayer bond orders and multicenter bonding. Electron localization function ELF-π isosurfaces above the molecular plane further suggested that effective longitudinal conduction is closely relevant with the bicyclic domain involving good electron delocalization and strong π-π stacking between layers. This work presents theoretical evidence for the cooperative effects of metal decoration and ring-expansion modifications on the electronic properties of the modified base pairs and also proves that the base pairs designed here could be competent building blocks for the DNA-based nanowires with improved electron activity and excellent conductivity.

DFT/TDDFT insights into the chemistry, biochemistry and photophysics of copper coordination compounds

Highlighting the recent progress in DFT/TDDFT application to coordination chemistry of copper.

DFT studies of structural and some spectral parameters of copper(II) complexes with N,N,N',N″-tetrakis (2-hydroxyethyl/propyl) ethylenediamine and tris(2-hydroxyethyl)amine

The structures and some spectral parameters of three copper(II) complexes; [Cu(THEEN)(H2O)](PIC)2 (1), [Cu(THPEN)] (PIC)2 C3H8O (2) and [Cu(TEAH3)(PIC)] (PIC)⋅(H2O) (3), previously synthesized and characterized by X-ray diffraction, are here computationally studied by using density functional theory (DFT) in its hybrid form B3LYP. In these complexes, THEEN is N,N,N',N″-tetrakis(2-hydroxyethyl) ethylenediamine and THPEN is N,N,N',N″-tetrakis(2-hydroxypropyl) ethylenediamine, tetrapodal ligands and TEAH3 is tris(2-hydroxyethyl)amine, a tripodal ligand. The primary coordination sphere of copper(II) ion in complexes (1), (2) and (3) are optimized, structural parameters are calculated, vibrational bands are assigned and energy gaps of frontier orbital (HOMO-LUMO) have been calculated with B3LYP/6-31G/LANL2DZ level of theory using DMSO as solvent. The calculated geometric and spectral results reproduced the experimental data with well agreement. Theoretical calculated molecular orbitals (HOMO-LUMO) and their energies have been calculated that suggest charge transfer occurs within the complexes.

Theoretical investigation on DNA/RNA base pairs mediated by copper, silver, and gold cations

B3LYP density functional based computations were performed in order to characterize the interactions present in some Cu(+), Ag(+), and Au(+) metal ion-mediated DNA and RNA base pairs from both structural and electronic points of view. Examined systems involve as ligands canonical Watson-Crick, Hoogsteen and Wobble base pairs. Two artificial Hoogsteen base pairs were also taken into account. Binding energy values indicate that complexes involving silver cations are less stable than those in which copper or gold are present, and propose a similar behaviour for these two latter ions. The nature of the bond linking metal ions and bases was described by the NBO analysis that suggests metal coordinative interactions to be covalent. An evaluation of the dispersion contributions for the investigated systems was performed with the B3LYP-D3 functional.

Development and validation of a ReaxFF reactive force field for Cu cation/water interactions and copper metal/metal oxide/metal hydroxide condensed phases

To enable large-scale reactive dynamic simulations of copper oxide/water and copper ion/water interactions we have extended the ReaxFF reactive force field framework to Cu/O/H interactions. To this end, we employed a multistage force field development strategy, where the initial training set (containing metal/metal oxide/metal hydroxide condensed phase data and [Cu(H(2)O)(n)](2+) cluster structures and energies) is augmented by single-point quantum mechanices (QM) energies from [Cu(H(2)O)(n)](2+) clusters abstracted from a ReaxFF molecular dynamics simulation. This provides a convenient strategy to both enrich the training set and to validate the final force field. To further validate the force field description we performed molecular dynamics simulations on Cu(2+)/water systems. We found good agreement between our results and earlier experimental and QM-based molecular dynamics work for the average Cu/water coordination, Jahn-Teller distortion, and inversion in [Cu(H(2)O)(6)](2+) clusters and first- and second-shell O-Cu-O angular distributions, indicating that this force field gives a satisfactory description of the Cu-cation/water interactions. We believe that this force field provides a computationally convenient method for studying the solution and surface chemistry of metal cations and metal oxides and, as such, has applications for studying protein/metal cation complexes, pH-dependent crystal growth/dissolution, and surface catalysis.

Linear energy relationships for the octahedral preference of Mg, Ca and transition metal ions

The geometry, atomic charges, force constants, and relative energies of the symmetric and distorted M(2+)(H(2)O)(4)(F(-))(2), M(3+)(H(2)O)(4)(F(-))(2), M(2+)(H(2)O)(3)(F(-))(2), and M(3+)(H(2)O)(3)(F(-))(2) metal complexes, M = Mg, Ca, Co, Cu, Fe, Mn, Ni, Zn, Cr, V, were calculated by using the B3LYP/TZVP density functional method in both gas phase and aqueous solution, modeled using the polarized continuum model. The deformation energy associated with moving one water ligand 12 degrees from the initial "octahedral" arrangement, in which all O-M-O, O-M-F, and F-M-F angles are either 90 degrees or 180 degrees, was calculated to examine the angular ligand flexibility. For all M(2+)(H(2)O)(4)(F(-))(2) complexes, this distortion increased the energy of the complex in proportion to the electrostatic potential-derived (ESP) charge of the metal, and in proportion to D(-10), where D is the distance from the distorted ligand to its closest neighbor. The octahedral stability was further examined by calculating the energies for the removal of a water ligand from the octahedral complex to form a square-pyramidal or trigonal-bipyramidal complex. The octahedral preference, defined as the negative of the corresponding binding energy of the ligand, was found to linearly correlate with the ESP charge of the metal in both the gas phase and aqueous solution. The obtained results indicate that quantum-mechanical covalent effects are of secondary importance for both the flexibility and the octahedral preference of M(2+)(H(2)O)(4)(F(-))(2) and M(3+)(H(2)O)(4)(F(-))(2) complexes. This conclusion and supporting data are important for the development of consistent molecular mechanical force fields of the studied metal ions.

Dissociations of copper(II)-containing complexes of aromatic amino acids: radical cations of tryptophan, tyrosine, and phenylalanine

The dissociations of two types of copper(II)-containing complexes of tryptophan (Trp), tyrosine (Tyr), or phenylalanine (Phe) are described. The first type is the bis-amino acid complex, [Cu(II)(M)(2)].(2+), where M = Trp, Tyr, or Phe; the second [Cu(II)(4Cl-tpy)(M)].(2+), where 4Cl-tpy is the tridendate ligand 4'-chloro-2,2':6',2''-terpyridine. Dissociations of the Cu(ii) bis-amino acid complexes produce abundant radical cation of the amino acid, M.(+), and/or its secondary products. By contrast, dissociations of the 4Cl-tpy-bearing ternary complexes give abundant M.(+) only for Trp. Density functional theory (DFT) calculations show that for Tyr and Phe, amino-acid displacement reactions by H(2)O and CH(3)OH (giving [Cu(II)(4Cl-tpy)(H(2)O)].(2+) and [Cu(II)(4Cl-tpy)(CH(3)OH)].(2+)) are energetically more favorable than dissociative electron transfer (giving M.(+) and [Cu(I)(4Cl-tpy)](+)). The fragmentation pathway common to all these [Cu(II)(4Cl-tpy)(M)].(2+) ions is the loss of NH(3). DFT calculations show that the loss of NH(3) proceeds via a "phenonium-type" intermediate. Dissociative electron transfer in [Cu(II)(4Cl-tpy)(M-NH(3))].(2+) results in [M-NH(3)].(+). The [Phe-NH(3)] (+) ion dissociates facilely by eliminating CO(2) and giving a metastable phenonium-type ion that rearranges readily into the styrene radical cation.

Molecular Dynamics Simulations of Cu(II) and the PHGGGWGQ Octapeptide

The interaction between Cu2+ and the copper-binding octapeptide region in the human prion protein has been investigated by molecular dynamics simulations. In total four different nonbonded and bonded models were used in the study. Charge sets containing atomic partial charges were developed for these models. Out of the considered models, the bonded model performed physically in the most correct way. The simulations with the bonded model showed that the water molecules in the axial position are very labile. The tryptophan indole ring can remain in a stable position on top of the equatorial coordination plane of copper without water mediation. Strong aromatic interaction was observed between the imidazole and indole rings. The nonbonded models showed a tendency for water-mediated interaction between the copper ion and different carbonyl oxygen atoms. In the case of the bonded model, a carbonyl group could also interact directly with the copper ion in one of the apical position.

Gas-phase reactivity of 2,7-dimethyl-[1,2,4]-triazepine thio derivatives toward Cu+ cation: a DFT study

The gas-phase interactions of 2,7-dimethyl-[1,2,4]-triazepine and its thio derivatives with Cu+ were studied through the use of high-level density functional theory (DFT) calculations. The structure of all possible tautomers and their conformers was optimized at the B3LYP/6-31G(d) level of theory. Final energies were obtained at the B3LYP/6-311+G(2df,2p) level. It has been found that the direct association of Cu+ occurs at the oxygen atom attached to position 3 in the case of the dioxo derivative and at the sulfur atom in all other cases. For the dithio derivatives, the global minimum of the PES corresponds to the structure in which the metal ion bridges between the heteroatom at position 3 and the nitrogen atom at position 4 of the corresponding enolic tautomer, forming a four-membered ring structure; for the dioxo derivative, this conformer competes with the ketone tautomer. Moreover, the isomerization processes leading from the most stable adduct to the other stable conformers were investigated. Among all the considered compounds, the 3,5-dithiotriazepines-Cu+ is found to be the one that associates Cu+ more tightly in the gas phase. The calculated Cu+ binding energies show a good correlation with the experimental proton affinities.

Localized Orbital Corrections for the Calculation of Ionization Potentials and Electron Affinities in Density Functional Theory

This paper describes the extension of a previously reported empirical localized orbital correction model to the correction of ionization potential energies (IP) and electron affinities (EA) for atoms and molecules of first and second row elements. The B3LYP localized orbital correction version of the model (B3LYP-LOC) uses 22 heuristically determined parameters that improve B3LYP DFT IP and EA energy calculations on the G2 data set of 134 molecules from a mean absolute deviation (MAD) from experiment of 0.137 to 0.039 eV. The method significantly reduces the number of outliers and overall MAD to error levels below that achieved with G2 wave function based theory; furthermore, the new model has zero additional computational cost beyond standard DFT calculations. Although the model is heuristic and is based on a multiple linear regression to experimental errors, each of the parameters is justified on physical grounds, and each provides insight into the fundamental limitations of DFT, most importantly the failure of current DFT methods to accurately account for nondynamical electron correlation.