Assessing therapeutic potential of molecules: molecular property diagnostic suite for tuberculosis $$(\mathbf{MPDS}^{\mathbf{TB}})$$ ( MPDS TB )

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/.

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 based dynamic consensus model for predicting blood-brain barrier permeability

The blood-brain barrier (BBB) is an important defence mechanism that restricts disease-causing pathogens and toxins to enter the brain from the bloodstream. In recent years, many in silico methods were proposed for predicting BBB permeability, however, the reliability of these models is questionable due to the smaller and class-imbalance dataset which subsequently leads to a very high false positive rate. In this study, machine learning and deep learning-based predictive models were built using XGboost, Random Forest, Extra-tree classifiers and deep neural network. A dataset of 8153 compounds comprising both the BBB permeable and BBB non-permeable was curated and subjected to calculations of molecular descriptors and fingerprints for generating the features for machine learning and deep learning models. Three balancing techniques were then applied to the dataset to address the class-imbalance issue. A comprehensive comparison among the models showed that the deep neural network model generated on the balanced MACCS fingerprint dataset outperformed with an accuracy of 97.8% and a ROC-AUC score of 0.98 among all the models. Additionally, a dynamic consensus model was prepared with the machine learning models and validated with a benchmark dataset for predicting BBB permeability with higher confidence scores.

Towards systematic exploration of chemical space: building the fragment library module in molecular property diagnostic suite

A fragment-based drug discovery (FBDD) approach has traditionally been of utmost significance in drug design studies. It allows the exploration of large chemical space to find novel scaffolds and chemotypes which can be improved into selective inhibitors with good affinity. In the current work, several public domain chemical libraries (ChEMBL, DrugCentral, PDB ligands, COCONUT, and SAVI) comprising bioactive and virtual molecules were retrieved to develop a fragment library. A systematic fragmentation method that breaks a given molecule into rings, linkers, and substituents was used to cleave the molecules and the fragments were analyzed. Further, only the ring framework was taken into the consideration to develop a fragment library that consists of a total number of 107,614 unique fragments. This set represents a rich diverse structure framework that covers a wide variety of yet-to-be-explored fragments for a wide range of small molecule-based applications. This fragment library is an integral part of the molecular property diagnostic suite (MPDS) suite that can be used with other modeling and informatics methods for FBDD approaches. The fragment library module of MPDS can be accessed at http://mpds.neist.res.in:8085.

Drug repositioning for anti-tuberculosis drugs: an in silico polypharmacology approach

