Building a virtual ligand screening pipeline using free software: a survey

Innovative computational approaches in drug discovery and design

In the current scenario of pandemics, drug discovery and design have undergone a significant transformation due to the integration of advanced computational methodologies. These methodologies utilize sophisticated algorithms, machine learning, artificial intelligence, and high-performance computing to expedite the drug development process, enhances accuracy, and reduces costs. Machine learning and AI have revolutionized predictive modeling, virtual screening, and de novo drug design, allowing for the identification and optimization of novel compounds with desirable properties. Molecular dynamics simulations provide a detailed insight into protein-ligand interactions and conformational changes, facilitating an understanding of drug efficacy at the atomic level. Quantum mechanics/molecular mechanics methods offer precise predictions of binding energies and reaction mechanisms, while structure-based drug design employs docking studies and fragment-based design to improve drug-receptor binding affinities. Network pharmacology and systems biology approaches analyze polypharmacology and biological networks to identify novel drug targets and understand complex interactions. Cheminformatics explores vast chemical spaces and employs data mining to find patterns in large datasets. Computational toxicology predicts adverse effects early in development, reducing reliance on animal testing. Bioinformatics integrates genomic, proteomic, and metabolomics data to discover biomarkers and understand genetic variations affecting drug response. Lastly, cloud computing and big data technologies facilitate high-throughput screening and comprehensive data analysis. Collectively, these computational innovations are driving a paradigm shift in drug discovery and design, making it more efficient, accurate, and cost-effective.

In silico prospection of receptors associated with the biological activity of U1-SCTRX-lg1a: an antimicrobial peptide isolated from the venom of Loxosceles gaucho

The emergence of antibiotic-resistant pathogens generates impairment to human health. U1-SCTRX-lg1a is a peptide isolated from a phospholipase D extracted from the spider venom of Loxosceles gaucho with antimicrobial activity against Gram-negative bacteria (between 1.15 and 4.6 μM). The aim of this study was to suggest potential receptors associated with the antimicrobial activity of U1-SCTRX-lg1a using in silico bioinformatics tools. The search for potential targets of U1-SCRTX-lg1a was performed using the PharmMapper server. Molecular docking between U1-SCRTX-lg1a and the receptor was performed using PatchDock software. The prediction of ligand sites for each receptor was conducted using the PDBSum server. Chimera 1.6 software was used to perform molecular dynamics simulations only for the best dock score receptor. In addition, U1-SCRTX-lg1a and native ligand interactions were compared using AutoDock Vina software. Finally, predicted interactions were compared with the ligand site previously described in the literature. The bioprospecting of U1-SCRTX-lg1a resulted in the identification of three hundred (300) diverse targets (Table S1), forty-nine (49) of which were intracellular proteins originating from Gram-negative microorganisms (Table S2). Docking results indicate Scores (10,702 to 6066), Areas (1498.70 to 728.40) and ACEs (417.90 to - 152.8) values. Among these, NAD + NH3-dependent synthetase (PDB ID: 1wxi) showed a dock score of 9742, area of 1223.6 and ACE of 38.38 in addition to presenting a Normalized Fit score of 8812 on PharmMapper server. Analysis of the interaction of ligands and receptors suggests that the peptide derived from brown spider venom can interact with residues SER48 and THR160. Furthermore, the C terminus (- 7.0 score) has greater affinity for the receptor than the N terminus (- 7.7 score). The molecular dynamics assay shown that free energy value for the protein complex of - 214,890.21 kJ/mol, whereas with rigid docking, this value was - 29.952.8 sugerindo that after the molecular dynamics simulation, the complex exhibits a more favorable energy value compared to the previous state. The in silico bioprospecting of receptors suggests that U1-SCRTX-lg1a may interfere with NAD + production in Escherichia coli, a Gram-negative bacterium, altering the homeostasis of the microorganism and impairing growth.

Supplementary Information

The online version contains supplementary material available at 10.1007/s40203-024-00190-8.

Introducing the Automated Ligand Searcher

The Automated Ligand Searcher (ALISE) is designed as an automated computational drug discovery tool. To approximate the binding free energy of ligands to a receptor, ALISE includes a three-stage workflow, with each stage involving an increasingly sophisticated computational method: molecular docking, molecular dynamics, and free energy perturbation, respectively. To narrow the number of potential ligands, poorly performing ligands are gradually segregated out. The performance and usability of ALISE are benchmarked for a case study containing known active ligands and decoys for the HIV protease. The example illustrates that ALISE filters the decoys successfully and demonstrates that the automation, comprehensiveness, and user-friendliness of the software make it a valuable tool for improved and faster drug development workflows.

A novel in silico scaffold-hopping method for drug repositioning in rare and intractable diseases

In the field of rare and intractable diseases, new drug development is difficult and drug repositioning (DR) is a key method to improve this situation. In this study, we present a new method for finding DR candidates utilizing virtual screening, which integrates amino acid interaction mapping into scaffold-hopping (AI-AAM). At first, we used a spleen associated tyrosine kinase inhibitor as a reference to evaluate the technique, and succeeded in scaffold-hopping maintaining the pharmacological activity. Then we applied this method to five drugs and obtained 144 compounds with diverse structures. Among these, 31 compounds were known to target the same proteins as their reference compounds and 113 compounds were known to target different proteins. We found that AI-AAM dominantly selected functionally similar compounds; thus, these selected compounds may represent improved alternatives to their reference compounds. Moreover, the latter compounds were presumed to bind to the targets of their references as well. This new "compound-target" information provided DR candidates that could be utilized for future drug development.

Identification of novel inhibitors of S-adenosyl-L-homocysteine hydrolase via structure-based virtual screening and molecular dynamics simulations

S-adenosyl-L-homocysteine hydrolase (SAHase) is an important regulator in the methylation reactions in many organisms and thus is crucial for numerous cellular functions. In recent years, SAHase has become one of the popular targets for drug design, and SAHase inhibitors have exhibited potent antiviral activity. In this study, we established the complex-based pharmacophore models based on the known crystal complex of SAHase (PDB ID: 1A7A) to screen the drug-likeness compounds of ChEMBL database. Then, three molecular docking programs were used to validate the reliability of compounds, involving Libdock, CDOCKER, and AutoDock Vina programs. The four promising hit compounds (CHEMBL420751, CHEMBL346387, CHEMBL1569958, and CHEMBL4206648) were performed molecular dynamics simulations and MM-PBSA calculations to evaluate their stability and binding-free energy in the binding site of SAHase. After screening and analyzing, the hit compounds CHEMBL420751 and CHEMBL346387 were suggested to further research to obtain novel potential SAHase inhibitors. A series of computer-aided drug design methods, including pharmacophore, molecular docking, molecular dynamics simulation and MM-PBSA calculations, were employed in this study to identity novel inhibitors of S-adenosyl-L-homocysteine hydrolase (SAHase). Some compounds from virtual screening could form various interactions with key residues of SAHase. Among them, compounds CHEMBL346387 and CHEMBL420751 exhibited potent binding affinity from molecular docking and MM-PBSA, and maintained good stability at the binding site during molecular dynamics simulations as well. All these results indicated that the selected compounds might have the potential to be novel SAHase inhibitors.

