Structure-based virtual screening: an overview

Targeting the glutamine-arginine-proline metabolism axis in cancer

Metabolic abnormalities are an important feature of tumours. The glutamine-arginine-proline axis is an important node of cancer metabolism and plays a major role in amino acid metabolism. This axis also acts as a scaffold for the synthesis of other nonessential amino acids and essential metabolites. In this paper, we briefly review (1) the glutamine addiction exhibited by tumour cells with accelerated glutamine transport and metabolism; (2) the methods regulating extracellular glutamine entry, intracellular glutamine synthesis and the fate of intracellular glutamine; (3) the glutamine, proline and arginine metabolic pathways and their interaction; and (4) the research progress in tumour therapy targeting the glutamine-arginine-proline metabolic system, with a focus on summarising the therapeutic research progress of strategies targeting of one of the key enzymes of this metabolic system, P5CS (ALDH18A1). This review provides a new basis for treatments targeting the metabolic characteristics of tumours.

Advances, opportunities, and challenges in methods for interrogating the structure activity relationships of natural products

Time span in literature: 1985-early 2024Natural products play a key role in drug discovery, both as a direct source of drugs and as a starting point for the development of synthetic compounds. Most natural products are not suitable to be used as drugs without further modification due to insufficient activity or poor pharmacokinetic properties. Choosing what modifications to make requires an understanding of the compound's structure-activity relationships. Use of structure-activity relationships is commonplace and essential in medicinal chemistry campaigns applied to human-designed synthetic compounds. Structure-activity relationships have also been used to improve the properties of natural products, but several challenges still limit these efforts. Here, we review methods for studying the structure-activity relationships of natural products and their limitations. Specifically, we will discuss how synthesis, including total synthesis, late-stage derivatization, chemoenzymatic synthetic pathways, and engineering and genome mining of biosynthetic pathways can be used to produce natural product analogs and discuss the challenges of each of these approaches. Finally, we will discuss computational methods including machine learning methods for analyzing the relationship between biosynthetic genes and product activity, computer aided drug design techniques, and interpretable artificial intelligence approaches towards elucidating structure-activity relationships from models trained to predict bioactivity from chemical structure. Our focus will be on these latter topics as their applications for natural products have not been extensively reviewed. We suggest that these methods are all complementary to each other, and that only collaborative efforts using a combination of these techniques will result in a full understanding of the structure-activity relationships of natural products.

Development of a novel representation of drug 3D structures and enhancement of the TSR-based method for probing drug and target interactions

Understanding the mechanisms underlying interactions between drugs and target proteins is critical for drug discovery. In our earlier studies, we introduced the Triangular Spatial Relationship (TSR)-based algorithm, which enables the representation of a protein's 3D structure as a vector of integers (TSR keys). These TSR keys correspond to substructures of the 3D structure of a protein and are computed based on the triangles constructed by all possible triples of Cα atoms within the protein. In this study, we report on a new TSR-based algorithm for probing drug and target interactions. Specifically, we have extended the previous algorithm in three novel directions: TSR keys for representing the 3D structure of a drug or a ligand, cross TSR keys between drugs and their targets and intra-residual TSR keys for phosphorylated amino acids. The outcomes illustrate the key contributions as follows: (i) The TSR-based method, which uses the TSR keys as features, is unique in its capability to interpret hierarchical relationships of drugs as well as drug - target complexes using common and specific TSR keys. (ii) The method can distinguish not only the binding sites from the rest of the protein structures, but also the binding sites of primary targets from those of off-targets. (iii) The method has the potential to correlate the 3D structures of drugs with their functions. (iv) Representation of 3D structures by TSR keys has its unique advantage in terms of ease of making searching for similar substructures across structure datasets easier. In summary, this study presents a novel computational methodology, with significant advantages, for providing insights into the mechanism underlying drug and target interactions.

EC-Conf: A ultra-fast diffusion model for molecular conformation generation with equivariant consistency

Despite recent advancement in 3D molecule conformation generation driven by diffusion models, its high computational cost in iterative diffusion/denoising process limits its application. Here, an equivariant consistency model (EC-Conf) was proposed as a fast diffusion method for low-energy conformation generation. In EC-Conf, a modified SE (3)-equivariant transformer model was directly used to encode the Cartesian molecular conformations and a highly efficient consistency diffusion process was carried out to generate molecular conformations. It was demonstrated that, with only one sampling step, it can already achieve comparable quality to other diffusion-based models running with thousands denoising steps. Its performance can be further improved with a few more sampling iterations. The performance of EC-Conf is evaluated on both GEOM-QM9 and GEOM-Drugs sets. Our results demonstrate that the efficiency of EC-Conf for learning the distribution of low energy molecular conformation is at least two magnitudes higher than current SOTA diffusion models and could potentially become a useful tool for conformation generation and sampling. SCIENTIFIC CONTRIBUTIONS: In this work, we proposed an equivariant consistency model that significantly improves the efficiency of conformation generation in diffusion-based models while maintaining high structural quality. This method serves as a general framework and can be further extended to more complex structure generation and prediction tasks, including those involving proteins, in future steps.

PIDiff: Physics informed diffusion model for protein pocket-specific 3D molecular generation

Designing drugs capable of binding to the structure of target proteins for treating diseases is essential in drug development. Recent remarkable advancements in geometric deep learning have led to unprecedented progress in three-dimensional (3D) generation of ligands that can bind to the protein pocket. However, most existing methods primarily focus on modeling the geometric information of ligands in 3D space. Consequently, these methods fail to consider that the binding of proteins and ligands is a phenomenon driven by intrinsic physicochemical principles. Motivated by this understanding, we propose PIDiff, a model for generating molecules by accounting in the physicochemical principles of protein-ligand binding. Our model learns not only the structural information of proteins and ligands but also to minimize the binding free energy between them. To evaluate the proposed model, we introduce an experimental framework that surpasses traditional assessment methods by encompassing various essential aspects for the practical application of generative models to actual drug development. The results confirm that our model outperforms baseline models on the CrossDocked2020 benchmark dataset, demonstrating its superiority. Through diverse experiments, we have illustrated the promising potential of the proposed model in practical drug development.

Artificial intelligence for small molecule anticancer drug discovery

INTRODUCTION

The transition from conventional cytotoxic chemotherapy to targeted cancer therapy with small-molecule anticancer drugs has enhanced treatment outcomes. This approach, which now dominates cancer treatment, has its advantages. Despite the regulatory approval of several targeted molecules for clinical use, challenges such as low response rates and drug resistance still persist. Conventional drug discovery methods are costly and time-consuming, necessitating more efficient approaches. The rise of artificial intelligence (AI) and access to large-scale datasets have revolutionized the field of small-molecule cancer drug discovery. Machine learning (ML), particularly deep learning (DL) techniques, enables the rapid identification and development of novel anticancer agents by analyzing vast amounts of genomic, proteomic, and imaging data to uncover hidden patterns and relationships.

AREA COVERED

In this review, the authors explore the important landmarks in the history of AI-driven drug discovery. They also highlight various applications in small-molecule cancer drug discovery, outline the challenges faced, and provide insights for future research.

EXPERT OPINION

The advent of big data has allowed AI to penetrate and enable innovations in almost every stage of medicine discovery, transforming the landscape of oncology research through the development of state-of-the-art algorithms and models. Despite challenges in data quality, model interpretability, and technical limitations, advancements promise breakthroughs in personalized and precision oncology, revolutionizing future cancer management.

