key: cord-0708798-19rsqrqb authors: Di Gaetano, Sonia; Pirone, Luciano; Galdadas, Ioannis; Traboni, Serena; Iadonisi, Alfonso; Pedone, Emilia; Saviano, Michele; Gervasio, Francesco Luigi; Capasso, Domenica title: Design, Synthesis, and Anticancer Activity of a Selenium-Containing Galectin-3 and Galectin-9N Inhibitor date: 2022-02-25 journal: Int J Mol Sci DOI: 10.3390/ijms23052581 sha: 288adef77caf7141aba2834256928db8b9c9cb52 doc_id: 708798 cord_uid: 19rsqrqb Galectins are soluble β-D-galactoside-binding proteins whose implication in cancer progression and disease outcome makes them prominent targets for therapeutic intervention. In this frame, the development of small inhibitors that block selectively the activity of galectins represents an important strategy for cancer therapy which is, however, still relatively underdeveloped. To this end, we designed here a rationally and efficiently novel diglycosylated compound, characterized by a selenoglycoside bond and the presence of a lipophilic benzyl group at both saccharide residues. The relatively high binding affinity of the new compound to the carbohydrate recognition domain of two galectins, galectin 3 and galectin 9, its good antiproliferative and anti-migration activity towards melanoma cells, as well as its anti-angiogenesis properties, pave the way for its further development as an anticancer agent. Galectins are a family of soluble β-D-galactoside-binding proteins characterized by a highly conserved carbohydrate recognition domain (CRD) [1] and are expressed in a broad range of animal species. These proteins are present in various cell types, including epithelial and endothelial cells from the skin, lung, and gut, as well as immune cells and fibroblasts [2] [3] [4] . Based on their molecular structure, galectins have been classified into three groups: (a) proto-galectins (galectin 1, -2, -5, -7, -10, -11, -13, -14 and -15) with a single CRD; (b) tandem repeat galectins (galectin 4, -6, -8, -9 and -12) with two distinct CRDs connected by a linker peptide; and (c) chimeric galectins represented by only one member, galectin 3, composed of a canonical domain connected to a long tail [5] . Galectins, synthesized on the cytoplasm, translocate rapidly to the extracellular space through non-classical exocytic pathways (non-ER/Golgi secretory route) [6] , where they get exposed to a large variety of glycans. The binding of galectins to several carbohydrates attached to extracellular proteins and lipids and other membrane components leads to the formations of the so-called "galectin-glycoprotein lattice" [7] , which controls the concentration, localization, and availability of glycosylated receptors on the cell surface. In that way, galectins modulate the transmembrane signaling events of diverse physiological and pathological processes, including among others cell-cell and cell-extracellular matrix interactions, cell adhesion, proliferation, differentiation, as well as angiogenesis [8, 9] , autoimmune processes [10] , pulmonary fibrosis [11] , cancer progression, and metastasis [12] . Moreover, galectins control infections by mediating the endocytosis of different pathogens through continuous surveying of the extracellular environment [13] [14] [15] . The role of galectins in all these cellular pathways would imply lack of specificity, but this is far from being the case, as the structural complexity and variability of the carbohydrates that galectins recognize result in specific interactions that allow galectins to selectively activate different pathways. Besides extracellular functions, intracellular galectins bind to several intracellular substrates to regulate different biological processes. As they do not interact with glycans, their activity is predominantly mediated by protein-protein interactions, although several carbohydrate molecules, such as lactose, are capable of inhibiting intracellular actions [16] . Galectins, in particular Gal-1 and Gal-3, have been shown to influence cell signaling through association with proteins such β-catenin and members of the RAS family [17] . Intracellular Gal-3 interacts with Bcl-2 and is able to block the intrinsic apoptotic pathway [18] . High intracellular Gal-3 levels are also known to affect the inflammatory response through various mechanisms [19] . As to cancer biology, galectins are highly involved in the regulation of apoptosis, tumor metastasis, tumor cell adhesion, migration, and angiogenesis [20] . Interestingly, the expression of galectins is altered in tumor cells compared to their normal counterparts. In many situations, this modified galectin expression correlates with the aggressiveness of tumors and the acquisition of the metastatic phenotype, indicating that galectins might modulate tumor progression and influence the disease outcome [21, 22] . In particular, current literature data highlight the importance of the role of Gal-3 in melanoma progression [23, 24] . Gordon-Alonso et al. [23] demonstrated that in melanoma with a dismal prognosis, Gal-3 binds the glycans of the extracellular matrix and of IFNγ, decreasing T-cell recruitment, which is the best predictor to date of patient survival [23] . A previous study reported that in melanoma cells, extracellular Gal-3 is able to induce the MMP-9 expression, activating the p38 MAPK pathway involved in cellular processes such as matrix degradation [24] . Therefore, the development of potent and selective small inhibitors that block the activity of galectins is a promising therapeutic approach for cancer. Accordingly, several natural or synthetic carbohydrates, able to bind galectins and interfere with their metabolism confirm the therapeutic potential of galectin inhibition. Simple lactose and N-acetyllactosamine (LacNAc) exhibit binding affinities in the high micromolar range, but their multivalent derivatives show greater affinity [25] . However, the use of these natural carbohydrates as drugs is impaired by their sensitivity to glucosidases, which reduce drastically their half-life in vivo. For this reason, specific derivatizations aimed at increasing their metabolic stability and specificity have been considered [26, 27] . Among these, the presence of an alternative sulfur bridge in β-D-galactopyranosyl-(1-1)-thio-β-Dgalactopyranoside (TDG) made the molecule stable to hydrolysis while retaining the same affinity to galectins as lactose and paved the way for the production of small symmetrical inhibitors with favorable properties [28] . TDG derivatives carrying different substitutions such as 4-aryl-1,2,3-triazolyl (TD139) [29] , 3-benzamido [28] , 4-amido-1,2,3-triazolyl [30] , or 3-O-coumarylmethyl [31] moieties at the 3-and 3'-positions are able to bind Gal-1 and Gal-3 with nanomolar affinity. In particular, TD139 has recently been approved by the FDA for the treatment of idiopathic pulmonary fibrosis [31] and, 3,3 -O-dicoumarylmethyl substituted thiodigalactoside appears to be the first molecule capable of binding Gal-3 with selectivity over Gal-1 [31] . Recently, our research group has evaluated the antiproliferative effect of a series of sulfurated diglycosylated compounds against different human tumor cell lines [32] . In an attempt to establish a structure-activity relationship, we focused on the effect of both the saccharide residues and the number of bridging sulfur atoms (sulfide or disulfide structures) [32] on the antiproliferative properties and affinity of the different compounds against galectins. Here, we designed, synthesized, and characterized a new diglycosylated compound as a potential antitumor agent, characterized by a selenoglycoside bond and a suitably placed aromatic benzyl group at C-3 of both sugar residues. The choice of selenium as the bridging atom in place of sulfur was based on its chemical properties [33] . Selenium is an essential element for human metabolism and is a component of selenoamino acids that are necessary for numerous cellular functions. Organoselenides always sparked scientific interest for their anticancer and anti-inflammatory functions due to the duality of selenium which, unlike the lighter chalcogen, sulfur, has both antioxidant and peroxidant properties [33] . Ebselen, a glutathione peroxidase mimetic [34] , was reported as anticancer [35, 36] , antiinflammatory, antidepressant, and antibacterial [37] , among others, and recently discovered as a potent inhibitor of the main protease from SARS-CoV-2 [38] . The poor solubility of the organoselenides that reduces their availability can be easily overcome by selenoglycosides [39, 40] , but the biological assessment of the latter has been limited to date by the need for expensive synthetic procedures. As a matter of fact, poor data are available on the bioactivity of selenoglycosides or 3,3 -substituted selenoglycosides derivatives covered by patent [41] , as essentially are focused only on the assessment of their affinity towards lectin receptors [42] rather than their biological effect A recently introduced approach provided us with a very useful tool for quick production of digalactosyl selenides starting from a cheap selenenating agent [43] , a method which is herein complemented by the availability of a simple synthetic strategy for the regioselective functionalization of carbohydrates with aromatic groups. The new compound herein described was evaluated for its ability to bind Gal-3 CRD and Gal-9N CRD as well as for its antiproliferative, anti-migration, and anti-angiogenesis activity. Its relatively high affinity, anti-migration and anti-angiogenesis properties pave the way for further development of such compounds as antitumor agents. Among the thiogalactosides described in the literature, we selected the simplest βgalactoside sugar, thiodigalactoside, and its selenated counterpart as lead compounds and we tried to further optimize their structures. Derivatization at position C3 has been proposed to improve the affinity for galectins through favorable interactions with charged residues [28] . Although derivatization at positions C4 and C6 is feasible, the hydrogen bond formed between the β-galactoside scaffold and galectins ( Figure 1 ) indicate that these interactions not only stabilize the β-galactoside in the binding site but are expected to be necessary for substrate recognition. Consequently, we have focused on potential modifications only at position C3. After comparison of several crystal structures of the CRD of Gal-3, and Gal-9N bound to different ligands, and given the presence of many arginines in the binding groove of the two galectins, we anticipated that arene substitutions in position C3 could enhance the affinity against specific galectins. Indeed, molecular docking of the digalactosyl selenide, SeDG, and its benzyl 3,3 -derivative, SeDG-Bn, to Gal-3, and Gal-9N indicates that the latter could bind more strongly to either galectin. In particular, in addition to the hydrogen bonds that the hydroxyl groups of the galactoside core form with the side chains of H158, N174, R162, E184, and W181, the presence of the phenyl rings in SeDG-Bn allows the formation of π-cation interactions with the side chains of R144 and R186 ( Figure 1) . Moreover, the relatively flexible benzyloxy group of SeDG-Bn is expected to be able to adopt an alternative conformation in both Gal-3 CRD , and Gal-9N CRD , in which the phenyl ring forms a π-π interaction with W181. The enhanced binding of SeDG-Bn to Gal-3 CRD , and Gal-9N CRD due to the additional interactions is reflected to its better docking score (−7.92 and −7.72 kcal/mol, respectively) with respect to SeDG (−6.75 and −6.72 kcal/mol, respec-tively). Although it is known that the docking score of a compound alone does not necessarily correlate with its activity, the predicted binding pose of SeDG-Bn with respect to SeDG in combination with its better docking score prompted us to synthesize and evaluate the in vitro activity of SeDG-Bn. After comparison of several crystal structures of the CRD of Gal-3, and Gal-9N bound to different ligands, and given the presence of many arginines in the binding groove of the two galectins, we anticipated that arene substitutions in position C3 could enhance the affinity against specific galectins. Indeed, molecular docking of the digalactosyl selenide, SeDG, and its benzyl 3,3'-derivative, SeDG-Bn, to Gal-3, and Gal-9N indicates that the The synthesis of SeDG 10 started from the conversion of pentacetyl galactose 1 into the aryl glycoside 2, through a recently described solvent-free glycosidation catalyzed by methansulfonic acid (Scheme 1) [44] . In this step, we introduced a para-methoxy-aryl as a suitable anomeric protecting group, which allowed the selective manipulation of the remaining reactive functionalities of the sugar. Compound 2 was smoothly deacetylated under standard Zemplén conditions (catalytic sodium methoxide in methanol) and the resulting tetrol 3 was submitted to a regioselective benzylation at O-3, following a recently established solvent-free methodology allowing the catalytic use of dibutyltin oxide [45] . The reaction provided 3-O-benzylated triol 4 in a good yield and excellent regioselectivity. After chromatographical purification, this compound was smoothly peracetylated under standard conditions, and the resulting derivative 5 was then deprotected at O-1 under oxidative conditions with cerium ammonium hexanitrate (CAN) in an acetonitrile/water mixture. The resulting anomeric hemiacetal 6 was acetylated with sodium acetate and acetic anhydride to yield compound 7 as an approximately equimolar anomeric mixture. The key step of the sequence was the synthesis of the diglycosyl selenide, which was performed by adapting a recently developed procedure, applied to date only to peracetylated precursors [43] . To this aim, compound 7 was first iodinated at the anomeric position with the iodine/silane combined reagent [45, 46] . It is worth noting that such a reaction is usually performed at 40 • C (refluxing dichloromethane) with a peracetylated precursor, but the presence of an activating benzyl group made the substrate expectedly more reactive and, thus, a lower temperature could be used. The use of relatively milder conditions allowed us to minimize the undesired occurrence of a HI-mediated benzyl cleavage [47] . Glycosyl iodide 8 was isolated as a crude product after an extractive work-up and directly treated (Scheme 2) with a mixture of selenium and sodium borohydride, preliminarily kept under stirring for 40 min to favor the reduction of elementary selenium. In this step, triethyl phosphite was added to guarantee the chemoselectivity of the reaction, by minimizing the generation of the corresponding digalactosyl diselenide. The proce-Scheme 1. Synthesis of the key-intermediate 8. Glycosyl iodide 8 was isolated as a crude product after an extractive work-up and directly treated (Scheme 2) with a mixture of selenium and sodium borohydride, preliminarily kept under stirring for 40 min to favor the reduction of elementary selenium. In this step, triethyl phosphite was added to guarantee the chemoselectivity of the reaction, by minimizing the generation of the corresponding digalactosyl diselenide. The procedure provided selenide 9 in a satisfying overall yield (33% from 7). After chromatographical purification, 9 was deacetylated under Zemplén conditions. Interestingly, application of the standard conditions (catalytic sodium methoxide at RT) led only to O-4 and O-6 deacetylation, and deprotection at the more congested O-2 site required more forced conditions in terms of both temperature and amount of NaOCH 3 . The deprotected selenide was eventually obtained in a satisfying yield (72%) after silica gel chromatography. All synthetic intermediates and the final product were characterized by NMR ( Figure S1 ). Scheme 1. Synthesis of the key-intermediate 8. Glycosyl iodide 8 was isolated as a crude product after an extractive work-up and directly treated (Scheme 2) with a mixture of selenium and sodium borohydride, preliminarily kept under stirring for 40 min to favor the reduction of elementary selenium. In this step, triethyl phosphite was added to guarantee the chemoselectivity of the reaction, by minimizing the generation of the corresponding digalactosyl diselenide. The procedure provided selenide 9 in a satisfying overall yield (33% from 7). After chromatographical purification, 9 was deacetylated under Zemplén conditions. Interestingly, application of the standard conditions (catalytic sodium methoxide at RT) led only to O-4 and O-6 deacetylation, and deprotection at the more congested O-2 site required more forced conditions in terms of both temperature and amount of NaOCH3. The deprotected selenide was eventually obtained in a satisfying yield (72%) after silica gel chromatography. All synthetic intermediates and the final product were characterized by NMR ( Figure S1 ). The genes encoding Gal-3 CRD (amino acid residues 112-250) and Gal-9N CRD (amino acid residues 1-148) were cloned into the vector pETM11. Escherichia coli BL21(DE3)GOLD cells were transformed with recombinant constructs and the expression of the recombinant proteins were induced with 0.5 mM isopropyl β-D-1-thiogalactopyranoside (IPTG) for 16 h at 22 • C. Cells were harvested by centrifugation, suspended in 20 mM Tris⁄HCl, 500 mM NaCl, and 1 mM dithiothretiol pH 7.5. The proteins were purified by a two-step procedure consisting of a His-trap affinity chromatography and a gel filtration. The proteins were concentrated using a Centricon-10 ultrafiltration unit (Amicon, Merck, Darmstadt, Germany). The final yield was 20 mg/mL and 2 mg/mL for Gal-3 CRD and Gal-9N CRD , respectively. Gal-3 CRD and Gal-9N CRD , once expressed as recombinant proteins in a prokaryotic system and purified to homogeneity, were used to evaluate the binding of lactose, which is their natural substrate, to verify their functional integrity. Both proteins were fully active as they were able to bind lactose (Supplementary Figure S2) . Both purified proteins were then probed for their binding affinity to SeDG and SeDG-Bn by ITC. The Kd values obtained for Gal-3 CRD and Gal-9N CRD to SeDG resulted to be of 20.2 ± 2. Table 1 ). On the contrary, the Kd values of SeDG-Bn were 2.2 ± 0.45 and 2.9 ± 0.3 µM for Gal-3 CRD and Gal-9N CRD , respectively, indicating the higher affinity of SeDG-Bn against both proteins ( Figure 3 and Table 1 ). fully active as they were able to bind lactose (Supplementary Figure S2) . Both purified proteins were then probed for their binding affinity to SeDG and SeDG-Bn by ITC. The Kd values obtained for Gal-3 CRD and Gal-9N CRD to SeDG resulted to be of 20.2 ± 2.2 µ M and 22.4 ± 4.4 μM, respectively (Figure 2A-C and Table 1 ). On the contrary, the Kd values of SeDG-Bn were 2.2 ± 0.45 and 2.9 ± 0.3 µ M for Gal-3 CRD and Gal-9N CRD , respectively, indicating the higher affinity of SeDG-Bn against both proteins ( Figure 3 and Table 1 ). Interaction studies between SeDG and Gal-3 CRD or Gal-9N CRD by ITC. Titration of (A) Gal-3 CRD and (B) Gal-9N CRD with SeDG. Titration of (C) buffer with SeDG-Bn (negative control). The upper panels of the graphs depict the injections, while the lower panels depict the integrations of the peaks according to the molar ratio of the molecules tested. The upper panels of the graphs depict the injections, while the lower panels depict the integrations of the peaks according to the molar ratio of the molecules tested. The resulting values of Kd, ∆G, ∆H, −T∆S and n (stoichiometry) are summarized in Table 1 . In all cases, the n-values indicate a 1:1 ratio of galectin/inhibitor. The negative ∆H suggests an exothermic reaction between the galectins and the ligands, while the positive ∆S suggests entropically driven reactions involving hydrophobic interactions, especially in the case of SeDG-Bn. The introduction of the benzyl group at position C3 of SeDG-Bn (Kd~2 µM) leads to a~10-fold enhancement of the binding affinity for both proteins compared to the SeDG (Kd~20 µM), underlying the importance of the substitution. The enhanced binding affinity of SeDG-Bn to Gal-3 CRD and Gal-9N CRD with respect to SeDG is reflected in a small change in the binding free energy (∆∆GSeDG-Bn/SeDG of −1.38 kcal/mol and -0.4 kcal/mol, respectively). The resulting values of Kd, ΔG, ΔH, −TΔS and n (stoichiometry) are summarized in Table 1 . In all cases, the n-values indicate a 1:1 ratio of galectin/inhibitor. The negative ΔH suggests an exothermic reaction between the galectins and the ligands, while the positive ΔS suggests entropically driven reactions involving hydrophobic interactions, especially in the case of SeDG-Bn. The introduction of the benzyl group at position C3 of SeDG-Bn (Kd ~ 2 μM) leads to a ~10-fold enhancement of the binding affinity for both proteins compared to the SeDG (Kd ~20 μM), underlying the importance of the substitution. The enhanced binding affinity of SeDG-Bn to Gal-3 CRD and Gal-9N CRD with respect to SeDG is reflected in a small change in the binding free energy (ΔΔGSeDG-Bn/SeDG of −1.38 kcal/mol and -0.4 kcal/mol, respectively). As previously reported, SeDG is unable to affect tumor cell proliferation even at 2 mM concentration [43] . In order to investigate the cellular effect of the aromatic derivatization at position C3 of both galactoside residues of SeDG-Bn, we performed a cytotoxicity assay on cell lines that differ in metastatic potential or histological origin. In particular, for this assay, we used a metastatic melanoma cell line (WM266), as well as low metastatic melanoma (A375), cervix carcinoma (HeLa), and human T lymphoblastoid cells (Jurkat). To verify the selectivity of action against cancer cells, we used the human dermal fibroblasts (HDF) as a healthy cellular model. First, we analyzed the cells to evaluate the presence of Gal-3 and Gal-9 and the results showed that all of them present significant As previously reported, SeDG is unable to affect tumor cell proliferation even at 2 mM concentration [43] . In order to investigate the cellular effect of the aromatic derivatization at position C3 of both galactoside residues of SeDG-Bn, we performed a cytotoxicity assay on cell lines that differ in metastatic potential or histological origin. In particular, for this assay, we used a metastatic melanoma cell line (WM266), as well as low metastatic melanoma (A375), cervix carcinoma (HeLa), and human T lymphoblastoid cells (Jurkat). To verify the selectivity of action against cancer cells, we used the human dermal fibroblasts (HDF) as a healthy cellular model. First, we analyzed the cells to evaluate the presence of Gal-3 and Gal-9 and the results showed that all of them present significant amounts of one or both proteins (Supplementary Figure S3 ) [32] . We then incubated the cells in the presence of the SeDG and SeDG-Bn for 48 h with concentrations ranging from 1 to 800 µM. The most interesting result has been obtained with WM266 cells, in which the dose-effect curve ( Figure 4 ) provides an IC50 of 430 µM for SeDG-Bn compared to 990 µM in healthy HDF cells demonstrating a selective activity towards cancer cells even though it ends up on the high micromolar range. Furthermore, very high values of IC50 (>900 µM) were obtained on other tested tumor cell lines treated with SeDG-Bn, demonstrating really scarce interference of SeDG-Bn with A375, HeLa, and Jurkat cell proliferation. The selective activity of SeDG-Bn towards WM266 cells suggests that SeDG-Bn is a good candidate for further optimization in melanoma therapy. Several studies report that non-cytotoxic galectin inhibitors are able to interfere with galectin-mediated cancer development [48, 49] . In in vivo experiments with nude mice, galectin inhibitors significantly reduced Gal-3 mediated lung metastasis [49] . Considering these previous data and encouraged by the bioactivity of SeDG-Bn on WM266 cells, we chose this tumor cell line to perform further biological assays. In particular, we evaluated the ability of SeDG and SeDG-Bn to inhibit WM266 cell migration and invasion, as well as angiogenesis progression carried out by means of endothelial cells [50, 51] . μM in healthy HDF cells demonstrating a selective activity towards cancer cells even though it ends up on the high micromolar range. Furthermore, very high values of IC50 (>900 µ M) were obtained on other tested tumor cell lines treated with SeDG-Bn, demonstrating really scarce interference of SeDG-Bn with A375, HeLa, and Jurkat cell proliferation. The selective activity of SeDG-Bn towards WM266 cells suggests that SeDG-Bn is a good candidate for further optimization in melanoma therapy. Several studies report that non-cytotoxic galectin inhibitors are able to interfere with galectin-mediated cancer development [48, 49] . In in vivo experiments with nude mice, galectin inhibitors significantly reduced Gal-3 mediated lung metastasis [49] . Considering these previous data and encouraged by the bioactivity of SeDG-Bn on WM266 cells, we chose this tumor cell line to perform further biological assays. In particular, we evaluated the ability of SeDG and SeDG-Bn to inhibit WM266 cell migration and invasion, as well as angiogenesis progression carried out by means of endothelial cells [50, 51] . Firstly, we performed an in vitro wound healing assay on WM266. To this aim, cell confluent monolayers were scratched and SeDG was added at 1 and 2.5 mM and SeDG-Bn at 100 and 200 µ M concentrations. To avoid the cytotoxic effect of SeDG-Bn, in motility and invasion assays, we used concentrations lower than the IC50. The results demonstrate that the wound healing of treated cells was significantly delayed compared to untreated cells. In particular, both compounds decreased the wound closure at all used concentrations and time intervals ( Figure 5A,B) . Notably, SeDG-Bn was able to inhibit cell migration at concentrations of an order of magnitude lower than SeDG did. Firstly, we performed an in vitro wound healing assay on WM266. To this aim, cell confluent monolayers were scratched and SeDG was added at 1 and 2.5 mM and SeDG-Bn at 100 and 200 µM concentrations. To avoid the cytotoxic effect of SeDG-Bn, in motility and invasion assays, we used concentrations lower than the IC50. The results demonstrate that the wound healing of treated cells was significantly delayed compared to untreated cells. In particular, both compounds decreased the wound closure at all used concentrations and time intervals (Figure 5A,B) . Notably, SeDG-Bn was able to inhibit cell migration at concentrations of an order of magnitude lower than SeDG did. Under the same experimental conditions, the specificity of SeDG-Bn for Gal-3 CRD was validated using Gal-3 CRD as a cell motility stimulator in the absence or presence of SeDG-Bn ( Figure 5C ). Given the reported implication of Gal-3 CRD in cell migration [52] , our results confirm that cells treated with 50 µM Gal-3 CRD migrate more than in the control. Noteworthy, Gal-3 CRD -mediated wound closure was completely inhibited by SeDG-Bn treatment ( Figure 5C ). Such observation provides experimental evidence that SeDG-Bn can reduce cell migration by inhibiting the stimulating effect of Gal-3 CRD . Then, we analyzed the ability of compounds to interfere with the invasion of WM266 cells throughout the extracellular matrix. The results obtained showed a marked inhibition of the WM266 invasion induced by the addition of SeDG-Bn acting in a dose-dependent mode ( Figure 6 ). Moreover, SeDG was able to reduce the cells invasion but only at millimolar concentrations. Taken together, these results demonstrate that SeDG-Bn is able to interfere specifically with essential Gal-3 CRD -driven cellular mechanisms involved in the metastatic cell diffusion, such as mobility and invasion of cancer cells. Galectins have been shown to protect the cells from oxidative stress due to oxygen lack within the tumor microenvironment and act as pro-angiogenic receptors stimulating the growth of vascular endothelial cells [53] [54] [55] . Considering that melanoma requires high angiogenic activity [56] , we evaluated the antiangiogenic activity of SeDG-Bn. To this aim, we used a common model of angiogenesis in which human umbilical vein cells (HuVeC) are layered on an appropriate matrix-forming cellular network (mesh-like structures). This model allowed us to verify the ability of digalactosyl selenides to inhibit cell vasculogenic capacity. Notably, SeDG-Bn was able to drastically prevent the formation of vessels at all of the evaluated concentrations in a dose-dependent manner, resulting in up to 60% reduction (Figure 7) . In contrast, we found that SeDG at 2.5 mM results in a non-statistically significant change of capillary tube formation. Under the same experimental conditions, the specificity of SeDG-Bn for Gal-3 CRD was validated using Gal-3 CRD as a cell motility stimulator in the absence or presence of SeDG-Bn ( Figure 5C ). Given the reported implication of Gal-3 CRD in cell migration [52] , our results confirm that cells treated with 50 µ M Gal-3 CRD migrate more than in the control. Noteworthy, Gal-3 CRD -mediated wound closure was completely inhibited by SeDG-Bn The distance between the limits of the scratch was quantified and the size of wound closure was determined as follows: % wound closure = 1 − (wound width tx/wound width t0) × 100. Bars depict mean ± SE of three independent experiments, ** p < 0.01 * p < 0.05. tion of the WM266 invasion induced by the addition of SeDG-Bn acting in a dose-dependent mode ( Figure 6 ). Moreover, SeDG was able to reduce the cells invasion but only at millimolar concentrations. Taken together, these results demonstrate that SeDG-Bn is able to interfere specifically with essential Gal-3 CRD -driven cellular mechanisms involved in the metastatic cell diffusion, such as mobility and invasion of cancer cells. Galectins have been shown to protect the cells from oxidative stress due to oxygen lack within the tumor microenvironment and act as pro-angiogenic receptors stimulating the growth of vascular endothelial cells [53] [54] [55] . Considering that melanoma requires high angiogenic activity [56] , we evaluated the antiangiogenic activity of SeDG-Bn. To this aim, are layered on an appropriate matrix-forming cellular network (mesh-like structures). This model allowed us to verify the ability of digalactosyl selenides to inhibit cell vasculogenic capacity. Notably, SeDG-Bn was able to drastically prevent the formation of vessels at all of the evaluated concentrations in a dose-dependent manner, resulting in up to 60% reduction (Figure 7) . In contrast, we found that SeDG at 2.5 mM results in a nonstatistically significant change of capillary tube formation. Prior to docking calculations, the crystal structures of Gal-3 CRD (PDB ID 5H9P) and Gal-9N CRD isoform (PDB ID 3WLU) were prepared using the Protein Preparation Wizard of the Schrödinger suite [57] in order to add bond orders and formal charges to the starting protein structures. The protonation state of each residue was determined using the PROPKA algorithm at pH 7.0 as implemented in Maestro. Crystallographic water molecules that formed less than two hydrogen bonds with non-water atoms were removed, and the structure of each protein was relaxed using the OPLS3 force field [58] . Specifically, in the case of Gal-3CRD, the water molecule that was found to mediate an interaction between TD139 and W181 of Gal-3 CRD in the crystal structure (PDB ID 5H9P), and which was present also in the same position in other Gal-3 CRD structures (PDB IDs 3AYE, 4JC1, 4JCK), was considered to have a structural role and was maintained. Similarly, in the case of Gal-9NCRD, the water molecule that was found to Prior to docking calculations, the crystal structures of Gal-3 CRD (PDB ID 5H9P) and Gal-9N CRD isoform (PDB ID 3WLU) were prepared using the Protein Preparation Wizard of the Schrödinger suite [57] in order to add bond orders and formal charges to the starting protein structures. The protonation state of each residue was determined using the PROPKA algorithm at pH 7.0 as implemented in Maestro. Crystallographic water molecules that formed less than two hydrogen bonds with nonwater atoms were removed, and the structure of each protein was relaxed using the OPLS3 force field [58] . Specifically, in the case of Gal-3CRD, the water molecule that was found to mediate an interaction between TD139 and W181 of Gal-3 CRD in the crystal structure (PDB ID 5H9P), and which was present also in the same position in other Gal-3 CRD structures (PDB IDs 3AYE, 4JC1, 4JCK), was considered to have a structural role and was maintained. Similarly, in the case of Gal-9NCRD, the water molecule that was found to mediate an interaction between R44, N63 and the co-crystallized ligand (PDB ID 3WLU) and that was found to mediate the same interaction in other Gal-9N CRD structures (PDB IDs 2D6N, 2EAK) was maintained as well. Se is not an element found frequently in drug-like molecules and, therefore, there are no standard parameters for it in many popular force fields, including OPLS3 which is the force filed used in the docking calculations of this study. Therefore, to get the necessary bonded and non-bonded parameters for the docking, we performed first principles calculations (with density functional) using the Force Field Builder tool of the Schrödinger suite, prior to the docking calculations. For instance, Force Field Builder ran density functional theory (DFT) calculations with the B3LYP/cc-pVTZ functional on SeDG and SeDG-Bn, which were then fitted to obtain the parameters for the torsions that include the Se atom, that are not part of the standard OPLS3 force field. The torsional energies were calculated by performing coordinate scans with Jaguar as implemented in the Schrödinger suite, followed by fitting the resulting energy curves to produce new force field parameters. To dock the SeDG and SeDG-Bn compounds to Gal-3CRD and Gal-9N CRD , we used the extra precision (XP) scoring function and docking protocol of Glide 7.6 (as implemented in the Schrödinger suite) [59] . Glide searches for favorable interactions between a ligand and a receptor. In our study, Glide was run with flexible ligand conformations, where multiple conformations of each ligand were generated before docking. Given the increased number of rotatable bonds in the proposed compounds and the conserved binding mode of the sugar rings that were found to be placed in the same position in all of the examined crystal structures for each system, we used positional constraints to enforce docked poses to place the sugar rings in a predefined volume of space relative to each receptor. Namely, we used the co-crystallized ligand of each receptor and defined a sphere of 1.8 Å into which the center of mass of each sugar ring would be forced to be placed. Moreover, although each receptor was treated as a rigid body during the docking process, the hydroxyl groups of S237 in Gal-3CRD, and S135 in Gal-9N CRD were allowed to adopt different orientations with different ligands to produce the most favorable interactions. The introduction of these assumptions (i.e., positional constraints on the sugar rings and rotatable hydroxyl groups) was necessary for Glide to be able to reproduce the crystallographic poses of the co-crystallized ligand for each system. Figure 1 was generated using the ChimeraX visualization software [60] . Acetic anhydride (1.5 mL) was added to a solution of compound 2 in (350 mg, 0.93 mmol) in pyridine (3 mL). The mixture was left overnight at RT and then methanol (ca 1 mL) was added to the mixture kept in a cool bath. The mixture was diluted with dichloromethane and transferred into a separatory funnel. The organic phase was washed with sat sodium carbonate and the aqueous phase re-extracted with dichloromethane. Combined organic phases were dried with anhydrous sodium sulfate, filtered, and concentrated in vacuo to yield pure compound 5 as an oil (600 mg, quantitative yield). 1 13 Ceric ammonium hexanitrate (CAN) (340 mg, 1.44 mmol) was added to a solution of compound 5 (312 mg, 0.48 mmol) in 4:1 acetonitrile/water (5 mL). The mixture was kept under stirring at RT for 30 min, then diluted with ethyl acetate and transferred into a separatory funnel, washing the vessel with additional amounts of water and ethyl acetate. The organic phase was separated and the aqueous phase re-extracted with further ethyl acetate. The combined organic phases were dried with anhydrous sodium sulfate, filtered, and concentrated in vacuo. The residue was purified by silica-gel flash chromatography (eluent: hexane/ethyl acetate from 1:1 to 0:1) to yield hemiacetal compound 6 as an oil (150 mg, yield 78%). Anomeric mixture α/β ca 5:1. 1 Acetic anhydride (0.7 mL) was added to a vessel containing hemiacetal 6 (140 mg, 0.36 mmol) and sodium acetate (10 mg, 0.060 mmol).The mixture was kept under stirring at 120 • C (oil bath) for 30 min, then diluted with dichloromethane and transferred into a separatory funnel. The organic phase was washed with sat aq carbonate, and the aqueous phase re-extracted with further dichloromethane. The combined organic phases were dried with anhydrous sodium sulfate, filtered, and concentrated in vacuo to yield compound 7 (152 mg, 97%) directly used in the following step. 1 To a solution of 7 (145 mg, 0.33 mmol) and I 2 (67 mg, 0.26 mmol) in anhydrous DCM (4 mL), triethylsilane (42 µL, 0.26 mmol) was added at 0 • C. After 5 min, TLC analysis of the mixture (eluent: hexane/ethyl acetate 3:2) proves the nearly complete conversion of 7 into a less polar and very UV-visible product. After adding solid NaHCO 3 (70 mg), the cool mixture was diluted with dichloromethane and the organic phase was washed with aq. sodium carbonate containing sodium thiosulfate (the latter is added portionwise as a solid until consumption of residual iodine indicated by discoloration of the organic phase upon shaking). The organic phase was then washed with water, dried, and concentrated. The crude residue was directly adopted for the subsequent selenoglycosidation steps. A mixture of elementary selenium (32 mg, 0.40 mmol) and sodium borohydride (45 mg, 1.2 mmol) was suspended in DMF (1 mL) and the suspension was kept under stirring at RT. After 40 min, triethyl phosphite (58 µL, 0.40 mmol) was added (with an instantaneous discoloration) and the mixture was poured into a vessel containing crude glycosyl iodide 8 (0.33 mmol). The vessel adopted for the first step was washed with portionwise DMF (1 mL overall). The mixture was kept under stirring at room temperature for one hour. The reaction was quenched with acetic acid (0.3 mL), and the mixture was transferred into a separatory funnel and diluted with dichloromethane. The organic phase was washed with aq sodium carbonate and the aqueous phase re-extracted with dichloromethane. Combined organic phases were dried with sodium sulfate, filtrated and concentrated. The residue was purified by silica-gel flash chromatography (eluent: hexane/ethyl acetate from 3:2 to 0:1) to yield compound 9 (46 mg, 33% yield). 1 14, 2.04, 2.00 (9H, 3 × s, 3 × -COCH 3 ). 13 To a solution of compound 9 (24 mg, 0.028 mmol) in MeOH (1 mL), a solution of sodium methoxide in methanol (5.6 µmol) was added, preliminarily prepared by adding a weighted amount of NaH (60% suspension) in methanol. The mixture was left under stirring at RT until the detection via TLC (eluent: chloroform/methanol 9:1) of a single product. The mixture was treated with Amberlyst H+ resin (preliminarily washed with methanol) until neutrality. The supernatant was transferred to another vessel, where it was concentrated in vacuo. NMR of the residue evidenced a large prevalence of a partially deprotected product still acetylated at O-2. The residue was then dissolved again with methanol and treated with sodium methoxide in methanol (0.028 mmol) 75 • C for 30 min. The mixture was treated with Amberlyst H+ resin (preliminarily washed with methanol) until neutrality. The supernatant was transferred to another vessel, where it was concentrated in vacuo. The residue was purified by silica-gel flash chromatography (eluent: chloroform/methanol from 10:1 to 8:1) to yield pure SeDG-Bn (12 mg, 72% yield). 1 Genes encoding the CRD of Gal-3 and Gal-9N, corresponding to residues 112-250 and 1-148, respectively, were purchased from Genewiz. The two proteins were cloned into pETM11 expression vector with a N-terminal polyHis-tag that is removable by TEV protease. The recombinant vectors pETM11-gal3 CRD and pETM11-gal9 CRD were used to transform E. coli BL21(DE3)GOLD cells and a single colony was grown overnight at 37 • C in 20 mL of LB medium. In total, 10 mL of the saturated cultures were transferred in 1 L of LB medium and grown at 37 • C to 0.8 OD 600 nm before induction with IPTG (0.5 mM). The inductions were kept for 16 and for 3 h at 22 • C in the expression of Gal-3 CRD and Gal-9N CRD , respectively. The cells were harvested by centrifugation, re-suspended in 30 mL of lysis buffer (20 mM NaP, 500 mM NaCl, 10 mM imidazole pH 7.0 for Gal-3 CRD and pH 7.5 Gal-9N CRD ), and sonicated. After centrifugation (30,000× g, 30 min, 4 • C) the soluble fraction was subjected to affinity chromatography onto an HisTrapTM column (5 mL, GE Healthcare) previously equilibrated in lysis buffer. Elution was carried out with 300 mM imidazole. Fractions containing the recombinant protein were, finally, applied onto a Superdex 75 10/300 column (GE Healthcare), connected to an AKTA system (GE Healthcare) pre-equilibrated in 20 mM NaP, 150 mM NaCl pH 7.0 for the purification of Gal-3 CRD and pH 7.5 for Gal-9N CRD . The purified proteins were used to perform interaction studies. The interaction between the molecule SeDG-Bn with Gal-3 CRD or Gal-9N CRD was tested by isothermal titration calorimetry (ITC) experiments, which were carried out on an ITC200 calorimeter (MicroCal/GE Healthcare) at 25 • C. SeDG or SeDG-Bn (500 µM) were titrated into a solution of Gal-3 CRD or Gal-9N CRD (20 µM) . Before all titrations, molecules were dissolved in 20 mM NaP, 150 mM NaCl pH 7.0 (when the molecules were assayed against Gal-3 CRD ) or pH 7.5 (for Gal-9N CRD ). The sample cell was filled with 280 µL of Gal-3 CRD or Gal-9N CRD . Molecules were injected in a volume of 0.2 µL for the first injection and 1 µL for the next 39 injections using an injection interval of 150 s. The control experiment, obtained by injection of the same amount of each molecule in buffer solution, was subtracted to correct for heat production upon molecules dilution. All data were analyzed and fitted using the Microcal Origin version 7.0 software package. Binding enthalpy, dissociation constants, and stoichiometry were determined by fitting the data using a one-set-of-site-binding model [61] . Human adenocarcinoma cell line (HeLa), human low metastatic melanoma (A375), human metastatic melanoma (WM266), and human T lymphoblastoid cell line (Jurkat) were from ATCC (Manassas, VA, USA); normal human dermal fibroblasts (HDF) were kindly provided by Dr Annalisa Tito (Arterra, Biosciences, Naples, Italy). HeLa, A375, and HDF were grown in DMEM supplemented with 10% fetal bovine serum (FBS), 1% glutamine, 100 U/mL penicillin, and 100 µg/mL streptomycin (Euroclone, Milano, Italy). Under the same experimental conditions, WM266 and Jurkat were grown in RPMI supplemented with heat-inactivated 10% fetal bovine serum (FBS), 2.5 mM glutamine, 100 U/mL penicillin, and 100 µg/mL streptomycin (Euroclone). Human umbilical vein endothelial cells (HuVec) were purchased from Lonza. They were grown in EBM2 basal medium enriched with Single Quot Supplements. All the cells were maintained in humidified air containing 5% CO 2 , at 37 • C. Cell lysates were obtained incubating human tumor cells in a lysis buffer (50 mM Hepes pH 7.4, 50 mM NaCl, 1% Triton) supplemented with protease inhibitor cocktail (Roche, Basilea Switzerland) for 30 min on ice and centrifuged at 13,000 rpm for 30 min. The protein concentrations were determined by Bradford method using Bio-Rad reagent (Bio-Rad Laboratories, Hercules, CA, USA). The obtained lysates were loaded on a SDSpolyacrylamide gel electrophoresis and transferred to a PVDF membrane (Bio-Rad). Successively, the membrane was incubated with primary antibodies (anti-Gal-3 and anti-Gal-9 from R&D, anti-actin from Sigma Aldrich, St. Louis, MO, USA) 1 h at RT. Actin was used as an internal standard. Bands were visualized with LiteAblot (Euroclone, Pero, Italy) chemiluminescence as substrate and the signal was acquired with ChemiDoc XRS System (Bio-Rad) and analyzed with the QuantityONE software. Cells were plated at a density of 1200 cells/well for HeLa, 1000 for A375, 10,000 for Jurkat and 2000 for WM266 and HDF in 96-well microplates (Thermofisher, Waltham, MA, USA). After 24 h incubation, cells were treated with increasing concentrations of synthetized compounds for 48 h. Cell proliferation was determined by using MTT (Sigma Aldrich, St. Louis, MO, USA) for adherent cells [62] . In the case of Jurkat that are suspension lymphoblasts, the proliferation was investigated by CCK-8 assay [63] . Plates were then analyzed by using a microplate reader (Enspire, Perkin Elmer, Waltham, MA, USA) at 570 nm for MTT and at 450 for CCK-8. The mean value of proliferating cells for each treatment was compared to untreated cells (control). All experiments were performed at least in triplicate and repeated at least 3 times. The IC50 values were calculated using Graph-Pad Prism (version 6.01, Graph-Pad software Inc., San Diego, CA, USA). Analysis of cell migration in vitro has been carried out by using the wound healing scratch assay. WM266 cells were grown until they reached a confluent monolayer. After, they were linearly scratched with a pipette tip to generate a cell-free space. The detached cells were removed with PBS and digalactosyl selenides at indicated concentrations or 50 µM Gal-3 CRD or 50 µM Gal-3 CRD pre-incubated with 200 µM SeDG-Bn for 5 min, were added to the adherent cells at different concentrations in RPMI, then cells were incubated at 37 • C and successively the scratch area was photographed at 0, 24, and 48 h [64] . The area related to each incision was photographed at the time 0, 24, and 48 h using a phase contrast optical microscope with a 10× magnification (Zeiss). The distance between the edges of the incisions was measured by AxioVision Rel. 4.8 software and its mean value was determined as follows: wound closure (%) = 1 − (wound width tx/wound width t0) × 100. Bars depict mean ± SE of three independent experiments. To analyze cell invasion, the surface of a filter membrane of a transwell insert was coated with extracellular matrix (ECL) (Millipore). Cells (5 × 10 4 ) suspended in 100 µL of serum-free medium were incubated with various concentrations of molecules and successively they were seeded onto the upper compartment of the transwell, whereas the lower chambers were filled with medium containing 10% FCS [65] . After 24 h, the cells were fixed with 10% formalin and stained with crystal violet. The remaining cells in the upper side of the transwell were removed and cells on the undersurface of the filters, the invaded cells, were visualized and counted from six randomly selected fields (10× magnification) under a phase-contrast microscope with AxioVision Rel. 4.8 software. Angiogenesis assay was performed by using "In Vitro Angiogenesis Assay Kit" (Merck, Burlington, MA, USA). Briefly, 100 µL of ECL were plated in a pre-cooled 96-well plate and incubated for 1 h at 37 • C. 10,000 cells (HuVec), suspended in EGM-2 [66] , incubated with digalactosyl selenides, and seeded onto solidified matrix, according to manufacturer's instructions. After 5-6 h of incubation at 37 • C, the tube formation was examined by means of a phase contrast optical microscope (Axiovert 200 Zeiss) at 10× magnification and the branch point number was calculated from 3 random fields per each treatment. Imagines were acquired by AxioVision Rel. 4.8 software. All data are presented as mean values ± SE. The statistical analysis was performed using Student's t-test, unpaired, two-sided, p < 0.05 was considered significant. Although the involvement of galectins in many diseases such as diabetes, obesity, atherosclerosis, autoimmune diseases, and cancer has been ascertained, many aspects of the complex metabolic networks in which they are involved are still unknown and their therapeutic potential has not been exploited. Recent studies have reported high expression of Gal-3 and Gal-9 in different tumors such as glioma [67, 68] , ovariancarcinoma [69] , pancreatic ductal adenocarcinoma [70] , chronic lymphocytic leukemia [71] , and melanoma [23] . Given the implication of Gal-3 and Gal-9 in metabolic pathways that regulate tumor progression, we designed and synthesized a new molecule, SeDG-Bn, able to bind Gal-3 CRD and Gal-9N CRD and affect tumor cell mobility, invasion, and angiogenesis. We started from a selectively functionalized digalactosyl selenide; thus, selecting a bridging selenium atom between the sugars in order to exploit the interesting in vivo properties of selenium. For this purpose, we took advantage of two synthetic procedures developed in our laboratories to carry out both the selenoglycosidation step by using elementary selenium as a cheap selenenating agent [43] , and the regioselective installation of the benzyl group at the O-3 position of the sugar [45] . SeDG-Bn proved capable of binding Gal-3 CRD and Gal-9N CRD with low micromolar affinity, demonstrating also a higher affinity for both proteins than the SeDG analogue, which lacks the benzyl groups. Encouraged by its good affinity, we performed a biological screening on a selected panel of human cancer cell lines and healthy cells expressing Gal-3 CRD or Gal-9N CRD . After we verified the non-toxicity of SeDG-Bn, we demonstrated its selective antiproliferative action towards melanoma cells without showing significant cytotoxicity towards normal cells. Then, we selected the WM266 cell line to evaluate the effect of our compounds on galectins-mediated activities that lead to cancer progressions [48] . Through a wound healing assay, we showed that cells treated with the SeDG-Bn migrate more slowly than the control and, remarkably, the cell migration induced by Gal-3 CRD is specifically abolished by the presence of the examined compound. Moreover, SeDG-Bn inhibits both the melanoma cell invasion of the extracellular matrix and the formation of pseudo-capillaries from endothelial cells, which are two processes that are important for metastatic progression. Our multidisciplinary approach, encompassing molecular modeling and design, synthetic, and biological chemistry as well as biological characterization, emphasizes the role of synthetic digalactosides as galectin inhibitors in cancer treatment and identified a new molecule with a specific biological activity. Moreover, our observations highlight the great potential of galectins as therapeutic targets for the design of selective molecules with enhanced activity. Supplementary Materials: The following supporting information can be downloaded at: https: //www.mdpi.com/article/10.3390/ijms23052581/s1. Galectin-Carbohydrate Interactions in Biomedicine and Biotechnology Galectinomics: Finding themes in complexity Phylogenetic analysis of the vertebrate galectin family Galectins at a glance Design-Functionality relationships for adhesion/growth-regulatory galectins Extracellular and intracellular small-molecule galectin-3 inhibitors Galectin-3: One Molecule for an Alphabet of Diseases, from A to Z Galectins as modulators of cell adhesion Control of Angiogenesis by Galectins Involves the Release of Platelet-Derived Proangiogenic Factors Galectin-3 in autoimmunity and autoimmune diseases Galactin-1, 3 and 9: Potential biomarkers in idiopathic pulmonary fibrosis and other interstitial lung diseases Galectins as modulators of tumour progression A Multi-Targeting Approach to Fight SARS-CoV-2 More is always better than one: The N-terminal domain of the spike protein as another emerging target for hampering the SARS-CoV-2 attachment to host cells Microbe-host interactions are positively and negatively regulated by galectin-glycan interactions Intracellular functions of galectins Galectin-3, a novel binding partner of beta-catenin Galectin-3 in apoptosis, a novel therapeutic target The role of galectin-3 in promotion of the inflammatory response Galectins: Structure, function and therapeutic potential Lactosamine-Based Derivatives as Tools to Delineate the Biological Functions of Galectins: Application to Skin Tissue Repair Synthetic 1,2,3-triazole-linked glycoconjugates bind with high affinity to human galectin Galectin-3 captures interferon-gamma in the tumor matrix reducing chemokine gradient production and T-cell tumor infiltration Extracellular galectin-3 induces MMP9 expression by activating p38 MAPK pathway via lysosome-associated membrane protein-1 (LAMP1) Tailored Multivalent Neo-Glycoproteins: Synthesis, Evaluation, and Application of a Library of Galectin-3-Binding Glycan Ligands Synthesis of a phenyl thio-[small beta]-d-galactopyranoside library from 1,5-difluoro-2,4-dinitrobenzene: Discovery of efficient and selective monosaccharide inhibitors of galectin-Org Regioselective 3-O-Substitution of Unprotected Thiodigalactosides: Direct Route to Galectin Inhibitors C2-symmetrical thiodigalactoside bis-benzamido derivatives as high-affinity inhibitors of galectin-3: Efficient lectin inhibition through double arginine-arene interactions Galectin-3-Binding Glycomimetics that Strongly Reduce Bleomycin-Induced Lung Fibrosis and Modulate Intracellular Glycan Recognition 1H-1,2,3-triazol-1-yl thiodigalactoside derivatives as high affinity galectin-3 inhibitors A Selective Galactose-Coumarin-Derived Galectin-3 Inhibitor Demonstrates Involvement of Galectin-3-glycan Interactions in a Pulmonary Fibrosis Model Synthesis of diglycosylated (di)sulfides and comparative evaluation of their antiproliferative effect against tumor cell lines: A focus on the nature of sugar-recognizing mediators involved Sweet Selenium: Synthesis and Properties of Selenium-Containing Sugars and Derivatives Selenium compounds, apoptosis and other types of cell death: An overview for cancer therapy Selenium compounds as therapeutic agents in cancer Chemistry of biologically important synthetic organoselenium compounds Organoselenium and organotellurium compounds: Toxicology and pharmacology Structure of Mpro from SARS-CoV-2 and discovery of its inhibitors Glutathione peroxidase-like antioxidant activity of diaryl diselenides: A Mechanistic Study Diferrocenyl diselenides: Excellent thiol peroxidase-like antioxidants Selenogalactoside Compounds for the Prevention and Treatment of Diseases Associated with Galectin and the Use Thereof Thio-and selenoglycosides as ligands for biomedically relevant lectins: Valencyactivity correlations for benzene-based dithiogalactoside clusters and first assessment for (di)selenodigalactosides Switchable synthesis of glycosyl selenides or diselenides with direct use of selenium as selenating agent Practical β-stereoselective O-glycosylation of phenols with penta-O-acetyl-β-D-Gglucopyranose Polymethylhydrosiloxane (PMHS): A convenient option for synthetic applications of the iodine/silane combined reagent. Straightforward entries to 2-hydroxyglycals and useful building-blocks of glucuronic acid and glucosamine Efficient and direct synthesis of saccharidic 1,2-ethylidenes, orthoesters, and glycals from peracetylated sugars via the in situ generation of glycosyl iodides with I2/Et3SiH An Easy and Versatile Approach for the Regioselective De-O-benzylation of Protected Sugars Based on the I2/Et3SiH Combined System Thiodigalactoside inhibits murine cancers by concurrently blocking effects of galectin-1 on immune dysregulation, angiogenesis and protection against oxidative stress Chemically modified, non-anticoagulant heparin derivatives are potent galectin-3 binding inhibitors and inhibit circulating galectin-3-promoted metastasis Unraveling galectin-1 as a novel therapeutic target for cancer Galectin-Inhibitory Thiodigalactoside Ester Derivatives Have Antimigratory Effects in Cultured Lung and Prostate Cancer Cells Monocyte migration driven by galectin-3 occurs through distinct mechanisms involving selective interactions with the extracellular matrix Galectin-1 and -9 in angiogenesis: A sweet couple Galectin-3 acts as an angiogenic switch to induce tumor angiogenesis via Jagged-1/Notch activation Galectin-3 induces endothelial cell morphogenesis and angiogenesis Protein and ligand preparation: Parameters, protocols, and influence on virtual screening enrichments OPLS3: A Force Field Providing Broad Coverage of Drug-like Small Molecules and Proteins Extra precision glide: Docking and scoring incorporating a model of hydrophobic enclosure for protein-ligand complexes Meeting modern challenges in visualization and analysis Investigating the properties of TBA variants with twin thrombin binding domains Benzodifurans for biomedical applications: BZ4, a selective anti-proliferative and anti-amyloid lead compound Extracellular vesiclemediated delivery of miR-101 inhibits lung metastasis in osteosarcoma Therapeutic Potential of a Novel αvβ 3 Antagonist to Hamper the Aggressiveness of Mesenchymal Triple Negative Breast Cancer Sub-Type Chemical Modification for Proteolytic Stabilization of the Selective αvβ3 Integrin RGDechi Peptide: In Vitro and in Vivo Activities on Malignant Melanoma Cells Unveiling a VEGF-mimetic peptide sequence in IQGAP1 protein Macropinocytosis requires Gal-3 in a subset of patient-derived glioblastoma stem cells Galectin-9: A predictive biomarker negatively regulating immune response in glioma patients Tissue and plasma levels of galectins in patients with high grade serous ovarian carcinoma as new predictive biomarkers Detection of pancreatic ductal adenocarcinoma with galectin-9 serum levels Increased serum levels of Galectin-9 in patients with chronic lymphocytic leukemia The docking data presented in this study will be openly available in our public repository upon acceptance. The work has been performed under the Project HPC-EUROPA3 (INFRAIA-2016-1-730897), with the support of the EC Research Innovation Action under the H2020 Programme; in particular, we gratefully acknowledge the computer resources and technical support provided by ARCHER UK National Supercomputing Service. We also acknowledge CSCS (pr126) for computer time. We are thankful for the use of molecular graphics and analyses UCSF ChimeraX, developed by the Resource for Biocomputing, Visualization, and Informatics at the University of California, San Francisco, with support from National Institutes of Health R01-GM129325 and the Office of Cyber Infrastructure and Computational Biology, National Institute of Allergy and Infectious Diseases. We thank also Alfonso Esposito for his valuable technical assistance. The authors declare no competing financial interest. Digalactosyl selenide SeDG-Bn 3-3 -O-dibenzyl Digalactosyl selenide with IPTG Isopropyl β-D-1-thiogalactopyranoside