Coordination and ligand exchange dynamics of solvated metal ions

Structures of Th4+ aqueous solutions: insights from AIMD and metadynamics simulations

Solution chemistry of actinide ions is critical to understanding the solvation behaviors and hydrolysis process. Using tetravalent thorium ion Th4+ as a representative example, we investigate the local structures and dynamic behaviors of hydrated Th4+ ions by ab initio molecular dynamics (AIMD) simulations using the recently developed norm-conserving pseudopotentials and basis sets optimized for actinides (J.-B. Lu et al., J. Chem. Theory Comput. 2021, 17, 3360-3371). AIMD simulations reveal two distinct solvation shells, with the first shell comprising 9 water molecules at approximately rTh-O = 2.50 Å and exhibiting a tricapped trigonal prism geometry. These conclusions are confirmed through metadynamics simulations and further structural analysis. AIMD simulations also show the slight effect of temperature and counterions on the structure of the solution. The structured solvation shells of the highly charged Th4+ ion with the specific geometry, distinct from the structure of liquid water, lead to corresponding structural changes in the hydrogen bond network in water. Additionally, beyond the solvent-shared ion pair (SIP) state observed in the unbiased AIMD simulations, the metadynamics simulations reconstruct a two-dimensional free energy surface that clearly indicates the potential stability of the contact ion pair (CIP) state in the system with Cl- as a counterion. The findings in this work provide insights into the solution chemistry of actinides and serve as a reference for studying other actinide solution systems.

Modelling ligand exchange in metal complexes with machine learning potentials

Metal ions are irreplaceable in many areas of chemistry, including (bio)catalysis, self-assembly and charge transfer processes. Yet, modelling their structural and dynamic properties in diverse chemical environments remains challenging for both force fields and ab initio methods. Here, we introduce a strategy to train machine learning potentials (MLPs) using MACE, an equivariant message-passing neural network, for metal-ligand complexes in explicit solvents. We explore the structure and ligand exchange dynamics of Mg2+ in water and Pd2+ in acetonitrile as two illustrative model systems. The trained potentials accurately reproduce equilibrium structures of the complexes in solution, including different coordination numbers and geometries. Furthermore, the MLPs can model structural changes between metal ions and ligands in the first coordination shell, and reproduce the free energy barriers for the corresponding ligand exchange. The strategy presented here provides a computationally efficient approach to model metal ions in solution, paving the way for modelling larger and more diverse metal complexes relevant to biomolecules and supramolecular assemblies.

Solvation energies of the ferrous ion in water and in ammonia at various temperatures

CONTEXT

The solvation of metal ions is crucial to understanding relevant properties in physics, chemistry, or biology. Therefore, we present solvation enthalpies and solvation free energies of the ferrous ion in water and ammonia. Our results agree well with the experimental reports for the hydration free energy and hydration enthalpy. We obtained [Formula: see text] kJ mol[Formula: see text] for the hydration free energy and [Formula: see text] kJ mol[Formula: see text] for the hydration enthalpy of ferrous ion in water at room temperature. At ambient temperature, we obtained [Formula: see text] kJ mol[Formula: see text] as the [Formula: see text] ammoniation free energy and [Formula: see text] kJ mol[Formula: see text] for the ammoniation enthalpy. In addition, the free energy of solvation is deeply affected when the temperature increases. This pattern can be attributed to the rise of entropy when the temperature rises. Besides, the temperature does not affect the ammoniation enthalpies and the hydration enthalpy of the [Formula: see text] ion.

METHOD

All the geometry optimizations are performed at the MP2 methods associated with the 6-31++g(d,p) basis set of Pople. solvated phase structures of [Formula: see text] ion in water or in ammonia are performed using the PCM model. The [Formula: see text] program suite was used to perform all the calculations. The program TEMPO was also used to evaluate the temperature sensitivity of the different obtained geometries.

Mechanistic insight into homogeneous catalytic crosslinking behavior between cellulose and epoxide by explicit solvent models

Inspired by recent advances on functional modification of cellulosic materials, the crosslinking behaviors of epoxide with cellulose under the catalysis of different homogeneous catalysts including H2O, Brønsted acid, Brønsted base, Lewis acid and neutral salt were systematically investigated using density functional theory (DFT) methods with hybrid micro-solvation-continuum approach. The results showed that catalytic activity, reaction mechanism and regioselectivity are determined by the combined effect of catalyst type, electronic effect and steric hindrance. All the homogeneous catalysts have catalytic activity for the crosslinking reaction, which decreases in the order of NaOH > HCl > NCl3 > MCl2 > CH3COOH > NaCl (N = Fe3+, Al3+; M = Zn2+, Ca2+). Upon the catalysis of NaOH, hydroxyl group of cellulose is firstly deprotonated to form a carbanion-like intermediate which will further attack the less sterically hindered C atom of epoxide showing excellent regioselectivity. Acidic catalysts readily cause epoxide protonated, which suffers from nucleophilic attack of cellulose and forms the carbocation-like intermediate. Brønsted acid exhibits poor regioselectivity, however, Lewis acid shows an interesting balance between catalytic activity and regioselectivity for the crosslinking reaction, which may be attributed to the unique catalysis and stabilization effects of its coordinated H2O on the transition state structure.