Efficient molecular conformation generation with quantum-inspired algorithm

CONTEXT

Conformation generation, also known as molecular unfolding (MU), is a crucial step in structure-based drug design, remaining a challenging combinatorial optimization problem. Quantum annealing (QA) has shown great potential for solving certain combinatorial optimization problems over traditional classical methods such as simulated annealing (SA). However, a recent study showed that a 2000-qubit QA hardware was still unable to outperform SA for the MU problem. Here, we propose the use of quantum-inspired algorithm to solve the MU problem, in order to go beyond traditional SA. We introduce a highly compact phase encoding method which can exponentially reduce the representation space, compared with the previous one-hot encoding method. For benchmarking, we tested this new approach on the public QM9 dataset generated by density functional theory (DFT). The root-mean-square deviation between the conformation determined by our approach and DFT is negligible (less than about 0.5Å), which underpins the validity of our approach. Furthermore, the median time-to-target metric can be reduced by a factor of five compared to SA. Additionally, we demonstrate a simulation experiment by MindQuantum using quantum approximate optimization algorithm (QAOA) to reach optimal results. These results indicate that quantum-inspired algorithms can be applied to solve practical problems even before quantum hardware becomes mature.

METHODS

The objective function of MU is defined as the sum of all internal distances between atoms in the molecule, which is a high-order unconstrained binary optimization (HUBO) problem. The degree of freedom of variables is discretized and encoded with binary variables by the phase encoding method. We employ the quantum-inspired simulated bifurcation algorithm for optimization. The public QM9 dataset is generated by DFT. The simulation experiment of quantum computation is implemented by MindQuantum using QAOA.

Ensemble-Based Virtual Screening Led to the Discovery of Novel Lead Molecules as Potential NMBAs

Neuromuscular blocking agents (NMBAs) are routinely used during anesthesia to relax skeletal muscle. Nicotinic acetylcholine receptors (nAChRs) are ligand-gated ion channels; NMBAs can induce muscle paralysis by preventing the neurotransmitter acetylcholine (ACh) from binding to nAChRs situated on the postsynaptic membranes. Despite widespread efforts, it is still a great challenge to find new NMBAs since the introduction of cisatracurium in 1995. In this work, an effective ensemble-based virtual screening method, including molecular property filters, 3D pharmacophore model, and molecular docking, was applied to discover potential NMBAs from the ZINC15 database. The results showed that screened hit compounds had better docking scores than the reference compound d-tubocurarine. In order to further investigate the binding modes between the hit compounds and nAChRs at simulated physiological conditions, the molecular dynamics simulation was performed. Deep analysis of the simulation results revealed that ZINC257459695 can stably bind to nAChRs' active sites and interact with the key residue Asp165. The binding free energies were also calculated for the obtained hits using the MM/GBSA method. In silico ADMET calculations were performed to assess the pharmacokinetic properties of hit compounds in the human body. Overall, the identified ZINC257459695 may be a promising lead compound for developing new NMBAs as an adjunct to general anesthesia, necessitating further investigations.

Hit discovery of potential CDK8 inhibitors and analysis of amino acid mutations for cancer therapy through computer-aided drug discovery

Cyclin-dependent kinase 8 (CDK8) has emerged as a promising target for inhibiting cancer cell function, intensifying efforts towards the development of CDK8 inhibitors as potential cancer therapeutics. Mutations in CDK8, a protein kinase, are also implicated as a primary factor associated with tumor formation. In this study, we identified potential inhibitors through virtual screening for CDK8 and single amino acid mutations in CDK8, namely D173A (Aspartate 173 mutate to Alanine), D189N (Aspartate 189 mutate to Asparagine), T196A (Threonine 196 mutate to Alanine) and T196D (Threonine 196 mutate to Aspartate). Four databases (CHEMBEL, ZINC, MCULE, and MolPort) containing 65,209,131 molecules have been searched to identify new inhibitors for CDK8 and its single mutations. In the first step, structure-based pharmacophore modeling in the Pharmit server was used to select the compounds to know the inhibitors. Then molecules with better predicted drug-like molecule properties were selected. The final filter used to select more effective inhibitors among the previously selected molecules was molecular docking. Finally, 13 hits for CDK8, 11 hits for D173A, 11 hits for D189N, 15 hits for T196A, and 12 hits for T196D were considered potential inhibitors. A majority of the virtual screening hits exhibited satisfactorily predict pharmacokinetic characteristics and toxicity properties.

Orphan GPR52 as an emerging neurotherapeutic target

GPR52 is a highly conserved, brain-enriched, Gs/olf-coupled orphan G protein-coupled receptor (GPCR) that controls various cyclic AMP (cAMP)-dependent physiological and pathological processes. Stimulation of GPR52 activity might be beneficial for the treatment of schizophrenia, psychiatric disorders and other human neurological diseases, whereas inhibition of its activity might provide a potential therapeutic approach for Huntington's disease. Excitingly, HTL0048149 (HTL'149), an orally available GPR52 agonist, has been advanced into phase I human clinical trials for the treatment of schizophrenia. In this concise review, we summarize the current understanding of GPR52 receptor distribution as well as its structure and functions, highlighting the recent advances in drug discovery efforts towards small-molecule GPR52 ligands. The opportunities and challenges presented by targeting GPR52 for novel therapeutics are also briefly discussed.

Molecular Gas-Phase Conformational Ensembles

Accurately determining the global minima of a molecular structure is important in diverse scientific fields, including drug design, materials science, and chemical synthesis. Conformational search engines serve as valuable tools for exploring the extensive conformational space of molecules and for identifying energetically favorable conformations. In this study, we present a comparison of Auto3D, CREST, Balloon, and ETKDG (from RDKit), which are freely available conformational search engines, to evaluate their effectiveness in locating global minima. These engines employ distinct methodologies, including machine learning (ML) potential-based, semiempirical, and force field-based approaches. To validate these methods, we propose the use of collisional cross-section (CCS) values obtained from ion mobility-mass spectrometry studies. We hypothesize that experimental gas-phase CCS values can provide experimental evidence that we likely have the global minimum for a given molecule. To facilitate this effort, we used our gas-phase conformation library (GPCL) which currently consists of the full ensembles of 20 small molecules and can be used by the community to validate any conformational search engine. Further members of the GPCL can be readily created for any molecule of interest using our standard workflow used to compute CCS values, expanding the ability of the GPCL in validation exercises. These innovative validation techniques enhance our understanding of the conformational landscape and provide valuable insights into the performance of conformational generation engines. Our findings shed light on the strengths and limitations of each search engine, enabling informed decisions for their utilization in various scientific fields, where accurate molecular structure determination is crucial for understanding biological activity and designing targeted interventions. By facilitating the identification of reliable conformations, this study significantly contributes to enhancing the efficiency and accuracy of molecular structure determination, with particular focus on metabolite structure elucidation. The findings of this research also provide valuable insights for developing effective workflows for predicting the structures of unknown compounds with high precision.

Integrated virtual screening and molecular dynamics simulation approaches revealed potential natural inhibitors for DNMT1 as therapeutic solution for triple negative breast cancer

