Comprehensive molecular docking and dynamic simulations for drug repurposing of clinical drugs against multiple cancer kinase targets.

Drug repurposing is a method to identify novel therapeutic agents from the existing drugs and clinical compounds. In the present comprehensive work, molecular docking, virtual screening and dynamics simulations were carried out for ten cancer types viz breast, colon, central nervous system, leukaemia, melanoma, ovarian, prostate, renal and lung (non-small and small cell) against validated eighteen kinase targets. The study aims to understand the action of chemotherapy drugs mechanism through binding interactions against selected targets via comparative docking simulations with the state-art molecular modelling suits such as MOE, Cresset-Flare, AutoDock Vina, GOLD and GLIDE. Chemotherapeutic drugs (n = 112) were shortlisted from standard drug databases with appropriate chemoinformatic filters. Based on docking studies it was revealed that leucovorin, nilotinib, ellence, thalomid and carfilzomib drugs possessed potential against other cancer targets. A library was built to enumerate novel molecules based on the scaffold and functional groups extracted from known drugs and clinical compounds. Twenty novel molecules were prioritised further based on drug-like attributes. These were cross docked against 1MQ4 Aurora-A Protein Kinase for prostate cancer and 4UYA Mitogen-activated protein kinase for renal cancer. All docking programs yielded similar results but interestingly AutoDock Vina yielded the lowest RMSD with the native ligand. To further validate the final docking results at atomistic level, molecular dynamics simulations were performed to ascertain the stability of the protein-ligand complex. The study enables repurposing of drugs and lead identification by employing a host of structure and ligand based virtual screening tools and techniques.Communicated by Ramaswamy H. Sarma.

A Review on Computational Analysis of Big Data in Breast Cancer for Predicting Potential Biomarkers.

Breast cancer is the most predominantly occurring cancer in the world. Several genes and proteins have been recently studied to predict biomarkers that enable early disease identification and monitor its recurrence. In the era of high-throughput technology, studies show several applications of big data for identifying potential biomarkers. The review aims to provide a comprehensive overview of big data analysis in breast cancer towards the prediction of biomarkers with emphasis on computational methods like text mining, network analysis, next-generation sequencing technology (NGS), machine learning (ML), deep learning (DL), and precision medicine. Integrating data from various computational approaches enables the stratification of cancer patients and the identification of molecular signatures in cancer and their subtypes. The computational methods and statistical analysis help expedite cancer prognosis and develop precision cancer medicine (PCM). As a part of case study in the present work, we constructed a large gene-drug interaction network to predict new biomarkers genes. The gene-drug network helped us to identify eight genes that could serve as novel potential biomarkers.

BESFA: bioinformatics based evolutionary, structural & functional analysis of prostrate, Placenta, Ovary, Testis, and Embryo (POTE) paralogs.

The POTE family comprises 14 paralogues and is primarily expressed in Prostrate, Placenta, Ovary, Testis, Embryo (POTE), and cancerous cells. The prospective function of the POTE protein family under physiological conditions is less understood. We systematically analyzed their cellular localization and molecular docking analysis to elucidate POTE proteins' structure, function, and Adaptive Divergence. Our results suggest that group three POTE paralogs (POTEE, POTEF, POTEI, POTEJ, and POTEKP (a pseudogene)) exhibits significant variation among other members could be because of their Adaptive Divergence. Furthermore, our molecular docking studies on POTE protein revealed the highest binding affinity with NCI-approved anticancer compounds. Additionally, POTEE, POTEF, POTEI, and POTEJ were subject to an explicit molecular dynamic simulation for 50ns. MM-GBSA and other essential electrostatics were calculated that showcased that only POTEE and POTEF have absolute binding affinities with minimum energy exploitation. Thus, this study's outcomes are expected to drive cancer research to successful utilization of POTE genes family as a new biomarker, which could pave the way for the discovery of new therapies.

Design of Novel Drug-like Molecules Using Informatics Rich Secondary Metabolites Analysis of Indian Medicinal and Aromatic Plants.

