Combining docking with pharmacophore filtering for improved virtual screening

In-vitro Evaluation of Triazine Scaffold for Anticancer Drug Development: A Review

INTRODUCTION

The widespread importance of the synthesis and modification of anticancer agents has given rise to many numbers of medicinal chemistry programs. In this regard, triazine derivatives have attracted attention due to their remarkable activity against a wide range of cancer cells. This evaluation covers work reports to define the anticancer activity, the most active synthesized compound for the target, the SAR and, when described, the probable MOA besides similarly considered to deliver complete and target-pointed data for the development of types of anti-tumour medicines of triazine derivatives. Triazine scaffold for the development of anticancer analogues. Triazine can also relate to numerous beneficial targets, and their analogues have auspicious in-vitro and in-vivo anti-tumour activity. Fused molecules can improve efficacy, and drug resistance and diminish side effects, and numerous hybrid molecules are beneath diverse stages of clinical trials, so hybrid derivatives of triazine may offer valuable therapeutic involvement for the dealing of tumours.

OBJECTIVE

The objective of the recent review was to summarize the recent reports on triazine as well as its analogues with respect to its anticancer therapeutic potential.

CONCLUSION

The content of the review would be helpful to update the researchers working towards the synthesis and designing of new molecules for the treatment of various types of cancer disease with the recent molecules that have been produced from the triazine scaffold. Triazine scaffolds based on 1,3,5-triazine considerably boost molecular diversity levels and enable covering chemical space in key medicinal chemistry fields.

Bioactivity of dihydropyrimidinone derivatives as inhibitors of cyclooxygenase-2 (COX-2): an in silico approach

Cyclooxygenase-2 (COX-2) is an enzyme involved in inflammation. The overexpression of COX-2 causes chronic inflammation, which can be prevented by COX-2 inhibitors. Generally, COX-2 inhibitors possess a carboxyl group and an aromatic ring in their molecular structure. These moieties are involved in the interaction with the active site of COX-2, thus playing a pivotal role in the inhibitory activity. Regarding the requisite molecular structure of COX-2 inhibitors, derivatives of dihydropyrimidinone (DHPM) are ideal candidates to be explored as COX-2 inhibitors, due to the ease of synthesis and their versatility to be transformed chemically. In this study, we prepared a novel small library consisting of 288 designed DHPM derivatives by varying the constituent components. The selection criteria of potential candidates for the COX-2 inhibitor of the data bank involve in silico studies via molecular docking investigations, prediction of ADMET and druglikeness, as well as molecular dynamics (MD) simulations. Molecular docking served as the initial step of selection, based on the comparison of grid score, docking pose, and interactions with those of lumiracoxib (LUR) as the original ligand of COX-2. The next criteria of selection were scores obtained from the ADMET and druglikeness by comparing the designed candidates with COX-2 inhibitors that were already marketed. Compound RDUE2 and SDT29 were the most potential candidates, which were further analyzed using the MD simulation. The results of the MD simulation indicated that RDUE2 and SDT29 interacted stably with amino acid residues on the active site of COX-2. The estimation of binding free energy indicated that SDT29 exhibited an inhibitory activity comparable to that of LUR, whereas RDUE2 showed a lower inhibitory activity than that of SDT29 and LUR.

Past, Present, and Future Perspectives on Computer-Aided Drug Design Methodologies

The application of computational approaches in drug discovery has been consolidated in the last decades. These families of techniques are usually grouped under the common name of "computer-aided drug design" (CADD), and they now constitute one of the pillars in the pharmaceutical discovery pipelines in many academic and industrial environments. Their implementation has been demonstrated to tremendously improve the speed of the early discovery steps, allowing for the proficient and rational choice of proper compounds for a desired therapeutic need among the extreme vastness of the drug-like chemical space. Moreover, the application of CADD approaches allows the rationalization of biochemical and interactive processes of pharmaceutical interest at the molecular level. Because of this, computational tools are now extensively used also in the field of rational 3D design and optimization of chemical entities starting from the structural information of the targets, which can be experimentally resolved or can also be obtained with other computer-based techniques. In this work, we revised the state-of-the-art computer-aided drug design methods, focusing on their application in different scenarios of pharmaceutical and biological interest, not only highlighting their great potential and their benefits, but also discussing their actual limitations and eventual weaknesses. This work can be considered a brief overview of computational methods for drug discovery.