Triple negative breast cancers (TNBC) are clinically heterogeneous but mostly aggressive malignancies devoid of expression of the estrogen, progesterone, and HER2 (ERBB2 or NEU) receptors. It accounts for 15-20% of all cases. Altered epigenetic regulation including DNA hypermethylation by DNA methyltransferase 1 (DNMT1) has been implicated as one of the causes of TNBC tumorigenesis. The antitumor effect of DNMT1 has also been explored in TNBC that currently lacks targeted therapies. However, the actual treatment for TNBC is yet to be discovered. This study is attributed to the identification of novel drug targets against TNBC. A comprehensive docking and simulation analysis was performed to optimize promising new compounds by estimating their binding affinity to the target protein. Molecular dynamics simulation of 500 ns well complemented the binding affinity of the compound and revealed strong stability of predicted compounds at the docked site. Calculation of binding free energies using MMPBSA and MMGBSA validated the strong binding affinity between compound and binding pockets of DNMT1. In a nutshell, our study uncovered that Beta-Mangostin, Gancaonin Z, 5-hydroxysophoranone, Sophoraflavanone L, and Dorsmanin H showed maximum binding affinity with the active sites of DNMT1. Furthermore, all of these compounds depict maximum drug-like properties. Therefore, the proposed compounds can be a potential candidate for patients with TNBC, but, experimental validation is needed to ensure their safety.Communicated by Ramaswamy H. Sarma.

Curation and cheminformatics analysis of a Ugi-reaction derived library (URDL) of synthetically tractable small molecules for virtual screening application

Virtual screening (VS) is an important approach in drug discovery and relies on the availability of a virtual library of synthetically tractable molecules. Ugi reaction (UR) represents an important multi-component reaction (MCR) that reliably produces a peptidomimetic scaffold. Recent literature shows that a tactically assembled Ugi adduct can be subjected to further chemical modifications to yield a variety of rings and scaffolds, thus, renewing the interest in this old reaction. Given the reliability and efficiency of UR, we collated an UR derived library (URDL) of small molecules (total = 5773) for VS. The synthesis of the majority of URDL molecules may be carried out in 1-2 pots in a time and cost-effective manner. The detailed analysis of the average property and chemical space of URDL was also carried out using the open-source Datawarrior program. The comparison with FDA-approved oral drugs and inhibitors of protein-protein interactions (iPPIs) suggests URDL molecules are 'clean', drug-like, and conform to a structurally distinct space from the other two categories. The average physicochemical properties of compounds in the URDL library lie closer to iPPI molecules than oral drugs thus suggesting that the URDL resource can be applied to discover novel iPPI molecules. The URDL molecules consist of diverse ring systems, many of which have not been exploited yet for drug design. Thus, URDL represents a small virtual library of drug-like molecules with unexplored chemical space designed for VS. The structures of all molecules of URDL, oral drugs, and iPPI compounds are being made freely accessible as supplementary information for broader application.

Computational-guided approach for identification of PI3K alpha inhibitor in the treatment of hepatocellular carcinoma by virtual screening and water map analysis

Hepatocellular carcinoma (HCC) is one of the most deadly disorders, with a relative survival rate of 36% in the last 5 years. After an extensive literature survey and pathophysiology analysis, PI3Kα was found to be a promising biological target as PIK3CA gene upregulation was observed in HCC, resulting in the loss of apoptosis of cells, which leads to uncontrollable growth and proliferation. Due to superior selectivity and promising therapeutic activity, the PI3K-targeted molecule library was selected, and the ligand preparation was executed. The study mainly focused on e-pharmacophore development, virtual screening and receptor-ligand docking analysis. Then, MMGBSA and ADME prediction analysis was performed with the top 10 molecules; for further analysis of ligand-receptor binding affinity at the catalytic binding site, induced fit docking was performed with the top two molecules. The analysis of quantum chemical stability descriptors, i.e., frontier molecular orbital analysis, was performed followed by molecular dynamics simulation of 100 ns to better understand the ligand-receptor binding. In this study, water map analysis played a significant role in the hit optimization and analysis of the thermodynamic properties of the receptor-ligand complex. The two hit molecules K894-1435 and K894-1045 represented superior docking scores, enhanced stability, and inhibitory action targeting Valine 851 amino acid residue at the catalytic binding site. Hence, the study has significance for the quest for selective PI3Kα inhibitors through the process of hit-to-lead optimization.Communicated by Ramaswamy H. Sarma.

Network-Based Methods and Their Applications in Drug Discovery

Drug discovery is time-consuming, expensive, and predominantly follows the "one drug → one target → one disease" paradigm. With the rapid development of systems biology and network pharmacology, a novel drug discovery paradigm, "multidrug → multitarget → multidisease", has emerged. This new holistic paradigm of drug discovery aligns well with the essence of networks, leading to the emergence of network-based methods in the field of drug discovery. In this Perspective, we initially introduce the concept and data sources of networks and highlight classical methodologies employed in network-based methods. Subsequently, we focus on the practical applications of network-based methods across various areas of drug discovery, such as target prediction, virtual screening, prediction of drug therapeutic effects or adverse drug events, and elucidation of molecular mechanisms. In addition, we provide representative web servers for researchers to use network-based methods in specific applications. Finally, we discuss several challenges of network-based methods and the directions for future development. In a word, network-based methods could serve as powerful tools to accelerate drug discovery.

Structure-based virtual screening and molecular docking approaches to identify potential inhibitors against KIF2C to combat glioma

Glioma, a kind of malignant brain tumor, is extremely lethal. Kinesin family member 2C (KIF2C) was found to have an aberrant expression in several cancer types, including lung cancer and glioma. KIF2C may therefore be a useful therapeutic target for the treatment of glioma. In the current study, new drug candidates that may function as KIF2C enzyme inhibitors were discovered. MTi OpenScreen was used to carry out the structure-based virtual screening of an inbuilt drug library containing 150,000 compounds. These compounds belong to different classes, such as natural product-based compounds (NP-lib), purchasable approved drugs (Drugs-lib), and food constituents compound collection (FOOD-lib). Based on their binding affinities, a total of 84 compounds were further pushed to calculate ADMET properties. The compounds (16) meeting the ADMET cutoff ranges were then further docked to the receptor to find their plausible binding modes using the Glide tool's standard precision (SP) technique. The docking results were examined using the Glide gscore, and the best binding compounds (Rimacalib and Sarizotan) were chosen to test their stability with KIF2C protein through molecular dynamics (MD) simulation. Similarly, Principal Component Analysis and cross-correlation matrix were also examined. The MM/GBSA binding free energies showed a considerable energy contribution in the binding of hits with the KIF2C. Collectively, these findings strongly suggest the potential of the lead compounds to inhibit the biological function of KIF2C, emphasizing the need for further investigation in this area.Communicated by Ramaswamy H. Sarma.

Exploring natural compound, Panicutine as leucine-rich repeat kinase 2 inhibitor against Parkinson's disease: a structure-guided approach