Continuous Manufacturing and Molecular Modeling of Pharmaceutical Amorphous Solid Dispersions

Amorphous solid dispersions enhance solubility and oral bioavailability of poorly water-soluble drugs. The escalating number of drugs with poor aqueous solubility, poor dissolution, and poor oral bioavailability is an unresolved problem that requires adequate interventions. This review article highlights recent solubility and bioavailability enhancement advances using amorphous solid dispersions (ASDs). The review also highlights the mechanism of enhanced dissolution and the challenges faced by ASD-based products, such as stability and scale-up. The role of process analytical technology (PAT) supporting continuous manufacturing is highlighted. Accurately predicting interactions between the drug and polymeric carrier requires long experimental screening methods, and this is a space where computational tools hold significant potential. Recent advancements in data science, computational tools, and easy access to high-end computation power are set to accelerate ASD-based research. Hence, particular emphasis has been given to molecular modeling techniques that can address some of the unsolved questions related to ASDs. With the advancement in PAT tools and artificial intelligence, there is an increasing interest in the continuous manufacturing of pharmaceuticals. ASDs are a suitable option for continuous manufacturing, as production of a drug product from an ASD by direct compression is a reality, where the addition of multiple excipients is easy to avoid. Significant attention is necessary for ongoing clinical studies based on ASDs, which is paving the way for the approval of many new ASDs and their introduction into the market.

The Current Perspectives in Clinical Research: Computer-Assisted Drug Designing, Ethics, and Good Clinical Practice

In the era of emerging microbial and non-communicable diseases and re-emerging microbial infections, the medical fraternity and the public are plagued by under-preparedness. It is evident by the severity of the Coronavirus disease (COVID-19) pandemic that novel microbial diseases are a challenge and are challenging to control. This is mainly attributed to the lack of complete knowledge of the novel microbe’s biology and pathogenesis and the unavailability of therapeutic drugs and vaccines to treat and control the disease. Clinical research is the only answer utilizing which can handle most of these circumstances. In this review, we highlight the importance of computer-assisted drug designing (CADD) and the aspects of molecular docking, molecular superimposition, 3D-pharmacophore technology, ethics, and good clinical practice (GCP) for the development of therapeutic drugs, devices, and vaccines.

An Overview of Antiviral Peptides and Rational Biodesign Considerations

Viral diseases have contributed significantly to worldwide morbidity and mortality throughout history. Despite the existence of therapeutic treatments for many viral infections, antiviral resistance and the threat posed by novel viruses highlight the need for an increased number of effective therapeutics. In addition to small molecule drugs and biologics, antimicrobial peptides (AMPs) represent an emerging class of potential antiviral therapeutics. While AMPs have traditionally been regarded in the context of their antibacterial activities, many AMPs are now known to be antiviral. These antiviral peptides (AVPs) have been shown to target and perturb viral membrane envelopes and inhibit various stages of the viral life cycle, from preattachment inhibition through viral release from infected host cells. Rational design of AMPs has also proven effective in identifying highly active and specific peptides and can aid in the discovery of lead peptides with high therapeutic selectivity. In this review, we highlight AVPs with strong antiviral activity largely curated from a publicly available AMP database. We then compile the sequences present in our AVP database to generate structural predictions of generic AVP motifs. Finally, we cover the rational design approaches available for AVPs taking into account approaches currently used for the rational design of AMPs.

Binding site identification of G protein-coupled receptors through a 3D Zernike polynomials-based method: application to C. elegans olfactory receptors

Studying the binding processes of G protein-coupled receptors (GPCRs) proteins is of particular interest both to better understand the molecular mechanisms that regulate the signaling between the extracellular and intracellular environment and for drug design purposes. In this study, we propose a new computational approach for the identification of the binding site for a specific ligand on a GPCR. The method is based on the Zernike polynomials and performs the ligand-GPCR association through a shape complementarity analysis of the local molecular surfaces. The method is parameter-free and it can distinguish, working on hundreds of experimentally GPCR-ligand complexes, binding pockets from randomly sampled regions on the receptor surface, obtaining an Area Under ROC curve of 0.77. Given its importance both as a model organism and in terms of applications, we thus investigated the olfactory receptors of the C. elegans, building a list of associations between 21 GPCRs belonging to its olfactory neurons and a set of possible ligands. Thus, we can not only carry out rapid and efficient screenings of drugs proposed for GPCRs, key targets in many pathologies, but also we laid the groundwork for computational mutagenesis processes, aimed at increasing or decreasing the binding affinity between ligands and receptors.

Screening S protein - ACE2 blockers from natural products: Strategies and advances in the discovery of potential inhibitors of COVID-19

The Coronavirus disease, 2019 (COVID-19) is caused by severe acute respiratory syndrome Coronavirus 2 (SARS-CoV-2), which poses a major threat to human life and health. Given its continued development, limiting the spread of COVID-19 in the population remains a challenging task. Currently, multiple therapies are being tried around the world to deal with SARS-CoV-2 infection, and a variety of studies have shown that natural products have a significant effect on COVID-19 patients. The combination of SARS-CoV-2 S protein with Angiotensin converting enzyme II(ACE2) of host cell to promote membrane fusion is an initial critical step for SARS-CoV-2 infection. Therefore, screening natural products that inhibit the binding of SARS-CoV-2 S protein and ACE2 also provides a feasible strategy for the treatment of COVID-19. Establishment of high throughput screening model is an important basis and key technology for screening S protein-ACE2 blockers. Based on this, the molecular structures of SARS-CoV-2 and ACE2 and their processes in the life cycle of SARS-CoV-2 and host cell infection were firstly reviewed in this paper, with emphasis on the methods and techniques of screening S protein-ACE2 blockers, including Virtual Screening (VS), Surface Plasmon Resonance (SPR), Biochromatography, Biotin-avidin with Enzyme-linked Immunosorbent assay and Gene Chip Technology. Furthermore, the technical principle, advantages and disadvantages and application scope were further elaborated. Combined with the application of the above screening technologies in S protein-ACE2 blockers, a variety of natural products, such as flavonoids, terpenoids, phenols, alkaloids, were summarized, which could be used as S protein-ACE2 blockers, in order to provide ideas for the efficient discovery of S protein-ACE2 blockers from natural sources and contribute to the development of broad-spectrum anti coronavirus drugs.

