Molecular Hosts for the Sensing and Separation of 99TcO4.

In recent years, European Union member states have hastened energy policy deliberations to address supply and sustainability concerns, placing a significant emphasis on nuclear energy as a means to achieve decarbonization goals. However, despite its significant role in power generation, nuclear energy faces significant challenges linked to fuel reprocessing and waste disposal, that hinder its broader expansion. In this context, the separation of technetium represents a concerning issue. Indeed, technetium's catalytic activity can impede the extraction of uranium, neptunium, and plutonium, affecting waste reprocessing efficiency. Additionally, the stable form of technetium in aerobic conditions, pertechnetate (TcO4 -), poses risks of groundwater contamination due to its mobility and solubility. Hence, sensing and separation of TcO4 - is imperative for both nuclear fuel processing and minimising radioactive contamination in the environment. However, the binding of TcO4 - and its separation from contaminated solutions present challenges due to the acidic (or basic) waste components and the high ionic strength in real matrices. Supramolecular chemists have addressed these issues by designing receptors inspired by molecular recognition principles. This article explores recent advancements and future directions in TcO4 - sensing and separation (using extraction and sorption) with a focus on molecular hosts. Metal-organic receptors will also be discussed.

Modeling Lanthanide Ions in Solution: A Versatile Force Field in Aqueous and Organic Solvents.

In this paper, we propose a new nonpolarizable force field for describing the Ln3+ (Ln = lanthanide) series based on a 12-6-4 Lennard-Jones potential. The development of the force field was performed in pure water by adjusting both the ion-oxygen distance and the hydration free energy. This force field accurately reproduces the Ln3+ hydration properties through the series, especially the coordination number that is hardly accessible using a nonpolarizable force field. Then, the validity and the transferability of the current force field were evaluated for two different systems containing Ln3+ in various solvents, namely, 0.1 mol L-1 La(NO3)3 salts in methanol and Eu(NO3)3 salts in solvent organic phases composed of DMDOHEMA molecules in n-heptane. The good agreement between our simulations and the data available in the literature confirms the accuracy of the force field for describing the lanthanide cations in both aqueous and nonaqueous media.

Investigation of the Plutonium(IV) Interactions with Two Variants of the EF-Hand Ca-Binding Site I of Calmodulin.

Due to its presence in the nuclear industry and its strong radiotoxicity, plutonium is an actinide of major interest in the event of internal contamination. To improve the understanding of its mechanisms of transport and accumulation in the body, the complexation of Pu(IV) to the most common protein calcium-binding motif found in cells, the EF-hand motif of calmodulin, was investigated. Visible and X-ray absorption spectroscopies (XAS) in solution made it possible to investigate the speciation of plutonium at physiological pH (pH 7.4) and pH 6 in two variants of the calmodulin Ca-binding site I and using Pu(IV) in different media: carbonate, chloride, or nitrate solutions. Three different species of Pu were identified in the samples, with formation of 1:1 Pu(IV):calmodulin peptide complexes, Pu(IV) reduction, and formation of peptide-mediated Pu(IV) hexanuclear cluster.

Optimized Coordination of Uranyl in Engineered Calmodulin Site 1 Provides a Subnanomolar Affinity for Uranyl and a Strong Uranyl versus Calcium Selectivity.

As an alpha emitter and chemical toxicant, uranium toxicity in living organisms is driven by its molecular interactions. It is therefore essential to identify main determinants of uranium affinity for proteins. Others and we showed that introducing a phosphoryl group in the coordination sphere of uranyl confers a strong affinity of proteins for uranyl. In this work, using calmodulin site 1 as a template, we modulate the structural organization of a metal-binding loop comprising carboxylate and/or carbonyl ligands and reach affinities for uranyl comparable to that provided by introducing a strong phosphoryl ligand. Shortening the metal binding loop of calmodulin site 1 from 12 to 10 amino acids in CaMΔ increases the uranyl-binding affinity by about 2 orders of magnitude to log KpH7 = 9.55 ± 0.11 (KdpH7 = 280 ± 60 pM). Structural analysis by FTIR, XAS, and molecular dynamics simulations suggests an optimized coordination of the CaMΔ-uranyl complex involving bidentate and monodentate carboxylate groups in the uranyl equatorial plane. The main role of this coordination sphere in reaching subnanomolar dissociation constants for uranyl is supported by similar uranyl affinities obtained in a cyclic peptide reproducing CaMΔ binding loop. In addition, CaMΔ presents a uranyl/calcium selectivity of 107 that is even higher in the cyclic peptide.

