Artificial intelligence in virtual screening: models versus experiments

A computational odyssey: uncovering classical β-lactamase inhibitors in dry fruits

In the antibacterial arsenal, β-lactams have held a prominent position, but increasing resistance due to unauthorized use and genetic factors requires new strategies. Combining β-lactamase inhibitors with broad-spectrum β-lactams proves effective in combating this resistance. ESBL producers demand new inhibitors, leading to the exploration of plant-derived secondary metabolites for potent β-lactam antibiotics or alternative inhibitors. Using virtual screening, molecular docking, ADMET analysis, and molecular dynamic simulation, this study actively analyzed the inhibitory activity of figs, cashews, walnuts, and peanuts against SHV-1, NDM-1, KPC-2, and OXA-48 β-lactamases. Using AutoDock Vina, the docking affinities of various compounds for target enzymes were initially screened, revealing 12 bioactive compounds with higher affinities for the target enzymes compared to Avibactam and Tazobactam. Top-scoring metabolites, including Oleanolic acid, Protocatechuic acid, and Tannin, were subjected to MD simulation studies to further analyze the stability of the docked complexes using WebGro. The simulation coordinates, in terms of RMSD, RMSF, SASA, Rg, and hydrogen bonds formed, showed that these phytocompounds are stable enough to retain in the active sites at various orientations. The PCA and FEL analysis also showed the stability of the dynamic motion of Cα residues of phytochemical-bound enzymes. The pharmacokinetic analysis of the top phytochemicals was performed to analyze their bioavailability and toxicity. This study provides new insights into the therapeutic potential of phytochemicals of selected dry fruits and contributes to future experimental studies to identify βL inhibitors from plants.Communicated by Ramaswamy H. Sarma.

Protein Structure Inspired Drug Discovery

Drug discovery starts with known function, either of a compound or a protein, in-turn prompting investigations to probe 3D structure of the compound-protein interface. As protein structure determines function, we hypothesized that unique 3D structural motifs represent primary information denoting unique function that can drive discovery of novel agents. Using a physics-based protein structure analysis platform developed by us, designed to conduct computationally intensive analysis at supercomputing speeds, we probed a high-resolution protein x-ray crystallographic library developed by us. We selected 3D structural motifs whose function was not otherwise established, that offered environments supporting binding of drug-like chemicals and were present on proteins that were not established therapeutic targets. For each of eight potential binding pockets on six different proteins we accessed a 60 million compound library and used our analysis platform to evaluate binding. Using eight-day colony formation assays acquired compounds were screened for efficacy against human breast, prostate, colon and lung cancer cells and toxicity against human bone marrow stem cells. Compounds selectively inhibiting cancer growth segregated to two pockets on separate proteins. The compound, Dxr2-017, exhibited selective activity against human melanoma cells in the NCI-60 cell line screen, had an IC50 of 19 nM against human melanoma M14 cells in our eight-day assay, while over 2100-fold higher concentrations inhibited stem cells by less than 30%. We show that Dxr2-017 induces anoikis, a unique form of programmed cell death in need of targeted therapeutics. The predicted target protein for Dxr2-017 is expressed in bacteria, not in humans. This supports our strategy of focusing on unique 3D structural motifs. It is known that functionally important 3D structures are evolutionarily conserved. Here we demonstrate proof-of-concept that protein structure represents high value primary data to support discovery of novel therapeutics. This approach is widely applicable.

Accelerating BRPF1b hit identification with BioPhysical and Active Learning Screening (BioPALS)

We report the development of BioPhysical and Active Learning Screening (BioPALS); a rapid and versatile hit identification protocol combining AI-powered virtual screening with a GCI-driven biophysical confirmation workflow. Its application to the BRPF1b bromodomain afforded a range of novel micromolar binders with favorable ADMET properties. In addition to the excellent in silico/in vitro confirmation rate demonstrated with BRPF1b, binding kinetics were determined, and binding topologies predicted for all hits. BioPALS is a lean, data-rich, and standardized approach to hit identification applicable to a wide range of biological targets.

