Protein-Ligand Docking in the Machine-Learning Era

SAnDReS 2.0: Development of machine-learning models to explore the scoring function space

Classical scoring functions may exhibit low accuracy in determining ligand binding affinity for proteins. The availability of both protein-ligand structures and affinity data make it possible to develop machine-learning models focused on specific protein systems with superior predictive performance. Here, we report a new methodology named SAnDReS that combines AutoDock Vina 1.2 with 54 regression methods available in Scikit-Learn to calculate binding affinity based on protein-ligand structures. This approach allows exploration of the scoring function space. SAnDReS generates machine-learning models based on crystal, docked, and AlphaFold-generated structures. As a proof of concept, we examine the performance of SAnDReS-generated models in three case studies. For all three cases, our models outperformed classical scoring functions. Also, SAnDReS-generated models showed predictive performance close to or better than other machine-learning models such as KDEEP, CSM-lig, and ΔVinaRF20. SAnDReS 2.0 is available to download at https://github.com/azevedolab/sandres.

Computational toolbox for the analysis of protein-glycan interactions

Protein-glycan interactions play pivotal roles in numerous biological processes, ranging from cellular recognition to immune response modulation. Understanding the intricate details of these interactions is crucial for deciphering the molecular mechanisms underlying various physiological and pathological conditions. Computational techniques have emerged as powerful tools that can help in drawing, building and visualising complex biomolecules and provide insights into their dynamic behaviour at atomic and molecular levels. This review provides an overview of the main computational tools useful for studying biomolecular systems, particularly glycans, both in free state and in complex with proteins, also with reference to the principles, methodologies, and applications of all-atom molecular dynamics simulations. Herein, we focused on the programs that are generally employed for preparing protein and glycan input files to execute molecular dynamics simulations and analyse the corresponding results. The presented computational toolbox represents a valuable resource for researchers studying protein-glycan interactions and incorporates advanced computational methods for building, visualising and predicting protein/glycan structures, modelling protein-ligand complexes, and analyse MD outcomes. Moreover, selected case studies have been reported to highlight the importance of computational tools in studying protein-glycan systems, revealing the capability of these tools to provide valuable insights into the binding kinetics, energetics, and structural determinants that govern specific molecular interactions.

SPRank─A Knowledge-Based Scoring Function for RNA-Ligand Pose Prediction and Virtual Screening

The growing interest in RNA-targeted drugs underscores the need for computational modeling of interactions between RNA molecules and small compounds. Having a reliable scoring function for RNA-ligand interactions is essential for effective computational drug screening. An ideal scoring function should not only predict the native pose for ligand binding but also rank the affinity of the binding for different ligands. However, existing scoring functions are primarily designed to predict the native binding modes for a given RNA-ligand pair and have not been thoroughly assessed for virtual screening purposes. In this paper, we introduce SPRank, a combination of machine-learning and knowledge-based scoring functions developed through a weighted iterative approach, specifically designed to tackle both binding mode prediction and virtual screening challenges. Our approach incorporates third-party docking software, such as rDock and AutoDock Vina, to sample flexible ligands against an ensemble of RNA structures, capturing the conformational flexibility of both the RNA and the ligand. Through rigorous testing, SPRank demonstrates improved performance compared to the tested scoring functions across four test sets comprising 122, 42, 55, and 71 nucleic acid-ligand complexes. Furthermore, SPRank exhibits improved performance in virtual screening tests targeting the HIV-1 TAR ensemble, which highlights its advantage in drug discovery. These results underscore the advantages of SPRank as a potentially promising tool for the RNA-targeted drug design. The source code of SPRank and the data sets are freely accessible at https://github.com/Vfold-RNA/SPRank.

Mechanistic prediction and validation of Brevilin A Therapeutic effects in Lung Cancer

BACKGROUND

Traditional Chinese medicine (TCM) has been found widespread application in neoplasm treatment, yielding promising therapeutic candidates. Previous studies have revealed the anti-cancer properties of Brevilin A, a naturally occurring sesquiterpene lactone derived from Centipeda minima (L.) A.Br. (C. minima), a TCM herb, specifically against lung cancer. However, the underlying mechanisms of its effects remain elusive. This study employs network pharmacology and experimental analyses to unravel the molecular mechanisms of Brevilin A in lung cancer.

METHODS

The Batman-TCM, Swiss Target Prediction, Pharmmapper, SuperPred, and BindingDB databases were screened to identify Brevilin A targets. Lung cancer-related targets were sourced from GEO, Genecards, OMIM, TTD, and Drugbank databases. Utilizing Cytoscape software, a protein-protein interaction (PPI) network was established. Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG), Gene set enrichment analysis (GSEA), and gene-pathway correlation analysis were conducted using R software. To validate network pharmacology results, molecular docking, molecular dynamics simulations, and in vitro experiments were performed.

RESULTS

