The Jahn-Teller effect in solutions of flexible molecules: a neutron diffraction study on the structure of a Cu2+solution in ethylene glycol

Hydration numbers of biologically relevant divalent metal cations from ab initio molecular dynamics and continuum solvation methods

Hydration and, in particular, the coordination number of a metal ion is of paramount importance as it defines many of its (bio)physicochemical properties. It is not only essential for understanding its behavior in aqueous solutions but also determines the metal ion reference state and its binding energy to (bio)molecules. In this paper, for divalent metal cations Ca2+, Cd2+, Cu2+, Fe2+, Hg2+, Mg2+, Ni2+, Pb2+, and Zn2+, we compare two approaches for predicting hydration numbers: (1) a mixed explicit/continuum DFT-D3//COSMO-RS solvation model and (2) density functional theory based ab initio molecular dynamics. The former approach is employed to calculate the Gibbs free energy change for the sequential hydration reactions, starting from [M(H2O)2]2+ aqua complexes to [M(H2O)9]2+, allowing explicit water molecules to bind in the first or second coordination sphere and determining the most stable [M(H2O)n]2+ structure. In the latter approach, the hydration number is obtained by integrating the ion-water radial distribution function. With a couple of exceptions, the metal ion hydration numbers predicted by the two approaches are in mutual agreement, as well as in agreement with the experimental data.

[Cu(aq)]2+ is structurally plastic and the axially elongated octahedron goes missing

High resolution (k = 18 Å^-1 or k = 17 Å^-1) copper K-edge EXAFS and MXAN (Minuit X-ray Absorption Near Edge) analyses have been used to investigate the structure of dissolved [Cu(aq)]²⁺ in 1,3-propanediol (1,3-P) or 1,5-pentanediol (1,5-P) aqueous frozen glasses. EXAFS analysis invariably found a single axially asymmetric 6-coordinate (CN6) site, with 4×O_eq = 1.97 Å, O_ax1 = 2.22 Å, and O_ax2 = 2.34 Å, plus a second-shell of 4×O_water = 3.6 Å. However, MXAN analysis revealed that [Cu(aq)]²⁺ occupies both square pyramidal (CN5) and axially asymmetric CN6 structures. The square pyramid included 4×H₂O = 1.95 Å and 1×H₂O = 2.23 Å. The CN6 sites included either a capped, near perfect, square pyramid with 5×H₂O = 1.94 ± 0.04 Å and H₂O_ax = 2.22 Å (in 1,3-P) or a split axial configuration with 4×H₂O = 1.94, H₂O_ax1 = 2.14 Å, and H₂O_ax2 = 2.28 Å (in 1,5-P). The CN6 sites also included an 8-H₂O second-shell near 3.7 Å, which was undetectable about the strictly pyramidal sites. Equatorial angles averaging 94° ± 5° indicated significant departures from tetragonal planarity. MXAN assessment of the solution structure of [Cu(aq)]²⁺ in 1,5-P prior to freezing revealed the same structures as previously found in aqueous 1M HClO₄, which have become axially compressed in the frozen glasses. [Cu(aq)]²⁺ in liquid and frozen solutions is dominated by a 5-coordinate square pyramid, but with split axial CN6 appearing in the frozen glasses. Among these phases, the Cu-O axial distances vary across 1 Å, and the equatorial angles depart significantly from the square plane. Although all these structures remove the d_x²-y² , d_z² degeneracy, no structure can be described as a Jahn-Teller (JT) axially elongated octahedron. The JT-octahedral description for dissolved [Cu(aq)]²⁺ should thus be abandoned in favor of square pyramidal [Cu(H₂O)₅]²⁺. The revised ligand environments have bearing on questions of the Cu(i)/Cu(ii) self-exchange rate and on the mechanism for ligand exchange with bulk water. The plasticity of dissolved Cu(ii) complex ions falsifies the foundational assumption of the rack-induced bonding theory of blue copper proteins and obviates any need for a thermodynamically implausible protein constraint.

Mechanisms of the polyol reduction of copper(ii) salts depending on the anion type and diol chain length

Systematic experiments were carried out to identify the main factors influencing the polyol reduction of copper(ii) compounds to elemental copper.

Copper Oxidation/Reduction in Water and Protein: Studies with DFTB3/MM and VALBOND Molecular Dynamics Simulations

