Pharmacophore-Based Virtual Screening,Quantum Mechanics Calculations,and Molecular Dynamics Simulation Approaches Identified Potential Natural Antiviral Drug Candidates against MERS-CoV S1-NTD

Middle East respiratory syndrome coronavirus (MERS-CoV) is a highly infectious zoonotic virus first reported into the human population in September 2012 on the Arabian Peninsula. The virus causes severe and often lethal respiratory illness in humans with an unusually high fatality rate. The N-terminal domain (NTD) of receptor-binding S1 subunit of coronavirus spike (S) proteins can recognize a variety of host protein and mediates entry into human host cells. Blocking the entry by targeting the S1-NTD of the virus can facilitate the development of effective antiviral drug candidates against the pathogen. Therefore, the study has been designed to identify effective antiviral drug candidates against the MERS-CoV by targeting S1-NTD. Initially, a structure-based pharmacophore model (SBPM) to the active site (AS) cavity of the S1-NTD has been generated, followed by pharmacophore-based virtual screening of 11,295 natural compounds. Hits generated through the pharmacophore-based virtual screening have re-ranked by molecular docking and further evaluated through the ADMET properties. The compounds with the best ADME and toxicity properties have been retrieved, and a quantum mechanical (QM) based density-functional theory (DFT) has been performed to optimize the geometry of the selected compounds. Three optimized natural compounds, namely Taiwanhomoflavone B (Amb23604132), 2,3-Dihydrohinokiflavone (Amb23604659), and Sophoricoside (Amb1153724), have exhibited substantial docking energy >−9.00 kcal/mol, where analysis of frontier molecular orbital (FMO) theory found the low chemical reactivity correspondence to the bioactivity of the compounds. Molecular dynamics (MD) simulation confirmed the stability of the selected natural compound to the binding site of the protein. Additionally, molecular mechanics generalized born surface area (MM/GBSA) predicted the good value of binding free energies (ΔG bind) of the compounds to the desired protein. Convincingly, all the results support the potentiality of the selected compounds as natural antiviral candidates against the MERS-CoV S1-NTD.     

Exploring Natural Product Activity and Species Source Candidates for Hunting ABCB1 Transporter Inhibitors: An In Silico Drug Discovery Study

The P-glycoprotein (P-gp/ABCB1) is responsible for a xenobiotic efflux pump that shackles intracellular drug accumulation. Additionally, it is included in the dud of considerable antiviral and anticancer chemotherapies because of the multidrug resistance (MDR) phenomenon. In the search for prospective anticancer drugs that inhibit the ABCB1 transporter, the Natural Product Activity and Species Source (NPASS) database, containing >35,000 molecules, was explored for identifying ABCB1 inhibitors. The performance of AutoDock4.2.6 software to anticipate ABCB1 docking score and pose was first assessed according to available experimental data. The docking scores of the NPASS molecules were predicted against the ABCB1 transporter. Molecular dynamics (MD) simulations were conducted for molecules with docking scores lower than taxol, a reference inhibitor, pursued by molecular mechanics-generalized Born surface area (MM-GBSA) binding energy estimations. On the basis of MM-GBSA calculations, five compounds revealed promising binding affinities as ABCB1 inhibitors with ΔG_binding < −105.0 kcal/mol. The binding affinity and stability of the identified inhibitors were compared to the chemotherapeutic agent. Structural and energetical analyses unveiled great steadiness of the investigated inhibitors within the ABCB1 active site throughout 100 ns MD simulations. Conclusively, these findings point out that NPC104372, NPC475164, NPC2313, NPC197736, and NPC477344 hold guarantees as potential ABCB1 drug candidates and warrant further in vitro/in vivo tests.     

Chemical Profiles and Pharmacological Properties with in Silico Studies on Elatostema papillosum Wedd