Leucine-rich repeat kinase 2 (LRRK2) is a promising drug target for the therapeutic management of Parkinson's disease (PD) and other neurodegenerative disorders. LRRK2 inhibitors have the potential to modulate neuroinflammation, reduce alpha-synuclein aggregation and improve motor symptoms in PD patients. Although LRRK2 inhibitors are still in the early stages of clinical development, the identification of potent and selective inhibitors through structure-guided approaches provides a promising avenue for the development of effective therapies for PD and other neurodegenerative disorders. In this study, natural compounds from the IMPPAT database were screened using a state-of-the-art computational virtual screening approach to identify potential inhibitors of LRRK2. We carried out a docking screening on a library of natural compounds and identified a few compounds with strong binding affinity, docking score and specificity towards LRRK2 as the top hits. These hits were then subjected to further analysis based on multiple parameters for the Pan-assay interference compounds and their physicochemical and pharmacokinetics evaluation followed by a detailed interaction analysis. After careful evaluation, one natural compound, Panicutine, was identified as a promising candidate for LRRK2 due to its significant affinity and specificity towards the LRRK2 binding pocket. Additionally, it exhibited drug-like properties with blood-brain barrier permeability as determined by ADMET properties. To gain a deeper understanding of the stability and conformational changes of the LRRK2-ligand complex, MD simulations were conducted for 100 nanoseconds under explicit solvent conditions followed by principal component analysis and free energy dynamics. The simulation results demonstrated that the LRRK2-Panicutine complex remained stable throughout the simulation trajectories. Based on these findings, it is concluded that Panicutine has the potential to act as a LRRK2 inhibitor against PD and other neurodegenerative disorders.Communicated by Ramaswamy H. Sarma.

High-Quality Conformer Generation with CONFORGE: Algorithm and Performance Assessment

Knowledge of the putative bound-state conformation of a molecule is an essential prerequisite for the successful application of many computer-aided drug design methods that aim to assess or predict its capability to bind to a particular target receptor. An established approach to predict bioactive conformers in the absence of receptor structure information is to sample the low-energy conformational space of the investigated molecules and derive representative conformer ensembles that can be expected to comprise members closely resembling possible bound-state ligand conformations. The high relevance of such conformer generation functionality led to the development of a wide panel of dedicated commercial and open-source software tools throughout the last decades. Several published benchmarking studies have shown that open-source tools usually lag behind their commercial competitors in many key aspects. In this work, we introduce the open-source conformer ensemble generator CONFORGE, which aims at delivering state-of-the-art performance for all types of organic molecules in drug-like chemical space. The ability of CONFORGE and several well-known commercial and open-source conformer ensemble generators to reproduce experimental 3D structures as well as their computational efficiency and robustness has been assessed thoroughly for both typical drug-like molecules and macrocyclic structures. For small molecules, CONFORGE clearly outperformed all other tested open-source conformer generators and performed at least equally well as the evaluated commercial generators in terms of both processing speed and accuracy. In the case of macrocyclic structures, CONFORGE achieved the best average accuracy among all benchmarked generators, with RDKit's generator coming close in second place.

Green Synthesis of Zinc Oxide Nanoparticles Using Cymbopogon citratus Extract and Its Antibacterial Activity

Excessive use of antimicrobial medications including antibiotics has led to the emerging menace of antimicrobial resistance, which, as per the World Health Organization (WHO), is among the top ten public health threats facing humanity, globally. This necessitates that innovative technologies be sought that can aid in the elimination of pathogens and hamper the spread of infections. Zinc oxide (ZnO) has multifunctionality owing to its extraordinary physico-chemical properties and functionality in a range of applications. In this research, ZnO nanoparticles (NPs) were synthesized from zinc nitrate hexahydrate, by a green synthesis approach using Cymbopogon citratus extract followed by characterization of the NPs. The obtained X-ray diffraction peaks of ZnO NPs matched with the standard JCPDS card (no. 89-510). The particles had a size of 20-24 nm, a wurtzite structure with a high crystallinity, and hexagonal rod-like shape. UV-Vis spectroscopy revealed absorption peaks between 369 and 374 nm of ZnO NPs synthesized from C. citratus extract confirming the formation of ZnO. Fourier transform infrared confirmed the ZnO NPs as strong absorption bands were observed in the range of 381-403 cm-1 corresponding to Zn-O bond stretching. Negative values of the highest occupied molecular orbital-lowest unoccupied molecular orbital for ZnO NPs indicated the good potential to form a stable ligand-protein complex. Docking results indicated favorable binding interaction between ZnO and DNA gyrase subunit b with a binding energy of -2.93 kcal/mol. ZnO NPs at various concentrations inhibited the growth of Escherichia coli and Staphylococcus aureus. Minimum inhibitory concentration values of ZnO NPs against E. coli and S. aureus were found to be 92.07 ± 0.13 and 88.13 ± 0.35 μg/mL, respectively, at a concentration of 2 mg/mL. AO/EB staining and fluorescence microscopy revealed the ability of ZnO NPs to kill E. coli and S. aureus cells. Through the findings of this study, it has been shown that C. citratus extract can be used in a green synthesis approach to generate ZnO NPs, which can be employed as alternatives to antibiotics and a tool to eliminate drug-resistant microbes in the future.

A Review on Artificial Intelligence Approaches and Rational Approaches in Drug Discovery

Artificial intelligence (AI) speeds up the drug development process and reduces its time, as well as the cost which is of enormous importance in outbreaks such as COVID-19. It uses a set of machine learning algorithms that collects the available data from resources, categorises, processes and develops novel learning methodologies. Virtual screening is a successful application of AI, which is used in screening huge drug-like databases and filtering to a small number of compounds. The brain's thinking of AI is its neural networking which uses techniques such as Convoluted Neural Network (CNN), Recursive Neural Network (RNN) or Generative Adversial Neural Network (GANN). The application ranges from small molecule drug discovery to the development of vaccines. In the present review article, we discussed various techniques of drug design, structure and ligand-based, pharmacokinetics and toxicity prediction using AI. The rapid phase of discovery is the need of the hour and AI is a targeted approach to achieve this.

The Impact of Supervised Learning Methods in Ultralarge High-Throughput Docking

Structure-based virtual screening methods are, nowadays, one of the key pillars of computational drug discovery. In recent years, a series of studies have reported docking-based virtual screening campaigns of large databases ranging from hundreds to thousands of millions compounds, further identifying novel hits after experimental validation. As these larg-scale efforts are not generally accessible, machine learning-based protocols have emerged to accelerate the identification of virtual hits within an ultralarge chemical space, reaching impressive reductions in computational time. Herein, we illustrate the motivation and the problem behind the screening of large databases, providing an overview of key concepts and essential applications of machine learning-accelerated protocols, specifically concerning supervised learning methods. We also discuss where the field stands with these novel developments, highlighting possible insights for future studies.

Design and Identification of a Novel Antiviral Affinity Peptide against Fowl Adenovirus Serotype 4 (FAdV-4) by Targeting Fiber2 Protein

Outbreaks of hydropericardium hepatitis syndrome caused by fowl adenovirus serotype 4 (FAdV-4) with a novel genotype have been reported in China since 2015, with significant economic losses to the poultry industry. Fiber2 is one of the important structural proteins on FAdV-4 virions. In this study, the C-terminal knob domain of the FAdV-4 Fiber2 protein was expressed and purified, and its trimer structure (PDB ID: 7W83) was determined for the first time. A series of affinity peptides targeting the knob domain of the Fiber2 protein were designed and synthesized on the basis of the crystal structure using computer virtual screening technology. A total of eight peptides were screened using an immunoperoxidase monolayer assay and RT-qPCR, and they exhibited strong binding affinities to the knob domain of the FAdV-4 Fiber2 protein in a surface plasmon resonance assay. Treatment with peptide number 15 (P15; WWHEKE) at different concentrations (10, 25, and 50 μM) significantly reduced the expression level of the Fiber2 protein and the viral titer during FAdV-4 infection. P15 was found to be an optimal peptide with antiviral activity against FAdV-4 in vitro with no cytotoxic effect on LMH cells up to 200 μM. This study led to the identification of a class of affinity peptides designed using computer virtual screening technology that targeted the knob domain of the FAdV-4 Fiber2 protein and may be developed as a novel potential and effective antiviral strategy in the prevention and control of FAdV-4.