We identified 599 Brevilin A-associated targets and 3864 lung cancer-related targets, with 155 overlapping genes considered as candidate targets for Brevilin A against lung cancer. The PPI network highlighted STAT3, TNF, HIF1A, PTEN, ESR1, and MTOR as potential therapeutic targets. GO and KEGG analyses revealed 2893 enriched GO terms and 157 enriched KEGG pathways, including the PI3K-Akt signaling pathway, FoxO signaling pathway, and HIF-1 signaling pathway. GSEA demonstrated a close association between hub genes and lung cancer. Gene-pathway correlation analysis indicated significant associations between hub genes and the cellular response to hypoxia pathway. Molecular docking and dynamics simulations confirmed Brevilin A's interaction with PTEN and HIF1A, respectively. In vitro experiments demonstrated Brevilin A-induced dose- and time-dependent cell death in A549 cells. Notably, Brevilin A treatment significantly reduced HIF-1α mRNA expression while increasing PTEN mRNA levels.

CONCLUSIONS

This study demonstrates that Brevilin A exerts anti-cancer effects in treating lung cancer through a multi-target and multi-pathway manner, with the HIF pathway potentially being involved. These results lay a theoretical foundation for the prospective clinical application of Brevilin A.

Computational investigation of 2, 4-Di Tert Butyl Phenol as alpha amylase inhibitor isolated from Coccinia grandis (L.) Voigt using molecular docking, and ADMET parameters

INTRODUCTION

Diabetes Mellitus is the metabolic disorder most prevalent globally, accounting for a substantial morbidity rate. The conventional drugs available for the management of diabetes are either expensive or lack the required efficacy. The purpose of this research is to isolate and characterize an active phytoconstituent from Coccinia grandis and assess its anti-diabetic properties.

METHODS AND MATERIALS

Stems of Coccinia grandis are subjected to successive extraction and isolation. The isolated compound by column chromatography was characterized by FTIR (fourier-transform infrared), 1 H NMR (proton nuclear magnetic resonance), and Mass spectroscopy. The antidiabetic potential of the isolated compound was evaluated by in-vitro alpha-amylase inhibitory activity. Further, the compound was subjected to molecular docking studies to study its interaction with the human pancreatic alpha-amylase (Molegro Virtual Docker) as well to determine the pharmacokinetic and toxicity profile using computational techniques (OSIRIS property explorer, Swiss ADME, pkCSM, and PreADMET).

RESULTS

The characterization of the compound suggests the structure to be 2,4-ditertiary butyl phenol. The in-vitro alpha-amylase inhibitory study indicated a concentration-dependent inhibition and the IC50 (median lethal dose) value of the isolated compound was found to be 64.36 μg/ml. The docking study with the A chain of receptor 5EMY yielded a favorable docking score of -81.48 Kcal mol-1, suggesting that the compound binds to the receptor with high affinity through electrostatic, hydrophobic, and hydrogen bonds. Furthermore, the silico ADME analysis of the compound revealed improved metabolism, a skin permeability of -3.87 cm/s, gastrointestinal absorption of 95.48 %, and a total clearance of 0.984 log ml min-1 kg-1. In silico toxicity analysis also predicted cutaneous irritations but no carcinogenicity, mutagenicity, or hepatotoxicity.

CONCLUSION

The data suggested that the isolated compound (2, 4-tertiary butyl phenol) has the potential to inhibit the alpha-amylase activity and possess optimal ADME properties as well as tolerable side effects.

Advances in machine-learning approaches to RNA-targeted drug design

RNA molecules play multifaceted functional and regulatory roles within cells and have garnered significant attention in recent years as promising therapeutic targets. With remarkable successes achieved by artificial intelligence (AI) in different fields such as computer vision and natural language processing, there is a growing imperative to harness AI's potential in computer-aided drug design (CADD) to discover novel drug compounds that target RNA. Although machine-learning (ML) approaches have been widely adopted in the discovery of small molecules targeting proteins, the application of ML approaches to model interactions between RNA and small molecule is still in its infancy. Compared to protein-targeted drug discovery, the major challenges in ML-based RNA-targeted drug discovery stem from the scarcity of available data resources. With the growing interest and the development of curated databases focusing on interactions between RNA and small molecule, the field anticipates a rapid growth and the opening of a new avenue for disease treatment. In this review, we aim to provide an overview of recent advancements in computationally modeling RNA-small molecule interactions within the context of RNA-targeted drug discovery, with a particular emphasis on methodologies employing ML techniques.

Geometry Optimization Algorithms in Conjunction with the Machine Learning Potential ANI-2x Facilitate the Structure-Based Virtual Screening and Binding Mode Prediction