High-Throughput Aqueous Electrolyte Structure Prediction Using IonSolvR and Equivariant Graph Neural Network Potentials

Neural network potentials have recently emerged as an efficient and accurate tool for accelerating ab initio molecular dynamics (AIMD) in order to simulate complex condensed phases such as electrolyte solutions. Their principal limitation, however, is their requirement for sufficiently large and accurate training sets, which are often composed of Kohn-Sham density functional theory (DFT) calculations. Here we examine the feasibility of using existing density functional tight-binding (DFTB) molecular dynamics trajectory data available in the IonSolvR database in order to accelerate the training of E(3)-equivariant graph neural network potentials. We show that the solvation structure of Na+ and Cl- in aqueous NaCl solutions can be accurately reproduced with remarkably small amounts of data (i.e., 100 MD frames). We further show that these predictions can be systematically improved further via an embarrassingly parallel resampling approach.

Molecular Mechanisms of Mercury-Sensitive Aquaporins

Aquaporins are transmembrane channels that allow for the passive permeation of water and other small molecules across biological membranes. Their channel activities are sensitive to mercury ions. Intriguingly, while most aquaporins are inhibited by mercury ions, several aquaporins are activated by mercury ions. The molecular basis of the opposing aquaporin regulation by mercury remains poorly understood. Herein, we investigated AqpZ inhibition and AQP6 activation upon binding of mercury ions using solid-state NMR (ssNMR) and molecular dynamics (MD) simulations. Based on the structure of the Hg-AqpZ complex constructed by MD simulations and ssNMR, we identified that the pore closure was caused by mercury-induced conformational changes of the key residue R189 in the selectivity filter region, while pore opening was caused by conformational changes of residues H181 and R196 in the selectivity filter region in AQP6. Both conformational changes were caused by the disruption of the H-bond network of R189/R196 by mercury. The molecular details provided a structural basis for mercury-mediated functional changes in aquaporins.

A quantum chemical molecular dynamics repository of solvated ions

The importance of ion-solvent interactions in predicting specific ion effects in contexts ranging from viral activity through to electrolyte viscosity cannot be underestimated. Moreover, investigations of specific ion effects in nonaqueous systems, highly relevant to battery technologies, biochemical systems and colloid science, are severely limited by data deficiency. Here, we report IonSolvR – a collection of more than 3,000 distinct nanosecond-scale ab initio molecular dynamics simulations of ions in aqueous and non-aqueous solvent environments at varying effective concentrations. Density functional tight binding (DFTB) is used to detail the solvation structure of up to 55 solutes in 28 different protic and aprotic solvents. DFTB is a fast quantum chemical method, and as such enables us to bridge the gap between efficient computational scaling and maintaining accuracy, while using an internally-consistent simulation technique. We validate the database against experimental data and provide guidance for accessing individual IonSolvR records.

Measurement(s)	solvation structure
Technology Type(s)	quantum chemistry computational method • Molecular Dynamics

Adsorption and solidification of cadmium by calcium sulfate dihydrate (gypsum) in an aqueous environment: a dispersion-corrected DFT and ab initio molecular dynamics study

Phosphogypsum is a by-product of phosphoric acid production by a wet-process and can be used to prepare adsorption materials to treat Cd(II) in sewage. Its main component is calcium sulfate dihydrate (gypsum). In the present study, incorporation and adsorption of cadmium on gypsum in an aqueous environment have been investigated using dispersion-corrected density functional theory, and the effect of different vacancy defects on the adsorption and doping of Cd on gypsum was also studied. The results show that Cd impurity defects can form more easily in gypsum crystals with Ca vacancy defects. The increase in both calcium vacancy defects and sulfate vacancy defects is beneficial to the incorporation of Cd in calcium sulfate dihydrate. In addition, the calcium vacancy defects on the gypsum surface can promote the adsorption of Cd(II), while sulfate vacancy defects on the gypsum surface inhibit the adsorption of Cd(II). Therefore, increasing the Ca vacancy defects of gypsum will help in improving the adsorption and curing properties of phosphogypsum for cadmium.

In Situ Characterization of Dehydration during Ion Transport in Polymeric Nanochannels

The transport of hydrated ions across nanochannels is central to biological systems and membrane-based applications, yet little is known about their hydrated structure during transport due to the absence of in situ characterization techniques. Herein, we report experimentally resolved ion dehydration during transmembrane transport using modified in situ liquid ToF-SIMS in combination with MD simulations for a mechanistic reasoning. Notably, complete dehydration was not necessary for transport to occur across membranes with sub-nanometer pores. Partial shedding of water molecules from ion solvation shells, observed as a decrease in the average hydration number, allowed the alkali-metal ions studied here (lithium, sodium, and potassium) to permeate membranes with pores smaller than their solvated size. We find that ions generally cannot hold more than two water molecules during this sterically limited transport. In nanopores larger than the size of the solvation shell, we show that ionic mobility governs the ion hydration number distribution. Viscous effects, such as interactions with carboxyl groups inside the membrane, preferentially hinder the transport of the mono- and dihydrates. Our novel technique for studying ion solvation in situ represents a significant technological leap for the nanofluidics field and may enable important advances in ion separation, biosensing, and battery applications.

