Combinatorial Chemistry & High Throughput Screening - Volume 18, Issue 7, 2015

Target based virtual screening has surpassed ligand based virtual screening methods in the recent past mainly as it provides more clues regarding intermolecular interactions and takes into consideration the flexible receptor as well. The current methodology describes a computational strategy of predicting Mycobacterium tuberculosis (M. tuberculosis) binders for five well studied targets representing M. tuberculosis proteome encompassing most of the known mechanisms of action. The diversity of the targets was affirmed by their active site analysis and structural studies. The current approach employed pharmacophore searching, docking and clustering techniques in tandem and was validated by enrichment studies using the available Schrödinger data set consisting of 1000 decoys. The application of this methodology was demonstrated by predicting potential molecular targets for fifty newly synthesized compounds. Cross docking studies on the targets were carried out with 4512 known inhibitors utilizing a high performance computing platform to reveal underlying affinity and promiscuity patterns. Optimum binding energy range for all targets as determined by high throughput docking was found to be -3 to -13 kcal/mol.

Every drug discovery research program involves synthesis of a novel and potential drug molecule utilizing atom efficient, economical and environment friendly synthetic strategies. The current work focuses on the role of the reactivity based fingerprints of compounds as filters for virtual screening using a tool ChemScore. A reactant-like (RLS) and a product- like (PLS) score can be predicted for a given compound using the binary fingerprints derived from the numerous known organic reactions which capture the molecule-molecule interactions in the form of addition, substitution, rearrangement, elimination and isomerization reactions. The reaction fingerprints were applied to large databases in biology and chemistry, namely ChEMBL, KEGG, HMDB, DSSTox, and the Drug Bank database. A large network of 1113 synthetic reactions was constructed to visualize and ascertain the reactant product mappings in the chemical reaction space. The cumulative reaction fingerprints were computed for 4000 molecules belonging to 29 therapeutic classes of compounds, and these were found capable of discriminating between the cognition disorder related and anti-allergy compounds with reasonable accuracy of 75% and AUC 0.8. In this study, the transition state based fingerprints were also developed and used effectively for virtual screening in drug related databases. The methodology presented here provides an efficient handle for the rapid scoring of molecular libraries for virtual screening.

The ligand-based virtual screening of combinatorial libraries employs a number of statistical modeling and machine learning methods. A comprehensive analysis of the application of these methods for the diversity oriented virtual screening of biological targets/drug classes is presented here. A number of classification models have been built using three types of inputs namely structure based descriptors, molecular fingerprints and therapeutic category for performing virtual screening. The activity and affinity descriptors of a set of inhibitors of four target classes DHFR, COX, LOX and NMDA have been utilized to train a total of six classifiers viz. Artificial Neural Network (ANN), k nearest neighbor (k- NN), Support Vector Machine (SVM), Naïve Bayes (NB), Decision Tree - (DT) and Random Forest - (RF). Among these classifiers, the ANN was found as the best classifier with an AUC of 0.9 irrespective of the target. New molecular fingerprints based on pharmacophore, toxicophore and chemophore (PTC), were used to build the ANN models for each dataset. A good accuracy of 87.27% was obtained using 296 chemophoric binary fingerprints for the COX-LOX inhibitors compared to pharmacophoric (67.82 %) and toxicophoric (70.64 %). The methodology was validated on the classical Ames mutagenecity dataset of 4337 molecules. To evaluate it further, selectivity and promiscuity of molecules from five drug classes viz. anti-anginal, anti-convulsant, anti-depressant, anti-arrhythmic and anti-diabetic were studied. The TPC fingerprints computed for each category were able to capture the drug-class specific features using the k-NN classifier. These models can be useful for selecting optimal molecules for drug design.

Natural products obtained from marine sources are considered to be a rich and diverse source of potential drugs. In the present work we demonstrate the use of chemoinformatics approach for the design of new molecules inspired by molecules from marine organisms. Accordingly we have assimilated information from two major scientific domains namely chemoinformatics and biodiversity informatics to develop an interactive marine database named MIMMO (Medicinally Important Molecules from Marine Organisms). The database can be queried for species, molecules, scaffolds, drugs, diseases and associated cumulative biological activity spectrum along with links to the literature resources. Molecular informatics analysis of the molecules obtained from MIMMO was performed to study their chemical space. The distinct skeletal features of the biologically active compounds isolated from marine species were identified. Scaffold molecules and species networks were created to identify common scaffolds from marine source and drug space. An analysis of the entire molecular data revealed a unique list of around 2000 molecules from which ten most frequently occurring distinct scaffolds were obtained.

The human 5-hydroxytryptamine receptor subtype 1A (5-HT1A) is highly expressed in the raphe nuclei region and limbic structures; for that reason 5-HT1A has served as a promising target for treating human mood disorders and neurodegenerative diseases. We have developed binary quantitative structure-activity relationship (QSAR) models for 5- HT1A binding using data retrieved from the WOMBAT database and the k-Nearest Neighbor (kNN) machine learning method. A rigorous QSAR modeling and screening workflow had been followed, with extensive internal and external validation processes. The models’ classification accuracies to discriminate 5-HT1A binders from the non-binders are as high as 96% for the external validation. These models were employed further to mine two major natural products screening libraries, i.e. TimTec Natural Product Library (NPL) and Natural Derivatives Library (NDL). In the end five screening hits were tested by radioligand binding assays with a success rate of 40%, and two Library compounds were confirmed to be binders at the μM concentration against the human 5-HT1A receptor. The combined application of rigorous QSAR modeling and model-based virtual screening presents a powerful means for profiling natural products compounds with important biomedical activities.

Histone deacetylases (HDACs) are part of a vast family of enzymes with crucial roles in numerous biological processes, largely through their repressive influence on transcription, with serious implications in a variety of human diseases. Among different isoforms, human HDAC2 in particular draws attention as a promising target for the treatment of cancer and memory deficits associated with neurodegenerative diseases. Now the challenge is to obtain a compound that is structurally novel and truly selective to HDAC2 because most current HDAC2 inhibitors do not show isoforms selectivity and suffer from metabolic instability. In order to identify novel, and isoform-selective inhibitors for human HDAC2, we designed a shape-based hybrid query from multiple scaffolds of known chemical classes and validated it to be a more effective approach to discover diverse scaffolds than single-molecule query. The hybrid query-based screening rendered a hit compound with the N-benzylaniline scaffold which showed moderate inhibitory activity against HDAC2, and its chemical structure is diverse compared to known HDAC2 inhibitors. Notably, this compound shows the selectivity against the HDAC6, a Class II enzyme, thus has the potential to further develop into the class- and isoform-selective inhibitors. Our present study supplies an useful approach to identifying novel HDAC2 inhibitors, and can be extended to the inquires of other important biomedical targets as well.