Peptide hemolytic activity analysis using visual data mining of similarity-based complex networks.

Peptides are promising drug development frameworks that have been hindered by intrinsic undesired properties including hemolytic activity. We aim to get a better insight into the chemical space of hemolytic peptides using a novel approach based on network science and data mining. Metadata networks (METNs) were useful to characterize and find general patterns associated with hemolytic peptides, whereas Half-Space Proximal Networks (HSPNs), represented the hemolytic peptide space. The best candidate HSPNs were used to extract various subsets of hemolytic peptides (scaffolds) considering network centrality and peptide similarity. These scaffolds have been proved to be useful in developing robust similarity-based model classifiers. Finally, using an alignment-free approach, we reported 47 putative hemolytic motifs, which can be used as toxic signatures when developing novel peptide-based drugs. We provided evidence that the number of hemolytic motifs in a sequence might be related to the likelihood of being hemolytic.

Unraveling druggable cancer-driving proteins and targeted drugs using artificial intelligence and multi-omics analyses.

The druggable proteome refers to proteins that can bind to small molecules with appropriate chemical affinity, inducing a favorable clinical response. Predicting druggable proteins through screening and in silico modeling is imperative for drug design. To contribute to this field, we developed an accurate predictive classifier for druggable cancer-driving proteins using amino acid composition descriptors of protein sequences and 13 machine learning linear and non-linear classifiers. The optimal classifier was achieved with the support vector machine method, utilizing 200 tri-amino acid composition descriptors. The high performance of the model is evident from an area under the receiver operating characteristics (AUROC) of 0.975 ± 0.003 and an accuracy of 0.929 ± 0.006 (threefold cross-validation). The machine learning prediction model was enhanced with multi-omics approaches, including the target-disease evidence score, the shortest pathways to cancer hallmarks, structure-based ligandability assessment, unfavorable prognostic protein analysis, and the oncogenic variome. Additionally, we performed a drug repurposing analysis to identify drugs with the highest affinity capable of targeting the best predicted proteins. As a result, we identified 79 key druggable cancer-driving proteins with the highest ligandability, and 23 of them demonstrated unfavorable prognostic significance across 16 TCGA PanCancer types: CDKN2A, BCL10, ACVR1, CASP8, JAG1, TSC1, NBN, PREX2, PPP2R1A, DNM2, VAV1, ASXL1, TPR, HRAS, BUB1B, ATG7, MARK3, SETD2, CCNE1, MUTYH, CDKN2C, RB1, and SMARCA4. Moreover, we prioritized 11 clinically relevant drugs targeting these proteins. This strategy effectively predicts and prioritizes biomarkers, therapeutic targets, and drugs for in-depth studies in clinical trials. Scripts are available at https://github.com/muntisa/machine-learning-for-druggable-proteins .

Multiquery Similarity Searching Models: An Alternative Approach for Predicting Hemolytic Activity from Peptide Sequence.

The desirable pharmacological properties and a broad number of therapeutic activities have made peptides promising drugs over small organic molecules and antibody drugs. Nevertheless, toxic effects, such as hemolysis, have hampered the development of such promising drugs. Hence, a reliable computational tool to predict peptide hemolytic toxicity is enormously useful before synthesis and experimental evaluation. Currently, four web servers that predict hemolytic activity using machine learning (ML) algorithms are available; however, they exhibit some limitations, such as the need for a reliable negative set and limited application domain. Hence, we developed a robust model based on a novel theoretical approach that combines network science and a multiquery similarity searching (MQSS) method. A total of 1152 initial models were constructed from 144 scaffolds generated in a previous report. These were evaluated on external data sets, and the best models were fused and improved. Our best MQSS model I1 outperformed all state-of-the-art ML-based models and was used to characterize the prevalence of hemolytic toxicity on therapeutic peptides. Based on our model's estimation, the number of hemolytic peptides might be 3.9-fold higher than the reported.

Exploring the dynamics of the ABCB1 membrane transporter P-glycoprotein in the presence of ATP and active/non-active compounds through molecular dynamics simulations.

