Molecular modeling study of checkpoint kinase 1 inhibitors by multiple docking strategies and prime/MM-GBSA calculation

Selective Inhibition of HDAC1 by Macrocyclic Polypeptide for the Treatment of Glioblastoma: A Binding Mechanistic Analysis Based on Molecular Dynamics

Glioblastoma (GBM) is the most common and aggressive intracranial malignant brain tumor, and the abnormal expression of HDAC1 is closely correlated to the progression, recurrence and metastasis of GBM cells, making selective inhibition of HDAC1 a promising strategy for GBM treatments. Among all available selective HDAC1 inhibitors, the macrocyclic peptides have gained great attention due to their remarkable inhibitory selectivity on HDAC1. However, the binding mechanism underlying this selectivity is still elusive, which increases the difficulty of designing and synthesizing the macrocyclic peptide-based anti-GBM drug. Herein, multiple computational approaches were employed to explore the binding behaviors of a typical macrocyclic peptide FK228 in both HDAC1 and HDAC6. Starting from the docking conformations of FK228 in the binding pockets of HDAC1&6, relatively long MD simulation (500 ns) shown that the hydrophobic interaction and hydrogen bonding of E91 and D92 in the Loop2 of HDAC1 with the Cap had a certain traction effect on FK228, and the sub-pocket formed by Loop1 and Loop2 in HDAC1 could better accommodate the Cap group, which had a positive effect on maintaining the active conformation of FK228. While the weakening of the interactions between FK228 and the residues in the Loop2 of HDAC6 during the MD simulation led to the large deflection of FK228 in the binding site, which also resulted in the decrease in the interactions between the Linker region of FK228 and the previously identified key amino acids (H134, F143, H174, and F203). Therefore, the residues located in Loop1 and Loop2 contributed in maintaining the active conformation of FK228, which would provide valuable hints for the discovery and design of novel macrocyclic polypeptide HDAC inhibitors.

Structure-guided screening of chemical database to identify NS3-NS2B inhibitors for effective therapeutic application in dengue infection

Dengue infection is the most common arthropod-borne disease caused by dengue viruses, predominantly affecting millions of human beings annually. To find out promising chemical entities for therapeutic application in Dengue, in the current research, a multi-step virtual screening effort was conceived to screen out the entire "screening library" of the Asinex database. Initially, through "Lipinski rule of five" filtration criterion almost 0.6 million compounds were collected and docked with NS3-NS2B protein. Thereby, the chemical space was reduced to about 3500 compounds through the analysis of binding affinity obtained from molecular docking study in AutoDock Vina. Further, the "Virtual Screening Workflow" (VSW) utility of Schrödinger suite was used, which follows a stepwise multiple docking programs such as - high-throughput virtual screening (HTVS), standard precision (SP), and extra precision (XP) docking, and in postprocessing analysis the MM-GBSA based free binding energy calculation. Finally, five potent molecules were proposed as potential inhibitors for the dengue NS3-NS2B protein based on the investigation of molecular interactions map and protein-ligand fingerprint analyses. Different pharmacokinetics and drug-likeness parameters were also checked, which favour the potentiality of selected molecules for being drug-like candidates. The molecular dynamics (MD) simulation analyses of protein-ligand complexes were explained that NS3-NS2B bound with proposed molecules quite stable in dynamic states as observed from the root means square deviation (RMSD) and root means square fluctuation (RMSF) parameters. The binding free energy was calculated using MM-GBSA method from the MD simulation trajectories revealed that all proposed molecules possess such a strong binding affinity towards the dengue NS3-NS2B protein. Therefore, proposed molecules may be potential chemical components for effective inhibition of dengue NS3-NS2B protein subjected to experimental validation.

Pharmacophore-based designing of putative ROS-1 targeting agents for NSCLC

Non-small cell lung cancer (NSCLC) is a fatal non-immunogenic malignancy, and proto-oncogene receptor tyrosine kinase (ROS-1) is one of its clinically relevant biomarkers. In this context, herein, we report a series of benzimidazol-2-amine derivatives which were synthesized on the basis of the pharmacophore of ROS-1 and evaluated for anti-proliferative activity. For this, the in silico receptor-ligand pharmacophore model of ROS-1, previously published by our own group, was utilized to screen out an in-house database of small molecule heterocycles. Docking analysis of the selected compounds was carried out within the active site of wild-type (WT) ROS-1 as well as Gly2032Arg mutant ROS-1 protein, which confirmed the retention of conserved interaction between selected molecules and hinge region amino acids Glu2027 and Met2029. Docking was followed by molecular dynamics simulations for the stability of the complexes and calculation of the MM-GBSA score for binding affinity. Finally, compounds were synthesized and the anti-proliferative potential of compounds was evaluated using the A549 cell line. Compounds 3a and 3b presented significant GI₅₀ values between 23.0 and 25.4 μM, among all the tested compounds.

Discovery of Ganoderma lucidum triterpenoids as potential inhibitors against Dengue virus NS2B-NS3 protease

Dengue virus (DENV) infection causes serious health problems in humans for which no drug is currently available. Recently, DENV NS2B-NS3 protease has been proposed as a primary target for anti-dengue drug discovery due to its important role in new virus particle formation by conducting DENV polyprotein cleavage. Triterpenoids from the medicinal fungus Ganoderma lucidum have been suggested as pharmacologically bioactive compounds and tested as anti-viral agents against various viral pathogens including human immunodeficiency virus. However, no reports are available concerning the anti-viral activity of triterpenoids from Ganoderma lucidum against DENV. Therefore, we employed a virtual screening approach to predict the functional triterpenoids from Ganoderma lucidum as potential inhibitors of DENV NS2B-NS3 protease, followed by an in vitro assay. From in silico analysis of twenty-two triterpenoids of Ganoderma lucidum, four triterpenoids, viz. Ganodermanontriol (−6.291 kcal/mol), Lucidumol A (−5.993 kcal/mol), Ganoderic acid C2 (−5.948 kcal/mol) and Ganosporeric acid A (−5.983 kcal/mol) were predicted to be viral protease inhibitors by comparison to reference inhibitor 1,8-Dihydroxy-4,5-dinitroanthraquinone (−5.377 kcal/mol). These results were further studied for binding affinity and stability using the molecular mechanics/generalized Born surface area method and Molecular Dynamics simulations, respectively. Also, in vitro viral infection inhibition suggested that Ganodermanontriol is a potent bioactive triterpenoid.

Comprehensive evaluation of the MM-GBSA method on bromodomain-inhibitor sets