Using Stochastic Computing for Virtual Screening Acceleration

Stochastic computing is an emerging scientific field pushed by the need for developing high-performance artificial intelligence systems in hardware to quickly solve complex data processing problems. This is the case of virtual screening, a computational task aimed at searching across huge molecular databases for new drug leads. In this work, we show a classification framework in which molecules are described by an energy-based vector. This vector is then processed by an ultra-fast artificial neural network implemented through FPGA by using stochastic computing techniques. Compared to other previously published virtual screening methods, this proposal provides similar or higher accuracy, while it improves processing speed by about two or three orders of magnitude.

Pharmacophore Model for SARS-CoV-2 3CLpro Small-Molecule Inhibitors and in Vitro Experimental Validation of Computationally Screened Inhibitors

Among the biomedical efforts in response to the current coronavirus (COVID-19) pandemic, pharmacological strategies to reduce viral load in patients with severe forms of the disease are being studied intensively. One of the main drug target proteins proposed so far is the SARS-CoV-2 viral protease 3CLpro (also called Mpro), an essential component for viral replication. Ongoing ligand- and receptor-based computational screening efforts would be facilitated by an improved understanding of the electrostatic, hydrophobic, and steric features that characterize small-molecule inhibitors binding stably to 3CLpro and by an extended collection of known binders. Here, we present combined virtual screening, molecular dynamics (MD) simulation, machine learning, and in vitro experimental validation analyses, which have led to the identification of small-molecule inhibitors of 3CLpro with micromolar activity and to a pharmacophore model that describes functional chemical groups associated with the molecular recognition of ligands by the 3CLpro binding pocket. Experimentally validated inhibitors using a ligand activity assay include natural compounds with the available prior knowledge on safety and bioavailability properties, such as the natural compound rottlerin (IC50 = 37 μM) and synthetic compounds previously not characterized (e.g., compound CID 46897844, IC50 = 31 μM). In combination with the developed pharmacophore model, these and other confirmed 3CLpro inhibitors may provide a basis for further similarity-based screening in independent compound databases and structural design optimization efforts to identify 3CLpro ligands with improved potency and selectivity. Overall, this study suggests that the integration of virtual screening, MD simulations, and machine learning can facilitate 3CLpro-targeted small-molecule screening investigations. Different receptor-, ligand-, and machine learning-based screening strategies provided complementary information, helping to increase the number and diversity of the identified active compounds. Finally, the resulting pharmacophore model and experimentally validated small-molecule inhibitors for 3CLpro provide resources to support follow-up computational screening efforts for this drug target.

Comprehensive virtual screening of 4.8 k flavonoids reveals novel insights into allosteric inhibition of SARS-CoV-2 MPRO

SARS-CoV-2 main protease is a common target for inhibition assays due to its high conservation among coronaviruses. Since flavonoids show antiviral activity, several in silico works have proposed them as potential SARS-CoV-2 main protease inhibitors. Nonetheless, there is reason to doubt certain results given the lack of consideration for flavonoid promiscuity or main protease plasticity, usage of short library sizes, absence of control molecules and/or the limitation of the methodology to a single target site. Here, we report a virtual screening study where dorsilurin E, euchrenone a11, sanggenol O and CHEMBL2171598 are proposed to inhibit main protease through different pathways. Remarkably, novel structural mechanisms were observed after sanggenol O and CHEMBL2171598 bound to experimentally proven allosteric sites. The former drastically affected the active site, while the latter triggered a hinge movement which has been previously reported for an inactive SARS-CoV main protease mutant. The use of a curated database of 4.8 k flavonoids, combining two well-known docking software (AutoDock Vina and AutoDock4.2), molecular dynamics and MMPBSA, guaranteed an adequate analysis and robust interpretation. These criteria can be considered for future screening campaigns against SARS-CoV-2 main protease.

PyRMD: A New Fully Automated AI-Powered Ligand-Based Virtual Screening Tool

Artificial intelligence (AI) algorithms are dramatically redefining the current drug discovery landscape by boosting the efficiency of its various steps. Still, their implementation often requires a certain level of expertise in AI paradigms and coding. This often prevents the use of these powerful methodologies by non-expert users involved in the design of new biologically active compounds. Here, the random matrix discriminant (RMD) algorithm, a high-performance AI method specifically tailored for the identification of new ligands, was implemented in a new fully automated tool, PyRMD. This ligand-based virtual screening tool can be trained using target bioactivity data directly downloaded from the ChEMBL repository without manual intervention. The software automatically splits the available training compounds into active and inactive sets and learns the distinctive chemical features responsible for the compounds' activity/inactivity. PyRMD was designed to easily screen millions of compounds in hours through an automated workflow and intuitive input files, allowing fine tuning of each parameter of the calculation. Additionally, PyRMD features a wealth of benchmark metrics, to accurately probe the model performance, which were used here to gauge its predictive potential and limitations. PyRMD is freely available on GitHub (https://github.com/cosconatilab/PyRMD) as an open-source tool.

Current trends in computer aided drug design and a highlight of drugs discovered via computational techniques: A review

Computer-aided drug design (CADD) is one of the pivotal approaches to contemporary pre-clinical drug discovery, and various computational techniques and software programs are typically used in combination, in a bid to achieve the desired outcome. Several approved drugs have been developed with the aid of CADD. On SciFinder®, we evaluated more than 600 publications through systematic searching and refining, using the terms, virtual screening; software methods; computational studies and publication year, in order to obtain data concerning particular aspects of CADD. The primary focus of this review was on the databases screened, virtual screening and/or molecular docking software program used. Furthermore, we evaluated the studies that subsequently performed molecular dynamics (MD) simulations and we reviewed the software programs applied, the application of density functional theory (DFT) calculations and experimental assays. To represent the latest trends, the most recent data obtained was between 2015 and 2020, consequently the most frequently employed techniques and software programs were recorded. Among these, the ZINC database was the most widely preferred with an average use of 31.2%. Structure-based virtual screening (SBVS) was the most prominently used type of virtual screening and it accounted for an average of 57.6%, with AutoDock being the preferred virtual screening/molecular docking program with 41.8% usage. Following the screening process, 38.5% of the studies performed MD simulations to complement the virtual screening and GROMACS with 39.3% usage, was the popular MD software program. Among the computational techniques, DFT was the least applied whereby it only accounts for 0.02% average use. An average of 36.5% of the studies included reports on experimental evaluations following virtual screening. Ultimately, since the inception and application of CADD in pre-clinical drug discovery, more than 70 approved drugs have been discovered, and this number is steadily increasing over time.

Docking Paradigm in Drug Design

