Raltegravir was the first HIV-1 integrase inhibitor that gained FDA approval for use in the treatment of HIV-1 infection. Because of the emergence of IN inhibitor-resistant viral strains, there is a need to identify innovative second-generation IN inhibitors. Previously, we identified 2-thioxo-4-thiazolidinone (rhodanine)-containing compounds as IN inhibitors. Herein, we report the design, synthesis and docking studies of a series of novel rhodanine derivatives as IN inhibitors. All these compounds were further tested against human apurinic/apyrimidinic endonuclease 1 (APE1) to determine their selectivity. Two compounds showed significant cytotoxicity in a panel of human cancer cell lines. Taken together, our results show that rhodanines are a promising class of compounds for developing drugs with antiviral and anticancer properties.
Raltegravir was the first HIV-1 integrase inhibitor that gained FDA approval for use in the treatment ofHIV-1 infection. Because of the emergence of IN inhibitor-resistant viral strains, there is a need to identify innovative second-generation IN inhibitors. Previously, we identified 2-thioxo-4-thiazolidinone (rhodanine)-containing compounds as IN inhibitors. Herein, we report the design, synthesis and docking studies of a series of novel rhodanine derivatives as IN inhibitors. All these compounds were further tested against humanapurinic/apyrimidinic endonuclease 1 (APE1) to determine their selectivity. Two compounds showed significant cytotoxicity in a panel ofhumancancercell lines. Taken together, our results show that rhodanines are a promising class ofcompounds for developing drugs with antiviral and anticancer properties.
Raltegravir (MK-0518) is the first HIV-1 integrase (IN) inhibitor approved for the treatment ofHIV-1 infection in treatment-experienced patients, and heralds a new era in HIV/AIDS treatment [1,2]. IN inhibitors differ from other classes of antiretroviral agents in that they target the integration stage in the viral life cycle. The integration step, catalyzed by IN, results in the insertion of proviral DNA into the host genome. This occurs as a two-stage reaction. First, the proviral DNA is primed for integration by an endonucleolyticcleavage at its 3’ end. This step is called 3’-processing and takes place in the cytoplasm of an infected cell. IN remains bound to the processed viral DNA and translocates into nucleus aided by several interacting cellular cofactors. Inside the nucleus, IN catalyzes a staggered cleavage in the host DNA. The nucleophilic hydroxyl group at the 3’ end of the viral DNA is then inserted into the host DNA by a trans-esterification reaction. This process is called the strand transfer reaction. These two steps catalyzed by IN are critical for HIV-1 replication and infection [3,4,5].Nearly two decades of intensive efforts towards developing clinically useful HIV-1 IN inhibitors has resulted in a rich database of structural information [6,7,8]. This is particularly useful in designing second generation inhibitors in order to overcome the challenges faced with emerging resistant strains against current IN inhibitors, including raltegravir. Among the various structural classes, β-diketoacid based pharmacophores have shown high potency and selectivity towards inhibition of the IN-catalyzed strand transfer [9,10]. S-1360, a β-diketoacid bioisostere, was the one offirst IN targeted antiretroviral agent to undergo humanclinical trials. Although it failed to show clinical efficacy, S-1360 has served as a good template for the design of several second-generation IN inhibitors [11,12,13].Previously, we identified rhodanine-containing compounds as IN inhibitors through a common feature pharmacophore search based on β−diketoacid bioisosteres (Figure 1) [14]. We reported that compounds containing both a rhodanine moiety and a salicylic acid substitution exhibited significant IN inhibition. The most potent compound 1 inhibited IN catalytic activities with IC50 values of 15 μM and 11 μM for 3’-processing and strand transfer, respectively. This prompted us to explore other modifications in the substructures ofcompound 1 (Figure 1) and study their effects on IN inhibition. Another earlier study also similarly applied a caffeic acid phenylester (CAPE)-based pharmacophore approach and identified rhodanine-containing compounds bearing a phenylsulfonamide moiety as IN inhibitors [15].
Figure 1
Design of new rhodanine derivatives and structural modifications for SAR study.
Several rhodanine-based small molecules have also known to inhibit various enzymes such as dual-specificity phosphatases, HCV NS3 protease, phosphodiesterases, etc. (Figure 2) [16,17,18,19]. The central rhodanine scaffold is commonly found in several antivirals, antimicrobials and antitumor agents [20,21]. In all these structures, modifications in the alkylidenechain and in the rhodanine moiety appear to influence biological activity. While there is diverse literature on rhodanine-containing compounds as various enzyme inhibitors [22], there is limited information regarding their use as HIV-1 IN inhibitors. Herein, we report the synthesis and structure-activity relationship (SAR) for a series ofrhodanine derivatives that inhibit IN in the low micromolar range. We further characterized the selectivity of these compounds for IN by counter-screening against APE1, another DNA binding enzyme, and analyzed their binding modes in the IN and APE1 active sites. Finally, we evaluated the antiproliferative effects of these rhodanines against a panel ofhumancancercell lines.
Figure 2
Structures of representative rhodanine derivatives with diverse biological activities [17,18,19,20]; (2) JSP-1 inhibitor; (3) PDE4 inhibitor; (4) HCV NS3 Protease inhibitor; (5) Antitumor agent (rhodanine moiety is highlighted in red).
Design of new rhodanine derivatives and structural modifications for SAR study.Structures of representative rhodanine derivatives with diverse biological activities [17,18,19,20]; (2) JSP-1 inhibitor; (3) PDE4 inhibitor; (4) HCV NS3 Protease inhibitor; (5) Antitumor agent (rhodanine moiety is highlighted in red).
2. Results and Discussion
2.1. Chemistry
In our earlier studies, we identified rhodanine-containing compounds through a pharmacophore-based high-throughput screening. We found that the presence of both a rhodaninecore and a salicylic acid group improved IN inhibitory effects [14]. Furthermore, substitutions in the aryl substructure and at the 4-position in the rhodanine nucleus (alkylidene substructure) also appeared to influence IN inhibition potency. Based on this, we synthesized several 2-thioxo-1,3-thiazolidin-4-ones bearing modifications in both the aryl and the alkylidene substructures to explore their effect on IN catalytic activities and understand their SAR. These rhodanines were synthesized according to previously a reported procedure (Scheme 1) [23]. All of these new analogues possessed a five-membered rhodanine moiety as a common structural unit and are linked to various substituted aromatic groups via alkylidene and four-carbon aliphaticthioxoamide linkers (Figure 1).
Scheme 1
Synthesis of compounds 6–54.
Synthesis ofcompounds 6–54.
2.2. In vitro IN inhibitory activities and SAR of rhodanine derivatives
Inhibition of IN catalytic activities, 3’-processing and strand transfer, was determined using an in vitro enzymatic assay. Compounds with a weakly electron-withdrawing chlorophenyl group in substructure B (e.g., 6–8) showed no IN inhibition at the highest concentration tested (Table 1). Interestingly, compound 9, which has a methoxy substitution in ring A, inhibited both 3’-processing and strand transfer with an equal potency (IC50 value of 33 μM). Substitution with the weakly basic and bulky indole group abolished any activity (compound 10). Compound 11, with a stronger electron-withdrawing substitution in substructure B, had a greater potency towards IN inhibition with an IC50 value of 58 ± 4 μM and 20 ± 6 μM for 3’-processing and strand transfer, respectively.
Table 1
Inhibition of HIV-1 IN and APE-1 catalytic activities by rhodanines 6–11.
Compound
Structure
Inhibition of IN catalytic activity IC50 (µM)
Inhibition of APE1 catalytic activityIC50 (μM)
3’-Processing
Strand Transfer
6
>100
>100
>100
7
>100
>100
45 ± 51
8
>100
>100
76 ± 5
9
33
33
47 ± 2
10
>100
>100
>100
11
58 ± 4
20 ± 6
>100
Inhibition ofHIV-1 IN and APE-1 catalytic activities by rhodanines 6–11.Compounds 12–21 have a relatively stronger electron-withdrawing nitrofuran ring in substructure B. Again, we observed that electron-withdrawing groups in ring A led to moderate inhibitory activity, while replacement by an electron-donating group such as a methoxyl improved activity further (Table 2).
Table 2
Inhibition of HIV-1 IN and APE1 catalytic activities by rhodanines 12–21.
Compound
Structure
Inhibition of IN catalytic activity IC50 (µM)
Inhibition of APE1 catalytic activity IC50 (μM)
3’-Processing
Strand Transfer
12
100
93
89 ± 13
13
68 ± 46
41 ± 27
93
14
58
58
>100
15
61 ± 35
56 ± 11
-
16
>100
74 ± 20
>100
17
77 ± 10
60
>100
18
34 ± 10
14 ± 5
>100
19
40 ± 28
25 ± 11
>100
20
>100
85 ± 21
>100
21
>100
>100
>100
Inhibition ofHIV-1 IN and APE1 catalytic activities by rhodanines 12–21.The presence of a basic moiety such as tetrazole ring at position 4 in ring A led to a loss of activity. Interestingly, introduction of electron-donating group in substructure B did not seem favorable for IN enzyme binding, as all of the compounds 22–26 were inactive (Table 3).
Table 3
Inhibition of HIV-1 IN and APE1 catalytic activities by rhodanines 22–26.
Compound
Structure
Inhibition of IN catalytic activity IC50 (µM)
Inhibition of APE1 catalytic activityIC50 (μM)
3’-Processing
Strand Transfer
22
>100
90
>100
23
>100
>100
83
24
>100
>100
>100
25
>100
>100
>100
26
>100
>100
-
Inhibition ofHIV-1 IN and APE1 catalytic activities by rhodanines 22–26.Since the presence of an electron-donating methoxy group in ring A seemed important for IN inhibitory activity, we synthesized analogs with trimethoxy substitution on ring A (Table 4). Substructure B was also modified using various substituted phenols. Electron-withdrawing substituents like Cl, Br and NO2 led to increased activity. Compound 35 with a salicylic acid had only a moderate effect on IN inhibition, while compound 34 with a 2-phenyldiazaphenol moiety inhibited 3’-processing and strand transfer, with IC50 values of 33 ± 19 μM and 26 ± 14, respectively. Aliphatic and basic substitutions had a negative effect on activity.
Table 4
Inhibition of HIV-1 IN and APE1 catalytic activities by rhodanines 27–40.
Compound
Structure
Inhibition of IN catalytic activity IC50 (µM)
Inhibition of APE1 catalytic activityIC50 (μM)
3’-Processing
Strand Transfer
27
83 ± 25
20
>100
28
44 ± 13
30 ± 19
>100
29
74 ± 18
72 ± 23
>100
30
47 ± 30
>100
>100
31
>100
>100
>100
32
82 ± 5
83 ± 3
>100
33
>100
>100
>100
34
33 ± 19
26 ± 14
>100
35
45 ± 24
61 ± 11
>100
36
46 ± 5
>100
>100
37
33 ± 7
35 ± 6
>100
38
51 ± 25
32 ± 11
65 ± 33
39
47 ± 20
31 ± 16
>100
40
>100
95
>100
Inhibition ofHIV-1 IN and APE1 catalytic activities by rhodanines 27–40.In order to further understand the effect of electron-donating groups in ring A on activity, we synthesized analogs with hydroxyl groups at position 4. Here again, strong electron-withdrawing and hydrophobic substitutions in substructure B appears to influence IN inhibition (Table 5). The most potent IN inhibition was seen with compound 53 having 3,5-diiodophenol substitution (IC50 value of 7 ± 3 μM and 3 ± 2 for 3’-processing and strand transfer, respectively). Compared to compounds bearing trimethoxy substitution (Table 4), some of these compounds are more potent, suggesting that electron-donating substituents in ring A, together with strong electron-withdrawing and hydrophobic substitution in substructure B, leads to favorable IN inhibitory activity. Interestingly, substitution with a hydroxyl group in ring A improved selectivity towards strand transfer inhibition (Figure 3A). This may be due to possible formation of stabilizing hydrogen bonds at the IN active site. Overall, these rhodanine derivatives with an aliphaticthioxoamide linker have an improved IN inhibitory activity over those previously reported with an aliphatic linker [14].
Table 5
Inhibition of HIV-1 IN and APE1 catalytic activities by rhodanines 41–54.
Compound
Structure
Inhibition of IN catalytic activity IC50 (µM)
Inhibition of APE1 catalytic activityIC50 (μM)
3’-Processing
Strand Transfer
41
>100
67
>100
42
>100
>100
>100
43
>100
88
>100
44
53 ± 8
45 ± 14
>100
45
>100
>100
>100
46
>100
100
>100
47
>100
>100
>100
48
17 ± 7
8 ± 2
>100
49
80
45
>100
50
60
55
>100
51
28 ± 1
21 ± 7
>100
52
38 ± 27
23 ± 9
>100
53
7 ± 3 (12.5)a
3 ± 2 (11)a
62 ± 3
54
>100
>100
>100
a Values in the parenthesis were obtained in the presence of Mg2+ as the cofactor.
Figure 3
(A) A representative gel showing inhibition of purified IN by selected compounds. Lanes 1: DNA alone; Lanes 2 and 15: DNA and IN alone with no drug; Lanes 3–14: DNA with IN and selected drug concentrations (lanes: 3–6: compound 48; lanes 7–10: compound 53; lanes 11–14: compound 51; at concentrations: 100, 33.3, 11.1 and 3.7 μM); (B) A representative gel showing inhibition of purified APE1 by selected compounds. Lanes 1: DNA alone; Lanes 2 and 15: DNA and APE1 alone with no drug; Lanes 3–26: DNA with APE1 and selected drug concentrations (lanes: 3–6: compound 48; lanes 7–10: compound 53; lanes 11–14: compound 51; lanes 15–18: compound 9; lanes 19–22: compound 7; lanes 23–26: compound 8; at concentrations: 100, 33.3, 11.1 and 3.7 μM).
Inhibition ofHIV-1 IN and APE1 catalytic activities by rhodanines 41–54.a Values in the parenthesis were obtained in the presence ofMg2+ as the cofactor.(A) A representative gel showing inhibition of purified IN by selected compounds. Lanes 1: DNA alone; Lanes 2 and 15: DNA and IN alone with no drug; Lanes 3–14: DNA with IN and selected drug concentrations (lanes: 3–6: compound 48; lanes 7–10: compound 53; lanes 11–14: compound 51; at concentrations: 100, 33.3, 11.1 and 3.7 μM); (B) A representative gel showing inhibition of purified APE1 by selected compounds. Lanes 1: DNA alone; Lanes 2 and 15: DNA and APE1 alone with no drug; Lanes 3–26: DNA with APE1 and selected drug concentrations (lanes: 3–6: compound 48; lanes 7–10: compound 53; lanes 11–14: compound 51; lanes 15–18: compound 9; lanes 19–22: compound 7; lanes 23–26: compound 8; at concentrations: 100, 33.3, 11.1 and 3.7 μM).
2.3. In vitro APE1 Inhibitory Activity
The specificity of these rhodanine derivatives against IN was assessed by counter-screening against humanapurinic/apyrimidinic endonuclease 1 (APE1), another DNA-binding enzyme. APE1 plays an important role in DNA repair by the base excision pathway. Like IN, APE1 possesses endonucleolytic activity and cleaves the phosphodiester DNA backbone 5’ to its recognition site. Like IN and other DNA-binding enzymes, optimal APE1catalytic activity requires the presence of a divalent metal ion. APE1 is also an attractive target for designing cancer therapeutics [24,25,26].Most compounds showed no inhibition against APE1 at the maximum tested concentration of 100 μM. Compound 53, with the most potent IN inhibitory activity, exhibited only a weak APE1 inhibition with an IC50 value of 62 ± 3 μM. Its inhibitory activity is 20-fold more selective for IN than APE1. Similarly, compounds 48 and 51, with promising IN inhibitory activity, did not show any significant APE1 inhibition at 100 μM. On the other hand, two compounds, 7 and 8, showed a moderate to weak APE1 inhibition with IC50 values of 45 ± 51 μM and 76 ± 5 μM, respectively, while they lacked any IN inhibitory activity. Similarly, compound 9 inhibited APE1 activity with an IC50 value of 47 ± 2 μM (Figure 3B). However, it also exhibited an equal inhibitory effect on IN (IC50 value of 33 μM), suggesting possible similarities in its mechanism of action.
2.4. Docking Studies
To predict the possible binding modes and enzyme inhibition mechanism, compounds 6–54 were docked onto the active sites of both HIV-1 IN and APE1, using GOLD 4.0, the automated docking program, and Glide, the grid-based ligand docking with energetics software from Schrödinger. The docking scores for these compounds are reported in Table 6. A plot of the pIC50 (ST activity) versus the predicted Glide docking score of these 49 compounds against IN is shown in Figure 4. Although a significant correlation between docking scores and inhibitory activities of the compounds was not routinely observed, there was a general trend where many of the active compounds scored high. Of the 32 active compounds with a ST IC50 value less than 100 µM, 24 compounds were correctly predicted as true positives. Eight compounds, which were moderately active or least active, were predicted as false negatives. Similarly, out of the 18 inactive compounds with ST IC50 values greater than 100 µM, 9 compounds were correctly predicted as true negatives, while 9 inactive compounds were over predicted as false positives. We found no correlation between the APE1 inhibitory activity and the docking scores.
Table 6
Comparison of docking scores of rhodanine-based compounds against HIV-1 IN and APE1.
Compound
HIV-1 IN ST IC50
APE1 IC50
Molecular Docking against HIV-IN
Molecular Docking against APE1
(μM)
(μM)
GOLD Fitness
Glide Score
GOLD Fitness
Glide Score
6
>100
>100
44.52
−3.26
54.65
−3.39
7
>100
45 ± 51
43.22
−3.42
56.72
−3.67
8
>100
76 ± 5
43.55
−3.24
50.35
−4.22
9
33
47 ± 2
42.9
−3.53
48.98
−3.60
10
>100
>100
44.44
−3.90
50.16
−5.63
11
20 ± 6
>100
41.15
−4.12
54.13
−5.31
12
93
89 ± 13
50.61
−3.80
55.33
−4.17
13
41 ± 27
93
56.35
−3.62
56.55
−4.44
14
58
>100
53.7
−3.78
61.16
−4.47
15
56 ± 11
−
50.97
−3.90
61.76
−3.90
16
74 ± 20
>100
49.48
−3.54
64.47
−3.21
17
60
>100
53.83
−3.59
51.32
−4.94
18
14 ± 5
>100
53.72
−4.19
48.2
−4.11
19
25 ± 11
>100
49.03
−4.20
48.36
−4.33
20
85 ± 21
>100
44.77
−3.44
57.11
−4.09
21
>100
>100
48.46
−3.71
53.61
−3.82
22
90
>100
36.49
−4.32
47.64
−4.67
23
>100
83
39.77
−3.44
52.31
−4.50
24
>100
>100
44.93
−3.96
50.75
−4.11
25
>100
>100
46.72
−2.82
49.01
−4.22
26
>100
−
47.85
−2.10
49.75
−4.46
27
20
>100
47.43
−5.03
46.7
−5.82
28
30 ± 19
>100
45.6
−5.37
47.93
−4.78
29
72 ± 23
>100
39.57
−5.45
49.09
−5.72
30
>100
>100
44.27
−5.33
49.84
−5.03
31
>100
>100
37.03
−5.24
47.14
−5.00
32
83 ± 3
>100
45.69
−5.33
47.69
−5.74
33
>100
>100
29.36
−5.29
39.69
−5.32
34
26 ± 14
>100
46.9
−4.98
52.47
−5.35
35
61 ± 11
>100
38.25
−4.97
47.61
−5.21
36
>100
>100
43.69
−5.35
46.27
−5.35
37
35 ± 6
>100
45.35
−5.17
49.67
−4.34
38
32 ± 11
65 ± 33
45.18
−5.10
49.09
−5.43
39
31 ± 16
>100
44.6
−5.08
50.75
−4.36
40
95
>100
21.11
−5.28
14.58
−3.79
41
67
>100
44.68
−5.35
52.5
−5.82
42
>100
>100
45.63
−5.29
49.24
−4.52
43
88
>100
50.44
−5.80
49.42
−6.18
44
45 ± 14
>100
44.75
−5.36
49.86
−5.08
45
>100
>100
42.79
−5.53
49.77
−5.20
46
100
>100
48.8
−5.80
52.75
−4.85
47
>100
>100
37.43
−5.54
45.29
−4.99
48
8 ± 2
>100
49.67
−5.16
59.85
−6.65
49
45
>100
44.68
−5.45
54.27
−7.48
50
55
>100
44.04
−5.74
51.34
−5.91
51
21 ± 7
>100
44.19
−5.64
54.26
−5.78
52
23 ± 9
>100
43.43
−5.73
50.25
−5.92
53
3 ± 2
62 ± 3
41.53
−5.72
55.04
−4.79
54
>100
>100
16.79
−5.40
17.8
−4.34
Figure 4
Plot of pIC50 (IN ST activity) versus the predicted Glide docking score.
Plot of pIC50 (IN ST activity) versus the predicted Glide docking score.Comparison of docking scores ofrhodanine-based compounds against HIV-1 IN and APE1.Compounds 6–8, with a weak electron-withdrawing group on ring B, did not show any IN inhibitory activity at the highest tested concentration (Table 1). However, the introduction of a carbonyl group in substructure B (11) significantly improved strand transfer activity. Interestingly, docking simulations with Glide also estimated a high fitness score for this compound. The observed activity of the compound 11 supports the predicted binding interactions inside the IN active side. As shown in Figure 5, the oxygen atom of the rhodaninecore and the oxygen atom of the substructure ring B are involved in interactions with both the Mg2+ ion and the active site residues D64 and D116. Additionally, the donornitrogen atom of 11 forms a hydrogen bond with E152. The importance of this additional carbonyl of 11 in substructure B can further be understood by comparing against the activity and binding mode ofcompound 9 (Figure 5). A possible explanation for the lower activity of 9, as compared to 11, might be the lack of such interactions with key active residues D64 and D116 in the substructure B. Compounds 6–8 show a similar binding mode at the HIV-1 IN active site.
Figure 5
The proposed binding modes of (a) compound 11, and (b) compound 6–8 in the HIV-1 IN active site. Close contact residues are shown as green sticks. Compounds 11 and 6–8 are shown in pink (colored by atom types) and Mg2+ is shown in red.
Compounds 12–21 with an electron-withdrawing nitrofuran ring in substructure B exhibit a moderate inhibition of strand IN transfer. Docking results predicted seven of these moderately active (12–17 and 20) compounds as false negatives, with their predicted activities lower than their actual activities. Compounds 27–40, with electron-donating trimethoxy group on ring A, have shown moderate IN inhibitory activity. Docking results predicted four inactive compounds (30–31, 33 and 36) as false positives, with their predicted activities being higher than their actual activities.The proposed binding modes of (a) compound 11, and (b) compound 6–8 in the HIV-1 IN active site. Close contact residues are shown as green sticks. Compounds 11 and 6–8 are shown in pink (colored by atom types) and Mg2+ is shown in red.The binding mode of the most active compound 53 is shown in Figure 6. The residues in close contact with the compound 53 were D64, C65, T66, H67, D116, I141, Q148, E152, and N155. As shown in Figure 6, the carbonyl oxygen atom of the rhodanine and the hydroxyl oxygen atom in ring B are involved in interactions with both the Mg2+ ion and the amino acid residues D64 and C65. The two donornitrogen atoms of 53 formed strong H-bonding interactions with the key active site residues D64 and E152. Additionally, the hydroxyl oxygen atom of the ring A was involved in a strong H-bonding interaction with the amino acid residue I141. The presence of these interactions with key active site residues explains the potency of 53. A close analogue to 53 is compound 39. The replacement of p-hydroxyl group (53) in ring A by a trimethoxy group (e.g., compound 39) resulted in a 10-fold decrease in potency. The binding mode ofcompound 39 is shown in Figure 6. A possible explanation for its lower activity, as compared to 53, might be the lack of H-bonding interactions with the amino acid residue I141. These docking results confirm the observations from the SAR studies. Indeed, the p-hydroxy group in ring A is able to form stabilizing H-bonds with I141. This accounts for the increased potency ofcompounds 41–54 as compared to their trimethoxy analogues 27–40.
Figure 6
The proposed binding modes of (a) compound 39; (b) compound 53; and (c) compound 54 in the HIV-1 IN active site. Close contact residues are shown as green sticks. Compounds 39, 53, and 54 are shown in pink (colored by atom types) and Mg2+ is shown in red.
Orientation of the rhodaninecore and substructure B are also important for IN inhibitory activity, contributing to the formation of key interactions with the metal ion and active site residues. Accordingly, replacement of the iodine groups on ring B by t-butyl groups (54) led to a drasticchange in the binding pattern of ring B and a loss of activity. Though Glide estimated similar scores for both 53 and 54, differences in activity can be clearly explained from the GOLD docking poses (Figure 6). With the introduction of the bulky t-butyl groups, the rhodanine ring becomes inverted and no longer interacts with the Mg2+ ion and the active site residues. This might be responsible for the loss of activity ofcompound 54.The proposed binding modes of (a) compound 39; (b) compound 53; and (c) compound 54 in the HIV-1 IN active site. Close contact residues are shown as green sticks. Compounds 39, 53, and 54 are shown in pink (colored by atom types) and Mg2+ is shown in red.
2.5. Drug-like Properties
Oral bioavailability is a desirable property ofcompounds under investigation in the drug discovery process. Lipinski’s rule-of-five is a simple model to forecast the absorption and intestinal permeability of a compound [27,28]. In the rule-of-five model, compounds are considered likely to be well absorbed when they possess these attributes—molecular weight < 500, cLogP < 5, number of H-bond donors < 5, number H-bond acceptors < 10, and number of rotatable bonds < 10. Further, we calculated the polar surface area, which is also used to predict drug absorption [29,30]. Earlier studies by Palm et al. [31,32,33] and Kelder et al. [34] suggest that the compounds with a polar surface area >140 Å2 would tend to show poor (<10%) absorption, whereas the compounds with polar surface area <60 Å2 could be predicted to show complete (>90%) absorption. The calculated atom-based Log P (S + logP) values of most of the active compounds ranged from 2.0 to 5.0, and the H-bond donor and acceptor counts are ≤5 and <10, respectively. High lipophilicity is often an issue with many rhodanine-containing structures impeding their cell permeability and in vivo activity [22]. However, most of the compounds presented here have polar surface areas <140 Å2 and log P value of <5. Therefore, most of these analogs have favorable drug-like properties and are potentially interesting for further optimization.Predicted physicochemical and ADMET properties ofrhodanine-based compounds.MWt: Molecular weight; MolVol: Molal volume at normal boiling point; MlogP: Moriguchi octanol-water partition coefficient; S + logP: Simulation plus octanol-water partition coefficient; S + logD: Simulation plus octanol-water distribution coefficient; S + Vd: Simulation plus volume of distribution; S + BBB: Simulation plus BBB permeability; PSA_2D: 2D-polar surface area; PolASA_3D: 3D-polar solvent accessible surface area; HBD: Hydrogen bond donors; HBA: Hydrogen bond acceptors.
2.6. Cytotoxicity and Antiviral Activity
Several rhodanine structures such as the rhodacyanine dyes have been found to possess anticancer and antitumor properties. Hence, we evaluated the ability of these rhodanine-based compounds to inhibit the proliferation of various cancercell lines. Compounds 13 and 17 were moderately active against a colon cancercell line (HCT 116) and a pancreaticcell line (Panc-1) (GI50 values provided in Table 8).
Table 8
In vitro cytotoxicity of rhodanine-based derivatives against a panel of cancer cell lines.
Compounds
GI50 (μM) a
MDA-MD435
Panc-1
HT 29
HCT 116p53+/+
HCT 116p53−/−
13
13.5 ± 5.2
8.6 ± 1.6
>10
>10
6.5 ± 1.2
17
9 ± 1.1
8.2 ± 1.2
9.6
5.1 ± 0.8
6.1 ± 2.2
a Values are reported as GI50 in μM, the concentration of the compound required to cause 50% inhibition of cell proliferation. GI50 values are calculated from at least three independent experiments at multiple doses.
In vitro cytotoxicity ofrhodanine-based derivatives against a panel ofcancercell lines.a Values are reported as GI50 in μM, the concentration of the compound required to cause 50% inhibition ofcell proliferation. GI50 values are calculated from at least three independent experiments at multiple doses.The remaining compounds lacked any significant cytotoxicity at the maximum tested concentration of 10 μM. We also tested the antiviral activity ofcompounds 48, 51 and 53 against HIV-1 inMT-4cell infectivity assay. Unfortunately, these compounds did not show any significant antiviral activity. Optimization studies are underway to design rhodanines with more favorable physiochemical properties and improved antiviral activity.
3. Experimental
3.1. General
The 1H-NMR spectra were recorded on a Bruker AC spectrometer (200 MHz) using DMSO-d6 as solvent. The data are given in parts per million (ppm) and are referenced to an internal standard oftetramethylsilane (TMS, δ 0.00 ppm). The spin-spin coupling constants (J) are given in Hz. Peak multiplicity is reported as s (singlet), d (doublet), dd (double doublet), t (triplet), m (multiplet), and br s (broad singlet). The mass spectra were obtained on a Kratos instrument using a direct inlet system; the ionization energy was 70 eV; the accelerating voltage was 1.75 kV. Melting points were measured on a Boetius hot-stage apparatus. 5-Benzylidene-4-oxo-2-thioxothiazolidines 6–54 were purified by silica gel flash chromatography column (eluent: 3:1 petroleum ether-ethyl acetate). 2-Thioxo-1,3-thiazolidin-4-ones were prepared according to a previously reported procedure (Scheme 1) [23].
3.2. General Procedures for the Synthesis of 5-Benzylidene-4-oxo-2-thioxothiazolidine (6–54)
All compounds were dissolved in DMSO to a final stock concentration of 10 mM and stored at −20 °C. Initial screening for IN inhibitory effect was done at 100 μM and dilutions were made in DMSO. Compounds showing at least 50% inhibition at 100 μM were selected for IC50 determinations. γ-[32P]-ATP was purchased from Perkin Elmer (Waltham, MA, USA). The expression system for wild-type IN was a generous gift of Dr. Robert Craigie, Laboratory of Molecular Biology, NIDDK, NIH (Bethesda, MD, USA). His-tagged APE1 was purified from an expression system kindly provided by Dr. TomCurran, Department of Developmental Neurobiology, St. Jude Children’s Research Hospital (Memphis, TN, USA).
3.3.2. Preparation of Oligonucleotide Substrates
21-mer oligonucleotides used in IN inhibition assay [21top (5’-GTGTGGAAAATCTCTAGCAGT-3’) and 21bot (5’-ACTGCTAGAGATTTTCCA CAC-3’)] were purchased from Integrated DNA Technologies (Coralville, IA, USA) and purified by UV shadowing on polyacrylamide gel. The oligonucleotides used in APE1 inhibition assay [top strand: 5′-ATTTCACCGGTACG(F)TCTAGAATCCG-3′ containing the tetrahydrofuran synthetic abasic residue (F), and bottom strand: 3′-TAAAGTGG CCATGC(C)AGATCTTAGGC-5′] were purchased from The Midland Certified Reagent Company (Midland, TX, USA). To analyze the extent of 3’-processing and strand transfer with 5’-end labeled substrates, 21top was 5’-end labeled using T4 polynucleotide kinase (Epicentre, Madison, WI, USA) and [γ-32P]-ATP (Perkin Elmer, Waltham, MA, USA). The kinase was heat-inactivated at 95 °C and 21bot was added in 1.5 M excess. The reaction mixture was allowed to cool slowly to room temperature, and purified through a spin-25 minicolumn (USA Scientific, Ocala, FL, USA) to separate annealed double-stranded oligonucleotidefrom unincorporated material. The oligonucleotides used in the APE1 inhibition assay were labeled using a similar protocol.
3.3.3. In Vitro IN Inhibition Assay
To determine the extent of 3’-processing and strand transfer, the test compounds, diluted in DMSO, were pre-incubated with purified recombinant IN (at a final concentration of 200 nM) in the reaction buffer [50 mM NaCl, 1 mM HEPES (pH 7.5), 50 μM EDTA, 50 μM dithiothreitol, 10% glycerol (w/v), 7.5 mM MnCl2, 0.1 mg mL–1 bovine serum albumin, 10 mM 2-mercaptoethanol, 10% dimethyl sulfoxide, and 25 mM MOPS (pH 7.2)] at 30 °Cfor 30 min. For Mg2+ assay, 7.5 mM MgCl2 was substituted for MnCl2 and an optimal concentration ofPEG was included in the reaction buffer. Following the pre-incubation, 20 nM of the 5’-end 32P-labeled linear oligonucleotide substrate was added to the reaction mixture and incubated for an additional 1 h. Reactions were quenched by addition of 50X loading dye (98% deionized formamide, 10 mM EDTA, 0.025% xylene cyanol, and 0.025 % bromophenol blue). An aliquot was subjected to electrophoresis on a denaturing polyacrylamide gel (0.09 M tris–borate, pH 8.3, 2 mM EDTA, 20% acrylamide, 8M urea). Gels were dried under vacuum, exposed in a PhosphorImager cassette and visualized with a Typhoon 8610 Variable Mode Imager (Amersham Biosciences). Quantification was done with Image Quant 5.2 software. Percent inhibition (%I) was calculated using the following equation:
%I = 100 × [1 − (D − C) / (N − C)]
where C, N, and D are the fractions of 21-mer substrate converted into 19-mer (product of 3’-processing) or strand-transfer products for DNA alone, DNA plus IN, and DNA, IN plus test compound, respectively. IC50 values were determined by plotting the logarithm of test compound concentrations as a function of %I to obtain the concentration that produced 50% inhibition.
3.3.4. APE1 Inhibition Assay
APE1 inhibition assay was performed as previously described [26]. To determine the extent of endonucleolyticcleavage of the abasic residue, the test compounds, diluted in DMSO, were preincubated with recombinant APE1 at a final concentration of 0.05 nM in the reaction buffer [50 mM NaCl, 1 mM HEPES, pH 7.5, 50 μM EDTA, 50 μM dithiothreitol, 10% glycerol (w/v), 7.5 mM MnCl2, 0.1 mg mL‑1 bovine serum albumin, 10 mM 2-mercaptoethanol, 10% DMSO, and 25 mM MOPS, pH 7.2] at 30 °Cfor 10 min. Following preincubation, 200 nM of the 5′-end 32P-labeled linear oligonucleotide substrate was added, and incubated for an additional 10 minutes. Reactions were quenched by the addition of an equal volume (16 μL) of loading dye (98% deionized formamide, 10 mM EDTA, 0.025% xylene cyanol, and 0.025% bromophenol blue). An aliquot (5 μL) was electrophoresed on a denaturing 20% polyacrylamide gel (0.09 M tris-borate pH 8.3, 2 mM EDTA, 20% acrylamide, 8 M urea). The gels were dried, exposed in a PhosphorImager cassette and analyzed using a Typhoon 8610 Variable Mode Imager (Amersham Biosciences). Quantification of enzyme inhibition was done using Image Quant 5.2 software. The percent inhibition (% I) was calculated using the following equation:
%I = 100 × [1 − (D − C) / (N − C)]
where C, N, and D are the fractions of 26-mer substrate converted to 13-mer incision products (or integration products) for DNA alone, DNA plus APE1, and DNA, APE1 plus test compound, respectively. IC50 values were determined by plotting the logarithm ofcompound concentration as a function of %I to obtain the concentration that produced 50% inhibition. An optimumconcentration ofAPE1 enzyme, at which there is complete conversion of the synthetic abasic-site containing oligonucleotide substrate to its cleaved product without subsequent exonuclease activity, was used in the assay.
3.4. Docking Studies
3.4.1. GOLD
HIV-1 IN. Structures of all the compounds 6–54 from Table 1–Table 5 were built and energy minimized using Catalyst (Accelrys, Inc., San Diego, CA, USA) [35]. Using the Poling algorithm and the best flexible conformation generation method in Catalyst, feasible and unique conformations for each compound was generated over a 20 kcal/mol range of energies. Three-dimensional coordinates of the crystallized structure of the catalytic domain ofHIV-1 IN were obtained from the Protein Data Bank at the Research Collaboratory for Structural Bioinformatics (RCSB PDB; PDB ID: 1BL3). All the water molecules present in protein were removed. Hydrogen atoms were added to the protein to assign appropriate ionization states to both the acidic and basic amino acid residues in the active site. Docking was performed using version 4.0 of the GOLD (Genetic Optimization for Ligand Docking; Cambridge Crystallographic Data Centre) software package [36]. A 20 Å radius, centered at the carboxylateoxygen atom (OD1) of the amino acid residue D64, within the IN active site was defined as the binding site. The docking simulations were carried out using the standard default settings—a population size of 100, a maximum number of 100,000 operations, and a mutation and crossover rate of 95. The GOLD Score fitness function is the default scoring function provided with GOLD. It has been optimized for the prediction of ligand binding positions and takes into account factors like H-bonding energy, Van der Waals energy and ligand torsion strain. On the other hand, the ChemScore fitness function estimates the total free energy change that occurs on ligand binding and is trained by regression against binding affinity data. All the compounds (6–54) were docked into the defined IN active site and the GOLD fitness scores were calculated. For each molecule, the bound conformation with a high fitness score was considered as the best bound conformation.APE1. Three-dimensional coordinates of the crystallized structure ofhumanAPE1 endonuclease with bound abasic DNA and Mn2+ ion were obtained from RCSB PDB (PDB ID: 1DE9). For the protein target, a 20 Å radius active site was defined using xyz coordinates of the 5’-phosphate fragment of the abasic DNA in the co-crystal structure. Prior to docking, the bound abasic DNA and chain B of the APE1 monomer were removed, and proper protonation states were assigned for the acidic and basic residues. Docking simulations with standard parameters and GOLD fitness score calculations were performed as described above. Based on the GOLD fitness scores for each molecule, a bound conformation with high fitness score was considered as the best bound-conformation.
3.4.2. Glide
Compounds 6–54 from Table 1-Table 5 were processed using the Schrödinger LigPrep utility (Schrödinger, LLC, USA). This program produces a number of low energy 3D structures, with various ionization states, tautomers, stereochemistries, and ring conformations, from each molecule input. For this study, a pH range of 6-8 was used. Both the neutral and the anionic states ofcarboxylic acid groups were generated. Three-dimensional coordinates of the crystal structures ofHIV-1 IN catalytic domain and humanAPE1 endonuclease, with bound abasic DNA and Mn2+ ion, were obtained from the RCSB PDB (PDB ID: 1BL3 and 1DE9). All the water molecules present in the protein were removed, and hydrogen atoms were added considering the appropriate ionization states for both the acidic and basic amino acid residues. All compounds (6–54) were docked onto the binding site in HIV-1 IN or humanAPE1 protein using version 5.0 of Glide (Grid-Based Ligand Docking With Energetics) software from Schrödinger [37,38,39]. A grid with a bounding box size of 20 Å was used. The coordinates of the enclosing box (1BL3: x = 22.23 Å; y = 1.33 Å; z = −19.16 Å and 1DE9: x = 47.866; y = 9.534; z = 40.274) were defined starting from the set of active site residues involved. The Glide algorithm is based on a systematic search of positions, orientations, and conformations of the ligand in the receptor binding site using funnel type approach. The search begins with a rough positioning and scoring phase. This significantly limits the search space, and reduces the number of poses to be selected for minimization on the precomputed OPLS-2005 Van der Waals and electrostatic grids for the protein. In order to obtain an accurate correlation between good poses and good scores, the Glide Extra-Precision (XP) Mode was used.
3.5. Drug-Like Property Prediction
The ligand structures were energy minimized using the Catalyst software package (Accelrys, Inc.) [35] and the lowest- energy conformation of each compound was exported to ADMET Predictor (Simulations Plus, Inc., Lancaster, CA, USA) [40] to calculate various ADME (absorption, distribution, metabolism and excretion) properties.
3.6. Biological Activity
3.6.1. Cell Lines and Chemicals
Humanmelanomacell line MDA-MB 435, pancreatic cancercell line Panc-1, and colon cancercell line HT 29 were purchased from the American Type Cell Culture (Manassas, VA, USA). HCT 116 p53+/+ and HCT 116 p53−/− were kindly provided by Dr. Bert Vogelstein, Johns Hopkins Medical Institutions (Baltimore, MD, USA). HCT 116 p53+/+, HCT116p53−/−, and HT 29 cells were maintained as monolayer cultures in RPMI 1640 supplemented with 10% fetal bovine serum (FBS; Gemini-Bioproducts, Woodland, CA, USA) and 2 mmol/L L-glutamine at 37 °C in a humidified atmosphere of 5% CO2. Similarly, MDA-MB 435 and Panc-1cells were maintained in Dulbecco’s modified Eagle’s medium (DMEM) supplemented with 10% fetal bovine serum as described above. To remove the adherent cells from the flask for passaging and counting, cells were washed with PBS without calcium or magnesium, incubated with a small volume of 0.25% trypsin-EDTA solution (Sigma-Aldrich, St. Louis, MO, USA) for 5 to 10 minutes and washed with culture medium and centrifuged. All experiments were done using cells in exponential cell growth.
3.6.2. Cytotoxicity Assay
Cytotoxicity was assessed by 3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide (MTT) assay as previously described [41,42]. Briefly, cells were seeded in 96-well microtitre plates (Panc-1 and HT 29 cells were seeded at 8,000 cells/well, and MDA-MB-435 and HCT116cells at 4,000 cells/well, respectively). After overnight attachment, cells were treated with the test compounds for 72 h. Initial screening in MTT assays was performed at 10 μM using sterile drug solutions prepared in suitable cell culture media (final DMSOconcentration in drug dilutions was less 0.001%). Compounds showing at least 50% growth inhibition at 10 μM were selected for IC50 determinations.After 72 h exposure to the test compounds, MTT solution (5 mg/mL; 20 μL) was added to each well and cells were incubated for 4 h at 37 °C. After incubation, media from each well was removed and the dark blue formazancrystals formed by live cells were dissolved in DMSO (150 μL/well). The absorbance intensity was measured at 570 nm against appropriate blank controls using a microplate reader (Molecular Devices, Sunnyvale, CA, USA). All assays were done in triplicate. The growth inhibitory concentration GI50 was determined by plotting the logarithm ofcompound concentration against percentage ofcells killed to obtain the concentration that produced 50% cell kill.
4. Conclusions
In this study, we have investigated the SAR ofrhodanine-based IN inhibitors by modifying the aryl and alkylidene sub-groups. Several of these compounds showed moderate to potent IN inhibition and appeared to be more selective for IN as compared to APE1. Results from this study confirm the importance of the rhodanine moiety as a useful structural scaffold with promising IN inhibitory activity and provides a framework for designing more potent and selective IN inhibitors. In addition, some of these compounds inhibited APE1 endonuclease activity and also showed antiproliferative effects against cancercell lines; thus making rhodanine-based derivatives suitable leads for antiviral and anticancer drug development.
Table 7
Predicted physicochemical and ADMET properties of rhodanine-based compounds.
Compound
MWt
MolVol
MlogP
S+logP
S+logD
S+Vd
S+BBB
PSA_2D
PolASA_3D
HBD
HBA
6
451.94
322
3.58
3.68
2.65
0.25
High
64
59
2
5
7
468.39
336
3.70
3.88
2.85
0.22
High
64
59
2
5
8
478.94
333
3.29
3.03
2.26
0.21
Low
106
130
2
7
9
463.97
347
2.92
3.53
2.46
0.27
High
73
55
2
6
10
438.54
322
2.25
2.93
2.13
0.26
Low
79
66
3
5
11
441.54
336
2.58
2.87
1.86
0.29
High
81
73
2
6
12
452.45
298
2.44
2.50
1.62
0.23
Low
119
127
2
8
13
468.90
312
2.56
2.67
1.93
0.23
Low
119
126
2
8
14
479.46
308
2.24
2.47
1.50
0.19
Low
162
236
2
10
15
502.46
326
2.91
2.90
1.95
0.25
Low
119
140
2
8
16
570.46
357
3.75
3.61
2.56
0.38
Low
119
140
2
8
17
464.49
322
1.83
2.37
1.43
0.18
Low
128
136
2
9
18
494.51
350
1.61
2.41
1.38
0.18
Low
137
137
2
10
19
524.54
378
1.12
2.33
1.25
0.18
Low
146
138
2
11
20
590.60
420
2.08
2.20
1.11
0.25
Low
176
174
2
13
21
516.52
343
2.54
1.97
1.02
0.21
Low
158
234
2
11
22
579.66
469
0.36
2.76
1.35
0.21
Low
117
58
2
11
23
645.72
511
1.34
2.82
1.19
0.19
Low
147
93
2
13
24
571.64
434
2.13
2.74
1.28
0.19
Low
129
126
2
11
25
557.58
417
2.45
3.73
2.24
0.29
Low
90
50
2
8
26
625.58
448
2.96
4.41
2.58
0.39
Low
90
48
2
8
27
540.02
410
1.39
3.41
1.86
0.21
Low
111
78
3
9
28
584.47
417
1.50
3.47
1.88
0.18
Low
111
94
3
9
29
550.58
406
1.10
2.99
1.14
0.19
Low
154
157
3
11
30
550.58
406
1.10
3.12
1.15
0.20
Low
154
161
3
11
31
549.63
441
0.87
2.98
1.53
0.19
Low
120
89
3
10
32
545.64
448
1.74
3.79
2.36
0.17
Low
111
93
3
9
33
576.70
490
1.29
3.41
2.31
0.23
Low
115
96
3
10
34
609.69
469
1.51
4.04
2.06
0.25
Low
134
90
3
11
35
549.59
413
0.27
3.10
0.60
0.20
Low
149
165
4
11
36
574.47
427
1.60
3.83
2.14
0.20
Low
111
73
3
9
37
663.37
441
1.81
3.88
1.96
0.19
Low
111
87
3
9
38
614.50
445
0.98
3.13
1.59
0.18
Low
120
92
3
10
39
757.37
455
2.02
4.17
2.16
0.20
Low
111
79
3
9
40
617.79
560
2.82
5.04
3.82
0.34
Low
111
73
3
9
41
465.94
333
2.21
2.97
1.70
0.23
Low
105
108
4
7
42
510.39
340
2.32
3.04
1.75
0.22
Low
105
107
4
7
43
476.50
329
1.85
2.84
1.22
0.18
Low
148
195
4
9
44
476.50
329
1.85
2.70
1.05
0.18
Low
148
206
4
9
45
475.55
364
1.68
2.89
1.69
0.21
Low
114
124
4
8
46
471.56
371
2.32
3.20
2.05
0.18
Low
105
122
4
7
47
502.62
413
2.13
3.23
2.33
0.26
Low
109
131
4
8
48
535.61
392
2.35
3.90
2.19
0.27
Low
128
121
4
9
49
475.51
336
1.05
2.62
0.49
0.22
Low
143
195
5
9
50
500.39
350
2.43
3.40
1.97
0.21
Low
105
102
4
7
51
589.29
364
2.66
3.46
1.89
0.20
Low
105
118
4
7
52
540.42
368
1.80
3.11
1.80
0.19
Low
114
127
4
8
53
683.29
378
2.87
3.67
1.98
0.22
Low
105
111
4
7
54
543.71
483
3.44
4.53
3.51
0.39
Low
105
104
4
7
MWt: Molecular weight; MolVol: Molal volume at normal boiling point; MlogP: Moriguchi octanol-water partition coefficient; S + logP: Simulation plus octanol-water partition coefficient; S + logD: Simulation plus octanol-water distribution coefficient; S + Vd: Simulation plus volume of distribution; S + BBB: Simulation plus BBB permeability; PSA_2D: 2D-polar surface area; PolASA_3D: 3D-polar solvent accessible surface area; HBD: Hydrogen bond donors; HBA: Hydrogen bond acceptors.
Authors: Thomas A Halgren; Robert B Murphy; Richard A Friesner; Hege S Beard; Leah L Frye; W Thomas Pollard; Jay L Banks Journal: J Med Chem Date: 2004-03-25 Impact factor: 7.446
Authors: Richard A Friesner; Jay L Banks; Robert B Murphy; Thomas A Halgren; Jasna J Klicic; Daniel T Mainz; Matthew P Repasky; Eric H Knoll; Mee Shelley; Jason K Perry; David E Shaw; Perry Francis; Peter S Shenkin Journal: J Med Chem Date: 2004-03-25 Impact factor: 7.446
Authors: Meihua Luo; Sarah Delaplane; Aihua Jiang; April Reed; Ying He; Melissa Fishel; Rodney L Nyland; Richard F Borch; Xiaoxi Qiao; Millie M Georgiadis; Mark R Kelley Journal: Antioxid Redox Signal Date: 2008-11 Impact factor: 8.401
Authors: Assem Barakat; Matti Haukka; Saied M Soliman; M Ali; Abdullah Mohammed Al-Majid; Ayman El-Faham; Luis R Domingo Journal: Molecules Date: 2021-11-30 Impact factor: 4.411