Photoinduced Ligand Exchange Dynamics of a Polypyridyl Ruthenium Complex in Aqueous Solution

The understanding of photoinduced ligand exchange mechanisms in polypyridyl ruthenium(II) complexes operating in aqueous solution is of crucial importance to rationalize their photoreactivity. Herein, we demonstrate that a synergetic use of ab initio molecular dynamics simulations and static calculations, both conducted at the DFT level, can provide a full understanding of photosubstitution mechanisms of a monodentate ligand by a solvent water molecule in archetypal ruthenium complexes in explicit water. The simulations show that the photoinduced loss of a monodentate ligand generates an unreactive 16-electron species in a hitherto undescribed pentacoordinated triplet excited state that converts, via an easily accessible crossing point, to a reactive 16-electron singlet ground state, which combines with a solvent water molecule to yield the experimentally observed aqua complex in less than 10 ps. This work paves the way for the rational design of novel photoactive metal complexes relevant for biological applications.

Modulation of High-Spin Co(II) in Li/Co-MOFs as Efficient Fenton-like Catalysts

Developing high-performance catalysts toward the Fenton reaction is important for environmental protection and sustainable development, yet it is still challenging. The high-spin states of first-row transition metal atoms with tetrahedral coordination provide a flexible electronic environment to activate the catalyst and elevate its catalytic activity. As a type of material with adjustable structures, metal-organic frameworks (MOFs) are excellent candidate catalysts as they can accurately regulate the coordination configurations of metal ions. In this paper, we investigate and summarize the direct formation of bimetallic carboxylate Li/Co-MOFs with tetrahedral coordination metal centers in a mixed H2O/polar organic solvent system. The induction of Li(I) ions is manifested in the generation of hydroxides during the dissociation of the Co(II) solvation structure to trigger the tetrahedral coordination behavior of Co(II). These Li/Co-MOFs containing high-spin Co(II) centers can serve as highly efficient Fenton-like catalysts for organics. This study provides a promising strategy for rational design of MOF-based catalysts with high-spin metal centers for application in environment governance.

Structures and relative stability of hydrated ferrous ion clusters and temperature effects

Structures of solvated ferrous ion clusters have been investigated in the singlet and quintet spin states of the ferrous ion. Relative stabilities of isomers are also discussed at different temperatures.

Cu2+ in liquid ammonia-The impact of solvent flexibility and electron correlation in ab initio quantum mechanical charge field molecular dynamics

The impact of solvent flexibility and electron correlation on the simulation results of Cu²⁺ in liquid ammonia has been investigated via an ab initio quantum mechanical charge field molecular dynamics (QMCF MD) simulation approach. To achieve this, three different simulation systems were considered in this study, namely Cu²⁺ in rigid and flexible ammonia at Hartree-Fock (HF) level of theory, as well as resolution of identity second order Møller-Plesset (MP2) perturbation theory in the rigid body case. In all cases, a stable octahedral [Cu(NH₃ )₆ ]²⁺ complex subject to dynamic Jahn-Teller distortions without the occurrence of ligand exchange was observed. The Cu²⁺  - NH₃ distance in the first shell agrees well with the experimental and other theoretical data. In all three cases, the structural data shows that the rigid-body ammonia model in conjunction with the HF level of theory provides accurate data for the first solvation shell, while at the same time, the computational demand and thus the achievable simulation time are much more beneficial. The vibrational analysis of the Cu²⁺  - NH₃ interaction yields similar force constants in the three investigated systems indicating that there is no distinct difference on the dynamical properties of the first solvation shell. In addition to the QMCF MD simulations, a number of natural bond orbital (NBO) analyses were carried out, confirming the strong electrostatic character of the Cu²⁺  - NH₃ interaction.

Density functional theory based molecular dynamics study of solution composition effects on the solvation shell of metal ions