Docking is in demand for the rational computer aided structure based drug design. A review of docking methods and programs is presented. Different types of docking programs are described. They include docking of non-covalent small ligands, protein-protein docking, supercomputer docking, quantum docking, the new generation of docking programs and the application of docking for covalent inhibitors discovery. Taking into account the threat of COVID-19, we present here a short review of docking applications to the discovery of inhibitors of SARS-CoV and SARS-CoV-2 target proteins, including our own result of the search for inhibitors of SARS-CoV-2 main protease using docking and quantum chemical post-processing. The conclusion is made that docking is extremely important in the fight against COVID-19 during the process of development of antivirus drugs having a direct action on SARS-CoV-2 target proteins.

Network bioinformatics analysis provides insight into drug repurposing for COVID-19

The COVID-19 disease caused by the SARS-CoV-2 virus is a health crisis worldwide. While developing novel drugs and vaccines is long, repurposing existing drugs against COVID-19 can yield treatments with known preclinical, pharmacokinetic, pharmacodynamic, and toxicity profiles, which can rapidly enter clinical trials. In this study, we present a novel network-based drug repurposing platform to identify candidates for the treatment of COVID-19. At the time of the initial outbreak, knowledge about SARS-CoV-2 was lacking, but based on its similarity with other viruses, we sought to identify repurposing candidates to be tested rapidly at the clinical or preclinical levels. We first analyzed the genome sequence of SARS-CoV-2 and confirmed SARS as the closest virus by genome similarity, followed by MERS and other human coronaviruses. Using text mining and database searches, we obtained 34 COVID-19-related genes to seed the construction of a molecular network where our module detection and drug prioritization algorithms identified 24 disease-related human pathways, five modules, and 78 drugs to repurpose. Based on clinical knowledge, we re-prioritized 30 potentially repurposable drugs against COVID-19 (including pseudoephedrine, andrographolide, chloroquine, abacavir, and thalidomide). Our work shows how in silico repurposing analyses can yield testable candidates to accelerate the response to novel disease outbreaks.

Survey of Similarity-Based Prediction of Drug-Protein Interactions

Therapeutic activity of a significant majority of drugs is determined by their interactions with proteins. Databases of drug-protein interactions (DPIs) primarily focus on the therapeutic protein targets while the knowledge of the off-targets is fragmented and partial. One way to bridge this knowledge gap is to employ computational methods to predict protein targets for a given drug molecule, or interacting drugs for given protein targets. We survey a comprehensive set of 35 methods that were published in high-impact venues and that predict DPIs based on similarity between drugs and similarity between protein targets. We analyze the internal databases of known PDIs that these methods utilize to compute similarities, and investigate how they are linked to the 12 publicly available source databases. We discuss contents, impact and relationships between these internal and source databases, and well as the timeline of their releases and publications. The 35 predictors exploit and often combine three types of similarities that consider drug structures, drug profiles, and target sequences. We review the predictive architectures of these methods, their impact, and we explain how their internal DPIs databases are linked to the source databases. We also include a detailed timeline of the development of these predictors and discuss the underlying limitations of the current resources and predictive tools. Finally, we provide several recommendations concerning the future development of the related databases and methods.

Mono-Alkylated Ligands Based on Pyrazole and Triazole Derivatives Tested Against Fusarium oxysporum f. sp. albedinis: Synthesis, Characterization, DFT, and Phytase Binding Site Identification Using Blind Docking/Virtual Screening for Potent Fophy Inhibitors

Twelve recent compounds, incorporating several heterocyclic moieties such as pyrazole, thiazole, triazole, and benzotriazole, made in excellent yield up to 37–99.6%. They were tested against Fusarium oxysporum f. sp. albedinis fungi (Bayoud disease), where the best results are for compounds 2, 4, and 5 with IC₅₀ = 18.8–54.4 μg/mL. Density functional theory (DFT) study presented their molecular reactivity, while the docking simulations to describe the synergies between the trained compounds of dataset containing all the tested compounds (57 molecules) and F. oxysporum phytase domain (Fophy) enzyme as biological target. By comparing the results of the docking studies for the Fophy protein, it is found that compound 5 has the best affinity followed by compounds 2 and 4, so there is good agreement with the experimental results where their IC₅₀ values are in the following order: 74.28 (5) < 150 (2) < 214.10 (4), using Blind docking/virtual screening of the homology modeled protein and two different tools as Autodock Vina and Dockthor web tool that gave us predicted sites for further antifungal drug design.

An ABSINTH-Based Protocol for Predicting Binding Affinities between Proteins and Small Molecules

The core task in computational drug discovery is to accurately predict binding free energies in receptor-ligand systems for large libraries of putative binders. Here, the ABSINTH implicit solvent model and force field are extended to describe small, organic molecules and their interactions with proteins. We show that an automatic pipeline based on partitioning arbitrary molecules into substructures corresponding to model compounds with known free energies of solvation can be combined with the CHARMM general force field into a method that is successful at the two important challenges a scoring function faces in virtual screening work flows: it ranks known binders with correlation values rivaling that of comparable state-of-the-art methods and it enriches true binders in a set of decoys. Our protocol introduces innovative modifications to common virtual screening workflows, notably the use of explicit ions as competitors and the integration over multiple protein and ligand species differing in their protonation states. We demonstrate the value of modifications to both the protocol and ABSINTH itself. We conclude by discussing the limitations of high-throughput implicit methods such as the one proposed here.

Deciphering inhibitory activity of flavonoids against tau protein kinases: a coupled molecular docking and quantum chemical study

Today, Alzheimer's disease (AD) is one of the most important neurodegenerative disorders that affected millions of people worldwide. Hundreds of academic investigations highlighted the potential roles of natural metabolites in the cornerstone of AD prevention. Nevertheless, alkaloids are only metabolites that successfully showed promising clinical therapeutic effects on the prevention of AD. In this regard, other plant metabolites such as flavonoids are also considered as promising substances in the improvement of AD complications. The lack of data on molecular mode of action of flavonoids inside brain tissues, and their potential to transport across the blood-brain barrier, a physical hindrance between bloodstream and brain tissues, limited the large-scale application of these compounds for AD therapy programs. Herein, a coupled docking and quantum study was applied to determine the binding mode of flavonoids and three protein kinases involved in the pathogenesis of AD. The results suggested that all docked metabolites showed considerable binding affinity to interact with target receptors, but some compounds possessed higher binding energy values. Because docking simulation cannot entirely reveal the potential roles of ligand substructures in the interaction with target residues, quantum chemical analyses (QCAs) were performed to cover this drawback. Accordingly, QCAs determined that distribution of molecular orbitals have a pivotal function in the determination of the type of reaction between ligands and receptors; therefore, using such quantum chemical descriptors may correct the results of virtual docking outcomes to highlight promising backbones for further developments.Communicated by Ramaswamy H. Sarma.

EasyAmber: A comprehensive toolbox to automate the molecular dynamics simulation of proteins