Inter-Site Cooperativity of Calmodulin N-Terminal Domain and Phosphorylation Synergistically Improve the Affinity and Selectivity for Uranyl.

Uranyl-protein interactions participate in uranyl trafficking or toxicity to cells. In addition to their qualitative identification, thermodynamic data are needed to predict predominant mechanisms that they mediate in vivo. We previously showed that uranyl can substitute calcium at the canonical EF-hand binding motif of calmodulin (CaM) site I. Here, we investigate thermodynamic properties of uranyl interaction with site II and with the whole CaM N-terminal domain by spectrofluorimetry and ITC. Site II has an affinity for uranyl about 10 times lower than site I. Uranyl binding at site I is exothermic with a large enthalpic contribution, while for site II, the enthalpic contribution to the Gibbs free energy of binding is about 10 times lower than the entropic term. For the N-terminal domain, macroscopic binding constants for uranyl are two to three orders of magnitude higher than for calcium. A positive cooperative process driven by entropy increases the second uranyl-binding event as compared with the first one, with ΔΔG = -2.0 ± 0.4 kJ mol-1, vs. ΔΔG = -6.1 ± 0.1 kJ mol-1 for calcium. Site I phosphorylation largely increases both site I and site II affinity for uranyl and uranyl-binding cooperativity. Combining site I phosphorylation and site II Thr7Trp mutation leads to picomolar dissociation constants Kd1 = 1.7 ± 0.3 pM and Kd2 = 196 ± 21 pM at pH 7. A structural model obtained by MD simulations suggests a structural role of site I phosphorylation in the affinity modulation.

Force Field Parameterization of Actinyl Molecular Cations Using the 12-6-4 Model.

In this work, a set of 12-6-4 force fields (FFs) parameters were developed for the actinyl molecular cations, AnO2n+ (n = 1, 2), from uranium to plutonium for classical molecular dynamics (MD) for four water models: TIP3P, SPC/E, OPC3, and TIP4Pew. Such a non-bonded potential model taking into account the induced dipole between the metallic center and the surrounding molecules has shown better performances for various cations than the classic 12-6 non-bonded potentials. The parametrization method proposed elsewhere for metallic cations has been extended to these molecular cations. In contrast to the actinyl 12-6 FFs from the literature, the new models reproduce correctly both solvation and thermodynamic properties, thanks to the inclusion of the induced dipole term (C4). The transferability of such force fields was assessed by performing MD simulations of carbonato actinyl species, which are highly implicated in actinide migration or actinide extraction from seawater. A highly satisfying agreement was found when comparing the EXAFS signals computed from our MD simulation to the experimental ones. The set of FFs developed here opens new possibilities for the study of actinide chemistry.

Effect of metal complexation on diglycolamide radiolysis: a comparison between ex situ gamma and in situ alpha irradiation.

Radiolytic degradation is an important aspect to consider when developing a ligand or a complexant for radionuclides. Diglycolamide extractants (DGAs) have been playing an important role in many partition processes for spent nuclear fuel. In particular, the extractant N,N,N'N'-tetraoctyl diglycolamide (TODGA) has been studied intensively for the purpose of solvent extraction processes such as ARTIST, i-Sanex, EURO-GANEX and EURO-EXAM, which have been developed around the TODGA extractant. For the first time, the radiolytic stability of TODGA was investigated both by in situ alpha irradiation using a macroconcentration of americum(III) and by ex situ gamma irradiation in the presence of a macroconcentration of neodymium(III). It was shown that metal ions complexed in the organic phase protect TODGA from degradation by irradiation and that the degradation was slower using in situ alpha irradiation compared to ex situ gamma irradiation. By comparison to gamma irradiation of Nd-TODGA solution, alpha irradiation of Am-TODGA solutions showed the presence of 2 additional compounds identified as a TODGA molecule with a CC bond and a TODGA molecule with the addition of a NO3 group. The major degradation products were identified and a degradation schema was proposed. The direct analysis of the solution containing Am(III) or Nd(III) showed that the degradation compounds retaining a diglycolamide skeleton are involved in heteroleptic complexes with TODGA, without a negative impact on An(III) or Ln(III) complexation.

Aggregation of Bifunctional Extractants Used for Uranium(VI) Separation.