P-glycoprotein (Pgp) is a member of the ATP-binding cassette family of transporters that confers multidrug resistance to cancer cells and is actively involved in the pharmacokinetics and toxicokinetics of a big variety of drugs. Extensive studies have provided insights into the binding of many compounds, but the precise mechanism of translocation across the membrane remains unknown; in this context, the major challenge has been to understand the basis for its polyspecificity. In this study, molecular dynamics (MD) simulations of human P-gp (hP-gp) in an explicit membrane-and-water environment were performed to investigate the dynamic behavior of the transporter in the presence of different compounds (active and inactive) in the binding pocket and ATP molecules within the nucleotide binding domains (NBDs). The complexes studied involve four compounds: cyclosporin A (CSA), amiodarone (AMI), pamidronate (APD), and valproic acid (VPA). While CSA and AMI are known to interact with P-gp, APD and VPA do not. The results highlighted how the presence of ATP notably contributed to increased flexibility of key residues in NBD1 of active systems, indicating potential conformational changes activating the translocation mechanism. MD simulations reveal how these domains adapt and respond to the presence of different substrates, as well as the influence of ATP binding on their flexibility. Furthermore, distinctive behavior was observed in the presence of active and inactive compounds, particularly in the arrangement of ATP between NBDs, supporting the proposed nucleotide sandwich dimer mechanism for ATP binding. This study provides comprehensive insights into P-gp behavior with various ligands and ATP, offering implications for drug development, toxicity assessment and demonstrating the validity of the results derived from the MD simulations.

Validation of venous thromboembolism predictive model in hematologic malignancies.

Although several scores stratify venous thromboembolism (VTE) risk in solid tumors, hematologic malignancies (HM) are underrepresented. To develop an internal and external validation of a logistic regression model to predict VTE risk in hospitalized HM patients. Validation of the existing VTE predictive model was performed through a prospective case-control study in 496 hospitalized HM patients between December 2010 and 2020 at the Arnaldo Milián University Hospital, Cuba. The predictive model designed with data from 285 patients includes 5 predictive factors: hypercholesterolemia, tumoral activity, use of thrombogenic drugs, diabetes mellitus, and immobilization. The model was internally validated using bootstrap analysis. External validation was realized in a prospective cohort of 211 HM patients. The predictive model had a 76.4% negative predictive value (NPV) and an 81.7% positive predictive value (PPV) in the bootstrapping validation. The area under curve (AUC) in the bootstrapping set was 0.838. Accuracy was 80.1% and 82.9% in the internal and external validation, respectively. In the external validation, the model produced 89.7% of NPV, 67.7% of PPV, 74.6% of sensitivity, and 86.2% of specificity. The AUC in the external validation was 0.900. VTE predictive model is a reproducible and simple tool with good accuracy and discrimination.

Molecular Modeling and In Vitro Evaluation of Piplartine Analogs against Oral Squamous Cell Carcinoma.

Cancer is a principal cause of death in the world, and providing a better quality of life and reducing mortality through effective pharmacological treatment remains a challenge. Among malignant tumor types, squamous cell carcinoma-esophageal cancer (EC) is usually located in the mouth, with approximately 90% located mainly on the tongue and floor of the mouth. Piplartine is an alkamide found in certain species of the genus Piper and presents many pharmacological properties including antitumor activity. In the present study, the cytotoxic potential of a collection of piplartine analogs against human oral SCC9 carcinoma cells was evaluated. The analogs were prepared via Fischer esterification reactions, alkyl and aryl halide esterification, and a coupling reaction with PyBOP using the natural compound 3,4,5-trimethoxybenzoic acid as a starting material. The products were structurally characterized using 1H and 13C nuclear magnetic resonance, infrared spectroscopy, and high-resolution mass spectrometry for the unpublished compounds. The compound 4-methoxy-benzyl 3,4,5-trimethoxybenzoate (9) presented an IC50 of 46.21 µM, high selectively (SI > 16), and caused apoptosis in SCC9 cancer cells. The molecular modeling study suggested a multi-target mechanism of action for the antitumor activity of compound 9 with CRM1 as the main target receptor.