Molecular Property Diagnostic Suite for COVID-19 (MPDSCOVID-19): an open-source disease-specific drug discovery portal

Molecular Property Diagnostic Suite (MPDS) was conceived and developed as an open-source disease-specific web portal based on Galaxy. MPDSCOVID-19 was developed for COVID-19 as a one-stop solution for drug discovery research. Galaxy platforms enable the creation of customized workflows connecting various modules in the web server. The architecture of MPDSCOVID-19 effectively employs Galaxy v22.04 features, which are ported on CentOS 7.8 and Python 3.7. MPDSCOVID-19 provides significant updates and the addition of several new tools updated after six years. Tools developed by our group in Perl/Python and open-source tools are collated and integrated into MPDSCOVID-19 using XML scripts. Our MPDS suite aims to facilitate transparent and open innovation. This approach significantly helps bring inclusiveness in the community while promoting free access and participation in software development.

Availability & Implementation

The MPDSCOVID-19 portal can be accessed at https://mpds.neist.res.in:8085/.

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.

Virtual Screening Process: A Guide in Modern Drug Designing

Due to its capacity to drastically cut the cost and time necessary for experimental screening of compounds, virtual screening (VS) has grown to be a crucial component of drug discovery and development. VS is a computational method used in drug design to identify potential drugs from enormous libraries of chemicals. This approach makes use of molecular modeling and docking simulations to assess the small molecule's ability to bind to the desired protein. Virtual screening has a bright future, as high computational power and modern techniques are likely to further enhance the accuracy and speed of the process.

ComplexEye: a multi-lens array microscope for high-throughput embedded immune cell migration analysis

Autonomous migration is essential for the function of immune cells such as neutrophils and plays an important role in numerous diseases. The ability to routinely measure or target it would offer a wealth of clinical applications. Video microscopy of live cells is ideal for migration analysis, but cannot be performed at sufficiently high-throughput (HT). Here we introduce ComplexEye, an array microscope with 16 independent aberration-corrected glass lenses spaced at the pitch of a 96-well plate to produce high-resolution movies of migrating cells. With the system, we enable HT migration analysis of immune cells in 96- and 384-well plates with very energy-efficient performance. We demonstrate that the system can measure multiple clinical samples simultaneously. Furthermore, we screen 1000 compounds and identify 17 modifiers of migration in human neutrophils in just 4 days, a task that requires 60-times longer with a conventional video microscope. ComplexEye thus opens the field of phenotypic HT migration screens and enables routine migration analysis for the clinical setting.

The latest perspectives of small molecules FMS kinase inhibitors

FMS kinase is a type III tyrosine kinase receptor that plays a central role in the pathophysiology and management of several diseases, including a range of cancer types, inflammatory disorders, neurodegenerative disorders, and bone disorders among others. In this review, the pathophysiological pathways of FMS kinase in different diseases and the recent developments of its monoclonal antibodies and inhibitors during the last five years are discussed. The biological and biochemical features of these inhibitors, including binding interactions, structure-activity relationships (SAR), selectivity, and potencies are discussed. The focus of this article is on the compounds that are promising leads and undergoing advanced clinical investigations, as well as on those that received FDA approval. In this article, we attempt to classify the reviewed FMS inhibitors according to their core chemical structure including pyridine, pyrrolopyridine, pyrazolopyridine, quinoline, and pyrimidine derivatives.

Artificial intelligence revolutionizing drug development: Exploring opportunities and challenges

