Prediction of pathogenic single amino acid substitutions using molecular fragment descriptors.

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.

TCR-Pred: A new web-application for prediction of epitope and MHC specificity for CDR3 TCR sequences using molecular fragment descriptors.

The search for the relationships between CDR3 TCR sequences and epitopes or MHC types is a challenging task in modern immunology. We propose a new approach to develop the classification models of structure-activity relationships (SAR) using molecular fragment descriptors MNA (Multilevel Neighbourhoods of Atoms) to represent CDR3 TCR sequences and the naïve Bayes classifier algorithm. We have created the freely available TCR-Pred web application (http://way2drug.com/TCR-pred/) to predict the interactions between α chain CDR3 TCR sequences and 116 epitopes or 25 MHC types, as well as the interactions between ß chain CDR3 TCR sequences and 202 epitopes or 28 MHC types. The TCR-Pred web application is based on the data (more 250 000 unique CDR3 TCR sequences) from VDJdb, McPAS-TCR, and IEDB databases and the proposed approach. The average AUC values of the prediction accuracy calculated using a 20-fold cross-validation procedure varies from 0.857 to 0.884. The created web application may be useful in studies related with T-cell profiling based on CDR3 TCR sequences.

Phyto4Health: Database of Phytocomponents from Russian Pharmacopoeia Plants.

Medicinal plants growing in Russia are a rich source of biologically active compounds. However, the evaluation of the hidden pharmacological potential of these compounds by in silico methods is complicated by the lack of specialized databases. We have created a database of 3128 phytocomponents from 268 medical plants included in the Russian Pharmacopoeia. The information about the compounds was supplemented with their physical-chemical characteristics and biological activity profiles estimated using the PASS software. Comparison with phytocomponents of medicinal plants from five other countries showed that the similarity of phytocomponents in our database is rather small. The uniqueness of the contents significantly enriches and provides easy access to the necessary information. The Phyto4Health data are freely available at http://www.way2drug.com/p4h/.

HGMMX: Host Gut Microbiota Metabolism Xenobiotics Database.

The metagenome of bacteria colonizing the human intestine is a set of genes that is almost 150 times greater than the set of host genes. Some of these genes encode enzymes whose functioning significantly expands the number of potential pathways for xenobiotic metabolism. The resulting metabolites can exhibit activity different from that of the parent compound. This can decrease the efficacy of pharmacotherapy as well as induce undesirable and potentially life-threatening side effects. Thus, analysis of the biotransformation of small drug-like compounds mediated by the gut microbiota is an important step in the development of new pharmaceutical agents and repurposing of the approved drugs. In vitro research, the interaction of drug-like compounds with the gut microbiota is a multistep and time-consuming process. Systematic testing of large sets of chemical structures is associated with a number of challenges, including the lack of standardized techniques and significant financial costs to identify the structure of the final metabolites. Estimation of the compounds' ability to be biotransformed by the gut microbiota and prediction of the structures of their metabolites are possible in silico. However, the development of computational approaches is limited by the lack of information about chemical structures metabolized by microbiota enzymes. The aim of this study is to create a database containing information on the metabolism of drug-like compounds by the gut microbiota. We created the data set containing information about 368 structures metabolized and 310 structures not metabolized by the human gut microbiota. The HGMMX database is freely available at https://www.way2drug.com/hgmmx. The information presented will be useful in the development of computational approaches for analyzing the impact of the human microbiota on metabolism of drug-like molecules.

SAV-Pred: A Freely Available Web Application for the Prediction of Pathogenic Amino Acid Substitutions for Monogenic Hereditary Diseases Studied in Newborn Screening.

Next Generation Sequencing (NGS) technologies are rapidly entering clinical practice. A promising area for their use lies in the field of newborn screening. The mass screening of newborns using NGS technology leads to the discovery of a large number of new missense variants that need to be assessed for association with the development of hereditary diseases. Currently, the primary analysis and identification of pathogenic variations is carried out using bioinformatic tools. Although extensive efforts have been made in the computational approach to variant interpretation, there is currently no generally accepted pathogenicity predictor. In this study, we used the sequence-structure-property relationships (SSPR) approach, based on the representation of protein fragments by molecular structural formula. The approach predicts the pathogenic effect of single amino acid substitutions in proteins related with twenty-five monogenic heritable diseases from the Uniform Screening Panel for Major Conditions recommended by the Advisory Committee on Hereditary Disorders in Newborns and Children. In order to create SSPR models of classification, we modified a piece of cheminformatics software, MultiPASS, that was originally developed for the prediction of activity spectra for drug-like substances. The created SSPR models were compared with traditional bioinformatic tools (SIFT 4G, Polyphen-2 HDIV, MutationAssessor, PROVEAN and FATHMM). The average AUC of our approach was 0.804 ± 0.040. Better quality scores were achieved for 15 from 25 proteins with a significantly higher accuracy for some proteins (IVD, HADHB, HBB). The best SSPR models of classification are freely available in the online resource SAV-Pred (Single Amino acid Variants Predictor).

Identification of Molecular Mechanisms Involved in Viral Infection Progression Based on Text Mining: Case Study for HIV Infection.

Viruses cause various infections that may affect human lifestyle for durations ranging from several days to for many years. Although preventative and therapeutic remedies are available for many viruses, they may still have a profound impact on human life. The human immunodeficiency virus type 1 is the most common cause of HIV infection, which represents one of the most dangerous and complex diseases since it affects the immune system and causes its disruption, leading to secondary complications and negatively influencing health-related quality of life. While highly active antiretroviral therapy may decrease the viral load and the velocity of HIV infection progression, some individual peculiarities may affect viral load control or the progression of T-cell malfunction induced by HIV. Our study is aimed at the text-based identification of molecular mechanisms that may be involved in viral infection progression, using HIV as a case study. Specifically, we identified human proteins and genes which commonly occurred, overexpressed or underexpressed, in the collections of publications relevant to (i) HIV infection progression and (ii) acute and chronic stages of HIV infection. Then, we considered biological processes that are controlled by the identified protein and genes. We verified the impact of the identified molecules in the associated clinical study.

CLC-Pred 2.0: A Freely Available Web Application for In Silico Prediction of Human Cell Line Cytotoxicity and Molecular Mechanisms of Action for Druglike Compounds.

In vitro cell-line cytotoxicity is widely used in the experimental studies of potential antineoplastic agents and evaluation of safety in drug discovery. In silico estimation of cytotoxicity against hundreds of tumor cell lines and dozens of normal cell lines considerably reduces the time and costs of drug development and the assessment of new pharmaceutical agent perspectives. In 2018, we developed the first freely available web application (CLC-Pred) for the qualitative prediction of cytotoxicity against 278 tumor and 27 normal cell lines based on structural formulas of 59,882 compounds. Here, we present a new version of this web application: CLC-Pred 2.0. It also employs the PASS (Prediction of Activity Spectra for Substance) approach based on substructural atom centric MNA descriptors and a Bayesian algorithm. CLC-Pred 2.0 provides three types of qualitative prediction: (1) cytotoxicity against 391 tumor and 47 normal human cell lines based on ChEMBL and PubChem data (128,545 structures) with a mean accuracy of prediction (AUC), calculated by the leave-one-out (LOO CV) and the 20-fold cross-validation (20F CV) procedures, of 0.925 and 0.923, respectively; (2) cytotoxicity against an NCI60 tumor cell-line panel based on the Developmental Therapeutics Program's NCI60 data (22,726 structures) with different thresholds of IG50 data (100, 10 and 1 nM) and a mean accuracy of prediction from 0.870 to 0.945 (LOO CV) and from 0.869 to 0.942 (20F CV), respectively; (3) 2170 molecular mechanisms of actions based on ChEMBL and PubChem data (656,011 structures) with a mean accuracy of prediction 0.979 (LOO CV) and 0.978 (20F CV). Therefore, CLC-Pred 2.0 is a significant extension of the capabilities of the initial web application.

Relationships between the Structure and Severe Drug-Induced Liver Injury for Low, Medium, and High Doses of Drugs.

Assessment of structure-activity relationships (SARs) for predicting severe drug-induced liver injury (DILI) is essential since in vivo and in vitro preclinical methods cannot detect many druglike compounds disrupting liver functions. To date, plenty of SAR models for the prediction of DILI have been developed; however, none of them considered the route of drug administration and daily dose, which may introduce significant bias into prediction results. We have created a dataset of 617 drugs with parenteral and oral administration routes and consistent information on DILI severity. We have found a clear relationship between route, dose, and DILI severity. According to SAR, nearly 40% of moderate- and non-DILI-causing drugs would cause severe DILI if they were administered at high oral doses. We have proposed the following approach to predict severe DILI. New compounds recommended to be used at low oral doses (<∼10 mg daily), or parenterally, can be considered not causing severe DILI. DILI for compounds administered at medium oral doses (∼10-100 mg daily; 22.2% of drugs under consideration) can be considered unpredictable because reasonable SAR models were not obtained due to the small size and heterogeneity of the corresponding dataset. The DILI potential of the compounds recommended to be used at high oral doses (more than ∼100 mg daily) can be estimated using SAR modeling. The balanced accuracy of the approach calculated by a 10-fold cross-validation procedure is 0.803. The developed approach can be used to estimate severe DILI for druglike compounds proposed to use at low and high oral doses or parenterally at the early stages of drug development.

Computational Prediction of Inhibitors and Inducers of the Major Isoforms of Cytochrome P450.

Human cytochrome P450 enzymes (CYPs) are heme-containing monooxygenases. This superfamily of drug-metabolizing enzymes is responsible for the metabolism of most drugs and other xenobiotics. The inhibition of CYPs may lead to drug-drug interactions and impair the biotransformation of drugs. CYP inducers may decrease the bioavailability and increase the clearance of drugs. Based on the freely available databases ChEMBL and PubChem, we have collected over 70,000 records containing the structures of inhibitors and inducers together with the IC50 values for the inhibitors of the five major human CYPs: 1A2, 3A4, 2D6, 2C9, and 2C19. Based on the collected data, we developed (Q)SAR models for predicting inhibitors and inducers of these CYPs using GUSAR and PASS software. The developed (Q)SAR models could be applied for assessment of the interaction of novel drug-like substances with the major human CYPs. The created (Q)SAR models demonstrated reasonable accuracy of prediction. They have been implemented in the web application P450-Analyzer that is freely available via the Internet.

AntiHIV-Pred: web-resource for in silico prediction of anti-HIV/AIDS activity.

MOTIVATION: Identification of new molecules promising for treatment of HIV-infection and HIV-associated disorders remains an important task in order to provide safer and more effective therapies. Utilization of prior knowledge by application of computer-aided drug discovery approaches reduces time and financial expenses and increases the chances of positive results in anti-HIV R&D. To provide the scientific community with a tool that allows estimating of potential agents for treatment of HIV-infection and its comorbidities, we have created a freely-available web-resource for prediction of relevant biological activities based on the structural formulae of drug-like molecules. RESULTS: Over 50 000 experimental records for anti-retroviral agents from ChEMBL database were extracted for creating the training sets. After careful examination, about seven thousand molecules inhibiting five HIV-1 proteins were used to develop regression and classification models with the GUSAR software. The average values of R2 = 0.95 and Q2 = 0.72 in validation procedure demonstrated the reasonable accuracy and predictivity of the obtained (Q)SAR models. Prediction of 81 biological activities associated with the treatment of HIV-associated comorbidities with 92% mean accuracy was realized using the PASS program. AVAILABILITY AND IMPLEMENTATION: Freely available on the web at http://www.way2drug.com/hiv/. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Extraction of Data on Parent Compounds and Their Metabolites from Texts of Scientific Abstracts.

The growing amount of experimental data on chemical objects includes properties of small molecules, results of studies of their interaction with human and animal proteins, and methods of synthesis of organic compounds (OCs). The data obtained can be used to identify the names of OCs automatically, including all possible synonyms and relevant data on the molecular properties and biological activity. Utilization of different synonymic names of chemical compounds allows researchers to increase the completeness of data on their properties available from publications. Enrichment of the data on the names of chemical compounds by information about their possible metabolites can help estimate the biological effects of parent compounds and their metabolites more thoroughly. Therefore, an attempt at automated extraction of the names of parent compounds and their metabolites from the texts is a rather important task. In our study, we aimed at developing a method that provides the extraction of the named entities (NEs) of parent compounds and their metabolites from abstracts of scientific publications. Based on the application of the conditional random fields' algorithm, we extracted the NEs of chemical compounds. We developed a set of rules allowing identification of parent compound NEs and their metabolites in the texts. We evaluated the possibility of extracting the names of potential metabolites based on cosine similarity between strings representing names of parent compounds and all other chemical NEs found in the text. Additionally, we used conditional random fields to fetch the names of parent compounds and their metabolites from the texts based on the corpus of texts labeled manually. Our computational experiments showed that usage of rules in combination with cosine similarity could increase the accuracy of recognition of the names of metabolites compared to the rule-based algorithm and application of a machine-learning algorithm (conditional random fields).

QSAR without borders.

Prediction of chemical bioactivity and physical properties has been one of the most important applications of statistical and more recently, machine learning and artificial intelligence methods in chemical sciences. This field of research, broadly known as quantitative structure-activity relationships (QSAR) modeling, has developed many important algorithms and has found a broad range of applications in physical organic and medicinal chemistry in the past 55+ years. This Perspective summarizes recent technological advances in QSAR modeling but it also highlights the applicability of algorithms, modeling methods, and validation practices developed in QSAR to a wide range of research areas outside of traditional QSAR boundaries including synthesis planning, nanotechnology, materials science, biomaterials, and clinical informatics. As modern research methods generate rapidly increasing amounts of data, the knowledge of robust data-driven modelling methods professed within the QSAR field can become essential for scientists working both within and outside of chemical research. We hope that this contribution highlighting the generalizable components of QSAR modeling will serve to address this challenge.

Correction: QSAR without borders.

Correction for 'QSAR without borders' by Eugene N. Muratov et al., Chem. Soc. Rev., 2020, DOI: 10.1039/d0cs00098a.

Assessment of the cardiovascular adverse effects of drug-drug interactions through a combined analysis of spontaneous reports and predicted drug-target interactions.

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.

Monoclinic α-Na2FePO4F with Strong Antisite Disorder and Enhanced Na+ Diffusion.

A new monoclinic α-polymorph of the Na2FePO4F fluoride-phosphate has been directly synthesized via a hydrothermal method for application in metal-ion batteries. The crystal structure of the as-prepared α-Na2FePO4F studied with powder X-ray and neutron diffraction (P21/c, a = 13.6753(10) Å, b = 5.2503(2) Å, c = 13.7202(8) Å, ß = 120.230(4)°) demonstrates strong antisite disorder between the Na and Fe atoms. As revealed with DFT-based calculations, α-Na2FePO4F has low migration barriers for Na+ along the main pathway parallel to the b axis, and an additional diffusion bypass allowing the Na+ cations to go around the Na/Fe antisite defects. These results corroborate with the extremely high experimental Na-ion diffusion coefficient of (1-5)·10-11 cm2·s-1, which is 2 orders of magnitude higher than that for the orthorhombic ß-polymorph ((5-10)·10-13 cm2·s-1). Being tested as a cathode material in Na- and Li-ion battery cells, monoclinic α-Na2FePO4F exhibits a reversible specific capacity of 90 and 80 mAh g-1, respectively.

Prediction of Protein-ligand Interaction Based on Sequence Similarity and Ligand Structural Features.

Computationally predicting the interaction of proteins and ligands presents three main directions: the search of new target proteins for ligands, the search of new ligands for targets, and predicting the interaction of new proteins and new ligands. We proposed an approach providing the fuzzy classification of protein sequences based on the ligand structural features to analyze the latter most complicated case. We tested our approach on five protein groups, which represented promised targets for drug-like ligands and differed in functional peculiarities. The training sets were built with the original procedure overcoming the data ambiguity. Our study showed the effective prediction of new targets for ligands with an average accuracy of 0.96. The prediction of new ligands for targets displayed the average accuracy 0.95; accuracy estimates were close to our previous results, comparable in accuracy to those of other methods or exceeded them. Using the fuzzy coefficients reflecting the target-to-ligand specificity, we provided predicting interactions for new proteins and new ligands; the obtained accuracy values from 0.89 to 0.99 were acceptable for such a sophisticated task. The protein kinase family case demonstrated the ability to account for subtle features of proteins and ligands required for the specificity of protein-ligand interaction.

Computer-Aided Estimation of Biological Activity Profiles of Drug-Like Compounds Taking into Account Their Metabolism in Human Body.

Most pharmaceutical substances interact with several or even many molecular targets in the organism, determining the complex profiles of their biological activity. Moreover, due to biotransformation in the human body, they form one or several metabolites with different biological activity profiles. Therefore, the development and rational use of novel drugs requires the analysis of their biological activity profiles, taking into account metabolism in the human body. In silico methods are currently widely used for estimating new drug-like compounds' interactions with pharmacological targets and predicting their metabolic transformations. In this study, we consider the estimation of the biological activity profiles of organic compounds, taking into account the action of both the parent molecule and its metabolites in the human body. We used an external dataset that consists of 864 parent compounds with known metabolites. It is shown that the complex assessment of active pharmaceutical ingredients' interactions with the human organism increases the quality of computer-aided estimates. The toxic and adverse effects showed the most significant difference: reaching 0.16 for recall and 0.14 for precision.

A Computational Approach for the Prediction of Treatment History and the Effectiveness or Failure of Antiretroviral Therapy.

Human Immunodeficiency Virus Type 1 (HIV-1) infection is associated with high mortality if no therapy is provided. Currently, the treatment of an HIV-1 positive patient requires that several drugs should be taken simultaneously. The resistance of the virus to an antiretroviral drug may lead to treatment failure. Our approach focuses on predicting the exposure of a particular viral variant to an antiretroviral drug or drug combination. It also aims at the prediction of drug treatment success or failure. We utilized nucleotide sequences of HIV-1 encoding protease and reverse transcriptase to perform such types of prediction. The PASS (Prediction of Activity Spectra for Substances) algorithm based on the naive Bayesian classifier was used to make a prediction. We calculated the probability of whether a sequence belonged (P1) or did not belong (P0) to the class associated with exposure of the viral sequence to the set of drugs that can be associated with resistance to the set of drugs. The accuracy calculated as the average Area Under the ROC (Receiver Operating Characteristic) Curve (AUC/ROC) for classifying exposure of the sequence to the HIV-1 protease inhibitors was 0.81 (±0.07), and for HIV-1 reverse transcriptase, it was 0.83 (±0.07). To predict cases of treatment effectiveness or failure, we used P1 and P0 values, obtained in PASS, along with the binary vector constructed based on short nucleotide descriptors and the applied random forest classifier. Average AUC/ROC prediction accuracy for the prediction of treatment effectiveness or failure for the combinations of HIV-1 protease inhibitors was 0.82 (±0.06) and of HIV-1 reverse transcriptase was 0.76 (±0.09).

Data and Text Mining Help Identify Key Proteins Involved in the Molecular Mechanisms Shared by SARS-CoV-2 and HIV-1.

Viruses can be spread from one person to another; therefore, they may cause disorders in many people, sometimes leading to epidemics and even pandemics. New, previously unstudied viruses and some specific mutant or recombinant variants of known viruses constantly appear. An example is a variant of coronaviruses (CoV) causing severe acute respiratory syndrome (SARS), named SARS-CoV-2. Some antiviral drugs, such as remdesivir as well as antiretroviral drugs including darunavir, lopinavir, and ritonavir are suggested to be effective in treating disorders caused by SARS-CoV-2. There are data on the utilization of antiretroviral drugs against SARS-CoV-2. Since there are many studies aimed at the identification of the molecular mechanisms of human immunodeficiency virus type 1 (HIV-1) infection and the development of novel therapeutic approaches against HIV-1, we used HIV-1 for our case study to identify possible molecular pathways shared by SARS-CoV-2 and HIV-1. We applied a text and data mining workflow and identified a list of 46 targets, which can be essential for the development of infections caused by SARS-CoV-2 and HIV-1. We show that SARS-CoV-2 and HIV-1 share some molecular pathways involved in inflammation, immune response, cell cycle regulation.

ROSC-Pred: web-service for rodent organ-specific carcinogenicity prediction.

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: alexey.lagunin@ibmc.msk.ru. Supplementary information: Supplementary data are available at Bioinformatics online.