Comparative Assessment of Docking Programs for Docking and Virtual Screening of Ribosomal Oxazolidinone Antibacterial Agents

Oxazolidinones are a broad-spectrum class of synthetic antibiotics that bind to the 50S ribosomal subunit of Gram-positive and Gram-negative bacteria. Many crystal structures of the ribosomes with oxazolidinone ligands have been reported in the literature, facilitating structure-based design using methods such as molecular docking. It would be of great interest to know in advance how well docking methods can reproduce the correct ligand binding modes and rank these correctly. We examined the performance of five molecular docking programs (AutoDock 4, AutoDock Vina, DOCK 6, rDock, and RLDock) for their ability to model ribosomal–ligand interactions with oxazolidinones. Eleven ribosomal crystal structures with oxazolidinones as the ligands were docked. The accuracy was evaluated by calculating the docked complexes’ root-mean-square deviation (RMSD) and the program’s internal scoring function. The rankings for each program based on the median RMSD between the native and predicted were DOCK 6 > AD4 > Vina > RDOCK >> RLDOCK. Results demonstrate that the top-performing program, DOCK 6, could accurately replicate the ligand binding in only four of the eleven ribosomes due to the poor electron density of said ribosomal structures. In this study, we have further benchmarked the performance of the DOCK 6 docking algorithm and scoring in improving virtual screening (VS) enrichment using the dataset of 285 oxazolidinone derivatives against oxazolidinone binding sites in the S. aureus ribosome. However, there was no clear trend between the structure and activity of the oxazolidinones in VS. Overall, the docking performance indicates that the RNA pocket’s high flexibility does not allow for accurate docking prediction, highlighting the need to validate VS. protocols for ligand-RNA before future use. Later, we developed a re-scoring method incorporating absolute docking scores and molecular descriptors, and the results indicate that the descriptors greatly improve the correlation of docking scores and pMIC values. Morgan fingerprint analysis was also used, suggesting that DOCK 6 underpredicted molecules with tail modifications with acetamide, n-methylacetamide, or n-ethylacetamide and over-predicted molecule derivatives with methylamino bits. Alternatively, a ligand-based approach similar to a field template was taken, indicating that each derivative’s tail groups have strong positive and negative electrostatic potential contributing to microbial activity. These results indicate that one should perform VS. campaigns of ribosomal antibiotics with care and that more comprehensive strategies, including molecular dynamics simulations and relative free energy calculations, might be necessary in conjunction with VS. and docking.

Lessons Learnt from COVID-19: Computational Strategies for Facing Present and Future Pandemics

Since its outbreak in December 2019, the COVID-19 pandemic has caused the death of more than 6.5 million people around the world. The high transmissibility of its causative agent, the SARS-CoV-2 virus, coupled with its potentially lethal outcome, provoked a profound global economic and social crisis. The urgency of finding suitable pharmacological tools to tame the pandemic shed light on the ever-increasing importance of computer simulations in rationalizing and speeding up the design of new drugs, further stressing the need for developing quick and reliable methods to identify novel active molecules and characterize their mechanism of action. In the present work, we aim at providing the reader with a general overview of the COVID-19 pandemic, discussing the hallmarks in its management, from the initial attempts at drug repurposing to the commercialization of Paxlovid, the first orally available COVID-19 drug. Furthermore, we analyze and discuss the role of computer-aided drug discovery (CADD) techniques, especially those that fall in the structure-based drug design (SBDD) category, in facing present and future pandemics, by showcasing several successful examples of drug discovery campaigns where commonly used methods such as docking and molecular dynamics have been employed in the rational design of effective therapeutic entities against COVID-19.

Implementing a Scoring Function Based on Interaction Fingerprint for Autogrow4: Protein Kinase CK1δ as a Case Study