DEHCNPB (butyl-N,N-di(2-ethylhexyl)carbamoyl-nonylphosphonate) is an amido-phosphonic acid that has remarkable properties for the separation of uranium from wet phosphoric acid. Despite previous studies, a detailed description of the DEHCNPB organic solutions at the supramolecular and molecular scales is missing. In the present work, we use classical Molecular Dynamics (MD) combined with SANS and SAXS experimental data in order to describe the aggregation of the bifunctional extractant DEHCNPB as well as the speciation of uranium(VI) in such systems. We provide a fine description of the molecular species in the organic solution and of the interactions within the aggregates formed, shedding light on solvent extraction mechanisms. Without uranium, the organic phase is highly composed of dimers and trimers H-bonded through phosphonate functions and without water molecules. With uranium, two to three extractant molecules coordinate directly the uranyl cation by their phosphonate groups. Uranyl is not fully dehydrated in this organic solution, and the amide groups of the extractants are found to form H-bonds with the water molecules bound to uranyl. These H-bond networks around the metallic cation stabilize the complexes and facilitate the extraction. These results underline the importance of considering weak interactions in the understanding of extraction processes and demonstrate how molecular simulations provide essential insights into such complex organic phase chemistry with a high number of species.

Thermodynamics of Malonamide Aggregation Deduced from Molecular Dynamics Simulations.

The aggregation of malonamide extractants diluted in an aliphatic solvent phase has been studied in the presence of water by molecular dynamics simulation. Using association criteria based on distances between molecules and graphs theory, the aggregate distribution has been computed and the corresponding Gibbs energy of aggregates and mass action law constants have been determined. Finally, a model allowing us to the compute critical micelle concentration and osmotic data for a variable concentration of extractants, with or without a correction of the organic phase activity, was developed. It appears however that the accurate depiction of the aggregation allows modeling the thermodynamics of the solution even without an explicit calculation of the activity: both models give results in good agreement with the experiments.

Syntheses and evaluation of new hydrophilic azacryptands used as masking agents of technetium in solvent extraction processes.

The extraction of technetium, present in nitric acid medium as pertechnetate anion, is an issue in solvent extraction processes used to recover uranium and plutonium. In the present study, a complexing agent is added in the aqueous nitric acid solution to bind selectively the pertechnetate anion and prevent its extraction into the organic phase or to back extract it in the aqueous phase. Several azacryptands with the addition of hydrophilic groups were synthesized to improve the solubility of the previously studied cage molecule in nitric acid medium. Solvent extraction tests reveal that all the synthesized ligands have a similar complexation strength towards pertechnetate and are able to maintain this anion in the aqueous phase (0.5 M HNO3). These ligands are able to overcome the Hofmeister bias and selectively bind technetium in nitric acid solution. The azacryptand concentration can be increased by a factor of three in the liquid-liquid extraction conditions compared to our previous work. Coordination studies using microcalorimetry, Single Crystal X-Ray Diffraction (SC-XRD), infrared and Raman spectroscopies show the formation of an inclusion complex with hydrogen bonds stabilizing the oxo-anion within the cavity. This solubility improvement is promising for the introduction of this kind of macrocyclic azacryptands in a solvent extraction process.

Liquid/liquid interface in periodic boundary condition.

We study how surface phenomena can change the interface geometry in liquid-liquid two-phase systems with periodic boundary conditions. Without any curvature effect on surface tension, planar (slab), cylindrical, and spherical structures are successively obtained as a function of the total composition and elongation of the box, in accordance with molecular dynamics simulations for a water/heptane system. The curvature effects described by Tolman relationship desymmetrize the phase diagram by stabilizing a concavity but it leads to inconsistencies with high curvature. Helfrich model partially resolves this and predicts the possible presence of shells reflecting a frustrated system.

Perrhenate and pertechnetate complexation by an azacryptand in nitric acid medium.

Technetium is present as the pertechnetate anion in spent nuclear fuel solutions, and its extraction by several extractant systems is a major problem for the liquid-liquid extraction processes used to separate uranium and plutonium. To prevent technetium extraction into the organic phase, a complexing agent may be added to the aqueous nitric acid phase to selectively bind the pertechnetate anion. In the present study, liquid-liquid extraction experiments reveal that technetium distribution ratios are considerably lowered with addition of an azacryptand, which is a good receptor for pertechnetate anion recognition. This ligand is able to overcome the Hofmeister bias and selectively bind techetium in nitric acid solution. Coordination studies using infrared and Raman spectoscopies and DFT calculations show the formation of an inclusion complex with hydrogen bonds stabilizing the oxo-anion within the cavity. For the first time, the cage molecules are studied for an extraction process.