The current study attempted, for the first time, to qualitatively and quantitatively determine the phytochemical components of Elatostema papillosum methanol extract and their biological activities. The present study represents an effort to correlate our previously reported biological activities with a computational study, including molecular docking, and ADME/T (absorption, distribution, metabolism, and excretion/toxicity) analyses, to identify the phytochemicals that are potentially responsible for the antioxidant, antidepressant, anxiolytic, analgesic, and anti-inflammatory activities of this plant. In the gas chromatography-mass spectroscopy analysis, a total of 24 compounds were identified, seven of which were documented as being bioactive based on their binding affinities. These seven were subjected to molecular docking studies that were correlated with the pharmacological outcomes. Additionally, the ADME/T properties of these compounds were evaluated to determine their drug-like properties and toxicity levels. The seven selected, isolated compounds displayed favorable binding affinities to potassium channels, human serotonin receptor, cyclooxygenase-1 (COX-1), COX-2, nuclear factor (NF)-κB, and human peroxiredoxin 5 receptor proteins. Phytol acetate, and terpene compounds identified in E. papillosum displayed strong predictive binding affinities towards the human serotonin receptor. Furthermore, 3-trifluoroacetoxypentadecane showed a significant binding affinity for the KcsA potassium channel. Eicosanal showed the highest predicted binding affinity towards the human peroxiredoxin 5 receptor. All of these findings support the observed in vivo antidepressant and anxiolytic effects and the in vitro antioxidant effects observed for this extract. The identified compounds from E. papillosum showed the lowest binding affinities towards COX-1, COX-2, and NF-κB receptors, which indicated the inconsequential impacts of this extract against the activities of these three proteins. Overall, E. papillosum appears to be bioactive and could represent a potential source for the development of alternative medicines; however, further analytical experiments remain necessary.     

Oxazinethione Derivatives as a Precursor to Pyrazolone and Pyrimidine Derivatives: Synthesis,Biological Activities,Molecular Modeling,ADME, and Molecular Dynamics Studies

In this study, we used oxazinethione as a perfect precursor to synthesize new pyrimidine and pyrazole derivatives with potent biological activities. Biological activities were determined for all compounds against A. flavus, E. coli, S. aureus, and F. moniliform. Compounds 3, 4a-b, and 5 exhibited higher activities toward A. flavus, E. coli, S. aureus, and F. moniliform; this was indicated through the MIC (minimum inhibitory concentration). At the same time, anticancer activities were determined through four cell lines, Ovcar-3, Hela, MCF-7, and LCC-MMk. The results obtained indicated that compound 5 was the most potent compound for both cell lines. Molecular docking was studied by the MOE (molecular operating environment). The in silico ADME of compounds 2 and 5 showed good pharmacokinetic properties. The present research strengthens the applicability of these compounds as encouraging anticancer and antibacterial drugs. Moreover, JAGUAR module MD simulations were carried out at about 100 ns. In addition, spectroscopic studies were carried out to establish the reactions of the synthesized structure derivatives.     

Identification of Potential Inhibitors of MurD Enzyme of Staphylococcus aureus from a Marine Natural Product Library

Staphylococcus aureus is an opportunistic pathogen that can cause fatal bacterial infections. MurD catalyzes the formation of peptide bond between UDP-N-acetylehyl-l-alanine and d-glutamic acid, which plays an important role in the synthesis of peptidoglycan and the formation of cell wall by S. aureus. Because S. aureus is resistant to most existing antibiotics, it is necessary to develop new inhibitors. In this study, Schrodinger 11.5 Prime homology modeling was selected to prepare the protein model of MurD enzyme, and its structure was optimized. We used a virtual screening program and similarity screening to screen 47163 compounds from three marine natural product libraries to explore new inhibitors of S. aureus. ADME provides analysis of the physicochemical properties of the best performing compounds during the screening process. To determine the stability of the docking effect, a 100 ns molecular dynamics was performed to verify how tightly the compound was bound to the protein. By docking analysis and molecular dynamics analysis, both 46604 and 46608 have strong interaction with the docking pocket, have good pharmacological properties, and maintain stable conformation with the target protein, so they have a chance to become drugs for S. aureus. Through virtual screening, similarity screening, ADME study and molecular dynamics simulation, 46604 and 46608 were selected as potential drug candidates for S. aureus.     

Synthesis,pharmacological evaluation,and molecular modeling studies of novel isatin hybrids as potential anticancer agents

A novel series of isatin hybrids 5a-g was designed, synthesized, and characterized spectroscopically. The synthesized compounds were evaluated for their cytotoxic activity against the human breast cancer cell line (MCF-7) by in vitro MTT assay. Amongst the tested compounds, 5e compound bearing benzyl moiety at N₄ piperazine was found to be the most active with the promising IC₅₀ (12.47 µM). Moreover, the active compounds 5e and 5g were subjected to antitumor evaluation (in vivo) against Dalton’s ascitic lymphoma (DAL) cell line and the results suggested that the best active compound 5e can normalize the blood picture in comparison to the standard drug. An in silico molecular docking study using the crystal structure of Hsp90 protein described the role of significant protein–ligand interactions and revealed more insights into the binding mode. The drug-likeliness of the compounds was predicted based on Lipinski's rule of five and pharmacokinetic ADME parameters. Hence, the synthesized isatin hybrids could be novel starting point anticancer lead compounds demonstrating drug-like properties which can be explored further for anticancer drug discovery.     