In the last 20 years, fragment-based drug discovery (FBDD) has become a popular and consolidated approach within the drug discovery pipeline, due to its ability to bring several drug candidates to clinical trials, some of them even being approved and introduced to the market. A class of targets that have proven to be particularly suitable for this method is represented by kinases, as demonstrated by the approval of BRAF inhibitor vemurafenib. Within this wide and diverse set of proteins, protein kinase CK1δ is a particularly interesting target for the treatment of several widespread neurodegenerative diseases, such as Alzheimer’s disease, Parkinson’s disease, and amyotrophic lateral sclerosis. Computational methodologies, such as molecular docking, are already routinely and successfully applied in FBDD campaigns alongside experimental techniques, both in the hit-discovery and in the hit-optimization stage. Concerning this, the open-source software Autogrow, developed by the Durrant lab, is a semi-automated computational protocol that exploits a combination between a genetic algorithm and a molecular docking software for de novo drug design and lead optimization. In the current work, we present and discuss a modified version of the Autogrow code that implements a custom scoring function based on the similarity between the interaction fingerprint of investigated compounds and a crystal reference. To validate its performance, we performed both a de novo and a lead-optimization run (as described in the original publication), evaluating the ability of our fingerprint-based protocol to generate compounds similar to known CK1δ inhibitors based on both the predicted binding mode and the electrostatic and shape similarity in comparison with the standard Autogrow protocol.

Discovery of S-217622, a Noncovalent Oral SARS-CoV-2 3CL Protease Inhibitor Clinical Candidate for Treating COVID-19

The coronavirus disease 2019 (COVID-19) pandemic, caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), has resulted in millions of deaths and threatens public health and safety. Despite the rapid global spread of COVID-19 vaccines, effective oral antiviral drugs are urgently needed. Here, we describe the discovery of S-217622, the first oral noncovalent, nonpeptidic SARS-CoV-2 3CL protease inhibitor clinical candidate. S-217622 was discovered via virtual screening followed by biological screening of an in-house compound library, and optimization of the hit compound using a structure-based drug design strategy. S-217622 exhibited antiviral activity in vitro against current outbreaking SARS-CoV-2 variants and showed favorable pharmacokinetic profiles in vivo for once-daily oral dosing. Furthermore, S-217622 dose-dependently inhibited intrapulmonary replication of SARS-CoV-2 in mice, indicating that this novel noncovalent inhibitor could be a potential oral agent for treating COVID-19.

The Discovery of Potent SHP2 Inhibitors with Anti-Proliferative Activity in Breast Cancer Cell Lines

Despite available treatments, breast cancer is the leading cause of cancer-related death. Knowing that the tyrosine phosphatase SHP2 is a regulator in tumorigenesis, developing inhibitors of SHP2 in breast cells is crucial. Our study investigated the effects of new compounds, purchased from NSC, on the phosphatase activity of SHP2 and the modulation of breast cancer cell lines’ proliferation and viability. A combined ligand-based and structure-based virtual screening protocol was validated, then performed, against SHP2 active site. Top ranked compounds were tested via SHP2 enzymatic assay, followed by measuring IC50 values. Subsequently, hits were tested for their anti-breast cancer viability and proliferative activity. Our experiments identified three compounds 13030, 24198, and 57774 as SHP2 inhibitors, with IC50 values in micromolar levels and considerable selectivity over the analogous enzyme SHP1. Long MD simulations of 500 ns showed a very promising binding mode in the SHP2 catalytic pocket. Furthermore, these compounds significantly reduced MCF-7 breast cancer cells’ proliferation and viability. Interestingly, two of our hits can have acridine or phenoxazine cyclic system known to intercalate in ds DNA. Therefore, our novel approach led to the discovery of SHP2 inhibitors, which could act as a starting point in the future for clinically useful anticancer agents.

Computational strategies to identify new drug candidates against neuroinflammation