We present an ab initio molecular dynamics study of the alkali metal ions Li+, Na+, K+ and Cs+, and of the alkaline earth metal ions Mg2+ and Ca2+ in both pure water and electrolyte solutions containing the counterions Cl- and SO42-. Simulations were conducted using different density functional theory methods (PBE, BLYP and revPBE), with and without the inclusion of dispersion interactions (-D3). Analysis of the ion-water structure and interaction strength, water exchange between the first and second hydration shell, and hydrogen bond network and low-frequency reorientation dynamics around the metal ions have been used to characterise the influence of solution composition on the ionic solvation shell. Counterions affect the properties of the hydration shell not only when they are directly coordinated to the metal ion, but also when they are at the second coordination shell. Chloride ions reduce the sodium hydration shell and expand the calcium hydration shell by stabilizing under-coordinated hydrated Na(H2O)5+ complexes and over-coordinated Ca(H2O)72+. The same behaviour is observed in CaSO4(aq), where Ca2+ and SO42- form almost exclusively solvent-shared ion pairs. Water exchange between the first and second hydration shell around Ca2+ in CaSO4(aq) is drastically decelerated compared with the simulations of the hydrated metal ion (single Ca2+, no counterions). Velocity autocorrelation function analysis, used to probe the strength of the local ion-water interaction, shows a smoother decay of Mg2+ in MgCl2(aq), which is a clear indication of a looser inter-hexahedral vibration in the presence of chloride ions located in the second coordination shell of Mg2+. The hydrogen bond statistics and orientational dynamics in the ionic solvation shell show that the influence on the water-water network cannot only be ascribed to the specific cation-water interaction, but also to the subtle interplay between the level of hydration of the ions, and the interactions between ions, especially those of opposite charge. As many reactive processes involving solvated metal ions occur in environments that are far from pure water but rich in ions, this computational study shows how the solution composition can result in significant differences in behaviour and function of the ionic solvation shell.

DFT-Assisted Spectroscopic Studies on the Coordination of Small Ligands to Palladium: From Isolated Ions to Nanoparticles

A combination of experimental spectroscopies (UV-vis and Fourier-transform infrared) and computational modeling was used to investigate the coordination of small ligands (aminopropanol and propanediol) to Pd species during the metal nanoparticle formation process. Differences emerged between O- (propanediol) and N-containing (aminopropanol) ligands. In particular, a strong interaction between the NH amino group and Pd²⁺ ions could be inferred on the basis of spectroscopic evidences, which was corroborated by theoretical simulations, which confirmed the preferential coordination of aminopropanol through the NH group. This interaction seems to potentially cause the aminopropanol ligand to control the particle shape through a selective blocking of Pd(100) facets, which promote the growth on the Pd(111) facets.

Reconsidering Calcium Dehydration as the Rate-Determining Step in Calcium Mineral Growth

The dehydration of cations is generally accepted as the rate-limiting step in many processes. Molecular dynamics (MD) can be used to investigate the dynamics of water molecules around cations, and two different methods exist to obtain trajectory-based water dehydration frequencies. Here, these two different post-processing methods (direct method versus survival function) have been implemented to obtain calcium dehydration frequencies from a series of trajectories obtained using a range of accepted force fields. None of the method combinations reproduced the commonly accepted experimental water exchange frequency of 10-8.2 s-1. Instead, our results suggest much faster water dynamics, comparable with more accurate ab initio MD simulations and with experimental values obtained using neutron scattering techniques. We obtained the best agreement using the survival function method to characterize the water dynamics, and we show that different method combinations significantly affect the outcome. Our work strongly suggests that the fast water exchange kinetics around the calcium ions is not rate-limiting for reactions involving dissolved/solvated calcium. Our results further suggest that, for alkali and most of the earth alkali metals, mechanistic rate laws for growth, dissolution, and adsorption, which are based on the principle of rate-limiting cation dehydration, need careful reconsideration.

Glass polymorphism and liquid-liquid phase transition in aqueous solutions: experiments and computer simulations

One of the most intriguing anomalies of water is its ability to exist as distinct amorphous ice forms (glass polymorphism or polyamorphism). This resonates well with the possible first-order liquid-liquid phase transition (LLPT) in the supercooled state, where ice is the stable phase. In this Perspective, we review experiments and computer simulations that search for LLPT and polyamorphism in aqueous solutions containing salts and alcohols. Most studies on ionic solutes are devoted to NaCl and LiCl; studies on alcohols have mainly focused on glycerol. Less attention has been paid to protein solutions and hydrophobic solutes, even though they reveal promising avenues. While all solutions show polyamorphism and an LLPT only in dilute, sub-eutectic mixtures, there are differences regarding the nature of the transition. Isocompositional transitions for varying mole fractions are observed in alcohol but not in ionic solutions. This is because water can surround alcohol molecules either in a low- or high-density configuration whereas for ionic solutes, the water ion hydration shell is forced into high-density structures. Consequently, the polyamorphic transition and the LLPT are prevented near the ions, but take place in patches of water within the solutions. We highlight discrepancies and different interpretations within the experimental community as well as the key challenges that need consideration when comparing experiments and simulations. We point out where reinterpretation of past studies helps to draw a unified, consistent picture. In addition to the literature review, we provide original experimental results. A list of eleven open questions that need further consideration is identified.

The Jahn-Teller effect in mixed aqueous solution: the solvation of Cu2+ in 18.6% aqueous ammonia obtained from ab initio quantum mechanical charge field molecular dynamics