Assay Development and Identification of the First Plasmodium falciparum 7,8-dihydro-6-hydroxymethylpterin-pyrophosphokinase Inhibitors

In the fight towards eradication of malaria, identifying compounds active against new drug targets constitutes a key approach. Plasmodium falciparum 7,8-dihydro-6-hydroxymethylpterin-pyrophosphokinase (PfHPPK) has been advanced as a promising target, as being part of the parasite essential folate biosynthesis pathway while having no orthologue in the human genome. However, no drug discovery efforts have been reported on this enzyme. In this study, we conducted a three-step screening of our in-house antifolate library against PfHPPK using a newly designed PfHPPK-GFP protein construct. Combining virtual screening, differential scanning fluorimetry and enzymatic assay, we identified 14 compounds active against PfHPPK. Compounds’ binding modes were investigated by molecular docking, suggesting competitive binding with the HMDP substrate. Cytotoxicity and in vitro ADME properties of hit compounds were also assessed, showing good metabolic stability and low toxicity. The most active compounds displayed low micromolar IC₅₀ against drug-resistant parasites. The reported hit compounds constitute a good starting point for inhibitor development against PfHPPK, as an alternative approach to tackle the malaria parasite.     

A Comprehensive In Silico Exploration of Pharmacological Properties,Bioactivities, Molecular Docking,and Anticancer Potential of Vieloplain F from Xylopia vielana Targeting B-Raf Kinase

Compounds derived from plants have several anticancer properties. In the current study, one guaiane-type sesquiterpene dimer, vieloplain F, isolated from Xylopia vielana species, was tested against B-Raf kinase protein (PDB: 3OG7), a potent target for melanoma. A comprehensive in silico analysis was conducted in this research to understand the pharmacological properties of a compound encompassing absorption, distribution, metabolism, excretion, and toxicity (ADMET), bioactivity score predictions, and molecular docking. During ADMET estimations, the FDA-approved medicine vemurafenib was hepatotoxic, cytochrome-inhibiting, and non-cardiotoxic compared to the vieloplain F. The bioactivity scores of vieloplain F were active for nuclear receptor ligand and enzyme inhibitor. During molecular docking experiments, the compound vieloplain F has displayed a higher binding potential with −11.8 kcal/mol energy than control vemurafenib −10.2 kcal/mol. It was shown that intermolecular interaction with the B-Raf complex and the enzyme’s active gorge through hydrogen bonding and hydrophobic contacts was very accurate for the compound vieloplain F, which was then examined for MD simulations. In addition, simulations using MM-GBSA showed that vieloplain F had the greatest propensity to bind to active site residues. The vieloplain F has predominantly represented a more robust profile compared to control vemurafenib, and these results opened the road for vieloplain F for its utilization as a plausible anti-melanoma agent and anticancer drug in the next era.     

Synthesis,anticancer and antimicrobial evaluation of new pyridyl and thiazolyl clubbed hydrazone scaffolds

Abstract

A series of new hydrazones bearing pyridyl and thiazolyl scaffolds have been synthesized and evaluated for their in vitro anticancer and antimicrobial activities. The anticancer activity was evaluated against the A549 lung cancer cell line. The eight hydrazone derivatives have shown better anticancer activity than positive control doxorubicin against the A549 lung cancer cell line. The antimicrobial activity was evaluated against bacterial and fungal pathogens by using well diffusion method. The four hydrazone derivatives have displayed good antimicrobial activities. Molecular docking studies of the synthesized hydrazone derivatives revealed good binding via hydrogen bond interactions with key residues on active sites as well as neighboring residues with an active site of Focal adhesion kinase (PDB ID 2JKO). A computational study for the prediction of absorption, distribution, metabolism, and excretion (ADME) properties of all compounds has also been performed.     

A pharmacoinformatic approach on Cannabinoid receptor 2 (CB2) and different small molecules: Homology modelling,molecular docking,MD simulations,drug designing and ADME analysis