Cytotoxic and Antifungal Amides Derived from Ferulic Acid: Molecular Docking and Mechanism of Action.

Amides derived from ferulic acid have a wide spectrum of pharmacological activities, including antitumor and antifungal activity. In the present study, a series of ten amides were obtained by coupling reactions using the reagents (benzotriazol-1-yloxy) tripyrrolidinophosphonium hexafluorophosphate (PyBOP) and N,N'-dicyclohexylcarbodiimide (DCC). All the compounds were identified on the basis of their IR, 1H- and 13C-NMR, HRMS data, and with yields ranging from 43.17% to 91.37%. The compounds were subjected to cytotoxic tests by the alamar blue technique and antifungal screening by the broth microdilution method to determine the minimum inhibitory concentration (MIC). The amides 10 and 11 displayed the best result in both biological evaluations, and compound 10 was the most potent and selective in HL-60 cancer cells, with no cytotoxicity on healthy cells. This amide had antifungal activity in all strains and had the lowest MIC against Candida albicans and Candida tropicalis. The possible mechanism of antifungal action occurs via the fungal cell wall. Molecular modeling suggested that compounds 10 and 11 interact with the enzymes GWT1 and GSC1, which are essential for the development of C. albicans. The findings of the present study demonstrated that compounds 10 and 11 may be used as a platform in drug development in the future.

Evolutionary algorithm-based generation of optimum peptide sequences with dengue virus inhibitory activity.

Background: There is currently no effective dengue virus (DENV) therapeutic. We aim to develop a genetic algorithm-based framework for the design of peptides with possible DENV inhibitory activity. Methods & results: A Python-based tool (denominated AutoPepGEN) based on a DENV support vector machine classifier as the objective function was implemented. AutoPepGEN was applied to the design of three- to seven-amino acid sequences and ten peptides were selected. Peptide-protease (DENV) docking and Molecular Mechanics-Generalized Born Surface Area calculations were performed for the selected sequences and favorable binding energies were observed. Conclusion: It is hoped that AutoPepGEN will serve as an in silico alternative to the experimental design of positional scanning combinatorial libraries, known to be prone to a combinatorial explosion. AutoPepGEN is available at: https://github.com/sjbarigye/AutoPepGEN.

A Multi-Objective Approach for Drug Repurposing in Preeclampsia.

Preeclampsia is a hypertensive disorder that occurs during pregnancy. It is a complex disease with unknown pathogenesis and the leading cause of fetal and maternal mortality during pregnancy. Using all drugs currently under clinical trial for preeclampsia, we extracted all their possible targets from the DrugBank and ChEMBL databases and labeled them as "targets". The proteins labeled as "off-targets" were extracted in the same way but while taking all antihypertensive drugs which are inhibitors of ACE and/or angiotensin receptor antagonist as query molecules. Classification models were obtained for each of the 55 total proteins (45 targets and 10 off-targets) using the TPOT pipeline optimization tool. The average accuracy of the models in predicting the external dataset for targets and off-targets was 0.830 and 0.850, respectively. The combinations of models maximizing their virtual screening performance were explored by combining the desirability function and genetic algorithms. The virtual screening performance metrics for the best model were: the Boltzmann-Enhanced Discrimination of ROC (BEDROC)α=160.9 = 0.258, the Enrichment Factor (EF)1% = 31.55 and the Area Under the Accumulation Curve (AUAC) = 0.831. The most relevant targets for preeclampsia were: AR, VDR, SLC6A2, NOS3 and CHRM4, while ABCG2, ERBB2, CES1 and REN led to the most relevant off-targets. A virtual screening of the DrugBank database identified estradiol, estriol, vitamins E and D, lynestrenol, mifrepristone, simvastatin, ambroxol, and some antibiotics and antiparasitics as drugs with potential application in the treatment of preeclampsia.

Structure-Function Relationships in the Human P-Glycoprotein (ABCB1): Insights from Molecular Dynamics Simulations.

