Computational evaluation of factors governing catalytic 2-keto acid decarboxylation

Lymphatic filariasis (LF) is a vector borne infectious disease caused by the nematode Wuchereria bancrofti, Brugia malayi, and Brugia timori. Over 120 million people are affected by LF in the world, of which two-thirds are in Asia. The infection restricts the normal flow of lymph from the infected area resulting in swelling of the extremities and causing permanent disability. As the available drugs for the treatment of LF are becoming ineffective due to the development of resistance, there is an urgent need to find new leads for drug development. In this study, asparaginyl-tRNA synthetase (AsnRS; PDB ID: 2XGT) essential for the protein bio-synthesis in the filarial nematode was used to carry out virtual screening (VS) of plant constituents from traditional Chinese medicine (TCM) database. Docking as well as E-pharmacophore based VS were carried out to identify the hits. The top scoring hits, Agri 1 (1,3,8-trihydroxy-4,5-dimethoxyxanthen-9-one-3-O-beta-D-glucopyranoside) and Agri 2 (5,7-dihydroxy-2-propylchromone 7-O-beta-D-glucopyranoside), constituents of Agrimonia pilosa, were selected for molecular dynamics (MD) simulation study for 10 ns. MD simulation showed that both the glycosides Agri 1 and Agri 2 were forming stable interactions with the target protein. Moreover, docking and MD simulation of the lead A (1,3,8-trihydroxy-4,5-dimethoxyxanthen-9-one; Mol. Wt.: 304.25; CLogP: 3.07) and lead B (5,7-dihydroxy-2-propylchromone; Mol. Wt.: 220.22; CLogP: 3.02), the aglycones of Agri 1 and Agri 2, respectively, were carried out with the target AsnRS. The in silico investigations of the aglycones suggest that the lead B could be a suitable fragment-like lead molecule for anti-filarial drug discovery.

Shape-based virtual screening,docking, and molecular dynamics simulations to identify Mtb-ASADH inhibitors

Aspartate β-semialdehyde dehydrogenase (ASADH) is a key enzyme for the biosynthesis of essential amino acids and several important metabolites in microbes. Inhibition of ASADH enzyme is a promising drug target strategy against Mycobacterium tuberculosis (Mtb). In this work, in silico approach was used to identify potent inhibitors of Mtb-ASADH. Aspartyl β-difluorophosphonate (β-AFP), a known lead compound, was used to understand the molecular recognition interactions (using molecular docking and molecular dynamics analysis). This analysis helped in validating the computational protocol and established the participation of Arg99, Glu224, Cys130, Arg249, and His256 amino acids as the key amino acids in stabilizing ligand–enzyme interactions for effective binding, an essential feature is H-bonding interactions with the two arginyl residues at the two ends of the ligand. Best binding conformation of β-AFP was selected as a template for shape-based virtual screening (ZINC and NCI databases) to identify compounds that competitively inhibit the Mtb-ASADH. The top rank hits were further subjected to ADME and toxicity filters. Final filter was based on molecular docking analysis. Each screened molecule carries the characteristics of the highly electronegative groups on both sides separated by an average distance of 6?Å. Finally, the best predicted 20 compounds exhibited minimum three H-bonding interactions with Arg99 and Arg249. These identified hits can be further used for designing the more potent inhibitors against ASADH family. MD simulations were also performed on two selected compounds (NSC4862 and ZINC02534243) for further validation. During the MD simulations, both compounds showed same H-bonding interactions and remained bound to key active residues of Mtb-ASADH.     

Origins of the specificity of inhibitor P218 toward wild-type and mutant PfDHFR: a molecular dynamics analysis

Molecular dynamics simulations were performed to evaluate the origin of the antimalarial effect of the lead compound P218. The simulations of the ligand in the cavities of wild-type, mutant Plasmodium falciparum Dihydrofolate Reductase (PfDHFR) and the human DHFR revealed the differences in the atomic-level interactions and also provided explanation for the specificity of this ligand toward PfDHFR. The binding free energy estimation using Molecular Mechanics Poisson-Boltzmann Surface Area method revealed that P218 has higher binding affinity (~ ?30 to ?35 kcal/mol) toward PfDHFR (both in wild-type and mutant forms) than human DHFR (~ ?22 kcal/mol), corroborating the experimental observations. Intermolecular hydrogen bonding analysis of the trajectories showed that P218 formed two stable hydrogen bonds with human DHFR (Ile7 and Glu30), wild-type and double-mutant PfDHFR’s (Asp54 and Arg122), while it formed three stable hydrogen bonds with quadruple-mutant PfDHFR (Asp54, Arg59, and Arg122). Additionally, P218 binding in PfDHFR is stabilized by hydrogen bonds with residues Ile14 and Ile164. It was found that mutant residues do not reduce the binding affinity of P218 to PfDHFR, in contrast, Cys59Arg mutation strongly favors inhibitor binding to quadruple-mutant PfDHFR. The atomistic-level details explored in this work will be highly useful for the design of non-resistant novel PfDHFR inhibitors as antimalarial agents.     