CB2 receptor belongs to the family of G-protein coupled receptors (GPCRs), which extensively controls a range of pointer transduction. CB2 plays an essential role in the immune system. It also associates in the pathology of different ailment conditions. In this scenario, the synthetic drugs are inducing side effects to the human beings after the drug use. Therefore, this study is seeking novel alternate drug molecules with least side effects than conventional drugs. The alternative drug molecules were chosen from the natural sources. These molecules were selected from cyanobacteria with the help of earlier research findings. The target and ligand molecules were obtained from recognized databases. The bioactive molecules are selected from various cyanobacterial species, which are selected by their biological and pharmacological properties, after, which we incorporated to the crucial findings such as homology modelling, molecular docking, MD simulations along with absorption, distribution, metabolism, and excretion (ADME) analysis. Initially, the homology modelling was performed to frame the target from unknown sequences of CB2, which revealed 44% of similarities and 66% of identities with the A2A receptor. Subsequently, the CB2 protein molecule has docked with already known and prepared bioactive molecules, agonists and antagonist complex. In the present study, the agonists (5) and antagonist (1) were also taken for comparing the results with natural molecules. At the end of the docking analysis, the cyanobacterial molecules and an antagonist TNC-201 are revealed better docking scores with well binding contacts than the agonists. Especially, the usneoidone shows better results than other cyanobacterial molecules, and it is very close docking scores with that of TCN-201. Therefore, the usneoidone has incorporated to MD simulation with Cannabinoid receptors 2 (CB2). In MD simulations, the complex (CB2 and usneoidone) reveals better stability in 30 ns. Based on the computational outcome, we concluded that usneoidone is an effectual and appropriate drug candidate for activating CB2 receptors and it will be serving as a better component for the complications of CB2. Moreover, these computational approaches can be motivated to discover novel drug candidates in the pharmacological and healthcare sectors.     

Synthesis and biological evaluation of theophylline acetohydrazide hydrazone derivatives as antituberculosis agents

A series of small molecules, theophylline acetohydrazide hydrazone derivatives were obtained via condensation of theophylline-7-acetohydrazide with different aromatic/heterocyclic aldehydes. The compounds were synthesized and characterized by using conventional methods. Further, the compounds and standard drugs were evaluated against Mycobacterium tuberculosis H37Rv strain, the activity obtained was varying depending on the functional group attached to theophylline acetohydrazide hydrazone compounds. Among these, Br-substituted compounds showed more potent against M. tuberculosis with MIC 3.6–4 μM and better than the reference drugs used. The molecular docking studies have shown the possible binding modes of the compounds with protein (PDB ID: 4RHX); the compound 4h has shown highest glide score and binding energy. For all compounds, ADME properties were predicted.     

吲哚衍生物类VEGFR-2酪氨酸激酶抑制剂的从头设计

以血管内皮生长因子受体-2(VEGFR-2)酪氨酸激酶的晶体结构为基础, 采用从头药物设计方法, 设计了一系列吲哚类化合物, 并用类药性和分子对接进行了筛选, 最后得到10个对接能量较低的化合物分子, 对具有最低结合能的化合物与VEGFR-2酪氨酸激酶的复合物进行了10 ns的分子动力学模拟, 并对其结合模式进行了分析. 这些化合物结构新颖, 可能作为抗肿瘤的先导化合物或候选药物. 本文结果为VEGFR-2酪氨酸激酶抑制剂的进一步改造、 设计及合成提供了理论基础, 并有助于开发高活性和高选择性的抗肿瘤药物.     

Novel and Potential Small Molecule Scaffolds as DYRK1A Inhibitors by Integrated Molecular Docking-Based Virtual Screening and Dynamics Simulation Study

The dual-specificity tyrosine phosphorylation-regulated kinase 1A (DYRK1A) is a novel, promising and emerging biological target for therapeutic intervention in neurodegenerative diseases, especially in Alzheimer’s disease (AD). The molMall database, comprising rare, diverse and unique compounds, was explored for molecular docking-based virtual screening against the DYRK1A protein, in order to find out potential inhibitors. Ligands exhibiting hydrogen bond interactions with key amino acid residues such as Ile165, Lys188 (catalytic), Glu239 (gk+1), Leu241 (gk+3), Ser242, Asn244, and Asp307, of the target protein, were considered potential ligands. Hydrogen bond interactions with Leu241 (gk+3) were considered key determinants for the selection. High scoring structures were also docked by Glide XP docking in the active sites of twelve DYRK1A related protein kinases, viz. DYRK1B, DYRK2, CDK5/p25, CK1, CLK1, CLK3, GSK3β, MAPK2, MAPK10, PIM1, PKA, and PKCα, in order to find selective DYRK1A inhibitors. MM/GBSA binding free energies of selected ligand–protein complexes were also calculated in order to remove false positive hits. Physicochemical and pharmacokinetic properties of the selected six hit ligands were also computed and related with the proposed limits for orally active CNS drugs. The computational toxicity webserver ProTox-II was used to predict the toxicity profile of selected six hits (molmall IDs 9539, 11352, 15938, 19037, 21830 and 21878). The selected six docked ligand–protein systems were exposed to 100 ns molecular dynamics (MD) simulations to validate their mechanism of interactions and stability in the ATP pocket of human DYRK1A kinase. All six ligands were found to be stable in the ATP binding pocket of DYRK1A kinase.     