Structure-based virtual screening utilizes molecular docking to explore and analyze ligand-macromolecule interactions, crucial for identifying and developing potential drug candidates. Although there is availability of several widely used docking programs, the accurate prediction of binding affinity and binding mode still presents challenges. In this study, we introduced a novel protocol that combines our in-house geometry optimization algorithm, the conjugate gradient with backtracking line search (CG-BS), which is capable of restraining and constraining rotatable torsional angles and other geometric parameters with a highly accurate machine learning potential, ANI-2x, renowned for its precise molecular energy predictions reassembling the wB97X/6-31G(d) model. By integrating this protocol with binding pose prediction using the Glide, we conducted additional structural optimization and potential energy prediction on 11 small molecule-macromolecule and 12 peptide-macromolecule systems. We observed that ANI-2x/CG-BS greatly improved the docking power, not only optimizing binding poses more effectively, particularly when the RMSD of the predicted binding pose by Glide exceeded around 5 Å, but also achieving a 26% higher success rate in identifying those native-like binding poses at the top rank compared to Glide docking. As for the scoring and ranking powers, ANI-2x/CG-BS demonstrated an enhanced performance in predicting and ranking hundreds or thousands of ligands over Glide docking. For example, Pearson's and Spearman's correlation coefficients remarkedly increased from 0.24 and 0.14 with Glide docking to 0.85 and 0.69, respectively, with the addition of ANI-2x/CG-BS for optimizing and ranking small molecules binding to the bacterial ribosomal aminoacyl-tRNA receptor. These results suggest that ANI-2x/CG-BS holds considerable potential for being integrated into virtual screening pipelines due to its enhanced docking performance.

Deciphering the Lexicon of Protein Targets: A Review on Multifaceted Drug Discovery in the Era of Artificial Intelligence

Understanding protein sequence and structure is essential for understanding protein-protein interactions (PPIs), which are essential for many biological processes and diseases. Targeting protein binding hot spots, which regulate signaling and growth, with rational drug design is promising. Rational drug design uses structural data and computational tools to study protein binding sites and protein interfaces to design inhibitors that can change these interactions, thereby potentially leading to therapeutic approaches. Artificial intelligence (AI), such as machine learning (ML) and deep learning (DL), has advanced drug discovery and design by providing computational resources and methods. Quantum chemistry is essential for drug reactivity, toxicology, drug screening, and quantitative structure-activity relationship (QSAR) properties. This review discusses the methodologies and challenges of identifying and characterizing hot spots and binding sites. It also explores the strategies and applications of artificial-intelligence-based rational drug design technologies that target proteins and protein-protein interaction (PPI) binding hot spots. It provides valuable insights for drug design with therapeutic implications. We have also demonstrated the pathological conditions of heat shock protein 27 (HSP27) and matrix metallopoproteinases (MMP2 and MMP9) and designed inhibitors of these proteins using the drug discovery paradigm in a case study on the discovery of drug molecules for cancer treatment. Additionally, the implications of benzothiazole derivatives for anticancer drug design and discovery are deliberated.

An ensemble-based approach to estimate confidence of predicted protein-ligand binding affinity values

When designing a machine learning-based scoring function, we access a limited number of protein-ligand complexes with experimentally determined binding affinity values, representing only a fraction of all possible protein-ligand complexes. Consequently, it is crucial to report a measure of confidence and quantify the uncertainty in the model's predictions during test time. Here, we adopt the conformal prediction technique to evaluate the confidence of a prediction for each member of the core set of the CASF 2016 benchmark. The conformal prediction technique requires a diverse ensemble of predictors for uncertainty estimation. To this end, we introduce ENS-Score as an ensemble predictor, which includes 30 models with different protein-ligand representation approaches and achieves Pearson's correlation of 0.842 on the core set of the CASF 2016 benchmark. Also, we comprehensively investigate the residual error of each data point to assess the normality behavior of the distribution of the residual errors and their correlation to the structural features of the ligands, such as hydrophobic interactions and halogen bonding. In the end, we provide a local host web application to facilitate the usage of ENS-Score. All codes to repeat results are provided at https://github.com/miladrayka/ENS_Score.

HIPPO: HIstogram-based Pseudo-POtential for scoring protein-ssRNA fragment-based docking poses

BACKGROUND

The RNA-Recognition motif (RRM) is a protein domain that binds single-stranded RNA (ssRNA) and is present in as much as 2% of the human genome. Despite this important role in biology, RRM-ssRNA interactions are very challenging to study on the structural level because of the remarkable flexibility of ssRNA. In the absence of atomic-level experimental data, the only method able to predict the 3D structure of protein-ssRNA complexes with any degree of accuracy is ssRNA'TTRACT, an ssRNA fragment-based docking approach using ATTRACT. However, since ATTRACT parameters are not ssRNA-specific and were determined in 2010, there is substantial opportunity for enhancement.

RESULTS