Conformational plasticity of the functionally important regions and binding sites in protein/enzyme structures is one of the key factors affecting their function and interaction with substrates/ligands. Molecular dynamics (MD) can address the challenge of accounting for protein flexibility by predicting the time-dependent behavior of a molecular system. It has a potential of becoming a particularly important tool in protein engineering and drug discovery, but requires specialized training and skills, what impedes practical use by many investigators. We have developed the easyAmber - a comprehensive set of programs to automate the molecular dynamics routines implemented in the Amber package. The toolbox can address a wide set of tasks in computational biology struggling to account for protein flexibility. The automated workflow includes a complete set of steps from the initial "static" molecular model to the MD "production run": the full-atom model building, optimization/equilibration of the molecular system, classical/conventional and accelerated molecular dynamics simulations. The easyAmber implements advanced MD protocols, but is highly automated and easy-to-operate to attract a broad audience. The toolbox can be used on a personal desktop station equipped with a compatible gaming GPU-accelerator, as well as help to manage huge workloads on a powerful supercomputer. The software provides an opportunity to operate multiple simulations of different proteins at the same time, thus significantly increasing work efficiency. The easyAmber takes the molecular dynamics to the next level in terms of usability for complex processing of large volumes of data, thus supporting the recent trend away from inefficient "static" approaches in biology toward a deeper understanding of the dynamics in protein structures. The software is freely available for download at https://biokinet.belozersky.msu.ru/easyAmber, no login required.

The Advent of Generative Chemistry

Generative adversarial networks (GANs), first published in 2014, are among the most important concepts in modern artificial intelligence (AI). Bridging deep learning and game theory, GANs are used to generate or "imagine" new objects with desired properties. Since 2016, multiple GANs with reinforcement learning (RL) have been successfully applied in pharmacology for de novo molecular design. Those techniques aim at a more efficient use of the data and a better exploration of the chemical space. We review recent advances for the generation of novel molecules with desired properties with a focus on the applications of GANs, RL, and related techniques. We also discuss the current limitations and challenges in the new growing field of generative chemistry.

Anti-cancer potential of (1,2-dihydronaphtho[2,1-b]furan-2-yl)methanone derivatives

A series of 1,2-dihydronaphtho[2,1-b]furan derivatives were synthesized by cyclizing 1-(aryl/alkyl(arylthio)methyl)-naphthalen-2-ol and pyridinium bromides in the presence of 1,8-diazabicyclo[5.4.0]undec-7-ene (DBU) in very good yield. The synthesized compounds were evaluated for their anti-proliferative potential against human triple negative MDA-MB-468 and MCF-7 breast cancer cells and non-cancerous WI-38 cells (lung fibroblast cell) using MTT experiments. Among 21 synthesized compounds, three compounds (3a, 3b and 3 s) showed promising anti-cancer potential and compound 3b was found to have best anti-proliferative activities based on the results of several biochemical and microscopic experiments.

The impact of compound library size on the performance of scoring functions for structure-based virtual screening

Larger training datasets have been shown to improve the accuracy of machine learning (ML)-based scoring functions (SFs) for structure-based virtual screening (SBVS). In addition, massive test sets for SBVS, known as ultra-large compound libraries, have been demonstrated to enable the fast discovery of selective drug leads with low-nanomolar potency. This proof-of-concept was carried out on two targets using a single docking tool along with its SF. It is thus unclear whether this high level of performance would generalise to other targets, docking tools and SFs. We found that screening a larger compound library results in more potent actives being identified in all six additional targets using a different docking tool along with its classical SF. Furthermore, we established that a way to improve the potency of the retrieved molecules further is to rank them with more accurate ML-based SFs (we found this to be true in four of the six targets; the difference was not significant in the remaining two targets). A 3-fold increase in average hit rate across targets was also achieved by the ML-based SFs. Lastly, we observed that classical and ML-based SFs often find different actives, which supports using both types of SFs on those targets.

In-Silico Modeling in Drug Metabolism and Interaction: Current Strategies of Lead Discovery

BACKGROUND

Drug metabolism is a complex mechanism of human body systems to detoxify foreign particles, chemicals, and drugs through bio alterations. It involves many biochemical reactions carried out by invivo enzyme systems present in the liver, kidney, intestine, lungs, and plasma. After drug administration, it crosses several biological membranes to reach into the target site for binding and produces the therapeutic response. After that, it may undergo detoxification and excretion to get rid of the biological systems. Most of the drugs and its metabolites are excreted through kidney via urination. Some drugs and their metabolites enter into intestinal mucosa and excrete through feces. Few of the drugs enter into hepatic circulation where they go into the intestinal tract. The drug leaves the liver via the bile duct and is excreted through feces. Therefore, the study of total methodology of drug biotransformation and interactions with various targets is costly.

METHODS

To minimize time and cost, in-silico algorithms have been utilized for lead-like drug discovery. Insilico modeling is the process where a computer model with a suitable algorithm is developed to perform a controlled experiment. It involves the combination of both in-vivo and in-vitro experimentation with virtual trials, eliminating the non-significant variables from a large number of variable parameters. Whereas, the major challenge for the experimenter is the selection and validation of the preferred model, as well as precise simulation in real physiological status.

RESULTS

The present review discussed the application of in-silico models to predict absorption, distribution, metabolism, and excretion (ADME) properties of drug molecules and also access the net rate of metabolism of a compound.

CONCLUSION

It helps with the identification of enzyme isoforms; which are likely to metabolize a compound, as well as the concentration dependence of metabolism and the identification of expected metabolites. In terms of drug-drug interactions (DDIs), models have been described for the inhibition of metabolism of one compound by another, and for the compound-dependent induction of drug-metabolizing enzymes.

Plectrabarbene, a New Abietane Diterpene from Plectranthus barbatus Aerial Parts

A new abietane diterpene namely plectrabarbene (2), together with two known compounds: sugiol (1) and 11,14-dihydroxy-8,11,13-abietatrien-7-one (3) have been isolated from the aerial parts of Plectranthus barbatus Andr. (Labiatae). The structures of these compounds were determined by various spectral techniques (e.g., UV, IR, NMR, and FAB) and by comparison with the literature data. A molecular docking study of the isolated diterpenes (1–3) was performed with AChE to gain an insight into their AChE inhibition mechanism. The results of docking experiments revealed that the all tested compounds showed binding affinity at the active site of AchE in comparison to donepezil.

Unraveling Plant Natural Chemical Diversity for Drug Discovery Purposes