Computational Approaches for the Design of Novel Anticancer Compounds Based on Pyrazolo[3,4-d]pyrimidine Derivatives as TRAP1 Inhibitor

In the present in-silico study, various computational techniques were applied to determine potent compounds against TRAP1 kinase. The pharmacophore hypothesis DHHRR_1 consists of important features required for activity. The 3D QSAR study showed a statistically significant model with R² = 0.96 and Q² = 0.57. Leave one out (LOO) cross-validation (R² CV = 0.58) was used to validate the QSAR model. The molecular docking study showed maximum XP docking scores (−11.265, −10.532, −10.422, −10.827, −10.753 kcal/mol) for potent pyrazole analogs (42, 46, 49, 56, 43), respectively, with significant interactions with amino acid residues (ASP 594, CYS 532, PHE 583, SER 536) against TRAP1 kinase receptors (PDB ID: 5Y3N). Furthermore, the docking results were validated using the 100 ns MD simulations performed for the selected five docked complexes. The selected inhibitors showed relatively higher binding affinities than the TRAP1 inhibitor molecules present in the literature. The ZINC database was used for a virtual screening study that screened ZINC05297837, ZINC05434822, and ZINC72286418, which showed similar binding interactions to those shown by potent ligands. Absorption, distribution, metabolism, and excretion (ADME) analysis showed noticeable results. The results of the study may be helpful for the further development of potent TRAP1 inhibitors     

Role of bridging water molecules in GSK3β-inhibitor complexes: insights from QM/MM, MD, and molecular docking studies

The role of water molecules is increasingly gaining interest in drug design, and several studies have highlighted their paramount contributions to the specificity and the affinity of ligand binding. In this study, we employ the two-layer ONIOM-based quantum mechanics/molecular mechanics (QM/MM) calculations, molecular dynamics (MD) simulations, and molecular docking studies to investigate the effect of bridging water molecules at the GSK3β-inhibitors interfaces. The results obtained from the ONIOM geometry optimization and AIM analysis corroborated the presence of bridging water molecules that form hydrogen bonds with protein side chain of Thr138 and/or backbone of Gln185, and mediate interactions with inhibitors in the 10 selected GSK3β-inhibitor complexes. Subsequently, MD simulations carried out on a representative system of 1R0E demonstrated that the bridging water molecule is stable at the GSK3β-inhibitor interface and appears to contribute to the stability of the protein-inhibitor interactions. Furthermore, molecular docking studies of GSK3β-inhibitor complexes indicated that the inhibitors can increase binding affinities and the better docked conformation of inhibitors can be obtained by inclusion of the bridging water molecules, especially for the flexible inhibitors, in docking experiments into individual protein conformations. Our results elucidate the importance of bridging water molecules at the GSK3β-inhibitor interfaces and suggest that they might prove useful in rational drug design.     

De novo design based identification of potential HIV-1 integrase inhibitors: A pharmacoinformatics study

In the present study, pharmacoinformatics paradigms include receptor-based de novo design, virtual screening through molecular docking and molecular dynamics (MD) simulation are implemented to identify novel and promising HIV-1 integrase inhibitors. The de novodrug/ligand/molecule design is a powerful and effective approach to design a large number of novel and structurally diverse compounds with the required pharmacological profiles. A crystal structure of HIV-1 integrase bound with standard inhibitor BI-224436 is used and a set of 80,000 compounds through the de novo approach in LigBuilder is designed. Initially, a number of criteria including molecular docking, in-silico toxicity and pharmacokinetics profile assessments are implied to reduce the chemical space. Finally, four de novo designed molecules are proposed as potential HIV-1 integrase inhibitors based on comparative analyses. Notably, strong binding interactions have been identified between a few newly identified catalytic amino acid residues and proposed HIV-1 integrase inhibitors. For evaluation of the dynamic stability of the protein-ligand complexes, a number of parameters are explored from the 100 ns MD simulation study. The MD simulation study suggested that proposed molecules efficiently retained their molecular interaction and structural integrity inside the HIV-1 integrase. The binding free energy is calculated through the Molecular Mechanics Poisson-Boltzmann Surface Area (MM-PBSA) approach for all complexes and it also explains their thermodynamic stability. Hence, proposed molecules through de novo design might be critical to inhibiting the HIV-1 integrase.     

