Design, Synthesis and Biological Evaluation of New Carbohydrate-Based Coumarin Derivatives as Selective Carbonic Anhydrase IX Inhibitors via “Click” Reaction

In this work, we designed a series of new carbohydrate-based coumarin carbonic anhydrase IX inhibitors by using 1,2,3-triazoles as linker. Next, these designed compounds were synthesized by the optimized one-pot click chemistry reaction condition. Subsequently, these target compounds were assayed for the inhibition of three carbonic anhydrase isoforms (CA I, CA II and CA IX). Intriguingly, all the compounds showed better CA IX inhibitory activity than initial coumarin fragments. Among them, compound 10a (IC50: 11 nM) possessed the most potent CA IX inhibitory activity, which was more potent than the reference drug acetazolamide (IC50: 30 nM). Notably, compound 10a showed 3018-fold, 1955-fold selectivity relative to CA I and CA II, respectively. Meanwhile, representative compounds could reduce tumor cell viability and the extracellular acidification in HT-29 and MDA-MB-231 cancer cell lines. Even more interestingly, our target compounds had no apparent cytotoxicity toward MCF-10A cell line. In addition, the in vitro stability assays also indicated our developed compounds possessed good liver microsomal metabolic stabilities and plasma stability. Furthermore, representative compounds revealed relatively low hERG cardiac toxicity and acute toxicity. Furthermore, docking studies were carried out to understand the interactions of our target compounds with the protein target CA IX. Collectively, our results suggest that compound 10a, as a selective CA IX inhibitor, could be an important lead compound for further optimization and development as an anticancer agent.


Introduction
Carbonic anhydrases (CAs, EC 4.2.1.1) are a very omnipresent zinc metalloenzymes, which are encoded by eight families, namely α-, β-, γ-, δ-, ζ-, η-, θ and ι. The α-CAs are represented by 15 different isoforms of human carbonic anhydrase (hCA) varies from CA I-CA XIV [1][2][3][4][5][6]. The expression of CAs is widely distributed in various tissues and they play an important role in numerous physiological and pathological processes, such as pH and CO 2 homeostasis, respiration and transport of CO 2 and HCO 3 − , biosynthetic reactions, tumorigenicity, bone resorption, calcification, etc., [7][8][9][10]. The transmembrane isoform CA IX has been characterized as biomarkers for several tumors. CA IX is not expressed in healthy tissues, but it is overexpressed in tumor microenvironment. CA IX is also one of the best markers for cellular hypoxia, and it promotes the acidification of the extracellular matrix, which is conducive to the growth and metastasis of tumor cells [11][12][13][14]. Therefore, CA IX is considered as a potential antitumor drug target. Among them, the natural product coumarin 1 acted as an effective CAI against CA I and CA II. In addition, it showed medium inhibitory activity against CA IX. The simple nonsubstituted derivative 2 is generally very poor CA IX inhibitor. As CAIs isosters of the simple coumarin 2, thiocoumarin, 2-thioxocoumarin, and dithiocoumarin were reported and these isosters showed improved inhibitory activity against CA IX. It is worth noting that these structures showed only micromolar level of inhibitory activity against CA IX. Interestingly, 7-hydroxy coumarin (umbelliferone) and 4-methyl-7-hydroxy coumarin (4methylumbelliferone) not only showed significantly selective for CA IX, but also possessed more potent inhibitory activity against CA IX [26][27][28]. Therefore, structural modification of these simple coumarins need to be performed to improve the inhibitory effect on hCA IX. For another example, glycosyl coumarins 8 and 9 showed low nanomolar inhibitory activity against CA IX [29].
It is widely known that coumarins as hydrophobic fragment could interact with the hydrophobic pocket at the entrance of the CA active site cavity. However, the neighboring hydrophilic subpocket cannot be fully occupied by these simple coumarins. Therefore, the CA IX inhibitory activity may be improved when some lipophilic bulky moieties was introduced into coumarins. The tail approaches have been extensively utilized in the design of CAIs, especially sugar tail approaches [30][31][32]. In our previous study, the introduction of the sugar moiety led to a significant enhancement of the CA inhibitory activity. As an extension of coumarins-based CAIs, we designed and synthesized a series of glycosyl conjugates of coumarins which were assessed for their inhibitory actions against three hCAs such as the hCA I, II, IX. Next, the inhibition of proliferation, the extracellular acidification Among them, the natural product coumarin 1 acted as an effective CAI against CA I and CA II. In addition, it showed medium inhibitory activity against CA IX. The simple nonsubstituted derivative 2 is generally very poor CA IX inhibitor. As CAIs isosters of the simple coumarin 2, thiocoumarin, 2-thioxocoumarin, and dithiocoumarin were reported and these isosters showed improved inhibitory activity against CA IX. It is worth noting that these structures showed only micromolar level of inhibitory activity against CA IX. Interestingly, 7-hydroxy coumarin (umbelliferone) and 4-methyl-7-hydroxy coumarin (4methylumbelliferone) not only showed significantly selective for CA IX, but also possessed more potent inhibitory activity against CA IX [26][27][28]. Therefore, structural modification of these simple coumarins need to be performed to improve the inhibitory effect on hCA IX. For another example, glycosyl coumarins 8 and 9 showed low nanomolar inhibitory activity against CA IX [29].
It is widely known that coumarins as hydrophobic fragment could interact with the hydrophobic pocket at the entrance of the CA active site cavity. However, the neighboring hydrophilic subpocket cannot be fully occupied by these simple coumarins. Therefore, the CA IX inhibitory activity may be improved when some lipophilic bulky moieties was introduced into coumarins. The tail approaches have been extensively utilized in the design of CAIs, especially sugar tail approaches [30][31][32]. In our previous study, the introduction of the sugar moiety led to a significant enhancement of the CA inhibitory activity. As an extension of coumarins-based CAIs, we designed and synthesized a series of glycosyl conjugates of coumarins which were assessed for their inhibitory actions against three hCAs such as the hCA I, II, IX. Next, the inhibition of proliferation, the extracellular acidification of tumor cells was explored. Meanwhile, the metabolic stability of target compounds was evaluated. Finally, the hERG cardiac toxicity and acute toxicity were also investigated. of tumor cells was explored. Meanwhile, the metabolic stability of target compounds was evaluated. Finally, the hERG cardiac toxicity and acute toxicity were also investigated.