Here we present HIPPO, a composite RRM-ssRNA scoring potential derived analytically from contact frequencies in near-native versus non-native docking models. HIPPO consists of a consensus of four distinct potentials, each extracted from a distinct reference pool of protein-trinucleotide docking decoys. To score a docking pose with one potential, for each pair of RNA-protein coarse-grained bead types, each contact is awarded or penalised according to the relative frequencies of this contact distance range among the correct and incorrect poses of the reference pool. Validated on a fragment-based docking benchmark of 57 experimentally solved RRM-ssRNA complexes, HIPPO achieved a threefold or higher enrichment for half of the fragments, versus only a quarter with the ATTRACT scoring function. In particular, HIPPO drastically improved the chance of very high enrichment (12-fold or higher), a scenario where the incremental modelling of entire ssRNA chains from fragments becomes viable. However, for the latter result, more research is needed to make it directly practically applicable. Regardless, our approach already improves upon the state of the art in RRM-ssRNA modelling and is in principle extendable to other types of protein-nucleic acid interactions.

Deconvoluting low yield from weak potency in direct-to-biology workflows with machine learning

High throughput and rapid biological evaluation of small molecules is an essential factor in drug discovery and development. Direct-to-biology (D2B), whereby compound purification is foregone, has emerged as a viable technique in time efficient screening, specifically for PROTAC design and biological evaluation. However, one notable limitation is the prerequisite of high yielding reactions to ensure the desired compound is indeed the compound responsible for biological activity. Herein, we report a machine learning based yield-assay deconfounder capable of deconvoluting low yield from low potency to identify false negatives. We validated this approach by identifying promising SARS-CoV-2 main protease inhibitors with nanomolar activity that rivaled potency observed from the standard D2B workflow. Furthermore, we show how our framework can be utilized in a broad, in silico screen to produce compounds of similar potency as a D2B assay.

Drugging the entire human proteome: Are we there yet?

Each of the ∼20,000 proteins in the human proteome is a potential target for compounds that bind to it and modify its function. The 3D structures of most of these proteins are now available. Here, we discuss the prospects for using these structures to perform proteome-wide virtual HTS (VHTS). We compare physics-based (docking) and AI VHTS approaches, some of which are now being applied with large databases of compounds to thousands of targets. Although preliminary proteome-wide screens are now within our grasp, further methodological developments are expected to improve the accuracy of the results.

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.

Augmenting a training dataset of the generative diffusion model for molecular docking with artificial binding pockets

This study introduces the PocketCFDM generative diffusion model, aimed at improving the prediction of small molecule poses in the protein binding pockets. The model utilizes a novel data augmentation technique, involving the creation of numerous artificial binding pockets that mimic the statistical patterns of non-bond interactions found in actual protein-ligand complexes. An algorithmic method was developed to assess and replicate these interaction patterns in the artificial binding pockets built around small molecule conformers. It is shown that the integration of artificial binding pockets into the training process significantly enhanced the model's performance. Notably, PocketCFDM surpassed DiffDock in terms of non-bond interaction and steric clash numbers, and the inference speed. Future developments and optimizations of the model are discussed. The inference code and final model weights of PocketCFDM are accessible publicly via the GitHub repository: https://github.com/vtarasv/pocket-cfdm.git.

Computer-Aided Drug Discovery and Design: Recent Advances and Future Prospects

Computer-aided drug discovery and design involve the use of information technologies to identify and develop, on a rational ground, chemical compounds that align a set of desired physicochemical and biological properties. In its most common form, it involves the identification and/or modification of an active scaffold (or the combination of known active scaffolds), although de novo drug design from scratch is also possible. Traditionally, the drug discovery and design processes have focused on the molecular determinants of the interactions between drug candidates and their known or intended pharmacological target(s). Nevertheless, in modern times, drug discovery and design are conceived as a particularly complex multiparameter optimization task, due to the complicated, often conflicting, property requirements.This chapter provides an updated overview of in silico approaches for identifying active scaffolds and guiding the subsequent optimization process. Recent groundbreaking advances in the field have also analyzed the integration of state-of-the-art machine learning approaches in every step of the drug discovery process (from prediction of target structure to customized molecular docking scoring functions), integration of multilevel omics data, and the use of a diversity of computational approaches to assist target validation and assess plausible binding pockets.

Designing drugs and chemical probes with the dualsteric approach

Traditionally, drugs are monovalent, targeting only one site on the protein surface. This includes orthosteric and allosteric drugs, which bind the protein at orthosteric and allosteric sites, respectively. Orthosteric drugs are good in potency, whereas allosteric drugs have better selectivity and are solutions to classically undruggable targets. However, it would be difficult to simultaneously reach high potency and selectivity when targeting only one site. Also, both kinds of monovalent drugs suffer from mutation-caused drug resistance. To overcome these obstacles, dualsteric modulators have been proposed in the past twenty years. Compared to orthosteric or allosteric drugs, dualsteric modulators are bivalent (or bitopic) with two pharmacophores. Each of the two pharmacophores bind the protein at the orthosteric and an allosteric site, which could bring the modulator with special properties beyond monovalent drugs. In this study, we comprehensively review the current development of dualsteric modulators. Our main effort reason and illustrate the aims to apply the dualsteric approach, including a "double win" of potency and selectivity, overcoming mutation-caused drug resistance, developments of function-biased modulators, and design of partial agonists. Moreover, the strengths of the dualsteric technique also led to its application outside pharmacy, including the design of highly sensitive fluorescent tracers and usage as molecular rulers. Besides, we also introduced drug targets, designing strategies, and validation methods of dualsteric modulators. Finally, we detail the conclusions and perspectives.

