Ground and Low-Lying States of Cu2+−H2O. A Difficult Case for Density Functional Methods

The Role of the CuCl Active Complex in the Stereoselectivity of the Salt-Induced Peptide Formation Reaction: Insights from Density Functional Theory Calculations

The salt-induced peptide formation (SIPF) reaction is a prebiotically plausible mechanism for the spontaneous polymerization of amino acids into peptides on early Earth. Experimental investigations of the SIPF reaction have found that in certain conditions, the l enantiomer is more reactive than the d enantiomer, indicating its potential role in the rise of biohomochirality. Previous work hypothesized that the distortion of the CuCl active complex toward a tetrahedral-like structure increases the central chirality on the Cu ion, which amplifies the inherent parity-violating energy differences between l- and d-amino acid enantiomers, leading to stereoselectivity. Computational evaluations of this theory have been limited to the protonated-neutral l + l forms of the CuCl active complex. Here, density functional theory methods were used to compare the energies and geometries of the homochiral (l + l and d + d) and heterochiral (l + d) CuCl-amino acid complexes for both the positive-neutral and neutral-neutral forms for alanine, valine, and proline. Significant energy differences were not observed between different chiral active complexes (i.e., d + d, l + l vs. l + d), and the distortions of active complexes between stereoselective systems and non-selective systems were not consistent, indicating that the geometry of the active complex is not the primary driver of the observed stereoselectivity of the SIPF reaction.

Proportional coexistence of okanin chalcone glycoside and okanin flavanone glycoside in Bidens pilosa leaves and theoretical investigation on the antioxidant properties of their aglycones

Bidens pilosa plant has been shown to produce okanin flavanone glycoside and its chalcone derivative. In most other plants, due to chalcone isomerase enzyme, the flavanone tends to exist in higher proportions than their chalcone precursors. Herein we have utilized liquid chromatography-mass spectrometry approach and shown that within the leaves of Bidens pilosa plant the two okanin glycosides exist in unusual equal proportional distribution, which indicates that Bidens pilosa plant is an alternative rich source of these highly sought-after antioxidant molecules. The aglycone okanin chalcone (ONC) and okanin flavanone (ONF) have experimentally been shown to exhibit antioxidant activity. However, experimental findings have not conclusively determined which of the two compounds is a more potent antiradical than the other. Herein, the density functional theory (DFT) method is utilized to establish, from structural and thermodynamic energetic considerations, the preferred antioxidant molecule between the two aglycone okanins. A theoretical study on the antioxidant properties of ONC and ONF has been performed by considering their radical scavenging and metal cation (Mⁿ⁺, where M = Cu(II) or Fe (III)) chelation ability. The study has been performed using B3LYP/6-31 + G(d,p) method. In the case of the metal chelation mechanism, the LANL2DZ pseudo-potential was selected to describe the selected Mⁿ⁺ cations. The results of the study suggest that ONC is a better radical scavenger than ONF because of the extended electron delocalization on its neutral radical, which is due to the presence of conjugation within the ONC neutral radical after hydrogen atom abstraction. In the metal chelation mechanism, it is noted that the binding energies depend on the media, the nature of the ligand and the cation and the cation coordination site on the ligand. The charge and the spin density on Mⁿ⁺ decrease on coordination to the ligand. The ability of the ligands to reduce Mⁿ⁺ cations, coupled with the strong Mⁿ⁺ binding properties, has significant implication on the antioxidant ability of both okanins. However, since ONC⋅⋅⋅M⁺ⁿ interaction results in higher binding energy than ONF⋅⋅⋅M⁺ⁿ interaction, the implication is that ONC is a preferred free metal ion chelator than ONF.

Cyclo[18]carbon: the smallest all-carbon electron acceptor

The recently synthesized C₁₈ molecule demonstrates strong electron acceptor properties similar to C₆₀.

Hydrogenated polycyclic aromatic hydrocarbons (HnPAHs) as catalysts for hydrogenation reactions in the interstellar medium: a quantum chemical model

Density functional studies show that neutral H_nPAHs are able to catalyze the formation of water with no activation barrier.

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.

Lanthanide and transition metal complexes of bioactive coumarins: molecular modeling and spectroscopic studies

The present paper summarizes theoretical and spectroscopic investigations on a series of active coumarins and their lanthanide and transition metal complexes with application in medicine and pharmacy. Molecular modeling as well as IR, Raman, NMR and electronic spectral simulations at different levels of theory were performed to obtain important molecular descriptors: total energy, formation energy, binding energy, stability, conformations, structural parameters, electron density distribution, molecular electrostatic potential, Fukui functions, atomic charges, and reactive indexes. The computations are performed both in gas phase and in solution with consideration of the solvent effect on the molecular structural and energetic parameters. The investigations have shown that the advanced computational methods are reliable for prediction of the metal-coumarin binding mode, electron density distribution, thermodynamic properties as well as the strength and nature of the metal-coumarin interaction (not experimentally accessible) and correctly interpret the experimental spectroscopic data. Known results from biological tests for cytotoxic, antimicrobial, anti-fungal, spasmolytic and anti-HIV activities on the studied metal complexes are reported and discussed.