The even more increasing application of computational approaches in these last decades has deeply modified the process of discovery and commercialization of new therapeutic entities. This is especially true in the field of neuroinflammation, in which both the peculiar anatomical localization and the presence of the blood-brain barrier makeit mandatory to finely tune the candidates' physicochemical properties from the early stages of the discovery pipeline. The aim of this review is therefore to provide a general overview to the readers about the topic of neuroinflammation, together with the most common computational strategies that can be exploited to discover and design small molecules controlling neuroinflammation, especially those based on the knowledge of the three-dimensional structure of the biological targets of therapeutic interest. The techniques used to describe the molecular recognition mechanisms, such as molecular docking and molecular dynamics, will therefore be eviscerated, highlighting their advantages and their limitations. Finally, we report several case studies in which computational methods have been applied in drug discovery on neuroinflammation, focusing on the last decade's research.

Combinatorial Virtual Screening Revealed a Novel Scaffold for TNKS Inhibition to Combat Colorectal Cancer

The abnormal Wnt signaling pathway leads to a high expression of β-catenin, which causes several types of cancer, particularly colorectal cancer (CRC). The inhibition of tankyrase (TNKS) activity can reduce cancer cell growth, invasion, and resistance to treatment by blocking the Wnt signaling pathway. A pharmacophore search and pharmacophore docking were performed to identify potential TNKS inhibitors in the training databases. The weighted MM/PBSA binding free energy of the docking model was calculated to rank the databases. The reranked results indicated that 26.98% of TNKS inhibitors that were present in the top 5% of compounds in the database and near an ideal value ranked 28.57%. The National Cancer Institute database was selected for formal virtual screening, and 11 potential TNKS inhibitors were identified. An enzyme-based experiment was performed to demonstrate that of the 11 potential TNKS inhibitors, NSC295092 and NSC319963 had the most potential. Finally, Wnt pathway analysis was performed through a cell-based assay, which indicated that NSC319963 is the most likely TNKS inhibitor (pIC50 = 5.59). The antiproliferation assay demonstrated that NSC319963 can decrease colorectal cancer cell growth; therefore, the proposed method successfully identified a novel TNKS inhibitor that can alleviate CRC.

Identify potent SARS-CoV-2 main protease inhibitors via accelerated free energy perturbation-based virtual screening of existing drugs

Drug repurposing effort for treatment of a new disease, such as COVID-19, usually starts from a virtual screening of existing drugs, followed by experimental validation, but the actual hit rate is generally rather low with traditional computational methods. It has been demonstrated that a new virtual screening approach with accelerated free energy perturbation-based absolute binding free energy (FEP-ABFE) predictions can reach an unprecedentedly high hit rate, leading to successful identification of 15 potent inhibitors of SARS-CoV-2 main protease (M^pro) from 25 computationally selected drugs under a threshold of K_i = 4 μM. The outcomes of this study are valuable for not only drug repurposing to treat COVID-19 but also demonstrating the promising potential of the FEP-ABFE prediction-based virtual screening approach.

The COVID-19 pandemic caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has become a global crisis. There is no therapeutic treatment specific for COVID-19. It is highly desirable to identify potential antiviral agents against SARS-CoV-2 from existing drugs available for other diseases and thus repurpose them for treatment of COVID-19. In general, a drug repurposing effort for treatment of a new disease, such as COVID-19, usually starts from a virtual screening of existing drugs, followed by experimental validation, but the actual hit rate is generally rather low with traditional computational methods. Here we report a virtual screening approach with accelerated free energy perturbation-based absolute binding free energy (FEP-ABFE) predictions and its use in identifying drugs targeting SARS-CoV-2 main protease (M^pro). The accurate FEP-ABFE predictions were based on the use of a restraint energy distribution (RED) function, making the practical FEP-ABFE−based virtual screening of the existing drug library possible. As a result, out of 25 drugs predicted, 15 were confirmed as potent inhibitors of SARS-CoV-2 M^pro. The most potent one is dipyridamole (inhibitory constant K_i = 0.04 µM) which has shown promising therapeutic effects in subsequently conducted clinical studies for treatment of patients with COVID-19. Additionally, hydroxychloroquine (K_i = 0.36 µM) and chloroquine (K_i = 0.56 µM) were also found to potently inhibit SARS-CoV-2 M^pro. We anticipate that the FEP-ABFE prediction-based virtual screening approach will be useful in many other drug repurposing or discovery efforts.