The solvation structure and dynamics of Cu2+ in 18.6 % aqueous ammonia have been investigated using an ab initio quantum mechanical charge field molecular dynamics (QMCF MD) simulation approach at the Hartree–Fock (HF) level of theory applying the LANL2DZ ECP and Dunning DZP basis sets for Cu2+, ammonia and water, respectively. During a simulation time of 20 ps, only NH3 molecules are observed within the first solvation shell of Cu2+, resulting in the formation of an octahedral [Cu(NH3)6]2+ complex. While no exchange of these ligands with the second solvation shell are observed along the simulation, the monitoring of the associated N-Ntrans distances highlight the dynamics of the associated Jahn-Teller distortions, showing on average 2 elongated axial (2.19 Å) and 4 equatorial Cu–N bonds (2.39 Å). The observed structural properties are found in excellent agreement with experimental studies. In addition, an NBO analysis was carried out, confirming the strong electrostatic character of the Cu2+–NH3 interaction.

Ab initio molecular dynamics studies of hydroxide coordination of alkaline earth metals and uranyl

Ab initio molecular dynamics (AIMD) simulations of the Mg2+, Ca2+, Sr2+ and UO22+ ions in either a pure aqueous environment or an environment containing two hydroxide ions have been carried out at the density functional level of theory, employing the generalised gradient approximation via the PBE exchange-correlation functional. Calculated mean M-O bond lengths in the first solvation shell of the aquo systems compared very well to existing experimental and computational literature, with bond lengths well within values measured previously and coordination numbers in line with previously calculated values. When applied to systems containing additional hydroxide ions, the methodology revealed increased bond lengths in all systems. Proton transfer events (PTEs) were recorded and were found to be most prevalent in the strontium hydroxide systems, likely due to the low charge density of the ion and the consequent lack of hydroxide coordination. For all alkaline earths, intrashell PTEs which occurred outside of the first solvation shell were most prevalent. Only three PTEs were identified in the entire simulation data of the uranium dihydroxide system, indicating the clear impact of the increased charge density of the hexavalent uranium ion on the strength of metal-oxygen bonds in aqueous solution. Broadly, systems containing more charge dense ions were found to exhibit fewer PTEs than those containing ions of lower charge density.

Structures of solvated ferrous ion clusters in ammonia and spin-crossover at various temperatures

The solvated ferrous ion in ammonia is hexa-coordinated, irrespective of the temperature.

ISOLDE: a physically realistic environment for model building into low-resolution electron-density maps

This paper introduces ISOLDE, a new software package designed to provide an intuitive environment for high-fidelity interactive remodelling/refinement of macromolecular models into electron-density maps. ISOLDE combines interactive molecular-dynamics flexible fitting with modern molecular-graphics visualization and established structural biology libraries to provide an immersive interface wherein the model constantly acts to maintain physically realistic conformations as the user interacts with it by directly tugging atoms with a mouse or haptic interface or applying/removing restraints. In addition, common validation tasks are accelerated and visualized in real time. Using the recently described 3.8 Å resolution cryo-EM structure of the eukaryotic minichromosome maintenance (MCM) helicase complex as a case study, it is demonstrated how ISOLDE can be used alongside other modern refinement tools to avoid common pitfalls of low-resolution modelling and improve the quality of the final model. A detailed analysis of changes between the initial and final model provides a somewhat sobering insight into the dangers of relying on a small number of validation metrics to judge the quality of a low-resolution model.

Real-time observation of the photoionization-induced water rearrangement dynamics in the 5-hydroxyindole-water cluster by time-resolved IR spectroscopy

Solvation plays an essential role in controlling the mechanism and dynamics of chemical reactions in solution. The present study reveals that changes in the local solute-solvent interaction have a great impact on the timescale of solvent rearrangement dynamics. Time-resolved IR spectroscopy has been applied to a hydration rearrangement reaction in the monohydrated 5-hydroxyindole-water cluster induced by photoionization of the solute molecule. The water molecule changes the stable hydration site from the indolic NH site to the substituent OH site, both of which provide a strongly attractive potential for hydration. The rearrangement time constant amounts to 8 ± 2 ns, and is further slowed down by a factor of more than five at lower excess energy. These rearrangement times are slower by about three orders of magnitude than those reported for related systems where the water molecule is repelled from a repulsive part of the interaction potential toward an attractive well. The excess energy dependence of the time constant is well reproduced by RRKM theory. Differences in the reaction mechanism are discussed on the basis of energy relaxation dynamics.

Towards a dissociative SPC-like water model II. The impact of Lennard-Jones and Buckingham non-coulombic forces

The dissociative water potential by Garofalini and coworkers has been re-formulated in the framework of the widely employed Lennard-Jones and Buckingham potentials, enhancing the transferability of the model to third party simulation programs.

Square planar or octahedral after all? The indistinct solvation of platinum(ii)

The solvation structures of Pd(ii) and Pt(ii) are typically reduced to the well-known square-planar structural motif, although it has been shown, in both experimental and theoretical investigations, that these solutes demonstrate the affinity to bind ligand molecules at elongated distance in axial coordination sites.

Nanosecond Dynamics at Protein Metal Sites: An Application of Perturbed Angular Correlation (PAC) of γ-Rays Spectroscopy