Biochemical characterization of Recombinase A from Wolbachia endosymbiont of filarial nematode Brugia malayi (wBmRecA)

Lymphatic filariasis is a debilitating disease that affects over 890 million people in 49 countries. A lack of vaccines, non-availability of adulticidal drugs, the threat of emerging drug resistance against available chemotherapeutics and an incomplete understanding of the immunobiology of the disease have sustained the problem. Characterization of Wolbachia proteins, the bacterial endosymbiont which helps in the growth and development of filarial worms, regulates fecundity in female worms and mediates immunopathogenesis of Lymphatic Filariasis, is an important approach to gain insights into the immunopathogenesis of the disease. In this study, we carried out extensive biochemical characterization of Recombinase A from Wolbachia of the filarial nematode Brugia malayi (wBmRecA) using an Electrophoretic Mobility Shift Assay, an ATP binding and hydrolysis assay, DNA strand exchange reactions, DAPI displacement assay and confocal microscopy, and evaluated anti-filarial activity of RecA inhibitors. Confocal studies showed that wBmRecA was expressed and localised within B. malayi microfilariae (Mf) and uteri and lateral chord of adult females. Recombinant wBmRecA was biochemically active and showed intrinsic binding capacity towards both single-stranded DNA and double-stranded DNA that were enhanced by ATP, suggesting ATP-induced cooperativity. wBmRecA promoted ATP hydrolysis and DNA strand exchange reactions in a concentration-dependent manner, and its binding to DNA was sensitive to temperature, pH and salt concentration. Importantly, the anti-parasitic drug Suramin, and Phthalocyanine tetrasulfonate (PcTs)-based inhibitors Fe-PcTs and 3,4-Cu-PcTs, inhibited wBmRecA activity and affected the motility and viability of Mf. The addition of Doxycycline further enhanced microfilaricidal activity of wBmRecA, suggesting potential synergism. Taken together, the omnipresence of wBmRecA in B. malayi life stages and the potent microfilaricidal activity of RecA inhibitors suggest an important role of wBmRecA in filarial pathogenesis.     

Identification of hotspot regions of MurB oxidoreductase enzyme using homology modeling,molecular dynamics and molecular docking techniques

Despite the availability of effective chemotherapy and a moderately protective vaccine, new anti-tuberculosis agents are urgently needed to decrease the global incidence of tuberculosis (TB) disease. The MurB gene belongs to the bacterial cell wall biosynthesis pathway and is an essential drug target in Mycobacterium tuberculosis (Mtb) that has no mammalian counterparts. Here, we present an integrated approach involving homology modeling, molecular dynamics and molecular docking studies on Mtb-MurB oxidoreductase enzyme. A homology model of Mtb-MurB enzyme was built for the first time in order to carry out structure-based inhibitor design. The accuracy of the model was validated using different techniques. The molecular docking study on this enzyme was undertaken using different classes of well known MurB inhibitors. Estimation of binding free energy by docking analysis indicated the importance of Tyr155, Arg156, Ser237, Asn241 and His304 residues within the Mtb-MurB binding pocket. Our computational analysis is in good agreement with experimental results of site-directed mutagenesis. The present study should therefore play a guiding role in the experimental design of Mtb-MurB inhibitors for in vitro/in vivo analysis.     

Exploration of the chlorpyrifos escape pathway from acylpeptide hydrolases using steered molecular dynamics simulations

