Predicting the stability constants of metal-ion complexes from first principles

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.

Comprehensive Mechanistic View of the Hydrolysis of Oxadiazole-Based Inhibitors by Histone Deacetylase 6 (HDAC6)

Histone deacetylase (HDAC) inhibitors used in the clinic typically contain a hydroxamate zinc-binding group (ZBG). However, more recent work has shown that the use of alternative ZBGs, and, in particular, the heterocyclic oxadiazoles, can confer higher isoenzyme selectivity and more favorable ADMET profiles. Herein, we report on the synthesis and biochemical, crystallographic, and computational characterization of a series of oxadiazole-based inhibitors selectively targeting the HDAC6 isoform. Surprisingly, but in line with a very recent finding reported in the literature, a crystal structure of the HDAC6/inhibitor complex revealed that hydrolysis of the oxadiazole ring transforms the parent oxadiazole into an acylhydrazide through a sequence of two hydrolytic steps. An identical cleavage pattern was also observed both in vitro using the purified HDAC6 enzyme as well as in cellular systems. By employing advanced quantum and molecular mechanics (QM/MM) and QM calculations, we elucidated the mechanistic details of the two hydrolytic steps to obtain a comprehensive mechanistic view of the double hydrolysis of the oxadiazole ring. This was achieved by fully characterizing the reaction coordinate, including identification of the structures of all intermediates and transition states, together with calculations of their respective activation (free) energies. In addition, we ruled out several (intuitively) competing pathways. The computed data (ΔG^‡ ≈ 21 kcal·mol^-1 for the rate-determining step of the overall dual hydrolysis) are in very good agreement with the experimentally determined rate constants, which a posteriori supports the proposed reaction mechanism. We also clearly (and quantitatively) explain the role of the -CF₃ or -CHF₂ substituent on the oxadiazole ring, which is a prerequisite for hydrolysis to occur. Overall, our data provide compelling evidence that the oxadiazole warheads can be efficiently transformed within the active sites of target metallohydrolases to afford reaction products possessing distinct selectivity and inhibition profiles.

Redox Potentials with COSMO-RS: Systematic Benchmarking with Different Databases

Recent techniques of computational electrochemistry can yield redox potentials with accuracy as good as 0.1 V. Yet, many such methods are not universal, easy to use, or computationally efficient. Herein, we provide a systematic benchmarking of a relatively cheap and straightforward computational approach for fairly accurate computations of redox potentials. It is based on a combination of the conductor-like screening model for real solvents (COSMO-RS) and the density functional theory (DFT). The benchmarking is done with databases covering diverse redox systems, including transition-metal aquacomplexes and various organic and inorganic compounds. In addition, we also present our own test set aiming at maximum chemical diversity and maximum range of redox potential values. We assess the performance of the fairly efficient computational protocol combining the COSMO-RS with the BP86 DFT functional. This is done by calibrating it against different high-level state-of-the-art techniques, in particular, polarizable continuum model (PCM) connected to composite G3(MP2,CC)(+) method, domain-based pair natural orbital implementation of coupled cluster theory, or complete basis set CBS-QB3 method. We also elaborate on the absolute reduction potential value of standard hydrogen electrode to be used with COSMO-RS, and we propose the value of approx. 4.4 V. The COSMO-RS/BP86-D3(BJ) combination outperforms other methods on a wide range of redox systems. However, we show that its accuracy is based on a balanced error cancelation and, therefore, it cannot be further systematically improved. As a result, the proposed procedure represents a pragmatic choice for large-scale screening, yet it could be combined with more advanced methods for detailed studies.

The role of tridentate ligands on the redox stability of anticancer gold(III) complexes