Metalloproteins are essential to numerous reactions in nature, and constitute approximately one-third of all known proteins. Molecular dynamics of proteins has been elucidated with great success both by experimental and theoretical methods, revealing atomic level details of function involving the organic constituents on a broad spectrum of time scales. However, the characterization of dynamics at biomolecular metal sites on nanosecond time scales is scarce in the literature. The aqua ions of many biologically relevant metal ions exhibit exchange of water molecules on the nanosecond time scale or faster, often defining their reactivity in aqueous solution, and this is presumably also a relevant time scale for the making and breaking of coordination bonds between metal ions and ligands at protein metal sites. Ligand exchange dynamics is critical for a variety of elementary steps of reactions in metallobiochemistry, for example, association and dissociation of metal bound water, association of substrate and dissociation of product in the catalytic cycle of metalloenzymes, at regulatory metal sites which require binding and dissociation of metal ions, as well as in the transport of metal ions across cell membranes or between proteins involved in metal ion homeostasis. In Perturbed Angular Correlation of γ-rays (PAC) spectroscopy, the correlation in time and space of two γ-rays emitted successively in a nuclear decay is recorded, reflecting the hyperfine interactions of the PAC probe nucleus with the surroundings. This allows for characterization of molecular and electronic structure as well as nanosecond dynamics at the PAC probe binding site. Herein, selected examples describing the application of PAC spectroscopy in probing the dynamics at protein metal sites are presented, including (1) exchange of Cd²⁺ bound water in de novo designed synthetic proteins, and the effect of remote mutations on metal site dynamics; (2) dynamics at the β-lactamase active site, where the metal ion appears to jump between the two adjacent sites; (3) structural relaxation in small blue copper proteins upon ¹¹¹Ag⁺ to ¹¹¹Cd²⁺ transformation in radioactive nuclear decay; (4) metal ion transfer between two HAH1 proteins with change in coordination number; and (5) metal ion sensor proteins with two coexisting metal site structures. With this Account, we hope to make our modest contribution to the field and perhaps spur additional interest in dynamics at protein metal sites, which we consider to be severely underexplored. Relatively little is known about detailed atomic motions at metal sites, for example, how ligand exchange processes affect protein function, and how the amino acid composition of the protein may control this facet of metal site characteristics. We also aim to provide the reader with a qualitative impression of the possibilities offered by PAC spectroscopy in bioinorganic chemistry, especially when elucidating dynamics at protein metal sites, and finally present data that may serve as benchmarks on a relevant time scale for development and tests of theoretical molecular dynamics methods applied to biomolecular metal sites.

Palladium(ii) in liquid ammonia: an investigation of structural and dynamical properties by applying quantum mechanical charge field molecular dynamics (QMCF-MD)

The structural and dynamical properties of Pd²⁺ in liquid ammonia have been investigated via quantum mechanical charge field molecular dynamics.

Towards a dissociative SPC-like water model – probing the impact of intramolecular Coulombic contributions

A modification of the dissociative Garofalini water model towards an SPC-like Coulombic formulation proved to enhance accuracy and transferability of this successful force field approach.

Modelling how incorporation of divalent cations affects calcite wettability-implications for biomineralisation and oil recovery

Using density functional theory and geochemical speciation modelling, we predicted how solid-fluid interfacial energy is changed, when divalent cations substitute into a calcite surface. The effect on wettability can be dramatic. Trace metal uptake can impact organic compound adsorption, with effects for example, on the ability of organisms to control crystal growth and our ability to predict the wettability of pore surfaces. Wettability influences how easily an organic phase can be removed from a surface, either organic compounds from contaminated soil or crude oil from a reservoir. In our simulations, transition metals substituted exothermically into calcite and more favourably into sites at the surface than in the bulk, meaning that surface properties are more strongly affected than results from bulk experiments imply. As a result of divalent cation substitution, calcite-fluid interfacial energy is significantly altered, enough to change macroscopic contact angle by tens of degrees. Substitution of Sr, Ba and Pb makes surfaces more hydrophobic. With substitution of Mg and the transition metals, calcite becomes more hydrophilic, weakening organic compound adsorption. For biomineralisation, this provides a switch for turning on and off the activity of organic crystal growth inhibitors, thereby controlling the shape of the associated mineral phase.

Equilibrium between 5- and 6-Fold Coordination in the First Hydration Shell of Cu(II)

The hydration structure dynamics of Cu(II) ion is characterized by a combination of classical molecular dynamics simulation and X-ray absorption near-edge spectroscopy. Previous experimental data have been analyzed on the basis of 5- or 6-fold first hydration structure, with a quite well-established equatorial structure. This 4-fold equatorial geometry has been our starting point to develop a simple but effective in silico model, which provides ab initio theoretical X-ray absorption spectra in very good agreement with the experimental data. Our results point out two equally populated 6- and 5-fold hydration structures with remarkable different water residence times of 5 and 98 ps, respectively, and a low free energy barrier between first and second hydration shell.

Improved radiolabeling of DOTATOC with trivalent radiometals for clinical application by addition of ethanol

Background