BACKGROUND: Several medicinal plants are being used in Indian medicine systems from ancient times. However, in most cases, the specific molecules or the active ingredients responsible for the medicinal or therapeutic properties are not yet known. OBJECTIVE: This study aimed to report a computational protocol as well as a tool for generating novel potential drug candidates from the bioactive molecules of Indian medicinal and aromatic plants through the chemoinformatics approach. METHODS: We built a database of the Indian medicinal and aromatic plants coupled with associated information (plant families, plant parts used for the medicinal purpose, structural information, therapeutic properties, etc.) We also developed a Java-based chemoinformatics open-source tool called DoMINE (Database of Medicinally Important Natural products from plantaE) for the generation of virtual library and screening of novel molecules from known medicinal plant molecules. We employed chemoinformatics approaches to in-silico screened metabolites from 104 Indian medicinal and aromatic plants and designed novel drug-like bioactive molecules. For this purpose, 1665 ring containing molecules were identified by text mining of literature related to the medicinal plant species, which were later used to extract 209 molecular scaffolds. Different scaffolds were further used to build a focused virtual library. Virtual screening was performed with cluster analysis to predict drug-like and lead-like molecules from these plant molecules in the context of drug discovery. The predicted drug-like and lead-like molecules were evaluated using chemoinformatics approaches and statistical parameters, and only the most significant molecules were proposed as the candidate molecules to develop new drugs. RESULTS AND CONCLUSION: The supra network of molecules and scaffolds identifies the relationship between the plant molecules and drugs. Cluster analysis of virtual library molecules showed that novel molecules had more pharmacophoric properties than toxicophoric and chemophoric properties. We also developed the DoMINE toolkit for the advancement of natural product-based drug discovery through chemoinformatics approaches. This study will be useful in developing new drug molecules from the known medicinal plant molecules. Hence, this work will encourage experimental organic chemists to synthesize these molecules based on the predicted values. These synthesized molecules need to be subjected to biological screening to identify potential molecules for drug discovery research.

Zebrafish acid ceramidase: Expression in Pichia pastoris GS115and biochemical characterization.

Acid ceramidase (N-acylsphingosine deacylase EC 3.5.1.23; AC) catalyzes the hydrolysis of ceramide into sphingosine (SPH) and free fatty acid. Zebrafish acid ceramidase (AC) has 60% homology with the human AC). Mutations in the human AC gene asah1 are known to cause Farber disease and spinal muscular atrophy with progressive myoclonic epilepsy. Zebrafish AC was overexpressed in Pichia pastoris by inserting asah1b gene into the genome. The majority of the overexpressed enzyme was secreted into the culture medium and purified to apparent homogeneity by stepwise chromatography. The recombinant protein was glycosylated precursor, that further undergoes limited autoproteolytic processing into two subunits (α and ß) which are visible in SDS-PAGE. The zebrafish AC is heterodimer associated with an inter-subunit disulfide bond. SDS-PAGE estimated the mass of native enzyme to be approximately 50 kDa & size exclusion chromatography estimated the mass of the active enzyme as approximately 100 kDa, suggesting the formation of a dimer of heterodimers. The protein was secreted as a mixture of processed and unprocessed forms in the culture media. A preliminary characterization of purified zebrafish AC was done by an enzyme assay. The zebrafish AC expressed in Pichia pastoris would be used for further structural and functional analysis.

Application of Genetic Programming (GP) Formalism for Building Disease Predictive Models from Protein-Protein Interactions (PPI) Data.

Protein-protein interactions (PPIs) play a vital role in the biological processes involved in the cell functions and disease pathways. The experimental methods known to predict PPIs require tremendous efforts and the results are often hindered by the presence of a large number of false positives. Herein, we demonstrate the use of a new Genetic Programming (GP) based Symbolic Regression (SR) approach for predicting PPIs related to a disease. In a case study, a dataset consisting of one hundred and thirty five PPI complexes related to cancer was used to construct a generic PPI predicting model with good PPI prediction accuracy and generalization ability. A high correlation coefficient(CC) of 0.893, low root mean square error (RMSE) and mean absolute percentage error (MAPE) values of 478.221 and 0.239, respectively were achieved for both the training and test set outputs. To validate the discriminatory nature of the model, it was applied on a dataset of diabetes complexes where it yielded significantly low CC values. Thus, the GP model developed here serves a dual purpose: (a)a predictor of the binding energy of cancer related PPI complexes, and (b)a classifier for discriminating PPI complexes related to cancer from those of other diseases.

ChemEngine: harvesting 3D chemical structures of supplementary data from PDF files.

