Exponential consensus ranking improves the outcome in docking and receptor ensemble docking

Building shape-focused pharmacophore models for effective docking screening

The performance of molecular docking can be improved by comparing the shape similarity of the flexibly sampled poses against the target proteins' inverted binding cavities. The effectiveness of these pseudo-ligands or negative image-based models in docking rescoring is boosted further by performing enrichment-driven optimization. Here, we introduce a novel shape-focused pharmacophore modeling algorithm O-LAP that generates a new class of cavity-filling models by clumping together overlapping atomic content via pairwise distance graph clustering. Top-ranked poses of flexibly docked active ligands were used as the modeling input and multiple alternative clustering settings were benchmark-tested thoroughly with five demanding drug targets using random training/test divisions. In docking rescoring, the O-LAP modeling typically improved massively on the default docking enrichment; furthermore, the results indicate that the clustered models work well in rigid docking. The C+ +/Qt5-based algorithm O-LAP is released under the GNU General Public License v3.0 via GitHub ( https://github.com/jvlehtonen/overlap-toolkit ). SCIENTIFIC CONTRIBUTION: This study introduces O-LAP, a C++/Qt5-based graph clustering software for generating new type of shape-focused pharmacophore models. In the O-LAP modeling, the target protein cavity is filled with flexibly docked active ligands, the overlapping ligand atoms are clustered, and the shape/electrostatic potential of the resulting model is compared against the flexibly sampled molecular docking poses. The O-LAP modeling is shown to ensure high enrichment in both docking rescoring and rigid docking based on comprehensive benchmark-testing.

From Proteome to Potential Drugs: Integration of Subtractive Proteomics and Ensemble Docking for Drug Repurposing against Pseudomonas aeruginosa RND Superfamily Proteins

Pseudomonas aeruginosa (P. aeruginosa) poses a significant threat as a nosocomial pathogen due to its robust resistance mechanisms and virulence factors. This study integrates subtractive proteomics and ensemble docking to identify and characterize essential proteins in P. aeruginosa, aiming to discover therapeutic targets and repurpose commercial existing drugs. Using subtractive proteomics, we refined the dataset to discard redundant proteins and minimize potential cross-interactions with human proteins and the microbiome proteins. We identified 12 key proteins, including a histidine kinase and members of the RND efflux pump family, known for their roles in antibiotic resistance, virulence, and antigenicity. Predictive modeling of the three-dimensional structures of these RND proteins and subsequent molecular ensemble-docking simulations led to the identification of MK-3207, R-428, and Suramin as promising inhibitor candidates. These compounds demonstrated high binding affinities and effective inhibition across multiple metrics. Further refinement using non-covalent interaction index methods provided deeper insights into the electronic effects in protein-ligand interactions, with Suramin exhibiting superior binding energies, suggesting its broad-spectrum inhibitory potential. Our findings confirm the critical role of RND efflux pumps in antibiotic resistance and suggest that MK-3207, R-428, and Suramin could be effectively repurposed to target these proteins. This approach highlights the potential of drug repurposing as a viable strategy to combat P. aeruginosa infections.

Docking Score ML: Target-Specific Machine Learning Models Improving Docking-Based Virtual Screening in 155 Targets

In drug discovery, molecular docking methods face challenges in accurately predicting energy. Scoring functions used in molecular docking often fail to simulate complex protein-ligand interactions fully and accurately leading to biases and inaccuracies in virtual screening and target predictions. We introduce the "Docking Score ML", developed from an analysis of over 200,000 docked complexes from 155 known targets for cancer treatments. The scoring functions used are founded on bioactivity data sourced from ChEMBL and have been fine-tuned using both supervised machine learning and deep learning techniques. We validated our approach extensively using multiple data sets such as validation of selectivity mechanism, the DUDE, DUD-AD, and LIT-PCBA data sets, and performed a multitarget analysis on drugs like sunitinib. To enhance prediction accuracy, feature fusion techniques were explored. By merging the capabilities of the Graph Convolutional Network (GCN) with multiple docking functions, our results indicated a clear superiority of our methodologies over conventional approaches. These advantages demonstrate that Docking Score ML is an efficient and accurate tool for virtual screening and reverse docking.

PROTACable Is an Integrative Computational Pipeline of 3-D Modeling and Deep Learning To Automate the De Novo Design of PROTACs

Proteolysis-targeting chimeras (PROTACs) that engage two biological targets at once are a promising technology in degrading clinically relevant protein targets. Since factors that influence the biological activities of PROTACs are more complex than those of a small molecule drug, we explored a combination of computational chemistry and deep learning strategies to forecast PROTAC activity and enable automated design. A new method named PROTACable was developed for the de novo design of PROTACs, which includes a robust 3-D modeling workflow to model PROTAC ternary complexes using a library of E3 ligase and linker and an SE(3)-equivariant graph transformer network to predict the activity of newly designed PROTACs. PROTACable is available at https://github.com/giaguaro/PROTACable/.

Benchmarking ANI potentials as a rescoring function and screening FDA drugs for SARS-CoV-2 Mpro

Here, we introduce the use of ANI-ML potentials as a rescoring function in the host-guest interaction in molecular docking. Our results show that the "docking power" of ANI potentials can compete with the current scoring functions at the same level of computational cost. Benchmarking studies on CASF-2016 dataset showed that ANI is ranked in the top 5 scoring functions among the other 34 tested. In particular, the ANI predicted interaction energies when used in conjunction with GOLD-PLP scoring function can boost the top ranked solution to be the closest to the x-ray structure. Rapid and accurate calculation of interaction energies between ligand and protein also enables screening of millions of drug candidates/docking poses. Using a unique protocol in which docking by GOLD-PLP, rescoring by ANI-ML potentials and extensive MD simulations along with end state free energy methods are combined, we have screened FDA approved drugs against the SARS-CoV-2 main protease (Mpro). The top six drug molecules suggested by the consensus of these free energy methods have already been in clinical trials or proposed as potential drug molecules in previous theoretical and experimental studies, approving the validity and the power of accuracy in our screening method.

ESSENCE-Dock: A Consensus-Based Approach to Enhance Virtual Screening Enrichment in Drug Discovery

Drug development is a complex, costly, and time-consuming endeavor. While high-throughput screening (HTS) plays a critical role in the discovery stage, it is one of many factors contributing to these challenges. In certain contexts, virtual screening can complement the HTS, potentially offering a more streamlined approach in the initial stages of drug discovery. Molecular docking is an example of a popular virtual screening technique that is often used for this purpose; however, its effectiveness can vary greatly. This has led to the use of consensus docking approaches that combine results from different docking methods to improve the identification of active compounds and reduce the occurrence of false positives. However, many of these methods do not fully leverage the latest advancements in molecular docking. In response, we present ESSENCE-Dock (Effective Structural Screening ENrichment ConsEnsus Dock), a new consensus docking workflow aimed at decreasing false positives and increasing the discovery of active compounds. By utilizing a combination of novel docking algorithms, we improve the selection process for potential active compounds. ESSENCE-Dock has been made to be user-friendly, requiring only a few simple commands to perform a complete screening while also being designed for use in high-performance computing (HPC) environments.

Consensus screening for a challenging target: the quest for P-glycoprotein inhibitors