Identification of Potent Small-Molecule PCSK9 Inhibitors Based on Quantitative Structure-Activity Relationship, Pharmacophore Modeling, and Molecular Docking Procedure

The leading cause of atherosclerotic cardiovascular disease (ASCVD) is elevated low-density lipoprotein cholesterol (LDL-C). Proprotein convertase subtilisin/kexin type 9 (PCSK9) attaches to the domain of LDL receptor (LDLR), diminishing LDL-C influx and LDLR cell surface presentation in hepatocytes, resulting in higher circulating LDL-C levels. PCSK9 dysfunction has been linked to lower levels of plasma LDLC and a decreased risk of coronary heart disease (CHD). Herein, using virtual screening tools, we aimed to identify a potent small-molecule PCSK9 inhibitor in compounds that are currently being studied in clinical trials. We first performed chemical absorption, distribution, metabolism, excretion, and toxicity (ADMET) filtering of 9800 clinical trial compounds obtained from the ZINC 15 database using Lipinski's rule of 5 and achieved 3853 compounds. Two-dimensional (2D) quantitative structure-activity relationship (QSAR) was initiated by computing molecular descriptors and selecting important descriptors of 23 PCSK9 inhibitors. Multivariate calibration was performed with the partial least square regression (PLS) method with 18 compounds for training to design the QSAR model and 5 compounds for the test set to assess the model. The best latent variables (LV) (LV=6) with the lowest value of Root-Mean-Square Error of Cross-Validation (RMSECV) of 0.48 and leave-one-out cross-validation correlation coefficient (R²CV) = 0.83 were obtained for the QSAR model. The low RMSEC (0.21) with high R²cal (0.966) indicates the probability of fit between the experimental data and the calibration model. Using QSAR analysis of 3853 compounds, 2635 had a pIC₅₀<1 and were considered for pharmacophore screening. The PHASE module (a complete package for pharmacophore modeling) designed the pharmacophore hypothesis through multiple ligands. The top 14 compounds (pIC₅₀>1) were defined as active, whereas 9 (pIC₅₀<1) were considered as an inactive set. Three five-point pharmacophore hypotheses achieved the highest score: DHHRR1, DHHRR2, and DHRRR1. The highest and best model with survival scores (5.365) was DHHRR1, comprising 1 hydrogen donor (D), 2 hydrophobic groups (H), and 2 rings of aromatic (R) features. We selected the molecules with a higher 1.5 fitness score (257 compounds) in pharmacophore screening (DHHRR1) for molecular docking screening. Molecular docking indicates that ZINC000051951669, with a binding affinity: of -13.2 kcal/mol and 2 H-bonds, has the highest binding to the PCSK9 protein. ZINC000011726230 with energy binding: -11.4 kcal/mol and 3 H-bonds, ZINC000068248147 with binding affinity: -10.7 kcal/mol and 1 H-bond, ZINC000029134440 with a binding affinity: -10.6 kcal/mol and 4 H-bonds were ranked next, respectively. To conclude, the archived molecules identified as inhibitory PCSK9 candidates, and especially ZINC000051951669 may therefore significantly inhibit PCSK9 and should be considered in the newly designed trials.

Application of computational methods for class A GPCR Ligand discovery

G protein-coupled receptors (GPCR) are integral membrane proteins of considerable interest as targets for drug development due to their role in transmitting cellular signals in a multitude of biological processes. Of the six classes categorizing GPCR (A, B, C, D, E, and F), class A contains the largest number of therapeutically relevant GPCR. Despite their importance as drug targets, many challenges exist for the discovery of novel class A GPCR ligands serving as drug precursors. Though knowledge of the structural and functional characteristics of GPCR has grown significantly over the past 20 years, a large portion of GPCR lack reported, experimentally determined structures. Furthermore, many GPCR have no known endogenous and/or synthetic ligands, limiting further exploration of their biochemical, cellular, and physiological roles. While many successes in GPCR ligand discovery have resulted from experimental high-throughput screening, computational methods have played an increasingly important role in GPCR ligand identification in the past decade. Here we discuss computational techniques applied to GPCR ligand discovery. This review summarizes class A GPCR structure/function and provides an overview of many obstacles currently faced in GPCR ligand discovery. Furthermore, we discuss applications and recent successes of computational techniques used to predict GPCR structure as well as present a summary of ligand- and structure-based methods used to identify potential GPCR ligands. Finally, we discuss computational hit list generation and refinement and provide comprehensive workflows for GPCR ligand identification.

Rational design of stapled antimicrobial peptides

The global increase in antimicrobial drug resistance has dramatically reduced the effectiveness of traditional antibiotics. Structurally diverse antibiotics are urgently needed to combat multiple-resistant bacterial infections. As part of innate immunity, antimicrobial peptides have been recognized as the most promising candidates because they comprise diverse sequences and mechanisms of action and have a relatively low induction rate of resistance. However, because of their low chemical stability, susceptibility to proteases, and high hemolytic effect, their usage is subject to many restrictions. Chemical modifications such as D-amino acid substitution, cyclization, and unnatural amino acid modification have been used to improve the stability of antimicrobial peptides for decades. Among them, a side-chain covalent bridge modification, the so-called stapled peptide, has attracted much attention. The stapled side-chain bridge stabilizes the secondary structure, induces protease resistance, and increases cell penetration and biological activity. Recent progress in computer-aided drug design and artificial intelligence methods has also been used in the design of stapled antimicrobial peptides and has led to the successful discovery of many prospective peptides. This article reviews the possible structure-activity relationships of stapled antimicrobial peptides, the physicochemical properties that influence their activity (such as net charge, hydrophobicity, helicity, and dipole moment), and computer-aided methods of stapled peptide design. Antimicrobial peptides under clinical trial: Pexiganan (NCT01594762, 2012-05-07). Omiganan (NCT02576847, 2015-10-13).

CADD, AI and ML in drug discovery: A comprehensive review

Computer-aided drug design (CADD) is an emerging field that has drawn a lot of interest because of its potential to expedite and lower the cost of the drug development process. Drug discovery research is expensive and time-consuming, and it frequently took 10-15 years for a drug to be commercially available. CADD has significantly impacted this area of research. Further, the combination of CADD with Artificial Intelligence (AI), Machine Learning (ML), and Deep Learning (DL) technologies to handle enormous amounts of biological data has reduced the time and cost associated with the drug development process. This review will discuss how CADD, AI, ML, and DL approaches help identify drug candidates and various other steps of the drug discovery process. It will also provide a detailed overview of the different in silico tools used and how these approaches interact.

Improving drug discovery through parallelism

Compound identification in ligand-based virtual screening is limited by two key issues: the quality and the time needed to obtain predictions. In this sense, we designed OptiPharm, an algorithm that obtained excellent results in improving the sequential methods in the literature. In this work, we go a step further and propose its parallelization. Specifically, we propose a two-layer parallelization. Firstly, an automation of the molecule distribution process between the available nodes in a cluster, and secondly, a parallelization of the internal methods (initialization, reproduction, selection and optimization). This new software, called pOptiPharm, aims to improve the quality of predictions and reduce experimentation time. As the results show, the performance of the proposed methods is good. It can find better solutions than the sequential OptiPharm, all while reducing its computation time almost proportionally to the number of processing units considered.