P-Glycoprotein (P-gp) is a transmembrane protein belonging to the ATP binding cassette superfamily of transporters, and it is a xenobiotic efflux pump that limits intracellular drug accumulation by pumping compounds out of cells. P-gp contributes to a reduction in toxicity, and has broad substrate specificity. It is involved in the failure of many cancer and antiviral chemotherapies due to the phenomenon of multidrug resistance (MDR), in which the membrane transporter removes chemotherapeutic drugs from target cells. Understanding the details of the ligand-P-gp interaction is therefore critical for the development of drugs that can overcome the MDR phenomenon, for the early identification of P-gp substrates that will help us to obtain a more effective prediction of toxicity, and for the subsequent outdesign of substrate properties if needed. In this work, a series of molecular dynamics (MD) simulations of human P-gp (hP-gp) in an explicit membrane-and-water environment were performed to investigate the effects of binding different compounds on the conformational dynamics of P-gp. The results revealed significant differences in the behaviour of P-gp in the presence of active and non-active compounds within the binding pocket, as different patterns of movement were identified that could be correlated with conformational changes leading to the activation of the translocation mechanism. The predicted ligand-P-gp interactions are in good agreement with the available experimental data, as well as the estimation of the binding-free energies of the studied complexes, demonstrating the validity of the results derived from the MD simulations.

A Multi-Objective Approach for Anti-Osteosarcoma Cancer Agents Discovery through Drug Repurposing.

Osteosarcoma is the most common type of primary malignant bone tumor. Although nowadays 5-year survival rates can reach up to 60-70%, acute complications and late effects of osteosarcoma therapy are two of the limiting factors in treatments. We developed a multi-objective algorithm for the repurposing of new anti-osteosarcoma drugs, based on the modeling of molecules with described activity for HOS, MG63, SAOS2, and U2OS cell lines in the ChEMBL database. Several predictive models were obtained for each cell line and those with accuracy greater than 0.8 were integrated into a desirability function for the final multi-objective model. An exhaustive exploration of model combinations was carried out to obtain the best multi-objective model in virtual screening. For the top 1% of the screened list, the final model showed a BEDROC = 0.562, EF = 27.6, and AUC = 0.653. The repositioning was performed on 2218 molecules described in DrugBank. Within the top-ranked drugs, we found: temsirolimus, paclitaxel, sirolimus, everolimus, and cabazitaxel, which are antineoplastic drugs described in clinical trials for cancer in general. Interestingly, we found several broad-spectrum antibiotics and antiretroviral agents. This powerful model predicts several drugs that should be studied in depth to find new chemotherapy regimens and to propose new strategies for osteosarcoma treatment.

Perturbation-Theory Machine Learning (PTML) Multilabel Model of the ChEMBL Dataset of Preclinical Assays for Antisarcoma Compounds.

Sarcomas are a group of malignant neoplasms of connective tissue with a different etiology than carcinomas. The efforts to discover new drugs with antisarcoma activity have generated large datasets of multiple preclinical assays with different experimental conditions. For instance, the ChEMBL database contains outcomes of 37,919 different antisarcoma assays with 34,955 different chemical compounds. Furthermore, the experimental conditions reported in this dataset include 157 types of biological activity parameters, 36 drug targets, 43 cell lines, and 17 assay organisms. Considering this information, we propose combining perturbation theory (PT) principles with machine learning (ML) to develop a PTML model to predict antisarcoma compounds. PTML models use one function of reference that measures the probability of a drug being active under certain conditions (protein, cell line, organism, etc.). In this paper, we used a linear discriminant analysis and neural network to train and compare PT and non-PT models. All the explored models have an accuracy of 89.19-95.25% for training and 89.22-95.46% in validation sets. PTML-based strategies have similar accuracy but generate simplest models. Therefore, they may become a versatile tool for predicting antisarcoma compounds.

Prediction of breast cancer proteins involved in immunotherapy, metastasis, and RNA-binding using molecular descriptors and artificial neural networks.