Computational Modeling on Aquaporin-3 as Skin Cancer Target: A Virtual Screening Study

Aquaporin-3 (AQP3) is one of the aquaglyceroporins, which is expressed in the basolateral layer of the skin membrane. Studies have reported that human skin squamous cell carcinoma overexpresses AQP3 and inhibition of its function may alleviate skin tumorigenesis. In the present study, we have applied a virtual screening method that encompasses filters for physicochemical properties and molecular docking to select potential hit compounds that bind to the Aquaporin-3 protein. Based on molecular docking results, the top 20 hit compounds were analyzed for stability in the binding pocket using unconstrained molecular dynamics simulations and further evaluated for binding free energy. Furthermore, examined the ligand-unbinding pathway of the inhibitor from its bound form to explore possible routes for inhibitor approach to the ligand-binding site. With a good docking score, stability in the binding pocket, and free energy of binding, these hit compounds can be developed as Aquaporin-3 inhibitors in the near future.

PHARMIP: An insilico method to predict genetics that underpin adverse drug reactions

Pharmacovigilance is the pharmacological science that focuses on the safe and appropriate use of drugs.Variability in response to drug therapy in both terms of safety and efficacy is highly related to patient's personal genomics. Hence, pharmacovigilance considers pharmacogenomics methodologies in the evaluation of medicinal products. The aim of this work is to introduce the pharmacovigilance/ pharmacogenomics insilico pipeline (PHARMIP) that uses the drug (or drug candidate) digital structure and the advances in bioinformatics tools and databases to figure-out the genetic factors underlying the drug reported adverse reactions (ADRs).PHARMIP uses user-friendly freely available bioinformatics resources to help pharmacovigilance and pharmacogenomics scientists with minimal bioinformatics experience to retrieve helpful information for their daily basis activities. Also, PHARMIP could help the advances in precision medicine in a drug-centric approach as it can be used to reveal genetic risk factors for certain drug ADRs. Domperidone was used as an example to the application of PHARMIP as the pipeline was initially developed during the insilico exploration of domperidone cardiotoxic ADRs.

Method is composed of 3 main steps:

•

Preparing the drug off-label targets (OLT) list.

•

Retrieving the related diseases/ adverse reactions (DA) list.

•

Analysis of DA list to get answers.

Discovery of GlyT2 Inhibitors Using Structure-Based Pharmacophore Screening and Selectivity Studies by FEP+ Calculations

In recent years, mammalian Glycine transporter 2 (GlyT2) has emerged as a promising target for the development of compounds against chronic pain states. In our current work, we discovered a new set of promising hits that inhibit the glycine transporter at nano- and micromolar activity and have excellent selectivity over GlyT1 (as shown by in vitro studies) using a newly designed virtual screening (VS) protocol that combines a structure-based pharmacophore and docking screens with a success rate of 75%. Furthermore, the free energy perturbation calculations and molecular dynamics (MD) studies revealed the GlyT2 amino acid residues critical for the binding and selectivity of both Glycine and our Hit1 compound. The FEP+ results well-matched with the available literature mutational data proving the quality of the generated GlyT2 structure. On the basis of these results, we propose that our hit compounds may lead to new chronic pain agents to address unmet and challenging clinical needs.

NMR-Fragment Based Virtual Screening: A Brief Overview

Fragment-based drug discovery (FBDD) using NMR has become a central approach over the last twenty years for development of small molecule inhibitors against biological macromolecules, to control a variety of cellular processes. Yet, several considerations should be taken into account for obtaining a therapeutically relevant agent. In this review, we aim to list the considerations that make NMR fragment screening a successful process for yielding potent inhibitors. Factors that may govern the competence of NMR in fragment based drug discovery are discussed, as well as later steps that involve optimization of hits obtained by NMR-FBDD.

Best Practices of Computer-Aided Drug Discovery: Lessons Learned from the Development of a Preclinical Candidate for Prostate Cancer with a New Mechanism of Action