MM-PB/GBSA methods represent a higher-level scoring theory than docking. This study reports an extensive testing of different MM-GBSA scoring schemes on two bromodomain (BRD) datasets. The first set is composed of 24 BRPF1 complexes, and the second one is a nonredundant set constructed from the PDBbind and composed of 28 diverse BRD complexes. A variety of MM-GBSA schemes were analyzed to evaluate the performance of four protocols with different numbers of minimization and MD steps, 10 different force fields and three different water models. Results showed that neither additional MD steps nor unfixing the receptor atoms improved scoring or ranking power. On the contrary, our results underscore the advantage of fixing receptor atoms or limiting the number of MD steps not only for a reduction in the computational costs but also for boosting the prediction accuracy. Among Amber force fields tested, ff14SB and its derivatives rather than ff94 or polarized force fields provided the most accurate scoring and ranking results. The TIP3P water model yielded the highest scoring and ranking power compared to the others. Posing power was further evaluated for the BRPF1 set. A slightly better posing power for the protocol which uses both minimization and MD steps with a fixed receptor than the one which uses only minimization with a fully flexible receptor-ligand system was observed. Overall, this study provides insights into the usage of the MM-GBSA methods for screening of BRD inhibitors, substantiating the benefits of shorter protocols and latest force fields and maintaining the crystal waters for accuracy.

Stability behaviour of antiretroviral drugs and their combinations. 11: Characterization of interaction products of zidovudine and efavirenz, and evaluation of their anti HIV-1 activity, and physiochemical and ADMET properties

Zidovudine (ZDV) and efavirenz (EFV), which belong to two separate classes of antiretroviral drugs, viz., NRTI and NNRTI, respectively, were subjected to different stability test conditions alone and in solid mixtures to evaluate possibility of interaction among them. The exposed samples were analyzed by high performance liquid chromatography (HPLC) using a C18 column and a PDA detector. Two new peaks were observed in the sample in which 50 μl CH₃CN was added to increase the contact among the drugs, and which was subjected in open beaker to accelerated stability test condition of 40 °C/75%RH for 15 d. Subsequently, liquid chromatography-high resolution mass spectrometric (LC-HRMS) studies were carried out to obtain their accurate mass. The products were also isolated, and subjected to ¹H, ¹³C, DEPT-135, COSY, HSQC and HMBC nuclear magnetic resonance (NMR) studies. The collective information allowed their structural characterization as isomeric cycloaddition products of the two drugs. As these were novel compounds, they were subjected to testing for cytotoxicity and in vitro anti-HIV-1 activity against primary isolates HIV-1_UG070 (X4, subtype D) and HIV-1_VB59 (R5, subtype C) in TZM-bl cell line. The two were found to show weak activity against the standard drugs. The reason was sought through molecular docking studies, which highlighted that it was perhaps their comparative bigger molecular size than the drugs of both classes used currently in HIV therapy. Being previously unknown molecules, their in silico physicochemical and ADMET properties were also evaluated using ADMET Predictor™ and TOPKAT software.

Structure-based identification of potential novel inhibitors targeting FAM3B (PANDER) causing type 2 diabetes mellitus through virtual screening

Type 2 diabetes mellitus is a metabolic disorder that requires potent therapeutic approaches. The FAM3B is a cytokine-like protein also referred to as PANcreatic-DERrived factor (PANDER) which mainly exists in pancreatic islets. In the process of identifying potential inhibitors with the aid of structure-based method PANDER protein is identified as a novel therapeutic target against type 2 diabetes mellitus as it involved in the development of type 2 diabetes by negatively regulating the pancreatic β-cell function and insulin sensitivity in the liver. In the present study, the 3d model of target protein FAM3B was generated by homology modeling technique using the MODELLER9.9 program. The assessment of the structural stability of the 3d model was established by energy minimization technique. The structural quality was evaluated with standard validating protocols. Binding regions of the target protein has been determined by literature and SiteMap tool. In the current study of research, the FAM3B model was subjected to molecular screening with the Asinex-elite database of 14849 output molecules using the Glide virtual screening module in the Schrodinger suite. The final XP descriptor output of 14 molecules was analyzed and prioritized based on molecular interactions at the FAM3B active site. The docking score, binding free energies (Prime MM/GBSA) and bioavailability were undertaken into the consideration to identify lead inhibitors. The identified lead compounds were checked for ADME properties all falling within the permeable ranges. The analysis of results gave the insight to develop the novel therapeutic strategies against type 2 diabetes mellitus.

Next-Generation Noncompetitive Nanosystems Based on Gambogic Acid: In silico Identification of Transferrin Receptor Binding Sites, Regulatory Shelf Stability, and Their Preliminary Safety in Healthy Rodents

A major challenge in drug delivery is to enhance the transport of drugs across biological barriers, such as the small intestine, the blood-brain barrier, and the blood-retinal/ocular barrier, and to effectively reach the site of action while minimizing the systemic impact. In recent years, piggybacking cell surface receptors have been considered a viable strategy for active drug delivery across the biological barriers. However, the ligands used to target drugs to plasma membrane receptors often have to compete against endogenous ligands, thereby limiting their binding to the cell surface and their transport across barriers. To address this problem, gambogic acid (GA) was identified as a noncompetitive ligand specific to the transferrin receptor (TfR), a receptor present on various barriers. However, the binding sites of the GA on TfR remain unknown, an essential step toward establishing structure-activity relationships. In silico binding site prediction tools, blind docking, and molecular docking simulation confirm that the GA binding site on the TfR is independent of the transferrin-bound iron binding sites. The GA-conjugated polyesters were processed into nanoparticles suitable for drug delivery applications that possess excellent storage stability under regulatory conditions. Traditionally, GA has been used as an anticancer compound that warrants safety assessment. The preliminary studies in healthy rodents on 10-repeated oral doses show no adverse effects. This work will generate paradigm shifting, new knowledge in the field of nanomedicines using unique noncompetitive nanosystems that do not compete with endogenous transferrin.

Computational aided mechanistic understanding of Camellia sinensis bioactive compounds against co-chaperone p23 as potential anticancer agent