Breast cancer (BC) is a heterogeneous disease where genomic alterations, protein expression deregulation, signaling pathway alterations, hormone disruption, ethnicity and environmental determinants are involved. Due to the complexity of BC, the prediction of proteins involved in this disease is a trending topic in drug design. This work is proposing accurate prediction classifier for BC proteins using six sets of protein sequence descriptors and 13 machine-learning methods. After using a univariate feature selection for the mix of five descriptor families, the best classifier was obtained using multilayer perceptron method (artificial neural network) and 300 features. The performance of the model is demonstrated by the area under the receiver operating characteristics (AUROC) of 0.980 ± 0.0037, and accuracy of 0.936 ± 0.0056 (3-fold cross-validation). Regarding the prediction of 4,504 cancer-associated proteins using this model, the best ranked cancer immunotherapy proteins related to BC were RPS27, SUPT4H1, CLPSL2, POLR2K, RPL38, AKT3, CDK3, RPS20, RASL11A and UBTD1; the best ranked metastasis driver proteins related to BC were S100A9, DDA1, TXN, PRNP, RPS27, S100A14, S100A7, MAPK1, AGR3 and NDUFA13; and the best ranked RNA-binding proteins related to BC were S100A9, TXN, RPS27L, RPS27, RPS27A, RPL38, MRPL54, PPAN, RPS20 and CSRP1. This powerful model predicts several BC-related proteins that should be deeply studied to find new biomarkers and better therapeutic targets. Scripts can be downloaded at https://github.com/muntisa/neural-networks-for-breast-cancer-proteins.

OncoOmics approaches to reveal essential genes in breast cancer: a panoramic view from pathogenesis to precision medicine.

Breast cancer (BC) is the leading cause of cancer-related death among women and the most commonly diagnosed cancer worldwide. Although in recent years large-scale efforts have focused on identifying new therapeutic targets, a better understanding of BC molecular processes is required. Here we focused on elucidating the molecular hallmarks of BC heterogeneity and the oncogenic mutations involved in precision medicine that remains poorly defined. To fill this gap, we established an OncoOmics strategy that consists of analyzing genomic alterations, signaling pathways, protein-protein interactome network, protein expression, dependency maps in cell lines and patient-derived xenografts in 230 previously prioritized genes to reveal essential genes in breast cancer. As results, the OncoOmics BC essential genes were rationally filtered to 140. mRNA up-regulation was the most prevalent genomic alteration. The most altered signaling pathways were associated with basal-like and Her2-enriched molecular subtypes. RAC1, AKT1, CCND1, PIK3CA, ERBB2, CDH1, MAPK14, TP53, MAPK1, SRC, RAC3, BCL2, CTNNB1, EGFR, CDK2, GRB2, MED1 and GATA3 were essential genes in at least three OncoOmics approaches. Drugs with the highest amount of clinical trials in phases 3 and 4 were paclitaxel, docetaxel, trastuzumab, tamoxifen and doxorubicin. Lastly, we collected ~3,500 somatic and germline oncogenic variants associated with 50 essential genes, which in turn had therapeutic connectivity with 73 drugs. In conclusion, the OncoOmics strategy reveals essential genes capable of accelerating the development of targeted therapies for precision oncology.

Gene Prioritization through Consensus Strategy, Enrichment Methodologies Analysis, and Networking for Osteosarcoma Pathogenesis.

Osteosarcoma is the most common subtype of primary bone cancer, affecting mostly adolescents. In recent years, several studies have focused on elucidating the molecular mechanisms of this sarcoma; however, its molecular etiology has still not been determined with precision. Therefore, we applied a consensus strategy with the use of several bioinformatics tools to prioritize genes involved in its pathogenesis. Subsequently, we assessed the physical interactions of the previously selected genes and applied a communality analysis to this protein-protein interaction network. The consensus strategy prioritized a total list of 553 genes. Our enrichment analysis validates several studies that describe the signaling pathways PI3K/AKT and MAPK/ERK as pathogenic. The gene ontology described TP53 as a principal signal transducer that chiefly mediates processes associated with cell cycle and DNA damage response It is interesting to note that the communality analysis clusters several members involved in metastasis events, such as MMP2 and MMP9, and genes associated with DNA repair complexes, like ATM, ATR, CHEK1, and RAD51. In this study, we have identified well-known pathogenic genes for osteosarcoma and prioritized genes that need to be further explored.