Effect of chemical environment on the radiation chemistry of N,N-di-(2-ethylhexyl)butyramide (DEHBA) and plutonium retention.

N,N-di-(2-ethylhexyl)butyramide (DEHBA) has been proposed as part of a hydro-reprocessing solvent extraction system for the co-extraction of uranium and plutonium from spent nuclear fuel, owing to its selectivity for hexavalent uranium and tetravalent plutonium. However, there is a critical lack of quantitative understanding regarding the impact of chemical environment on the radiation chemistry of DEHBA, and how this would affect process performance. Here we present a systematic investigation into the radiolytic degradation of DEHBA in a range of n-dodecane solvent system formulations, where we subject DEHBA to gamma irradiation, measure reaction kinetics, ligand integrity, degradation product formation, and investigate solvent system performance through uranium and plutonium extraction and strip distribution ratios. The rate of DEHBA degradation in n-dodecane was found to be slow (G = -0.31 ± 0.02 µmol J-1) but enhanced upon contact with the oxidizing conditions of the investigated solvent systems (organic-only, or in contact with either 0.1 or 3.0 M aqueous nitric acid). Two major degradation products were identified in the organic phase, bis-2-ethylhexylamine (b2EHA) and N-(2-ethylhexyl)butyramide (MEHBA), resulting from the cleavage of C-N bonds, and could account for the total loss of DEHBA up to ∼300 kGy for organic-only conditions. Both b2EHA and MEHBA were also found to be susceptible to radiolytic degradation, having G-values of -0.12 ± 0.01 and -0.08 ± 0.01 µmol J-1, respectively. Solvent extraction studies showed: (i) negligible change in uranium extraction and stripping with increasing absorbed dose; and (ii) plutonium extraction and retention exhibits complex dependencies on absorbed dose and chemical environment. Organic-only conditions afforded enhanced plutonium extraction and retention attributed to b2EHA, while acid contacts inhibited this effect and promoted significant plutonium retention for the highest acidity. Overall it has been demonstrated that chemical environment during irradiation has a significant influence on the extent of DEHBA degradation and plutonium retention.

Probing the existence of uranyl trisulfate structures in the AMEX solvent extraction process.

Knowledge of the complex microstructure in solvent extraction phases is mandatory for a full comprehension of ionic separation. Coupling EXAFS with MD simulations for uranyl extraction in sulfuric media with tertiary amine extractants enabled unravelling of the unprecedented uranyl tri-sulfate structure.

Insight of the Metal-Ligand Interaction in f-Element Complexes by Paramagnetic NMR Spectroscopy.

The magnetic properties of LnIII and AnIII complexes formed with dipicolinate ligands have been studied by NMR spectroscopy. To know precisely the geometries of these complexes, a crystallographic study by single-crystal X-ray diffraction (XRD) and extended X-ray absorption fine structure (EXAFS) in solution was performed. Several methods to separate the paramagnetic shifts observed in the NMR spectra were applied to these complexes. Methods using a number of nuclei of the dipicolinate ligands revealed an abrupt change in the geometries of the complexes and a metal-ligand interaction in the middle of the lanthanide series. A study of the variation of the paramagnetic shifts with temperature demonstrated that higher-order terms of the dipolar and contact contributions are required, especially for the lightest LnIII and almost all the studied AnIII . Bleaney's parameters a and C a D relating to the contact and dipolar terms, respectively, were deduced from experimental data and compared with the results of ab initio calculations. Quite a good agreement was found for the temperature dependencies of a and C a D . However, the C a D values obtained from cation magnetic anisotropy calculations showed some discrepancies with the values derived from Bleaney's equation defined for LnIII . Other parameters, such as the crystal field parameter and the hyperfine constants Fi obtained from the experimental data of the [An(ethyl-dpa)3 ]3- complexes (ethyl-dpa=4-ethyl-2,6-dipicolinic acid), are at odds with the assumptions underlying Bleaney's theory.

UO22+ structure in solvent extraction phases resolved at molecular and supramolecular scales: a combined molecular dynamics, EXAFS and SWAXS approach.