Co-chaperon p23 has been well established as molecular chaperon for the heat shock protein 90 (Hsp90) that further leads to immorality in cancer cells by providing defense against Hsp90 inhibitors, and as stimulating agent for generating overexpressed antiapoptotic proteins, that is, Hsp70 and Hsp27. The natural compounds such as catechins from Camellia sinensis (green tea) are also well known for inhibition activity against various cancer. However, molecular interaction profile and potential lead bioactive compounds against co-chaperon p23 from green tea are not yet reported. To this context, we study the various secondary metabolites of green tea against co-chaperon p23 using structure-based virtual screening from Traditional Chinese Medicine (TCM) database. Following 26 compounds were obtained from TCM database and further studied for extra precision molecular docking that showed binding score between -10.221 and -2.276 kcal/mol with co-chaperon p23. However, relative docking score to known inhibitors, that is, ailanthone (-4.54 kcal/mol) and gedunin ( 3.60 kcal/mol) along with ADME profile analysis concluded epicatechin (-7.013 kcal/mol) and cis-theaspirone (-4.495 kcal/mol) as potential lead inhibitors from green tea against co-chaperone p23. Furthermore, molecular dynamics simulation and molecular mechanics generalized born surface area calculations validated that epicatechin and cis-theaspirone have significantly occupied the active region of co-chaperone p23 by hydrogen and hydrophobic interactions with various residues including most substantial amino acids, that is, Thr90, Ala94, and Lys95. Hence, these results supported the fact that green tea contained potential compounds with an ability to inhibit the cancer by disrupting the co-chaperon p23 activity.

Hierarchical-Clustering, Scaffold-Mining Exercises and Dynamics Simulations for Effectual Inhibitors Against Lipid-A Biosynthesis of Helicobacter pylori

Introduction

Treatment failures of standard regimens and new strains egression are due to the augmented drug resistance conundrum. These confounding factors now became the drug designers spotlight to implement therapeutics against Helicobacter pylori strains and to safeguard infected victims with devoid of adverse drug reactions. Thereby, to navigate the chemical space for medicine, paramount vital drug target opting considerations should be imperative. The study is therefore aimed to develop potent therapeutic variants against an insightful extrapolative, common target LpxC as a follow-up to previous studies.

Methods

We explored the relationships between existing inhibitors and novel leads at the scaffold level in an appropriate conformational plasticity for lead-optimization campaign. Hierarchical-clustering and shape-based screening against an in-house library of > 21 million compounds resulted in panel of 11,000 compounds. Rigid-receptor docking through virtual screening cascade, quantum-polarized-ligand, induced-fit dockings, post-docking processes and system stability assessments were performed.

Results

After docking experiments, an enrichment performance unveiled seven ranked actives better binding efficiencies with Zinc-binding potency than substrate and in-actives (decoy-set) with ROC (1.0) and area under accumulation curve (0.90) metrics. Physics-based membrane permeability accompanied ADME/T predictions and long-range dynamic simulations of 250 ns chemical time have depicted good passive diffusion with no toxicity of leads and sustained consistency of lead1-LpxC in the physiological milieu respectively.

Conclusions

In the study, as these static outcomes obtained from this approach competed with the substrate and existing ligands in binding affinity estimations as well as positively correlated from different aspects of predictions, which could facilitate promiscuous new chemical entities against H. pylori.

Insights into the substrate binding specificity of quorum-quenching acylase PvdQ

Quorum sensing is a cell to cell signaling mechanism that enables them to coordinate their behaviors in a density-dependent manner mediated by small diffusible signaling molecules, which can control the virulence and biofilm gene expression in many Gram-negative and positive bacteria. N-acyl homoserine lactone acylase PvdQ from human opportunistic pathogen Pseudomonas aeruginosa is a quorum-quenching enzyme that can hydrolyze the amide bond of the quorum signaling N-acyl homoserine lactones (AHLs) thereby degrading the signaling molecules, turning off the biofilm phenotype and resulting in a reduction of bacterial virulence. Previous studies demonstrated that PvdQ has different preferences for N-acyl substrates with different acyl chain lengths and substituents. However, the substrate binding specificity determinants of the quorum-quenching enzyme PvdQ with the different bacterial ligands are unknown and unintuitive. Further, elucidation of these determinants can lead to mutants with efficiency and broader substrate promiscuity. To investigate this question, a computational study was carried out combining multiple molecular docking methods, molecular dynamics simulations, residue interaction network analysis, and binding free energy calculations. The main findings are: firstly, the results from pKa predictions support that the pKa of the N-terminus of Serβ1 was depressed due to the surrounding residues. Multiple molecular docking studies provide useful information about the detailed binding modes and binding affinities. Secondly, 300 ns molecular dynamics simulations were carried out to analyze the overall molecular motions of substrate-bound and substrate-free PvdQ. The specific interactions between the active site of PvdQ and different ligands revealed the determinants for the preference among the ligands. A systematic comparison and analysis of the protein dynamic fingerprint of each complex demonstrated that binding of the most favorable ligand, C12-homoserine lactone (C12-HSL), reduced the global motions of the complex and maintained the correct arrangement of the catalytic site. Further, the residue interaction network analysis of each system illustrated that there are more communication contacts and pathways between the residues in the C12-HSL complex as compared to complexes with the other ligands. The binding of the C12-HSL ligand facilitates structural communication between the two knobs and the active site. While the binding of the other ligands tend to impair specific communication pathways between the two knobs and the active site, and lead to a catalytically inefficient state. Finally, simulation results from free energy landscape and binding free energy analysis revealed that the C12-HSL ligand has the lowest binding free energy and greater stability than the less favored ligands. Each of the following residues: Serβ1, Hisβ23, Pheβ24, Metβ30, Pheβ32, Leuβ50, Asnβ57, Thrβ69, Valβ70, Trpβ162, Trpβ186, Asnβ269, Argβ297 and Leuα146, play different roles in substrate binding specificity. This is the first computational study that provides molecular information for structure-dynamic-function relationships of PvdQ with different ligands and demonstrates determinants of bacterial substrate binding specificity.

Structure-guided approach identifies a novel class of HIV-1 ribonuclease H inhibitors: binding mode insights through magnesium complexation and site-directed mutagenesis studies

Persistent HIV infection requires lifelong treatment and among the 2.1 million new HIV infections that occur every year there is an increased rate of transmitted drug-resistant mutations. This fact requires a constant and timely effort in order to identify and develop new HIV inhibitors with innovative mechanisms. The HIV-1 reverse transcriptase (RT) associated ribonuclease H (RNase H) is the only viral encoded enzyme that still lacks an efficient inhibitor despite the fact that it is a well-validated target whose functional abrogation compromises viral infectivity. Identification of new drugs is a long and expensive process that can be speeded up by in silico methods. In the present study, a structure-guided screening is coupled with a similarity-based search on the Specs database to identify a new class of HIV-1 RNase H inhibitors. Out of the 45 compounds selected for experimental testing, 15 inhibited the RNase H function below 100 μM with three hits exhibiting IC50 values <10 μM. The most active compound, AA, inhibits HIV-1 RNase H with an IC50 of 5.1 μM and exhibits a Mg-independent mode of inhibition. Site-directed mutagenesis studies provide valuable insight into the binding mode of newly identified compounds; for instance, compound AA involves extensive interactions with a lipophilic pocket formed by Ala502, Lys503, and Trp (406, 426 and 535) and polar interactions with Arg557 and the highly conserved RNase H primer-grip residue Asn474. The structural insights obtained from this work provide the bases for further lead optimization.