NFBTA: A Potent Cytotoxic Agent against Glioblastoma.

Piplartine (PPL), also known as piperlongumine, is a biologically active alkaloid extracted from the Piper genus which has been found to have highly effective anticancer activity against several tumor cell lines. This study investigates in detail the antitumoral potential of a PPL analogue; (E)-N-(4-fluorobenzyl)-3-(3,4,5-trimethoxyphenyl) acrylamide (NFBTA). The anticancer potential of NFBTA on the glioblastoma multiforme (GBM) cell line (U87MG) was determined by 3-(4,5-dimethyl-2-thia-zolyl)-2, 5-diphenyl-2H-tetrazolium bromide (MTT), and lactate dehydrogenase (LDH) release analysis, and the selectivity index (SI) was calculated. To detect cell apoptosis, fluorescent staining via flow cytometry and Hoechst 33258 staining were performed. Oxidative alterations were assessed via colorimetric measurement methods. Alterations in expressions of key genes related to carcinogenesis were determined. Additionally, in terms of NFBTA cytotoxic, oxidative, and genotoxic damage potential, the biosafety of this novel agent was evaluated in cultured human whole blood cells. Cell viability analyses revealed that NFBTA exhibited strong cytotoxic activity in cultured U87MG cells, with high selectivity and inhibitory activity in apoptotic processes, as well as potential for altering the principal molecular genetic responses in U87MG cell growth. Molecular docking studies strongly suggested a plausible anti-proliferative mechanism for NBFTA. The results of the experimental in vitro human glioblastoma model and computational approach revealed promising cytotoxic activity for NFBTA, helping to orient further studies evaluating its antitumor profile for safe and effective therapeutic applications.

Gene prioritization, communality analysis, networking and metabolic integrated pathway to better understand breast cancer pathogenesis.

Consensus strategy was proved to be highly efficient in the recognition of gene-disease association. Therefore, the main objective of this study was to apply theoretical approaches to explore genes and communities directly involved in breast cancer (BC) pathogenesis. We evaluated the consensus between 8 prioritization strategies for the early recognition of pathogenic genes. A communality analysis in the protein-protein interaction (PPi) network of previously selected genes was enriched with gene ontology, metabolic pathways, as well as oncogenomics validation with the OncoPPi and DRIVE projects. The consensus genes were rationally filtered to 1842 genes. The communality analysis showed an enrichment of 14 communities specially connected with ERBB, PI3K-AKT, mTOR, FOXO, p53, HIF-1, VEGF, MAPK and prolactin signaling pathways. Genes with highest ranking were TP53, ESR1, BRCA2, BRCA1 and ERBB2. Genes with highest connectivity degree were TP53, AKT1, SRC, CREBBP and EP300. The connectivity degree allowed to establish a significant correlation between the OncoPPi network and our BC integrated network conformed by 51 genes and 62 PPi. In addition, CCND1, RAD51, CDC42, YAP1 and RPA1 were functional genes with significant sensitivity score in BC cell lines. In conclusion, the consensus strategy identifies both well-known pathogenic genes and prioritized genes that need to be further explored.

A desirability-based multi objective approach for the virtual screening discovery of broad-spectrum anti-gastric cancer agents.