Typically, metal-based radiopharmaceuticals are synthesized in aqueous solutions with no or low ethanol content. Labeling yields are defined by temperature, period of labeling, amount of precursor, pH etc. As recently observed, radiolabeling yields (RCY) seem to increase in the presence of non-aqueous solvents. Consequently, this effect was investigated systematically using ethanol as non-aqueous solvent (n-as), which is widely utilized in medicine, and DOTATOC as model compound.

Methods

To determine the impact of ethanol on the radiolabeling efficacy, "standard" labeling conditions of ⁶⁸Ga-DOTATOC (95-100 °C, 10-15 min, 20-50 μg DOTATOC, aqueous solution), i.e. 10 nmol (2.9 μM, 14.2 μg), were modified in terms of lower temperature (70 °C) to achieve lower RCY (<75 %). From those lower RCY, positive effects of increasing amounts of ethanol (0-40 vol%) could directly be observed. Labeling parameters were finally evaluated in terms of shorter reaction time and lower amount of precursor. To investigate whether the effects observed are also true for other trivalent radiometals, labeling was also performed with ⁴⁴Sc.

Results

For increasing amounts of ethanol, ⁶⁸Ga-DOTATOC RCY at 70 °C improved significantly. RCY of ~95 % can be achieved within 10 min using 30 vol% ethanol compared to 46 % in the pure aqueous system. If "standard" temperatures of 95 °C are applied, high RCY of 89 % can be achieved within 5 min with much lower amounts of precursor, i.e. even at 0.93 nmol (0.3 μM, 1.3 μg). This also represents significantly increased specific activities. Similar behavior was observed for ⁴⁴Sc where RCY increase successively with increasing amounts of ethanol.

Conclusion

There is clear experimental evidence, that adding more than 20 vol% ethanol to the reaction mixtures significantly improve labeling efficacies. This could be demonstrated for ⁶⁸Ga-DOTATOC and ⁴⁴Sc-DOTATOC in terms of temperature, time and concentration of required precursor. Whether this is a principal phenomenon with practical impact on the radiopharmaceutical chemistry of trivalent metals and whether this applies to other non-aqueous solvents as well - and what the physico-chemical reasons are, remains to be studied in more detail. Nevertheless, the effect observed here will improve ⁶⁸Ga-DOTATOC labeling and may save at least half of the usually applied amount of precursor.

3D-RISM-MP2 Approach to Hydration Structure of Pt(II) and Pd(II) Complexes: Unusual H-Ahead Mode vs Usual O-Ahead One

Solvation of transition metal complexes with water has been one of the fundamental topics in physical and coordination chemistry. In particular, Pt(II) complexes have recently attracted considerable interest for their relation to anticancer activity in cisplatin and its analogues, yet the interaction of the water molecule and the metal center has been obscured. The challenge from a theoretical perspective remains that both the microscopic solvation effect and the dynamical electron correlation (DEC) effect have to be treated simultaneously in a reasonable manner. In this work we derive the analytical gradient for the three-dimensional reference interaction site model Møller-Plesset second order (3D-RISM-MP2) free energy. On the basis of the three-regions 3D-RISM self-consistent field (SCF) method recently proposed by us, we apply a new layer of the Z-vector method to the CP-RISM equation as well as point-charge approximation to the derivatives with respect to the density matrix elements in the RISM-CPHF equation to remarkably reduce the computational cost. This method is applied to study the interaction of H2O with the d(8) square planar transition metal complexes in aqueous solution, trans-[Pt(II)Cl2(NH3)(glycine)] (1a), [Pt(II)(NH3)4](2+) (1b), [Pt(II)(CN)4](2-) (1c), and their Pd(II) analogues 2a, 2b, and 2c, respectively, to elucidate whether the usual H2O interaction through O atom (O-ahead mode) or unusual one through H atom (H-ahead mode) is stable in these complexes. We find that the interaction energy of the coordinating water and the transition metal complex changes little when switching from gas to aqueous phase, but the solvation free energy differs remarkably between the two interaction modes, thereby affecting the relative stability of the H-ahead and O-ahead modes. Particularly, in contrast to the expectation that the O-ahead mode is preferred due to the presence of positive charges in 1b, the H-ahead mode is also found to be more stable. The O-ahead mode is found to be more stable than the H-ahead one only in 2b. The energy decomposition analysis (EDA) at the 3D-RISM-MP2 level revealed that the O-ahead mode is stabilized by the electrostatic (ES) interaction, whereas the H-ahead one is mainly stabilized by the DEC effect. The ES interaction is also responsible for the difference between the Pd(II) and Pt(II) complexes; because the electrostatic potential is more negative along the z-axis in the Pt(II) complex than in the Pd(II) one, the O-ahead mode prefers the Pd(II) complexes, whereas the H-ahead becomes predominant in the Pt(II) complexes.

Mechanism of Formation of Copper(II) Chloro Complexes Revealed by Transient Absorption Spectroscopy and DFT/TDDFT Calculations