Vitamin D receptor (VDR) non-synonymous single nucleotide polymorphisms (nsSNPs) affect the calcitriol drug response - A theoretical insight

Pharmacogenetics and pharmacogenomics have become presumptive with advancements in next-generation sequencing technology. In complex diseases, distinguishing the feasibility of pathogenic and neutral disease-causing variants is a time consuming and expensive process. Recent drug research and development processes mainly rely on the relationship between the genotype and phenotype through Single nucleotide polymorphisms (SNPs). The SNPs play an indispensable role in elucidating the individual's vulnerability to disease and drug response. The understanding of the interplay between these leads to the establishment of personalized medicine. In order to address this issue, we developed a computational pipeline of vitamin D receptor (VDR) for SNP centered study by application of elegant molecular docking and molecular dynamics simulation approaches. In a few SNPs the volume of the binding cavities has increased in mutant structures when compared to the wild type, indicating a weakening in interaction (699.1 ų in wild type Vs. 738.8 in Leu230Val, 820.7 ų in Arg247Leu). This also differently reflected in the H-bond interactions and binding free energies -169.93 kcal/mol (wild type) Vs -156.43 kcal/mol (R154W), -105.49 kcal/mol (R274L) in Leu230Val and Arg247Leu respectively. Although we could not find noteworthy changes in the binding free energies and binding pocket in the remaining mutations, the H-bond interactions made these SNPs deleterious. Thus, we further analyzed the H-bond interactions and distances using molecular dynamics (MD) simulation studies.

Revealing vilazodone's binding mechanism underlying its partial agonism to the 5-HT1A receptor in the treatment of major depressive disorder

It has been estimated that major depressive disorder (MDD) will become the second largest global burden among all diseases by 2030. Various types of drugs, including selective serotonin reuptake inhibitors (SSRIs), serotonin-norepinephrine reuptake inhibitors (SNRIs), and serotonin receptor partial agonist/reuptake inhibitors (SPARIs), have been approved and become the primary or first-line medications prescribed for MDD. SPARI was expected to demonstrate more enhanced drug efficacy and a rapid onset of action as compared to SSRI and SNRI. As one of the most famous SPARIs, vilazodone was approved by the FDA for the treatment of MDD. Because of the great clinical importance of vilazodone, its binding mechanism underlying its partial agonism to the 5-HT_1A receptor (5-HT_1AR) could provide valuable information to SPARIs' drug-like properties. However, this mechanism has not been reported to date; consequently, the rational design of new efficacious SPARI-based MDD drugs is severely hampered. To explore the molecular mechanism of vilazodone, an integrated computational strategy was adopted in this study to reveal its binding mechanism and prospective structural feature at the agonist binding site of 5-HT_1AR. As a result, 22 residues of this receptor were identified as hotspots, consistently favoring the binding of vilazodone and its analogues, and a common binding mechanism underlying their partial agonism to 5-HT_1AR was, therefore, discovered. Moreover, three main interaction features between vilazodone and 5-HT_1AR have been revealed and schematically summarized. In summary, this newly identified binding mechanism will provide valuable information for medicinal chemists working in the field of rational design of novel SPARIs for MDD treatment.

A Mitochondrial-targeted purine-based HSP90 antagonist for leukemia therapy

Reprogramming of mitochondrial functions sustains tumor growth and may provide therapeutic opportunities. Here, we targeted the protein folding environment in mitochondria by coupling a purine-based inhibitor of the molecular chaperone Heat Shock Protein-90 (Hsp90), PU-H71 to the mitochondrial-targeting moiety, triphenylphosphonium (TPP). Binding of PU-H71-TPP to ADP-Hsp90, Hsp90 co-chaperone complex or mitochondrial Hsp90 homolog, TRAP1 involved hydrogen bonds, π-π stacking, cation-π contacts and hydrophobic interactions with the surrounding amino acids in the active site. PU-H71-TPP selectively accumulated in mitochondria of tumor cells (17-fold increase in mitochondria/cytosol ratio), whereas unmodified PU-H71 showed minimal mitochondrial localization. Treatment of tumor cells with PU-H71-TPP dissipated mitochondrial membrane potential, inhibited oxidative phosphorylation in sensitive cell types, and reduced ATP production, resulting in apoptosis and tumor cell killing. Unmodified PU-H71 had no effect. Bioinformatics analysis identified a “mitochondrial Hsp90” signature in Acute Myeloid Leukemia (AML), which correlates with worse disease outcome. Accordingly, inhibition of mitochondrial Hsp90s killed primary and cultured AML cells, with minimal effects on normal peripheral blood mononuclear cells. These data demonstrate that directing Hsp90 inhibitors with different chemical scaffolds to mitochondria is feasible and confers improved anticancer activity. A potential “addiction” to mitochondrial Hsp90s may provide a new therapeutic target in AML.

Exploring putative inhibitors of Death Associated Protein Kinase 1 (DAPK1) via targeting Gly-Glu-Leu (GEL) and Pro-Glu-Asn (PEN) substrate recognition motifs

Recently, a new signaling complex Death Associated Protein Kinase 1 (DAPK1) ̶ N-methyl-D-aspartate receptor subtype 2B (NMDAR2B or NR2B) engaged in the neuronal death cascade was identified and it was found that after stroke injury, N-methyl-D-aspartate glutamate (NMDA) receptors interact with DAPK1 through NR2B subunit and lead to excitotoxicity via over-activation of NMDA receptors. An acute brain injury, such as stroke, is a serious life-threatening medical condition which occurs due to poor blood supply to the brain and further leads to neuronal cell death. During a stroke, activated DAPK1 migrates towards the extra-synaptic site and binds to NR2B subunit of NMDA receptor. It is this DAPK1-NR2B interaction that arbitrates the pathological processes like apoptosis, necrosis, and autophagy of neuronal cells observed in stroke injury, hence we aimed to inhibit this vital interaction to prevent neuronal damage. In the present study, using PubChem database, we applied an integrative approach of virtual screening and molecular dynamic simulations and identified a potential lead compound 11 that interrupts DAPK1-NR2B interaction by competing with both ATP and substrate for their binding sites on DAPK1. This inhibitor was found potent and considerably selective to DAPK1 as it made direct contact with the ATP binding sites as well as substrate recognition motifs: Gly-Glu-Leu (GEL) and Pro-Glu-Asn (PEN). Further in vitro and in vivo experiments are demanded to validate the efficacy of compound 11 nevertheless, it can be considered as suitable starting point for designing DAPK1 inhibitors.

