MOTIVATION: Next Generation Sequencing technologies make it possible to detect rare genetic variants in individual patients. Currently, more than a dozen software and web services have been created to predict the pathogenicity of variants related with changing of amino acid residues. Despite considerable efforts in this area, at the moment there is no ideal method to classify pathogenic and harmless variants, and the assessment of the pathogenicity is often contradictory. In this article, we propose to use peptides structural formulas of proteins as an amino acid residues substitutions description, rather than a single-letter code. This allowed us to investigate the effectiveness of chemoinformatics approach to assess the pathogenicity of variants associated with amino acid substitutions. RESULTS: The structure-activity relationships analysis relying on protein-specific data and atom centric substructural multilevel neighborhoods of atoms (MNA) descriptors of molecular fragments appeared to be suitable for predicting the pathogenic effect of single amino acid variants. MNA-based Naïve Bayes classifier algorithm, ClinVar and humsavar data were used for the creation of structure-activity relationships models for 10 proteins. The performance of the models was compared with 11 different predicting tools: 8 individual (SIFT 4G, Polyphen2 HDIV, MutationAssessor, PROVEAN, FATHMM, MVP, LIST-S2, MutPred) and 3 consensus (M-CAP, MetaSVM, MetaLR). The accuracy of MNA-based method varies for the proteins (AUC: 0.631-0.993; MCC: 0.191-0.891). It was similar for both the results of comparisons with the other individual predictors and third-party protein-specific predictors. For several proteins (BRCA1, BRCA2, COL1A2, and RYR1), the performance of the MNA-based method was outstanding, capable of capturing the pathogenic effect of structural changes in amino acid substitutions. AVAILABILITY AND IMPLEMENTATION: The datasets are available as supplemental data at Bioinformatics online. A python script to convert amino acid and nucleotide sequences from single-letter codes to SD files is available at https://github.com/SmirnygaTotoshka/SequenceToSDF. The authors provide trial licenses for MultiPASS software to interested readers upon request.

Adverse drug effects (ADEs) are one of the leading causes of death in developed countries and are the main reason for drug recalls from the market, whereas the ADEs that are associated with action on the cardiovascular system are the most dangerous and widespread. The treatment of human diseases often requires the intake of several drugs, which can lead to undesirable drug-drug interactions (DDIs), thus causing an increase in the frequency and severity of ADEs. An evaluation of DDI-induced ADEs is a nontrivial task and requires numerous experimental and clinical studies. Therefore, we developed a computational approach to assess the cardiovascular ADEs of DDIs. This approach is based on the combined analysis of spontaneous reports (SRs) and predicted drug-target interactions to estimate the five cardiovascular ADEs that are induced by DDIs, namely, myocardial infarction, ischemic stroke, ventricular tachycardia, cardiac failure, and arterial hypertension. We applied a method based on least absolute shrinkage and selection operator (LASSO) logistic regression to SRs for the identification of interacting pairs of drugs causing corresponding ADEs, as well as noninteracting pairs of drugs. As a result, five datasets containing, on average, 3100 potentially ADE-causing and non-ADE-causing drug pairs were created. The obtained data, along with information on the interaction of drugs with 1553 human targets predicted by PASS Targets software, were used to create five classification models using the Random Forest method. The average area under the ROC curve of the obtained models, sensitivity, specificity and balanced accuracy were 0.837, 0.764, 0.754 and 0.759, respectively. The predicted drug targets were also used to hypothesize the potential mechanisms of DDI-induced ventricular tachycardia for the top-scoring drug pairs. The created five classification models can be used for the identification of drug combinations that are potentially the most or least dangerous for the cardiovascular system.