Development of potential antitubercular molecules is a challenging task due to the rapidly emerging drug-resistant strains of Mycobacterium tuberculosis (M.tb). Structure-based approaches hold greater benefit in identifying compounds/drugs with desired polypharmacological profiles. These methods can be employed based on the knowledge of protein binding sites to identify the complementary ligands. In this study, polypharmacology guided computational drug repurposing approach was applied to identify potential antitubercular drugs. 20 important druggable protein targets in M.tb were considered from the target library of Molecular Property Diagnostic Suite-Tuberculosis (MPDS^TB- http://mpds.neist.res.in:8084 ) for virtual screening. FDA approved drugs were collected, preprocessed and docked in the active sites of the 20 M.tb targets. The top 300 drug molecules from each target (20 × 300) were filtered-in and subsequently screened for possible antitubercular and antimycobacterial activity using PASS tool. Using this approach, 34 drugs with predicted antitubercular and anti-mycobacterial activity were identified along with good binding affinity against multiple M.tb targets. Interestingly, 21 out of the 34 identified drugs are antibiotics while 4 drug molecules (nitrofural, stavudine, quinine and quinidine) are non-antibiotics showing promising predicted antitubercular activity. Most of these molecules have the similar privileged antimycobacterial drugs scaffold. Further drug likeness properties were calculated to get deeper insights to M.tb lead molecules. Interestingly, it was also observed that the drugs identified from the study are under different stages of drug discovery (i.e., in vitro, clinical trials) for the effective treatment of various diseases including cancer, degenerative diseases, dengue virus infection, tuberculosis, etc. Krasavin et al., 2017 synthesized nitrofuran analogues with appreciable MICs (22-23 µM) against M.tb H37Rv. These experiments further add to the credibility of the drugs identified in this study (TB).

Mycobacterium tuberculosis Cell Wall Permeability Model Generation Using Chemoinformatics and Machine Learning Approaches

The drug-resistant strains of Mycobacterium tuberculosis (M.tb) are evolving at an alarming rate, and this indicates the urgent need for the development of novel antitubercular drugs. However, genetic mutations, complex cell wall system of M.tb, and influx-efflux transporter systems are the major permeability barriers that significantly affect the M.tb drugs activity. Thus, most of the small molecules are ineffective to arrest the M.tb cell growth, even though they are effective at the cellular level. To address the permeability issue, different machine learning models that effectively distinguish permeable and impermeable compounds were developed. The enzyme-based (IC₅₀) and cell-based (minimal inhibitory concentration) data were considered for the classification of M.tb permeable and impermeable compounds. It was assumed that the compounds that have high activity in both enzyme-based and cell-based assays possess the required M.tb cell wall permeability. The XGBoost model was outperformed when compared to the other models generated from different algorithms such as random forest, support vector machine, and naïve Bayes. The XGBoost model was further validated using the validation data set (21 permeable and 19 impermeable compounds). The obtained machine learning models suggested that various descriptors such as molecular weight, atom type, electrotopological state, hydrogen bond donor/acceptor counts, and extended topochemical atoms of molecules are the major determining factors for both M.tb cell permeability and inhibitory activity. Furthermore, potential antimycobacterial drugs were identified using computational drug repurposing. All the approved drugs from DrugBank were collected and screened using the developed permeability model. The screened compounds were given as input in the PASS server for the identification of possible antimycobacterial compounds. The drugs that were retained after two filters were docked to the active site of 10 different potential antimycobacterial drug targets. The results obtained from this study may improve the understanding of M.tb permeability and activity that may aid in the development of novel antimycobacterial drugs.

Repurposing of approved drug molecules for viral infectious diseases: a molecular modelling approach

The identification of new viral drugs has become a task of paramount significance due to the frequent occurrence of viral infections and especially during the current pandemic. Despite the recent advancements, the development of antiviral drugs has not made parallel progress. Reduction of time frame and cost of the drug development process is the major advantage of drug repurposing. Therefore, in this study, a drug repurposing strategy using molecular modelling techniques, i.e. biological activity prediction, virtual screening, and molecular dynamics simulation was employed to find promising repurposing candidates for viral infectious diseases. The biological activities of non-redundant (4171) drug molecules were predicted using PASS analysis, and 1401 drug molecules were selected which showed antiviral activities in the analysis. These drug molecules were subjected to virtual screening against the selected non-structural viral proteins. A series of filters, i.e. top 10 drug molecules based on binding affinity, mean value of binding affinity, visual inspection of protein-drug complexes, and number of H-bond between protein and drug molecules were used to narrow down the drug molecules. Molecular dynamics simulation analysis was carried out to validate the intrinsic atomic interactions and binding conformations of protein-drug complexes. The binding free energies of drug molecules were assessed by employing MMPBSA analysis. Finally, nine drug molecules were prioritized, as promising repurposing candidates with the potential to inhibit the selected non-structural viral proteins.Communicated by Ramaswamy H. Sarma.

Hybrid Dynamic Pharmacophore Models as Effective Tools to Identify Novel Chemotypes for Anti-TB Inhibitor Design: A Case Study With Mtb-DapB

Antimicrobial resistance (AMR) is one of the most serious global public health threats as it compromises the successful treatment of deadly infectious diseases like tuberculosis. New therapeutics are constantly needed but it takes a long time and is expensive to explore new biochemical space. One way to address this issue is to repurpose the validated targets and identify novel chemotypes that can simultaneously bind to multiple binding pockets of these targets as a new lead generation strategy. This study reports such a strategy, dynamic hybrid pharmacophore model (DHPM), which represents the combined interaction features of different binding pockets contrary to the conventional approaches, where pharmacophore models are generated from single binding sites. We have considered Mtb-DapB, a validated mycobacterial drug target, as our model system to explore the effectiveness of DHPMs to screen novel unexplored compounds. Mtb-DapB has a cofactor binding site (CBS) and an adjacent substrate binding site (SBS). Four different model systems of Mtb-DapB were designed where, either NADPH/NADH occupies CBS in presence/absence of an inhibitor 2, 6-PDC in the adjacent SBS. Two more model systems were designed, where 2, 6-PDC was linked to NADPH and NADH to form hybrid molecules. The six model systems were subjected to 200 ns molecular dynamics simulations and trajectories were analyzed to identify stable ligand-receptor interaction features. Based on these interactions, conventional pharmacophore models (CPM) were generated from the individual binding sites while DHPMs were created from hybrid-molecules occupying both binding sites. A huge library of 1,563,764 publicly available molecules were screened by CPMs and DHPMs. The screened hits obtained from both types of models were compared based on their Hashed binary molecular fingerprints and 4-point pharmacophore fingerprints using Tanimoto, Cosine, Dice and Tversky similarity matrices. Molecules screened by DHPM exhibited significant structural diversity, better binding strength and drug like properties as compared to the compounds screened by CPMs indicating the efficiency of DHPM to explore new chemical space for anti-TB drug discovery. The idea of DHPM can be applied for a wide range of mycobacterial or other pathogen targets to venture into unexplored chemical space.

Molecular property diagnostic suite (MPDS): Development of disease-specific open source web portals for drug discovery

Molecular property diagnostic suite (MPDS) is a Galaxy-based open source drug discovery and development platform. MPDS web portals are designed for several diseases, such as tuberculosis, diabetes mellitus, and other metabolic disorders, specifically aimed to evaluate and estimate the drug-likeness of a given molecule. MPDS consists of three modules, namely data libraries, data processing, and data analysis tools which are configured and interconnected to assist drug discovery for specific diseases. The data library module encompasses vast information on chemical space, wherein the MPDS compound library comprises 110.31 million unique molecules generated from public domain databases. Every molecule is assigned with a unique ID and card, which provides complete information for the molecule. Some of the modules in the MPDS are specific to the diseases, while others are non-specific. Importantly, a suitably altered protocol can be effectively generated for another disease-specific MPDS web portal by modifying some of the modules. Thus, the MPDS suite of web portals shows great promise to emerge as disease-specific portals of great value, integrating chemoinformatics, bioinformatics, molecular modelling, and structure- and analogue-based drug discovery approaches.

A QSAR and molecular modelling study towards new lead finding: polypharmacological approach to Mycobacterium tuberculosis

Developing effective inhibitors against Mycobacterium tuberculosis (Mtb) is a challenging task, primarily due to the emergence of resistant strains. In this study, we have proposed and implemented an in silico guided polypharmacological approach, which is expected to be effective against resistant strains by simultaneously inhibiting several potential Mtb drug targets. A combination of pharmacophore and QSAR based virtual screening strategy taking three key targets such as InhA (enoyl-acyl-carrier-protein reductase), GlmU (N-acetyl-glucosamine-1-phosphate uridyltransferase) and DapB (dihydrodipicolinate reductase) have resulted in initial 784 hits from Asinex database of 435,000 compounds. These hits were further subjected to docking with 33 Mtb druggable targets. About 110 potential polypharmacological hits were taken by integrating the aforementioned screening protocols. Further screening was conducted by taking various parameters and properties such as cell permeability, drug-likeness, drug-induced phospholipidosisand structural alerts. A consensus analysis has yielded 59 potential hits that pass through all the filters and can be prioritized for effective drug-resistant tuberculosis. This study proposes about nine potential hits which are expected to be promising molecules, having not only drug-like properties, but also being effective against multiple Mtb targets.