ATP-binding cassette (ABC) transporters are a large family of proteins involved in membrane transport of a wide variety of substrates. Among them, ABCB1, also known as MDR-1 or P-glycoprotein (P-gp), is the most characterized. By exporting xenobiotics out of the cell, P-gp activity can affect the ADME properties of several drugs. Moreover, P-gp has been found to mediate multidrug resistance in cancer cells. Thus, the inhibition of P-gp activity may lead to increased absorption and/or intracellular accumulation of co-administered drugs, enhancing their effectiveness. Using the human-mouse chimeric cryoEM 3D structure of the P-gp in the inhibitor-bound intermediate form (PDBID: 6qee), approximately 200 000 commercially available natural compounds from the ZINC database were virtually screened. To build a model able to discriminate between substrate and inhibitors, two datasets of compounds with known activity, including P-gp inhibitors, substrates, and inactive molecules were also docked. The best docking pose of selected substrates and inhibitors were used to generate 3D common feature pharmacophoric models that were combined with the Autodock Vina binding energy values to prioritize compounds for visual inspection. With this consensus approach, 13 potential candidates were identified and then tested for their ability to inhibit P-gp, using zosuquidar, a third generation P-gp inhibitor, as a reference drug. Eight compounds were found to be active with 6 of them having an IC50 lower than 5 μM in a membrane-based ATPase activity assay. Moreover, the P-gp inhibitory activity was also confirmed by two different cell-based in vitro methods. Both retrospective and prospective results demonstrate the ability of the combined structure-based pharmacophore modeling and docking-based virtual screening approach to predict novel hit compounds with inhibitory activity toward P-gp. The resulting chemical scaffolds could serve as inspiration for the optimization of novel and more potent P-gp inhibitors.

Susceptibility of Tetranychus urticae to the Alkaloidal Extract of Zanthoxylum schreberi Bark: Phenotypic and Biochemical Insights for Biotechnological Exploitation

Tetranychus urticae Koch, a phytophagous mite, is one of the most significant crop pests globally. The primary method employed for controlling T. urticae involves chemical means, utilizing synthesized products, posing the risk of developing resistance. The urgency for novel strategies integrated into pest management programs to combat this mite is becoming increasingly imperative. Botanical pesticides emerge as a promising tool to forestall arthropod resistance. Among these, extracts from Rutaceae plants, abundant in bioactive specialized metabolites, have demonstrated potential as insecticides and miticides. In this study, various concentrations of alkaloidal extracts sourced from the bark of Zanthoxylum schreberi J.F.Gmel. (Rutaceae) were evaluated against T. urticae adult females. Furthermore, the extract's combination with three distinct commercial acaricides (i.e., chlorfenapyr, cyflumetofen, and abamectin) was also assessed for this mite. Chemical characterization of the extract via LC-MS allowed for the annotation of various compounds related to ten benzylisoquinoline-derived alkaloids. The extract, both alone and in combination with commercial insecticides, yielded varying responses, inducing over 40% mortality at 2% w/w, demonstrating a 90% repellency rate at the same concentration, and exerting a moderate impact on fecundity. These treatments extended beyond phenotypic responses, delving into the biochemical effects on treated T. urticae females through an exploration of the impact on four enzymes, i.e., acetylcholinesterase (AChE), glutathione S-transferase (GST), esterases (GE), and P450-like monooxygenases (PMO). Employing consensus docking studies and in vitro enzymatic evaluations, it was discovered that the Z. schreberi-derived extract and its constituents significantly affected two key enzymes, AChE and GST (IC50 < 6 µM), which were associated with the phenotypic observations of T. urticae females. The evaluation of alkaloid-rich botanicals showcases promising potential as a relevant biotechnological strategy in addressing mite-related concerns, offering a pathway toward innovative and sustainable pest management solutions.

Integrating QSAR modelling and deep learning in drug discovery: the emergence of deep QSAR

Quantitative structure-activity relationship (QSAR) modelling, an approach that was introduced 60 years ago, is widely used in computer-aided drug design. In recent years, progress in artificial intelligence techniques, such as deep learning, the rapid growth of databases of molecules for virtual screening and dramatic improvements in computational power have supported the emergence of a new field of QSAR applications that we term 'deep QSAR'. Marking a decade from the pioneering applications of deep QSAR to tasks involved in small-molecule drug discovery, we herein describe key advances in the field, including deep generative and reinforcement learning approaches in molecular design, deep learning models for synthetic planning and the application of deep QSAR models in structure-based virtual screening. We also reflect on the emergence of quantum computing, which promises to further accelerate deep QSAR applications and the need for open-source and democratized resources to support computer-aided drug design.

In Search of Novel SGLT2 Inhibitors by High-throughput Virtual Screening

BACKGROUND

Type 2 diabetes mellitus constitutes approximately 90% of all reported forms of diabetes mellitus. Insulin resistance characterizes this manifestation of diabetes. The prevalence of this condition is commonly observed in patients aged 45 and above; however, there is an emerging pattern of younger cohorts receiving diagnoses primarily attributed to lifestyle-related variables, including obesity, sedentary behavior, and poor dietary choices. The enzyme SGLT2 exerts a negative regulatory effect on insulin signaling pathways, resulting in the development of insulin resistance and subsequent elevation of blood glucose levels. The maintenance of glucose homeostasis relies on the proper functioning of insulin signaling pathways, while disruptions in insulin signaling can contribute to the development of type 2 diabetes.

OBJECTIVE

Our study aimed to identify novel SGLT2 inhibitors by high-throughput virtual Screening.

METHODS

We screened the May bridge Hit Discover database to identify potent hits followed by druglikeness, synthetic accessibility, PAINS alert, toxicity estimation, ADME assessment, and consensus molecular docking.

RESULTS

The screening process led to the identification of three molecules that demonstrated significant binding affinity, favorable drug-like properties, effective ADME, and minimal toxicity.

CONCLUSION

The identified molecules could manage T2DM effectively by inhibiting SGLT2, providing a promising avenue for future therapeutic strategies.

Synthesis of Pyrrolo[3,4-b]pyridin-5-ones via Ugi-Zhu Reaction and In Vitro-In Silico Studies against Breast Carcinoma

An Ugi-Zhu three-component reaction (UZ-3CR) coupled in a one-pot manner to a cascade process (N-acylation/aza Diels-Alder cycloaddition/decarboxylation/dehydration) was performed to synthesize a series of pyrrolo[3,4-b]pyridin-5-ones in 20% to 92% overall yields using ytterbium triflate as a catalyst, toluene as a solvent, and microwaves as a heat source. The synthesized molecules were evaluated in vitro against breast cancer cell lines MDA-MB-231 and MCF-7, finding that compound 1f, at a concentration of 6.25 μM, exhibited a potential cytotoxic effect. Then, to understand the interactions between synthesized compounds and the main proteins related to the cancer cell lines, docking studies were performed on the serine/threonine kinase 1 (AKT1) and Orexetine type 2 receptor (Ox2R), finding moderate to strong binding energies, which matched accurately with the in vitro results. Additionally, molecular dynamics were performed between proteins related to the studied cell lines and the three best ligands.

ChemFlow_py: a flexible toolkit for docking and rescoring