A new polarizable force field for describing the solvation of the uranyl (UO22+) cation in solvent extraction phases has been developed for molecular dynamics simulations. The validity of the polarizable force field has been established by comparison with EXAFS and SWAXS experiments. This new force field allows for describing both the UO22+ hydration and solvation properties in good agreement with the experiments. In aqueous phases we demonstrated that the UO22+ force field has been improved from the previous one we developed. Indeed, the UO22+ structural and dynamics properties, i.e., the dynamics of the water molecules in the vicinity of the uranyl cation, calculated from molecular dynamics are in better agreement with the EXAFS experiments. Furthermore, the transferability of the UO22+ force field proposed here has been validated on typical solvent extraction phases containing uranyl nitrate salts with extractant molecules, namely DMDOHEMA molecules, in n-heptane. The good agreements observed between the theoretical (MD simulations) and experimental UO22+ structures at the molecular (EXAFS) and supramolecular (SWAXS) scales prove the accuracy of the UO22+ force field developed and proposed in the present paper.

Activity Coefficients of Aqueous Sodium, Calcium, and Europium Nitrate Solutions from Osmotic Equilibrium MD Simulations.

Osmotic and activity coefficients of three aqueous electrolyte solutions with cations of similar ionic radius, but different charges, are described by molecular dynamics with the help of the osmotic equilibrium method using polarizable force fields up to high concentration. Simulations of vapor-liquid interfaces of aqueous solutions of NaNO3, Ca(NO3)2, and Eu(NO3)3 at different concentrations and at 298.15 K provide time-averaged number density profiles and consequently the quantity of solvent molecules in the vapor phase. These three cations of similar ionic radii exhibit an increasing amount of water in their first coordination sphere due to their increasing charge. The solvent activity is directly determined by the vapor phase density at different salt concentrations with respect to the vapor phase density of the pure solvent. The obtained densities of the liquid bulk and the osmotic and activity coefficients for the three different nitrate salts are in good agreement with the experimental results. Time-averaged concentration profiles and the interpretation of radial distribution functions are used to explain the role of coordination on the thermodynamic properties of aqueous electrolyte solutions.

Structural Analysis of Uranyl Complexation by the EF-Hand Motif of Calmodulin: Effect of Phosphorylation.

Better understanding of uranyl-protein interactions is a prerequisite to predict uranium chemical toxicity in cells. The EF-hand motif of the calmodulin site I is about thousand times more affine for uranyl than for calcium, and threonine phosphorylation increases the uranyl affinity by two orders of magnitude at pH 7. In this study, we confront X-ray absorption spectroscopy with Fourier transform infrared (FTIR) spectroscopy, time-resolved laser-induced fluorescence spectroscopy (TRLFS), and structural models obtained by molecular dynamics simulations to analyze the uranyl coordination in the native and phosphorylated calmodulin site I. For the native site I, extended X-ray absorption fine structure (EXAFS) data evidence a short U-Oeq distance, in addition to distances compatible with mono- and bidentate coordination by carboxylate groups. Further analysis of uranyl speciation by TRLFS and thorough investigation of the fluorescence decay kinetics strongly support the presence of a hydroxide uranyl ligand. For a phosphorylated site I, the EXAFS and FTIR data support a monodentate uranyl coordination by the phosphoryl group and strong interaction with mono- and bidentate carboxylate ligands. This study confirms the important role of a phosphoryl ligand in the stability of uranyl-protein interactions. By evidencing a hydroxide uranyl ligand in calmodulin site I, this study also highlights the possible role of less studied ligands as water or hydroxide ions in the stability of protein-uranyl complexes.

Simulating Osmotic Equilibria: A New Tool for Calculating Activity Coefficients in Concentrated Aqueous Salt Solutions.

Herein, a new theoretical method is presented for predicting osmotic equilibria and activities, where a bulk liquid and its corresponding vapor phase are simulated by means of molecular dynamics using explicit polarization. Calculated time-averaged number density profiles provide the amount of evaporated molecules present in the vapor phase and consequently the vapor-phase density. The activity of the solvent and the corresponding osmotic coefficient are determined by the vapor density at different solute concentrations with respect to the reference vapor density of the pure solvent. With the extended Debye-Hückel equation for the activity coefficient along with the corresponding Gibbs-Duhem relation, the activity coefficients of the solutes are calculated by fitting the osmotic coefficients. A simple model based on the combination of Poisson processes and Maxwell-Boltzmann velocity distributions is introduced to interpret statistical phenomena observed during the simulations, which are related to evaporation and recondensation. This method is applied to aqueous dysprosium nitrate [Dy(NO3)3] solutions at different concentrations. The obtained densities of the liquid bulk and the osmotic and activity coefficients are in good agreement with the experimental results for concentrated and saturated solutions. Density profiles of the liquid-vapor interface at different concentrations provide detailed insight into the spatial distributions of all compounds.