Lu
Liu
a,
Seyfeddine
Rahali
bc,
Rémi
Maurice
a,
Cecilia
Gomez Pech
ab,
Gilles
Montavon
a,
Jean-Yves
Le Questel
b,
Jérôme
Graton
b,
Julie
Champion
*a and
Nicolas
Galland
*b
aSUBATECH UMR 6457, CNRS, IMT Atlantique, Université de Nantes, 4 Rue Alfred Kastler, 44307 Nantes, France. E-mail: julie.champion@subatech.in2p3.fr
bUniversité de Nantes, CNRS, CEISAM UMR 6230, 44000 Nantes, France. E-mail: nicolas.galland@univ-nantes.fr
cDepartment of Chemistry, College of Science and Arts, Qassim University, Ar Rass, Saudi Arabia
First published on 12th July 2021
As a non-covalent interaction, halogen bonding is now acknowledged to be useful in all fields where the control of intermolecular recognition plays a pivotal role. Halogen-bond basicity scales allow quantification of the halogen bonding of referential donors with organic functional groups from a thermodynamic point of view. Herein we present the pKBAtI basicity scale to provide the community an overview of halogen-bond acceptor strength towards astatine, the most potent halogen-bond donor element. This experimental scale is erected on the basis of complexation constants measured between astatine monoiodide (AtI) and sixteen selected Lewis bases. It spans over 6 log units and culminates with a value of 5.69 ± 0.32 for N,N,N′,N′-tetramethylthiourea. On this scale, the carbon π-bases are the weakest acceptors, the oxygen derivatives cover almost two-thirds of the scale, and sulphur bases exhibit the highest AtI basicity. Regarding the applications of 211At in targeted radionuclide therapy, stronger labelling of carrier agents could be envisaged on the basis of the pKBAtI scale.
Astatine (At, Z = 85) is the heaviest naturally occurring halogen element. By combining liquid/liquid competition experiments and quantum mechanical calculations, we have recently evidenced the first XBs involving At in complexes between astatine monoiodide (AtI) and a few Lewis bases. Among the halogens, we confirmed the highest donating ability of astatine, especially compared to iodine.17,18 There are fundamental interests to get more experimental data, not only in order to confirm previous findings, but also to increase our understanding of halogen bonding from the perspective of the most potent halogen element, in particular for the refinement of theoretical models.19,20 On the other hand, the 211 radioisotope of astatine presents potential use in targeted therapy of cancers.21–23 Practical interests notably lie in the improvement of the radiolabelling strategy for these applications.24–28 For example, an At-mediated XB can be assumed to explain, in particular, the efficient radiolabelling of bis-(nido-carboranylmethyl)benzene derivatives,29 the latter having been considered as pendant groups for attaching 211At to carrier biomolecules. Therefore, we proposed to investigate on a more extended and diversified chemical sample the various features (structural, thermodynamic) of XBs involving astatine. One may expect to extend the pKBI2 scale to astatine species. However, some AtI-mediated complexes show some different behaviors compared to the trends depicted by pKBI2 (diethyl ether being for instance a significant stronger base compared to hexamethyl benzene according to ref. 17). Moreover, the metallic character of At30,31 makes its similarity with other halogen elements questionable. In this work, we aim therefore to establish an appropriate basicity scale of At-mediated XBs gathering representative chemical functionalities of Lewis bases, and to look for the strongest XB interactions involving At.
It is worth noting that astatine is considered as an “invisible” element.21,32 Indeed, it is a short-lived radioelement (t1/2 ≤ 8.1 h). Only minute quantities can be produced artificially, leading to astatine solutions of concentration typically below 10−10 mol L−1. The conventional spectroscopic tools are therefore inapplicable and most information is obtained by indirect methods. In particular, we successfully evaluated the reactivity of At-compounds by studying the distribution of 211At radionuclide in biphasic liquid systems.17,18,30,31,33 However, a careful selection of the XB donor and acceptors was mandatory to guarantee the occurrence of XB interactions. Among the dihalogens with strong potent donor ability,34 the At2 and AtF species are still hypothetical. The narrow predominance domains of AtBr and AtCl in aqueous phases and their apparent immiscibility in organic solvents35 prohibit so far any measurements involving these XB donors. As a result, AtI imposes itself as the referential XB donor. Furthermore, the choice of the XB acceptors is severely limited. They must be stable under acidic and oxidizing conditions defining the predominance domain of the AtI species. Therefore, most of the nitrogen Lewis bases, such as amines, pyridines, anilines, guanidines or imines, are experimentally inaccessible owing to their protonation state under the applied conditions. The experimental constraints also imply that the ligands are soluble in organic solvents to promote, upon the formation of XBs, the transfer of astatine species from the aqueous phase to the organic one. Accordingly, we have considered eight new ligands with regard to our previous studies,17,18 which introduce the following chemical functions: ketone, thioketone, amide, thioamide, urea, thiourea, selenoether and nitrile (Fig. 1).
![]() | (1) |
The distribution ratio of At (D) is calculated as the ratio of the volumetric activity of 211At between the organic phase and the aqueous one. In a series of experiments, the ligand concentration in the organic phase is increased, which is supposed to influence the At speciation and thus induce a variation of D.
Upon increase of the ligand concentration, the distribution of At can display two typical behaviors. Fig. 2a shows the results for the Lewis base 3, which is also representative of the evolution observed for 1 and 4 (see Fig. S1 in the ESI†). In the low ligand concentration range, the distribution of At is actually ruled by the ratio between AtI and AtI2− in the aqueous phase and the extraction of AtI in the organic phase; D is thus constant. Then, a sharp increase of D can be observed, which implies the complexation of AtI with ligands in the organic phase, leading to a transfer of At species from the aqueous phase to the organic one. However, as shown in Fig. 2b and S1,† an additional feature appears in the cases of 5, 6 and 7: an inflection point of the distribution vs. concentration curve at the highest ligand concentrations. This behaviour could be rationalized by the hypothesis that B⋯AtI or an additional astatine species is present in the aqueous phase. Because the former possibility is the most probable one under the experimental conditions, as previously shown for other Lewis bases,17 we here consider that a 1:
1 complex is also formed in the aqueous phase.
The changes of D were reproduced quantitatively from two thermodynamic models of the biphasic system to derive the value of the KBAtI complexation constant. In this indirect approach, it is noteworthy that we limit as much as possible the number of thermodynamic equilibria involved in the models, i.e., the number of parameters to be adjusted during the fitting. The considered reactions, together with the corresponding analytical expressions of D used to fit the experimental data, are presented in Section 3.3 in the ESI.† The relevance of the chosen models (number of species, equilibrium constants) is notably supported by the results of quantum mechanical calculations (vide infra). For the ligands 2 and 8, the D values remain constant over the whole investigated concentration range, preventing any determination of the complexation constants (see Fig. S2 in the ESI†). Finally, the whole sample of KBAtI equilibrium constants is reported in Table 1, together with previously determined XB complexation constants.17,18
Lewis base | No. | Solvent | log![]() |
||
---|---|---|---|---|---|
Experimenta | PW6B95/TZVPD | B3LYP/TZVPD | |||
a The values in parenthesis correspond to two standard deviations. | |||||
Cyclohexanone | 1 | Cyclohexane | 1.61(37) | 1.47 | 1.66 |
N,N,N′,N′-Tetramethylurea | 3 | Heptane | 1.76(50) | 2.59 | 2.52 |
(1R)-(−)-Thiocamphor | 4 | Heptane | 3.25(17) | 3.58 | 3.46 |
N,N-Dimethyl thioacetamide | 5 | Heptane | 3.91(86) | 4.15 | 4.09 |
N,N,N′,N′-Tetramethylthiourea | 6 | Cyclohexane | 5.69(32) | 5.28 | 5.02 |
Dimethyl selenide | 7 | Cyclohexane | 3.60(44) | 3.69 | 3.79 |
Toluene | 9 | Cyclohexane | −0.67(24)17 | −0.31 | −0.62 |
Ethyl butanoate | 10 | Cyclohexane | 0.46(56)17 | 0.56 | 0.87 |
Hexamethylbenzene | 11 | Cyclohexane | 0.67(64)17 | 0.92 | 0.02 |
Diethyl ether | 12 | Cyclohexane | 1.53(46)17 | 0.34 | 0.15 |
Diethyl methylphosphonate | 13 | Cyclohexane | 1.75(44)17 | 2.41 | 3.18 |
Tributyl phosphate | 14 | Cyclohexane | 2.84(26)17 | 2.30 | 3.26 |
Triphenylphosphine sulfide | 15 | Cyclohexane | 3.41(76)17 | 3.81 | 3.90 |
Dibutyl sulfoxide | 16 | Cyclohexane | 3.78(40)17 | 3.41 | 3.17 |
Diethyl sulfide | 17 | Cyclohexane | 4.01(62)17 | 3.09 | 2.87 |
Tributylphosphine oxide | 18 | Cyclohexane | 4.24(35)18 | 4.53 | 4.51 |
This work also reveals an important advantage of our approach based on biphasic systems. The interactions between the studied compounds can possibly be quantified in the two solvents, and simultaneously. Indeed, for the three strongest Lewis bases, 5, 6 and 7, their interactions with AtI and with water molecules in the aqueous phase become competitive. Therefore, the complexation constants between AtI and these ligands can also be determined in the aqueous phase (KwBAtI), as detailed in Section 3.2 in the ESI.† The obtained values can be larger (KwBAtI = 103.82 for 7) or smaller (KwBAtI = 103.80 for 5 and 104.40 for 6) than the corresponding KBAtI values. A comprehensive analysis of the influence of the solvent (organic, water) on the strength of the XB interactions becomes feasible, but it goes beyond the scope of this work.36 However, some information can be learned from crude calculations based on an implicit model of water (see Section 2 in the ESI†). The sum of the solvation free energies of AtI and of ligand 5 or 6 exceeds that of the corresponding adduct. Compared to the organic conditions, the adduct in the aqueous phase is therefore destabilized with respect to the reactants and the complexation constant is reduced. The solvation free energy in water of 7 is much weaker than the ones of 5 and 6, which (i) can be justified by the presence of nitrogen functions in 5 and 6, and (ii) is corroborated by a partition coefficient between organic and aqueous phases much bigger for 7 (101.40±0.16, see Section 1 in the ESI†) than for 5 and 6 (10−0.13±0.01 and 10−0.37±0.04, respectively). The weakest solvation free energy of 7 may explain the relative stabilization of its adduct with AtI when solvated in water, and the increased complexation constant. The KwBAtI constants can also be compared to constants measured in aqueous phases for the formation of other adducts of AtI, e.g. IAtBr− as the heaviest-known trihalogen species. This species can be viewed as an XB adduct between AtI and Br−,37 and the corresponding equilibrium constant in water (101.4±0.3)31 is weaker by more than two orders of magnitude. This outcome is probably the consequence of the much stronger solvation free energy for the small spherical Br− anion compared to neutral ligands 5–7, making IAtBr− relatively destabilized with respect to the other adducts.
As an interhalogen, AtI may form two types of XB interaction, one mediated by the astatine atom and the other one through the iodine atom. A descriptor commonly used to characterize the XB donating ability of a given donor is VS,max, that is a local maximum value of the electrostatic potential at the molecular surface (molecular electrostatic potential, MEP).2 Fig. S3a in the ESI† displays the calculated MEP for AtI at the PW6B95/TZVPD level of theory. Two positive regions are observed, one on the astatine side and the other one on the iodine side. Although astatine presents a significantly higher VS,max value (180.4 kJ mol−1) than iodine (71.3 kJ mol−1), we have nevertheless investigated the two types of interaction with the selected ligands. Fig. 3 displays the most-stable structures computed for the interactions between AtI and the ligands 1, 3–7, while Fig. S4 in the ESI† presents the structures obtained with the previously investigated ligands (9–18). Whatever the B⋯AtI system, the interaction between the two molecular fragments (AtI and B) is mediated by the astatine atom in the most stable structure. The structures stabilized by an interaction with the iodine atom are significantly less stable, for instance representing less than 7% of the whole population of complexes formed with toluene according to the Boltzmann distribution calculated at the PW6B95/TZVPD level of theory. Halogen bonding with the iodine atom of AtI is a minor phenomenon. Focusing on the At-mediated interactions with the ligands 1 and 3–7, they are all predicted with distances shorter than the sum of the van der Waals radii of the two involved atoms,40 by 22% for 7 up to 25% for 3 (see normalized interaction distances rXB from Table S3 in the ESI†). Furthermore, the angle formed between each ligand and AtI is close to 180°, the largest deviation being 3.7° for 7. These structural features are typical of halogen bonding. Hence, the theoretical calculations suggest that the interactions between AtI and the newly studied Lewis bases are stabilized by At-mediated XBs, as it was previously established for 9–18.17,18
It has previously been shown that the equilibrium constants KBAtI can be accurately determined,17,18 provided that the relativistic DFT calculations are carried out on isodesmic-like reactions as detailed in Section 2.2 in the ESI.† If 12 was chosen as the reference ligand, the set of the obtained complexation constants for the 16 B⋯AtI complexes is however translated according to the experimental values; the resulting values are given in Table 1. The PW6B95/TZVPD results are compared with the experimental ones in Fig. 4a. A strong relationship is established. The linear regression slope is very close to one (0.988), which constitutes a first clue that the property calculated corresponds to that measured (i.e., KBAtI equilibrium constants). The associated coefficient of determination R2 (0.976) is also quite good. Furthermore, the mean absolute deviation (MAD) between the experimental and computed logKBAtI values (0.45) is smaller than the average experimental uncertainty (0.47). The calculated values at the B3LYP/TZVPD level of theory also fairly match the experimental data (R2 = 0.968 and MAD = 0.55, see Fig. 4b). Thus, concerning the interactions between AtI and the ligands 1, 3–7, 9–18, the nature of the predicted most stable complexes and the agreement between the calculated KBAtI constants and the measured property allow us to conclude that the experimentally studied species are At-mediated XB complexes.
pKBAtI = log10(KBAtI) | (2) |
Gathering all the chemical functionalities available within the experimental limitations, the pKBAtI scale spans over 6 log units, as illustrated in Fig. 5. 6 is at the top of this scale with a pKBAtI value of 5.69. It must be stressed that nitrogen-based Lewis bases are unfortunately absent from this experimental dataset; even 8 which shows a low pKa value41 did not lead to any variations of the astatine distribution under our experimental conditions. Scanning the established scale, the carbon π-bases exhibit the weakest XB accepting ability, with pKBAtI from −0.67 to 0.67. The oxygen acceptor atomic sites are generally stronger XB bases, except the ester derivative 10. Reaching a pKBAtI value of 4.24 with the phosphine oxide compound 18, the oxygen family, therefore, covers almost 60% of this basicity scale. The sulphur family gathers even stronger XB acceptors, the scale culminating with the thiourea 6. Within the ketone (1vs.4), urea (3vs.6) and ether (12vs.17) families, the sulphur derivatives are systematically found to be stronger acceptors than the oxygen compounds with several orders of magnitude. These trends are fully supported by the computed logKBAtI values, clearly indicating the much better affinity of sulphur than oxygen for the At-containing XB donor. Conversely, selenoether 7 shows a weaker XB basicity (3.60) than the corresponding thioether 17 (4.01), which is not in line with the DFT calculations.
Note that the trends which emerge from the pKBAtI basicity scale can hardly be anticipated from the calculated structural parameters. No correlation exists in particular with the lengths of the interaction distances (R2 = 0.019), or even with the normalized interaction distances (rXB in Table S3†). However, the two carbon π-bases exhibit XBs with the largest deviations from linearity (6.3° and 7.8°) in line with their weak XB accepting ability. In addition, the lengthening of the At–I bond upon complexation is the most important for sulphur bases bearing an sp2 S atom, with a maximum of 0.092 Å for 6, in agreement with the thermodynamic trends.
Finally, the pKBAtI scale of astatine basicity can be compared to the pKBI2 scale of iodine basicity. As shown in Fig. 6, there is an overall linearity between both sets, with R2 = 0.916, supporting again the interactions with AtI through halogen bonding. For 15 Lewis bases among 16, the pKBAtI value is larger than the pKBI2 value, which can only be explained by considering a stronger donor ability of At among the halogen elements. Indeed, the propensity of a halogen atom to form XB interactions is commonly assumed to increase with increasing atom polarizability and decreasing atom electronegativity.2,42 This is supported in addition by the calculated MEP for AtI and I2. Fig. S3† shows that the value of VS,max on the At side of AtI, 180.4 kJ mol−1, is significantly higher than that calculated for diiodine, 129.1 kJ mol−1. The comparison to the I2 donor may also be extended to structural aspects. For instance, the interaction distances here calculated with the AtI donor (Table S3†) can be compared with those known between diiodine and the same or representative XB acceptors. Some interaction distances measured in crystallographic XB complexes with I2 and observed in the Cambridge structural database43 are gathered in Table S4 in the ESI.† Since the iodine atom is smaller than astatine, it is not unexpected that the lengths of the interactions involving iodine are shorter. This is more surprising if we consider, for the same type of acceptor site, the normalized interaction distances. While distances are weighted by the sum of the van der Waals radii of the two involved atoms, rXB is 0.722 for the interaction between I2 and 15, and 0.783 when AtI mediates the interaction. This trend holds for other sulphur bases as for the selenoether 7. Such finding might appear at first sight in contradiction with the stronger donor ability of AtI, but keep in mind that in vacuo calculated values are assessed with respect to crystal distances that endure the packing effects.
![]() | ||
Fig. 6 Relationship between the pKBAtI and pKBI2 scales for the 16 selected Lewis bases. The numbers in parenthesis in the analytical expression correspond to one standard deviation. |
Other interhalogen donors can be indirectly compared with the AtI donor via the pKBI2 scale. For instance, the pKBICl and pKBIBr scales compared to the pKBI2 scale give linear relationships expressed as pKBICl = 1.79 × pKBI2 + 1.42 and pKBIBr = 1.53 × pKBI2 + 0.87,15,16 which lie above the line formed by the pKBAtI scale (Fig. 6). We can deduce that the donor ability among interhalogens follows the order: I2 < AtI < IBr < ICl. The ordering between the iodine donors, I2, IBr and ICl, is in line with the general consensus that more electron-withdrawing groups R, bound to the electrophilic halogen atom X, will lead to more stable complexes of R–X with Lewis bases.2,42 Conversely, the position of AtI vs. IBr and ICl could not have been anticipated since the donor atom X and the electron-withdrawing group R are both different. Even the calculated values of VS,max for these XB donors are useless (see Fig. S3 in the ESI†). One may wonder if the pKBAtI scale encodes different information from that provided by other Lewis basicity scales that are already known, and especially hydrogen bonding scales. However, it has been shown that the pKBI2 scale, which is in overall linearly correlated with the pKBAtI scale, as seen above, is orthogonal or quasi-orthogonal to these scales.16 For instance, the 4-fluorophenol pKBHX hydrogen-bond basicity scale explains only 36% (on a sample of 265 points) of the variance of the pKBI2 scale. Although correlations with pKBAtI could be established within given families of compounds (e.g., oxygen functionalities, or their thio-derivatives), a comprehensive understanding of this behaviour would need specific investigations beyond the scope of this work.
Footnote |
† Electronic supplementary information (ESI) available: A description of the experimental and the computational methodologies, and many tables and figures reporting additional data. See DOI: 10.1039/d1sc02133h |
This journal is © The Royal Society of Chemistry 2021 |