We apply two recently developed computational methods, DFTB3 and VALBOND, to study copper oxidation/reduction processes in solution and protein. The properties of interest include the coordination structure of copper in different oxidation states in water or in a protein (plastocyanin) active site, the reduction potential of the copper ion in different environments, and the environmental response to copper oxidation. The DFTB3/MM and VALBOND simulation results are compared to DFT/MM simulations and experimental results whenever possible. For a copper ion in aqueous solution, DFTB3/MM results are generally close to B3LYP/MM with a medium basis, including both solvation structure and reduction potential for Cu(II); for Cu(I), however, DFTB3/MM finds a two-water coordination, similar to previous Born-Oppenheimer molecular dynamics simulations using BLYP and HSE, whereas B3LYP/MM leads to a tetrahedron coordination. For a tetraammonia copper complex in aqueous solution, VALBOND and DFTB3/MM are consistent in terms of both structural and dynamical properties of solvent near copper for both oxidation states. For copper reduction in plastocyanin, DFTB3/MM simulations capture the key properties of the active site, and the computed reduction potential and reorganization energy are in fair agreement with experiment, especially when the periodic boundary condition is used. Overall, the study supports the value of VALBOND and DFTB3(/MM) for the analysis of fundamental copper redox chemistry in water and protein, and the results also help highlight areas where further improvements in these methods are desirable.

Structural observations of time dependent oscillatory behavior of CuIICl2 solutions measured via extended X-ray absorption fine structure

Cell surface ECTO-NOX proteins exhibit a clock-related, temperature-independent entrainable pattern of periodic (24 min) oscillations in the rate of oxidation of NAD(P)H. Aqueous solutions of copper salts also oxidize NAD(P)H with a similar temperature-independent pattern. For both, five maxima are observed, two of which are separated by 6 min and the remaining three are separated by 4.5 min. In D2O, the pattern is retained but the period length is proportionately increased to 30 min in direct relationship to the 30 h circadian day observed with D2O-grown organisms. With copper solutions, periodic changes in redox potential correlate precisely with the periodic changes in the rates of NAD(P)H oxidation. Consequently, the local environment of the Cu2+ ion in copper chloride solutions was investigated by X-ray absorption spectroscopy. Detailed extended X-ray absorption fine structure (EXAFS) analyses revealed a pattern of oscillations closely resembling those of the copper-catalyzed oxidation of NADH. With CuCl2 in D2O, a pattern with a period length of 30 min was observed. The findings suggest a regular pattern of distortion in the axial and/or equatorial oxygen atoms of the coordinated water molecules which correlate with redox potential changes sufficient to oxidize NADH. A metastable equilibrium condition in the ratio of ortho to para nuclear spin orientation of the water associated hydrogen atoms would be kinetically consistent with a 24-30 min timeframe. The temperature independence of the biological clock can thus be understood as the consequence of a physical rather than a chemical basis for the timing events.

New Insights into the Jahn-Teller Effect through ab initio Quantum-Mechanical/Molecular-Mechanical Molecular Dynamics Simulations of CuII in Water

The CuII hydration shell structure has been studied by means of classical molecular dynamics (MD) simulations including three-body corrections and hybrid quantum-mechanical/molecular-mechanical (QM/MM) molecular dynamics (MD) simulations at the Hartree-Fock level. The copper(II) ion is found to be six-fold coordinated and [Cu(H2O)6]2+ exhibits a distorted octahedral structure. The QM/MM MD approach reproduces correctly the experimentally observed Jahn-Teller effect but exhibits faster inversions (< 200 fs) and a more complex behaviour than expected from experimental data. The dynamic Jahn-Teller effect causes the high lability of [Cu(H2O)6]2+ with a ligand-exchange rate constant some orders or magnitude higher than its neighbouring ions NiII and ZnII. Nevertheless, no first-shell water exchange occurred during a 30-ps simulation. The structure of the hydrated ion is discussed in terms of radial distribution functions, coordination numbers, and various angular distributions and the dynamical properties as librational and vibrational motions and reorientational times were evaluated, which lead to detailed information about the first hydration shell. Second-shell water-exchange processes could be observed within the simulation time scale and yielded a mean ligand residence time of approximelty 20 ps.

First solvation shell of the Cu(II) aqua ion: evidence for fivefold coordination

We determined the structure of the hydrated Cu(II) complex by both neutron diffraction and first-principles molecular dynamics. In contrast with the generally accepted picture, which assumes an octahedrally solvated Cu(II) ion, our experimental and theoretical results favor fivefold coordination. The simulation reveals that the solvated complex undergoes frequent transformations between square pyramidal and trigonal bipyramidal configurations. We argue that this picture is also consistent with experimental data obtained previously by visible near-infrared absorption, x-ray absorption near-edge structure, and nuclear magnetic resonance methods. The preference of the Cu(II) ion for fivefold instead of sixfold coordination, which occurs for other cations of comparable charge and size, results from a Jahn-Teller destabilization of the octahedral complex.