2-Phenylindole derivatives as anticancer agents: synthesis and screening against murine melanoma,human lung and breast cancer cell lines

Indole derivatives have attractive anticancer properties and may be a future hope for better anticancer drug(s) of low toxicity and high potency. In this paper, syntheses of 2-phenylindole derivatives have been described via Fischer indole synthesis through a one-pot solvent-free method. The synthesized compounds were screened for anticancer potential in vitro against murine melanoma (B16F10), human lung cancer (A549), and human breast cancer (MDA-MB-231) cell lines. The results highlighted that 2-phenylindole derivatives are also promising anticancer agents in case of melanoma and lung cancer along with the breast cancer. Molecular docking analyses with possible targets for melanoma (NEDD4-1) and lung cancer (EGFR) were also performed to understand specific interactions of 2-phenylindole derivatives with the amino acid residues of the receptors.     

Molecular Docking,Molecular Dynamics Simulations,Computational Screening to Design Quorum Sensing Inhibitors Targeting LuxP of <Emphasis Type="Italic">Vibrio harveyi</Emphasis> and Its Biological Evaluation

Quorum sensing (QS) plays an important role in the biofilm formation, production of virulence factors and stress responses in Vibrio harveyi. Therefore, interrupting QS is a possible approach to modulate bacterial behavior. In the present study, three docking protocols, such as Rigid Receptor Docking (RRD), Induced Fit Docking (IFD), and Quantum Polarized Ligand Docking (QPLD) were used to elucidate the binding mode of boronic acid derivatives into the binding pocket of LuxP protein in V. harveyi. Among the three docking protocols, IFD accurately predicted the correct binding mode of the studied inhibitors. Molecular dynamics (MD) simulations of the protein-ligand complexes indicates that the inter-molecular hydrogen bonds formed between the protein and ligand complex remains stable during the simulation time. Pharmacophore and shape-based virtual screening were performed to find selective and potent compounds from ChemBridge database. Five hit compounds were selected and subjected to IFD and MD simulations to validate the binding mode. In addition, enrichment calculation was performed to discriminate and separate active compounds from the inactive compounds. Based on the computational studies, the potent Bicyclo [2.2.1] hept-5-ene-2,3-dicarboxylic acid-2,6-dimethylpyridine 1-oxide (ChemBridge_5144368) was selected for in vitro assays. The compound exhibited dose dependent inhibition in bioluminescence and also inhibits biofilm formation in V. harveyi to the level of 64.25 %. The result from the study suggests that ChemBridge_5144368 could serve as an anti-quorum sensing molecule for V. harveyi.     

Evaluation of inhibitory effects of benzothiazole and 3-amino-benzothiazolium derivatives on DNA topoisomerase II by molecular modeling studies

There has been considerable interest in DNA topoisomerases over the last decade, as they have been shown to be one of the major cellular targets in anticancer drug development. Previously we synthesized some benzothiazole derivatives and corresponding benzothiazolium forms, and tested their DNA inhibitory activity to develop novel antitumor agents. Among the 12 prepared compounds, compound BM3 (3-aminobenzothiazole-3-ium 4-methylbenzene sulfonate) exhibited extreme topoisomerase II inhibitory activity compared with the reference drug etoposide. We also tried to determine the DNA and enzyme binding abilities of BM3 and found that BM3 acted on topoisomerase II first at low doses, while it had also showed DNA minor groove binding properties at higher doses. In this study the interactions between DNA topoisomerase II and the compounds were examined in detail by molecular modelling studies such as molecular docking and pharmacophore analysis performed using Discovery Studio 3.5. As a result, it was found that benzothiazolium compounds exhibited a totally different mechanism than benzothiazoles by binding to the different amino acids at the active site of the protein molecule. 3-Aminobenzothiazoliums are worthy of carrying onto anticancer studies; BM3 especially would be a good anticancer candidate for preclinical studies.