Digital access to chemical journals resulted in a vast array of molecular information that is now available in the supplementary material files in PDF format. However, extracting this molecular information, generally from a PDF document format is a daunting task. Here we present an approach to harvest 3D molecular data from the supporting information of scientific research articles that are normally available from publisher's resources. In order to demonstrate the feasibility of extracting truly computable molecules from PDF file formats in a fast and efficient manner, we have developed a Java based application, namely ChemEngine. This program recognizes textual patterns from the supplementary data and generates standard molecular structure data (bond matrix, atomic coordinates) that can be subjected to a multitude of computational processes automatically. The methodology has been demonstrated via several case studies on different formats of coordinates data stored in supplementary information files, wherein ChemEngine selectively harvested the atomic coordinates and interpreted them as molecules with high accuracy. The reusability of extracted molecular coordinate data was demonstrated by computing Single Point Energies that were in close agreement with the original computed data provided with the articles. It is envisaged that the methodology will enable large scale conversion of molecular information from supplementary files available in the PDF format into a collection of ready- to- compute molecular data to create an automated workflow for advanced computational processes. Software along with source codes and instructions available at https://sourceforge.net/projects/chemengine/files/?source=navbar.Graphical abstract.

Building and analysis of protein-protein interactions related to diabetes mellitus using support vector machine, biomedical text mining and network analysis.

In order to understand the molecular mechanism underlying any disease, knowledge about the interacting proteins in the disease pathway is essential. The number of revealed protein-protein interactions (PPI) is still very limited compared to the available protein sequences of different organisms. Experiment based high-throughput technologies though provide some data about these interactions, those are often fairly noisy. Computational techniques for predicting protein-protein interactions therefore assume significance. 1296 binary fingerprints that encode a combination of structural and geometric properties were developed using the crystallographic data of 15,000 protein complexes in the pdb server. In a case study, these fingerprints were created for proteins implicated in the Type 2 diabetes mellitus disease. The fingerprints were input into a SVM based model for discriminating disease proteins from non disease proteins yielding a classification accuracy of 78.2% (AUC value of 0.78) on an external data set composed of proteins retrieved via text mining of diabetes related literature. A PPI network was constructed and analysed to explore new disease targets. The integrated approach exemplified here has a potential for identifying disease related proteins, functional annotation and other proteomics studies.

Role of Open Source Tools and Resources in Virtual Screening for Drug Discovery.

Advancement in chemoinformatics research in parallel with availability of high performance computing platform has made handling of large scale multi-dimensional scientific data for high throughput drug discovery easier. In this study we have explored publicly available molecular databases with the help of open-source based integrated in-house molecular informatics tools for virtual screening. The virtual screening literature for past decade has been extensively investigated and thoroughly analyzed to reveal interesting patterns with respect to the drug, target, scaffold and disease space. The review also focuses on the integrated chemoinformatics tools that are capable of harvesting chemical data from textual literature information and transform them into truly computable chemical structures, identification of unique fragments and scaffolds from a class of compounds, automatic generation of focused virtual libraries, computation of molecular descriptors for structure-activity relationship studies, application of conventional filters used in lead discovery along with in-house developed exhaustive PTC (Pharmacophore, Toxicophores and Chemophores) filters and machine learning tools for the design of potential disease specific inhibitors. A case study on kinase inhibitors is provided as an example.

Protein Ligand Complex Guided Approach for Virtual Screening.

The target ligand association data is a rich source of information which is not exploited enough for drug design efforts in virtual screening. A java based open-source toolkit for Protein Ligand Network Extraction (J-ProLiNE) focused on protein-ligand complex analysis with several features integrated in a distributed computing network has been developed. Sequence alignment and similarity search components have been automated to yield local, global alignment scores along with similarity and distance scores. 10000 proteins with co-crystallized ligands from pdb and MOAD databases were extracted and analyzed for revealing relationships between targets, ligands and scaffolds. Through this analysis, we could generate a protein ligand network to identify the promiscuous and selective scaffolds for multiple classes of proteins targets. Using J-ProLiNE we created a 507 x 507 matrix of protein targets and native ligands belonging to six enzyme classes and analyzed the results to elucidate the protein-protein, protein-ligand and ligand-ligand interactions. In yet another application of the J-ProLiNE software, we were able to process kinase related information stored in US patents to construct disease-gene-ligand-scaffold networks. It is hoped that the studies presented here will enable target ligand knowledge based virtual screening for inhibitor design.

ChemScreener: A Distributed Computing Tool for Scaffold based Virtual Screening.

In this work we present ChemScreener, a Java-based application to perform virtual library generation combined with virtual screening in a platform-independent distributed computing environment. ChemScreener comprises a scaffold identifier, a distinct scaffold extractor, an interactive virtual library generator as well as a virtual screening module for subsequently selecting putative bioactive molecules. The virtual libraries are annotated with chemophore-, pharmacophore- and toxicophore-based information for compound prioritization. The hits selected can then be further processed using QSAR, docking and other in silico approaches which can all be interfaced within the ChemScreener framework. As a sample application, in this work scaffold selectivity, diversity, connectivity and promiscuity towards six important therapeutic classes have been studied. In order to illustrate the computational power of the application, 55 scaffolds extracted from 161 anti-psychotic compounds were enumerated to produce a virtual library comprising 118 million compounds (17 GB) and annotated with chemophore, pharmacophore and toxicophore based features in a single step which would be non-trivial to perform with many standard software tools today on libraries of this size.

