scieee AI-readable full text Open interactive document viewer

5-Аrylidene derivatives of 3-(benzo[d]thiazol-2-ylamino)-2-thioxothiazolidin-4-one as potential antineoplastic agents – in silico evaluation

Mosula, L.M.; Zhmura, I.O.; Mosula, V.S.

Abstract

Computer modeling and in silico prediction of molecular properties are rational approaches in early Drug Development, facilitating the targeted synthesis. One of the promising scaffold for the Drug Design is the rhodanine cycle due to the possibility of introducing substituents in positions 3 and 5. The aim of the work is to study the spectrum of biological activity of a series of 5-arylidene derivatives of 3-(benzo[d]thiazol-2-ylamino)-2-thioxothiazolidin-4-one, predict the affinity for therapeutic targets and establish their possible affiliation to certain codes according to the Anatomical Therapeutic Chemical Classification System of drugs, and predict possible toxicity of hit compounds via the online services. The investigated series of compounds consists of 21 derivatives. In silico prediction was carried out with the use of the web resources SuperPred 3.0, ProTox 3.0. A wide spectrum of biological activity of the derivatives of the series with a predominant antitumour effect is demonstrated. The group structural similarity of compounds to antineoplastic and immunomodulating agents was established. A high probability of group-level efficacy is predicted for the antitumor therapy targets: Krüppel-like factor 5, Nuclear factor erythroid 2-related factor 2, and DNA-(apurinic/apyrimidinic site) lyase. Compounds 3 ((ethyl (E)-4-(2-(2-((3-(benzo[d]thiazol-2-ylamino)-4-oxo-2-thioxothiazolidin-5-ylidene)methyl)-4-chlorophenoxy)acetamido)benzoate)), 4 ((E)-2-(2-((3-(benzo[d]thiazol-2-yla­mino)-4-oxo-2-thioxothiazolidin-5-ylidene)methyl)-4-chlorophenoxy)-N-(4-sulfamoylphenyl)acetamide) and 9 ((E)-2-(2-((3-(benzo[d]thiazol-2-ylamino)-4-oxo-2-thioxothiazolidin-5-ylidene)methyl)-4,6-dimethylphenoxy)acetamide) can be conside­red as hit compounds. With high predictive model accuracy, their potential impact on the specified targets is 99.43%, 95.88%, and 96.92%, respectively. Compound 3 is a potential multi-hitter. The obtained results confirmed the viability of further in vitro, in vivo research.

Full text