Virtual Screening-Based Drug Development for the Treatment of Nervous System Diseases

The incidence rate of nervous system diseases has increased in recent years. Nerve injury or neurodegenerative diseases usually cause neuronal loss and neuronal circuit damage, which seriously affect motor nerve and autonomic nervous function. Therefore, safe and effective treatment is needed. As traditional drug research becomes slower and more expensive, it is vital to enlist the help of cutting- edge technology. Virtual screening (VS) is an attractive option for the identification and development of promising new compounds with high efficiency and low cost. With the assistance of computer- aided drug design (CADD), VS is becoming more and more popular in new drug development and research. In recent years, it has become a reality to transform non-neuronal cells into functional neurons through small molecular compounds, which provides a broader application prospect than transcription factor-mediated neuronal reprogramming. This review mainly summarizes related theory and technology of VS and the drug research and development using VS technology in nervous system diseases in recent years, and focuses more on the potential application of VS technology in neuronal reprogramming, thus facilitating new drug design for both prevention and treatment of nervous system diseases.

Identifying the natural compound Catechin from tropical mangrove plants as a potential lead candidate against 3CLpro from SARS-CoV-2: An integrated in silico approach

SARS-CoV-2, a member of beta coronaviruses, is a single-stranded, positive-sense RNA virus responsible for the COVID-19 pandemic. With global fatalities of the pandemic exceeding 4.57 million, it becomes crucial to identify effective therapeutics against the virus. A protease, 3CL^pro, is responsible for the proteolysis of viral polypeptides into functional proteins, which is essential for viral pathogenesis. This indispensable activity of 3CL^pro makes it an attractive target for inhibition studies. The current study aimed to identify potential lead molecules against 3CL^pro of SARS-CoV-2 using a manually curated in-house library of antiviral compounds from mangrove plants. This study employed the structure-based virtual screening technique to evaluate an in-house library of antiviral compounds against 3CL^pro of SARS-CoV-2. The library was comprised of thirty-three experimentally proven antiviral molecules extracted from different species of tropical mangrove plants. The molecules in the library were virtually screened using AutoDock Vina, and subsequently, the top five promising 3CL^pro-ligand complexes along with 3CL^pro-N3 (control molecule) complex were subjected to MD simulations to comprehend their dynamic behaviour and structural stabilities. Finally, the MM/PBSA approach was used to calculate the binding free energies of 3CLpro complexes. Among all the studied compounds, Catechin achieved the most significant binding free energy (-40.3 ± 3.1 kcal/mol), and was closest to the control molecule (-42.8 ± 5.1 kcal/mol), and its complex with 3CL^pro exhibited the highest structural stability. Through extensive computational investigations, we propose Catechin as a potential therapeutic agent against SARS-CoV-2. Communicated by Ramaswamy H. Sarma.

A Hybrid Docking and Machine Learning Approach to Enhance the Performance of Virtual Screening Carried out on Protein-Protein Interfaces

The modulation of protein-protein interactions (PPIs) by small chemical compounds is challenging. PPIs play a critical role in most cellular processes and are involved in numerous disease pathways. As such, novel strategies that assist the design of PPI inhibitors are of major importance. We previously reported that the knowledge-based DLIGAND2 scoring tool was the best-rescoring function for improving receptor-based virtual screening (VS) performed with the Surflex docking engine applied to several PPI targets with experimentally known active and inactive compounds. Here, we extend our investigation by assessing the vs. potential of other types of scoring functions with an emphasis on docking-pose derived solvent accessible surface area (SASA) descriptors, with or without the use of machine learning (ML) classifiers. First, we explored rescoring strategies of Surflex-generated docking poses with five GOLD scoring functions (GoldScore, ChemScore, ASP, ChemPLP, ChemScore with Receptor Depth Scaling) and with consensus scoring. The top-ranked poses were post-processed to derive a set of protein and ligand SASA descriptors in the bound and unbound states, which were combined to derive descriptors of the docked protein-ligand complexes. Further, eight ML models (tree, bagged forest, random forest, Bayesian, support vector machine, logistic regression, neural network, and neural network with bagging) were trained using the derivatized SASA descriptors and validated on test sets. The results show that many SASA descriptors are better than Surflex and GOLD scoring functions in terms of overall performance and early recovery success on the used dataset. The ML models were superior to all scoring functions and rescoring approaches for most targets yielding up to a seven-fold increase in enrichment factors at 1% of the screened collections. In particular, the neural networks and random forest-based ML emerged as the best techniques for this PPI dataset, making them robust and attractive vs. tools for hit-finding efforts. The presented results suggest that exploring further docking-pose derived SASA descriptors could be valuable for structure-based virtual screening projects, and in the present case, to assist the rational design of small-molecule PPI inhibitors.

Application of Computational Biology and Artificial Intelligence in Drug Design

Traditional drug design requires a great amount of research time and developmental expense. Booming computational approaches, including computational biology, computer-aided drug design, and artificial intelligence, have the potential to expedite the efficiency of drug discovery by minimizing the time and financial cost. In recent years, computational approaches are being widely used to improve the efficacy and effectiveness of drug discovery and pipeline, leading to the approval of plenty of new drugs for marketing. The present review emphasizes on the applications of these indispensable computational approaches in aiding target identification, lead discovery, and lead optimization. Some challenges of using these approaches for drug design are also discussed. Moreover, we propose a methodology for integrating various computational techniques into new drug discovery and design.

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

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

Identification of a new and diverse set of Mycobacterium tuberculosstais uracil-DNA glycosylase (MtUng) inhibitors using structure-based virtual screening: experimental validation and molecular dynamics studies

Mycobacterium tuberculosis uracil-DNA glycosylase (MtUng), a key DNA repair enzyme, represents an attractive target for the design of new antimycobacterial agents. However, only a limited number of weak MtUng inhibitors are reported, primarily based on the uracil ring, and hence, lack diversity. We report the first structure-based virtual screening (SBVS) using three separate libraries consisting of uracil and non-uracil small molecules, together with the FDA-approved drugs. Twenty diverse virtual hits with the highest predicted binding were procured and screened using a fluorescence-based assay to evaluate their potential to inhibit MtUng. Several of these molecules were found to inhibit MtUng activity at low mM and µM levels, comparable to or better than several other reported Ung inhibitors. Thus, these molecules represent a diverse set of scaffolds for developing next-generation MtUng inhibitors. The most active uracil-based compound 5 (IC₅₀ = 0.14 mM) was found to be ∼15-fold more potent than the positive control, uracil. The binding stability and conformation of compound 5 in complex with the enzyme were further confirmed using molecular dynamics simulation.

DENVIS: Scalable and High-Throughput Virtual Screening Using Graph Neural Networks with Atomic and Surface Protein Pocket Features