Pharmacoinformatic Study on the Selective Inhibition of the Protozoan Dihydrofolate Reductase Enzymes

Dihydrofolate reductase (DHFR) is an essential enzyme of the folate metabolic pathway in protozoa and it is a validated, potential drug target in many infectious diseases. Information about unique conserved residues of the DHFR enzyme is required to understand residual selectivity of the protozoan DHFR enzyme. The three dimensional crystal structures are not available for all the protozoan DHFR enzymes. Enzyme-substrate/inhibitor interaction information is required for the binding mode characterization in protozoan DHFR for selective inhibitor design. In this work, multiple sequence analysis was carried out in all the studied species. Homology models were built for protozoan DHFR enzymes, for which 3D structures are not available in PDB. The molecular docking and Prime-MMGBSA calculations of the natural substrate (dihydrofolate, DHF) and classical DHFR inhibitor (methotrexate, MTX) were performed in protozoan DHFR enzymes. Comparative sequence analysis showed that an overall sequence identity between the studied species ranging from 22.94 % (CfDHFR-BgDHFR) to 94.61 % (LdDHFR-LmDHFR). Interestingly, it was observed that most of the active site residues were conserved in all the cases and all the enzymes exhibit similar key binding interactions with DHF and MTX in molecular docking analysis, but there are a few key binding residues which differ in protozoan species that makes it suitable for target selectivity. This information can be used to design selective and potent protozoan DHFR enzyme inhibitors.

Discovering key residues of dengue virus NS2b-NS3-protease: New binding sites for antiviral inhibitors design

The NS2B-NS3 protease is essential for the Dengue Virus (DENV) replication process. This complex constitutes a target for efficient antiviral discovery because a drug could inhibit the viral polyprotein processing. Furthermore, since the protease is highly conserved between the four Dengue virus serotypes, it is probable that a drug would be equally effective against all of them. In this article, a strategy is reported that allowed us to identify influential residues on the function of the Dengue NS2b-NS3 Protease. Moreover, this is a strategy that could be applied to virtually any protein for the search of alternative influential residues, and for non-competitive inhibitor development. First, we incorporated several features derived from computational alanine scanning mutagenesis, sequence, structure conservation, and other structure-based characteristics. Second, these features were used as variables to obtain a multilayer perceptron model to identify defined groups (clusters) of key residues as possible candidate pockets for binding sites of new leads on the DENV protease. The identified residues included: i) amino acids close to the beta sheet-loop-beta sheet known to be important in its closed conformation for NS2b ii) residues close to the active site, iii) several residues evenly spread on the NS2b-NS3 contact surface, and iv) some inner residues most likely related to the overall stability of the protease. In addition, we found concordance on our list of residues with previously identified amino acids part of a highly conserved peptide studied for vaccine development.

The evaluation of QM/MM-driven molecular docking combined with MM/GBSA calculations as a halogen-bond scoring strategy

The combination of quantum mechanics/molecular mechanics-driven (QM/MM) molecular docking with binding free-energy calculations was successfully used to reproduce the X-ray geometries of protein–ligand complexes with halogen bonding. The procedure involves quantum-polarized ligand docking (QPLD) to obtain the QM-derived ligand atomic charges in the protein environment at the B3PW91/cc-pVTZ level and the MM/GBSA (generalized-Born/surface area) algorithm to calculate the binding free energies of resultant complexes. The performance was validated using a set of 106 X-ray complexes and compared with the Glide and AutoDock VinaXB scoring functions in terms of RMSD and the reconstruction of halogen-bond geometry (distance and σ-hole angle). The results revealed that docking and scoring using the QPLD–GBSA procedure outperformed the remaining scoring functions in the majority of instances. Additionally, a comparison of the orientation of the top ranked binding poses calculated using the fixed atomic charges of ligands obtained from force-field parameterization and by QM calculations in the protein environment provides strong evidence that the use of QM-derived charges is significant.

Structure and sequence based functional annotation of Zika virus NS2b protein: Computational insights

While Zika virus (ZIKV) outbreaks are a growing concern for global health, a deep understanding about the virus is lacking. Here we report a contribution to the basic science on the virus- a detailed computational analysis of the non structural protein NS2b. This protein acts as a cofactor for the NS3 protease (NS3Pro) domain that is important on the viral life cycle, and is an interesting target for drug development. We found that ZIKV NS2b cofactor is highly similar to other virus within the Flavivirus genus, especially to West Nile Virus, suggesting that it is completely necessary for the protease complex activity. Furthermore, the ZIKV NS2b has an important role to the function and stability of the ZIKV NS3 protease domain even when presents a low conservation score. In addition, ZIKV NS2b is mostly rigid, which could imply a non dynamic nature in substrate recognition. Finally, by performing a computational alanine scanning mutagenesis, we found that residues Gly 52 and Asp 83 in the NS2b could be important in substrate recognition.

A highly efficient sorbitol dehydrogenase from Gluconobacter oxydans G624 and improvement of its stability through immobilization

A sorbitol dehydrogenase (GoSLDH) from Gluconobacter oxydans G624 (G. oxydans G624) was expressed in Escherichia coli BL21(DE3)-CodonPlus RIL. The complete 1455-bp codon-optimized gene was amplified, expressed, and thoroughly characterized for the first time. GoSLDH exhibited K_m and k_cat values of 38.9 mM and 3820 s⁻¹ toward L-sorbitol, respectively. The enzyme exhibited high preference for NADP⁺ (vs. only 2.5% relative activity with NAD⁺). GoSLDH sequencing, structure analyses, and biochemical studies, suggested that it belongs to the NADP⁺-dependent polyol-specific long-chain sorbitol dehydrogenase family. GoSLDH is the first fully characterized SLDH to date, and it is distinguished from other L-sorbose-producing enzymes by its high activity and substrate specificity. Isothermal titration calorimetry showed that the protein binds more strongly to D-sorbitol than other L-sorbose-producing enzymes, and substrate docking analysis confirmed a higher turnover rate. The high oxidation potential of GoSLDH for D-sorbitol was confirmed by cyclovoltametric analysis. Further, stability of GoSLDH significantly improved (up to 13.6-fold) after cross-linking of immobilized enzyme on silica nanoparticles and retained 62.8% residual activity after 10 cycles of reuse. Therefore, immobilized GoSLDH may be useful for L-sorbose production from D-sorbitol.