The screening and testing of extracts against a variety of pharmacological targets in order to benefit from the immense natural chemical diversity is a concern in many laboratories worldwide. And several successes have been recorded in finding new actives in natural products, some of which have become new drugs or new sources of inspiration for drugs. But in view of the vast amount of research on the subject, it is surprising that not more drug candidates were found. In our view, it is fundamental to reflect upon the approaches of such drug discovery programs and the technical processes that are used, along with their inherent difficulties and biases. Based on an extensive survey of recent publications, we discuss the origin and the variety of natural chemical diversity as well as the strategies to having the potential to embrace this diversity. It seemed to us that some of the difficulties of the area could be related with the technical approaches that are used, so the present review begins with synthetizing some of the more used discovery strategies, exemplifying some key points, in order to address some of their limitations. It appears that one of the challenges of natural product-based drug discovery programs should be an easier access to renewable sources of plant-derived products. Maximizing the use of the data together with the exploration of chemical diversity while working on reasonable supply of natural product-based entities could be a way to answer this challenge. We suggested alternative ways to access and explore part of this chemical diversity with in vitro cultures. We also reinforced how important it was organizing and making available this worldwide knowledge in an "inventory" of natural products and their sources. And finally, we focused on strategies based on synthetic biology and syntheses that allow reaching industrial scale supply. Approaches based on the opportunities lying in untapped natural plant chemical diversity are also considered.

Molecular Functionality of Plant Proteins from Low- to High-Solid Systems with Ligand and Co-Solute

In the food industry, proteins are regarded as multifunctional systems whose bioactive hetero-polymeric properties are affected by physicochemical interactions with the surrounding components in formulations. Due to their nutritional value, plant proteins are increasingly considered by the new product developer to provide three-dimensional assemblies of required structure, texture, solubility and interfacial/bulk stability with physical, chemical or enzymatic treatment. This molecular flexibility allows them to form systems for the preservation of fresh food, retention of good nutrition and interaction with a range of microconstituents. While, animal- and milk-based proteins have been widely discussed in the literature, the role of plant proteins in the development of functional foods with enhanced nutritional profile and targeted physiological effects can be further explored. This review aims to look into the molecular functionality of plant proteins in relation to the transport of bioactive ingredients and interaction with other ligands and proteins. In doing so, it will consider preparations from low- to high-solids and the effect of structural transformation via gelation, phase separation and vitrification on protein functionality as a delivery vehicle or heterologous complex. Applications for the design of novel functional foods and nutraceuticals will also be discussed.

Discovery of potent apoptosis signal-regulating kinase 1 inhibitors via integrated computational strategy and biological evaluation

Apoptosis signal-regulating Kinase 1 (ASK1) has been confirmed as a potential therapeutic target for the treatment of non-alcoholic steatohepatitis (NASH) disorder and the discovery of ASK1 inhibitors has attracted increasing attention. In this work, a series of in silico methods including pharmacophore screening, docking binding site analysis, protein-ligand interaction fingerprint (PLIF) similarity investigation and molecular docking were applied to find the potential hits from commercial compound databases. Five compounds with potential inhibitory activity were purchased and submitted to biological activity validation. Thus, one hit compound was discovered with micromolar IC₅₀ value (10.59 μM) against ASK1. Results demonstrated that the integration of computation methods and biological test was quite reliable for the discovery of potent ASK1 inhibitors and the strategy could be extended to other similar targets of interest.

Supercomputer docking with a large number of degrees of freedom

Docking represents one of the most popular computational approaches in drug design. It has reached popularity owing to capability of identifying correct conformations of a ligand within an active site of the target-protein and of estimating the binding affinity of a ligand that is immensely helpful in prediction of compound activity. Despite many success stories, there are challenges, in particular, handling with a large number of degrees of freedom in solving the docking problem. Here, we show that SOL-P, the docking program based on the new Tensor Train algorithm, is capable to dock successfully oligopeptides having up to 25 torsions. To make the study comparative we have performed docking of the same oligopeptides with the SOL program which uses the same force field as that utilized by SOL-P and has common features of many docking programs: the genetic algorithm of the global optimization and the grid approximation. SOL has managed to dock only one oligopeptide. Moreover, we present the results of docking with SOL-P ligands into proteins with moveable atoms. Relying on visual observations we have determined the common protein atom groups displaced after docking which seem to be crucial for successful prediction of experimental conformations of ligands.

Recent applications of deep learning and machine intelligence on in silico drug discovery: methods, tools and databases

The identification of interactions between drugs/compounds and their targets is crucial for the development of new drugs. In vitro screening experiments (i.e. bioassays) are frequently used for this purpose; however, experimental approaches are insufficient to explore novel drug-target interactions, mainly because of feasibility problems, as they are labour intensive, costly and time consuming. A computational field known as 'virtual screening' (VS) has emerged in the past decades to aid experimental drug discovery studies by statistically estimating unknown bio-interactions between compounds and biological targets. These methods use the physico-chemical and structural properties of compounds and/or target proteins along with the experimentally verified bio-interaction information to generate predictive models. Lately, sophisticated machine learning techniques are applied in VS to elevate the predictive performance. The objective of this study is to examine and discuss the recent applications of machine learning techniques in VS, including deep learning, which became highly popular after giving rise to epochal developments in the fields of computer vision and natural language processing. The past 3 years have witnessed an unprecedented amount of research studies considering the application of deep learning in biomedicine, including computational drug discovery. In this review, we first describe the main instruments of VS methods, including compound and protein features (i.e. representations and descriptors), frequently used libraries and toolkits for VS, bioactivity databases and gold-standard data sets for system training and benchmarking. We subsequently review recent VS studies with a strong emphasis on deep learning applications. Finally, we discuss the present state of the field, including the current challenges and suggest future directions. We believe that this survey will provide insight to the researchers working in the field of computational drug discovery in terms of comprehending and developing novel bio-prediction methods.

Application of the SwissDrugDesign Online Resources in Virtual Screening

SwissDrugDesign is an important initiative led by the Molecular Modeling Group of the SIB Swiss Institute of Bioinformatics. This project provides a collection of freely available online tools for computer-aided drug design. Some of these web-based methods, i.e., SwissSimilarity and SwissTargetPrediction, were especially developed to perform virtual screening, while others such as SwissADME, SwissDock, SwissParam and SwissBioisostere can find applications in related activities. The present review aims at providing a short description of these methods together with examples of their application in virtual screening, where SwissDrugDesign tools successfully supported the discovery of bioactive small molecules.

Natural Alkaloid Berberine Activity against Pseudomonas aeruginosa MexXY-Mediated Aminoglycoside Resistance: In Silico and in Vitro Studies