Design of New Carbohydrate-Based Coumarin Derivatives
The target compounds comprised three main elements: (i) a sugar-tail scaffold, (ii) a coumarin moiety, and (iii) a 1,2,3-triazole linker. Due to possessing the CA IX selectivity, the coumarin moieties was chosen as initial fragments. In order to form more hydrogen bonds with the enzymatic cavity, different monosaccharide scaffolds were chosen as the hydrophilic groups for drawing favorable interaction with the hydrophilic halves of the CA active site. The 1,2,3-triazole ring is a bioisostere of amide group, endowed with a moderate dipole character, hydrogen bonding capability, rigidity, stability in the in vivo environment, and shows good tolerance to metabolic processes as well as to pH fluctuations [33][34][35][36][37]. Therefore, the 1,2,3-triazole played the role as a biocompatible covalent linker between the sugar and coumarin moiety. The approach is shown in Figure 2.

Chemistry
In order to prepare the designed compounds, we decided to develop an one-pot synthesis of 1,2,3-triazoles from benzoyl glycosyl bromides such as 13a, sodium azide, and coumarin-derived alkynes. Among the most common monosaccharides, D-glucose 11a was utilized as the starting material to synthesize model substrate 13a. As described in Scheme 1, the intermediate 12a was produced via benzoylation of 11a. The intermediate