Molecular Characterization and Enhancement of Anticancer Activity of Caffeic Acid Phenethyl Ester by γ Cyclodextrin

Caffeic Acid Phenethyl Ester (CAPE) is a key component in New Zealand propolis, known for a variety of health promoting and therapeutic potentials. We investigated the molecular mechanism of anticancer and anti-metastasis activities of CAPE. cDNA array performed on the control and CAPE-treated breast cancer cells revealed activation of DNA damage signaling involving upregulation of GADD45α and p53 tumor suppressor proteins. Molecular docking analysis revealed that CAPE is capable of disrupting mortalin-p53 complexes. We provide experimental evidence and demonstrate that CAPE induced disruption of mortalin-p53 complexes led to nuclear translocation and activation of p53 resulting in growth arrest in cancer cells. Furthermore, CAPE-treated cells exhibited downregulation of mortalin and several other key regulators of cell migration accountable for its anti-metastasis activity. Of note, we found that whereas CAPE was unstable in the culture medium (as it gets degraded into caffeic acid by secreted esterases), its complex with gamma cyclodextrin (γCD) showed high efficacy in anti-tumor and anti-metastasis assays in vitro and in vivo (when administered through either intraperitoneal or oral route). The data proposes that CAPE-γCD complex is a potent anti-cancer and anti-metastasis reagent.

Insights into mechanism of pyrido[2,3-d]pyrimidines as DYRK1A inhibitors based on molecular dynamic simulations

DYRK1A is characterized by the early development and regulation of neuronal proliferation, and its over expression gives rise to neurological abnormalities. As the promising DYRK1A inhibitors, the binding mechanism between DYRK1A and pyrido[2,3-d]pyrimidines derivatives at molecular level are still veiled. In this article, it was achieved to get the structural insights into pyrido[2,3-d]pyrimidines derivatives as DYRK1A inhibitors by means of comprehensive computational approaches involving molecular docking, molecular dynamics simulation, free energy calculation, and energy decomposition analysis. The calculated energy values were highly consistent with the experimental activities. Based on the individual energy terms analysis, the van der Waals interaction was the major leading force in the DYRK1A-ligand interaction. Lys188 was the important residue that formed the hydrogen bond, which improved the inhibitory activity. Furthermore, four novel inhibitors with higher predicted activity were designed based on the obtained findings and confirmed by molecular simulations. Our study is expected to provide significant drug design strategy for the development of more promising DYRK1A inhibitors. Proteins 2016; 84:1108-1123. © 2016 Wiley Periodicals, Inc.

A computational study on role of 6-(hydroxymethyl)-3-[3,4,5-trihydroxy-6-[(3,4,5-trihydroxyoxan-2-yl)oxymethyl]oxan-2-yl]oxyoxane-2,4,5-triol in the regulation of blood glucose level

6-(hydroxymethyl)-3-[3,4,5-trihydroxy-6-[(3,4,5-trihydroxyoxan-2-yl)oxymethyl]oxan-2-yl]oxyoxane-2,4,5-triol (SID 242078875) was isolated from the fruits of Syzygium densiflorum Wall. ex Wight & Arn (Myrtaceae), which has been traditionally used in the treatment of diabetes by the tribes of The Nilgiris, Tamil Nadu, India. In this study, reverse pharmacophore mapping approach and text-based database search identified the dipeptidyl peptidase-IV, protein-tyrosine phosphatase 1B, phosphoenolpyruvate carboxykinase, glycogen synthase kinase-3β and glucokinase as potential targets of SID 242078875 in diabetes management. Further, molecular docking was performed to predict the binding pose of SID 242078875 in the active site region of the target protein. In addition, dynamic behaviour and stability of protein-ligand complexes were observed for a period of 50 ns through molecular dynamics simulation.

The Binding Mode Prediction and Similar Ligand Potency in the Active Site of Vitamin D Receptor with QM/MM Interaction, MESP, and MD Simulation

Non-secosteroidal ligands are well-known vitamin D receptor (VDR) agonists. In this study, we described a combined QM/MM to define the protein-ligand interaction energy a strong positive correlation in both QM-MM interaction energy and binding free energy against the biological activity. The molecular dynamics simulation study was performed, and specific interactions were extensively studied. The molecular docking results and surface analysis shed light on steric and electrostatic complementarities of these non-secosteroidal ligands to VDR. Finally, the drug likeness properties were also calculated and found within the acceptable range. The results show that bulky group substitutions in side chain decrease the VDR activity, whereas a small substitution increased it. Functional analyses of H393A and H301A mutations substantiate their roles in the VDR agonistic and antagonistic activities. Apart from the His393 and His301, two other amino acids in the hinge region viz. Ser233 and Arg270 acted as an electron donor/acceptor specific to the agonist in the distinct ligand potency. The results from this study disclose the binding mechanism of VDR agonists and structural modifications required to improve the selectivity.

Inhibitor design against JNK1 through e-pharmacophore modeling docking and molecular dynamics simulations

c-Jun-NH2 terminal kinases (JNKs) come under a class of serine/threonine protein kinases and are encoded by three genes, namely JNK1, JNK2 and JNK3. Human JNK1 is a cytosolic kinase belonging to mitogen-activated protein kinase (MAPK) family, which plays a major role in intracrinal signal transduction cascade mechanism. Overexpressed human JNK1, a key kinase interacts with other kinases involved in the etiology of many cancers, such as skin cancer, liver cancer, breast cancer, brain tumors, leukemia, multiple myeloma and lymphoma. Thus, to unveil a novel human JNK1 antagonist, receptor-based pharmacophore modeling was performed with the available eighteen cocrystal structures of JNK1 in the protein data bank. Eighteen e-pharmacophores were generated from the 18 cocrystal structures. Four common e-pharmacophores were developed from the 18 e-pharmacophores, which were used as three-dimensional (3D) query for shape-based similarity screening against more than one million small molecules to generate a JNK1 ligand library. Rigid receptor docking (RRD) performed using GLIDE v6.3 for the 1683 compounds from in-house library and 18 cocrystal ligands with human JNK1 from lower stringency to higher stringency revealed 17 leads. Further to derive the best leads, dock complexes obtained from RRD were studied further with quantum-polarized ligand docking (QPLD), induced fit docking (IFD) and molecular mechanics/generalized Born surface area (MM-GBSA). Four leads have showed lesser binding free energy and better binding affinity towards JNK1 compared to 18 cocrystal ligands. Additionally, JNK1-lead1 complex interaction stability was reasserted using 50 ns MD simulations run and also compared with the best resolute cocrystal structure using Desmond v3.8. Thus, the results obtained from RRD, QPLD, IFD and MD simulations indicated that lead1 might be used as a potent antagonist toward human JNK1 in cancer therapeutics.