The design of accurate virtual screening tools is an open challenge in drug discovery. Several structure-based methods have been developed at different levels of approximation. Among them, molecular docking is an established technique with high efficiency, but typically low accuracy. Moreover, docking performances are known to be target-dependent, which makes the choice of the docking program and corresponding scoring function critical when approaching a new protein target. To compare the performances of different docking protocols, we developed ChemFlow_py, an automated tool to perform docking and rescoring. Using four protein systems extracted from DUD-E with 100 known active compounds and 3000 decoys per target, we compared the performances of several rescoring strategies including consensus scoring. We found that the average docking results can be improved by consensus ranking, which emphasizes the relevance of consensus scoring when little or no chemical information is available for a given target. ChemFlow_py is a free toolkit to optimize the performances of virtual high-throughput screening (vHTS). The software is publicly available at https://github.com/IFMlab/ChemFlow_py .

Repurposing Navitoclax to block SARS-CoV-2 fusion and entry by targeting heptapeptide repeat sequence 1 in S2 protein

Along with the long pandemic of COVID-19 caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has come the dilemma of emerging viral variants of concern (VOC), particularly Omicron and its subvariants, able to deftly escape immune surveillance and the otherwise protective effect of current vaccines and antibody drugs. We previously identified a peptide-based pan-CoV fusion inhibitor, termed as EK1, able to bind the HR1 region in viral spike (S) protein S2 subunit. This effectively blocked formation of the six-helix bundle (6-HB) fusion core and, thus, showed efficacy against all human coronaviruses (HCoVs). EK1 is now in phase 3 clinical trials. However, the peptide drug generally lacks oral availability. Therefore, we herein performed a structure-based virtual screening of the libraries of biologically active molecules and identified nine candidate compounds. One is Navitoclax, an orally active anticancer drug by inhibition of Bcl-2. Like EK1 peptide, it could bind HR1 and block 6-HB formation, efficiently inhibiting fusion and infection of all SARS-CoV-2 variants tested, as well as SARS-CoV and MERS-CoV, with IC50 values ranging from 0.5 to 3.7 μM. These findings suggest that Navitoclax is a promising repurposed drug candidate for development as a safe and orally available broad-spectrum antiviral drug to combat the current SARS-CoV-2 and its variants, as well as other HCoVs.

Consolidated knowledge-guided computational pipeline for therapeutic intervention against bacterial biofilms - a review

Biofilm-associated bacterial infections attributed to multifactorial antimicrobial resistance have caused worldwide challenges in formulating successful treatment strategies. In search of accelerated yet cost-effective therapeutics, several researchers have opted for bioinformatics-based protocols to systemize targeted therapies against biofilm-producing strains. The present review investigated the up-to-date computational databases and servers dedicated to anti-biofilm research to design/screen novel biofilm inhibitors (antimicrobial peptides/phytocompounds/synthetic compounds) and predict their biofilm-inhibition efficacy. Scrutinizing the contemporary in silico methods, a consolidated approach has been highlighted, referred to as a knowledge-guided computational pipeline for biofilm-targeted therapy. The proposed pipeline has amalgamated prominently employed methodologies in genomics, transcriptomics, interactomics and proteomics to identify potential target proteins and their complementary anti-biofilm compounds for effective functional inhibition of biofilm-linked pathways. This review can pave the way for new portals to formulate successful therapeutic interventions against biofilm-producing pathogens.

Novel Inhibitors of androgen receptor's DNA binding domain identified using an ultra-large virtual screening

Androgen receptor (AR) inhibition remains the primary strategy to combat the progression of prostate cancer (PC). However, all clinically used AR inhibitors target the ligand-binding domain (LBD), which is highly susceptible to truncations through splicing or mutations that confer drug resistance. Thus, there exists an urgent need for AR inhibitors with novel modes of action. We thus launched a virtual screening of an ultra-large chemical library to find novel inhibitors of the AR DNA-binding domain (DBD) at two sites: protein-DNA interface (P-box) and dimerization site (D-box). The compounds selected through vigorous computational filtering were then experimentally validated. We identified several novel chemotypes that effectively suppress transcriptional activity of AR and its splice variant V7. The identified compounds represent previously unexplored chemical scaffolds with a mechanism of action that evades the conventional drug resistance manifested through LBD mutations. Additionally, we describe the binding features required to inhibit AR DBD at both P-box and D-box target sites.

AI-Driven De Novo Design and Molecular Modeling for Discovery of Small-Molecule Compounds as Potential Drug Candidates Targeting SARS-CoV-2 Main Protease

Over the past three years, significant progress has been made in the development of novel promising drug candidates against COVID-19. However, SARS-CoV-2 mutations resulting in the emergence of new viral strains that can be resistant to the drugs used currently in the clinic necessitate the development of novel potent and broad therapeutic agents targeting different vulnerable spots of the viral proteins. In this study, two deep learning generative models were developed and used in combination with molecular modeling tools for de novo design of small molecule compounds that can inhibit the catalytic activity of SARS-CoV-2 main protease (Mpro), an enzyme critically important for mediating viral replication and transcription. As a result, the seven best scoring compounds that exhibited low values of binding free energy comparable with those calculated for two potent inhibitors of Mpro, via the same computational protocol, were selected as the most probable inhibitors of the enzyme catalytic site. In light of the data obtained, the identified compounds are assumed to present promising scaffolds for the development of new potent and broad-spectrum drugs inhibiting SARS-CoV-2 Mpro, an attractive therapeutic target for anti-COVID-19 agents.

Hic Sunt Dracones: Molecular Docking in Uncharted Territories with Structures from AlphaFold2 and RoseTTAfold

AlphaFold2 and RoseTTAfold impress with their high accuracy in protein structure prediction. However, for structure-based virtual screenings, not only the overall structure but especially the binding sites need to be accurately predicted. In this work, the docking performance for 66 targets with known ligands but without experimental structures available in the protein data bank was elucidated. The results suggest that using an experimental surrogate-ligand complex is often superior over homology models, and only at low sequence identity to the closest homologue AlphaFold2 structures show an equal performance. The generally high fluctuation of receiver operating characteristic area under the curve values obtained for different homology models suggests that multiple combinations of docking programs and homology models should be tested prior to prospective virtual screenings, and in some cases post-processing of crude models might be necessary.

Advances in the discovery of new chemotypes through ultra-large library docking

INTRODUCTION

The size and complexity of virtual screening libraries in drug discovery have skyrocketed in recent years, reaching up to multiple billions of accessible compounds. However, virtual screening of such ultra-large libraries poses several challenges associated with preparing the libraries, sampling, and pre-selection of suitable compounds. The utilization of artificial intelligence (AI)-assisted screening approaches, such as deep learning, poses a promising countermeasure to deal with this rapidly expanding chemical space. For example, various AI-driven methods were recently successfully used to identify novel small molecule inhibitors of the SARS-CoV-2 main protease (M^pro).

AREAS COVERED

This review focuses on presenting various kinds of virtual screening methods suitable for dealing with ultra-large libraries. Challenges associated with these computational methodologies are discussed, and recent advances are highlighted in the example of the discovery of novel M^pro inhibitors targeting the SARS-CoV-2 virus.

EXPERT OPINION

With the rapid expansion of the virtual chemical space, the methodologies for docking and screening such quantities of molecules need to keep pace. Employment of AI-driven screening compounds has already been shown to be effective in a range from a few thousand to multiple billion compounds, furthered by de novo generation of drug-like molecules without human interference.

How good are AlphaFold models for docking-based virtual screening?