In Silico Screening of Natural Flavonoids against 3-Chymotrypsin-like Protease of SARS-CoV-2 Using Machine Learning and Molecular Modeling

The "Long-COVID syndrome" has posed significant challenges due to a lack of validated therapeutic options. We developed a novel multi-step virtual screening strategy to reliably identify inhibitors against 3-chymotrypsin-like protease of SARS-CoV-2 from abundant flavonoids, which represents a promising source of antiviral and immune-boosting nutrients. We identified 57 interacting residues as contributors to the protein-ligand binding pocket. Their energy interaction profiles constituted the input features for Machine Learning (ML) models. The consensus of 25 classifiers trained using various ML algorithms attained 93.9% accuracy and a 6.4% false-positive-rate. The consensus of 10 regression models for binding energy prediction also achieved a low root-mean-square error of 1.18 kcal/mol. We screened out 120 flavonoid hits first and retained 50 drug-like hits after predefined ADMET filtering to ensure bioavailability and safety profiles. Furthermore, molecular dynamics simulations prioritized nine bioactive flavonoids as promising anti-SARS-CoV-2 agents exhibiting both high structural stability (root-mean-square deviation < 5 Å for 218 ns) and low MM/PBSA binding free energy (<-6 kcal/mol). Among them, KB-2 (PubChem-CID, 14630497) and 9-O-Methylglyceofuran (PubChem-CID, 44257401) displayed excellent binding affinity and desirable pharmacokinetic capabilities. These compounds have great potential to serve as oral nutraceuticals with therapeutic and prophylactic properties as care strategies for patients with long-COVID syndrome.

The discovery of subunit-selective GluN1/GluN2B NMDAR antagonist via pharmacophere-based virtual screening

The incidence and mortality rates of neurodegenerative diseases, such as Alzheimer's disease and Parkinson's disease, are gradually increasing worldwide. Numerous studies have demonstrated that N-methyl-D-aspartic acid receptor (NMDAR)-mediated excitotoxicity contributes to neurodegenerative diseases. Ifenprodil, a subtype-selective NMDAR antagonist, showed strong therapeutic potential. However, it suffers from low oral bioavailability and off-target side effects. In this study, natural compounds were identified for selective inhibition of GluN1/GluN2B NMDAR of human. We obtained a set of natural compounds (n = 81,366) from COCONUT, TIPdb, PAMDB, CMNPD, YMDB, and NPAtlas databases, and then virtually screened by an ifenprodil-structure-based pharmacophore model and molecular docking. The top 100 compounds were selected for binding affinity prediction via batch drug-target affinity (BatchDTA). Then, the top 50 compounds were analyzed by absorption, distribution, metabolism, excretion, toxicity (ADMET). Molecular dynamics involving molecular mechanics/position-Boltzmann surface area (MM-PBSA) calculation were performed to further screening. The top 15 compounds with strong binding affinity and ifenprodil, a proven GluN2B-selective NMDAR antagonist, were subjected to molecular dynamic simulations (100 ns), root-mean-square deviation (RMSD), root-mean-square fluctuation (RMSF), radius of gyration (Rg), H-bonds, solvent accessible surface area (SASA), principal component analysis (PCA), and binding free energy calculations. Based on these analyses, one possible lead compound carrying positive charges (CNP0099440) was identified, with great binding affinity and less off-target activity by contrast to ifenprodil. CNP0099440 has great potential to be a GluN1/GluN2B NMDAR antagonist candidate and can be further detected via in vitro and in vivo experiments.

Integrated Molecular Modeling and Machine Learning for Drug Design

Modern therapeutic development often involves several stages that are interconnected, and multiple iterations are usually required to bring a new drug to the market. Computational approaches have increasingly become an indispensable part of helping reduce the time and cost of the research and development of new drugs. In this Perspective, we summarize our recent efforts on integrating molecular modeling and machine learning to develop computational tools for modulator design, including a pocket-guided rational design approach based on AlphaSpace to target protein-protein interactions, delta machine learning scoring functions for protein-ligand docking as well as virtual screening, and state-of-the-art deep learning models to predict calculated and experimental molecular properties based on molecular mechanics optimized geometries. Meanwhile, we discuss remaining challenges and promising directions for further development and use a retrospective example of FDA approved kinase inhibitor Erlotinib to demonstrate the use of these newly developed computational tools.

Anti-obesity effects of olivetol in adult zebrafish model induced by short-term high-fat diet