Gastric cancer is the third leading cause of cancer-related mortality worldwide and despite advances in prevention, diagnosis and therapy, it is still regarded as a global health concern. The efficacy of the therapies for gastric cancer is limited by a poor response to currently available therapeutic regimens. One of the reasons that may explain these poor clinical outcomes is the highly heterogeneous nature of this disease. In this sense, it is essential to discover new molecular agents capable of targeting various gastric cancer subtypes simultaneously. Here, we present a multi-objective approach for the ligand-based virtual screening discovery of chemical compounds simultaneously active against the gastric cancer cell lines AGS, NCI-N87 and SNU-1. The proposed approach relays in a novel methodology based on the development of ensemble models for the bioactivity prediction against each individual gastric cancer cell line. The methodology includes the aggregation of one ensemble per cell line using a desirability-based algorithm into virtual screening protocols. Our research leads to the proposal of a multi-targeted virtual screening protocol able to achieve high enrichment of known chemicals with anti-gastric cancer activity. Specifically, our results indicate that, using the proposed protocol, it is possible to retrieve almost 20 more times multi-targeted compounds in the first 1% of the ranked list than what is expected from a uniform distribution of the active ones in the virtual screening database. More importantly, the proposed protocol attains an outstanding initial enrichment of known multi-targeted anti-gastric cancer agents.

Chemoinformatics Profiling of the Chromone Nucleus as a MAO-B/A2AAR Dual Binding Scaffold.

BACKGROUND: In the context of the current drug discovery efforts to find disease modifying therapies for Parkinson's disease (PD) the current single target strategy has proved inefficient. Consequently, the search for multi-potent agents is attracting more and more attention due to the multiple pathogenetic factors implicated in PD. Multiple evidences points to the dual inhibition of the monoamine oxidase B (MAO-B), as well as adenosine A2A receptor (A2AAR) blockade, as a promising approach to prevent the neurodegeneration involved in PD. Currently, only two chemical scaffolds has been proposed as potential dual MAO-B inhibitors/A2AAR antagonists (caffeine derivatives and benzothiazinones). METHODS: In this study, we conduct a series of chemoinformatics analysis in order to evaluate and advance the potential of the chromone nucleus as a MAO-B/A2AAR dual binding scaffold. RESULTS: The information provided by SAR data mining analysis based on network similarity graphs and molecular docking studies support the suitability of the chromone nucleus as a potential MAOB/ A2AAR dual binding scaffold. Additionally, a virtual screening tool based on a group fusion similarity search approach was developed for the prioritization of potential MAO-B/A2AAR dual binder candidates. Among several data fusion schemes evaluated, the MEAN-SIM and MIN-RANK GFSS approaches demonstrated to be efficient virtual screening tools. Then, a combinatorial library potentially enriched with MAO-B/A2AAR dual binding chromone derivatives was assembled and sorted by using the MIN-RANK and then the MEAN-SIM GFSS VS approaches. CONCLUSION: The information and tools provided in this work represent valuable decision making elements in the search of novel chromone derivatives with a favorable dual binding profile as MAOB inhibitors and A2AAR antagonists with the potential to act as a disease-modifying therapeutic for Parkinson's disease.

Fusing Docking Scoring Functions Improves the Virtual Screening Performance for Discovering Parkinson's Disease Dual Target Ligands.

BACKGROUND: Virtual methodologies have become essential components of the drug discovery pipeline. Specifically, structure-based drug design methodologies exploit the 3D structure of molecular targets to discover new drug candidates through molecular docking. Recently, dual target ligands of the Adenosine A2A Receptor and Monoamine Oxidase B enzyme have been proposed as effective therapies for the treatment of Parkinson's disease. METHODS: In this paper we propose a structure-based methodology, which is extensively validated, for the discovery of dual Adenosine A2A Receptor/Monoamine Oxidase B ligands. This methodology involves molecular docking studies against both receptors and the evaluation of different scoring functions fusion strategies for maximizing the initial virtual screening enrichment of known dual ligands. RESULTS: The developed methodology provides high values of enrichment of known ligands, which outperform that of the individual scoring functions. At the same time, the obtained ensemble can be translated in a sequence of steps that should be followed to maximize the enrichment of dual target Adenosine A2A Receptor antagonists and Monoamine Oxidase B inhibitors. CONCLUSION: Information relative to docking scores to both targets have to be combined for achieving high dual ligands enrichment. Combining the rankings derived from different scoring functions proved to be a valuable strategy for improving the enrichment relative to single scoring function in virtual screening experiments.