Copper(II) complexes are extremely labile with typical ligand exchange rate constants on the order of 10(6)-10(9) M(-1) s(-1). As a result, it is often difficult to identify the actual formation mechanism of these complexes. In this work, using UV-vis transient absorption when probing in a broad time range (20 ps to 8 μs) in conjunction with DFT/TDDFT calculations, we studied the dynamics and underlying reaction mechanisms of the formation of extremely labile copper(II) CuCl4(2-) chloro complexes from copper(II) CuCl3(-) trichloro complexes and chloride ions. These two species, produced via photochemical dissociation of CuCl4(2-) upon 420 nm excitation into the ligand-to-metal-charge-transfer electronic state, are found to recombine into parent complexes with bimolecular rate constants of (9.0 ± 0.1) × 10(7) and (5.3 ± 0.4) × 10(8) M(-1) s(-1) in acetonitrile and dichloromethane, respectively. In dichloromethane, recombination occurs via a simple one-step addition. In acetonitrile, where [CuCl3](-) reacts with the solvent to form a [CuCl3CH3CN](-) complex in less than 20 ps, recombination takes place via ligand exchange described by the associative interchange mechanism that involves a [CuCl4CH3CN](2-) intermediate. In both solvents, the recombination reaction is potential energy controlled.

A high-resolution XAS study of aqueous Cu(II) in liquid and frozen solutions: pyramidal, polymorphic, and non-centrosymmetric

High-resolution EXAFS (k = 18 Å(-1)) and MXAN XAS analyses show that axially elongated square pyramidal [Cu(H2O)5](2+) dominates the structure of Cu(II) in aqueous solution, rather than 6-coordinate JT-octahedral [Cu(H2O)6](2+). Freezing produced a shoulder at 8989.6 eV on the rising XAS edge and an altered EXAFS spectrum, while 1s → 3d transitions remained invariant in energy position and intensity. Core square pyramidal [Cu(H2O)5](2+) also dominates frozen solution. Solvation shells were found at ∼3.6 Å (EXAFS) or ∼3.8 Å (MXAN) in both liquid and frozen phases. However, MXAN analysis revealed that about half the time in liquid solution, [Cu(H2O)5](2+) associates with an axially non-bonding 2.9 Å water molecule. This distant water apparently organizes the solvation shell. When the 2.9 Å water molecule is absent, the second shell is undetectable to MXAN. The two structural arrangements may represent energetic minima of fluxional dissolved aqueous [Cu(H2O)5](2+). The 2.9 Å trans-axial water resolves an apparent conflict of the [Cu(H2O)5](2+) core model with a dissociational exchange mechanism. In frozen solution, [Cu(H2O)5](2+) is associated with either a 3.0 Å axial non-bonded water molecule or an axial ClO4(-) at 3.2 Å. Both structures are again of approximately equal presence. When the axial ClO4(-) is present, Cu(II) is ∼0.5 Å above the mean O4 plane. This study establishes [Cu(H2O)5](2+) as the dominant core structure for Cu(II) in water solution, and is the first to both empirically resolve multiple extended solution structures for fluxional [Cu(H2O)5](2+) and to provide direct evidence for second shell dynamics.

Structure and water exchange dynamics of hydrated oxo halo ions in aqueous solution using QMCF MD simulation, large angle X-ray scattering and EXAFS

Theoretical ab initio quantum mechanical charge field molecular dynamics (QMCF MD) has been applied in conjunction with experimental large angle X-ray scattering (LAXS) and EXAFS measurements to study structure and dynamics of the hydrated oxo chloro anions chlorite, ClO2(-), chlorate, ClO3(-), and perchlorate, ClO4(-). In addition, the structures of the hydrated hypochlorite, ClO(-), bromate, BrO3(-), iodate, IO3(-) and metaperiodate, IO4(-), ions have been determined in aqueous solution by means of LAXS. The structures of the bromate, metaperiodate, and orthoperiodate, H2IO6(3-), ions have been determined by EXAFS as solid sodium salts and in aqueous solution as well. The results show clearly that the only form of periodate present in aqueous solution is metaperiodate. The Cl-O bond distances in the hydrated oxo chloro anions as determined by LAXS and obtained in the QMCF MD simulations are in excellent agreement, being 0.01-0.02 Å longer than in solid anhydrous salts due to hydration through hydrogen bonding to water molecules. The oxo halo anions, all with unit negative charge, have low charge density making them typical structure breakers, thus the hydrogen bonds formed to the hydrating water molecules are weaker and more short-lived than those between water molecules in pure water. The water exchange mechanism of the oxo chloro anions resembles those of the oxo sulfur anions with a direct exchange at the oxygen atoms for perchlorate and sulfate. The water exchange rate for the perchlorate ion is significantly faster, τ0.5 = 1.4 ps, compared to the hydrated sulfate ion and pure water, τ0.5 = 2.6 and 1.7 ps, respectively. The angular radial distribution functions show that the chlorate and sulfite ions have a more complex water exchange mechanism. As the chlorite and chlorate ions are more weakly hydrated than the sulfite ion the spatial occupancy is less well-defined and it is not possible to follow any well-defined migration pattern as it is difficult to distinguish between hydrating water molecules and bulk water in the region close to the ions.