The multidrug efflux system MexXY-OprM, inside the resistance-nodulation-division family, is a major determinant of aminoglycoside resistance in Pseudomonas aeruginosa. In the fight aimed to identify potential efflux pump inhibitors among natural compounds, the alkaloid berberine emerged as a putative inhibitor of MexXY-OprM. In this work, we elucidated its interaction with the extrusor protein MexY and assessed its synergistic activity with aminoglycosides. In particular, we built an in silico model for the MexY protein in its trimeric association using both AcrB (E. coli) and MexB (P. aeruginosa) as 3D templates. This model has been stabilized in the bacterial cytoplasmic membrane using a molecular dynamics approach and used for ensemble docking to obtain the binding site mapping. Then, through dynamic docking, we assessed its binding affinity and its synergism with aminoglycosides focusing on tobramycin, which is widely used in the treatment of pulmonary infections. In vitro assays validated the data obtained: the results showed a 2-fold increase of the inhibitory activity and 2-4 log increase of the killing activity of the association berberine-tobramycin compared to those of tobramycin alone against 13/28 tested P. aeruginosa clinical isolates. From hemolytic assays, we preliminarily assessed berberine's low toxicity.

Discovery of an intrasubunit nicotinic acetylcholine receptor-binding site for the positive allosteric modulator Br-PBTC

Nicotinic acetylcholine receptor (nAChR) ligands that lack agonist activity but enhance activation in the presence of an agonist are called positive allosteric modulators (PAMs). nAChR PAMs have therapeutic potential for the treatment of nicotine addiction and several neuropsychiatric disorders. PAMs need to be selectively targeted toward certain nAChR subtypes to tap this potential. We previously discovered a novel PAM, (R)-7-bromo-N-(piperidin-3-yl)benzo[b]thiophene-2-carboxamide (Br-PBTC), which selectively potentiates the opening of α4β2*, α2β2*, α2β4*, and (α4β4)₂α4 nAChRs and reactivates some of these subtypes when desensitized (* indicates the presence of other subunits). We located the Br-PBTC-binding site through mutagenesis and docking in α4. The amino acids Glu-282 and Phe-286 near the extracellular domain on the third transmembrane helix were found to be crucial for Br-PBTC's PAM effect. E282Q abolishes Br-PBTC potentiation. Using (α4^E282Qβ2)₂α5 nAChRs, we discovered that the trifluoromethylated derivatives of Br-PBTC can potentiate channel opening of α5-containing nAChRs. Mutating Tyr-430 in the α5 M4 domain changed α5-selectivity among Br-PBTC derivatives. There are two kinds of α4 subunits in α4β2 nAChRs. Primary α4 forms an agonist-binding site with another β2 subunit. Accessory α4 forms an agonist-binding site with another α4 subunit. The pharmacological effect of Br-PBTC depends both on its own and agonists' occupancy of primary and accessory α4 subunits. Br-PBTC reactivates desensitized (α4β2)₂α4 nAChRs. Its full efficacy requires intact Br-PBTC sites in at least one accessory and one primary α4 subunit. PAM potency increases with higher occupancy of the agonist sites. Br-PBTC and its derivatives should prove useful as α subunit-selective nAChR PAMs.

Accelerating Drug Discovery by Early Protein Drug Target Prediction Based on a Multi-Fingerprint Similarity Search

In this continuing work, we have updated our recently proposed Multi-fingerprint Similarity Search algorithm (MuSSel) by enabling the generation of dominant ionized species at a physiological pH and the exploration of a larger data domain, which included more than half a million high-quality small molecules extracted from the latest release of ChEMBL (version 24.1, at the time of writing). Provided with a high biological assay confidence score, these selected compounds explored up to 2822 protein drug targets. To improve the data accuracy, samples marked as prodrugs or with equivocal biological annotations were not considered. Notably, MuSSel performances were overall improved by using an object-relational database management system based on PostgreSQL. In order to challenge the real effectiveness of MuSSel in predicting relevant therapeutic drug targets, we analyzed a pool of 36 external bioactive compounds published in the Journal of Medicinal Chemistry from October to December 2018. This study demonstrates that the use of highly curated chemical and biological experimental data on one side, and a powerful multi-fingerprint search algorithm on the other, can be of the utmost importance in addressing the fate of newly conceived small molecules, by strongly reducing the attrition of early phases of drug discovery programs.

AMMOS2: a web server for protein-ligand-water complexes refinement via molecular mechanics

AMMOS2 is an interactive web server for efficient computational refinement of protein-small organic molecule complexes. The AMMOS2 protocol employs atomic-level energy minimization of a large number of experimental or modeled protein-ligand complexes. The web server is based on the previously developed standalone software AMMOS (Automatic Molecular Mechanics Optimization for in silico Screening). AMMOS utilizes the physics-based force field AMMP sp4 and performs optimization of protein-ligand interactions at five levels of flexibility of the protein receptor. The new version 2 of AMMOS implemented in the AMMOS2 web server allows the users to include explicit water molecules and individual metal ions in the protein-ligand complexes during minimization. The web server provides comprehensive analysis of computed energies and interactive visualization of refined protein-ligand complexes. The ligands are ranked by the minimized binding energies allowing the users to perform additional analysis for drug discovery or chemical biology projects. The web server has been extensively tested on 21 diverse protein-ligand complexes. AMMOS2 minimization shows consistent improvement over the initial complex structures in terms of minimized protein-ligand binding energies and water positions optimization. The AMMOS2 web server is freely available without any registration requirement at the URL: http://drugmod.rpbs.univ-paris-diderot.fr/ammosHome.php.

Chiral Graphs: Reduced Representations of Ligand Scaffolds for Stereoselective Biomolecular Recognition, Drug Design, and Enhanced Exploration of Chemical Structure Space

Rational structure-based drug design relies on a detailed, atomic-level understanding of protein-ligand interactions. The chiral nature of drug binding sites in proteins has led to the discovery of predominantly chiral drugs. A mechanistic understanding of stereoselectivity (which governs how one stereoisomer of a drug might bind stronger than the others to a protein) depends on the topology of stereocenters in the chiral molecule. Chiral graphs and reduced chiral graphs, introduced here, are new topological representations of chiral ligands using graph theory, to facilitate a detailed understanding of chiral recognition of ligands/drugs by proteins. These representations are demonstrated by application to all ≈14 000+ chiral ligands in the Protein Data Bank (PDB), which will facilitate an understanding of protein-ligand stereoselectivity mechanisms. Ligand modifications during drug development can be easily incorporated into these chiral graphs. In addition, these chiral graphs present an efficient tool for a deep dive into the enormous chemical structure space to enable sampling of unexplored structural scaffolds.

Development and Testing of Druglike Screening Libraries