Small-molecule drug design is a complex and iterative decision-making process relying on pre-existing knowledge and driven by experimental data. Low-molecular-weight chemicals represent an attractive therapeutic option, as they are readily accessible to organic synthesis and can easily be characterized.1 Their potency as well as pharmacokinetic and pharmacodynamic properties can be systematically and rationally investigated and ultimately optimized via expert science behind medicinal chemistry and methods of computer-aided drug design (CADD). In recent years, significant advances in molecular modeling techniques have afforded a variety of tools to effectively identify potential binding pockets on prospective targets, to map key interactions between ligands and their binding sites, to construct and assess energetics of the resulting complexes, to predict ADMET properties of candidate compounds, and to systematically analyze experimental and computational data to derive meaningful structure-activity relationships leading to the creation of a drug candidate. This Perspective describes a real case of a drug discovery campaign accomplished in a relatively short time with limited resources. The study integrated an arsenal of available molecular modeling techniques with an array of experimental tools to successfully develop a novel class of potent and selective androgen receptor inhibitors with a novel mode of action. It resulted in the largest academic licensing deal in Canadian history, totaling $142M. This project exemplifies the importance of team science, an integrative approach to drug discovery, and the use of best practices in CADD. We posit that the lessons learned and best practices for executing an effective CADD project can be applied, with similar success, to many drug discovery projects in both academia and industry.

Azathioprine, Mercaptopurine, and 5-Aminosalicylic Acid Affect the Growth of IBD-Associated Campylobacter Species and Other Enteric Microbes

Campylobacter concisus is a bacterium that is associated with inflammatory bowel disease (IBD). Immunosuppressive drugs including azathioprine (AZA) and mercaptopurine (MP), and anti-inflammatory drug such as 5-aminosalicylic acid (5-ASA) are commonly used to treat patients with IBD. This study aimed to examine the effects of AZA, MP, and 5-ASA on the growth of IBD-associated bacterial species and to identify bacterial enzymes involved in immunosuppressive drug metabolism. A total of 15 bacterial strains of five species including 11 C. concisus strains, Bacteroides fragilis, Bacteroides vulgatus, Enterococcus faecalis, and Escherichia coli were examined. The impact of AZA, MP, and 5-ASA on the growth of these bacterial species was examined quantitatively using a plate counting method. The presence of enzymes involved in AZA and MP metabolism in these bacterial species was identified using bioinformatics tools. AZA and MP significantly inhibited the growth of all 11 C. concisus strains. C. concisus strains were more sensitive to AZA than MP. 5-ASA showed inhibitory effects to some C. concisus strains, while it promoted the growth of other C. concisus strains. AZA and MP also significantly inhibited the growth of B. fragilis and B. vulgatus. The growth of E. coli was significantly inhibited by 200 μg/ml of AZA as well as 100 and 200 μg/ml of 5-ASA. Bacterial enzymes related to AZA and MP metabolism were found, which varied in different bacterial species. In conclusion, AZA and MP have inhibitory effects to IBD-associated C. concisus and other enteric microbes, suggesting an additional therapeutic mechanism of these drugs in the treatment of IBD. The strain dependent differential impact of 5-ASA on the growth of C. concisus may also have clinical implication given that in some cases 5-ASA medications were found to cause exacerbations of colitis.

Common Hits Approach: Combining Pharmacophore Modeling and Molecular Dynamics Simulations