Computational methods for virtual screening can dramatically accelerate early-stage drug discovery by identifying potential hits for a specified target. Docking algorithms traditionally use physics-based simulations to address this challenge by estimating the binding orientation of a query protein-ligand pair and a corresponding binding affinity score. Over the recent years, classical and modern machine learning architectures have shown potential for outperforming traditional docking algorithms. Nevertheless, most learning-based algorithms still rely on the availability of the protein-ligand complex binding pose, typically estimated via docking simulations, which leads to a severe slowdown of the overall virtual screening process. A family of algorithms processing target information at the amino acid sequence level avoid this requirement, however, at the cost of processing protein data at a higher representation level. We introduce deep neural virtual screening (DENVIS), an end-to-end pipeline for virtual screening using graph neural networks (GNNs). By performing experiments on two benchmark databases, we show that our method performs competitively to several docking-based, machine learning-based, and hybrid docking/machine learning-based algorithms. By avoiding the intermediate docking step, DENVIS exhibits several orders of magnitude faster screening times (i.e., higher throughput) than both docking-based and hybrid models. When compared to an amino acid sequence-based machine learning model with comparable screening times, DENVIS achieves dramatically better performance. Some key elements of our approach include protein pocket modeling using a combination of atomic and surface features, the use of model ensembles, and data augmentation via artificial negative sampling during model training. In summary, DENVIS achieves competitive to state-of-the-art virtual screening performance, while offering the potential to scale to billions of molecules using minimal computational resources.

MO-MEMES: A method for accelerating virtual screening using multi-objective Bayesian optimization

The pursuit of potential inhibitors for novel targets has become a very important problem especially over the last 2 years with the world in the midst of the COVID-19 pandemic. This entails performing high throughput screening exercises on drug libraries to identify potential “hits”. These hits are identified using analysis of their physical properties like binding affinity to the target receptor, octanol-water partition coefficient (LogP) and more. However, drug libraries can be extremely large and it is infeasible to calculate and analyze the physical properties for each of those molecules within acceptable time and moreover, each molecule must possess a multitude of properties apart from just the binding affinity. To address this problem, in this study, we propose an extension to the Machine learning framework for Enhanced MolEcular Screening (MEMES) framework for multi-objective Bayesian optimization. This approach is capable of identifying over 90% of the most desirable molecules with respect to all required properties while explicitly calculating the values of each of those properties on only 6% of the entire drug library. This framework would provide an immense boost in identifying potential hits that possess all properties required for a drug molecules.

Small Molecules as Toll-like Receptor 4 Modulators Drug and In-House Computational Repurposing

The innate immunity toll-like receptor 4 (TLR4) system is a receptor of paramount importance as a therapeutic target. Virtual screening following a “computer-aided drug repurposing” approach was applied to the discovery of novel TLR4 modulators with a non-lipopolysaccharide-like structure. We screened almost 29,000 approved drugs and drug-like molecules from commercial, public, and in-house academia chemical libraries and, after biological assays, identified several compounds with TLR4 antagonist activity. Our computational protocol showed to be a robust approach for the identification of hits with drug-like scaffolds as possible inhibitors of the TLR4 innate immune pathways. Our collaborative work broadens the chemical diversity for inspiration of new classes of TLR4 modulators.

Fragment-Based Drug Discovery against Mycobacteria: The Success and Challenges

The emergence of drug-resistant mycobacteria, including Mycobacterium tuberculosis (Mtb) and non-tuberculous mycobacteria (NTM), poses an increasing global threat that urgently demands the development of new potent anti-mycobacterial drugs. One of the approaches toward the identification of new drugs is fragment-based drug discovery (FBDD), which is the most ingenious among other drug discovery models, such as structure-based drug design (SBDD) and high-throughput screening. Specialized techniques, such as X-ray crystallography, nuclear magnetic resonance spectroscopy, and many others, are part of the drug discovery approach to combat the Mtb and NTM global menaces. Moreover, the primary drawbacks of traditional methods, such as the limited measurement of biomolecular toxicity and uncertain bioavailability evaluation, are successfully overcome by the FBDD approach. The current review focuses on the recognition of fragment-based drug discovery as a popular approach using virtual, computational, and biophysical methods to identify potent fragment molecules. FBDD focuses on designing optimal inhibitors against potential therapeutic targets of NTM and Mtb (PurC, ArgB, MmpL3, and TrmD). Additionally, we have elaborated on the challenges associated with the FBDD approach in the identification and development of novel compounds. Insights into the applications and overcoming the challenges of FBDD approaches will aid in the identification of potential therapeutic compounds to treat drug-sensitive and drug-resistant NTMs and Mtb infections.

Precise Assembly of Multiple Antigens on Nanoparticles with Specially Designed Affinity Peptides

Antigen proteins, assembled on nanoparticles, can be recognized by antigen-presenting cells effectively to enhance antigen immunogenicity. The ability to simultaneously display multiantigens on the same nanoparticle could have numerous applications but remained technical challenges. Here, we described a method for precise assembly of multiple antigens on nanoparticles with specially designed affinity peptides. First, we designed and screened affinity peptides with high affinity and specificity, which could respectively target the key amino acid residues of classical swine fever virus (CSFV) E2 protein or porcine circovirus type 2 capsid protein (PCV2 Cap) accurately. Then, we conjugated the antigen proteins to poly(lactic acid-glycolic acid) copolymer (PLGA) and Gram-positive enhancer matrix (GEM) nanoparticles through the peptides and perfectly assembled two kinds of multiantigen display nanoparticles with different particle sizes. Subsequently, the immunological properties of the assembled nanoparticles were tested. The results showed that the antigen display nanoparticles could promote the maturation, phagocytosis, and proinflammatory effects of antigen-presenting cells (APCs). Besides, compared with the antigen proteins, multiantigen display nanoparticles could induce much higher levels of antibodies and neutralizing antibodies in mice. This strategy may provide a technical support for the study of protein structure and the research and development of polyvalent vaccines.

Computational Methods in Cooperation with Experimental Approaches to Design Protein Tyrosine Phosphatase 1B Inhibitors in Type 2 Diabetes Drug Design: A Review of the Achievements of This Century

Protein tyrosine phosphatase 1B (PTP1B) dephosphorylates phosphotyrosine residues and is an important regulator of several signaling pathways, such as insulin, leptin, and the ErbB signaling network, among others. Therefore, this enzyme is considered an attractive target to design new drugs against type 2 diabetes, obesity, and cancer. To date, a wide variety of PTP1B inhibitors that have been developed by experimental and computational approaches. In this review, we summarize the achievements with respect to PTP1B inhibitors discovered by applying computer-assisted drug design methodologies (virtual screening, molecular docking, pharmacophore modeling, and quantitative structure–activity relationships (QSAR)) as the principal strategy, in cooperation with experimental approaches, covering articles published from the beginning of the century until the time this review was submitted, with a focus on studies conducted with the aim of discovering new drugs against type 2 diabetes. This review encourages the use of computational techniques and includes helpful information that increases the knowledge generated to date about PTP1B inhibition, with a positive impact on the route toward obtaining a new drug against type 2 diabetes with PTP1B as a molecular target.

New solid phase methodology for the synthesis of biscoumarin derivatives: experimental and in silico approaches

The simple and greener one-pot approach for the synthesis of biscoumarin derivatives using catalytic amounts of nano-MoO₃ catalyst under mortar-pestle grinding was described. The use of non-toxic and mild catalyst, cost-effectiveness, ordinary grinding, and good to the excellent yield of the final product makes this procedure a more attractive pathway for the synthesis of biologically remarkable pharmacophores. Accordingly, biscoumarin derivatives were successfully extended in the developed protocols. Next, a computational investigation was performed to identify the potential biological targets of this set of compounds. In this case, first, a similarity search on different virtual libraries was performed to find an ideal biological target for these derivatives. Results showed that the synthesized derivatives can be α-glucosidase inhibitors. In another step, molecular docking studies were carried out against human lysosomal acid-alpha-glucosidase (PDB ID: 5NN8) to determine the detailed binding modes and critical interactions with the proposed target. In silico assessments showed the gold score value in the range of 17.56 to 29.49. Additionally, molecular dynamic simulations and the MM-GBSA method of the most active derivative against α-glucosidase were conducted to study the behavior of selected compounds in the biological system. Ligand 1 stabilized after around 30 ns and participated in various interactions with Trp481, Asp518, Asp616, His674, Phe649, and Leu677 residues.