By harnessing artificial intelligence (AI) algorithms and machine learning techniques, the entire drug discovery process stands to undergo a profound transformation, offering a myriad of advantages. Foremost among these is the ability of AI to conduct swift and efficient screenings of expansive compound libraries, significantly augmenting the identification of potential drug candidates. Moreover, AI algorithms can prove instrumental in predicting the efficacy and safety profiles of candidate compounds, thus endowing invaluable insights and reducing reliance on extensive preclinical and clinical testing. This predictive capacity of AI has the potential to streamline the drug development pipeline and enhance the success rate of clinical trials, ultimately resulting in the emergence of more efficacious and safer therapeutic agents. However, the deployment of AI in drug discovery introduces certain challenges that warrant attention. A primary hurdle entails the imperative acquisition of high-quality and diverse data. Furthermore, ensuring the interpretability of AI models assumes critical importance in securing regulatory endorsement and cultivating trust within scientific and medical communities. Addressing ethical considerations, including data privacy and mitigating bias, represents an additional momentous challenge, requiring assiduous navigation. In this review, we provide an intricate and comprehensive overview of the multifaceted challenges intrinsic to conventional drug development paradigms, while simultaneously interrogating the efficacy of AI in effectively surmounting these formidable obstacles.

Machine Learning Approaches Identify Chemical Features for Stage-Specific Antimalarial Compounds

Efficacy data from diverse chemical libraries, screened against the various stages of the malaria parasite Plasmodium falciparum, including asexual blood stage (ABS) parasites and transmissible gametocytes, serve as a valuable reservoir of information on the chemical space of compounds that are either active (or not) against the parasite. We postulated that this data can be mined to define chemical features associated with the sole ABS activity and/or those that provide additional life cycle activity profiles like gametocytocidal activity. Additionally, this information could provide chemical features associated with inactive compounds, which could eliminate any future unnecessary screening of similar chemical analogs. Therefore, we aimed to use machine learning to identify the chemical space associated with stage-specific antimalarial activity. We collected data from various chemical libraries that were screened against the asexual (126 374 compounds) and sexual (gametocyte) stages of the parasite (93 941 compounds), calculated the compounds' molecular fingerprints, and trained machine learning models to recognize stage-specific active and inactive compounds. We were able to build several models that predict compound activity against ABS and dual activity against ABS and gametocytes, with Support Vector Machines (SVM) showing superior abilities with high recall (90 and 66%) and low false-positive predictions (15 and 1%). This allowed the identification of chemical features enriched in active and inactive populations, an important outcome that could be mined for essential chemical features to streamline hit-to-lead optimization strategies of antimalarial candidates. The predictive capabilities of the models held true in diverse chemical spaces, indicating that the ML models are therefore robust and can serve as a prioritization tool to drive and guide phenotypic screening and medicinal chemistry programs.

Comparative Study of Single-stranded Oligonucleotides Secondary Structure Prediction Tools

BACKGROUND

Single-stranded nucleic acids (ssNAs) have important biological roles and a high biotechnological potential linked to their ability to bind to numerous molecular targets. This depends on the different spatial conformations they can assume. The first level of ssNAs spatial organisation corresponds to their base pairs pattern, i.e. their secondary structure. Many computational tools have been developed to predict the ssNAs secondary structures, making the choice of the appropriate tool difficult, and an up-to-date guide on the limits and applicability of current secondary structure prediction tools is missing. Therefore, we performed a comparative study of the performances of 9 freely available tools (mfold, RNAfold, CentroidFold, CONTRAfold, MC-Fold, LinearFold, UFold, SPOT-RNA, and MXfold2) on a dataset of 538 ssNAs with known experimental secondary structure.

RESULTS

The minimum free energy-based tools, namely mfold and RNAfold, and some tools based on artificial intelligence, namely CONTRAfold and MXfold2, provided the best results, with [Formula: see text] of exact predictions, whilst MC-fold seemed to be the worst performing tool, with only [Formula: see text] of exact predictions. In addition, UFold and SPOT-RNA are the only options for pseudoknots prediction. Including in the analysis of mfold and RNAfold results 5-10 suboptimal solutions further improved the performances of these tools. Nevertheless, we could observe issues in predicting particular motifs, such as multiple-ways junctions and mini-dumbbells, or the ssNAs whose structure has been determined in complex with a protein. In addition, our benchmark shows that some effort has to be paid for ssDNA secondary structure predictions.