We present a new approach that incorporates flexibility based on extensive MD simulations of protein-ligand complexes into structure-based pharmacophore modeling and virtual screening. The approach uses the multiple coordinate sets saved during the MD simulations and generates for each frame a pharmacophore model. Pharmacophore models with the same pharmacophore features are pooled. In this way the high number of pharmacophore models that results from the MD simulation is reduced to only a few hundred representative pharmacophore models. Virtual screening runs are performed with every representative pharmacophore model; the screening results are combined and rescored to generate a single hit-list. The score for a particular molecule is calculated based on the number of representative pharmacophore models which classified it as active. Hence, the method is called common hits approach (CHA). The steps between the MD simulation and the final hit-list are performed automatically and without user interaction. We test the performance of CHA for virtual screening using screening databases with active and inactive compounds for 40 protein-ligand systems. The results of the CHA are compared to the (i) median screening performance of all representative pharmacophore models of protein-ligand systems, as well as to the virtual screening performance of (ii) a random classifier, (iii) the pharmacophore model derived from the experimental structure in the PDB, and (iv) the representative pharmacophore model appearing most frequently during the MD simulation. For the 34 (out of 40) protein-ligand complexes, for which at least one of the approaches was able to perform better than a random classifier, the highest enrichment was achieved using CHA in 68% of the cases, compared to 12% for the PDB pharmacophore model and 20% for the representative pharmacophore model appearing most frequently. The availabilithy of diverse sets of different pharmacophore models is utilized to analyze some additional questions of interest in 3D pharmacophore-based virtual screening.

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.

In silico screening and study of novel ERK2 inhibitors using 3D QSAR, docking and molecular dynamics

ERK2 is a dual specificity protein kinase, part of the Ras/Raf/MEK/ERK signal transduction cascade. It forms an interesting target for inhibition based on its relationship with cell proliferation and oncogenesis. A 3D QSAR pharmacophore model (Hypo1) with high correlation (r=0.938) was developed for ERK2 ATP site on the basis of experimentally known inhibitors. The model included three hydrogen bonds, and one hydrophobic site. Assessment of Hypo1 through Fisher randomization, cost analysis, leave one out method and decoy test suggested that the model can reliably detect ERK2 inhibitors. Hypo1 has been used for virtual screening of potential inhibitors from ZINC, Drug Bank, NCI, Maybridge and Chembank databases. Using Hypo1 as a query, databases have been interrogated for compounds who meet the pharmacophore features. The resulting hit compounds were subject to docking and analysis. Docking and molecular dynamics analysis showed that in order to achieve a higher potency compounds have to interact with catalytic site, glycine rich loop, Hinge region, Gatekeeper region and ATP site entrance residues. We also identified catalytic site and Glycine rich loop as important regions to bind by molecules for better potency and selectivity.

An integrated approach to knowledge-driven structure-based virtual screening

In many practical applications of structure-based virtual screening (VS) ligands are already known. This circumstance requires that the obtained hits need to satisfy initial made expectations i.e., they have to fulfill a predefined binding pattern and/or lie within a predefined physico-chemical property range. Based on the RApid Index-based Screening Engine (RAISE) approach, we introduce CRAISE-a user-controllable structure-based VS method. It efficiently realizes pharmacophore-guided protein-ligand docking to assess the library content but thereby concentrates only on molecules that have a chance to fulfill the given binding pattern. In order to focus only on hits satisfying given molecular properties, library profiles can be utilized to simultaneously filter compounds. CRAISE was evaluated on a range of strict to rather relaxed hypotheses with respect to its capability to guide binding-mode predictions and VS runs. The results reveal insights into a guided VS process. If a pharmacophore model is chosen appropriately, a binding mode below 2 Å is successfully reproduced for 85% of well-prepared structures, enrichment is increased up to median AUC of 73%, and the selectivity of the screening process is significantly enhanced leading up to seven times accelerated runtimes. In general, CRAISE supports a versatile structure-based VS approach allowing to assess hypotheses about putative ligands on a large scale.

Software and resources for computational medicinal chemistry

Computer-aided drug design plays a vital role in drug discovery and development and has become an indispensable tool in the pharmaceutical industry. Computational medicinal chemists can take advantage of all kinds of software and resources in the computer-aided drug design field for the purposes of discovering and optimizing biologically active compounds. This article reviews software and other resources related to computer-aided drug design approaches, putting particular emphasis on structure-based drug design, ligand-based drug design, chemical databases and chemoinformatics tools.

Reviewing ligand-based rational drug design: the search for an ATP synthase inhibitor

