A Comprehensive Survey of Prospective Structure-Based Virtual Screening for Early Drug Discovery in the Past Fifteen Years

DiffBindFR: an SE(3) equivariant network for flexible protein-ligand docking

Molecular docking, a key technique in structure-based drug design, plays pivotal roles in protein-ligand interaction modeling, hit identification and optimization, in which accurate prediction of protein-ligand binding mode is essential. Conventional docking approaches perform well in redocking tasks with known protein binding pocket conformation in the complex state. However, in real-world docking scenario without knowing the protein binding conformation for a new ligand, accurately modeling the binding complex structure remains challenging as flexible docking is computationally expensive and inaccurate. Typical deep learning-based docking methods do not explicitly consider protein side chain conformations and fail to ensure the physical plausibility and detailed atomic interactions. In this study, we present DiffBindFR, a full-atom diffusion-based flexible docking model that operates over the product space of ligand overall movements and flexibility and pocket side chain torsion changes. We show that DiffBindFR has higher accuracy in producing native-like binding structures with physically plausible and detailed interactions than available docking methods. Furthermore, in the Apo and AlphaFold2 modeled structures, DiffBindFR demonstrates superior advantages in accurate ligand binding pose and protein binding conformation prediction, making it suitable for Apo and AlphaFold2 structure-based drug design. DiffBindFR provides a powerful flexible docking tool for modeling accurate protein-ligand binding structures.

Molecular docking based comparative study of antiviral compounds on SARS-CoV-2 spike protein

The urgent need for effective therapeutic interventions against SARS-CoV-2 has prompted extensive exploration of potential drug candidates. Among the viral proteins, the spike (S) protein presents an attractive target due to its critical role in viral entry and infection. In this study, we employed molecular docking techniques to investigate the binding affinities and interaction profiles of a panel of active compounds against the SARS-CoV-2 spike protein. Utilising computational simulations, we assessed the binding properties of these compounds within the receptor-binding domain (RBD) and other key regions of the spike protein. Our comparative analysis elucidates the differential binding patterns and identifies promising lead compounds with high binding affinity and favourable interaction profiles. Furthermore, we discuss the implications of these findings for the development of potential therapeutics targeting the SARS-CoV-2 spike protein. Using molecular docking and the Lipinski five rule, this study illustrates possible compounds with strong binding affinities, their molecular interactions, for both naturally occurring and man-made drugs. Computational approach is applied, and it is concluded that, drugs like Withanolide, Dihydroergocristine, Fenebrutinib, and Ergotamine shows binding energies between -8.3 and -9.1 kcal/mol, and are possible candidate for anti covid drug.

Interaction of Norsecurinine-Type Oligomeric Alkaloids with α-Tubulin: A Molecular Docking Study

The medicinal plant Securinega virosa (Roxb ex. Willd) Baill., also known as Flueggea virosa (Roxb. ex Willd.) Royle, is commonly used in traditional medicine in Africa and Asia for the management of diverse pathologies, such as parasite infections, diabetes, and gastrointestinal diseases. Numerous alkaloids have been isolated from the twigs and leaves of the plant, notably a variety of oligomeric indolizidine alkaloids derived from the monomers securinine and norsecurinine which both display anticancer properties. The recent discovery that securinine can bind to tubulin and inhibit microtubule assembly prompted us to investigate the potential binding of two series of alkaloids, fluevirosines A-H and fluevirosinine A-J, with the tubulin dimer by means of molecular modeling. These natural products are rare high-order alkaloids with tri-, tetra-, and pentameric norsecurinine motifs. Despite their large size (up to 2500 Å3), these alkaloids can bind easily to the large drug-binding cavity (about 4800 Å3) on α-tubulin facing the β-tubulin unit. The molecular docking analysis suggests that these hydrophobic macro-alkaloids can form stable complexes with α/β-tubulin. The tubulin-binding capacity varies depending on the alkaloid size and structure. Structure-binding relationships are discussed. The docking analysis identifies the trimer fluevirosine D, tetramer fluevirosinine D, and pentamer fluevirosinine H as the most interesting tubulin ligands in the series. This study is the first to propose a molecular target for these atypical oligomeric Securinega alkaloids.