CONCLUSIONS

In general, Mfold, RNAfold, and MXfold2 seem to currently be the best choice for the ssNAs secondary structure prediction, although they still show some limits linked to specific structural motifs. Nevertheless, actual trends suggest that artificial intelligence has a high potential to overcome these remaining issues, for example the recently developed UFold and SPOT-RNA have a high success rate in predicting pseudoknots.

GPCNDTA: Prediction of drug-target binding affinity through cross-attention networks augmented with graph features and pharmacophores

Drug-target affinity prediction is a challenging task in drug discovery. The latest computational models have limitations in mining edge information in molecule graphs, accessing to knowledge in pharmacophores, integrating multimodal data of the same biomolecule and realizing effective interactions between two different biomolecules. To solve these problems, we proposed a method called Graph features and Pharmacophores augmented Cross-attention Networks based Drug-Target binding Affinity prediction (GPCNDTA). First, we utilized the GNN module, the linear projection unit and self-attention layer to correspondingly extract features of drugs and proteins. Second, we devised intramolecular and intermolecular cross-attention to respectively fuse and interact features of drugs and proteins. Finally, the linear projection unit was applied to gain final features of drugs and proteins, and the Multi-Layer Perceptron was employed to predict drug-target binding affinity. Three major innovations of GPCNDTA are as follows: (i) developing the residual CensNet and the residual EW-GCN to correspondingly extract features of drug and protein graphs, (ii) regarding pharmacophores as a new type of priors to heighten drug-target affinity prediction performance, and (iii) devising intramolecular and intermolecular cross-attention, in which the intramolecular cross-attention realizes the effective fusion of different modal data related to the same biomolecule, and the intermolecular cross-attention fulfills the information interaction between two different biomolecules in attention space. The test results on five benchmark datasets imply that GPCNDTA achieves the best performance compared with state-of-the-art computational models. Besides, relying on ablation experiments, we proved effectiveness of GNN modules, pharmacophores and two cross-attention strategies in improving the prediction accuracy, stability and reliability of GPCNDA. In case studies, we applied GPCNDTA to predict binding affinities between 3C-like proteinase and 185 drugs, and observed that most binding affinities predicted by GPCNDTA are close to corresponding experimental measurements.

"You've got the Body I've got the Brains" - Could the current AI-based tools replace the human ingenuity for designing new drug candidates?

The emergence of artificial intelligence (AI) tools has transformed the landscape of drug discovery, providing unprecedented speed, efficiency, and cost-effectiveness in the search for new therapeutics. From target identification to drug formulation and delivery, AI-driven algorithms have revolutionized various aspects of medicinal chemistry, significantly accelerating the drug design process. Despite the transformative power of AI, this perspective article emphasizes the limitations of AI tools in drug discovery, requiring inventive skills of medicinal chemists. However, the article highlighted that there is a need for a harmonious integration of AI-based tools and human expertise in drug discovery. Such a synergistic approach promises to lead to groundbreaking therapies that address unmet medical needs and benefit humankind. As the world evolves technologically, the question remains: When will AI tools effectively design and develop drugs? The answer may lie in the seamless collaboration between AI and human researchers, unlocking transformative therapies that combat diseases effectively.

Molecular Property Diagnostic Suite Compound Library (MPDS-CL): a structure-based classification of the chemical space

Molecular Property Diagnostic Suite Compound Library (MPDS-CL) is an open-source Galaxy-based cheminformatics web portal which presents a structure-based classification of the molecules. A structure-based classification of nearly 150 million unique compounds, obtained from 42 publicly available databases and curated for redundancy removal through 97 hierarchically well-defined atom composition-based portions, has been done. These are further subjected to 56-bit fingerprint-based classification algorithm which led to the formation of 56 structurally well-defined classes. The classes thus obtained were further divided into clusters based on their molecular weight. Thus, the entire set of molecules was put into 56 different classes and 625 clusters. This led to the assignment of a unique ID, named as MPDS-AadharID, for each of these 149,169,443 molecules. MPDS-AadharID is akin to the unique number given to citizens in India (similar to SSN in the US and NINO in the UK). The unique features of MPDS-CL are (a) several search options, such as exact structure search, substructure search, property-based search, fingerprint-based search, using SMILES, InChIKey and key-in; (b) automatic generation of information for the processing for MPDS and other galaxy tools; (c) providing the class and cluster of a molecule which makes it easier and fast to search for similar molecules and (d) information related to the presence of the molecules in multiple databases. The MPDS-CL can be accessed at https://mpds.neist.res.in:8086/ .