A stable, mononuclear, cationic Pt(iii) complex stabilised by bulky N-heterocyclic carbenes

The thermally stable, mononuclear, paramagnetic Pt(iii) complex [PtI₂(IPr)₂][BAr^F] exhibiting a square-planar structure has been isolated and fully characterised.

Strong interactions between copper halides and unsaturated systems: new metallocycles? Or the importance of deformation

The complexes formed by CuF with CC double and triple bonds have been studied at the MP2 and CCSD(T) computational levels. The interaction of CuF with acetylene, ethylene and their fluoro derivatives is very strong, with interaction energies close to those of conventional covalent bonds. Hence, these complexes could be actually viewed as a new kind of metallocycles, with significantly strong Cu-C linkages. All electronic indexes analyzed by means of the AIM, ELF and NBO formalisms, indicate that the strength of the interaction should increase with the number of fluorine substituents in both series of compounds. Surprisingly, however, although both series of compounds exhibit the same bonding arrangements, they follow opposite stability trends and the expected increase of the interaction energies with the number of fluorine substituents is only observed in the acetylene series. The reason for this unexpected behavior is once more associated with the effects triggered by the distortion of the interacting subunits. Deformation not only has a direct energetic cost but dramatically affects the intrinsic properties of the interacting systems.

A quantum-chemical study of the binding ability of βXaaHisGlyHis towards copper(II) ion

The present study analyzed binding of Cu²⁺ to tetrapeptides in water solution at several levels of theoretical approximation. The methods used to study the energetic and structural properties of the complexes in question include semiempirical hamiltonians, density functional theory as well as ab initio approaches including electron correlation effects. In order to shed light on the character of interactions between Cu²⁺ and peptides, which are expected to be mainly electrostatic in nature, decomposition of interaction energy into physically meaningful components was applied.

Three dimensional models of Cu(2+)-Aβ(1-16) complexes from computational approaches

Elucidation of the coordination of metal ions to Aβ is essential to understand their role in its aggregation and to rationally design new chelators with potential therapeutic applications in Alzheimer disease. Because of that, in the last 10 years several studies have focused their attention in determining the coordination properties of Cu(2+) interacting with Aβ. However, more important than characterizing the first coordination sphere of the metal is the determination of the whole Cu(2+)-Aβ structure. In this study, we combine homology modeling (HM) techniques with quantum mechanics based approaches (QM) to determine plausible three-dimensional models for Cu(2+)-Aβ(1-16) with three histidines in their coordination sphere. We considered both ε and δ coordination of histidines 6, 13, and 14 as well as the coordination of different possible candidates containing oxygen as fourth ligand (Asp1, Glu3, Asp7, Glu11, and CO(Ala2)). Among the 32 models that enclose COO(-), the lowest energy structures correspond to [O(E3),N(δ)(H6),N(ε)(H13),N(ε)(H14)] (1), [O(E3),N(δ)(H6),N(δ)(H13),N(δ)(H14)] (2), and [O(D7),N(ε)(H6),N(δ)(H13),N(δ)(H14)] (3). The most stable model containing CO(Ala2) as fourth ligand in the Cu(2+) coordination sphere is [O(c)(A2),N(ε)(H6),N(δ)(H13),N(ε)(H14)] (4). An estimation of the relative stability between Glu3 (1) and CO(Ala2) (4) coordinated complexes seems to indicate that the preference for the latter coordination may be due to solvent effects. The present results also show the relationship between the peptidic and metallic moieties in defining the overall geometry of the complex and illustrate that the final stability of the complexes results from a balance between the metal coordination site and amyloid folding upon complexation.

Modeling the interactions of phthalocyanines in water: from the Cu(II)-tetrasulphonate to the metal-free phthalocyanine

A quantum and statistical study on the effects of the ions Cu(2+) and SO(3)(-) in the solvent structure around the metal-free phthalocyanine (H(2)Pc) is presented. We developed an ab initio interaction potential for the system CuPc-H(2)O based on quantum chemical calculations and studied its transferability to the H(2)Pc-H(2)O and [CuPc(SO(3))(4)](4-)-H(2)O interactions. The use of the molecular dynamics technique allows the determination of energetic and structural properties of CuPc, H(2)Pc, and [CuPc(SO(3))(4)](4-) in water and the understanding of the keys for the different behaviors of the three phthalocyanine (Pc) derivatives in water. The inclusion of the Cu(2+) cation in the Pc structure reinforces the appearance of two axial water molecules and second-shell water molecules in the solvent structure, whereas the presence of SO(3)(-) anions implies a well defined hydration shell of about eight water molecules around them making the macrocycle soluble in water. Debye-Waller factors for axial water molecules have been obtained in order to examine the potential sensitivity of the extended x-ray absorption fine structure technique to detect the axial water molecules.

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.