A High-Quality Data Set of Protein-Ligand Binding Interactions Via Comparative Complex Structure Modeling

High-quality protein-ligand complex structures provide the basis for understanding the nature of noncovalent binding interactions at the atomic level and enable structure-based drug design. However, experimentally determined complex structures are scarce compared with the vast chemical space. In this study, we addressed this issue by constructing the BindingNet data set via comparative complex structure modeling, which contains 69,816 modeled high-quality protein-ligand complex structures with experimental binding affinity data. BindingNet provides valuable insights into investigating protein-ligand interactions, allowing visual inspection and interpretation of structural analogues' structure-activity relationships. It can also be used for evaluating machine-learning-based scoring functions. Our results indicate that machine learning models trained on BindingNet could reduce the bias caused by buried solvent-accessible surface area, as we previously found for models trained on the PDBbind data set. We also discussed strategies to improve BindingNet and its potential utilization for benchmarking the molecular docking methods and ligand binding free energy calculation approaches. The BindingNet complements PDBbind in constructing a sufficient and unbiased protein-ligand binding data set and is freely available at http://bindingnet.huanglab.org.cn.

Computational Methods in the Design of Anticancer Drugs

In recent years, continuous progress has been made in the development of new anticancer drugs, and several compounds (small molecules, engineered antibodies, immunomodulators, etc [...].

DockOpt: A Tool for Automatic Optimization of Docking Models

Molecular docking is a widely used technique for leveraging protein structure for ligand discovery, but it remains difficult to utilize due to limitations that have not been adequately addressed. Despite some progress toward automation, docking still requires expert guidance, hindering its adoption by a broader range of investigators. To make docking more accessible, we developed a new utility called DockOpt, which automates the creation, evaluation, and optimization of docking models prior to their deployment in large-scale prospective screens. DockOpt outperforms our previous automated pipeline across all 43 targets in the DUDE-Z benchmark data set, and the generated models for 84% of targets demonstrate sufficient enrichment to warrant their use in prospective screens, with normalized LogAUC values of at least 15%. DockOpt is available as part of the Python package Pydock3 included in the UCSF DOCK 3.8 distribution, which is available for free to academic researchers at https://dock.compbio.ucsf.edu and free for everyone upon registration at https://tldr.docking.org.

Structure-Based Identification of Organoruthenium Compounds as Nanomolar Antagonists of Cannabinoid Receptors

Metal complexes exhibit a diverse range of coordination geometries, representing novel privileged scaffolds with convenient click types of preparation inaccessible for typical carbon-centered organic compounds. Herein, we explored the opportunity to identify biologically active organometallic complexes by reverse docking of a rigid, minimum-size octahedral organoruthenium scaffold against thousands of protein-binding pockets. Interestingly, cannabinoid receptor type 1 (CB1) was identified based on the docking scores and the degree of overlap between the docked organoruthenium scaffold and the hydrophobic scaffold of the cocrystallized ligand. Further structure-based optimization led to the discovery of organoruthenium complexes with nanomolar binding affinities and high selectivity toward CB2. Our work indicates that octahedral organoruthenium scaffolds may be advantageous for targeting the large and hydrophobic binding pockets and that the reverse docking approach may facilitate the discovery of novel privileged scaffolds, such as organometallic complexes, for exploring chemical space in lead discovery.

Targeting Allosteric Site of PCSK9 Enzyme for the Identification of Small Molecule Inhibitors: An In Silico Drug Repurposing Study

The primary cause of atherosclerotic cardiovascular disease (ASCVD) is elevated levels of low-density lipoprotein cholesterol (LDL-C). Proprotein convertase subtilisin/kexin type 9 (PCSK9) plays a crucial role in this process by binding to the LDL receptor (LDL-R) domain, leading to reduced influx of LDL-C and decreased LDL-R cell surface presentation on hepatocytes, resulting higher circulating levels of LDL-C. As a consequence, PCSK9 has been identified as a crucial target for drug development against dyslipidemia and hypercholesterolemia, aiming to lower plasma LDL-C levels. This research endeavors to identify promising inhibitory candidates that target the allosteric site of PCSK9 through an in silico approach. To start with, the FDA-approved Drug Library from Selleckchem was selected and virtually screened by docking studies using Glide extra-precision (XP) docking mode and Smina software (Version 1.1.2). Subsequently, rescoring of 100 drug compounds showing good average docking scores were performed using Gnina software (Version 1.0) to generate CNN Score and CNN binding affinity. Among the drug compounds, amikacin, bestatin, and natamycin were found to exhibit higher docking scores and CNN affinities against the PCSK9 enzyme. Molecular dynamics simulations further confirmed that these drug molecules established the stable protein-ligand complexes when compared to the apo structure of PCSK9 and the complex with the co-crystallized ligand structure. Moreover, the MM-GBSA calculations revealed binding free energy values ranging from -84.22 to -76.39 kcal/mol, which were found comparable to those obtained for the co-crystallized ligand structure. In conclusion, these identified drug molecules have the potential to serve as inhibitors PCSK9 enzyme and these finding could pave the way for the development of new PCSK9 inhibitory drugs in future in vitro research.

DOCK 6: Incorporating hierarchical traversal through precomputed ligand conformations to enable large-scale docking

To allow DOCK 6 access to unprecedented chemical space for screening billions of small molecules, we have implemented features from DOCK 3.7 into DOCK 6, including a search routine that traverses precomputed ligand conformations stored in a hierarchical database. We tested them on the DUDE-Z and SB2012 test sets. The hierarchical database search routine is 16 times faster than anchor-and-grow. However, the ability of hierarchical database search to reproduce the experimental pose is 16% worse than that of anchor-and-grow. The enrichment performance is on average similar, but DOCK 3.7 has better enrichment than DOCK 6, and DOCK 6 is on average 1.7 times slower. However, with post-docking torsion minimization, DOCK 6 surpasses DOCK 3.7. A large-scale virtual screen is performed with DOCK 6 on 23 million fragment molecules. We use current features in DOCK 6 to complement hierarchical database calculations, including torsion minimization, which is not available in DOCK 3.7.

A flexible data-free framework for structure-based de novo drug design with reinforcement learning

Contemporary structure-based molecular generative methods have demonstrated their potential to model the geometric and energetic complementarity between ligands and receptors, thereby facilitating the design of molecules with favorable binding affinity and target specificity. Despite the introduction of deep generative models for molecular generation, the atom-wise generation paradigm that partially contradicts chemical intuition limits the validity and synthetic accessibility of the generated molecules. Additionally, the dependence of deep learning models on large-scale structural data has hindered their adaptability across different targets. To overcome these challenges, we present a novel search-based framework, 3D-MCTS, for structure-based de novo drug design. Distinct from prevailing atom-centric methods, 3D-MCTS employs a fragment-based molecular editing strategy. The fragments decomposed from small-molecule drugs are recombined under predefined retrosynthetic rules, offering improved drug-likeness and synthesizability, overcoming the inherent limitations of atom-based approaches. Leveraging multi-threaded parallel simulations combined with a real-time energy constraint-based pruning strategy, 3D-MCTS achieves remarkable efficiency. At a fixed computational cost, it outperforms other state-of-the-art (SOTA) methods by producing molecules with enhanced binding affinity. Furthermore, its fragment-based approach ensures the generation of more dependable binding conformations, exhibiting a success rate 43.6% higher than that of other SOTAs. This advantage becomes even more pronounced when handling targets that significantly deviate from the training dataset. 3D-MCTS is capable of achieving thirty times more hits with high binding affinity than traditional virtual screening methods, which demonstrates the superior ability of 3D-MCTS to explore chemical space. Moreover, the flexibility of our framework makes it easy to incorporate domain knowledge during the process, thereby enabling the generation of molecules with desirable pharmacophores and enhanced binding affinity. The adaptability of 3D-MCTS is further showcased in metalloprotein applications, highlighting its potential across various drug design scenarios.

Computational and biophysical methods for the discovery and optimization of covalent drugs

Drugs that act by covalently attaching to their targets have been used to treat human diseases for over a hundred years. However, the deliberate design of covalent drugs was discouraged due to concerns of toxicity and off-target effects. Recent successes in covalent drug discovery have sparked fresh interest in this field. New screening and testing methods aimed at covalent inhibitors can play pivotal roles in facilitating the discovery process. This feature article focuses on computational and biophysical advances originating from our labs over the past decade and how these approaches have contributed to the design of prolyl oligopeptidase (POP) and SARS-CoV-2 3CLpro covalent inhibitors.

Discovery of Novel SIRT1/2 Inhibitors with Effective Cytotoxicity against Human Leukemia Cells

The sirtuin enzyme family members, SIRT1 and SIRT2, play both tumor-promoting and tumor-suppressing roles, depending on the context and experimental conditions. Compounds that inhibit either SIRT1 or SIRT2 show promising antitumor effects in several types of cancer models, both in vitro and in vivo. The simultaneous inhibition of SIRT1 and SIRT2 is helpful in treating cancer by completely blocking p53 deacetylation, leading to cell death. However, only a few SIRT1/2 dual inhibitors have been developed. Here, we report the discovery of a novel series of SIRT1/2 dual inhibitors via a rational drug design that involved virtual screening and a substructure search. Eleven of the derived compounds exhibited high inhibitory activities, with IC50 < 5 μM and high specificity for both SIRT1 and SIRT2. Compounds hsa55 and PS9 strongly induced apoptosis and showed antiproliferative effects against human leukemia cell lines, which could be due to their ability to increase of p53 and α-tubulin acetylation, as we observed in MOLM-13 cells. Therefore, the new scaffolds of these compounds and their efficacy in leukemia cell lines provide important clues for the further development of novel anti-leukemia drugs.

Structure-based discovery of novel cruzain inhibitors with distinct trypanocidal activity profiles

Over 110 years after the first formal description of Chagas disease, the trypanocidal drugs thus far available have limited efficacy and several side effects. This encourages the search for novel treatments that inhibit T. cruzi targets. One of the most studied anti-T. cruzi targets is the cysteine protease cruzain; it is associated with metacyclogenesis, replication, and invasion of the host cells. We used computational techniques to identify novel molecular scaffolds that act as cruzain inhibitors. First, with a docking-based virtual screening, we identified compound 8, a competitive cruzain inhibitor with a K_i of 4.6 μM. Then, aided by molecular dynamics simulations, cheminformatics, and docking, we identified the analog compound 22 with a K_i of 27 μM. Surprisingly, despite sharing the same isoquinoline scaffold, compound 8 presented higher trypanocidal activity against the epimastigote forms, while compound 22, against the trypomastigotes and amastigotes. Taken together, compounds 8 and 22 represent a promising scaffold for further development of trypanocidal compounds as drug candidates for treating Chagas disease.

In Silico Methodologies to Improve Antioxidants' Characterization from Marine Organisms

Marine organisms have been reported to be valuable sources of bioactive molecules that have found applications in different industrial fields. From organism sampling to the identification and bioactivity characterization of a specific compound, different steps are necessary, which are time- and cost-consuming. Thanks to the advent of the -omic era, numerous genome, metagenome, transcriptome, metatranscriptome, proteome and microbiome data have been reported and deposited in public databases. These advancements have been fundamental for the development of in silico strategies for basic and applied research. In silico studies represent a convenient and efficient approach to the bioactivity prediction of known and newly identified marine molecules, reducing the time and costs of "wet-lab" experiments. This review focuses on in silico approaches applied to bioactive molecule discoveries from marine organisms. When available, validation studies reporting a bioactivity assay to confirm the presence of an antioxidant molecule or enzyme are reported, as well. Overall, this review suggests that in silico approaches can offer a valuable alternative to most expensive approaches and proposes them as a little explored field in which to invest.