Acylpeptide hydrolases (APH) catalyze the removal of an N-acylated amino acid from blocked peptides. APH is significantly more sensitive than acetylcholinesterase, a target of Alzheimer’s disease, to inhibition by organophosphorus (OP) compounds. Thus, OP compounds can be used as a tool to probe the physiological functions of APH. Here, we report the results of a computational study of molecular dynamics simulations of APH bound to the OP compounds and an exploration of the chlorpyrifos escape pathway using steered molecular dynamics (SMD) simulations. In addition, we apply SMD simulations to identify potential escape routes of chlorpyrifos from hydrolase hydrophobic cavities in the APH-inhibitor complex. Two previously proposed APH pathways were reliably identified by CAVER 3.0, with the estimated relative importance of P1 > P2 for its size. We identify the major pathway, P2, using SMD simulations, and Arg526, Glu88, Gly86, and Asn65 are identified as important residues for the ligand leaving via P2. These results may help in the design of APH-targeting drugs with improved efficacy, as well as in understanding APH selectivity of the inhibitor binding in the prolyl oligopeptidase family.     

Differential reactivity of filarial antigens with human sera from bancroftian filariasis endemic zone

The reacting pattern of circulating filarial antigen fraction-2 fromWuchereria bancrofti and soluble antigen from adultBrugia malayi with bancroftian filarial sera were analysed by immunoblotting technique and enzyme linked immunosorbent assay. Microfilaraemic sera reacted specifically with proteins of molecular weight 200, 120, 97, 56, 54, 43, 26 and 17 kDa of circulating Filarial antigen fraction-2 and 44, 40, 38, 31, 22 and 18 kDa ofBrugia malayi adult soluble antigen. Clinical filarial sera identified protein molecules of 56, 54 and 42 kDa of circulating filarial antigen fraction-2 and 19, 16 and 14 kDa ofBrugia malayi adult soluble antigen. Some components of both the antigen preparation were also identified by endemic normal serai.e.proteins 120, 97, 62, 43 and 33 kDa of circulating filarial antigen fraction-2 and 170, 120, 43, 31 and 12 kDa ofBrugia malayi adult soluble antigen. One of the sodium dodecyl sulphate-polyacrylamide gel electropherosis fractions of circulating filarial antigen fraction-2 (CFA₂-8) andBrugia malayi adult soluble antigen fraction-6 when used in enzyme linked immunosorbent assay could differentiate microfilaraemic sera from endemic normal and clinical filarial sera. The other antigen fractions (CFA₂-2, 6 and 7 andBmA-2) showed a high geometric mean titre of filarial immunoglobulin G antibodies in endemic normal sera when compared to microfilaraemia and clinical filarial sera. These proteins need to be further studied to assess their involvement in protecting from filarial infection in an endemic area.     

Virtual screening to identify <Emphasis Type="Italic">Leishmania braziliensis N</Emphasis>-myristoyltransferase inhibitors: pharmacophore models,docking, and molecular dynamics

Leishmaniasis is caused by several protozoa species belonging to genus Leishmania that are hosted by humans and other mammals. Millions of new cases are recorded every year and the drugs available on the market do not show satisfactory efficacy and safety. A hierarchical virtual screening approach based on the pharmacophore model, molecular docking, and molecular dynamics was conducted to identify possible Leishmania braziliensis N-misristoyltransferase (LbNMT) inhibitors. The adopted pharmacophore model had three main features: four hydrophobic centers, four hydrogen-bond acceptor atoms, and one positive nitrogen center. The molecules (n=15,000) were submitted to alignment with the pharmacophore model and only 27 molecules aligned to model. Six molecules were submitted to molecular docking, using receptor PDB ID 5A27. After docking, the ZINC35426134 was a top-ranked molecule (? 64.61 kcal/mol). The molecule ZINC35426134 shows hydrophobic interactions with Phe82, Tyr209, Val370, and Leu391 and hydrogen bonds with Asn159, Tyr318, and Val370. Molecular dynamics simulations were performed with the protein in its APO and HOLO forms for 37 ns in order to assess the stability of the protein–ligand complex. Results showed that the HOLO form was more stable than the APO one, and it suggests that the ZINC35426134 binding stabilizes the enzyme. Therefore, the selected molecule has the potential to meet the herein proposed target.     

In Vitro,In Silico and In Vivo Studies of Ursolic Acid as an Anti-Filarial Agent