Motivation: Identification of rodent carcinogens is an important task in risk assessment of chemicals. SAR methods were proposed to reduce the number of animal experiments. Most of these methods ignore information about organ-specificity of tumorigenesis. Our study was aimed at the creation of classification models and a freely available online service for prediction of rodent carcinogens considering the species (rats, mice), sex and tissue-specificity from structural formula of compounds. Results: The data from Carcinogenic Potency Database for 1011 organic compounds evaluated on the standard two-year rodent carcinogenicity bioassay was used for the creation of training sets. Structure-activity relationships models for prediction of rodent organ-specific carcinogenicity were created by PASS software, which was based on Bayesian-like approach and Multilevel Neighborhoods of Atoms descriptors. The average prediction accuracy for training sets calculated by leave-one-out and 10-fold cross-validation was 79 and 78.2%, respectively. Availability and implementation: Freely available on the web at http://www.way2drug.com/ROSC. Contact: [email protected]. Supplementary information: Supplementary data are available at Bioinformatics online.

UNLABELLED: A new freely available web server site of metabolism predictor to predict the sites of metabolism (SOM) based on the structural formula of chemicals has been developed. It is based on the analyses of 'structure-SOM' relationships using a Bayesian approach and labelled multilevel neighbourhoods of atoms descriptors to represent the structures of over 1000 metabolized xenobiotics. The server allows predicting SOMs that are catalysed by 1A2, 2C9, 2C19, 2D6 and 3A4 isoforms of cytochrome P450 and enzymes of the UDP-glucuronosyltransferase family. The average invariant accuracy of prediction that was calculated for the training sets (using leave-one-out cross-validation) and evaluation sets is 0.9 and 0.95, respectively. AVAILABILITY AND IMPLEMENTATION: Freely available on the web at http://www.way2drug.com/SOMP.

SUMMARY: Experimentally found gene expression profiles are used to solve different problems in pharmaceutical studies, such as drug repositioning, resistance, toxicity and drug-drug interactions. A special web service, DIGEP-Pred, for prediction of drug-induced changes of gene expression profiles based on structural formulae of chemicals has been developed. Structure-activity relationships for prediction of drug-induced gene expression profiles were determined by Prediction of Activity Spectra for Substances (PASS) software. Comparative Toxicogenomics Database with data on the known drug-induced gene expression profiles of chemicals was used to create mRNA- and protein-based training sets. An average prediction accuracy for the training sets (ROC AUC) calculated by leave-one-out cross-validation on the basis of mRNA data (1385 compounds, 952 genes, 500 up- and 475 down-regulations) and protein data (1451 compounds, 139 genes, 93 up- and 55 down-regulations) exceeded 0.85. AVAILABILITY: Freely available on the web at http://www.way2drug.com/GE.

BACKGROUND: The knowledge about proteins with specific interaction capacity to the protein partners is very important for the modeling of cell signaling networks. However, the experimentally-derived data are sufficiently not complete for the reconstruction of signaling pathways. This problem can be solved by the network enrichment with predicted protein interactions. The previously published in silico method PAAS was applied for prediction of interactions between protein kinases and their substrates. RESULTS: We used the method for recognition of the protein classes defined by the interaction with the same protein partners. 1021 protein kinase substrates classified by 45 kinases were extracted from the Phospho.ELM database and used as a training set. The reasonable accuracy of prediction calculated by leave-one-out cross validation procedure was observed in the majority of kinase-specificity classes. The random multiple splitting of the studied set onto the test and training set had also led to satisfactory results. The kinase substrate specificity for 186 proteins extracted from TRANSPATH database was predicted by PAAS method. Several kinase-substrate interactions described in this database were correctly predicted. Using the previously developed ExPlain system for the reconstruction of signal transduction pathways, we showed that addition of the newly predicted interactions enabled us to find the possible path between signal trigger, TNF-alpha, and its target genes in the cell. CONCLUSIONS: It was shown that the predictions of protein kinase substrates by PAAS were suitable for the enrichment of signaling pathway networks and identification of the novel signaling pathways. The on-line version of PAAS for prediction of protein kinase substrates is freely available at http://www.ibmc.msk.ru/PAAS/.

The concept of the biological activity spectrum was introduced to describe the properties of biologically active substances. The PASS (prediction of activity spectra for substances) software product, which predicts more than 300 pharmacological effects and biochemical mechanisms on the basis of the structural formula of a substance, may be efficiently used to find new targets (mechanisms) for some ligands and, conversely, to reveal new ligands for some biological targets. We have developed a WWW interface for the PASS software. A WWW server for the on-line prediction of the biological activity spectra of substances has been constructed.