Virtual screening and molecular dynamics investigations using natural compounds against autotaxin for the treatment of chronic pain

Chronic pain is a common and debilitating condition with a huge social and economic burden worldwide. Currently, available drugs in clinics are not adequately effective and possess a variety of severe side effects leading to treatment withdrawal and poor quality of life. Recent findings highlight the potential role of autotaxin (ATX) as a promising novel target for chronic pain management, extending beyond its previously established involvement in arthritis and other neurological disorders, such as Alzheimer's disease. In the present study, we used a virtual screening strategy by targeting ATX against commercially available natural compounds (enamine- phenotypic screening library) to identify the potential inhibitors for the treatment of chronic pain. After initial identification using molecular docking based virtual screening, molecular mechanics (MM/GBSA), ADMET profiling and molecular dynamics simulation were performed to verify top hits. The computational screening resulted in the identification of fifteen top scoring structurally diverse hits that have free energy of binding (ΔG) values in the range of -25.792 (for compound Enamine_1850) to -74.722 Kcal/mol (for compound Enamine_1687). Moreover, the top-scoring hits have favourable ADME properties as calculated using in-silico algorithms. Additionally, the molecular dynamics simulation revealed the stable nature of protein-ligand interaction and provided information about amino acid residues involved in binding. This study led to the identification of potential autotaxin inhibitors with favourable pharmacokinetic properties. Identified hits may further be investigated for their safety and efficacy potential using in-vitro and in-vivo models of chronic pain.Communicated by Ramaswamy H. Sarma.

Structure-based virtual screening and molecular dynamics simulation for the identification of sphingosine kinase-2 inhibitors as potential analgesics

Neuropathic pain is due to an injury or disease of the somatosensory nervous system, which accounts for a significant economical and health burden to society. Due to poor understanding of their underlying mechanisms, the available treatments merely provide symptomatic relief and precipitates a variety of adverse effects. This suggests that there is an unmet medical need that must be addressed with effective strategies for the development of novel therapeutics. Sphingosine kinase 2 (SphK2) is an oncogenic lipid kinase that has emerged as a promising target for chronic pain and other diseases. In the present study, we have explored the structure-based virtual high-throughput screening of the Nuclei of Bioassays, Ecophysiology, and Biosynthesis of Natural Products Database (NuBBE) to identify potent natural products as inhibitors of SphK2. A molecular docking study was performed to calculate binding affinities and specificity to identify potential leads against SphK2. Initially, hits were selected by the implementation of absorption, distribution, metabolism, excretion and toxicity properties, Lipinski rule, and PAINS filters. The top-scoring hits also exhibiting an optimal ADMET profile were subjected to MM/GBSA free binding free energy calculation and molecular dynamics simulation. The results from molecular dynamics simulation revealed a stable ligand -SphK2 complex with protein and ligand RMSD within reasonable limits. Overall, we identified compounds, NuBBE_972 and NuBBE_1107 as potential inhibitors of SphK2 with optimal pharmacokinetic properties which have the potential to be developed as novel therapeutics for the management of chronic pain.Communicated by Ramaswamy H. Sarma.

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

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

Current trends in computer aided drug design and a highlight of drugs discovered via computational techniques: A review

Computer-aided drug design (CADD) is one of the pivotal approaches to contemporary pre-clinical drug discovery, and various computational techniques and software programs are typically used in combination, in a bid to achieve the desired outcome. Several approved drugs have been developed with the aid of CADD. On SciFinder®, we evaluated more than 600 publications through systematic searching and refining, using the terms, virtual screening; software methods; computational studies and publication year, in order to obtain data concerning particular aspects of CADD. The primary focus of this review was on the databases screened, virtual screening and/or molecular docking software program used. Furthermore, we evaluated the studies that subsequently performed molecular dynamics (MD) simulations and we reviewed the software programs applied, the application of density functional theory (DFT) calculations and experimental assays. To represent the latest trends, the most recent data obtained was between 2015 and 2020, consequently the most frequently employed techniques and software programs were recorded. Among these, the ZINC database was the most widely preferred with an average use of 31.2%. Structure-based virtual screening (SBVS) was the most prominently used type of virtual screening and it accounted for an average of 57.6%, with AutoDock being the preferred virtual screening/molecular docking program with 41.8% usage. Following the screening process, 38.5% of the studies performed MD simulations to complement the virtual screening and GROMACS with 39.3% usage, was the popular MD software program. Among the computational techniques, DFT was the least applied whereby it only accounts for 0.02% average use. An average of 36.5% of the studies included reports on experimental evaluations following virtual screening. Ultimately, since the inception and application of CADD in pre-clinical drug discovery, more than 70 approved drugs have been discovered, and this number is steadily increasing over time.

Chemical similarity methods for analyzing secondary metabolite structures

The relation that exists between the structure of a compound and its function is an integral part of chemoinformatics. The similarity principle states that “structurally similar molecules tend to have similar properties and similar molecules exert similar biological activities”. The similarity of the molecules can either be studied at the structure level or at the descriptor level (properties level). Generally, the objective of chemical similarity measures is to enhance prediction of the biological activities of molecules. In this article, an overview of various methods used to compare the similarity between metabolite structures has been provided, including two-dimensional (2D) and three-dimensional (3D) approaches. The focus has been on methods description; e.g. fingerprint-based similarity in which the molecules under study are first fragmented and their fingerprints are computed, 2D structural similarity by comparing the Tanimoto coefficients and Euclidean distances, as well as the use of physiochemical properties descriptor-based similarity methods. The similarity between molecules could also be measured by using data mining (clustering) techniques, e.g. by using virtual screening (VS)-based similarity methods. In this approach, the molecules with the desired descriptors or /and structures are screened from large databases. Lastly, SMILES-based chemical similarity search is an important method for studying the exact structure search, substructure search and also descriptor similarity. The use of a particular method depends upon the requirements of the researcher.

Cheminformatics to Characterize Pharmacologically Active Natural Products

Natural products have a significant role in drug discovery. Natural products have distinctive chemical structures that have contributed to identifying and developing drugs for different therapeutic areas. Moreover, natural products are significant sources of inspiration or starting points to develop new therapeutic agents. Natural products such as peptides and macrocycles, and other compounds with unique features represent attractive sources to address complex diseases. Computational approaches that use chemoinformatics and molecular modeling methods contribute to speed up natural product-based drug discovery. Several research groups have recently used computational methodologies to organize data, interpret results, generate and test hypotheses, filter large chemical databases before the experimental screening, and design experiments. This review discusses a broad range of chemoinformatics applications to support natural product-based drug discovery. We emphasize profiling natural product data sets in terms of diversity; complexity; acid/base; absorption, distribution, metabolism, excretion, and toxicity (ADME/Tox) properties; and fragment analysis. Novel techniques for the visual representation of the chemical space are also discussed.