4
O. Dada et al. / Tetrahedron Letters xxx (2018) xxx–xxx
of S-glycopeptides and lanthionine derivatives [27]. A mixture of
10 and 8 in EtOAc was treated with an aqueous solution of
K2CO3, and as expected, the desired complex 3 was produced in
very high yield (Scheme 1). Both conditions were thus applied to
the formation of compounds 4 and 5, and fortunately, all reactions
worked well with the highest-yielding shown in Scheme 1. Thio-
glycosylation of 8 with thiol 11 under the action of Et3N gave com-
pound 4 in 91% yield, while with thiol 12 the reaction provided
compound 5 in 86% yield under the biphasic conditions. All of
the NHC-Au(I) complexes 2–5 were isolated as off-white solids
and fully characterized by spectroscopic techniques and micro-
analysis. The formation of complexes 2–4 could be readily deter-
mined by comparing the 1H NMR spectra of the complexes and
their corresponding glycosyl thiols 9–11. A triplet peak corre-
sponding to the anomeric protons of the thiols disappeared and
was replaced with a doublet peak in the complexes. The chemical
shifts of these protons were also clearly shifted. Also, the benzylic
protons observed as a singlet peak in the precursor chloride 8,
became diastereotopic and showed doublet of doublet peaks in
complexes 3 and 4, which was in accordance with the previous
report [12]. Interestingly, the benzylic protons remain a singlet
It is worth noting that the aim of this study was not to perform
a thorough investigation on the relative stability of all the possible
conformations but to evaluate the C-Au and Au-S distances and
benchmark the different functionals against the crystal structure
in order to select the best computational level to be used in larger
systems. All the functionals utilised provided reliable parameters
when compared with the crystal structure (see ESI for details).
The calculated C-Au distances range between 2.012 and 2.030 Å.
The worst case corresponds to B3LYP while MN15 and WB97XD
presented a good agreement with respect to the actual crystal
structure. Similar results are found for the Au-S distances in which
MN15 and WB97XD provided the closest values with respect to the
crystal structure in both conformers 1 and 2. Values of the Au-S-C
angle are very close to the crystal structure with deviations of 1°.
Fig. 4 shows the overlap between the calculated models and the
crystal structure (blue). As observed, all the functionals present a
good overlap with no significant deviations. As previously
observed, B3LYP presented the largest variations (Fig. 3, red) with
a total overall RMSD of the position of each atom with respect to
the crystal structure of 0.674 Å. The total overall RMSD for the rest
of the functionals were 0.428 (M06-2X, orange), 0.501 (MN15,
green) and 0.381 Å (WB97XD, cyan). Considering these results,
MN15 and WB97XD provided the most reliable results (less total
RMSD and distance deviation), and reproduce the crystal structure
more accurately. However, M06-2X could also describe the struc-
ture within a reasonably good accuracy.
A similar benchmark study was also carried out using the crys-
tal structure for p-methoxythiophenolate conjugated to the NHC⁄-
Au complex [28]. As observed in the case of NHC-Au-SCN, WB97XD
and MN15 functional provide the most accurate results, i.e. struc-
tural parameters such as bonding distances and angles from those
functionals are the closest to the crystal structure (see ESI). Once
the functionals are benchmarked and the computational level
selected, we obtained a model structure for compound 4. Fig. 4
shows the optimised structure and all the selected structural
parameters were gathered in Table 3 (see ESI). It is shown that
the Au-S bond distance ranges from 2.404 to 2.339 Å with an aver-
age Au-S bond distance of 2.372 0.027 Å. The C-Au distance is less
sensitive to the method with an average distance across the
method of 2.027 0.007 Å. Finally, the Au-S-C angle is very stable
showing an average value of 94.5 0.4°. Of course, there is a wide
range of different conformers to be considered since the OH groups
can have free rotation. But based on the previous crystal structures
and the benchmarks, we believe that this can be a very reliable and
possible structure.
peak for complexes 2 and 5 which is likely due to the a-anomeric
configuration of the sugar in 2 and the extra three-carbon linker in
5. The most conspicuous change in the 13C NMR spectra was the
migration of the chemical shift of the carbene carbon (NCN). The
peak, observed at around 170 ppm for chloride 8, was shifted
downfield to around 180 ppm for the NHC-Au(I) complexes 2–5,
indicating the different chemical environment of the carbon before
and after the reaction.
The in vitro cytotoxicity of the four synthesized complexes 2–5
were subsequently carried out at the National Cancer Institute
(NCI), and the results are given in the ESI. In all tested cell lines,
complex 2 had better activity in comparison with the other tested
compounds and is therefore a potential new anticancer drug candidate.
Computational results
Computational calculations were carried out in order to provide
a theoretical structure of 4.
Structures were optimized by DFT using different functionals
(B3LYP, M06-2X, WB97XD and MN15) with the 6–31 + G(d) basis
set for all the atoms and pseudopotentials LANL2DZ on the heavy
Au atom as described in the computational methods. Firstly, the
performance of each functional was benchmarked on the crystal
structure of NHC-Au-SCN [28]. For such purpose, the crystal struc-
ture was used as an input (Isomer 1) and optimised without any
constraints. Additionally, a second isomer was considered to
explore the possible influence of the different arrangements on
the structural parameters, i.e., Au-S distances, angles, etc, within
both isomers (Fig. 2).
In summary, as an extension of our previous project, four new
monomeric NHC-Au(I)-sugar complexes 2–5 were synthesized in
very high to excellent yields by reacting the appropriate pre-syn-
thesized acetylated thiosugars with NHC-Au(I) chloride. Two
different conditions were attempted for the reaction and the
Fig. 2. Two different conformers evaluated using DFT. Conformer
1
(Conf.1)
Fig. 3. Overlap between the crystal structure (blue), and those calculated using
corresponds to the crystal structure.
B3LYP (red), M06-2X (orange), MN15 (green) and WB97XD (cyan) functionals.