Design and Development of ChemInfoCloud: An Integrated Cloud Enabled Platform for Virtual Screening.

The power of cloud computing and distributed computing has been harnessed to handle vast and heterogeneous data required to be processed in any virtual screening protocol. A cloud computing platorm ChemInfoCloud was built and integrated with several chemoinformatics and bioinformatics tools. The robust engine performs the core chemoinformatics tasks of lead generation, lead optimisation and property prediction in a fast and efficient manner. It has also been provided with some of the bioinformatics functionalities including sequence alignment, active site pose prediction and protein ligand docking. Text mining, NMR chemical shift (1H, 13C) prediction and reaction fingerprint generation modules for efficient lead discovery are also implemented in this platform. We have developed an integrated problem solving cloud environment for virtual screening studies that also provides workflow management, better usability and interaction with end users using container based virtualization, OpenVz.

MegaMiner: A Tool for Lead Identification Through Text Mining Using Chemoinformatics Tools and Cloud Computing Environment.

Virtual screening is an indispensable tool to cope with the massive amount of data being tossed by the high throughput omics technologies. With the objective of enhancing the automation capability of virtual screening process a robust portal termed MegaMiner has been built using the cloud computing platform wherein the user submits a text query and directly accesses the proposed lead molecules along with their drug-like, lead-like and docking scores. Textual chemical structural data representation is fraught with ambiguity in the absence of a global identifier. We have used a combination of statistical models, chemical dictionary and regular expression for building a disease specific dictionary. To demonstrate the effectiveness of this approach, a case study on malaria has been carried out in the present work. MegaMiner offered superior results compared to other text mining search engines, as established by F score analysis. A single query term 'malaria' in the portlet led to retrieval of related PubMed records, protein classes, drug classes and 8000 scaffolds which were internally processed and filtered to suggest new molecules as potential anti-malarials. The results obtained were validated by docking the virtual molecules into relevant protein targets. It is hoped that MegaMiner will serve as an indispensable tool for not only identifying hidden relationships between various biological and chemical entities but also for building better corpus and ontologies.

Prediction of Bioactive Compounds Using Computed NMR Chemical Shifts.

NMR based chemical shifts are an important diagnostic parameter for structure elucidation as they capture rich information related to conformational, electronic and stereochemical arrangement of functional groups in a molecule which is responsible for its activity towards any biological target. The present work discusses the importance of computing NMR chemical shifts from molecular structures. The NMR chemical shift data (experimental or computed) was used to generate fingerprints in binary formats for mapping molecular fragments (as descriptors) and correlating with the bioactivity classes. For this study, chemical shift data derived binary fingerprints were computed for 149 classes and 4800 bioactive molecules. The sensitivity and selectivity of fingerprints in discriminating molecules belonging to different therapeutic categories was assessed using a LibSVM based classifier. An accuracy of 82% for proton and 94% for carbon NMR fingerprints were obtained for anti-psoriatic and anti-psychotic molecules demonstrating the effectiveness of this approach for virtual screening.

Pharmacophore and Docking Based Virtual Screening of Validated Mycobacterium tuberculosis Targets.

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.

Chemoinformatics Approach for Building Molecular Networks from Marine Organisms.

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.

Role of Chemical Reactivity and Transition State Modeling for Virtual Screening.

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.

A Study of Applications of Machine Learning Based Classification Methods for Virtual Screening of Lead Molecules.

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.

Chemical structure representations and applications in computational toxicity.

Efficient storage and retrieval of chemical structures is one of the most important prerequisite for solving any computational-based problem in life sciences. Several resources including research publications, text books, and articles are available on chemical structure representation. Chemical substances that have same molecular formula but several structural formulae, conformations, and skeleton framework/scaffold/functional groups of the molecule convey various characteristics of the molecule. Today with the aid of sophisticated mathematical models and informatics tools, it is possible to design a molecule of interest with specified characteristics based on their applications in pharmaceuticals, agrochemicals, biotechnology, nanomaterials, petrochemicals, and polymers. This chapter discusses both traditional and current state of art representation of chemical structures and their applications in chemical information management, bioactivity- and toxicity-based predictive studies.