Gold(III) complexes are promising compounds for cancer chemotherapy, whose action depends on their redox stability. In this context, the choice of ligands is crucial to adjust their reactivity and biological response. The present study addressed the effect of the gold coordination sphere on the reduction potential (E^o) for ten gold(III) complexes containing five or six-membered rings tridentate ligands - [Au^III(trident)Cl]³⁺ⁿ (trident = N^N^N, C^N^N, C^C^N, C^N^C, and N^C^N). The calculated E^o covered a broad range of 2500 mV with the most stable complexes containing two AuC bonds (E^o = -1.85 V for [Au^III(C^C^N)Cl] - f). For complexes with one AuC bond, the N^C^N ligands stabilize the gold(III) complex more efficiently than N^N^C; however, the inclusion of the non-innocent ligand bipy (2,2'-bipyridine) in N^N portion provides an extra stabilization effect. Among the derivatives with one AuC bond, [Au^III(N^N^C)Cl]⁺ (N^N = bipy) (a) showed E^o = -1.20 V. For the complexes with N^N^N ligands, E^o was positive and almost constant (+0.60 V). Furthermore, the kinetics for ligand exchange reactions (Cl^-/H₂O, H₂O/Cys and Cl^-/Cys) were monitored for the most stable compounds and the energy profiles compared to the reduction pathways.

Two Ligands of Interest in Recovering Uranium from the Oceans: The Correct Formation Constants of the Uranyl(VI) Cation with 2,2'-Bipyridyl-6,6'-dicarboxylic Acid and 1,10-Phenanthroline-2,9-dicarboxylic Acid

The ligands BDA (2,2'-bipyridyl-6,6'-dicarboxylic acid) and PDA (1,10-phenanthroline-2,9-dicarboxylic acid) are of interest as functional group types for ion-exchange materials for extracting uranium from the oceans, reported in a previous paper for PDA Lashley, M. A. ( Inorg. Chem. 2016 55 10818 10829). Yang, Y. ( Inorg. Chem. 2019, 58, 6064 6074) have published what they claim to be a more accurate result for the formation of the UO₂²⁺/PDA complex of log K₁ = 22.84 compared with our reported value of log K₁ = 16.5, as well as log K₁ = 21.52 for the BDA complex. The determination of log K₁ for the PDA and BDA complexes with the UO₂²⁺ cation was carried out by Yang et al. using a competition reaction between DTPA (diethylenetriamine pentaacetic acid) and BDA or PDA, monitoring the absorbance due to the BDA and PDA ligands. This competition method using absorbance versus pH titrations was developed for determining the formation constants of the complexes of several polypyridyl ligands plus PDA complexes of metal ions, which were too stable for log K determination by competition with protons. A key feature of such titrations is that in the competition reaction, the displacement of the pyridyl donor ligand (e.g., PDA) by the competing ligand (e.g., DTPA), the absorbance spectrum of the displaced pyridyl donor ligand should be observed. Competing ligands used to date have been EDTA (ethylenediamine tetraacetic acid), DTPA, or the hydroxide ion. In the study of Yang et al., no such displaced PDA or BDA was apparent in the absorbance spectra in their titrations so that their reported log K₁ values have no validity. Their log K₁ values are so much higher than log K₁ for the uranyl DTPA complex (∼13.6) that DTPA could not possibly displace BDA or PDA from the uranyl cation, and a competition reaction could not possibly occur. We report the correct value of log K₁ = 15.4 (ionic strength = zero) for the uranyl BDA complex, to illustrate the correct determination of such a constant by a competition reaction between BDA and hydroxide, showing how the characteristic absorbance spectrum for a BDA complex, here the UO₂²⁺ complex, disappears, and the distinctive absorbance spectrum of the free nonprotonated BDA ligand appears as the pH is increased, and BDA is displaced by the hydroxide ion.

Cyclic Octapeptides Composed of Two Glutathione Units Outperform the Monomer in Lead Detoxification

A rationally-designed scaffold of cyclic octapeptides composed of two units of the natural tripeptide glutathione (GSH) was optimized to strongly and selectively capture toxic lead ions (Pb(II)). Using state-of-the-art computational tools, a list of eleven plausible peptides was shortened to five analogs based on their calculated affinity to Pb(II) ions. We then synthesized and investigated them for their abilities to recover Pb-poisoned human cells. A clear pattern was observed from the in vitro detoxification results, indicating the importance of cavity size and polar moieties to enhance metal capturing. These, together with the apparent benefit in cyclizing the peptides, improved the detoxification of the two lead peptides by approximately two folds compared to GSH and the benchmark chelating agents against Pb poisoning. Moreover, the two peptides did not show any toxicity and, therefore, were thoroughly investigated to determine their potential as next-generation remedies for Pb poisoning.

Thiolation and Carboxylation of Glutathione Synergistically Enhance Its Lead-Detoxification Capabilities

The natural tripeptide glutathione (GSH) is a ubiquitous compound harboring various biological tasks, among them interacting with essential and toxic metal ions. Yet, although weakly binding the poisonous metal lead (Pb), GSH poorly detoxifies it. β-Mercaptoaspartic acid is a new-to-nature novel amino acid that was found to enhance the Pb-detoxification capability of a synthetic cyclic tetrapeptide. Aiming to explore the advantages of noncanonical amino acids (ncAAs) of this nature, we studied the detoxification capabilities of GSH and three analogue peptides, each of which contains at least one ncAA that harbors both free carboxylate and thiolate groups. A thorough investigation that includes in vitro detoxification and mechanistic evaluations, metal-binding affinity, metal selectivity, and computational studies shows that these ncAAs are highly beneficial in additively enhancing Pb binding and reveals the importance of both high affinity and metal selectivity in synergistically reducing Pb toxicity in cells. Hence, such ncAAs join the chemical toolbox against Pb poisoning and pollution, enabling peptides to strongly and selectively bind the toxic metal ion.

Potent Cyclic Tetrapeptide for Lead Detoxification

Lead (Pb) is a ubiquitous poisonous metal, affecting the health of vast populations worldwide. Medications to treat Pb poisoning suffer from various limitations and are often toxic owing to insufficient metal selectivity. Here, we report a cyclic tetrapeptide that selectively binds Pb and eradicates its toxic effect on the cellular level, with superior potency than state-of-the-art drugs. The Pb-peptide complex is remarkably strong and was characterized experimentally and computationally. Accompanied by the lack of toxicity and enhanced stability of this peptide, these qualities indicate its merit as a potential remedy for Pb poisoning.

Computational predictions of metal-macrocycle stability constants require accurate treatments of local solvent and pH effects

Rational design of molecular chelating agents requires a detailed understanding of physicochemical ligand-metal interactions in solvent phase. Computational quantum chemistry methods should be able to provide this, but computational reports have shown poor accuracy when determining absolute binding constants for many chelating molecules. To understand why, we compare and benchmark static- and dynamics-based computational procedures for a range of monovalent and divalent cations binding to a conventional cryptand molecule: 2.2.2-cryptand ([2.2.2]). The benchmarking comparison shows that dynamics simulations using standard OPLS-AA classical potentials can reasonably predict binding constants for monovalent cations, but these procedures fail for divalent cations. We also consider computationally efficient static procedure using Kohn-Sham density functional theory (DFT) and cluster-continuum modeling that accounts for local microsolvation and pH effects. This approach accurately predicts binding energies for monovalent and divalent cations with an average error of 3.2 kcal mol-1 compared to experiment. This static procedure thus should be useful for future molecular screening efforts, and high absolute errors in the literature may be due to inadequate modeling of local solvent and pH effects.

Phosphonate Chelators for Medicinal Metal Ions

A family of phosphonate-bearing chelators was synthesized to study their potential in metal-based (radio)pharmaceuticals. Three ligands (H₆phospa, H₆dipedpa, H₆eppy; structures illustrated in manuscript) were fully characterized, including X-ray crystallographic structures of H₆phospa and H₆dipedpa. NMR spectroscopy techniques were used to confirm the complexation of each ligand with selected trivalent metal ions. These methods were particularly useful in discerning structural information for Sc³⁺ and La³⁺ complexes. Solution studies were conducted to evaluate the complex stability of 15 metal complexes. As a general trend, H₆phospa was noted to form the most stable complexes, and H₆eppy associated with the least stable complexes. Moreover, In³⁺ complexes were determined to be the most stable, and complexes with La³⁺ were the least stable, across all metals. Density functional theory (DFT) was employed to calculate structures of H₆phospa and H₆dipedpa complexes with La³⁺ and Sc³⁺. A comparison of experimental ¹H NMR spectra with calculated ¹H NMR spectra using DFT-optimized structures was used as a method of structure validation. It was noted that theoretical NMR spectra were very sensitive to a number of variables, such as ligand configuration, protonation state, and the number/orientation of explicit water molecules. In general, the inclusion of an explicit second shell of water molecules qualitatively improved the agreement between theoretical and experimental NMR spectra versus a polarizable continuum solvent model alone. Formation constants were also calculated from DFT results using potential-energy optimized structures. Strong dependence of molecular free energies on explicit water molecule number, water molecule configuration, and protonation state was observed, highlighting the need for dynamic data in accurate first-principles calculations of metal-ligand stability constants.

Synthesis, Characterisation and In Vitro Anticancer Activity of Catalytically Active Indole-Based Half-Sandwich Complexes

The synthesis, characterisation and evaluation of the in vitro cytotoxicity of four indole-based half-sandwich metal complexes towards two ovarian cancer cell lines (A2780 and A2780cisR) and one normal prostate cell line (PNT2) are presented herein. Although capable of inducing catalytic oxidation of NADH and able to reduce NAD+ with high turnover frequencies, in cells and in the presence of sodium formate, these complexes also strongly interact with biomolecules such as glutathione. This work highlights that efficient out-of-cells catalytic activity might lead to higher reactivity towards biomolecules, thus inhibiting the in-cells catalytic processes.

A coumarin derivative-Cu2+ complex-based fluorescent chemosensor for detection of biothiols

Herein, a novel fluorescent sensor has been developed for the detection of biothiols based on theoretical calculations of the stability constant of the complex between a Cu²⁺ ion and (E)-3-((2-(benzo[d]thiazol-2-yl)hydrazono)methyl)-7-(diethylamino) coumarin (BDC) as a fluorescent ligand. In this study, on the basis of density functional theory method, the Gibbs free energy of ligand-exchange reaction and the solvation model were carried out using thermodynamic cycles. The obtained results are in good agreement with the experimental data. The BDC–Cu²⁺ complex can be used as a fluorescent sensor for the detection of biothiols in the presence of non-thiol containing amino acids, with a detection limit for cysteine at 0.3 μM. Moreover, theoretical calculations of excited states were used to elucidate variations in the fluorescence properties. The computed results show that the excited doublet states D₂ and D₁ are dark doublet states, which quench the fluorescence of the complex.

This study opens a new research direction towards the utilization of complexes between metal ions and fluorescent ligands for the detection of biothiols based on the theoretical calculations of stability constants.

Predicting the Thermodynamic Stability of Zirconium Radiotracers

The thermodynamic stability of a metal-ligand complex, as measured by the formation constant (log β), is one of the most important parameters that determines metal ion selectivity and potential applications in, for example, radiopharmaceutical science. The stable coordination chemistry of radioactive ⁸⁹Zr⁴⁺ in an aqueous environment is of paramount importance when developing positron-emitting radiotracers based on proteins (usually antibodies) for use with positron emission tomography. Desferrioxamine B (DFO) remains the chelate of choice for clinical applications of ⁸⁹Zr-labeled proteins, but the coordination of DFO to Zr⁴⁺ ions is suboptimal. Many alternative ligands have been reported, but the challenges in measuring very high log β values with metal ions such as Zr⁴⁺ that tend to hydrolyze mean that accurate thermodynamic data are scarce. In this work, density functional theory (DFT) calculations were used to predict the reaction energetics for metal ion complexation. Computed values of pseudoformation constants (log β') are correlated with experimental data and showed an excellent linear relationship (R² = 0.97). The model was then used to estimate the absolute and relative formation constants of 23 different Zr⁴⁺ complexes using a total of 17 different ligands, including many of the alternative bifunctional chelates that have been reported recently for use in ⁸⁹Zr⁴⁺ radiochemistry. In addition, detailed computational studies were performed on the geometric isomerism and hydration state of Zr-desferrioxamine. Collectively, the results offer new insights into Zr⁴⁺ coordination chemistry that will help guide the synthesis of future ligands. The computational model developed here is straightforward and reproducible and can be readily applied in the design of other metal coordination compounds.

Predicting Stability Constants for Terbium(III) Complexes with Dipicolinic Acid and 4-Substituted Dipicolinic Acid Analogues using Density Functional Theory

The relative stability constants of Tb(III) complexes exhibiting binding to a series of 4-substituted analogues of dipicolinic acid (2,6-pyridinedicarboxylic acid) (DPA) were calculated using density functional theory (DFT) with the standard thermodynamic cycle. DFT calculations showed that the strengths of the stability constants were modified by the substituents in the following (decreasing) order: -NH2 > -OH ∼ -CH2OH > -imidazole ∼ -Cl ∼ -Br ∼ -H > -F > -I, with the differences among them falling within one to two log units except for -NH2. Through population and structural analysis, we observed that the -NH2, -OH, -CH2OH, and halide substituents can donate electrons via resonance effect to the pyridine ring of DPA while inductively withdrawing electrons with different strengths, thus resulting in the different binding strengths of the 4-substituted DPAs to the Tb(III) ions. We believe that these observations possess utility not only in the ongoing development of luminescent probes for bioanalytical studies but also for more recent cross-industrial efforts to enhance reservoir surveillance capabilities using chemical tracers within the oil and gas sector.

A new approach to biomining: Bioengineering surfaces for metal recovery from aqueous solutions

Electronics waste production has been fueled by economic growth and the demand for faster, more efficient consumer electronics. The glass and metals in end-of-life electronics components can be reused or recycled; however, conventional extraction methods rely on energy-intensive processes that are inefficient when applied to recycling e-waste that contains mixed materials and small amounts of metals. To make e-waste recycling economically viable and competitive with obtaining raw materials, recovery methods that lower the cost of metal reclamation and minimize environmental impact need to be developed. Microbial surface adsorption can aid in metal recovery with lower costs and energy requirements than traditional metal-extraction approaches. We introduce a novel method for metal recovery by utilizing metal-binding peptides to functionalize fungal mycelia and enhance metal recovery from aqueous solutions such as those found in bioremediation or biomining processes. Using copper-binding as a proof-of-concept, we compared binding parameters between natural motifs and those derived in silico, and found comparable binding affinity and specificity for Cu. We then combined metal-binding peptides with chitin-binding domains to functionalize a mycelium-based filter to enhance metal recovery from a Cu-rich solution. This finding suggests that engineered peptides could be used to functionalize biological surfaces to recover metals of economic interest and allow for metal recovery from metal-rich effluent with a low environmental footprint, at ambient temperatures, and under circumneutral pH.

The AQUA-MER databases and aqueous speciation server: A web resource for multiscale modeling of mercury speciation

To assess the chemical reactivity, toxicity, and mobility of pollutants in the environment, knowledge of their species distributions is critical. Because their direct measurement is often infeasible, speciation modeling is widely adopted. Mercury (Hg) is a representative pollutant for which study of its speciation benefits from modeling. However, Hg speciation modeling is often hindered by a lack of reliable thermodynamic constants. Although computational chemistry (e.g., density functional theory [DFT]) can generate these constants, methods for directly coupling DFT and speciation modeling are not available. Here, we combine computational chemistry and continuum-scale modeling with curated online databases to ameliorate the problem of unreliable inputs to Hg speciation modeling. Our AQUA-MER databases and web server (https://aquamer.ornl.gov) provides direct speciation results by combining web-based interfaces to a speciation calculator, databases of thermodynamic constants, and a computational chemistry toolkit to estimate missing constants. Although Hg is presented as a concrete use case, AQUA-MER can also be readily applied to other elements. © 2019 Wiley Periodicals, Inc.

Precious metal-free molecular machines for solar thermal energy storage

Four benzothiazolium crown ether-containing styryl dyes were prepared through an optimized synthetic procedure. Two of the dyes (4b and 4d) having substituents in the 5-position of the benzothiazole ring are newly synthesized compounds. They demonstrated a higher degree of trans–cis photoisomerization and a longer life time of the higher energy forms in comparison with the known analogs. The chemical structures of all dyes in the series were characterized by NMR, UV–vis, IR spectroscopy and elemental analysis. The steady-state photophysical properties of the dyes were elucidated. The stability constants of metal complexes were determined and are in good agreement with the literature data for reference dyes. The temporal evolution of trans-to-cis isomerization was observed in a real-time regime. The dyes demonstrated a low intrinsic fluorescence of their Ba²⁺ complexes and high yield of E/Z photoisomerization with lifetimes of the higher energy form longer than 500 seconds. Density functional theory (DFT) calculations at the B3LYP/6-31+G(d,p) level were performed in order to predict the enthalpies (H) of the cis and trans isomers and the storage energies (ΔH) for the systems studied.

Selective heavy metal removal and water purification by microfluidically-generated chitosan microspheres: Characteristics, modeling and application

Many industrial wastewater streams contain heavy metals, posing serious and irreversible damage to humans and living organisms, even at low concentrations due to their high toxicity and persistence in the environment. In this study, high-performance monodispersed chitosan (CS) microspheres were prepared using a simple microfluidic method and evaluated for metal removal from contaminated water. Batch experiments were carried out to evaluate the adsorption characteristics for the removal of copper ions, one representative heavy metal, from aqueous solutions. The inherent advantages of microfluidics enabled a precise control of particle size (CV = 2.3%), while exhibiting outstanding selectivity towards target ions (adsorption capacity 75.52 mg g^-1) and fair regeneration (re-adsorption efficiency 74% after 5 cycles). An integrated adsorption mechanism analytic system was developed based on different adsorption kinetics and isotherms models, providing an excellent adsorption prediction model with pseudo-second order kinetics (R² = 0.999), while the isotherm was fitted best to the Langmuir model (R² = 0.998). The multi-step adsorption process was revealed via quantitative measurements and schematically described. Selective adsorption performance of CS microspheres in the present of other competitive metal ions with different valence states has been demonstrated and studied by both experimental and density functional theory (DFT) analysis.

Computational Tools for Calculating log β Values of Geochemically Relevant Uranium Organometallic Complexes

Uranium (U^VI) interacts with organic ligands, subsequently controlling its aqueous chemistry. It is therefore imperative to assess the binding ability of natural organic molecules. We evidence that density functional theory (DFT) can be used as a practical protocol for predicting the stability of U^VI organic ligand complexes, allowing for the development of a relative stability series for organic complexes with limited experimental data. Solvation methods and DFT settings were benchmarked to suggest a suitable off-the-shelf solution. The results indicate that the IEFPCM solvation method should be employed. A mixed solvation approach improves the accuracy of the calculated stability constant (log β); however, the calculated log β are approximately five times more favorable than experimental data. Different basis sets, functionals, and effective core potentials were tested to check that there were no major changes in molecular geometries and Δ_r G. The recommended method employed is the B3LYP functional, aug-cc-pVDZ basis set for ligands, MDF60 ECP and basis set for U^VI, and the IEFPCM solvation model. Using the fitting approach employed in the literature with these updated DFT settings allows fitting of 1:1 U^VI complexes with root-mean-square deviation of 1.38 log β units. Fitting multiple bound carboxylate ligands indicates a second, separate fitting for 1:2 and 1:3 complexes.

Metal leaching from refinery waste hydroprocessing catalyst

The present study aims to develop an eco-friendly methodology for the recovery of nickel (Ni), molybdenum (Mo), and vanadium (V) from the refinery waste spent hydroprocessing catalyst. The proposed process has two stages: the first stage is to separate alumina, while the second stage involves the separation of metal compounds. The effectiveness of leaching agents, such as NH₄OH, (NH₄)₂CO₃, and (NH₄)₂S₂O₈, for the extraction of Mo, V, Ni, and Al from the refinery spent catalyst has been reported as a function of reagent concentration (0.5 to 2.0 molar), leaching time (1 to 6 h), and temperature (35 to 60°C). The optimal leaching conditions were achieved to obtain the maximum recovery of Mo, Ni, and V metals. The effect of the mixture of multi-ammonium salts on the metal extraction was also studied, which showed an adverse effect for Ni and V, while marginal improvement was observed for Mo leaching. The ammonium salts can form soluble metal complexes, in which stability or solubility depends on the nature of ammonium salt and the reaction conditions. The extracted metals and support can be reused to synthesize a fresh hydroprocessing catalyst. The process will reduce the refinery waste and recover the expensive metals. Therefore, the process is not only important from an environmental point of view but also vital from an economic perspective.

Post-synthetic modifications of cadmium-based knots and links

Three topologically non-trivial cadmium(ii)-based complexes-Cd-[2]C, Cd-TK and Cd-SL-were simultaneously self-assembled in a dynamic library, individually isolated and fully characterized using solid-state, gas-phase and solution-phase techniques. Post-synthetic modifications, including reduction and transmetalation, were subsequently achieved. Imine bond reduction followed by demetallation led to the isolation of the corresponding organic molecules [2]C, TK and SL. Transmetalation of Cd-TK and Cd-SL with the zinc(ii) cation resulted in isolation of the corresponding zinc(ii)-containing complexes Zn-TK and Zn-SL.

Tuning the affinity of catechols and salicylic acids towards Al(iii): characterization of Al–chelator interactions

Aluminum is a non-essential element in the human body with unclear harmful effects; therefore, the design and tuning of new and efficient Al(iii) chelating agents is a subject of paramount importance nowadays.

Pseudo electron-deficient organometallics: limited reactivity towards electron-donating ligands

This work presents the unusual reactivity of a family of electron-deficient half-sandwich metal complexes.

Toward Quantitatively Accurate Calculation of the Redox-Associated Acid-Base and Ligand Binding Equilibria of Aquacobalamin

Redox processes in complex transition metal-containing species are often intimately associated with changes in ligand protonation states and metal coordination number. A major challenge is therefore to develop consistent computational approaches for computing pH-dependent redox and ligand dissociation properties of organometallic species. Reduction of the Co center in the vitamin B12 derivative aquacobalamin can be accompanied by ligand dissociation, protonation, or both, making these properties difficult to compute accurately. We examine this challenge here by using density functional theory and continuum solvation to compute Co-ligand binding equilibrium constants (Kon/off), pKas, and reduction potentials for models of aquacobalamin in aqueous solution. We consider two models for cobalamin ligand coordination: the first follows the hexa, penta, tetra coordination scheme for Co(III), Co(II), and Co(I) species, respectively, and the second model features saturation of each vacant axial coordination site on Co(II) and Co(I) species with a single, explicit water molecule to maintain six directly interacting ligands or water molecules in each oxidation state. Comparing these two coordination schemes in combination with five dispersion-corrected density functionals, we find that the accuracy of the computed properties is largely independent of the scheme used, but including only a continuum representation of the solvent yields marginally better results than saturating the first solvation shell around Co throughout. PBE performs best, displaying balanced accuracy and superior performance overall, with RMS errors of 80 mV for seven reduction potentials, 2.0 log units for five pKas and 2.3 log units for two log Kon/off values for the aquacobalamin system. Furthermore, we find that the BP86 functional commonly used in corrinoid studies suffers from erratic behavior and inaccurate descriptions of Co-axial ligand binding, leading to substantial errors in predicted pKas and Kon/off values. These findings demonstrate the effectiveness of the present approach for computing electrochemical and thermodynamic properties of a complex transition metal-containing cofactor.

Theoretical prediction of coordination environments and stability constants of lanthanum lactate complexes in solution

Using Density Functional Theory calculations in combination with explicit solvent and a continuum solvent model, this work sets out to understand the coordination environment and relevant thermodynamics of La(iii)-lactate complexes.

The synthesis and unexpected solution chemistry of thermochromic carborane-containing osmium half-sandwich complexes

The functionalisation of the 16-electron complex [Os(η⁶-p-cymene)(1,2-dicarba-closo-dodecarborane-1,2-dithiolato)] (1) with a series of Lewis bases to give the corresponding 18-electron complexes is reported.

Theoretical prediction of Am(iii)/Eu(iii) selectivity to aid the design of actinide-lanthanide separation agents

Selective extraction of minor actinides from lanthanides is a critical step in the reduction of radiotoxicity of spent nuclear fuels. However, the design of suitable ligands for separating chemically similar 4f- and 5f-block trivalent metal ions poses a significant challenge. First-principles calculations should play an important role in the design of new separation agents, but their ability to predict metal ion selectivity has not been systematically evaluated. In this work, we examine the ability of several density functional theory methods to predict selectivity of Am(iii) and Eu(iii) with oxygen, mixed oxygen-nitrogen, and sulfur donor ligands. The results establish a computational method capable of predicting the correct order of selectivities obtained from liquid-liquid extraction and aqueous phase complexation studies. To allow reasonably accurate predictions, it was critical to employ sufficiently flexible basis sets and provide proper account of solvation effects. The approach is utilized to estimate the selectivity of novel amide-functionalized diazine and 1,2,3-triazole ligands.

How simple is too simple? Computational perspective on importance of second-shell environment for metal-ion selectivity

The metal-ion selectivity in biomolecules represents one of the most important phenomena in bioinorganic chemistry. The open question to what extent is the selectivity in the complex bioinorganic structures such as metallopeptides determined by the first-shell ligands of the metal ion is answered herein using six model peptides complexed with the set of divalent metal ions (Mn(2+), Fe(2+), Co(2+), Ni(2+), Cu(2+), Zn(2+), Cd(2+), and Hg(2+)) and their various first-shell representations. By calculating the differences among the free energies of complexation of metal ions in these peptides and their model (truncated) systems it is quantitatively shown that the definition of the first shell is paramount to this discussion and revolves around the chemical nature of the binding site. Despite the vast conceivable diversity of peptidic structures, that suggest certain fluidity of this definition, major contributing factors are identified and assessed based on their importance for capturing metal-ion selectivity. These factors include soft/hard character of ligands and various non-covalent interactions in the vicinity of the binding site. The relative importance of these factors is considered and specific suggestions for effective construction of the models are made. The relationship of first-shell models and their corresponding parent peptides is discussed thoroughly, both with respect to their chemical similarity and potential disparity introduced by generally "non-alignable" conformational flexibility of the two systems. It is concluded that, in special cases, this disparity can be negligible and that heeding the chemical factors contributing to selectivity during construction of the model can successfully result in models that retain the affinity profile for various metal ions with high fidelity.