Simulation and Machine Learning Methods for Ion-Channel Structure Determination, Mechanistic Studies and Drug Design

Ion channels are expressed in almost all living cells, controlling the in-and-out communications, making them ideal drug targets, especially for central nervous system diseases. However, owing to their dynamic nature and the presence of a membrane environment, ion channels remain difficult targets for the past decades. Recent advancement in cryo-electron microscopy and computational methods has shed light on this issue. An explosion in high-resolution ion channel structures paved way for structure-based rational drug design and the state-of-the-art simulation and machine learning techniques dramatically improved the efficiency and effectiveness of computer-aided drug design. Here we present an overview of how simulation and machine learning-based methods fundamentally changed the ion channel-related drug design at different levels, as well as the emerging trends in the field.

Insights into the Pharmacological Effects of Flavonoids: The Systematic Review of Computer Modeling

Computer modeling is a method that is widely used in scientific investigations to predict the biological activity, toxicity, pharmacokinetics, and synthesis strategy of compounds based on the structure of the molecule. This work is a systematic review of articles performed in accordance with the recommendations of PRISMA and contains information on computer modeling of the interaction of classical flavonoids with different biological targets. The review of used computational approaches is presented. Furthermore, the affinities of flavonoids to different targets that are associated with the infection, cardiovascular, and oncological diseases are discussed. Additionally, the methodology of bias risks in molecular docking research based on principles of evidentiary medicine was suggested and discussed. Based on this data, the most active groups of flavonoids and lead compounds for different targets were determined. It was concluded that flavonoids are a promising object for drug development and further research of pharmacology by in vitro, ex vivo, and in vivo models is required.

Identifying SARS-CoV-2 main protease inhibitors by applying the computer screening of a large database of molecules

The outbreak of coronavirus disease 2019 (COVID-19) at the end of 2019 affected global health. Its infection agent was called severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Wearing a mask, maintaining social distance, and vaccination are effective ways to prevent infection of SARS-CoV-2, but none of them help infected people. Targeting the enzymes of SARS-CoV-2 is an effective way to stop the replication of the virus in infected people and treat COVID-19 patients. SARS-CoV-2 main protease is a therapeutic target which the inhibition of its enzymatic activity prevents from the replication of SARS-CoV-2. A large database of molecules has been searched to identify new inhibitors for SARS-CoV-2 main protease enzyme. At the first step, ligand screening based on similarity search was used to select similar compounds to known SARS-CoV-2 main protease inhibitors. Then molecules with better predicted pharmacokinetic properties were selected. Structure-based virtual screening based on the application of molecular docking and molecular dynamics simulation methods was used to select more effective inhibitors among selected molecules in previous step. Finally two compounds were considered as SARS-CoV-2 main protease inhibitors.

Development of New Inhibitors of HDAC1-3 Enzymes Aided by In Silico Design Strategies

Histone deacetylases (HDACs) are overexpressed in cancer, and their inhibition shows promising results in cancer therapy. In particular, selective class I HDAC inhibitors such as entinostat are proposed to be more beneficial in breast cancer treatment. Computational drug design is an inevitable part of today's drug discovery projects because of its unequivocal role in saving time and cost. Using three HDAC inhibitors trichostatin, vorinostat, and entinostat as template structures and a diverse fragment library, all synthetically accessible compounds thereof (∼3200) were generated virtually and filtered based on similarity against the templates and PAINS removal. The 298 selected structures were docked into the active site of HDAC I and ranked using a calculated binding affinity. Top-ranking structures were inspected manually, and, considering the ease of synthesis and drug-likeness, two new structures (3a and 3b) were proposed for synthesis and biological evaluation. The synthesized compounds were purified to a degree of more than 95% and structurally verified using various methods. The designed compounds 3a and 3b showed 65-80 and 5% inhibition on HDAC 1, 2, and 3 isoforms at a concentration of 10 μM, respectively. The novel compound 3a may be used as a lead structure for designing new HDAC inhibitors.

Highly Potent, Selective, Biostable, and Cell-Permeable Cyclic d-Peptide for Dual-Targeting Therapy of Lung Cancer

The application of peptide drugs in cancer therapy is impeded by their poor biostability and weak cell permeability. Therefore, it is imperative to find biostable and cell-permeable peptide drugs for cancer treatment. Here, we identified a potent, selective, biostable, and cell-permeable cyclic d-peptide, NKTP-3, that targets NRP1 and KRAS^G12D using structure-based virtual screening. NKTP-3 exhibited strong biostability and cellular uptake ability. Importantly, it significantly inhibited the growth of A427 cells with the KRAS^G12D mutation. Moreover, NKTP-3 showed strong antitumor activity against A427 cell-derived xenograft and KRAS^G12D-driven primary lung cancer models without obvious toxicity. This study demonstrates that the dual NRP1/KRAS^G12D-targeting cyclic d-peptide NKTP-3 may be used as a potential chemotherapeutic agent for KRAS^G12D-driven lung cancer treatment.

Yuel: Improving the Generalizability of Structure-Free Compound-Protein Interaction Prediction

Predicting binding affinities between small molecules and the protein target is at the core of computational drug screening and drug target identification. Deep learning-based approaches have recently been adapted to predict binding affinities and they claim to achieve high prediction accuracy in their tests; we show that these approaches do not generalize, that is, they fail to predict interactions between unknown proteins and unknown small molecules. To address these shortcomings, we develop a new compound-protein interaction predictor, Yuel, which predicts compound-protein interactions with a higher generalizability than the existing methods. Upon comprehensive tests on various data sets, we find that out of all the deep-learning approaches surveyed, Yuel manifests the best ability to predict interactions between unknown compounds and unknown proteins.

wSDTNBI: a novel network-based inference method for virtual screening

In recent years, the rapid development of network-based methods for the prediction of drug-target interactions (DTIs) provides an opportunity for the emergence of a new type of virtual screening (VS), namely, network-based VS. Herein, we reported a novel network-based inference method named wSDTNBI. Compared with previous network-based methods that use unweighted DTI networks, wSDTNBI uses weighted DTI networks whose edge weights are correlated with binding affinities. A two-pronged approach based on weighted DTI and drug-substructure association networks was employed to calculate prediction scores. To show the practical value of wSDTNBI, we performed network-based VS on retinoid-related orphan receptor γt (RORγt), and purchased 72 compounds for experimental validation. Seven of the purchased compounds were confirmed to be novel RORγt inverse agonists by in vitro experiments, including ursonic acid and oleanonic acid with IC50 values of 10 nM and 0.28 μM, respectively. Moreover, the direct contact between ursonic acid and RORγt was confirmed using the X-ray crystal structure, and in vivo experiments demonstrated that ursonic acid and oleanonic acid have therapeutic effects on multiple sclerosis. These results indicate that wSDTNBI might be a powerful tool for network-based VS in drug discovery.