Chemistry
In order to prepare the designed compounds, we decided to develop an one-pot synthesis of 1,2,3-triazoles from benzoyl glycosyl bromides such as 13a, sodium azide, and coumarin-derived alkynes. Among the most common monosaccharides, D-glucose 11a was utilized as the starting material to synthesize model substrate 13a. As  The model reaction, as shown in Table 1, was performed in the presence of different types of solvents, copper sulfate at different reaction temperature and time to find optimum conditions. Notably, this reaction did not provide the desired product 15a when t-BuOH, DMSO, DMF were chosen as solvent (Table 1, entries 1-3). When H2O was used as solvent, the desired product was observed (Table 1, entry 4). These preliminary results encouraged us to do further reaction conditions screening to increase the yield of 15a. Subsequently, increasing the reaction temperature and prolonging the reaction time could not improve the yield (Table 1, entry 5-9). It is worth noting that the hydrolysate of compound 13a was found due to the decomposition of benzoyl glycosyl bromide by water. Therefore, H2O is not suitable as a solvent alone. Next, the use of mixed solvents can significantly increase the yield such as H2O/t-BuOH, H2O/DMF and H2O/DMSO (Table 1, entry 10-12). Among these three mixed solvents, H2O/t-BuOH provides the highest yield (42%). However, at present, the obtained yields are not satisfactory to us due to the partial hydrolysis of benzoyl glycosyl bromide 13a. Interestingly, the yield can be further improved by adjusting the ratio of t-BuOH and H2O (Table 1, entry [13][14][15]. The best results were obtained when using H2O/t-BuOH (1:3). Several other catalysts, including CuBr and CuSO4·5H2O/VcNa were next evaluated and CuSO4·5H2O/VcNa was found to have higher catalytic activity (Table 1, entry [16][17]. Subsequently, the yield of the product 15a was increased to 90% when the reaction temperature was increased to 70 °C (Table 1, entry [18][19][20]. In addition, the reaction still has a high yield even if the reaction time is reduced to 30 min. Therefore, the optimal reaction condition entailed the use of CuSO4·5H2O/VcNa as a catalyst in H2O/t-BuOH (1:3) at 70 °C. The model reaction, as shown in Table 1, was performed in the presence of different types of solvents, copper sulfate at different reaction temperature and time to find optimum conditions. Notably, this reaction did not provide the desired product 15a when t-BuOH, DMSO, DMF were chosen as solvent (Table 1, entries 1-3). When H2O was used as solvent, the desired product was observed (Table 1, entry 4). These preliminary results encouraged us to do further reaction conditions screening to increase the yield of 15a. Subsequently, increasing the reaction temperature and prolonging the reaction time could not improve the yield (Table 1, entry 5-9). It is worth noting that the hydrolysate of compound 13a was found due to the decomposition of benzoyl glycosyl bromide by water. Therefore, H2O is not suitable as a solvent alone. Next, the use of mixed solvents can significantly increase the yield such as H2O/t-BuOH, H2O/DMF and H2O/DMSO (Table 1, entry 10-12). Among these three mixed solvents, H2O/t-BuOH provides the highest yield (42%). However, at present, the obtained yields are not satisfactory to us due to the partial hydrolysis of benzoyl glycosyl bromide 13a. Interestingly, the yield can be further improved by adjusting the ratio of t-BuOH and H2O (Table 1, entry [13][14][15]. The best results were obtained when using H2O/t-BuOH (1:3). Several other catalysts, including CuBr and CuSO4·5H2O/VcNa were next evaluated and CuSO4·5H2O/VcNa was found to have higher catalytic activity (Table 1, entry [16][17]. Subsequently, the yield of the product 15a was increased to 90% when the reaction temperature was increased to 70 °C (Table 1, entry [18][19][20]. In addition, the reaction still has a high yield even if the reaction time is reduced to 30 min. Therefore, the optimal reaction condition entailed the use of CuSO4·5H2O/VcNa as a catalyst in H2O/t-BuOH (1:3) at 70 °C. The model reaction, as shown in Table 1, was performed in the presence of different types of solvents, copper sulfate at different reaction temperature and time to find optimum conditions. Notably, this reaction did not provide the desired product 15a when t-BuOH, DMSO, DMF were chosen as solvent (  [16][17]. Subsequently, the yield of the product 15a was increased to 90% when the reaction temperature was increased to 70 • C (Table 1, entry [18][19][20]. In addition, the reaction still has a high yield even if the reaction time is reduced to 30 min. Therefore, the optimal reaction condition entailed the use of CuSO 4 ·5H 2 O/VcNa as a catalyst in To study the scope of the above reaction, a wide range of diversely benzoyl glycosyl bromides were tried. To our delight, under the conditions mentioned above, two coumarinderived alkynes could react smoothly with the in situ generated organic azides from benzoyl glycosyl bromides and NaN 3 to give the corresponding products in excellent yields ( Table 2). CuSO4·5H2O/VcNa were next evaluated and CuSO4·5H2O/VcNa was found to have higher catalytic activity (Table 1, entry [16][17]. Subsequently, the yield of the product 15a was increased to 90% when the reaction temperature was increased to 70 °C (Table 1, entry [18][19][20]. In addition, the reaction still has a high yield even if the reaction time is reduced to 30 min. Therefore, the optimal reaction condition entailed the use of CuSO4·5H2O/VcNa as a catalyst in H2O/t-BuOH (1:3) at 70 °C.  To study the scope of the above reaction, a wide range of diversely benzoyl glycosyl bromides were tried. To our delight, under the conditions mentioned above, two coumarin-derived alkynes could react smoothly with the in situ generated organic azides from benzoyl glycosyl bromides and NaN3 to give the corresponding products in excellent yields ( Table 2). To study the scope of the above reaction, a wide range of diversely benzoyl glycosyl bromides were tried. To our delight, under the conditions mentioned above, two coumarin-derived alkynes could react smoothly with the in situ generated organic azides from benzoyl glycosyl bromides and NaN3 to give the corresponding products in excellent yields ( Table 2). Table 2. Synthesis of triazoles 15a-15p using the optimal reaction condition a .

Entry Benzoyl Glycosyl Bromides
Products Yield b (%) To study the scope of the above reaction, a wide range of diversely benzoyl glycosyl bromides were tried. To our delight, under the conditions mentioned above, two coumarin-derived alkynes could react smoothly with the in situ generated organic azides from benzoyl glycosyl bromides and NaN3 to give the corresponding products in excellent yields ( Table 2). Table 2. Synthesis of triazoles 15a-15p using the optimal reaction condition a .

Entry Benzoyl Glycosyl Bromides
Products Yield b (%) To study the scope of the above reaction, a wide range of diversely benzoyl glycosyl bromides were tried. To our delight, under the conditions mentioned above, two coumarin-derived alkynes could react smoothly with the in situ generated organic azides from benzoyl glycosyl bromides and NaN3 to give the corresponding products in excellent yields ( Table 2). Table 2. Synthesis of triazoles 15a-15p using the optimal reaction condition a .

Entry Benzoyl Glycosyl Bromides
Products Yield b (%) To study the scope of the above reaction, a wide range of diversely benzoyl glycosyl bromides were tried. To our delight, under the conditions mentioned above, two coumarin-derived alkynes could react smoothly with the in situ generated organic azides from benzoyl glycosyl bromides and NaN3 to give the corresponding products in excellent yields ( Table 2). Table 2. Synthesis of triazoles 15a-15p using the optimal reaction condition a .

Entry Benzoyl Glycosyl Bromides
Products Yield b (%) To study the scope of the above reaction, a wide range of diversely benzoyl glycosyl bromides were tried. To our delight, under the conditions mentioned above, two coumarin-derived alkynes could react smoothly with the in situ generated organic azides from benzoyl glycosyl bromides and NaN3 to give the corresponding products in excellent yields ( Table 2). To study the scope of the above reaction, a wide range of diversely benzoyl glycosyl bromides were tried. To our delight, under the conditions mentioned above, two coumarin-derived alkynes could react smoothly with the in situ generated organic azides from benzoyl glycosyl bromides and NaN3 to give the corresponding products in excellent yields ( Table 2). To study the scope of the above reaction, a wide range of diversely benzoyl glycosyl bromides were tried. To our delight, under the conditions mentioned above, two coumarin-derived alkynes could react smoothly with the in situ generated organic azides from benzoyl glycosyl bromides and NaN3 to give the corresponding products in excellent yields ( Table 2). To study the scope of the above reaction, a wide range of diversely benzoyl glycosyl bromides were tried. To our delight, under the conditions mentioned above, two coumarin-derived alkynes could react smoothly with the in situ generated organic azides from benzoyl glycosyl bromides and NaN3 to give the corresponding products in excellent yields ( Table 2). To study the scope of the above reaction, a wide range of diversely benzoyl glycosyl bromides were tried. To our delight, under the conditions mentioned above, two coumarin-derived alkynes could react smoothly with the in situ generated organic azides from benzoyl glycosyl bromides and NaN3 to give the corresponding products in excellent yields ( Table 2).  The target compound 10a was obtained by deprotection of benzoyl groups using Na-OMe/MeOH in high yield (Scheme 3) [38]. All other target compounds 10b-10p were also prepared by this procedure. Since there are two possible configurations (α and β) of the C1-hydrogen, the configuration of representative compound 10a was studied by the data of 1 H-NMR(see Supplementary Materials) in this paper. The data of 1 H-NMR showed that the chemical shift of C1-hydrogen (double peak) appears in 5.57 ppm, the coupling constant is 9.3 Hz. These results are consistent with the β configuration. The target compound 10a was obtained by deprotection of benzoyl groups using Na-OMe/MeOH in high yield (Scheme 3) [38]. All other target compounds 10b-10p were also prepared by this procedure. Since there are two possible configurations (α and β) of the C1-hydrogen, the configuration of representative compound 10a was studied by the data of 1 H-NMR(see Supplementary Materials) in this paper. The data of 1 H-NMR showed that the chemical shift of C1-hydrogen (double peak) appears in 5.57 ppm, the coupling constant is 9.3 Hz. These results are consistent with the β configuration.  The target compound 10a was obtained by deprotection of benzoyl groups using Na-OMe/MeOH in high yield (Scheme 3) [38]. All other target compounds 10b-10p were also prepared by this procedure. Since there are two possible configurations (α and β) of the C1-hydrogen, the configuration of representative compound 10a was studied by the data of 1 H-NMR(see Supplementary Materials) in this paper. The data of 1 H-NMR showed that the chemical shift of C1-hydrogen (double peak) appears in 5.57 ppm, the coupling constant is 9.3 Hz. These results are consistent with the β configuration. The target compound 10a was obtained by deprotection of benzoyl groups using Na-OMe/MeOH in high yield (Scheme 3) [38]. All other target compounds 10b-10p were also prepared by this procedure. Since there are two possible configurations (α and β) of the C1-hydrogen, the configuration of representative compound 10a was studied by the data of 1 H-NMR(see Supplementary Materials) in this paper. The data of 1 H-NMR showed that the chemical shift of C1-hydrogen (double peak) appears in 5.57 ppm, the coupling constant is 9.3 Hz. These results are consistent with the β configuration.  The target compound 10a was obtained by deprotection of benzoyl groups using Na-OMe/MeOH in high yield (Scheme 3) [38]. All other target compounds 10b-10p were also prepared by this procedure. Since there are two possible configurations (α and β) of the C1-hydrogen, the configuration of representative compound 10a was studied by the data of 1 H-NMR(see Supplementary Materials) in this paper. The data of 1 H-NMR showed that the chemical shift of C1-hydrogen (double peak) appears in 5.57 ppm, the coupling constant is 9.3 Hz. These results are consistent with the β configuration. The target compound 10a was obtained by deprotection of benzoyl groups using Na-OMe/MeOH in high yield (Scheme 3) [38]. All other target compounds 10b-10p were also prepared by this procedure. Since there are two possible configurations (α and β) of the C1-hydrogen, the configuration of representative compound 10a was studied by the data of 1 H-NMR(see Supplementary Materials) in this paper. The data of 1 H-NMR showed that the chemical shift of C1-hydrogen (double peak) appears in 5.57 ppm, the coupling constant is 9.3 Hz. These results are consistent with the β configuration.  The target compound 10a was obtained by deprotection of benzoyl groups using Na-OMe/MeOH in high yield (Scheme 3) [38]. All other target compounds 10b-10p were also prepared by this procedure. Since there are two possible configurations (α and β) of the C1-hydrogen, the configuration of representative compound 10a was studied by the data of 1 H-NMR(see Supplementary Materials) in this paper. The data of 1 H-NMR showed that the chemical shift of C1-hydrogen (double peak) appears in 5.57 ppm, the coupling constant is 9.3 Hz. These results are consistent with the β configuration. The target compound 10a was obtained by deprotection of benzoyl groups using Na-OMe/MeOH in high yield (Scheme 3) [38]. All other target compounds 10b-10p were also prepared by this procedure. Since there are two possible configurations (α and β) of the C1-hydrogen, the configuration of representative compound 10a was studied by the data of 1 H-NMR(see Supplementary Materials) in this paper. The data of 1 H-NMR showed that the chemical shift of C1-hydrogen (double peak) appears in 5.57 ppm, the coupling constant is 9.3 Hz. These results are consistent with the β configuration.  The target compound 10a was obtained by deprotection of benzoyl groups using Na-OMe/MeOH in high yield (Scheme 3) [38]. All other target compounds 10b-10p were also prepared by this procedure. Since there are two possible configurations (α and β) of the C1-hydrogen, the configuration of representative compound 10a was studied by the data of 1 H-NMR(see Supplementary Materials) in this paper. The data of 1 H-NMR showed that the chemical shift of C1-hydrogen (double peak) appears in 5.57 ppm, the coupling constant is 9.3 Hz. These results are consistent with the β configuration. The target compound 10a was obtained by deprotection of benzoyl groups using Na-OMe/MeOH in high yield (Scheme 3) [38]. All other target compounds 10b-10p were also prepared by this procedure. Since there are two possible configurations (α and β) of the C1-hydrogen, the configuration of representative compound 10a was studied by the data of 1 H-NMR(see Supplementary Materials) in this paper. The data of 1 H-NMR showed that the chemical shift of C1-hydrogen (double peak) appears in 5.57 ppm, the coupling constant is 9.3 Hz. These results are consistent with the β configuration.  The target compound 10a was obtained by deprotection of benzoyl groups using Na-OMe/MeOH in high yield (Scheme 3) [38]. All other target compounds 10b-10p were also prepared by this procedure. Since there are two possible configurations (α and β) of the C1-hydrogen, the configuration of representative compound 10a was studied by the data of 1 H-NMR(see Supplementary Materials) in this paper. The data of 1 H-NMR showed that the chemical shift of C1-hydrogen (double peak) appears in 5.57 ppm, the coupling constant is 9.3 Hz. These results are consistent with the β configuration. The target compound 10a was obtained by deprotection of benzoyl groups using Na-OMe/MeOH in high yield (Scheme 3) [38]. All other target compounds 10b-10p were also prepared by this procedure. Since there are two possible configurations (α and β) of the C1-hydrogen, the configuration of representative compound 10a was studied by the data of 1 H-NMR(see Supplementary Materials) in this paper. The data of 1 H-NMR showed that the chemical shift of C1-hydrogen (double peak) appears in 5.57 ppm, the coupling constant is 9.3 Hz. These results are consistent with the β configuration.  The target compound 10a was obtained by deprotection of benzoyl groups using Na-OMe/MeOH in high yield (Scheme 3) [38]. All other target compounds 10b-10p were also prepared by this procedure. Since there are two possible configurations (α and β) of the C1-hydrogen, the configuration of representative compound 10a was studied by the data of 1 H-NMR(see Supplementary Materials) in this paper. The data of 1 H-NMR showed that the chemical shift of C1-hydrogen (double peak) appears in 5.57 ppm, the coupling constant is 9.3 Hz. These results are consistent with the β configuration. The target compound 10a was obtained by deprotection of benzoyl groups using Na-OMe/MeOH in high yield (Scheme 3) [38]. All other target compounds 10b-10p were also prepared by this procedure. Since there are two possible configurations (α and β) of the C1-hydrogen, the configuration of representative compound 10a was studied by the data of 1 H-NMR(see Supplementary Materials) in this paper. The data of 1 H-NMR showed that the chemical shift of C1-hydrogen (double peak) appears in 5.57 ppm, the coupling constant is 9.3 Hz. These results are consistent with the β configuration.  The target compound 10a was obtained by deprotection of benzoyl groups using Na-OMe/MeOH in high yield (Scheme 3) [38]. All other target compounds 10b-10p were also prepared by this procedure. Since there are two possible configurations (α and β) of the C1-hydrogen, the configuration of representative compound 10a was studied by the data of 1 H-NMR(see Supplementary Materials) in this paper. The data of 1 H-NMR showed that the chemical shift of C1-hydrogen (double peak) appears in 5.57 ppm, the coupling constant is 9.3 Hz. These results are consistent with the β configuration. The target compound 10a was obtained by deprotection of benzoyl groups using Na-OMe/MeOH in high yield (Scheme 3) [38]. All other target compounds 10b-10p were also prepared by this procedure. Since there are two possible configurations (α and β) of the C1-hydrogen, the configuration of representative compound 10a was studied by the data of 1 H-NMR(see Supplementary Materials) in this paper. The data of 1 H-NMR showed that the chemical shift of C1-hydrogen (double peak) appears in 5.57 ppm, the coupling constant is 9.3 Hz. These results are consistent with the β configuration.  The target compound 10a was obtained by deprotection of benzoyl groups using Na-OMe/MeOH in high yield (Scheme 3) [38]. All other target compounds 10b-10p were also prepared by this procedure. Since there are two possible configurations (α and β) of the C1-hydrogen, the configuration of representative compound 10a was studied by the data of 1 H-NMR(see Supplementary Materials) in this paper. The data of 1 H-NMR showed that the chemical shift of C1-hydrogen (double peak) appears in 5.57 ppm, the coupling constant is 9.3 Hz. These results are consistent with the β configuration. The target compound 10a was obtained by deprotection of benzoyl groups using Na-OMe/MeOH in high yield (Scheme 3) [38]. All other target compounds 10b-10p were also prepared by this procedure. Since there are two possible configurations (α and β) of the C1-hydrogen, the configuration of representative compound 10a was studied by the data of 1 H-NMR(see Supplementary Materials) in this paper. The data of 1 H-NMR showed that the chemical shift of C1-hydrogen (double peak) appears in 5.57 ppm, the coupling constant is 9.3 Hz. These results are consistent with the β configuration.  The target compound 10a was obtained by deprotection of benzoyl groups using NaOMe/MeOH in high yield (Scheme 3) [38]. All other target compounds 10b-10p were also prepared by this procedure. Since there are two possible configurations (α and β) of the C1-hydrogen, the configuration of representative compound 10a was studied by the data of 1 H-NMR(see Supplementary Materials) in this paper. The data of 1 H-NMR showed that the chemical shift of C1-hydrogen (double peak) appears in 5.57 ppm, the coupling constant is 9.3 Hz. These results are consistent with the β configuration.
The target compound 10a was obtained by deprotection of benzoyl groups using Na-OMe/MeOH in high yield (Scheme 3) [38]. All other target compounds 10b-10p were also prepared by this procedure. Since there are two possible configurations (α and β) of the C1-hydrogen, the configuration of representative compound 10a was studied by the data of 1 H-NMR(see Supplementary Materials) in this paper. The data of 1 H-NMR showed that the chemical shift of C1-hydrogen (double peak) appears in 5.57 ppm, the coupling constant is 9.3 Hz. These results are consistent with the β configuration.

Carbonic Anhydrase Inhibition
The CA inhibitory activities of compounds 10a-10p were measured against three isoforms hCA I, II, and IX by the esterase assay. Meanwhile, acetazolamide (AZA), umbelliferone and 4-methylumbelliferone served as standard inhibitors were tested. The CA Scheme 3. Preparation of representative target compound 10a. Reagents and conditions: (a) NaOMe, MeOH, r.t., 99%.

Carbonic Anhydrase Inhibition
The CA inhibitory activities of compounds 10a-10p were measured against three isoforms hCA I, II, and IX by the esterase assay. Meanwhile, acetazolamide (AZA), umbelliferone and 4-methylumbelliferone served as standard inhibitors were tested. The CA inhibition data are shown in Table 3. The selectivity ratios for CA IX inhibition over CA I and CA II are also presented in Table 3. (i) All new compounds exhibited weak inhibitory activity against the cytosolic isoform hCA I, with IC 50 in the micromolar range (from 17.5 to 65.6 µM). In addition, this isoform was not at all inhibited by the parent coumarins umbelliferone and 4-methylumbelliferone.
(ii) The physiologically dominant isoform hCA II was similarly poorly inhibited by all the compounds, which showed IC 50 values in the range of 12.2-32.8 µM. The unsubstituted coumarin showed no inhibitory power against hCA II again. As CA I and CA II are considered important off-target CA isozyme, the poor inhibition against CA I and CA II is an interesting feature for compounds designed to target the tumor-associated enzymes.
(iii) The tumor-associated isoform hCA IX was weakly inhibited by parent coumarins umbelliferone and 4-methylumbelliferone. Notably, the introduction of sugar moieties into the coumarin ring at position 7 lead to a significant increase of the hCA IX inhibitory potency with IC 50 in the range of 0.011-0.132 µM. Among them, the glucose derivative 10a and mannose derivative 10c showed most potent hCA IX inhibitory activities with IC 50 values of 11 nM in the case of 10a and 15 nM in the case of 10c. The increase in inhibitory potency against this isoform may be due to the matching the hydrophilic half of the active site with the hydrophilic glucosyl moiety, and in addition to interacting with hydrophobic half of the active site with coumarins.
(iv) Next, the isoform selectivities of compounds 10a-10p against CAI, CAII and CAIX were investigated. The results indicated that the target compounds showed good CA IX isoform selectivity (up to 3018-fold, 1955-fold selectivity toward the hCA I and CA II, respectively). The selective inhibition of carbonic anhydrase IX of coumarin-derivatives might decrease the side effects, and these target compounds were expected to be interesting lead compounds for treatment of cancer.

Effect of the Synthetic Compounds on the Cell Viability of Human Cancer Cell Lines
In the present study, the potential cytotoxicity of compounds 10a-10p was evaluated in vitro by MTT assay. As a matter of fact, HT-29 cells possessed a high CA IX expression under ambient air, while MDA-MB-231 cells showed enhanced CA IX expression only upon hypoxia. So, half of each cell line were subjected to hypoxic conditions (0.5% O 2 , 5% CO 2 and 94.5% N 2 ) for inducing the expression of CA IX and the other half were incubated under normal conditions (normoxia). Subsequently, all synthesized derivatives were screened at a concentration of 100 µM to test their cytotoxicity in mentioned above cells.
Results in Table 4 indicated all the compounds were capable of reducing tumor cell viability in these two cancer cell lines under both of the conditions. At normal oxygen content, the inhibitory activity of compounds against HT-29 cancer cells was slightly higher than that in MDA-MB-231 cells. Affected by hypoxic conditions which leading to overexpression of CA IX, the inhibitory activity of compounds against the same tumor cells was higher than that under normal oxygen. The inhibitory activity of compounds against HT-29 under hypoxia conditions was similar to that of MDA-MB-231 cells. What's more, note that new carbohydrate-based coumarins designed in this study exhibited comparatively potent inhibitions at the enzyme level. Whereas weak cellular antitumor activities (high µM) of these compounds reported in this study were also observed, which revealed that hCA IX inhibitors are not potent enough in anti-tumor monotherapy only depending on the regulation of tumor microenvironment. Therefore, more pharmacological effects of CA IX inhibitors as antitumor agents need to be further explored. Since CA IX can regulate the extracellular pH of tumor cells, we next examined whether our inhibitors had any effect on the extracellular pH of tumor cells.

Extracellular pH Measurement of Human Cancer Cell Lines in the Presence of Compounds 10a and 10h
One important feature of CA IX is regulating the pH of tumor cells. The expression of CA IX in solid tumors will make the extracellular pH become more acidic to promote tumor growth and metastasis. Therefore, the effect on extracellular pH of the representative compounds 10a and 10h was tested in hypoxia and normoxic conditions. The results showed that most potent compound 10a had an important impact on the extracellular acidification of cancer cells in two different concentrations (0.1 and 0.5 mM) (Figure 3). Especially, the extracellular acidification of these two cancer cells could be slightly reversed by this compound under normoxic condition. In addition, the reversal of extracellular acidification is more significant under hypoxia condition. Furthermore, a dose-dependent effect was observed. In addition, compared with the most potent compound 10a, the compound 10h with weak enzymatic activity showed weaker reversal effect of extracellular acidification on both cell lines in hypoxia and normoxic conditions. One important feature of CA IX is regulating the pH of tumor cells. The expression of CA IX in solid tumors will make the extracellular pH become more acidic to promote tumor growth and metastasis. Therefore, the effect on extracellular pH of the representative compounds 10a and 10h was tested in hypoxia and normoxic conditions. The results showed that most potent compound 10a had an important impact on the extracellular acidification of cancer cells in two different concentrations (0.1 and 0.5 mM) (Figure 3). Especially, the extracellular acidification of these two cancer cells could be slightly reversed by this compound under normoxic condition. In addition, the reversal of extracellular acidification is more significant under hypoxia condition. Furthermore, a dose-dependent effect was observed. In addition, compared with the most potent compound 10a, the compound 10h with weak enzymatic activity showed weaker reversal effect of extracellular acidification on both cell lines in hypoxia and normoxic conditions.

In Vitro Cytotoxicity Studies on Normal Mammalian Cells
To obtain the preliminary safety profiles of these carbohydrate-based coumarins in normal human cells, compounds 10a and 10c were chosen to evaluated for cytotoxicity against human normal cell line MCF-10A by measuring its IC 50 in MTT assay. The data (Table 5.) suggested the carbohydrate-based coumarins 10a and 10c had no apparent cytotoxicity toward MCF-10A cell line (IC 50 >100 µM). A perfect drug should be stable enough to reach its pharmacological target. Therefore, target compounds were subjected to evaluate the in vitro metabolic stability. As presented in Table 6, all target compounds exhibited excellent metabolic stability with half-life values range from 1153 min to 1843 min. Therefore, target compounds are worthy of further evaluation. In fact, most drugs reach their targets to show their pharmacological activity via the circulation system, therefore, the plasma stability of the target compound is a crucial consideration for the evaluation of its druggability. Therefore, the in vitro plasma stability assay of the representative compound was also performed in this study. The results revealed that compound 10a possessed high plasma stability, with a half-life value greater than 1230 min (Table 7).

Inhibition Evaluation on hERG Channel
The hERG channel is a potassium ion channel that plays an important role in the regulation of cardiac repolarization. Blockage of the hERG channel can cause serious cardiac side effects, such as prolongation of the QT interval. For this reason, the most potent compounds 10a and 10c was further assessed for inhibition of the hERG channel. As shown in Table 8, compounds 10a and 10c showed weak inhibition (IC 50 > 50 µM). The results indicated that compounds 10a and 10c were less likely to produce cardiotoxicity.

Preliminary Assessment of the Acute Toxicity
Finally, acute toxicity of representative compounds was conducted in mice to preliminarily evaluate the risk of target compounds to mammals. The compounds 10a and 10c were given by intragastric administration from 500 mg/kg to 1500 mg/kg, respectively. Animals were observed for lethality for 14 days. The obtained preliminary data (Table 9) allowed us to conclude that LD 50 values for these compounds were expected to exceed 1500 mg/kg, because after administering this dose none of the animals died. In order to predict the interact modes between target compounds and hCA IX, molecular docking studies were employed in this study. At first, a redocking experiment for AAZ with CA IX (PDB:3IAI) was processed by using AutoDock 4.2.6 (Developed and maintained by the Olson Laboratory at the Scripps Research Institute, La Jolla, CA, USA) with the AutoDock4Zn force field. The redocking results showed that the crystal pose of AAZ was reproduced fairly well. Next, the same models were applied to the previous reported sugar-containing coumarins 8 and 9. As seen from Figure 4, the hydrophilic and hydrophobic region of active site were displayed with blue and red background, respectively. The obtained results indicated that the coumarin groups of the compounds 8 and 9 could bind to the hydrophobic part of the active site while the sugar fragments failed to penetrate deep into the hydrophilic part to form sufficient binding because they were far away from the hydrophilic part. The most potent compounds 10a and 10c occupied the hydrophilic and hydrophobic region of active site simultaneously. More specifically, the sugar groups and 1,2,3-trizaole group of the compounds could bind to the hydrophilic part while the coumarin groups could bind to the hydrophobic part.

Chemistry
Unless otherwise noted, all solvents and reagents were obtained from commercial suppliers and used without further purification. 1 H-NMR and 13 C-NMR spectra were measured with an AV-600 spectrometer (Bruker Bioscience, Billerica, MA, USA), with tet-

Chemistry
Unless otherwise noted, all solvents and reagents were obtained from commercial suppliers and used without further purification. 1 H-NMR and 13 C-NMR spectra were measured with an AV-600 spectrometer (Bruker Bioscience, Billerica, MA, USA), with tetramethyl silane as an internal standard. The 200-300 mesh silica gel (Qingdao Haiyang Chemical, Qingdao, China) was used for Column chromatography. TLC analysis was performed on silica gel plates GF254 (Qingdao Haiyang Chemical, Qingdao, China). ESI-MS were obtained on Agilent ESI-QTOF instrument (Agilent, Santa Clara, CA, USA). High resolution mass spectra (HRMS) were recorded on an Agilent Accurate-Mass Q-TOF 6530 (Agilent, Santa Clara, CA, USA) in ESI mode.

General Procedure for the Synthesis of Target Compound 10a
To a solution of 11a (5 g, 27.8 mmol) in pyridine (50 mL) was added benzoyl chloride (16.1 mL, 139 mmol) under the ice bath and stirred for 5 h at room temperature. The H 2 O (100 mL) and dichloromethane (DCM, 100 mL) were added and the mixture was stirred for 20 min. Next, the mixture was extracted by CH 2 Cl 2 and the organic layer was washed with 0.1 M HCl (20 mL), saturated NaHCO 3 (30 mL). In addition, the organic layer was washed with water, then dried over anhydrous Na 2 SO 4 and concentrated under reduced pressure to give 12a. To a solution of 12a in dichloromethane (50 mL) was added HBr-HOAc (35 mL). The reaction mixture was stirred for 5 h. Next, H 2 O (50 mL) was added and the reaction mixture was vigorously stirred for 20 min. The organic layer was washed with saturated NaHCO 3 (30 mL), and then dried over anhydrous Na 2 SO 4 and concentrated under reduced pressure. The crude products were purified by silica gel column chromatography eluting with 20% v/v EtOAc/n-hexane to give 13a as a white solid (16.8 g, 92%). To a solution of 13a (0.66 g, 1 mmol) in H 2 O/t-BuOH (1:3) (20 mL) was added NaN 3 (0.13 g, 2 mmol), 7-(Prop-2-ynyloxy)-2H-chromen-2-one (0.2 g, 1 mmol), CuSO 4 ·5H 2 O (25.0 mg, 0.1 mmol) and sodium ascorbate (59.4 mg, 0.3 mmol). The reaction mixture was vigorously stirred at 70 • C for 30 min. The solvent t-BuOH was removed under vacuo and the rest of the reaction mixture was extracted by using ethyl acetate. Furthermore, the organic layer was washed with water, then dried over anhydrous Na 2 SO 4 and concentrated under reduced pressure. The crude products were purified by silica gel column chromatography eluting with 25% v/v acetone/n-hexane to obtain 15a (0.76 g, 0.93 mmol). Next, the benzoyl group was removed by sodium methoxide to acquire a crude residue which was purified by column chromatography eluting with 10% v/v MeOH/CH 2 Cl 2 to obtain target compound 10a (0.37 g, 99%) as a white solid. The compounds 10b-10p were synthesized by the above synthetic method.  13 13  The esterase assay was used for testing the inhibitory activities of target compounds against carbonic anhydrases [39]. Firstly, the substrate (4-nitrophenyl acetate) was dissolved in DMSO and the concentration of solution of was 2 mM at 25 • C. Next, the inhibitor was also dissolved in DMSO and the concentration of solution of was 200 mM at 25 • C. Next, the stock solution of inhibitor was diluted to the specified concentration by using analysis buffer. In this case, 10 different inhibitor concentrations have been prepared, varying between 5 mM and 0.5 nM. In addition, the purchased CA was diluted with analysis buffer and added to 96-well plates. Next, different concentrations of inhibitors were also added to 96-well plates. To form the inhibitors and enzymes complex, E-I solutions were preincubated together for 4 h before the test at room temperature. Meanwhile, the nonenzymatic hydrolysis rates were tested by the blank control experiments. The absorbances were determined by using multifunctional enzyme marking instrument.

Cell Viability Assay
The two human cancer cell lines were cultivated for 24 h at 37 • C in a 96-well plate. Subsequently, different concentrations of inhibitors were added to 96-well plate and then was incubated for 96 hours. After the MTT solution was added to per well, cancer cells were further incubated for 4 h at 37 • C. Finally, formazan crystals were dissolved by using DMSO. The cancer cells were from American Type Culture Collection.

Metabolic Stabilities Assay
The metabolic stability of target compounds was tested by predicting intrinsic clearance rate of mouse liver microsomes. The liver microsomes were extracted from mouse. The test compounds and control working solutions were prepared. Meanwhile, the NADPH cofactors were prepared by using NADPH powder. Next, the empty plates were pre-warmed for 10 min minutes and the liver microsomes were diluted with 100 mM phosphate buffer to 0.56 mg/mL. Next, transfer microsome working solutions into pre-warmed incubation plates. In addition, the incubation plates were pre-incubated for 10 min at 40 • C with constant shaking. Then the blank plates were added liver microsomes, NAPDH cofactor and the stop solutions to quench the blank plates. Compound working solution containing microsomes was added to the incubation plates. Furthermore, the NCF60 plates were added buffer and mixed thoroughly 3 times. Then shake plates and incubate for 60 min at 40 • C. The quenching solutions were added to quenching plate T0 and incubation plate (T60) was added NAPDH cofactor, mixing 3 times thoroughly, and immediately remove mixture to quenching plate. Then shake plates and incubate for 60 min at 40 • C. The final each component in incubation medium include liver microsome protein, 1 µM test compound, 1 µM control compound, 0.99% MeOH and 0.01% DMSO. At 5, 10, 20, 30, and 60 min, quenching solution was added to plates. All sampling plates were shaken and centrifuged. The supernatant was transferred into HPLC water. Each plate was sealed and shaken for 10 min prior to LC-MS/MS analysis.

Plasma Stability Assay
Rat plasma was collected from adult male SD rats. The in vitro stability of 10a was determined in rat plasma and in PBS at 37 • C, at the meantime, samples were taken at time points 0, 10, 30, 60, and 120, and tested by HPLC. The rest of the time was incubated in a constant temperature shaker to compare different time points. Compare the peak areas of 10a in PBS or in rat plasma (pH = 7.4, 37 • C) before and after constant temperature incubation at different time points. The compounds in rat plasma were extracted with methanol, centrifuged, and repeated three times, and then detected with a preparative liquid phase at 254 wavelength.

Inhibition Evaluation on hERG K + Channel
HEK 293 cells were stably transfected with human Ether-a-gogo related gene (hERG) channel. The voltage-gated hERG potassium channel current was recorded at room temperature from randomly selected transfected cells using whole-cell recording technique with an EPC10 USB (HEKA) or Multiclamp 700 B amplifier (Molecular Devices) while electrical data was digitalized by Digidata1440 A with acquisition rate of 10 kHz and signals filtered at 2.5 KHz using Patchmaster or pClamp 10, respectively. Dofetilide was included as a positive control to ensure the accuracy and sensitivity of the test system. All experiments were performed in 3 times.

Acute Toxicity Assay
Groups of mice (5 per group) were treated intraperitoneally with target compounds. After the compound was administered, the latency of death was observed and recorded over a span of 14 days. The experimental protocols were evaluated and approved by the ethics committee of the Shenyang Pharmaceutical University.

Experimental Protocol for Docking Study
The molecular docking studies were performed using crystallographic structure of CA IX complex (PDB: 3IAI). The crystallographic protein was obtained by removing all bound water and ligands. In this study, the molecular docking was performed by inserting target compounds into the active pocket of CA IX. The interaction models between protein with target compounds were analyzed after the end of molecular docking.

Conclusions
In this article, different sugars were connected to coumarins using a 1,2,3-triazole linker for the first time. Thus, 16 new carbohydrate-based coumarin derivatives were synthesized by our optimized one-pot click chemistry reaction. The in vitro CA inhibitory activities of the new derivatives had been tested and the result showed that almost all carbohydrate-based coumarin derivatives showed potent inhibitory activity against hCA IX, with the values of IC 50 ranging from 11 to 132 nM. Intriguingly, all the compounds showed better CA IX inhibitory activity than initial coumarin fragments. Among all newly synthesized compounds, compound 10a was found to be the most potent and selective CA IX inhibitor (selectivity index in favor of CA I and CA II equal to 3018, 1955, respectively). Meanwhile, they efficiently reduce tumor cell viability and the extracellular acidification in two different cell lines (HT-29 and MDA-MB-231). In addition, representative compounds were not toxic to normal cell line. In particularly, the target compounds showed good metabolic stability and plasma stability in vitro. Moreover, representative compounds showed low hERG inhibition and low acute toxicity. Furthermore, docking studies were carried out to understand the interactions of our target compounds with the protein target CA IX. The obtained results allowed for the consideration of this compound 10a as an interesting lead for the development of a new class of selective CA IX inhibitors for the treatment of cancer.