Obesity is a complex disease caused by various factors, and synthetic drugs used to treat it can have side effects. Natural compounds, such as olivetol, could be a promising alternative. Olivetol is a substance found in certain lichen species and has anti-inflammatory and anti-cancer properties. In this study, researchers conducted in-silico molecular docking studies and found that olivetol had significant binding affinity with receptors involved in obesity. They also investigated the effects of olivetol on a diet-induced obese zebrafish model and found that high doses of olivetol reduced excessive fat accumulation and triglyceride and lipid accumulation. The low dose of olivetol showed a significant reduction in liver enzymes' levels. However, the high dose of olivetol resulted in a significant increase in HMG-CoA levels. These results suggest that olivetol may be a promising anti-obesity agent for the treatment of hyperlipidemia-related disorders, but further research is necessary to understand its full effects on the body.

Multi-shelled ECIF: improved extended connectivity interaction features for accurate binding affinity prediction

Motivation

Extended connectivity interaction features (ECIF) is a method developed to predict protein-ligand binding affinity, allowing for detailed atomic representation. It performed very well in terms of Comparative Assessment of Scoring Functions 2016 (CASF-2016) scoring power. However, ECIF has the limitation of not being able to adequately account for interatomic distances.

Results

To investigate what kind of distance representation is effective for P-L binding affinity prediction, we have developed two algorithms that improved ECIF's feature extraction method to take distance into account. One is multi-shelled ECIF, which takes into account the distance between atoms by dividing the distance between atoms into multiple layers. The other is weighted ECIF, which weights the importance of interactions according to the distance between atoms. A comparison of these two methods shows that multi-shelled ECIF outperforms weighted ECIF and the original ECIF, achieving a CASF-2016 scoring power Pearson correlation coefficient of 0.877.

Availability and implementation