Machine Learning Advances in Predicting Peptide/Protein-Protein Interactions Based on Sequence Information for Lead Peptides Discovery

Peptides have shown increasing advantages and significant clinical value in drug discovery and development. With the development of high-throughput technologies and artificial intelligence (AI), machine learning (ML) methods for discovering new lead peptides have been expanded and incorporated into rational drug design. Predictions of peptide-protein interactions (PepPIs) and protein-protein interactions (PPIs) are both opportunities and challenges in computational biology, which will help to better understand the mechanisms of disease and provide the impetus for the discovery of lead peptides. This paper comprehensively reviews computational models for PepPI and PPI predictions. It begins with an introduction of various databases of peptide ligands and target proteins. Then it discusses data formats and feature representations for proteins and peptides. Furthermore, classical ML methods and emerging deep learning (DL) methods that can be used to train prediction models of PepPI and PPI are classified into four categories, and their advantages and disadvantages are analyzed. To assess the relative performance of different models, different validation protocols and evaluation indexes are discussed. The goal of this review is to help researchers quickly get started to develop computational frameworks using these integrated resources and eventually promote the discovery of lead peptides.

Novel Benzothiazole Derivatives Synthesis and its Analysis as Diuretic Agents

Benzothiazoles, an anticonvulsant, antiviral, antihypertensive, and cancer-fighting medication of the heterocyclic scaffold family, also acts as antibacterial and antiviral agents. There is much interest in this chemical’s production because of the strong and vital biological action it possesses. Substituted aromatic aldehydes were combined with 2-amino-benzothiazole-6-sulfonic acid amides, or Schiff base derivatives, to create Schiff base derivatives. Recrystallized, characterized, and tested for diuretic efficacy in vivo using online tools, m.p. (melting point), Rf, FTIR (Fourier transform infrared), 1H-NMR (proton nuclear magnetic resonance) data The molecular characteristics of all the substances created were estimated using Lipinski’s rule of 5, OSIRIS (software) molecular property explorer, Molsoft, and Autodock 4.0 docking software. Male Wistar rats were used to make all the compounds traditionally in order to test for diuretic activity. Neither the elemental nor the spectral information for the synthesized compounds disagreed. There were five different methods used to evaluate these compounds: Lipinski rule of five, Molsoft to determine molecular characteristics, PASS (prediction of activity spectra for substances) values to determine the diuretic effect, and OSIRIS software to determine toxicology. In order to investigate the diuretic effects of the selected drugs, docking analysis was used. Acetazolamide was shown to have a diuretic effect that was superior to that of compounds IIIb and IIIe, whereas 2-{(E)-[(3-hydroxyphenyl)methylidene]amino}-1,3-benzothiazole-6-sulfonamide (IIIb) was found to be the most promising potential.

Artificial intelligence assisted identification of potential tau aggregation inhibitors: ligand- and structure-based virtual screening, in silico ADME, and molecular dynamics study