Structural insights on identification of potential lead compounds targeting WbpP in Vibrio vulnificus through structure-based approaches

WbpP encoding UDP-GlcNAC C4 epimerase is responsible for the activation of virulence factor in marine pathogen Vibrio vulnificus (V. vulnificus) and it is linked to many aquatic diseases, thus making it a potential therapeutic target. There are few reported compounds that include several natural products and synthetic compounds targeting Vibrio sp, but specific inhibitor targeting WbpP are unavailable. Here, we performed structure-based virtual screening using chemical libraries such as Binding, TOSLab and Maybridge to identify small molecule inhibitors of WbpP with better drug-like properties. Deficient structural information forced to model the structure and the stable protein structure was obtained through 30 ns of MD simulations. Druggability regions are focused for new lead compounds and our screening protocol provides fast docking of entire small molecule library with screening criteria of ADME/Lipinski filter/Docking followed by re-docking of top hits using a method that incorporates both ligand and protein flexibility. Docking conformations of lead molecules interface displays strong H-bond interactions with the key residues Gly101, Ser102, Val195, Tyr165, Arg298, Val209, Ser142, Arg233 and Gln200. Subsequently, the top-ranking compounds were prioritized using the molecular dynamics simulation-based conformation and stability studies. Our study suggests that the proposed compounds may aid as a starting point for the rational design of novel therapeutic agents.

E-pharmacophore-based virtual screening to identify GSK-3β inhibitors

Glycogen synthase kinase-3β (GSK-3β) is a serine/threonine kinase which has attracted significant attention during recent years in drug design studies. The deregulation of GSK-3β increased the loss of hippocampal neurons by triggering apoptosis-mediating production of neurofibrillary tangles and alleviates memory deficits in Alzheimer's disease (AD). Given its role in the formation of neurofibrillary tangles leading to AD, it has been a major therapeutic target for intervention in AD, hence was targeted in the present study. Twenty crystal structures were refined to generate pharmacophore models based on energy involvement in binding co-crystal ligands. Four common e-pharmacophore models were optimized from the 20 pharmacophore models. Shape-based screening of four e-pharmacophore models against nine established small molecule databases using Phase v3.9 had resulted in 1800 compounds having similar pharmacophore features. Rigid receptor docking (RRD) was performed for 1800 compounds and 20 co-crystal ligands with GSK-3β to generate dock complexes. Interactions of the best scoring lead obtained through RRD were further studied with quantum polarized ligand docking (QPLD), induced fit docking (IFD) and molecular mechanics/generalized Born surface area. Comparing the obtained leads to 20 co-crystal ligands resulted in 18 leads among them, lead1 had the lowest docking score, lower binding free energy and better binding orientation toward GSK-3β. The 50 ns MD simulations run confirmed the stable nature of GSK-3β-lead1 docking complex. The results from RRD, QPLD, IFD and MD simulations confirmed that lead1 might be used as a potent antagonist for GSK-3β.

Computational Structure-Based De Novo Design of Hypothetical Inhibitors against the Anti- Inflammatory Target COX-2

Cyclooxygenase-2 (COX-2) produces prostaglandins in inflamed tissues and hence has been considered as an important target for the development of anti-inflammatory drugs since long. Administration of traditional non-steroidal anti-inflammatory drugs (NSAIDs) and other COX-2 selective inhibitors (COXIBS) for the treat of inflammation has been found to be associated with side effects, which mainly includes gastro-intestinal (GI) toxicity. The present study involves developing a virtual library of novel molecules with high druglikeliness using structure-based de novo drug designing and 2D fingerprinting approach. A library of 2657 drug like molecules was generated. 2D fingerprinting based screening of the designed library gave a unique set of compounds. Molecular docking approach was then used to identify two compounds highly specific for COX-2 isoform. Molecular dynamics simulations of protein-ligand complexes revealed that the candidate ligands were dynamically stable within the cyclooxygenase binding site of COX-2. The ligands were further analyzed for their druglikeliness, ADMET properties and synthetic accessibility using knowledge based set of rules. The results revealed that the molecules are predicted to selectively bind to COX-2 enzyme thereby potentially overcoming the limitations posed by the drugs in clinical use.

An SH2 domain model of STAT5 in complex with phospho-peptides define "STAT5 Binding Signatures"

The signal transducer and activator of transcription 5 (STAT5) is a member of the STAT family of proteins, implicated in cell growth and differentiation. STAT activation is regulated by phosphorylation of protein monomers at conserved tyrosine residues, followed by binding to phospho-peptide pockets and subsequent dimerization. STAT5 is implicated in the development of severe pathological conditions, including many cancer forms. However, nowadays a few STAT5 inhibitors are known, and only one crystal structure of the inactive STAT5 dimer is publicly available. With a view to enabling structure-based drug design, we have: (1) analyzed phospho-peptide binding pockets on SH2 domains of STAT5, STAT1 and STAT3; (2) generated a model of STAT5 bound to phospho-peptides; (3) assessed our model by docking against a class of known STAT5 inhibitors (Müller et al. in ChemBioChem 9:723-727, 2008); (4) used molecular dynamics simulations to optimize the molecular determinants responsible for binding and (5) proposed unique "Binding Signatures" of STAT5. Our results put in place the foundations to address STAT5 as a target for rational drug design, from sequence, structural and functional perspectives.

Mechanistic insights of SrtA–LPXTG blockers targeting the transpeptidase mechanism in Streptococcus mutans

The SrtA–LPXTG interaction plays a key role in transpeptidation reaction, cell wall and biofilm formations. This study explains the blocking of LEU interactions with SrtA will results as SrtA inhibitors through MD simulation and energy calculations methods.

Protease activated receptor-2 (PAR2): possible target of phytochemicals