As part of our drug discovery program for anti-filarial agents from Indian medicinal plants, leaves of Eucalyptus tereticornis were chemically investigated, which resulted in the isolation and characterization of an anti-filarial agent, ursolic acid (UA) as a major constituent. Antifilarial activity of UA against the human lymphatic filarial parasite Brugia malayi using in vitro and in vivo assays, and in silico docking search on glutathione-s-transferase (GST) parasitic enzyme were carried out. The UA was lethal to microfilariae (mf; LC₁₀₀: 50; IC₅₀: 8.84 µM) and female adult worms (LC₁₀₀: 100; IC50: 35.36 µM) as observed by motility assay; it exerted 86% inhibition in MTT reduction potential of the adult parasites. The selectivity index (SI) of UA for the parasites was found safe. This was supported by the molecular docking studies, which showed adequate docking (LibDock) scores for UA (−8.6) with respect to the standard antifilarial drugs, ivermectin (IVM −8.4) and diethylcarbamazine (DEC-C −4.6) on glutathione-s-transferase enzyme. Further, in silico pharmacokinetic and drug-likeness studies showed that UA possesses drug-like properties. Furthermore, UA was evaluated in vivo in B. malayi-M. coucha model (natural infection), which showed 54% macrofilaricidal activity, 56% female worm sterility and almost unchanged microfilaraemia maintained throughout observation period with no adverse effect on the host. Thus, in conclusion in vitro, in silico and in vivo results indicate that UA is a promising, inexpensive, widely available natural lead, which can be designed and developed into a macrofilaricidal drug. To the best of our knowledge this is the first ever report on the anti-filarial potential of UA from E. tereticornis, which is in full agreement with the Thomson Reuter''s ‘Metadrug’ tool screening predictions.     

Exploring the structure and dynamics of nano-confined water molecules using molecular dynamics simulations

Confinement effects can lead to drastic changes in the structural and dynamical properties of water molecules. In this work, we have performed classical molecular dynamics simulations of endohedral fullerenes of type (H₂O)_n@C_m (n = 1, 12, 21, 62, 108 and m = 60, 180, 240, 500 and 720) to explore the effects of spherical confinement on water properties. It is shown that these confined water molecules can form distinct solvation pattern depending upon the available space inside the fullerene cavity. For the systems with smaller diameter, cage-like structure is predominant whereas bulk-like structure is observed for larger fullerenes. The orientational relaxation of these confined water molecules showed slower relaxation as the cavity diameter increases except for the (H₂O)₂₁@C₂₄₀. In this case, stable cage-like structure hinders the overall dynamics of the trapped water molecules. Finally, we have calculated the hydrogen bond lifetimes from the hydrogen bond time correlation functions and compared with that of bulk water.     

Computational prediction of essential genes in an unculturable endosymbiotic bacterium, Wolbachia of Brugia malayi

Background  

Wolbachia (wBm) is an obligate endosymbiotic bacterium of Brugia malayi, a parasitic filarial nematode of humans and one of the causative agents of lymphatic filariasis. There is a pressing need for new drugs against filarial parasites, such as B. malayi. As wBm is required for B. malayi development and fertility, targeting wBm is a promising approach. However, the lifecycle of neither B. malayi nor wBm can be maintained in vitro. To facilitate selection of potential drug targets we computationally ranked the wBm genome based on confidence that a particular gene is essential for the survival of the bacterium.     

161 Discovery of potent KdsA inhibitors of Leptospira interrogans through homology modeling,docking, and molecular dynamics simulations

Leptospira interrogans is the foremost cause of human leptospirosis. Discovery of novel lead molecules for common drug targets of more than 250 Leptospira serovars is of significant research interest. Lipopolysaccharide (LPS) layer prevent entry of hydrophobic agents into the cell and protect structural integrity of the bacterium. KDO-8-phosphate synthase (KdsA) catalyzes the first step of KDO biosynthesis that leads to formation of inner core of LPS. KdsA was identified as a potential drug target against Leptospira interrogans through subtractive genomic approach, metabolic pathway analysis, and comparative analysis (Amineni et al., 2010). The present study rationalizes a systematic implementation of homology modeling, docking, and molecular dynamics simulations to discover potent KdsA inhibitors (Pradhan et al., 2013; Umamaheswari et al., 2010). A reliable tertiary structure of KdsA in complex with substrate PEP was constructed based on co-crystal structure of Aquifex aeolicus KdsA synthase with PEP using Modeller9v10. Geometry-based analog search for PEP was performed from LigandInfo database to generate an in house library of 352 ligands. The ligand data-set was docked into KdsA active site through three-stage docking technique (HTVS, SP, and XP) using Glidev5.7. Thirteen lead molecules were found to have better binding affinity compared to PEP (XP Gscore?=??7.38?kcal/mol; Figure 1). The best lead molecule (KdsA- lead1 docking complex) showed XP Gscore of ?10.26?kcal/mol and the binding interactions (Figure 2) were correlated favorably with PEP–KdsA interactions (Figure 1). Molecular dynamics simulations of KdsA– lead1 docking complex for 10?ns had revealed that the complex (Figure 3) remained stable in closer to physiological environmental condition. The predicted pharmacological properties of lead1 were well within the range of a drug molecule with good ADME profile, hence, would be intriguing towards development of potent inhibitor molecule against KdsA of Leptospira.     