Alzheimer's disease (AD) is a severe, growing, multifactorial disorder affecting millions of people worldwide characterized by cognitive decline and neurodegeneration. The accumulation of tau protein into paired helical filaments is one of the major pathological hallmarks of AD and has gained the interest of researchers as a potential drug target to treat AD. Lately, Artificial Intelligence (AI) has revolutionized the drug discovery process by speeding it up and reducing the overall cost. As a part of our continuous effort to identify potential tau aggregation inhibitors, and leveraging the power of AI, in this study, we used a fully automated AI-assisted ligand-based virtual screening tool, PyRMD to screen a library of 12 million compounds from the ZINC database to identify potential tau aggregation inhibitors. The preliminary hits from virtual screening were filtered for similar compounds and pan-assay interference compounds (the compounds containing reactive functional groups which can interfere with the assays) using RDKit. Further, the selected compounds were prioritized based on their molecular docking score with the binding pocket of tau where the binding pockets were identified using replica exchange molecular dynamics simulation. Thirty-three compounds showing good docking scores for all the tau clusters were selected and were further subjected to in silico pharmacokinetic prediction. Finally, top 10 compounds were selected for molecular dynamics simulation and MMPBSA binding free energy calculations resulting in the identification of UNK_175, UNK_1027, UNK_1172, UNK_1173, UNK_1237, UNK_1518, and UNK_2181 as potential tau aggregation inhibitors.

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.

ZINC-22─A Free Multi-Billion-Scale Database of Tangible Compounds for Ligand Discovery

Purchasable chemical space has grown rapidly into the tens of billions of molecules, providing unprecedented opportunities for ligand discovery but straining the tools that might exploit these molecules at scale. We have therefore developed ZINC-22, a database of commercially accessible small molecules derived from multi-billion-scale make-on-demand libraries. The new database and tools enable analog searching in this vast new space via a facile GUI, CartBlanche, drawing on similarity methods that scale sublinearly in the number of molecules. The new library also uses data organization methods, enabling rapid lookup of molecules and their physical properties, including conformations, partial atomic charges, c Log P values, and solvation energies, all crucial for molecule docking, which had become slow with older database organizations in previous versions of ZINC. As the libraries have continued to grow, we have been interested in finding whether molecular diversity has suffered, for instance, because certain scaffolds have come to dominate via easy analoging. This has not occurred thus far, and chemical diversity continues to grow with database size, with a log increase in Bemis-Murcko scaffolds for every two-log unit increase in database size. Most new scaffolds come from compounds with the highest heavy atom count. Finally, we consider the implications for databases like ZINC as the libraries grow toward and beyond the trillion-molecule range. ZINC is freely available to everyone and may be accessed at cartblanche22.docking.org, via Globus, and in the Amazon AWS and Oracle OCI clouds.

From Traditional Ethnopharmacology to Modern Natural Drug Discovery: A Methodology Discussion and Specific Examples

Ethnopharmacology, through the description of the beneficial effects of plants, has provided an early framework for the therapeutic use of natural compounds. Natural products, either in their native form or after crude extraction of their active ingredients, have long been used by different populations and explored as invaluable sources for drug design. The transition from traditional ethnopharmacology to drug discovery has followed a straightforward path, assisted by the evolution of isolation and characterization methods, the increase in computational power, and the development of specific chemoinformatic methods. The deriving extensive exploitation of the natural product chemical space has led to the discovery of novel compounds with pharmaceutical properties, although this was not followed by an analogous increase in novel drugs. In this work, we discuss the evolution of ideas and methods, from traditional ethnopharmacology to in silico drug discovery, applied to natural products. We point out that, in the past, the starting point was the plant itself, identified by sustained ethnopharmacological research, with the active compound deriving after extensive analysis and testing. In contrast, in recent years, the active substance has been pinpointed by computational methods (in silico docking and molecular dynamics, network pharmacology), followed by the identification of the plant(s) containing the active ingredient, identified by existing or putative ethnopharmacological information. We further stress the potential pitfalls of recent in silico methods and discuss the absolute need for in vitro and in vivo validation as an absolute requirement. Finally, we present our contribution to natural products’ drug discovery by discussing specific examples, applying the whole continuum of this rapidly evolving field. In detail, we report the isolation of novel antiviral compounds, based on natural products active against influenza and SARS-CoV-2 and novel substances active on a specific GPCR, OXER1.