Although significant advances in experimental high throughput screening (HTS) have been made for drug lead identification, in silico virtual screening (VS) is indispensable owing to its unique advantage over experimental HTS, target-focused, cheap, and efficient, albeit its disadvantage of producing false positive hits. For both experimental HTS and VS, the quality of screening libraries is crucial and determines the outcome of those studies. In this paper, we first reviewed the recent progress on screening library construction. We realized the urgent need for compiling high-quality screening libraries in drug discovery. Then we compiled a set of screening libraries from about 20 million druglike ZINC molecules by running fingerprint-based similarity searches against known drug molecules. Lastly, the screening libraries were objectively evaluated using 5847 external actives covering more than 2000 drug targets. The result of the assessment is very encouraging. For example, with the Tanimoto coefficient being set to 0.75, 36% of external actives were retrieved and the enrichment factor was 13. Additionally, drug target family specific screening libraries were also constructed and evaluated. The druglike screening libraries are available for download from https://mulan.pharmacy.pitt.edu .

A review of ligand-based virtual screening web tools and screening algorithms in large molecular databases in the age of big data

Virtual screening has become a widely used technique for helping in drug discovery processes. The key to this success is its ability to aid in the identification of novel bioactive compounds by screening large molecular databases. Several web servers have emerged in the last few years supplying platforms to guide users in screening publicly accessible chemical databases in a reasonable time. In this review, we discuss a representative set of online virtual screening servers and their underlying similarity algorithms. Other related topics, such as molecular representation or freely accessible databases are also treated. The most relevant contributions to this review arise from critical discussions concerning the pros and cons of servers and algorithms, and the challenges that future works must solve in a virtual screening framework.

Protein structure and computational drug discovery

The first protein structures revealed a complex web of weak interactions stabilising the three-dimensional shape of the molecule. Small molecule ligands were then found to exploit these same weak binding events to modulate protein function or act as substrates in enzymatic reactions. As the understanding of ligand–protein binding grew, it became possible to firstly predict how and where a particular small molecule might interact with a protein, and then to identify putative ligands for a specific protein site. Computer-aided drug discovery, based on the structure of target proteins, is now a well-established technique that has produced several marketed drugs. We present here an overview of the various methodologies being used for structure-based computer-aided drug discovery and comment on possible future developments in the field.

Review and comparative assessment of similarity-based methods for prediction of drug–protein interactions in the druggable human proteome

Drug–protein interactions (DPIs) underlie the desired therapeutic actions and the adverse side effects of a significant majority of drugs. Computational prediction of DPIs facilitates research in drug discovery, characterization and repurposing. Similarity-based methods that do not require knowledge of protein structures are particularly suitable for druggable genome-wide predictions of DPIs. We review 35 high-impact similarity-based predictors that were published in the past decade. We group them based on three types of similarities and their combinations that they use. We discuss and compare key aspects of these methods including source databases, internal databases and their predictive models. Using our novel benchmark database, we perform comparative empirical analysis of predictive performance of seven types of representative predictors that utilize each type of similarity individually and all possible combinations of similarities. We assess predictive quality at the database-wide DPI level and we are the first to also include evaluation over individual drugs. Our comprehensive analysis shows that predictors that use more similarity types outperform methods that employ fewer similarities, and that the model combining all three types of similarities secures area under the receiver operating characteristic curve of 0.93. We offer a comprehensive analysis of sensitivity of predictive performance to intrinsic and extrinsic characteristics of the considered predictors. We find that predictive performance is sensitive to low levels of similarities between sequences of the drug targets and several extrinsic properties of the input drug structures, drug profiles and drug targets. The benchmark database and a webserver for the seven predictors are freely available at http://biomine.cs.vcu.edu/servers/CONNECTOR/.

Ligand-Based Pharmacophore Screening Strategy: a Pragmatic Approach for Targeting HER Proteins

Targeting ErbB family of receptors is an important therapeutic option, because of its essential role in the broad spectrum of human cancers, including non-small cell lung cancer (NSCLC). Therefore, in the present work, considerable effort has been made to develop an inhibitor against HER family proteins, by combining the use of pharmacophore modelling, docking scoring functions, and ADME property analysis. Initially, a five-point pharmacophore model was developed using known HER family inhibitors. The generated model was then used as a query to screen a total of 468,880 compounds of three databases namely ZINC, ASINEX, and DrugBank. Subsequently, docking analysis was carried out to obtain hit molecules that could inhibit the HER receptors. Further, analysis of GLIDE scores and ADME properties resulted in one hit namely BAS01025917 with higher glide scores, increased CNS involvement, and good pharmaceutically relevant properties than reference ligand, afatinib. Furthermore, the inhibitory activity of the lead compounds was validated by performing molecular dynamic simulations. Of note, BAS01025917 was found to possess scaffolds with a broad spectrum of antitumor activity. We believe that this novel hit molecule can be further exploited for the development of a pan-HER inhibitor with low toxicity and greater potential.

Strategies of Virtual Screening in Medicinal Chemistry

Virtual screening represents an effective computational strategy to rise-up the chances of finding new bioactive compounds by accelerating the time needed to move from an initial intuition to market. Classically, the most pursued approaches rely on ligand- and structure-based studies, the former employed when structural data information about the target is missing while the latter employed when X-ray/NMR solved or homology models are instead available for the target. The authors will focus on the most advanced techniques applied in this area. In particular, they will survey the key concepts of virtual screening by discussing how to properly select chemical libraries, how to make database curation, how to applying and- and structure-based techniques, how to wisely use post-processing methods. Emphasis will be also given to the most meaningful databases used in VS protocols. For the ease of discussion several examples will be presented.

Evaluation of Protein-Ligand Docking by Cyscore

Protein-ligand docking is a powerful method in drug discovery. The reliability of docking can be quantified by RMSD between a docking structure and an experimentally determined one. However, most experimentally determined structures are not available in practice. Evaluation by scoring functions is an alternative for assessing protein-ligand docking results. This chapter first provides a brief introduction to scoring methods used in docking. Then details are provided on how to use Cyscore programs. Finally it describes a case study for evaluation of protein-ligand docking.

Synthesis, SAR, and Docking Studies Disclose 2-Arylfuran-1,4-naphthoquinones as In Vitro Antiplasmodial Hits

A total of 28 lapachol-related naphthoquinones with four different scaffolds were synthesized and spectroscopically characterized. In vitro antiplasmodial activity was assayed against the chloroquine-resistant Plasmodium falciparum W2 strain by the parasite lactate dehydrogenase (pLDH) method. Cytotoxicity against Hep G2A16 cell was determined by the MTT assay. All compounds disclosed higher in vitro antiplasmodial activity than lapachol. Ortho- and para-naphthoquinones with a furan ring fused to the quinonoid moiety were more potent than 2-hydroxy-3-(1′-alkenyl)-1,4-naphthoquinones, while ortho-furanonaphthoquinones were more cytotoxic. Molecular docking to Plasmodium targets Pfcyt bc1 complex and PfDHOD enzyme showed that five out of the 28 naphthoquinones disclosed favorable binding energies. Furanonaphthoquinones endowed with an aryl moiety linked to the furan ring are highlighted as new in vitro antiplasmodial lead compounds and warrant further investigation.

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.