Regulatory phosphorylation of cyclin-dependent kinase 2: insights from molecular dynamics simulations

The structures of fully active cyclin-dependent kinase-2 (CDK2) complexed with ATP and peptide substrate, CDK2 after the catalytic reaction, and CDK2 inhibited by phosphorylation at Thr14/Tyr15 were studied using molecular dynamics (MD) simulations. The structural details of the CDK2 catalytic site and CDK2 substrate binding box were described. Comparison of MD simulations of inhibited complexes of CDK2 was used to help understand the role of inhibitory phosphorylation at Thr14/Tyr15. Phosphorylation at Thr14/Tyr15 causes ATP misalignment for the phosphate-group transfer, changes in the Mg²⁺ coordination sphere, and changes in the H-bond network formed by CDK2 catalytic residues (Asp127, Lys129, Asn132). The inhibitory phosphorylation causes the G-loop to shift from the ATP binding site, which leads to opening of the CDK2 substrate binding box, thus probably weakening substrate binding. All these effects explain the decrease in kinase activity observed after inhibitory phosphorylation at Thr14/Tyr15 in the G-loop. Interaction of the peptide substrate, and the phosphorylated peptide product, with CDK2 was also studied and compared. These results broaden hypotheses drawn from our previous MD studies as to why a basic residue (Arg/Lys) is preferred at the P₊₂ substrate position. Figure View of the substrate binding site of the fully active cyclin-dependent kinase-2 (CDK2) (pT160-CDK2/cyclin A/ATP). The pThr160 activation site is located in the T-loop (yellow secondary structure). The G-loop, which partly forms the ATP binding site, is shown in blue. The Thr14 and Tyr15 inhibitory phosphorylation sites located in the G-loop are shown in licorice representation     

Identification of natural compound inhibitors against PfDXR: A hybrid structure-based molecular modeling approach and molecular dynamics simulation studies

In the present contribution, multicomplex-based pharmacophore studies were carried out on the structural proteome of Plasmodium falciparum 1-deoxy-D -xylulose-5-phosphate reductoisomerase. Among the constructed models, a representative model with complementary features, accountable for the inhibition was used as a primary filter for the screening of database molecules. Auxiliary evaluations of the screened molecules were performed via drug-likeness and molecular docking studies. Subsequently, the stability of the docked inhibitors was envisioned by molecular dynamics simulations, principle component analysis, and molecular mechanics-Poisson-Boltzmann surface area-based free binding energy calculations. The stability assessment of the hits was done by comparing with the reference (beta-substituted fosmidomycin analog, LC5) to prioritize more potent candidates. All the complexes showed stable dynamic behavior while three of them displayed higher binding free energy compared with the reference. The work resulted in the identification of the compounds with diverse scaffolds, which could be used as initial leads for the design of novel PfDXR inhibitors.     

Insights into molecular interactions between CaM and its inhibitors from molecular dynamics simulations and experimental data

In order to contribute to the structural basis for rational design of calmodulin (CaM) inhibitors, we analyzed the interaction of CaM with 14 classic antagonists and two compounds that do not affect CaM, using docking and molecular dynamics (MD) simulations, and the data were compared to available experimental data. The Ca²⁺-CaM-Ligands complexes were simulated 20 ns, with CaM starting in the “open” and “closed” conformations. The analysis of the MD simulations provided insight into the conformational changes undergone by CaM during its interaction with these ligands. These simulations were used to predict the binding free energies (ΔG) from contributions ΔH and ΔS, giving useful information about CaM ligand binding thermodynamics. The ΔG predicted for the CaM’s inhibitors correlated well with available experimental data as the r² obtained was 0.76 and 0.82 for the group of xanthones. Additionally, valuable information is presented here: I) CaM has two preferred ligand binding sites in the open conformation known as site 1 and 4, II) CaM can bind ligands of diverse structural nature, III) the flexibility of CaM is reduced by the union of its ligands, leading to a reduction in the Ca²⁺-CaM entropy, IV) enthalpy dominates the molecular recognition process in the system Ca²⁺-CaM-Ligand, and V) the ligands making more extensive contact with the protein have higher affinity for Ca²⁺-CaM. Despite their limitations, docking and MD simulations in combination with experimental data continue to be excellent tools for research in pharmacology, toward a rational design of new drugs.     

