In Silico ADMET and Molecular Docking Studies of Thiazolidinone Derivatives as Potential Antitubercular Agents
1. Introduction
Tuberculosis (TB) is a chronic infectious disease caused primarily by Mycobacterium tuberculosis, which mainly affects the lungs but may also involve other organs of the body. According to the World Health Organization (WHO), tuberculosis remains one of the leading causes of death from infectious diseases worldwide, with millions of new cases reported annually. Despite the availability of several therapeutic agents, tuberculosis continues to be a major cause of morbidity and mortality, particularly in developing countries. The emergence of multidrug-resistant tuberculosis (MDR-TB) and extensively drug-resistant tuberculosis (XDR-TB) has become a major global health challenge, reducing the effectiveness of existing therapies and emphasizing the urgent need for the discovery of novel antitubercular agents with improved efficacy, reduced toxicity, and better pharmacokinetic properties [1-7]. Heterocyclic compounds play an essential role in medicinal chemistry because of their broad spectrum of biological activities. Among these, thiazolidinone derivatives have attracted considerable attention due to their diverse pharmacological properties, including antimicrobial, antitubercular, anti-inflammatory, anticancer, antiviral, antifungal, and antioxidant activities. The thiazolidinone nucleus is considered an important pharmacophore capable of interacting effectively with biological targets through hydrogen bonding and hydrophobic interactions.
Pyridine-containing compounds are also well known for their significant therapeutic importance. Incorporation of pyridine moieties into heterocyclic frameworks often enhances biological activity, improves solubility, and increases binding affinity toward target proteins. Therefore, the combination of thiazolidinone and pyridine scaffolds may provide promising lead compounds for the development of effective antitubercular agents [8-14]. In recent years, computational approaches such as molecular docking and ADMET prediction studies have become valuable tools in modern drug discovery and development. Molecular docking helps predict ligand–protein interactions and binding affinity, while ADMET analysis provides important information regarding absorption, distribution, metabolism, excretion, and toxicity profiles of drug candidates. These computational approaches significantly reduce the time, cost, and failure rate associated with conventional drug discovery processes and assist in identifying promising lead molecules with favourable pharmacokinetic properties.
Based on these considerations, the present study was undertaken to synthesize a series of novel N-(4-oxo-3-aryl-1,2-thiazolidin-2-yl) pyridine-3-carboxamide derivatives and evaluate their antitubercular activity against Mycobacterium tuberculosis H37Rv strain. The synthesized compounds were characterized using FT-IR and ¹H NMR spectroscopy. Furthermore, in silico drug-likeness evaluation, bioactivity prediction, ADMET studies, and molecular docking analyses were performed to investigate their pharmacological potential and binding interactions with the target protein. The study aims to identify promising thiazolidinone derivatives that could serve as potential lead molecules for future antitubercular drug development.
4. Antitubercular Activity Evaluation
The antitubercular activity of the synthesized compounds was evaluated against Mycobacterium tuberculosis H37Rv using the Resazurin Microtiter Assay (REMA) plate method, a rapid and reliable technique for determining the minimum inhibitory concentration (MIC). The assay was performed to assess the inhibitory potential of the synthesized compounds against the growth of M. tuberculosis under in vitro conditions. Briefly, 100 µL of Middlebrook 7H9 broth medium containing the test microorganism was dispensed into each well of a sterile 96-well microtiter plate. The synthesized compounds were prepared at concentrations ranging from 1.25 to 10.0 µg/mL and added to the respective wells. To minimize evaporation during incubation, the peripheral wells of the plate were filled with sterile distilled water. A standardized bacterial suspension equivalent to McFarland standard No. 1 was prepared and further diluted (1:20) in 7H9 broth. Subsequently, 100 µL of the diluted inoculum was added to each test well. Wells containing inoculated medium without the test compound served as growth controls, while wells containing only sterile medium were maintained as negative controls. The plates were sealed and incubated at 37°C for seven days. Following incubation, 25 µL of freshly prepared 0.02% resazurin solution was added to each well, and the plates were incubated for an additional 48 hours. Resazurin acts as an oxidation–reduction indicator, changing from blue to pink in the presence of metabolically active bacterial cells. The MIC was defined as the lowest concentration of the test compound that prevented the color change from blue to pink, indicating complete inhibition of bacterial growth. The antitubercular susceptibility testing was carried out at AR Lab and Research Foundation, Salem, Tamil Nadu, India, following standard laboratory protocols for M. tuberculosis H37Rv evaluation. This method provided a quantitative assessment of the antimycobacterial potential of the synthesized compounds and facilitated comparison with established antitubercular agents.
Anti-tubercular activity of synthesized compounds (3a-3e) against Mycobacterium tuberculosis H37RV
5. COMPUTATIONAL METHODOLOGY
5.1 Molinspiration Physicochemical and Bioactivity Prediction
The physicochemical properties and bioactivity scores of the synthesized compounds were evaluated using the Molinspiration Cheminformatics online server. Parameters, including molecular weight (MW), logarithm of partition coefficient (logP), topological polar surface area (TPSA), number of hydrogen bond donors and acceptors, number of rotatable bonds, and Lipinski’s Rule of Five violations, were calculated to assess drug-likeness properties. Bioactivity scores for GPCR ligand, ion channel modulator, kinase inhibitor, nuclear receptor ligand, protease inhibitor, and enzyme inhibitor activities were also predicted using Molinspiration software.
5.2 ADMET Prediction Studies
The ADMET (Absorption, Distribution, Metabolism, Excretion, and Toxicity) properties of the synthesized compounds were predicted using the pkCSM online pharmacokinetic prediction server. The canonical SMILES of the synthesized compounds were used as input for the prediction of pharmacokinetic parameters, including water solubility, intestinal absorption, Caco-2 permeability, blood–brain barrier permeability, cytochrome P450 interactions, total clearance, AMES toxicity, hepatotoxicity, hERG inhibition, and skin sensitization. The obtained results were analyzed to evaluate the drug-likeness and safety profiles of the synthesized compounds.
5.3 Molecular Docking Studies
Molecular docking studies were carried out using AutoDock Vina to evaluate the binding affinity and interaction of the synthesized compounds with the Mycobacterium tuberculosis enoyl-acyl carrier protein reductase (InhA) target protein. The three-dimensional crystal structure of the target protein was obtained from the Protein Data Bank (PDB).
Protein preparation involved the removal of water molecules and co-crystallized ligands, followed by the addition of polar hydrogen atoms and Kollman charges before docking analysis. Energy minimization of the protein structure was performed to obtain a stable conformation for docking studies. The ligand structures were drawn, optimized, and converted into appropriate PDBQT file formats using molecular modeling tools.
Grid box dimensions were defined around the active binding site of the InhA protein to ensure proper accommodation of ligands during docking simulations. Docking calculations were performed to predict the most favourable binding conformations based on minimum binding energy values. The docking interactions, hydrogen bonding, hydrophobic interactions, and binding poses of the synthesized compounds within the active site were analyzed using molecular visualization software.
Computational Techniques
MW = molecular weight;
logP= logarithm of the partition coefficient between octanoland water;
TPSA = topological polar surface area;
n atom =number of atoms present in a given molecule (This value is often used to assess the potential for hydrogen bonding, acidity/basicity, and other functional group characteristics in the molecule);
nON = number of hydrogen bond acceptors (N and O atoms);
nOHNH = the number of hydroxyl (OH) and amine (NH) groups in a given molecule;
a_don = number of hydrogen bond donors;
a_acc = number of hydrogen bond acceptors;
Nrotb = number of rotatable bonds;
nViolation = number of Lipinski’s Rule of Five violations (nViolation value of 0 would indicate no structural issues or violations, meaning the molecule follows standard bonding and stereochemical conventions compounds 3c, 3d, and 3e were predicted to exhibit hepatotoxic potential.
5. RESULTS AND DISCUSSION
5.1 Spectral Characterization
FT-IR Spectroscopy
The synthesized N-(4-oxo-3-aryl-1,2-thiazolidin-2-yl) pyridine-3-carboxamide derivatives 3(a–e) were characterized by FT-IR spectroscopy, and the obtained spectral data confirmed the successful formation of the target compounds. Characteristic absorption bands corresponding to important functional groups were observed in all synthesized derivatives.
Broad absorption bands observed in the region 3481.62–3648.90 cm⁻¹ confirmed the presence of hydroxyl (–OH) and amide (>N–H) stretching vibrations. Aromatic C–H stretching vibrations appeared in the range of 3023.48–3100.14 cm⁻¹, indicating the presence of aromatic moieties in all compounds. Strong absorption bands observed between 1682.94–1749.08 cm⁻¹ corresponded to carbonyl (C=O) stretching vibrations of the thiazolidinone and amide functionalities.
The aromatic C=C stretching vibrations were observed around 1421.19–1472.81 cm⁻¹, while aromatic C–N stretching vibrations appeared in the range of 1303.22–1395.59 cm⁻¹. Characteristic C–S stretching vibrations of the thiazolidinone ring were observed at 1161.62–1188.94 cm⁻¹, confirming cyclization and formation of the thiazolidinone nucleus. In compound 3d, the presence of a C–Cl stretching vibration at 669.12 cm⁻¹ further confirmed substitution of the chlorophenyl group. Similarly, compounds containing methoxy and furyl substituents exhibited characteristic C–O–C stretching vibrations in the range of 1244.23–1264.23 cm⁻¹.
Overall, the FT-IR spectral data strongly supported the successful synthesis of the designed thiazolidinone derivatives.
¹H NMR Spectroscopy
The ¹H NMR spectral analysis further confirmed the structures of the synthesized compounds 3(a–e). The spectra showed characteristic signals corresponding to the thiazolidinone ring, aromatic protons, and substituted functional groups.
The methylene protons (CH₂) of the thiazolidinone ring appeared as doublets around δ 4.03–4.77 ppm, confirming formation of the heterocyclic ring system. The methine proton attached to the thiazolidinone ring appeared as a singlet at δ 5.05–5.34 ppm. Aromatic protons were observed as multiplets in the range of δ 6.22–8.55 ppm, corresponding to substituted phenyl, furyl, and pyridine ring protons.
The amide NH proton appeared as a singlet around δ 8.91 ppm in all compounds, indicating the presence of the carboxamide linkage. In compound 3e, methoxy group protons appeared as a singlet at δ 3.75 ppm, confirming the presence of the methoxy substituent. Compound 3a exhibited characteristic signals attributable to the hydroxyl-substituted aromatic ring.
Thus, the obtained ¹H NMR spectral data were in good agreement with the proposed structures of the synthesized thiazolidinone derivatives.
5.2 Antitubercular Activity
The synthesized compounds 3(a–e) were evaluated for antitubercular activity against Mycobacterium tuberculosis H37Rv strain using the Resazurin Microtiter Assay (REMA) plate method. The minimum inhibitory concentration (MIC) values indicated that all synthesized compounds exhibited appreciable antitubercular activity against the tested organism.
Among the synthesized derivatives, compounds 3a and 3d demonstrated significant antitubercular activity, suggesting that hydroxy and chloro substituents on the aromatic ring enhance biological activity. Compound 3b exhibited comparatively better activity among the tested derivatives, while compound 3c showed moderate inhibitory activity against Mycobacterium tuberculosis. Compound 3e displayed moderate activity, possibly due to the electron-donating methoxy substituent influencing interaction with the biological target.
The observed activity may be attributed to the presence of the thiazolidinone nucleus combined with the pyridine carboxamide moiety, which are known pharmacophores for antimycobacterial activity. Variations in aromatic substituents significantly influenced the biological activity of the synthesized compounds.
5.3 Physicochemical Evaluation for Drug-Likeness
The physicochemical properties of compounds 3(a–e) were evaluated using Molinspiration software, and the calculated parameters indicated favourable drug-like characteristics for the synthesized compounds.
The molecular weights of the synthesized compounds ranged from 289.30 to 333.79 g/mol, which fall within the acceptable range suggested by Lipinski’s Rule of Five, indicating good oral drug-likeness. The calculated logP values ranged from 0.27 to 2.25, demonstrating moderate lipophilicity and suggesting favourable membrane permeability and absorption characteristics.
The TPSA values for all synthesized compounds were found within the acceptable range of 62.30–82.53 Ų, suggesting good oral bioavailability and adequate hydrogen bonding capacity. All compounds possessed acceptable numbers of hydrogen bond donors and acceptors and exhibited zero Lipinski violations, indicating compliance with standard drug-likeness criteria.
The number of rotatable bonds ranged from 3 to 5, suggesting moderate molecular flexibility, which may facilitate better fitting within biological target active sites while maintaining structural stability.
Overall, the physicochemical evaluation demonstrated that the synthesized compounds possess favourable drug-like properties and may serve as potential candidates for further biological investigation.
5.4 Molinspiration Bioactivity Score Analysis
The Molinspiration bioactivity score analysis was performed to predict possible biological interactions of the synthesized compounds with various pharmacological targets, including GPCR ligands, ion channel modulators, kinase inhibitors, nuclear receptor ligands, protease inhibitors, and enzyme inhibitors.
Among the synthesized derivatives, compound 3e exhibited comparatively better bioactivity scores, particularly as a protease inhibitor and enzyme inhibitor, suggesting potential interaction with enzymatic targets. Compound 3a also showed moderate activity as an enzyme inhibitor and a protease inhibitor.
Compounds 3b and 3d displayed moderate negative scores across most categories, indicating comparatively lower biological interactions. Compound 3c showed the least favourable bioactivity scores among the synthesized derivatives, suggesting weaker interactions with the studied biological targets.
The bioactivity results indicate that electron-donating substituents such as methoxy and hydroxy groups may enhance biological interactions, whereas heterocyclic substitutions such as furyl rings may reduce predicted activity.
5.5 ADMET Prediction Studies
The ADMET properties of compounds 3(a–e) were predicted using pkCSM software to evaluate their pharmacokinetic and toxicity profiles.
Absorption Properties
All synthesized compounds demonstrated high predicted intestinal absorption values (>93%), indicating favourable oral bioavailability. Compounds 3b and 3e exhibited relatively higher Caco-2 permeability values, suggesting enhanced membrane permeability and absorption potential.
Water solubility analysis indicated that compounds 3a and 3c possessed comparatively better aqueous solubility, while compounds 3d and 3e were less soluble due to increased lipophilicity.
Most compounds were predicted to be non-substrates of P-glycoprotein, except compounds 3b and 3e, which may undergo efflux-mediated transport.
Distribution Properties
The predicted VDss values suggested low tissue distribution for all compounds. The BBB permeability values indicated poor penetration across the blood–brain barrier, suggesting limited central nervous system exposure. Similarly, the CNS permeability values also confirmed poor CNS distribution.
Metabolism
The metabolic predictions revealed that compounds 3b, 3d, and 3e are substrates for CYP3A4 enzymes, indicating possible hepatic metabolism through cytochrome-mediated pathways. Compounds 3b and 3d were predicted to inhibit CYP1A2 and CYP2C19 enzymes, while compound 3e additionally inhibited CYP2C9.
Excretion
The total clearance values suggested slow elimination of the compounds from the body, which may contribute to prolonged biological activity. None of the synthesized derivatives was predicted to be substrates of renal OCT2 transporters.
Toxicity Studies
The toxicity predictions demonstrated that all synthesized compounds were non-mutagenic according to the AMES toxicity test. None of the compounds was predicted to inhibit hERG channels, indicating low cardiotoxic potential.
Compounds 3c, 3d, and 3e showed mild hepatotoxicity predictions, whereas compounds 3a and 3b were predicted to be non-hepatotoxic. Compound 3e exhibited possible skin sensitization, while other compounds were predicted to be non-sensitizers.
Overall, the ADMET analysis suggested that the synthesized compounds possess acceptable pharmacokinetic profiles with relatively low toxicity risks.
5.6 Molecular Docking Studies
Molecular docking studies were performed using AutoDock Vina to evaluate binding interactions between the synthesized compounds and the selected Mycobacterium tuberculosis enoyl-acyl carrier protein reductase (InhA) target protein. Before docking, the target protein structure was prepared by removing water molecules and co-crystallized ligands, followed by the addition of polar hydrogen atoms and Kollman charges. The prepared protein and ligand structures were converted into PDBQT format for docking analysis.
The docking scores indicated favourable binding affinities toward the protein active site. Among the synthesized compounds, 3a and 3b demonstrated the best docking scores (–8.6 kcal/mol) within cavity C1, indicating strong interactions with the active site of the target protein. Compound 3e also exhibited good binding affinity with a docking score of –8.3 kcal/mol.
The cavity volumes ranged from 617 to 1108 ų, suggesting adequate space for ligand accommodation within the binding pockets. The docking dimensions and cavity centers confirmed stable fitting of the synthesized compounds within the receptor active sites.
The enhanced binding affinity observed for compounds containing hydroxy and phenyl substituents may be attributed to hydrogen bonding interactions and hydrophobic contacts with amino acid residues present in the active site. The docking studies supported the experimental antitubercular activity results and indicated that the synthesized thiazolidinone derivatives may act as promising antitubercular agents.
5.7 Structure–Activity Relationship (SAR)
. Electron-withdrawing substituents such as chloro groups enhanced activity, possibly due to increased hydrophobic interactions within the target binding pocket. Hydroxy substitution may improve activity through hydrogen bonding interactions with active-site residues, thereby enhancing ligand binding affinity.
Compounds containing unsubstituted phenyl rings also demonstrated good activity, suggesting that aromatic hydrophobic interactions contribute positively to receptor binding. In contrast, the methoxy-substituted derivative exhibited comparatively moderate activity, indicating that strong electron-donating groups may reduce optimal ligand–target interactions.
The furyl-substituted derivative displayed moderate biological activity, which may be attributed to altered electronic distribution and weaker hydrophobic interactions compared with phenyl-substituted analogues.
Overall, the SAR findings suggest that both electronic and hydrophobic properties of aromatic substituents play crucial roles in determining the antitubercular activity of thiazolidinone derivatives.
CONCLUSION
The present study successfully synthesized a series of novel N-(4-oxo-3-aryl-1,2-thiazolidin-2-yl) pyridine-3-carboxamide derivatives and characterized them using FT-IR and ¹H NMR spectral analyses. The spectral data confirmed the successful formation of the desired thiazolidinone framework.
Biological evaluation revealed promising antitubercular activity for all synthesized compounds against Mycobacterium tuberculosis H37Rv strain. Compounds containing hydroxy and chloro substituents demonstrated comparatively better activity, indicating the significant influence of substituent effects on biological performance.
The in-silico drug-likeness, ADMET, and molecular docking studies further supported the potential of these compounds as promising antitubercular candidates. Most compounds exhibited favourable pharmacokinetic properties, acceptable toxicity profiles, and strong binding affinities toward the Mycobacterium tuberculosis enoyl-acyl carrier protein reductase (InhA) target protein. Among the synthesized derivatives, compounds 3a, 3b, and 3d emerged as the most promising candidates for further optimization and biological investigation as potential antitubercular agents.
Overall, the combined experimental and computational findings suggest that these thiazolidinone derivatives may serve as promising lead compounds for the development of novel antitubercular drugs. Further structural optimization, detailed mechanistic studies, and in vivo investigations are warranted to establish their therapeutic potential and safety profiles.
References
- Bontha Venkata Subramanya Lokesh et al., (2019). “Synthesis, Biological Evaluation and Molecular Docking Studies of New Pyrazolines as Anti-Tubercular and Cytotoxic Agents.” Bioorganic Chemistry, pp. 1–12.
- Katharigatta N. Venugopala et al., (2021). “Molecular Docking and Anti-Tubercular Evaluation of Designed Compounds.” Journal of Molecular Structure, pp. 129245–129256.
- N. C. Desai et al., (2016). “Synthesis, Biological Evaluation and Molecular Docking Study of Novel Indole and Pyridine Based 1,3,4-Oxadiazole Derivatives as Potential Anti-Tubercular Agents.” Bioorganic & Medicinal Chemistry Letters, pp. 1776–1783.
- Sarvesh Kumar Pandey et al., (2020). “Synthesis, Biological Evaluation and Molecular Docking Studies of Novel Quinazolinones as Anti-Tubercular and Anti-Malarial Agents.” European Journal of Medicinal Chemistry, pp. 112867–112879.
- Mayursinh Zala et al., (2023). “Synthesis, Characterization, Anti-Tubercular Activity and Molecular Docking Studies of Pyrazolylpyrazoline-Clubbed Triazole and Tetrazole Hybrids.” Journal of Molecular Structure, pp. 135430–135438.
- Mubarak H. Shaikh et al., (2015). “1,2,3-Triazole Derivatives as Anti-Tubercular Agents: Synthesis, Biological Evaluation and Molecular Docking Study.” Medicinal Chemistry Research, pp. 149–158.
- Esra Tatar et al., (2016). “Design, Synthesis and Molecular Docking Studies of a Conjugated Thiadiazole-Thiourea Scaffold as Anti-Tubercular Agents.” Bioorganic Chemistry, pp. 1–9.
- Shankar Thapa et al., (2024). “Synthesis, Molecular Docking, Molecular Dynamic Simulation Studies and Anti-Tubercular Activity Evaluation of Substituted Benzimidazole Derivatives.” Journal of Biomolecular Structure and Dynamics, pp. 1–15.
- Pimonluck Sittikornpaiboon et al., (2016). “Molecular Docking Study of Mycobacterium tuberculosis Dihydrofolate Reductase in Complex with 2,4-Diaminopyrimidine Analogues.” Journal of Molecular Graphics and Modelling, pp. 223–230.
- Pran Kishore Deb et al., (2021). “In Vitro Anti-TB Properties, In Silico Target Validation, Molecular Docking and Dynamics Studies of Substituted 1,2,4-Oxadiazole Analogues against Mycobacterium tuberculosis.” Computers in Biology and Medicine, pp. 104984–104995.
- Nilam H. Lalavani et al., (2022). “Synthesis, Pharmacokinetic and Molecular Docking Studies of New Benzohydrazide Derivatives Possessing Anti-Tuberculosis Activity.” Journal of Molecular Structure, pp. 132569–132577.
- M. J. Zala et al., (2022). “Synthesis and Molecular Docking Study of Arylsulfanyl Pyrazolylpyrazoline Derivatives as Anti-Tubercular Agents.” Synthetic Communications, pp. 2438–2450.
- S. Bodige et al., (2020). “Design, Synthesis, Anti-Tubercular and Antibacterial Activities of 1,3,5-Triazinyl Carboxamide Derivatives and In Silico Docking Studies.” Chemical Biology & Drug Design, pp. 601–612.
- B. Srinu et al., (2019). “Synthesis, Anti-Tubercular Activity and Molecular Docking Studies of Novel Pyrimidine-5-Carboxamides.” Medicinal Chemistry Research, pp. 1200–1210.
- H. Surya Prakash Rao et al., (2023). “Design, Synthesis, Molecular Docking and Biological Activity of Pyrazolo[3,4-b] Pyridines against M. tuberculosis.” Journal of Heterocyclic Chemistry, pp. 1450–1462.
- Sravanthi Siliveri et al., (2024). “Design, Synthesis, In Silico Molecular Docking, ADMET Studies and Biological Evaluation of Novel Diphenyl Pyrazole Incorporated Triazolyl Benzene Sulfonamides.” Results in Chemistry, pp. 101234–101245.
- Sushil K. Kashaw et al. (2011). “Synthesis and Anti-Tubercular Evaluation of Novel Thiazolidinones and Brominated Derivatives.” European Journal of Medicinal Chemistry, pp. 384–390.
- Nazar Trotsko et al., (2021). “Anti-Tubercular Properties of Thiazolidin-4-ones.” Bioorganic Chemistry, pp. 105321–105335.
- Himaja Malipeddi et al., (2012). “Synthesis and Anti-Tubercular Activity of Novel Thiazolidinone Derivatives.” International Journal of Pharmaceutical Sciences and Research, pp. 221–227.
- Haresh Oza et al., (1998). “Synthesis and Anti-Tubercular Activity of Thiazolidinone Derivatives.” Indian Journal of Chemistry, pp. 1045–1049.