All the codes and data are available on GitHub (https://github.com/koji11235/MSECIFv2).

SubMDTA: drug target affinity prediction based on substructure extraction and multi-scale features

BACKGROUND

Drug-target affinity (DTA) prediction is a critical step in the field of drug discovery. In recent years, deep learning-based methods have emerged for DTA prediction. In order to solve the problem of fusion of substructure information of drug molecular graphs and utilize multi-scale information of protein, a self-supervised pre-training model based on substructure extraction and multi-scale features is proposed in this paper.

RESULTS

For drug molecules, the model obtains substructure information through the method of probability matrix, and the contrastive learning method is implemented on the graph-level representation and subgraph-level representation to pre-train the graph encoder for downstream tasks. For targets, a BiLSTM method that integrates multi-scale features is used to capture long-distance relationships in the amino acid sequence. The experimental results showed that our model achieved better performance for DTA prediction.

CONCLUSIONS

The proposed model improves the performance of the DTA prediction, which provides a novel strategy based on substructure extraction and multi-scale features.

Recent Developments in Ultralarge and Structure-Based Virtual Screening Approaches

Drug development is a wide scientific field that faces many challenges these days. Among them are extremely high development costs, long development times, and a small number of new drugs that are approved each year. New and innovative technologies are needed to solve these problems that make the drug discovery process of small molecules more time and cost efficient, and that allow previously undruggable receptor classes to be targeted, such as protein-protein interactions. Structure-based virtual screenings (SBVSs) have become a leading contender in this context. In this review, we give an introduction to the foundations of SBVSs and survey their progress in the past few years with a focus on ultralarge virtual screenings (ULVSs). We outline key principles of SBVSs, recent success stories, new screening techniques, available deep learning-based docking methods, and promising future research directions. ULVSs have an enormous potential for the development of new small-molecule drugs and are already starting to transform early-stage drug discovery.

Molecular docking in organic, inorganic, and hybrid systems: a tutorial review

Molecular docking simulation is a very popular and well-established computational approach and has been extensively used to understand molecular interactions between a natural organic molecule (ideally taken as a receptor) such as an enzyme, protein, DNA, RNA and a natural or synthetic organic/inorganic molecule (considered as a ligand). But the implementation of docking ideas to synthetic organic, inorganic, or hybrid systems is very limited with respect to their use as a receptor despite their huge popularity in different experimental systems. In this context, molecular docking can be an efficient computational tool for understanding the role of intermolecular interactions in hybrid systems that can help in designing materials on mesoscale for different applications. The current review focuses on the implementation of the docking method in organic, inorganic, and hybrid systems along with examples from different case studies. We describe different resources, including databases and tools required in the docking study and applications. The concept of docking techniques, types of docking models, and the role of different intermolecular interactions involved in the docking process to understand the binding mechanisms are explained. Finally, the challenges and limitations of dockings are also discussed in this review.

Graphical abstract

Recent Advances in Computer-Aided Structure-Based Drug Design on Ion Channels

Ion channels play important roles in fundamental biological processes, such as electric signaling in cells, muscle contraction, hormone secretion, and regulation of the immune response. Targeting ion channels with drugs represents a treatment option for neurological and cardiovascular diseases, muscular degradation disorders, and pathologies related to disturbed pain sensation. While there are more than 300 different ion channels in the human organism, drugs have been developed only for some of them and currently available drugs lack selectivity. Computational approaches are an indispensable tool for drug discovery and can speed up, especially, the early development stages of lead identification and optimization. The number of molecular structures of ion channels has considerably increased over the last ten years, providing new opportunities for structure-based drug development. This review summarizes important knowledge about ion channel classification, structure, mechanisms, and pathology with the main focus on recent developments in the field of computer-aided, structure-based drug design on ion channels. We highlight studies that link structural data with modeling and chemoinformatic approaches for the identification and characterization of new molecules targeting ion channels. These approaches hold great potential to advance research on ion channel drugs in the future.

Design and synthesis of multi-functional small-molecule based inhibitors of amyloid-β aggregation: Molecular modeling and in vitro evaluation

Amyloid-β1-42 (Aβ42) peptide aggregate formation in the brain plays a crucial role in the onset and progression of Alzheimer's disease. According to published research, the Aβ monomer's amino acid residues KLVFF (16-20) self-associate to create antiparallel β-sheet fibrils. Small compounds can prevent self-assembly and destroy Aβ fibrils by attaching to the Aβ16-20 regions of Aβ42. To enhance biological characteristics and binding affinity to the amyloid beta peptide, β-sheet breaker small molecules can be developed and modified with various scaffolds. In the current study, a novel series of 2,3-disubstitutedbenzofuran derivatives was designed and created by fusing the benzofuran core of a known iron chelator, neuroprotective, and neurorestorative agent, like VK-28, with a motif found in the structure of a known muscarinic inhibitor and amyloid binding agent, like SKF-64346. Measurements of the binding affinity and in vitro aggregation inhibition of the Aβ42 peptide were made using the thioflavin T (ThT) test. Using AutoDock 4.2 software, molecular docking studies of the synthesized compounds were performed on the monomer and fibril structures of amyloid beta peptide. The compounds 8a-8g exhibited strong binding energy and affinity to Aβ fibrils as well as a 50%-67% reduction of the growth of Aβ aggregation. Finally, the positive traits of our recently synthesized compounds make them excellent candidates for additional in vivo testing as a "β-sheet breaking agent."

Vaccinomics Approach for Multi-Epitope Vaccine Design against Group A Rotavirus Using VP4 and VP7 Proteins

Rotavirus A is the most common cause of Acute Gastroenteritis globally among children <5 years of age. Due to a segmented genome, there is a high frequency of genetic reassortment and interspecies transmission which has resulted in the emergence of novel genotypes. There are concerns that monovalent (Rotarix: GlaxoSmithKline Biologicals, Rixensart, Belgium) and pentavalent (RotaTeq: MERCK & Co., Inc., Kenilworth, NJ, USA) vaccines may be less effective against non-vaccine strains, which clearly shows the demand for the design of a vaccine that is equally effective against all circulating genotypes. In the present study, a multivalent vaccine was designed from VP4 and VP7 proteins of RVA. Epitopes were screened for antigenicity, allergenicity, homology with humans and anti-inflammatory properties. The vaccine contains four B-cell, three CTL and three HTL epitopes joined via linkers and an N-terminal RGD motif adjuvant. The 3D structure was predicted and refined preceding its docking with integrin. Immune simulation displayed promising results both in Asia and worldwide. In the MD simulation, the RMSD value varied from 0.2 to 1.6 nm while the minimum integrin amino acid fluctuation (0.05-0.1 nm) was observed with its respective ligand. Codon optimization was performed with an adenovirus vector in a mammalian expression system. The population coverage analysis showed 99.0% and 98.47% in South Asia and worldwide, respectively. These computational findings show potential against all RVA genotypes; however, in-vitro/in-vivo screening is essential to devise a meticulous conclusion.

Pharmacological Chaperones and Protein Conformational Diseases: Approaches of Computational Structural Biology

Whenever a protein fails to fold into its native structure, a profound detrimental effect is likely to occur, and a disease is often developed. Protein conformational disorders arise when proteins adopt abnormal conformations due to a pathological gene variant that turns into gain/loss of function or improper localization/degradation. Pharmacological chaperones are small molecules restoring the correct folding of a protein suitable for treating conformational diseases. Small molecules like these bind poorly folded proteins similarly to physiological chaperones, bridging non-covalent interactions (hydrogen bonds, electrostatic interactions, and van der Waals contacts) loosened or lost due to mutations. Pharmacological chaperone development involves, among other things, structural biology investigation of the target protein and its misfolding and refolding. Such research can take advantage of computational methods at many stages. Here, we present an up-to-date review of the computational structural biology tools and approaches regarding protein stability evaluation, binding pocket discovery and druggability, drug repurposing, and virtual ligand screening. The tools are presented as organized in an ideal workflow oriented at pharmacological chaperones' rational design, also with the treatment of rare diseases in mind.

GB-score: Minimally designed machine learning scoring function based on distance-weighted interatomic contact features

In recent years, thanks to advances in computer hardware and dataset availability, data-driven approaches (like machine learning) have become one of the essential parts of the drug design framework to accelerate drug discovery procedures. Constructing a new scoring function, a function that can predict the binding score for a generated protein-ligand pose during docking procedure or a crystal complex, based on machine and deep learning has become an active research area in computer-aided drug design. GB-Score is a state-of-the-art machine learning-based scoring function that utilizes distance-weighted interatomic contact features, PDBbind-v2019 general set, and Gradient Boosting Trees algorithm to the binding affinity prediction. The distance-weighted interatomic contact featurization method used the distance between different ligand and protein atom types for numerical representation of the protein-ligand complex. GB-Score attains Pearson's correlation 0.862 and RMSE 1.190 on the CASF-2016 benchmark test in the scoring power metric. GB-Score's codes are freely available on the web at https://github.com/miladrayka/GB_Score.

In silico screening of phenylethanoid glycosides, a class of pharmacologically active compounds as natural inhibitors of SARS-CoV-2 proteases

Since the advent of Covid-19, several natural products have been investigated regarding their in silico interactions with SARS-CoV-2 proteases - 3CL^pro and PL^pro, two of the most important pharmacological targets for antiviral development. Phenylethanoid glycosides (PG) are a class of natural products present in important medicinal plants and a drug containing this group of active ingredients has been successfully used in the treatment of Covid-19 in China. Thus, a dataset with 567 derivatives of this class was built from reviews published between 1994 and 2020, and their interaction against both SARS-CoV-2 proteases was investigated. The virtual screening was performed by filtering the PGs through the evaluation of scores based on the AutoDock Vina, GOLD/ChemPLP, and GOLD/GoldScore evaluation functions. The bRO5 pharmacokinetic parameters of the PGs ranked in the previous step were analyzed and their interaction with key amino acid residues of the 3CL^pro and PL^pro enzymes was evaluated. Ninety-eight compounds were identified by computational approaches against PL^pro and 80 PGs against 3CL^pro. Of these, four interacted with key catalytic residues of PL^pro, which is an indicative of inhibitory activity, and three compounds interacted with catalytic key residues of 3CL^pro. Of these, five PGs occur in plants of the Traditional Chinese Medicine (TCM), while two are components of plants/formulations currently used in the Covid-19 protocols in China. The data presented here show the potential of PGs as selective inhibitors of SARS-CoV-2 3CL^pro and PL^pro.

Biomanufacturing by In Vitro Biotransformation (ivBT) Using Purified Cascade Multi-enzymes

In vitro biotransformation (ivBT) refers to the use of an artificial biological reaction system that employs purified enzymes for the one-pot conversion of low-cost materials into biocommodities such as ethanol, organic acids, and amino acids. Unshackled from cell growth and metabolism, ivBT exhibits distinct advantages compared with metabolic engineering, including but not limited to high engineering flexibility, ease of operation, fast reaction rate, high product yields, and good scalability. These characteristics position ivBT as a promising next-generation biomanufacturing platform. Nevertheless, challenges persist in the enhancement of bulk enzyme preparation methods, the acquisition of enzymes with superior catalytic properties, and the development of sophisticated approaches for pathway design and system optimization. In alignment with the workflow of ivBT development, this chapter presents a systematic introduction to pathway design, enzyme mining and engineering, system construction, and system optimization. The chapter also proffers perspectives on ivBT development.

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.

Assessing How Residual Errors of Scoring Functions Correlate to Ligand Structural Features

Scoring functions (SFs) are ubiquitous tools for early stage drug discovery. However, their accuracy currently remains quite moderate. Despite a number of successful target-specific SFs appearing recently, up until now, no ideas on how to systematically improve the general scope of SFs have been formulated. In this work, we hypothesized that the specific features of ligands, corresponding to interactions well appreciated by medicinal chemists (e.g., hydrogen bonds, hydrophobic and aromatic interactions), might be responsible, in part, for the remaining SF errors. The latter provides direction to efforts aimed at the rational and systematic improvement of SF accuracy. In this proof-of-concept work, we took a CASF-2016 coreset of 285 ligands as a basis for comparison and calculated the values of scores for a representative panel of SFs (including AutoDock 4.2, AutoDock Vina, X-Score, NNScore2.0, ΔVina RF20, and DSX). The residual error of linear correlation of each SF value, with the experimental values of affinity and activity, was then analyzed in terms of its correlation with the presence of the fragments responsible for certain medicinal chemistry defined interactions. We showed that, despite the fact that SFs generally perform reasonably, there is room for improvement in terms of better parameterization of interactions involving certain fragments in ligands. Thus, this approach opens a potential way for the systematic improvement of SFs without their significant complication. However, the straightforward application of the proposed approach is limited by the scarcity of reliable available data for ligand-receptor complexes, which is a common problem in the field.