The use of phytochemicals either singly or in combination with other anticancer drugs comes with an advantage of less toxicity and minimal side effects. Signaling pathways play central role in cell cycle, cell growth, metabolism, etc. Thus, the identification of phytochemicals with promising antagonistic effect on the receptor/s playing key role in single transduction may have better therapeutic application. With this background, phytochemicals were screened against protease-activated receptor 2 (PAR2). PAR2 belongs to the superfamily of GPCRs and is an important target for breast cancer. Using in silico methods, this study was able to identify the phytochemicals with promising binding affinity suggesting their therapeutic potential in the treatment of breast cancer. The findings from this study acquires importance as the information on the possible agonists and antagonists of PAR2 is limited due its unique mechanism of activation.

Exploring the prominent performance of CX-4945 derivatives as protein kinase CK2 inhibitors by a combined computational study

Protein kinase CK2, also known as casein kinase II, is related to various cellular events and is a potential target for numerous cancers. In this study, we attempted to gain more insight into the inhibition process of CK2 by a series of CX-4945 derivatives through an integrated computational study that combines molecular docking, molecular dynamics (MD) simulations, and binding free energy calculations. Based on the binding poses predicted by molecular docking, the MD simulations were performed to explore the dynamic binding processes for ten selected inhibitors. Then, both Molecular Mechanics/Poisson Boltzmann Surface Area (MM/PBSA) and Molecular Mechanics/Generalized Born Surface Area (MM/GBSA) techniques were employed to predict the binding affinities of the studied systems. The predicted binding energies of the selected inhibitors correlate well with their experimental activities (r(2) = 0.78). The van der Waals term is the most favorable component for the total energies. The free energy decomposition on a per residue basis reveals that the residue K68 is essential for the electrostatic interactions between CK2 and the studied inhibitors and numerous residues, including L45, V53, V66, F113, M163 and I174, play critical roles in forming van der Waals interactions with the inhibitors. Finally, a number of new derivatives were designed and the binding affinity and the predicted binding free energies of each designed molecule were obtained on the basis of molecular docking and MM/PBSA. It is expected that our research will benefit the future rational design of novel and potent inhibitors of CK2.

Computational studies on the binding mechanism between triazolone inhibitors and Chk1 by molecular docking and molecular dynamics

Chk1, a serine/threonine protein kinase that participates in transducing DNA damage signals, is an attractive target due to its involvement in tumor initiation and progression. As a novel Chk1 inhibitor, the triazolone's bioactivity mechanism is not clear. In this study, we carried out an integrated computational study that combines molecular docking, molecular dynamics (MD) simulations, and binding free energy calculations to identify the key factors necessary for the bioactivities. With the aim of discerning the structural features that affect the inhibitory activity of triazolones, MK-8776, a Chk1 inhibitor that reached the clinical stage, was also used as a reference for simulations. A comparative analysis of the triazolone inhibitors at the molecular level offers valuable insight into the structural and energetic properties. A general feature is that all the studied inhibitors bind in the pocket characterized by residues Leu14, Val22, Ala35, Glu84, Tyr85, Cys86, and Leu136 of Chk1. Moreover, introducing hydrophobic groups into triazolone inhibitors is favorable for binding to Chk1, which is corroborated by residue Leu136 with a relatively large difference in the contribution between MK-8776 and five triazolones to the total binding free energies. A hydrogen bond between the polar hydrogen atoms at R1 and Cys86 can facilitate proper placement of the inhibitor in the binding pocket of Chk1 that favors binding. However, the introduction of hydrophilic groups into the R2 position diminishes binding affinity. The information provided by this research is of benefit for further rational design of novel promising inhibitors of Chk1.

Interaction investigations of crustacean β-GBP recognition toward pathogenic microbial cell membrane and stimulate upon prophenoloxidase activation

In invertebrates, crustaceans' immune system consists of pattern recognition receptors (PRRs) instead of immunoglobulin's, which involves in the microbial recognition and initiates the protein-ligand interaction between hosts and pathogens. In the present study, PRRs namely β-1,3 glucan binding protein (β-GBP) from mangrove crab Episesarma tetragonum and its interactions with the pathogens such as bacterial and fungal outer membrane proteins (OMP) were investigated through microbial aggregation and computational interaction studies. Molecular recognition and microbial aggregation results of Episesarma tetragonum β-GBP showed the specific binding affinity toward the fungal β-1,3 glucan molecule when compared to other bacterial ligands. Because of this microbial recognition, prophenoloxidase activity was enhanced and triggers the innate immunity inside the host animal. Our findings disclose the role of β-GBP in molecular recognition, host-pathogen interaction through microbial aggregation, and docking analysis. In vitro results were concurred with the in silico docking, and molecular dynamics simulation analysis. This study would be helpful to understand the molecular mechanism of β-GBP and update the current knowledge on the PRRs of crustaceans.

Identification of natural compound inhibitors for multidrug efflux pumps of Escherichia coli and Pseudomonas aeruginosa using in silico high-throughput virtual screening and in vitro validation

Pseudomonas aeruginosa and Escherichia coli are resistant to wide range of antibiotics rendering the treatment of infections very difficult. A main mechanism attributed to the resistance is the function of efflux pumps. MexAB-OprM and AcrAB-TolC are the tripartite efflux pump assemblies, responsible for multidrug resistance in P. aeruginosa and E. coli respectively. Substrates that are more susceptible for efflux are predicted to have a common pharmacophore feature map. In this study, a new criterion of excluding compounds with efflux substrate-like features was used, thereby refining the selection process and enriching the inhibitor identification process. An in-house database of phytochemicals was created and screened using high-throughput virtual screening against AcrB and MexB proteins and filtered by matching with the common pharmacophore models (AADHR, ADHNR, AAHNR, AADHN, AADNR, AAADN, AAADR, AAANR, AAAHN, AAADD and AAADH) generated using known efflux substrates. Phytochemical hits that matched with any one or more of the efflux substrate models were excluded from the study. Hits that do not have features similar to the efflux substrate models were docked using XP docking against the AcrB and MexB proteins. The best hits of the XP docking were validated by checkerboard synergy assay and ethidium bromide accumulation assay for their efflux inhibition potency. Lanatoside C and diadzein were filtered based on the synergistic potential and validated for their efflux inhibition potency using ethidium bromide accumulation study. These compounds exhibited the ability to increase the accumulation of ethidium bromide inside the bacterial cell as evidenced by these increase in fluorescence in the presence of the compounds. With this good correlation between in silico screening and positive efflux inhibitory activity in vitro, the two compounds, lanatoside C and diadzein could be promising efflux pump inhibitors and effective to use in combination therapy against drug resistant strains of P. aeruginosa and E. coli.