Effect of metal ions (Li+, Na+, K+, Mg2+, Ca2+, Ni2+, Cu2+ and Zn2+) and water coordination on the structure and properties of L-histidine and zwitterionic L-histidine

Interactions between metal ions and amino acids are common both in solution and in the gas phase. The effect of metal ions and water on the structure of L-histidine is examined. The effect of metal ions (Li+, Na+, K+, Mg2+, Ca2+, Ni2+, Cu2+ and Zn2+) and water on structures of His·M(H2O)m, m=0.1 complexes have been determined theoretically employing density functional theories using extended basis sets. Of the five stable complexes investigated the relative stability of the gas-phase complexes computed with DFT methods (with one exception of K+ systems) suggest metallic complexes of the neutral L-histidine to be the most stable species. The calculations of monohydrated systems show that even one water molecule has a profound effect on the relative stability of individual complexes. Proton dissociation enthalpies and Gibbs energies of L-histidine in the presence of the metal cations Li+, Na+, K+, Mg2+, Ca2+, Ni2+, Cu2+ and Zn2+ were also computed. Its gas-phase acidity considerably increases upon chelation. Of the Lewis acids investigated, the strongest affinity to L-histidine is exhibited by the Cu2+ cation. The computed Gibbs energies ΔG are negative, span a rather broad energy interval (from -130 to -1,300 kJ/mol), and upon hydration are appreciably lowered.

Density functional theory for transition metals and transition metal chemistry

We introduce density functional theory and review recent progress in its application to transition metal chemistry. Topics covered include local, meta, hybrid, hybrid meta, and range-separated functionals, band theory, software, validation tests, and applications to spin states, magnetic exchange coupling, spectra, structure, reactivity, and catalysis, including molecules, clusters, nanoparticles, surfaces, and solids.

Effect of metal ions (Li+, Na+, K+, Mg2+, Ca2+, Ni2+, Cu2+, and Zn2+) and water coordination on the structure and properties of L-arginine and zwitterionic L-arginine

Interactions between metal ions and amino acids are common both in solution and in the gas phase. The effect of metal ions and water on the structure of L-arginine is examined. The effects of metal ions (Li(+), Na(+), K(+), Mg(2+), Ca(2+), Ni(2+), Cu(2+), and Zn(2+)) and water on structures of Arg x M(H2O)m , m = 0, 1 complexes have been determined theoretically by employing the density functional theories (DFT) and using extended basis sets. Of the three stable complexes investigated, the relative stability of the gas-phase complexes computed with DFT methods (with the exception of K(+) systems) suggests metallic complexes of the neutral L-arginine to be the most stable species. The calculations of monohydrated systems show that even one water molecule has a profound effect on the relative stability of individual complexes. Proton dissociation enthalpies and Gibbs energies of arginine in the presence of the metal cations Li(+), Na(+), K(+), Mg(2+), Ca(2+), Ni(2+), Cu(2+), and Zn(2+) were also computed. Its gas-phase acidity considerably increases upon chelation. Of the Lewis acids investigated, the strongest affinity to arginine is exhibited by the Cu(2+) cation. The computed Gibbs energies DeltaG(o) are negative, span a rather broad energy interval (from -150 to -1500 kJ/mol), and are appreciably lowered upon hydration.

Cu2+/+ cation coordination to adenine--thymine base pair. Effects on intermolecular proton-transfer processes

Intermolecular proton-transfer processes in the Watson & Crick adenine-thymine Cu+ and Cu2+ cationized base pairs have been studied using the density functional theory (DFT) methods. Cationized systems subject to study are those resulting from cation coordination to the main basic sites of the base pair, N7 and N3 of adenine and O2 of thymine. For Cu+ coordinated to N7 or N3 of adenine, only the double proton-transferred product is found to be stable, similarly to the neutral system. However, when Cu+ interacts with thymine, through the O2 carbonyl atom, the single proton transfer from thymine to adenine becomes thermodynamically spontaneous, and thus rare forms of the DNA bases may spontaneously appear. For Cu2+ cation, important effects on proton-transfer processes appear due to oxidation of the base pair, which stabilizes the different single proton-transfer products. Results for hydrated systems show that the presence of the water molecules interacting with the metal cation (and their mode of coordination) can strongly influence the ability of Cu2+ to induce oxidation on the base pair.

Reactivity of the non stoichiometric Ni3O4 phase supported at the Pd(100) surface: interaction with Au and other transition metal atoms

In this work we have studied, employing ab initio periodic calculations, the interaction between the non-stoichiometric Ni(3)O(4) monolayer (a rhombic distribution of vacancies hereafter referred to as RH-Ni(3)O(4)) supported on the Pd(100) surface and several transition metal atoms (Ni, Cu, Pd, Pt, Ag, Au). The interaction produces a regular array of metal centers uniformly distributed in size and shape. According to the size of the atom and the ionization potential, the nature of the interaction ranges from an essentially electrostatic one to a polar-covalent one. The chemical reactivity versus the CO molecule of the overlayer resulting from the saturation with Au atoms of the Ni vacancies has been investigated.