The solution structure of the forkhead box‐O DNA binding domain of Brugia malayi DAF‐16a

Brugia malayi is a parasitic nematode that causes lymphatic filariasis in humans. Here the solution structure of the forkhead DNA binding domain of Brugia malayi DAF‐16a, a putative ortholog of Caenorhabditis elegans DAF‐16, is reported. It is believed to be the first structure of a forkhead or winged helix domain from an invertebrate. C. elegans DAF‐16 is involved in the insulin/IGF‐I signaling pathway and helps control metabolism, longevity, and development. Conservation of sequence and structure with human FOXO proteins suggests that B. malayi DAF‐16a is a member of the FOXO family of forkhead proteins. Proteins 2014; 82:3490–3496. © 2014 Wiley Periodicals, Inc.     

Single-channel recording from adult Brugia malayi

Lymphatic filariasis is a significant cause of morbidity in humans. One of the causative agents is Brugia malayi a clade III nematode. Current therapeutic agents are effective against the microfilaria but less so against the adults residing in the host lymphatics. A large number of anthelmintics act on nematode ion channels including the nicotinic receptors found on nematode somatic muscle. The purpose of this study was to develop a preparation from adult B. malayi that was amenable to patch-clamp recording to facilitate the study of the ion channels present in this organism. We also present a preliminary characterization of the single-channel properties of nicotinic receptors from the adult musculature.     

A hybrid structural model of the complete Brugia malayi cytoplasmic asparaginyl-tRNA synthetase

Aminoacyl-tRNA synthetases are validated molecular targets for anti-infective drug discovery because of their essentiality in protein synthesis. Thanks to genome sequencing, it is now possible to systematically study aminoacyl-tRNA synthetases from human eukaryotic parasites as putative targets for novel drug discovery. As part of a program targeting class IIb asparaginyl-tRNA synthetases (AsnRS) from the parasitic nematode Brugia malayi for anti-filarial drugs, we report the complete structure of a eukaryotic AsnRS. Metazoan and fungal AsnRS differ from their bacterial homologues by the addition of a conserved N-terminal extension of about 110 residues whose structure we have determined by solution NMR for the B. malayi enzyme. In addition, we solved by X-ray crystallography a series of structures of the catalytically active N-terminally truncated enzyme (residues 112-548), allowing the structural basis for the mechanism of asparagine activation to be elucidated. The N-terminal domain contains a structured region with a novel fold featuring a lysine-rich helix that is shown by NMR to interact with tRNA. This is connected by an unstructured tether to the remainder of the enzyme, which is highly similar to the known structure of bacterial AsnRS. These data enable a model of the complete AsnRS-tRNA complex to be constructed.     

Pharmacophore generation,atom-based 3D-QSAR,molecular docking and molecular dynamics simulation studies on benzamide analogues as FtsZ inhibitors

FtsZ is an appealing target for the design of antimicrobial agent that can be used to defeat the multidrug-resistant bacterial pathogens. Pharmacophore modelling, molecular docking and molecular dynamics (MD) simulation studies were performed on a series of three-substituted benzamide derivatives. In the present study a five-featured pharmacophore model with one hydrogen bond acceptors, one hydrogen bond donors, one hydrophobic and two aromatic rings was developed using 97 molecules having MIC values ranging from .07 to 957 μM. A statistically significant 3D-QSAR model was obtained using this pharmacophore hypothesis with a good correlation coefficient (R² = .8319), cross validated coefficient (Q² = .6213) and a high Fisher ratio (F = 103.9) with three component PLS factor. A good correlation between experimental and predicted activity of the training (R² = .83) and test set (R² = .67) molecules were displayed by ADHRR.1682 model. The generated model was further validated by enrichment studies using the decoy test and MAE-based criteria to measure the efficiency of the model. The docking studies of all selected inhibitors in the active site of FtsZ protein showed crucial hydrogen bond interactions with Val 207, Asn 263, Leu 209, Gly 205 and Asn-299 residues. The binding free energies of these inhibitors were calculated by the molecular mechanics/generalized born surface area VSGB 2.0 method. Finally, a 15 ns MD simulation was done to confirm the stability of the 4DXD–ligand complex. On a wider scope, the prospect of present work provides insight in designing molecules with better selective FtsZ inhibitory potential.