A crucial component in structure-based drug discovery is the availability of high-quality three-dimensional structures of the protein target. Whenever experimental structures were not available, homology modeling has been, so far, the method of choice. Recently, AlphaFold (AF), an artificial-intelligence-based protein structure prediction method, has shown impressive results in terms of model accuracy. This outstanding success prompted us to evaluate how accurate AF models are from the perspective of docking-based drug discovery. We compared the high-throughput docking (HTD) performance of AF models with their corresponding experimental PDB structures using a benchmark set of 22 targets. The AF models showed consistently worse performance using four docking programs and two consensus techniques. Although AlphaFold shows a remarkable ability to predict protein architecture, this might not be enough to guarantee that AF models can be reliably used for HTD, and post-modeling refinement strategies might be key to increase the chances of success.

Apo2ph4: A Versatile Workflow for the Generation of Receptor-based Pharmacophore Models for Virtual Screening

Pharmacophore models are widely used as efficient virtual screening (VS) filters for the target-directed enrichment of large compound libraries. However, the generation of pharmacophore models that have the power to discriminate between active and inactive molecules traditionally requires structural information about ligand-target complexes or at the very least knowledge of one active ligand. The fact that the discovery of the first known active ligand of a newly investigated target represents a major hurdle at the beginning of every drug discovery project underscores the need for methods that are able to derive high-quality pharmacophore models even without the prior knowledge of any active ligand structures. In this work, we introduce a novel workflow, called apo2ph4, that enables the rapid derivation of pharmacophore models solely from the three-dimensional structure of the target receptor. The utility of this workflow is demonstrated retrospectively for the generation of a pharmacophore model for the M2 muscarinic acetylcholine receptor. Furthermore, in order to show the general applicability of apo2ph4, the workflow was employed for all 15 targets of the recently published LIT-PCBA dataset. Pharmacophore-based VS runs using the apo2ph4-derived models achieved a significant enrichment of actives for 13 targets. In the last presented example, a pharmacophore model derived from the etomidate site of the α1β2γ2 GABA_A receptor was used in VS campaigns. Subsequent in vitro testing of selected hits revealed that 19 out of 20 (95%) tested compounds were able to significantly enhance GABA currents, which impressively demonstrates the applicability of apo2ph4 for real-world drug design projects.

Structure-Based Discovery of Small-Molecule Inhibitors of the Autocatalytic Proliferation of α-Synuclein Aggregates

The presence of amyloid fibrils of α-synuclein is closely associated with Parkinson's disease and related synucleinopathies. It is still very challenging, however, to systematically discover small molecules that prevent the formation of these aberrant aggregates. Here, we describe a structure-based approach to identify small molecules that specifically inhibit the surface-catalyzed secondary nucleation step in the aggregation of α-synuclein by binding to the surface of the amyloid fibrils. The resulting small molecules are screened using a range of kinetic and thermodynamic assays for their ability to bind α-synuclein fibrils and prevent the further generation of α-synuclein oligomers. This study demonstrates that the combination of structure-based and kinetic-based drug discovery methods can lead to the identification of small molecules that selectively inhibit the autocatalytic proliferation of α-synuclein aggregates.

Towards Effective Consensus Scoring in Structure-Based Virtual Screening