211 25/Том XXX/3 UDC 615.277.3.015.4:547.789.13:004.94 https://doi.org/10.26641/2307-0404.2025.3.340765 L.M. Mosula * , I.O. Zhmura, V.S. Mosula 5-АRYLIDENE DERIVATIVES OF 3-(BENZO[d]THIAZOL-2-YLAMINO)-2THIOXOTHIAZOLIDIN-4-ONE AS POTENTIAL ANTINEOPLASTIC AGENTS – IN SILICO EVALUATION I. Horbachevsky Ternopil National Medical University of Ministry of Health of Ukraine Maidan Voli, 1, Ternopil, 46001, Ukraine Тернопільський національний медичний університет ім. І. Я. Горбачевського МОЗ України Майдан Волі, 1, Тернопіль, 46001, Україна *e-mail: [email protected] Цитування: Медичні перспективи. 2025. Т. 30, № 3. С. 211-219 Cited: Medicni perspektivi. 2025;30(3):211-219 Key words: rhodanine, benzothiazole, 5-arylidene substituents, biological activity, drug development, in silico prediction, antineoplastic agents, SAR analysis Ключові слова: роданін, бензотіазол, 5-ариліденпохідні, біологічна активність, розробка ліків, in silico прогноз, антинеопластичні засоби, SAR-аналіз Abstract. 5-Arylidene derivatives of 3-(benzo[ d ]thiazol-2-ylamino)-2-thioxothiazolidin-4-one as potential antineoplastic agents – in silico evaluation. Mosula L.M., Zhmura I.O., Mosula V.S. Computer modeling and in silico prediction of molecular properties are rational approaches in early Drug Development, facilitating the targeted synthesis. One of the promising scaffold for the Drug Design is the rhodanine cycle due to the possibility of introducing substituents in positions 3 and 5. The aim of the work is to study the spectrum of biological activity of a series of 5-arylidene derivatives of 3-(benzo[d]thiazol-2-ylamino)-2-thioxothiazolidin-4-one, predict the affinity for therapeutic targets and establish their possible affiliation to certain codes according to the Anatomical Therapeutic Chemical Classification System of drugs, and predict possible toxicity of hit compounds via the online services. The investigated series of compounds consists of 21 derivatives. In silico prediction was carried out with the use of the web resources SuperPred 3.0, ProTox 3.0. A wide spectrum of biological activity of the derivatives of the series with a predominant antitumour effect is demonstrated. The group structural similarity of compounds to antineoplastic and immunomodulating agents was established. A high probability of group-level efficacy is predicted for the antitumor therapy targets: Krüppel-like factor 5, Nuclear factor erythroid 2-related factor 2, and DNA-(apurinic/apyrimidinic site) lyase. Compounds 3 ((ethyl (E)-4-(2-(2-((3-(benzo[d]thiazol-2-ylamino)-4oxo-2-thioxothiazolidin-5-ylidene)methyl)-4-chlorophenoxy)acetamido)benzoate)), 4 ((E)-2-(2-((3-(benzo[d]thiazol-2-ylamino)-4-oxo-2-thioxothiazolidin-5-ylidene)methyl)-4-chlorophenoxy)-N-(4-sulfamoylphenyl)acetamide) and 9 ((E)-2-(2-((3- (benzo[d]thiazol-2-ylamino)-4-oxo-2-thioxothiazolidin-5-ylidene)methyl)-4,6-dimethylphenoxy)acetamide) can be considered as hit compounds. With high predictive model accuracy, their potential impact on the specified targets is 99.43%, 95.88%, and 96.92%, respectively. Compound 3 is a potential multi-hitter. The obtained results confirmed the viability of further in vitro, in vivo research. Реферат. 5-Ариліденпохідні 3-(бензо[ d ]тіазол-2-іламіно)-2-тіоксотіазолідин-4-ону як потенційні протипухлинні агенти - оцінка in silico. Мосула Л.М., Жмура І.О, Мосула В.С. На початкових етапах створення ліків раціональним є використання комп’ютерного моделювання та in silico прогнозування властивостей молекул, що сприяє цілеспрямованому синтезу. Одним з перспективних структурних «каркасів» для дизайну ліків є роданіновий цикл завдяки можливості введення замісників у 3 і 5 положення. Мета роботи – дослідити спектр біологічної активності ряду 5-ариліденпохідних 3-(бензо[d]тіазол-2-іламіно)-2-тіоксотіазолідин-4-ону, спрогнозувати афінність до терапевтичних мішеней та встановити їх можливу приналежність до певних кодів згідно з анатомо-терапевтичною хімічною класифікацією лікарських засобів, а також передбачити можливу токсичність «сполук-хітів» за допомогою онлайн-сервісів. Досліджуваний ряд сполук містить 21 похідне. In silico прогнозування було проведено за допомогою веб-ресурсів SuperPred 3.0, ProTox 3.0. Продемонстровано потенційно широкий спектр біологічної активності серії похідних з переважаючим протипухлинним ефектом. Установлено групову структурну подібність серії похідних до антинеопластичних та імуномодулюючих лікарських засобів. З високою ймовірністю передбачається групова ефективність стосовно мішеней протухлинної терапії: Krüppel-like factor 5, Nuclear factor erythroid 2-related factor 2 і DNA-(apurinic or apyrimidinic site) lyase. Сполуки 3 (етил (E)-4-(2-(2-((3-(бензо[d]тіазол-2-іламіно)-4-оксо-2-тіоксотіазолідин-5іліден)метил)-4-хлорофеноксі)ацетамідо)бензоат), 4 ((E)-2-(2-((3-(бензо[d]тіазол-2-іламіно)-4-оксо-2- ФАРМАЦІЯ 212 Н а умовах ліцензі ї C C BY 4.0 тіоксотіазолідин-5-іліден)метил)-4-хлорофенокси)-N-(4-сульфамоїлфеніл)ацетамід) і 9 ((Е)-2-(2-((3-(бензо[d]тіазол-2-іламіно)-4-оксо-2-тіоксотіазолідин-5-іліден)метил)-4,6-диметилфеноксі)ацетамід) можна розглядати як «сполуки-хіти». При високій точності моделей прогнозування їх потенційний вплив на зазначені мішені становить 99,43%, 95,88%, 96,92%, відповідно. Сполука 3 є потенційним мультихітером. Одержані результати підтвердили доцільність подальших in vitro, in vivo досліджень. Increasing the efficiency of the search for new biologically active compounds (BACs), decreasing the terms of their implementation is only possible if effective pre-experimental studies are carried out. In recent years, researchers have been focused on the use of new approaches to Drug Development, such as exploring the possibilities of the use of targeted synthesis of new substances with potential pharmacological properties on the basis of the modern in silico technologies based on the relationship between chemical structure and biological activity. The “structure – activity”relationship (SAR analysis) is a great method for Drug Development, from the evaluation of drug targets to the improvement of the properties of a molecule [1]. Rhodanine derivatives with a benzothiazole moiety in the molecules are convenient objects of study. They are characterised by antitumour, antiviral and anti-tuberculosis activity [2].This leads to the logical conclusion of the viability of their use for the targeted design of highly active substances. The aim of the work is to study the spectrum of biological activity of a series of 5-arylidene derivatives of 3-(benzo[d]thiazol-2-ylamino)-2-thioxothiazolidin-4-one, predict the affinity for therapeutic targets and establish their possible affiliation to certain codes according to the Anatomical Therapeutic Chemical Classification System (ATC) of drugs, and predict possible toxicity of hit compounds via the online services. MATERIALS AND METHODS OF RESEARCH For the in silico study we used some earlier synthesized (compounds 1-6) with previously confirmed in vitro antitumour activity [2] and virtually modelled (compounds 7-21) derivatives on the basis of 3-(benzo[d]thiazol-2-ylamino)-2-thioxothiazolidin-4-one (Fig. 1). Fig. 1. General structural formula of differentially substituted 5-arylidene derivatives of 3-(benzo[d]thiazol-2-ylamino)-2-thioxothiazolidin-4-one Previous in silico studies of the physicochemical and pharmacokinetic properties of the molecules confirmed the prospects of searching for BACs among rhodanine derivatives with a benzothiazole fragment in the molecules. Prediction of ADME properties of compounds in a series, performed on the МЕДИЧНІ ПЕРСПЕКТИВИ / MEDICNI PERSPEKTIVI 213 25/Том XXX/3 Swiss ADME online platform, described by us earlier [3]. In this study, we describe the impact of variously substituted arylidene fragments introduced at position 5 of the rhodanine ring on the biological activity of the compounds, and present the results of in silico toxicity prediction, which represents an important component of the Drug Development process. The antitumour activity of some derivatives (compounds 1-6) has been proven in vitro and described earlier. These derivatives exhibited moderate antitumour activity with high selectivity of action against certain lines of cancer cells. Among them, a highly active compound 2((E)-2-(2-((3-(benzo[d]thiazol-2ylamino)-4-oxo-2-thioxothiazolidin-5-ylidene)methyl)-4-chlorophenoxy)-N-(4-methoxyphenyl)acetamide), was identified which was found to be nontoxic (LD50>1000 mg/kg) in vivo [2]. Chemistry The chemical composition and structure of the core heterocycle and compounds 1-6 were determined with the use of elemental analysis and 1Н NMR-spectroscopy [2]. This paper discusses the potential of the new hybrid structures (compounds 721), that were virtually modeled on the basis of 3- (benzo[d]thiazol-2-ylamino)-2-thioxothiazolidin-4one, and the earlier synthesized compounds 1-6. In silico experiment For our set goal we used the capabilities of the updated version of the web tool SuperPred 3.0 [4], that allows us with high reliability to compare the inputted molecule to all available targets at once and to predict their possible ATC codes with accuracy of up to 80.5%. According to the definition provided by SuperPred 3.0, a compound is considered a potential multi-hitter if it shows a high probability (>1%) of belonging to specific ATC classes and a very high possibility of binding to a biotarget (>80%). These thresholds may appear simple, yet they are underpinned by complex logic involving machine learning, statistical validity, and practical interpretation within the pharmaceutical context. They are designed to identify compounds capable of interacting with multiple targets and exhibiting therapeutic activity across various pharmacological classes. Small parts of the molecule are enough to detect structural similarity to finished pharmaceutical products (FPPs) [5].Using this web-platform it is possible to calculate quantitative parameters of the probable binding of ligands to targets (in %) and the overall accuracy of the prediction model (in %) with a high accuracy [6]. Additionally, to reduce the number of animal experiments in future research, we performed computational toxicity estimations for the main hit compounds using the online software ProTox 3.0 [7]. The prediction results were obtained with the use of in silico tools SuperPred 3.0 and ProTox 3.0, therefore they do not entail the use of statistical methods and research involving humans or animals. Independent bioethics committee approval and written consent were not obtained. RESULTS AND DISCUSSION We studied the structural similarity of 5-arylidene derivatives to the known FPPs and their ability to bind with biomolecular targets and their probable affiliation to certain ATC classes. According to the prognosis of SuperPred 3.0, all the derivatives of the series have more than one indication, so they belong to several ATC codes. There is a group structural similarity of the compounds to antineoplastic and immunomodulating agents, which are assigned code (L01XE – protein kinase inhibitors). The highest possibility of a structural similarity is predicted for compound 5 which is 35.80%. Additionally, for compounds 8, 12, 14, 16–21 a similarity to drugs with a code L04АА (selective immunosuppressants) is predicted. Compound 20 is the most structurally similar to them (28.60%). SuperPred 3.0 found structural similarity of the studied compounds to not just antitumour MPs, but also other drugs, however a group similarity to other FPPs with a different pharmacological activity is not predicted. Compounds 1, 3, 7 and 9, as well as the core heterocycle, are most structurally similar to antifungals for topical use (D01AE). For those derivatives the predicted values range from 92.52% to 15.17%, notably the highest value is predicted for compound 3. For other derivatives (compounds 1319 and 21) structural similarity to antithrombotic agents (platelet aggregation inhibitors, except for heparin B01AC) is predicted with a probability in the range of 39.12-8.74%, and compound 2 is most similar to FPPs from the group of direct factors Xa inhibitors (B01AF). With a probability of 16.19% and 17.38% compounds 10 and 11 are predicted to belong to the group of hypnotics and sedatives (N05CM). The highest value of structural similarity to antiglaucoma preparations and miotics (S01EC) is predicted for compound 4 (21.15%), and for compound 12, the predicted similarity to cardiovascular system drugs (simple angiotensin II receptor blockers plain C09CA) is 18.69%. Among the predicted targets, those associated with antitumor therapy predominated. We considered targets common to all compounds in the series that demonstrated the highest binding probabilities and high predictive model accuracy (>85%) (Fig. 2). Fig. 2. The highest probability of binding of the core heterocycle and compounds 1–21 with targets for cancer therapy (in %) 50.00% 55.00% 60.00% 65.00% 70.00% 75.00% 80.00% 85.00% 90.00% 95.00% 100.00% Core heterocycle Compound 1 Compound 2 Compound 3 Compound 4 Compound 5 Compound 6 Compound 7 Compound 8 Compound 9 Compound 10 Compound 11 Compound 12 Compound 13 Compound 14 Compound 15 Compound 16 Compound 17 Compound 18 Compound 19 Compound 20 Compound 21 90.90% 98.79% 98.74% 99.43% 97.64% 98.02% 98.32% 97.83% 98.75% 98.34% 95.79% 98.65% 99.01% 98.10% 99.11% 98.25% 99.11% 98.56% 98.76% 98.54% 98.76% 97.17% 72.74% 92.57% 95.10% 93.49% 95.88% 94.24% 94.65% 90.90% 91.07% 91.18% 94.80% 91.33% 91.07% 94.15% 91.87% 93.39% 92.34% 93.77% 91.75% 94.19% 91.75% 89.35% 89.89% 91.11% 94.11% 54.65% 83.60% 90.13% 94.39% 90.77% 88.63% 96.92% 71.61% 95.33% 89.47% 96.41% 90.83% 90.22% 92.08% 93.64% 94.44% 93.38% 95.29% 89.77% Krüppel-like factor 5 Nuclear factor erythroid 2-related factor 2 DNA-(apurinic or apyrimidinic site) lyase Probability of compound-target binding МЕДИЧНІ ПЕРСПЕКТИВИ / MEDICNI PERSPEKTIVI 215 25/Том XXX/3 The highest level of binding of the compounds was determined for three targets: Krüppel-like factor 5 (KLF5), Nuclear factor erythroid 2-related factor 2 (NRF2) and DNA-(apurinic or apyrimidinic site) lyase (DNA AP lyase). Compound 3 exhibited the highest predicted interaction probability with KLF5, compound 4 – with NRF2, and compound 9 – with DNA AP lyase, which supports their classification as hit compounds, i.e., promising candidates for further lead optimization. For these compounds, we performed in silico oral toxicity prediction using the ProTox 3.0 web server. The prediction results are presented in Table 1. Table 1 Predicted oral toxicity of hit compounds Compound 3 Compound 4 Compound 9 Hepatotoxicity + (62%) + (54%) - Carcinogenicity - - + (51%) Immunotoxicity - - - Mutagenicity - - - Cytotoxicity - - - LD50 (mg/kg) 1400 2000 1400 Class IV IV IV Compounds classified under toxicity class IV (harmful if swallowed; 300<LD₅₀≤2000 mg/kg) are not considered critically toxic, yet they require careful attention during Drug Development. Only the probability of binding compounds to therapeutic targets and model accuracy that exceeded 80% was taken into account. As evident from the results of the virtual prediction, all the studied compounds are likely to exhibit a wide range of biological activity, and a high probability of their interaction with therapeutic targets is predicted with a marked prediction accuracy. Based on the list of predicted targets, 5-arylidene derivatives are similar to the core heterocycle, and compound 3 (R1 = 1-OCH2CONHC6H4-п-COOC2H5; R2 = 4-Cl) is a potential multi-hitter with a predominant structural similarity to drugs intended to treat dermatological diseases (D01AE). The in silico prediction confirmed the in vitro antitumor activity of compound 2, with binding probabilities to the specified targets of 98.74%, 95.10%, and 94.11%, respectively. Compounds 3, 4, and 9 were also identified as potential hits (Fig. 3). Compound 3 ethyl (Е)-4-(2-(2-((3-(benzo[d]thiazol-2-ylamino)-4-oxo-2-thioxothiazolidin-5-ylidene)methyl)- 4-chloro-phenoxy)acetamido)benzoate (99.43% to KLF5) ФАРМАЦІЯ 216 Н а умовах ліцензі ї C C BY 4.0 Compound 4 (Е)-2-(2-((3-(benzo[d]thiazol-2-ylamino)-4-oxo-2-thioxothiazolidin-5-ylidene)methyl)-4-chlorophenoxy)-N-(4sulfamoylphenyl)acetamide (95.88% to NRF2) Compound 9 (Е)-2-(2-((3-(benzo[d]thiazol-2-ylamino)-4-oxo-2-thioxothiazolidin-5-ylidene)methyl)-4,6dimethylphenoxy)acetamide (96.92% to DNA AP lyase) Fig. 3. Chemical structures of the compounds with the highest predicted affinity values to anti-tumour targets KLF5 is a zinc-finger-containing transcription factor that has a significant role in the regulation of the expression of various genes, which impacts diverse cellular functions [8]. KLF5 is highly expressed in various types of cancer. KLF5 knockdown reduces the growth of cancer cells. KLF5 is involved in a variety of cellular functions, such as differentiation, proliferation, migration, apoptosis and regulation of cancer stem cells. The expression and activity of KLF5 is distorted in many types of cancer, including ovarian, cervical, breast, colorectal, lung, pancreatic and thyroid cancers [9]. NRF2, also known as nuclear factor erythroidderived 2-like 2, is a transcription factor encoded by the NFE2L2 gene in humans. It has a critical role in the regulation of the expression of genes involved in МЕДИЧНІ ПЕРСПЕКТИВИ / MEDICNI PERSPEKTIVI 217 25/Том XXX/3 the oxidative stress response and the maintenance of redox homeostasis [10]. It was found that in the context of cancer, NRF2 plays a dual role. On one hand, it might protect normal cells from oxidative damage and reduce the risk of cancer through the suppression of oxidative stress and inflammation, which contribute to tumour development. On the other hand, many types of cancer show an increased activity of NRF2, which may lead to a poor prognosis and therapeutic resistance [11]. The activation of NRF2 is linked to the upregulation of the genes that facilitate cell survival and proliferation, which makes it a key player in the development of resistance to anticancer therapy. Targeting NRF2 and its pathways downstream has become an area of interest in the development of new therapies for cancer [10]. DNA AP lyase is an enzyme that has a critical role for repairing DNA damage. It is specifically targeted at apurinic/apyrimidinic (AP) sites, which are places in DNA that lack a purine or pyrimidine base. These sites can occur spontaneously or as a result of DNA damage caused by a variety of factors, such as radiation, chemicals or enzymatic activity. Deficiencies in these pathways often lead to mutagenesis and instability of the genome, which is known to be linked to human diseases such as cancer and neurodegenerative disorders [12]. In the context of antitumour therapy, DNA AP lyase is important as it helps maintain genome integrity, repairing DNA damage. A targeted action on DNA repair enzymes, in particular on DNA AP lyase, can increase the effectiveness of antitumour therapy by preventing cancer cells from repairing the damage caused by chemotherapy and radiotherapy [13]. The latest advances in antitumour therapy are focused on various strategies, such as immune checkpoint inhibitors, therapeutic vaccines and adoptive Tcell therapy. These approaches aim to harness the body's immune system to more effectively detect and destroy cancer cells. Using DNA repair inhibitors in combination with these immunotherapies, researchers are hoping for improved treatment outcomes and reduced chances of recurrence of cancer [14, 15, 16]. We performed a comparative analysis of the predicted parameters for the identified hit compounds and known inhibitors of KLF5, NRF2, and DNA AP. According to SuperPred 3.0, the binding probability of established inhibitors KLF5 (ML264, Tolfenamic acid, Mithramycin A, JQ1, Kenpaullone), NRF2 (ML385, Brusatol, Luteolin, Trigonelline, Halofuginone, Keap1), and DNA AP (CRT0044876, Lucanthone, E3330 (APX3330), Methoxyamine, RN8) ranges from 70% to 95%, depending on the specific compound structure and the predictive model’s accuracy. For ML264, one of the most well-known KLF5 inhibitors, the predicted binding probability does not exceed 90%, which is consistent with its experimentally confirmed activity [17]. The predicted binding probability for compound 3 is higher (approximately 100%) at a model accuracy of 86.33%, indicating the reliability of the prediction. This supports the consideration of compound 3 as a potential candidate for further docking analyses and biological testing. For example, compound 4 demonstrates a higher predicted probability of binding to NRF2 (95.88%) at high model accuracy (96%) compared to ML385 (not exceeding 90%), one of the most well-known selective NRF2 inhibitors that acts through binding to the DLG motif of KEAP1 [18]. This indicates high selectivity and the potential of the novel compound as an NRF2 inhibitor. The in silico predicted interaction probabilities for compound 9 and the known non-specific inhibitor CRT0044876 with DNA AP lyase are comparably high (above 90%). The inhibitory activity of CRT0044876 against DNA AP lyase has been experimentally confirmed and described previously [19]. Subsequent studies demonstrated that CRT0044876 binds to a pocket distant from the active site of APE1 and forms colloidal aggregates, which may lead to non-specific inhibition [20]. The high predicted binding probability of compound 9 to the target, combined with the strong model accuracy (91.11%), indicates its considerable potential. The obtained results justify further investigation of these compounds as promising candidates for antitumor therapy. Aside from antitumour activity, the derivatives of the series are predicted to have a high probability of antidiabetic, antifungal, antiarteriosclerotic, antihypertensive and other kinds of activities. Analysis of the “chemical structure – biological activity” relationship in the series of derivatives allows us to state that arylidene fragment at position 5 has a significant influence on the expression of antitumour activity. The increase in binding probability of the compounds to their respective targets depends on the introduced substituents (R1, R2 and R3). The size of the molecule and branching also play a significant role. The most optimal in this context is the presence of fragments -C6H4-п-COOC2H5 (compound 3), -C6H4-п-SO2NH2(compound 4), which are connected to phenoxyacetamide moiety in molecules. Compound 9 is characterized by high branching: R1=1-OCH2CONH2, R2=4-CH3, R3=6-CH3 (Fig. 3). The introduction of arylidene fragments facilitated the increase in the structural similarity of the derivatives to antineoplastic agents and positively influenced their potential affinity for their aforementioned ФАРМАЦІЯ 218 Н а умовах ліцензі ї C C BY 4.0 targets. This can be clearly traced while analyzing the values obtained for 5-arylidene derivatives and the unsubstituted core heterocycle. The predicted toxicity class IV of the hit compounds, according to the ProTox 3.0 classification, indicates moderate toxicity and does not warrant their exclusion, particularly in light of their anticipated bioactivity. Thorough optimization of the dosage and administration route is advisable, as the predicted toxicity is associated with oral administration. Several antitumor agents (e.g., Methotrexate, Docetaxel, Vincristine) exhibit similar toxicological profiles yet are successfully used in clinical practice due to their proven efficacy and controlled dosing, which enables safe therapeutic application. CONCLUSIONS 1. The results of in silico prediction indicate a potentially broad therapeutic spectrum and predominantly antitumor activity of the 5-arylidene derivatives of 3-(benzo[d]thiazol-2-ylamino)-2-thioxothiazolidin-4-one. Structural group similarity to known antineoplastic and immunomodulatory agents was observed, and compound 3 is predicted to possess multifunctional properties. 2. Several hit compounds were identified within the series. Compounds 3, 4, and 9 demonstrated higher predicted binding affinities to key targets (Krüppel-like factor 5, Nuclear factor erythroid 2related factor 2, DNA-(apurinic or apyrimidinic site) lyase) compared to compound 2, whose antitumor activity has been previously confirmed in vitro. 3. Despite the predicted toxicity class IV, the identified candidates for further optimization exhibit a favorable balance between potential efficacy and acceptable toxicological profiles (absence of mutagenicity or/and carcinogenicity). The presence of shared pharmacophoric fragments between the hit compounds and known inhibitors of the respective targets supports the rationale for in-depth investigation. 4. The obtained results require further experimental validation. Acknowledgements The authors gratefully acknowledge the valuable knowledge and experience acquired during the lecture part of the II European Chemistry School for Ukrainians (website: https://acmin.agh.edu.pl/en/detail/s/iieuropean-chemistry-school-for-ukrainians). Contributors: Mosula L.M. – conceptualization, data curation, methodology, formal analysis, investigation, writing – original draft, writing – review and editing the manuscript, manuscript translation into English, supervision, project administration; Zhmura I.O. – formal analysis, investigation, resources search, writing – original draft, visualization. Mosula V.S. – formal analysis, investigation, resources search, writing – original draft. Funding. This research received no external funding. Conflict of interests. The authors declare no conflict of interest. REFERENCES 1. Akash S, Arefin F, Aovi FI. In silico investigation of potential therapeutic medication for the inhibition of dengue virus (DENV NS2B/NS3 and NS1) by modification of polycyclic quaternary alkaloid (sanguinarine derivatives) with different computational approaches. Biointerface Res Appl Chem. 2023;13(5):403. doi: https://doi.org/10.33263/BRIAC135.403 2. Mosula L, Zimenkovsky B, Havrylyuk D, Missir A-V, Chirita IC, Lesyk R. Synthesis and antitumor activity of novel 2-thioxo-4-thiazolidinones with benzothiazole moieties. [Internet]. Farmacia. 2009 [cited 2025 Jun 29];57(3):321-30. Available from: https://farmaciajournal.com/arhiva/20093/issue32009art08.pdf 3. Mosula LM, Kolodych IO, Mosula VS. [Virtual screening of the properties of a series of 5-arylidene derivatives of 3-(benzo[d]thiazol-2-ylamino)-2-thioxothiazolidin-4-one]. Annaly Mechnykovskoho Instytutu. 2025;1:73-80. Ukrainian. doi: https://doi.org/10.5281/zenodo.15011532 4. SuperPred 3.0 [Internet]. [cited 2025 Jun 29]. Available from: https://prediction.charite.de/index.php/ 5. Kopak NA. Searching of biological activity of Sesters 4acetylaminobenzenethiosulfoacid using methods of chemoinformatics. Chem Technol Appl Subst. 2023;6(2):76-86. doi: https://doi.org/10.23939/ctas2023.02.076 6. Gallo K, Goede A, Preissner R, Gohlke BO. SuperPred 3.0: drug classification and target prediction-a machine learning approach. Nucleic Acids Res. 2022;50(W1):W726-W31. doi: https://doi.org/10.1093/nar/gkac297 7. ProTox 3.0 [Internet]. [cited 2025 Jun 29]. Available from: https://tox.charite.de 8. Gao Y, Ding Y, Haiying C, Haijun C, Zhou J. Targeting Krüppel-Like Factor 5 (KLF5) for Cancer Therapy. Curr Top Med Chem. 2015;15(8):699-13. doi: https://doi.org/10.2174/1568026615666150302105052 9. Siraj AK, Pratheeshkumar P, Divya SP, et al. Krupple-Like Factor 5 is a Potential Therapeutic Target and Prognostic Marker in Epithelial Ovarian Cancer. Front Pharmacol. 2020;11:598880. doi: https://doi.org/10.3389/fphar.2020.598880 МЕДИЧНІ ПЕРСПЕКТИВИ / MEDICNI PERSPEKTIVI 219 25/Том XXX/3 10. Paramasivan P, Kankia IH, Langdon SP, Deeni YY. Emerging role of nuclear factor erythroid 2-related factor 2 in the mechanism of action and resistance to anticancer therapies. Cancer Drug Resist. 2019;2:490-15. doi: http://dx.doi.org/10.20517/cdr.2019.57 11. He F, Antonucci L, Karin M. NRF2 as a regulator of cell metabolism and inflammation in cancer. Carcinogenesis. 2020;41(4):405-16. doi: https://doi.org/10.1093/carcin/bgaa039 12. Lin Y, McMahon A, Driscoll G, Bullock S, Zhao J, Yan S. Function and molecular mechanisms of APE2 in genome and epigenome integrity. Mutat Res Rev Mutat Res. 2021;787:108347. doi: https://doi.org/10.1016/j.mrrev.2020.108347 13. Hans F, Senarisoy M, Bhaskar Naidu C, Timmins J. Focus on DNA Glycosylases-A Set of Tightly Regulated Enzymes with a High Potential as Anticancer Drug Targets. Int J Mol Sci. 2020;21(23):9226. doi: https://doi.org/10.3390/ijms21239226 14. Zhang A, Fan T, Liu Y, Yu G, Li C, Jiang Z. Regulatory T cells in immune checkpoint blockade antitumor therapy. Mol Cancer. 2024;23(1):251. doi: https://doi.org/10.1186/s12943-024-02156-y 15. Dembic Z. Antitumor Drugs and Their Targets. Molecules. 2020;25(23):5776. doi: https://doi.org/10.3390/molecules25235776 16. Albarran V, San Roman M, Pozas J, et al. Adoptive T cell therapy for solid tumors: current landscape and future challenges. Front Immunol. 2024;15:1352805. doi: https://doi.org/10.3389/fimmu.2024.1352805 17. Shen X, Zhang Y, Xu Z, et al. KLF5 inhibition overcomes oxaliplatin resistance in patient-derived colorectal cancer organoids by restoring apoptotic response. Cell Death Dis. 2022;13:303. doi: https://doi.org/10.1038/s41419-022-04773-1 18. Singh A, Venkannagari S, Oh KH, et al. Small Molecule Inhibitor of NRF2 Selectively Intervenes Therapeutic Resistance in KEAP1-Deficient NSCLC Tumors. ACS Chemical Biology. 2016;11(11):3214-25. doi: https://doi.org/10.1021/acschembio.6b00651 19. Madhusudan S, Smart F, Shrimpton P, et al. Isolation of a small molecule inhibitor of DNA base excision repair. Nucleic Acids Research. 2005;33(15):4711-24. doi: https://doi.org/10.1093/nar/gki781 20. Pidugu LS, Servius HW, Sevdalis SE, et al. Characterizing inhibitors of human AP endonuclease 1. PloS ONE. 2023;18(1):e0280526. doi: https://doi.org/10.1371/journal.pone.0280526 Стаття надійшла до редакції 04.08.2025; затверджена до публікації 01.09.2025