Following major advances in the field of medicinal chemistry, novel drugs can now be designed systematically, instead of relying on old trial and error approaches. Current drug design strategies can be classified as being either ligand- or structure-based depending on the design process. In this paper, by describing the search for an ATP synthase inhibitor, we review two frequently used approaches in ligand-based drug design: The pharmacophore model and the quantitative structure-activity relationship (QSAR) method. Moreover, since ATP synthase ligands are potentially useful drugs in cancer therapy, pharmacophore models were constructed to pave the way for novel inhibitor designs.

Recent trends and observations in the design of high-quality screening collections

The design of a high-quality screening collection is of utmost importance for the early drug-discovery process and provides, in combination with high-quality assay systems, the foundation of future discoveries. Herein, we review recent trends and observations to successfully expand the access to bioactive chemical space, including the feedback from hit assessment interviews of high-throughput screening campaigns; recent successes with chemogenomics target family approaches, the identification of new relevant target/domain families, diversity-oriented synthesis and new emerging compound classes, and non-classical approaches, such as fragment-based screening and DNA-encoded chemical libraries. The role of in silico library design approaches are emphasized.

Inhibitors of Helicobacter pylori protease HtrA found by 'virtual ligand' screening combat bacterial invasion of epithelia

Background

The human pathogen Helicobacter pylori (H. pylori) is a main cause for gastric inflammation and cancer. Increasing bacterial resistance against antibiotics demands for innovative strategies for therapeutic intervention.

Methodology/Principal Findings

We present a method for structure-based virtual screening that is based on the comprehensive prediction of ligand binding sites on a protein model and automated construction of a ligand-receptor interaction map. Pharmacophoric features of the map are clustered and transformed in a correlation vector (‘virtual ligand’) for rapid virtual screening of compound databases. This computer-based technique was validated for 18 different targets of pharmaceutical interest in a retrospective screening experiment. Prospective screening for inhibitory agents was performed for the protease HtrA from the human pathogen H. pylori using a homology model of the target protein. Among 22 tested compounds six block E-cadherin cleavage by HtrA in vitro and result in reduced scattering and wound healing of gastric epithelial cells, thereby preventing bacterial infiltration of the epithelium.

Conclusions/Significance

This study demonstrates that receptor-based virtual screening with a permissive (‘fuzzy’) pharmacophore model can help identify small bioactive agents for combating bacterial infection.

Virtual screening of selective multitarget kinase inhibitors by combinatorial support vector machines

Multitarget agents have been increasingly explored for enhancing efficacy and reducing countertarget activities and toxicities. Efficient virtual screening (VS) tools for searching selective multitarget agents are desired. Combinatorial support vector machines (C-SVM) were tested as VS tools for searching dual-inhibitors of 11 combinations of 9 anticancer kinase targets (EGFR, VEGFR, PDGFR, Src, FGFR, Lck, CDK1, CDK2, GSK3). C-SVM trained on 233-1,316 non-dual-inhibitors correctly identified 26.8%-57.3% (majority >36%) of the 56-230 intra-kinase-group dual-inhibitors (equivalent to the 50-70% yields of two independent individual target VS tools), and 12.2% of the 41 inter-kinase-group dual-inhibitors. C-SVM were fairly selective in misidentifying as dual-inhibitors 3.7%-48.1% (majority <20%) of the 233-1,316 non-dual-inhibitors of the same kinase pairs and 0.98%-4.77% of the 3,971-5,180 inhibitors of other kinases. C-SVM produced low false-hit rates in misidentifying as dual-inhibitors 1,746-4,817 (0.013%-0.036%) of the 13.56 M PubChem compounds, 12-175 (0.007%-0.104%) of the 168 K MDDR compounds, and 0-84 (0.0%-2.9%) of the 19,495-38,483 MDDR compounds similar to the known dual-inhibitors. C-SVM was compared to other VS methods Surflex-Dock, DOCK Blaster, kNN and PNN against the same sets of kinase inhibitors and the full set or subset of the 1.02 M Zinc clean-leads data set. C-SVM produced comparable dual-inhibitor yields, slightly better false-hit rates for kinase inhibitors, and significantly lower false-hit rates for the Zinc clean-leads data set. Combinatorial SVM showed promising potential for searching selective multitarget agents against intra-kinase-group kinases without explicit knowledge of multitarget agents.