Virtual screening (VS) is a computational strategy that uses in silico automated protein docking inter alia to rank potential ligands, or by extension rank protein-ligand pairs, identifying potential drug candidates. Most docking methods use preferred sets of physicochemical descriptors (PCDs) to model the interactions between host and guest molecules. Thus, conventional VS is often data-specific, method-dependent and with demonstrably differing utility in identifying candidate drugs. This study proposes four universality classes of novel consensus scoring (CS) algorithms that combine docking scores, derived from ten docking programs (ADFR, DOCK, Gemdock, Ledock, PLANTS, PSOVina, QuickVina2, Smina, Autodock Vina and VinaXB), using decoys from the DUD-E repository ( http://dude.docking.org/ ) against 29 MRSA-oriented targets to create a general VS formulation that can identify active ligands for any suitable protein target. Our results demonstrate that CS provides improved ligand-protein docking fidelity when compared to individual docking platforms. This approach requires only a small number of docking combinations and can serve as a viable and parsimonious alternative to more computationally expensive docking approaches. Predictions from our CS algorithm are compared against independent machine learning evaluations using the same docking data, complementing the CS outcomes. Our method is a reliable approach for identifying protein targets and high-affinity ligands that can be tested as high-probability candidates for drug repositioning.

Fragmented blind docking: a novel protein-ligand binding prediction protocol

In the present paper we propose a novel blind docking protocol based on Autodock-Vina. The developed docking protocol can provide binding site identification and binding pose prediction at the same time, by a systematical exploration of the protein volume performed with several preliminary docking calculations. In our opinion, this protocol can be successfully applied during the first steps of the virtual screening pipeline, because it provides binding site identification and binding pose prediction at the same time without visual evaluation of the binding site. After the binding pose prediction, MM/GBSA re-scoring rescoring procedures has been applied to improve the accuracy of the protein-ligand bound state. The FRAD protocol has been tested on 116 protein-ligand complexes of the Heat Shock Protein 90 - alpha, on 176 of Human Immunodeficiency virus protease 1, and on more than 100 protein-ligand system taken from the PDBbind dataset. Overall, the FRAD approach combined to MM/GBSA re-scoring can be considered as a powerful tool to increase the accuracy and efficiency with respect to other standard docking approaches when the ligand-binding site is unknown.Communicated by Ramaswamy H. Sarma.

Comprehensive Survey of Consensus Docking for High-Throughput Virtual Screening

The rapid advances of 3D techniques for the structural determination of proteins and the development of numerous computational methods and strategies have led to identifying highly active compounds in computer drug design. Molecular docking is a method widely used in high-throughput virtual screening campaigns to filter potential ligands targeted to proteins. A great variety of docking programs are currently available, which differ in the algorithms and approaches used to predict the binding mode and the affinity of the ligand. All programs heavily rely on scoring functions to accurately predict ligand binding affinity, and despite differences in performance, none of these docking programs is preferable to the others. To overcome this problem, consensus scoring methods improve the outcome of virtual screening by averaging the rank or score of individual molecules obtained from different docking programs. The successful application of consensus docking in high-throughput virtual screening highlights the need to optimize the predictive power of molecular docking methods.

Molecular docking-based virtual screening: Challenges in hits identification for Anti-SARS-Cov-2 activity

The emergence of severe acute respiratory syndrome coronavirus 2 (SARS-Cov-2) requires finding new drugs or repurposing drugs for clinical use. Molecular docking belongs to structure-based drug design providing a fast method for identifying the hit compounds with antiviral activity against SARS-Cov-2. However, the weakness of the docking method is compounded by the limited crystallographic information and comparison drugs due to the novelty of this virus can present challenges in identifying hits of anti-SARS-Cov-2. In the current review, we highlighted several aspects, especially those related to the target structure, docking validation, and virtual hit selection, that need to be considered to obtain reliable docking results. Here, we discussed several cases pertaining to the issue highlighted and approaches that could be used to solve them.

Pharmacophore model-aided virtual screening combined with comparative molecular docking and molecular dynamics for identification of marine natural products as SARS-CoV-2 papain-like protease inhibitors

Targeting SARS-CoV-2 papain-like protease using inhibitors is a suitable approach for inhibition of virus replication and dysregulation of host anti-viral immunity. Engaging all five binding sites far from the catalytic site of PLpro is essential for developing a potent inhibitor. We developed and validated a structure-based pharmacophore model with 9 features of a potent PLpro inhibitor. The pharmacophore model-aided virtual screening of the comprehensive marine natural product database predicted 66 initial hits. This hit library was downsized by filtration through a molecular weight filter of ≤ 500 g/mol. The 50 resultant hits were screened by comparative molecular docking using AutoDock and AutoDock Vina. Comparative molecular docking enables benchmarking docking and relieves the disparities in the search and scoring functions of docking engines. Both docking engines retrieved 3 same compounds at different positions in the top 1 % rank, hence consensus scoring was applied, through which CMNPD28766, aspergillipeptide F emerged as the best PLpro inhibitor. Aspergillipeptide F topped the 50-hit library with a pharmacophore-fit score of 75.916. Favorable binding interactions were predicted between aspergillipeptide F and PLpro similar to the native ligand XR8-24. Aspergillipeptide F was able to engage all the 5 binding sites including the newly discovered BL2 groove, site V. Molecular dynamics for quantification of Cα-atom movements of PLpro after ligand binding indicated that it exhibits highly correlated domain movements contributing to the low free energy of binding and a stable conformation. Thus, aspergillipeptide F is a promising candidate for pharmaceutical and clinical development as a potent SARS-CoV-2 PLpro inhibitor.

MILCDock: Machine Learning Enhanced Consensus Docking for Virtual Screening in Drug Discovery

Molecular docking tools are regularly used to computationally identify new molecules in virtual screening for drug discovery. However, docking tools suffer from inaccurate scoring functions with widely varying performance on different proteins. To enable more accurate ranking of active over inactive ligands in virtual screening, we created a machine learning consensus docking tool, MILCDock, that uses predictions from five traditional molecular docking tools to predict the probability a ligand binds to a protein. MILCDock was trained and tested on data from both the DUD-E and LIT-PCBA docking datasets and shows improved performance over traditional molecular docking tools and other consensus docking methods on the DUD-E dataset. LIT-PCBA targets proved to be difficult for all methods tested. We also find that DUD-E data, although biased, can be effective in training machine learning tools if care is taken to avoid DUD-E's biases during training.

Boosting the Full Potential of PyMOL with Structural Biology Plugins

Over the past few decades, the number of available structural bioinformatics pipelines, libraries, plugins, web resources and software has increased exponentially and become accessible to the broad realm of life scientists. This expansion has shaped the field as a tangled network of methods, algorithms and user interfaces. In recent years PyMOL, widely used software for biomolecules visualization and analysis, has started to play a key role in providing an open platform for the successful implementation of expert knowledge into an easy-to-use molecular graphics tool. This review outlines the plugins and features that make PyMOL an eligible environment for supporting structural bioinformatics analyses.

Systematic Investigation of Docking Failures in Large-Scale Structure-Based Virtual Screening

In recent years, large-scale structure-based virtual screening has attracted increasing levels of interest for identification of novel compounds corresponding to potential drug targets. It is critical to understand the strengths and weaknesses of docking algorithms to increase the success rate in practical applications. Here, we systematically investigated the docking successes and failures of two representative docking programs: UCSF DOCK 3.7 and AutoDock Vina. DOCK 3.7 performed better in early enrichment on the Directory of Useful Decoys: Enhanced (DUD-E) data set, although both docking methods were roughly comparable in overall enrichment performance. DOCK 3.7 also showed superior computational efficiency. Intriguingly, the Vina scoring function showed a bias toward compounds with higher molecular weights. Both the tested docking approaches yielded incorrectly predicted ligand binding poses caused by the limitations of torsion sampling. Based on a careful analysis of docking results from six representative cases, we propose the reasons underlying docking failures; furthermore, we provide a few solutions, representing practical guidance for large-scale virtual screening campaigns and future docking algorithm development.

Mono-n-butyl Malate-Derived Compounds from Camu-camu (Myrciaria dubia) Malic Acid: The Alkyl-Dependent Antihyperglycemic-Related Activity

Malic acid derivatives from camu-camu (Myrciaria dubia) fruit exhibited a strong in vitro inhibitory activity toward pancreatic α-amylase and α-glucosidase enzymes. During a bioguided chromatographic fractionation process of the whole fruit (pulp and peelings) polar extract, isomers (S)-4-butoxy-2-hydroxy-4-oxobutanoic acid (1) and (S)-4-butoxy-3-hydroxy-4-oxobutanoic acid (2) (84:16) were isolated and identified as a potent inhibitor of α-amylase (IC₅₀= 11.69 ± 1.75 μg/mL) and α-glucosidase (IC₅₀ = 102.69 ± 4.16 μg/mL). The chemical structures were confirmed by HPLC-ESIMS and ¹H and ¹³C NMR (one- and two-dimensional) analyses. The structure-based virtual screening demonstrated that the aliphatic moiety plays a significant role in the binding mode of the test alkyl malate esters. Compound 1 exhibited the best interaction profile to bind both enzymes, having key structural features to form relevant contacts by involving adequate enzyme-ligand complex stabilization and compactness over time.

Ginsenoside Rg5 allosterically interacts with P2RY12 and ameliorates deep venous thrombosis by counteracting neutrophil NETosis and inflammatory response

Background

Deep venous thrombosis (DVT) highly occurs in patients with severe COVID-19 and probably accounted for their high mortality. DVT formation is a time-dependent inflammatory process in which NETosis plays an important role. However, whether ginsenoside Rg5 from species of Panax genus could alleviate DVT and its underlying mechanism has not been elucidated.

Methods

The interaction between Rg5 and P2RY₁₂ was studied by molecular docking, molecular dynamics, surface plasmon resonance (SPR), and molecular biology assays. The preventive effect of Rg5 on DVT was evaluated in inferior vena cava stasis–induced mice, and immunocytochemistry, Western blot, and calcium flux assay were performed in neutrophils from bone marrow to explore the mechanism of Rg5 in NETosis via P2RY₁₂.

Results

Rg5 allosterically interacted with P2RY₁₂, formed stable complex, and antagonized its activity via residue E188 and R265. Rg5 ameliorated the formation of thrombus in DVT mice; accompanied by decreased release of Interleukin (IL)-6, IL-1β, and tumor necrosis factor-α in plasma; and suppressed neutrophil infiltration and neutrophil extracellular trap (NET) release. In lipopolysaccharide- and platelet-activating factor–induced neutrophils, Rg5 reduced inflammatory responses via inhibiting the activation of ERK/NF-κB signaling pathway while decreasing cellular Ca²⁺ concentration, thus reducing the activity and expression of peptidyl arginine deiminase 4 to prevent NETosis. The inhibitory effect on neutrophil activity was dependent on P2RY₁₂.

Conclusions

Rg5 could attenuate experimental DVT by counteracting NETosis and inflammatory response in neutrophils via P2RY₁₂, which may pave the road for its clinical application in the prevention of DVT-related disorders.

3pHLA-score improves structure-based peptide-HLA binding affinity prediction

Binding of peptides to Human Leukocyte Antigen (HLA) receptors is a prerequisite for triggering immune response. Estimating peptide-HLA (pHLA) binding is crucial for peptide vaccine target identification and epitope discovery pipelines. Computational methods for binding affinity prediction can accelerate these pipelines. Currently, most of those computational methods rely exclusively on sequence-based data, which leads to inherent limitations. Recent studies have shown that structure-based data can address some of these limitations. In this work we propose a novel machine learning (ML) structure-based protocol to predict binding affinity of peptides to HLA receptors. For that, we engineer the input features for ML models by decoupling energy contributions at different residue positions in peptides, which leads to our novel per-peptide-position protocol. Using Rosetta's ref2015 scoring function as a baseline we use this protocol to develop 3pHLA-score. Our per-peptide-position protocol outperforms the standard training protocol and leads to an increase from 0.82 to 0.99 of the area under the precision-recall curve. 3pHLA-score outperforms widely used scoring functions (AutoDock4, Vina, Dope, Vinardo, FoldX, GradDock) in a structural virtual screening task. Overall, this work brings structure-based methods one step closer to epitope discovery pipelines and could help advance the development of cancer and viral vaccines.

Consensus scoring evaluated using the GPCR-Bench dataset: Reconsidering the role of MM/GBSA

The recent availability of large numbers of GPCR crystal structures has provided an unprecedented opportunity to evaluate their performance in virtual screening protocols using established benchmarking datasets. In this study, we evaluated the ability of MM/GBSA in consensus scoring-based virtual screening enrichment together with nine classical scoring functions, using the GPCR-Bench dataset consisting of 24 GPCR crystal structures and 254,646 actives and decoys. While the performance of consensus scoring was modest overall, combinations which included MM/GBSA performed relatively well compared to combinations of classical scoring functions. Combinations of MM/GBSA and good-performing scoring functions provided the highest proportion of improvements, with improvements observed in 32% and 19% of all combinations across all targets at the EF_1% and EF_5% levels respectively. Combinations of MM/GBSA and poor-performing scoring functions still outperformed classical scoring functions, with improvements observed in 26% and 17% of all combinations at the EF_1% and EF_5% levels. In comparison, only 14-22% and 6-11% of combinations of classical scoring functions produced improvements at EF_1% and EF_5% respectively. Efforts to improve performance by increasing the number of scoring functions in consensus scoring to three were mostly ineffective. We also observed that consensus scoring performed better for individual scoring functions possessing initially low enrichment factors, potentially implying their benefits are more relevant in such scenarios. Overall, this study demonstrated the first implementation of MM/GBSA in consensus scoring using the GPCR-Bench dataset and could provide a valuable benchmark of the performance of MM/GBSA in comparison to classical scoring functions in consensus scoring for GPCRs.

In Silico Drug Repositioning to Target the SARS-CoV-2 Main Protease as Covalent Inhibitors Employing a Combined Structure-Based Virtual Screening Strategy of Pharmacophore Models and Covalent Docking

The epidemic caused by the SARS-CoV-2 coronavirus, which has spread rapidly throughout the world, requires urgent and effective treatments considering that the appearance of viral variants limits the efficacy of vaccines. The main protease of SARS-CoV-2 (M^pro) is a highly conserved cysteine proteinase, fundamental for the replication of the coronavirus and with a specific cleavage mechanism that positions it as an attractive therapeutic target for the proposal of irreversible inhibitors. A structure-based strategy combining 3D pharmacophoric modeling, virtual screening, and covalent docking was employed to identify the interactions required for molecular recognition, as well as the spatial orientation of the electrophilic warhead, of various drugs, to achieve a covalent interaction with Cys145 of M^pro. The virtual screening on the structure-based pharmacophoric map of the SARS-CoV-2 M^pro in complex with an inhibitor N3 (reference compound) provided high efficiency by identifying 53 drugs (FDA and DrugBank databases) with probabilities of covalent binding, including N3 (Michael acceptor) and others with a variety of electrophilic warheads. Adding the energy contributions of affinity for non-covalent and covalent docking, 16 promising drugs were obtained. Our findings suggest that the FDA-approved drugs Vaborbactam, Cimetidine, Ixazomib, Scopolamine, and Bicalutamide, as well as the other investigational peptide-like drugs (DB04234, DB03456, DB07224, DB7252, and CMX-2043) are potential covalent inhibitors of SARS-CoV-2 M^pro.

A human-based multi-gene signature enables quantitative drug repurposing for metabolic disease

Insulin resistance (IR) contributes to the pathophysiology of diabetes, dementia, viral infection, and cardiovascular disease. Drug repurposing (DR) may identify treatments for IR; however, barriers include uncertainty whether in vitro transcriptomic assays yield quantitative pharmacological data, or how to optimise assay design to best reflect in vivo human disease. We developed a clinical-based human tissue IR signature by combining lifestyle-mediated treatment responses (>500 human adipose and muscle biopsies) with biomarkers of disease status (fasting IR from >1200 biopsies). The assay identified a chemically diverse set of >130 positively acting compounds, highly enriched in true positives, that targeted 73 proteins regulating IR pathways. Our multi-gene RNA assay score reflected the quantitative pharmacological properties of a set of epidermal growth factor receptor-related tyrosine kinase inhibitors, providing insight into drug target specificity; an observation supported by deep learning-based genome-wide predicted pharmacology. Several drugs identified are suitable for evaluation in patients, particularly those with either acute or severe chronic IR.

Consensus combining outcomes of multiple ensemble dockings: examples using dDAT crystalized complexes

Docking using different programs provides more reliable information about the interaction of molecules than data obtained in a single program. An exponential consensus ranking (ECR) was developed to combine scoring functions across docking programs differing in efficiencies and scales of measurements. The ECR method was adapted to merge results of re- and cross-dockings (i.e., ensemble docking) made in multiple docking programs. Adapted ECR consisted of four consecutive steps: 1- determination of scoring functions for a ligand with a series of macromolecules in multiple docking programs; 2- ranking of the scoring functions per macromolecule in each program; 3- combining the ranking across the programs creating a ranking per macromolecule; 4- averaging the ranking per macromolecule creating a final ranking. This last step incorporated the heterogeneity of the macromolecule conformations in the consensual score. The final ranking based on the adapted ECR represents relative affinity of a series of ligands to a macromolecule on average. As an example, a ranking of the average affinity of antidepressants and other ligands to the Drosophila melanogaster dopamine transporter (dDAT) was presented. Adapted ECR generated a ranking similar to that based on the affinity constant of each ligand obtained from the literature.

• A final ranking of the average relative affinity of different ligands to the dDAT.

• A consensus method combining multiple ensemble dockings.

• A complete protocol to make re-docking and cross-docking using Autodock Vina, Gold and DockThor.

Computational Study on Potential Novel Anti-Ebola Virus Protein VP35 Natural Compounds

Ebola virus (EBOV) is one of the most lethal pathogens that can infect humans. The Ebola viral protein VP35 (EBOV VP35) inhibits host IFN-α/β production by interfering with host immune responses to viral invasion and is thus considered as a plausible drug target. The aim of this study was to identify potential novel lead compounds against EBOV VP35 using computational techniques in drug discovery. The 3D structure of the EBOV VP35 with PDB ID: 3FKE was used for molecular docking studies. An integrated library of 7675 African natural product was pre-filtered using ADMET risk, with a threshold of 7 and, as a result, 1470 ligands were obtained for the downstream molecular docking using AutoDock Vina, after an energy minimization of the protein via GROMACS. Five known inhibitors, namely, amodiaquine, chloroquine, gossypetin, taxifolin and EGCG were used as standard control compounds for this study. The area under the curve (AUC) value, evaluating the docking protocol obtained from the receiver operating characteristic (ROC) curve, generated was 0.72, which was considered to be acceptable. The four identified potential lead compounds of NANPDB4048, NANPDB2412, ZINC000095486250 and NANPDB2476 had binding affinities of −8.2, −8.2, −8.1 and −8.0 kcal/mol, respectively, and were predicted to possess desirable antiviral activity including the inhibition of RNA synthesis and membrane permeability, with the probable activity (Pa) being greater than the probable inactivity (Pi) values. The predicted anti-EBOV inhibition efficiency values (IC₅₀), found using a random forest classifier, ranged from 3.35 to 11.99 μM, while the Ki values ranged from 0.97 to 1.37 μM. The compounds NANPDB4048 and NANPDB2412 had the lowest binding energy of −8.2 kcal/mol, implying a higher binding affinity to EBOV VP35 which was greater than those of the known inhibitors. The compounds were predicted to possess a low toxicity risk and to possess reasonably good pharmacological profiles. Molecular dynamics (MD) simulations of the protein–ligand complexes, lasting 50 ns, and molecular mechanisms Poisson-Boltzmann surface area (MM-PBSA) calculations corroborated the binding affinities of the identified compounds and identified novel critical interacting residues. The antiviral potential of the molecules could be confirmed experimentally, while the scaffolds could be optimized for the design of future novel anti-EBOV chemotherapeutics.

Target identification for repurposed drugs active against SARS-CoV-2 via high-throughput inverse docking

Screening already approved drugs for activity against a novel pathogen can be an important part of global rapid-response strategies in pandemics. Such high-throughput repurposing screens have already identified several existing drugs with potential to combat SARS-CoV-2. However, moving these hits forward for possible development into drugs specifically against this pathogen requires unambiguous identification of their corresponding targets, something the high-throughput screens are not typically designed to reveal. We present here a new computational inverse-docking protocol that uses all-atom protein structures and a combination of docking methods to rank-order targets for each of several existing drugs for which a plurality of recent high-throughput screens detected anti-SARS-CoV-2 activity. We demonstrate validation of this method with known drug-target pairs, including both non-antiviral and antiviral compounds. We subjected 152 distinct drugs potentially suitable for repurposing to the inverse docking procedure. The most common preferential targets were the human enzymes TMPRSS2 and PIKfyve, followed by the viral enzymes Helicase and PLpro. All compounds that selected TMPRSS2 are known serine protease inhibitors, and those that selected PIKfyve are known tyrosine kinase inhibitors. Detailed structural analysis of the docking poses revealed important insights into why these selections arose, and could potentially lead to more rational design of new drugs against these targets.

An integrated method for optimized identification of effective natural inhibitors against SARS-CoV-2 3CLpro

The current severe situation of coronavirus disease 2019 (COVID-19) caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has not been reversed and posed great threats to global health. Therefore, there is an urgent need to find out effective antiviral drugs. The 3-chymotrypsin-like protease (3CLpro) in SARS-CoV-2 serve as a promising anti-virus target due to its essential role in the regulation of virus reproduction. Here, we report an improved integrated approach to identify effective 3CLpro inhibitors from effective Chinese herbal formulas. With this approach, we identified the 5 natural products (NPs) including narcissoside, kaempferol-3-O-gentiobioside, rutin, vicenin-2 and isoschaftoside as potential anti-SARS-CoV-2 candidates. Subsequent molecular dynamics simulation additionally revealed that these molecules can be tightly bound to 3CLpro and confirmed effectiveness against COVID-19. Moreover, kaempferol-3-o-gentiobioside, vicenin-2 and isoschaftoside were first reported to have SARS-CoV-2 3CLpro inhibitory activity. In summary, this optimized integrated strategy for drug screening can be utilized in the discovery of antiviral drugs to achieve rapid acquisition of drugs with specific effects on antiviral targets.

Improving Structure-Based Virtual Screening with Ensemble Docking and Machine Learning

One of the main challenges of structure-based virtual screening (SBVS) is the incorporation of the receptor's flexibility, as its explicit representation in every docking run implies a high computational cost. Therefore, a common alternative to include the receptor's flexibility is the approach known as ensemble docking. Ensemble docking consists of using a set of receptor conformations and performing the docking assays over each of them. However, there is still no agreement on how to combine the ensemble docking results to obtain the final ligand ranking. A common choice is to use consensus strategies to aggregate the ensemble docking scores, but these strategies exhibit slight improvement regarding the single-structure approach. Here, we claim that using machine learning (ML) methodologies over the ensemble docking results could improve the predictive power of SBVS. To test this hypothesis, four proteins were selected as study cases: CDK2, FXa, EGFR, and HSP90. Protein conformational ensembles were built from crystallographic structures, whereas the evaluated compound library comprised up to three benchmarking data sets (DUD, DEKOIS 2.0, and CSAR-2012) and cocrystallized molecules. Ensemble docking results were processed through 30 repetitions of 4-fold cross-validation to train and validate two ML classifiers: logistic regression and gradient boosting trees. Our results indicate that the ML classifiers significantly outperform traditional consensus strategies and even the best performance case achieved with single-structure docking. We provide statistical evidence that supports the effectiveness of ML to improve the ensemble docking performance.

Combination of pose and rank consensus in docking-based virtual screening: the best of both worlds

The use of high-throughput docking (HTD) in the drug discovery pipeline is today widely established. In spite of methodological improvements in docking accuracy (pose prediction), scoring power, ranking power, and screening power in HTD remain challenging. In fact, pose prediction is of critical importance in view of the pose-dependent scoring process, since incorrect poses will necessarily decrease the ranking power of scoring functions. The combination of results from different docking programs (consensus scoring) has been shown to improve the performance of HTD. Moreover, it has been also shown that a pose consensus approach might also result in database enrichment. We present a new methodology named Pose/Ranking Consensus (PRC) that combines both pose and ranking consensus approaches, to overcome the limitations of each stand-alone strategy. This approach has been developed using four docking programs (ICM, rDock, Auto Dock 4, and PLANTS; the first one is commercial, the other three are free). We undertook a thorough analysis for the best way of combining pose and rank strategies, and applied the PRC to a wide range of 34 targets sampling different protein families and binding site properties. Our approach exhibits an improved systematic performance in terms of enrichment factor and hit rate with respect to either pose consensus or consensus ranking alone strategies at a lower computational cost, while always ensuring the recovery of a suitable number of ligands. An analysis using four free docking programs (replacing ICM by Auto Dock Vina) displayed comparable results.

The new methodology named Pose/Ranking Consensus (PRC) combines both pose and ranking consensus strategies. It displays an enhanced performance in terms of enrichment factor and hit rate, ensuring the recovery of a suitable number of ligands.

Synthesis and insecticidal activity of novel 1,2,4-triazole derivatives containing trifluoroacetyl moieties

A series of compounds containing trifluoroacetyl groups were synthesized, and their insecticidal activity against Nilaparvata lugens and Aphis craccivora was evaluated. The compound structure was identified by NMR, HRMS, and single-crystal diffraction. The bioassay results indicated that compound 4-1 (R¹ is chloropyridine, R² is H), 4-2 (R¹ is chlorothiazole, R² is H) and 4-19 (R¹ is benzyl, R² is isopropyl) had the best activity against Nilaparvata lugens (93.5%, 94.1% and 95.5%) at 100 mg/L concentration. The effect of different substituents of R¹ or R² on the activity was studied through the structure-activity relationship. Molecular docking of compounds 4-1 and 4-2 was discussed.

dockECR: Open consensus docking and ranking protocol for virtual screening of small molecules

The development of open computational pipelines to accelerate the discovery of treatments for emerging diseases allows finding novel solutions in shorter periods of time. Consensus molecular docking is one of these approaches, and its main purpose is to increase the detection of real actives within virtual screening campaigns. Here we present dockECR, an open consensus docking and ranking protocol that implements the exponential consensus ranking method to prioritize molecular candidates. The protocol uses four open source molecular docking programs: AutoDock Vina, Smina, LeDock and rDock, to rank the molecules. In addition, we introduce a scoring strategy based on the average RMSD obtained from comparing the best poses from each single program to complement the consensus ranking with information about the predicted poses. The protocol was benchmarked using 15 relevant protein targets with known actives and decoys, and applied using the main protease of the SARS-CoV-2 virus. For the application, different crystal structures of the protease, and frames obtained from molecular dynamics simulations were used to dock a library of 79 molecules derived from previously co-crystallized fragments. The ranking obtained with dockECR was used to prioritize eight candidates, which were evaluated in terms of the interactions generated with key residues from the protease. The protocol can be implemented in any virtual screening campaign involving proteins as molecular targets. The dockECR code is publicly available at: https://github.com/rochoa85/dockECR.

Strengths and Weaknesses of Docking Simulations in the SARS-CoV-2 Era: the Main Protease (Mpro) Case Study

The scientific community is working against the clock to arrive at therapeutic interventions to treat patients with COVID-19. Among the strategies for drug discovery, virtual screening approaches have the capacity to search potential hits within millions of chemical structures in days, with the appropriate computing infrastructure. In this article, we first analyzed the published research targeting the inhibition of the main protease (Mpro), one of the most studied targets of SARS-CoV-2, by docking-based methods. An alarming finding was the lack of an adequate validation of the docking protocols (i.e., pose prediction and virtual screening accuracy) before applying them in virtual screening campaigns. The performance of the docking protocols was tested at some level in 57.7% of the 168 investigations analyzed. However, we found only three examples of a complete retrospective analysis of the scoring functions to quantify the virtual screening accuracy of the methods. Moreover, only two publications reported some experimental evaluation of the proposed hits until preparing this manuscript. All of these findings led us to carry out a retrospective performance validation of three different docking protocols, through the analysis of their pose prediction and screening accuracy. Surprisingly, we found that even though all tested docking protocols have a good pose prediction, their screening accuracy is quite limited as they fail to correctly rank a test set of compounds. These results highlight the importance of conducting an adequate validation of the docking protocols before carrying out virtual screening campaigns, and to experimentally confirm the predictions made by the models before drawing bold conclusions. Finally, successful structure-based drug discovery investigations published during the redaction of this manuscript allow us to propose the inclusion of target flexibility and consensus scoring as alternatives to improve the accuracy of the methods.

Computational Evaluation of Interaction Between Curcumin Derivatives and Amyloid-β Monomers and Fibrils: Relevance to Alzheimer's Disease

BACKGROUND

The most important hallmark in the neuropathology of Alzheimer's disease (AD) is the formation of amyloid-β (Aβ) fibrils due to the misfolding/aggregation of the Aβ peptide. Preventing or reverting the aggregation process has been an active area of research. Naturally occurring products are a potential source of molecules that may be able to inhibit Aβ42 peptide aggregation. Recently, we and others reported the anti-aggregating properties of curcumin and some of its derivatives in vitro, presenting an important therapeutic avenue by enhancing these properties.

OBJECTIVE

To computationally assess the interaction between Aβ peptide and a set of curcumin derivatives previously explored in experimental assays.

METHODS

The interactions of ten ligands with Aβ monomers were studied by combining molecular dynamics and molecular docking simulations. We present the in silico evaluation of the interaction between these derivatives and the Aβ42 peptide, both in the monomeric and fibril forms.

RESULTS

The results show that a single substitution in curcumin could significantly enhance the interaction between the derivatives and the Aβ42 monomers when compared to a double substitution. In addition, the molecular docking simulations showed that the interaction between the curcumin derivatives and the Aβ42 monomers occur in a region critical for peptide aggregation.

CONCLUSION

Results showed that a single substitution in curcumin improved the interaction of the ligands with the Aβ monomer more so than a double substitution. Our molecular docking studies thus provide important insights for further developing/validating novel curcumin-derived molecules with high therapeutic potential for AD.

Use of molecular docking computational tools in drug discovery

Molecular docking has become an important component of the drug discovery process. Since first being developed in the 1980s, advancements in the power of computer hardware and the increasing number of and ease of access to small molecule and protein structures have contributed to the development of improved methods, making docking more popular in both industrial and academic settings. Over the years, the modalities by which docking is used to assist the different tasks of drug discovery have changed. Although initially developed and used as a standalone method, docking is now mostly employed in combination with other computational approaches within integrated workflows. Despite its invaluable contribution to the drug discovery process, molecular docking is still far from perfect. In this chapter we will provide an introduction to molecular docking and to the different docking procedures with a focus on several considerations and protocols, including protonation states, active site waters and consensus, that can greatly improve the docking results.

Discovery of a Potent and Selective Chikungunya Virus Envelope Protein Inhibitor through Computer-Aided Drug Design

The worldwide expansion of chikungunya virus (CHIKV) into tropical and subtropical areas in the last 15 years has posed a currently unmet need for vaccines and therapeutics. The E2-E1 envelope glycoprotein complex binds receptors on the host cell and promotes membrane fusion during CHIKV entry, thus constituting an attractive target for the development of antiviral drugs. In order to identify CHIKV antivirals acting through inhibition of the envelope glycoprotein complex function, our first approach was to search for amenable druggable sites within the E2-E1 heterodimer. We identified a pocket located in the interface between E2 and E1 around the fusion loop. Then, via a structure-based virtual screening approach and in vitro assay of antiviral activity, we identified compound 7 as a specific inhibitor of CHIKV. Through a lead optimization process, we obtained compound 11 that demonstrated increased antiviral activity and low cytotoxicity (EC₅₀ 1.6 μM, CC₅₀ 56.0 μM). Molecular dynamics simulations were carried out and described a possible interaction pattern of compound 11 and the E1-E2 dimer that could be useful for further optimization. As expected from target site selection, compound 11 inhibited virus internalization during CHIKV entry. In addition, virus populations resistant to compound 11 included mutation E2-P173S, which mapped to the proposed binding pocket, and second site mutation E1-Y24H. Construction of recombinant viruses showed that these mutations conferred antiviral resistance in the parental background. Finally, compound 11 presents acceptable solubility values and is chemically and enzymatically stable in different media. Altogether, these findings uncover a suitable pocket for the design of CHIKV entry inhibitors with promising antiviral activity and pharmacological profiles.

Combination of consensus and ensemble docking strategies for the discovery of human dihydroorotate dehydrogenase inhibitors

The inconsistencies in the performance of the virtual screening (VS) process, depending on the used software and structural conformation of the protein, is a challenging issue in the drug design and discovery field. Varying performance, especially in terms of early recognition of the potential hit compounds, negatively affects the whole process and leads to unnecessary waste of the time and resources. Appropriate application of the ensemble docking and consensus-scoring approaches can significantly increase reliability of the VS results. Dihydroorotate dehydrogenase (DHODH) is a key enzyme in the pyrimidine biosynthesis pathway. It is considered as a valuable therapeutic target in cancer, autoimmune and viral diseases. Based on the conducted benchmark study and analysis of the effect of different combinations of the applied methods and approaches, here we suggested a structure-based virtual screening (SBVS) workflow that can be used to increase the reliability of VS.