Fragment-based in silico modeling of multi-target inhibitors against breast cancer-related proteins

Modeling and Interpretability Study of the Structure-Activity Relationship for Multigeneration EGFR Inhibitors

The fourth-generation EGFR inhibitors targeting L858R/T790M/C797S mutations are in clinical trials mostly, and it is necessary to develop new inhibitors. In this study, an internal data set containing 2302 multitarget EGFR inhibitors targeting the wild type (83%) and the L858R (92%), L858R/T790M (96%), and L858R/T790M/C797S (60%) mutations was collected. We established a structure-activity relationship model for predicting the bioactivities of multigeneration EGFR inhibitors by a multitask deep neural network (MT-DNN). We also constructed four single-task models on 1384 L858R/T790M/C797S (60%) mutation inhibitors by support vector machine (SVM), random forest (RF), XGBoost (XGB), and single-target neural network (ST-DNN), respectively. The MT-DNN model significantly outperformed single-task models on the external data set of 304 fourth-generation EGFR inhibitors. Furthermore, the combined application of MT-DNN and SHAP/delta-SHAP value interpretability analysis offers rigorous structural information from a global perspective. With SHAP/delta-SHAP methods, the MT-DNN model can mine the core scaffold and important fragments of multigeneration EGFR inhibitors and provide valuable information from a structure-activity relationship perspective to address the resistant mutation problem.

Computer-aided Drug Discovery Approaches in the Identification of Anticancer Drugs from Natural Products: A Review

Natural plant sources are essential in the development of several anticancer drugs, such as vincristine, vinblastine, vinorelbine, docetaxel, paclitaxel, camptothecin, etoposide, and teniposide. However, various chemotherapies fail due to adverse reactions, drug resistance, and target specificity. Researchers are now focusing on developing drugs that use natural compounds to overcome these issues. These drugs can affect multiple targets, have reduced adverse effects, and are effective against several cancer types. Developing a new drug is a highly complex, expensive, and time-consuming process. Traditional drug discovery methods take up to 15 years for a new medicine to enter the market and cost more than one billion USD. However, recent Computer Aided Drug Discovery (CADD) advancements have changed this situation. This paper aims to comprehensively describe the different CADD approaches in identifying anticancer drugs from natural products. Data from various sources, including Science Direct, Elsevier, NCBI, and Web of Science, are used in this review. In-silico techniques and optimization algorithms can provide versatile solutions in drug discovery ventures. The structure-based drug design technique is widely used to understand chemical constituents' molecular-level interactions and identify hit leads. This review will discuss the concept of CADD, in-silico tools, virtual screening in drug discovery, and the concept of natural products as anticancer therapies. Representative examples of molecules identified will also be provided.

AISMPred: A Machine Learning Approach for Predicting Anti-Inflammatory Small Molecules

Background/Objectives: Inflammation serves as a vital response to diverse harmful stimuli like infections, toxins, or tissue injuries, aiding in the elimination of pathogens and tissue repair. However, persistent inflammation can lead to chronic diseases. Peptide therapeutics have gained attention for their specificity in targeting cells, yet their development remains costly and time-consuming. Therefore, small molecules, with their stability, low immunogenicity, and oral bioavailability, have become a focal point for predicting anti-inflammatory small molecules (AISMs). Methods: In this study, we introduce a computational method called AISMPred, designed to classify AISMs and non-AISMs. To develop this approach, we constructed a dataset comprising 1750 AISMs and non-AISMs, each annotated with IC50 values sourced from the PubChem BioAssay database. We computed two distinct types of molecular descriptors using PaDEL and Mordred tools. Subsequently, these descriptors were concatenated to form a hybrid feature set. The SVC-L1 regularization method was implemented for the optimum feature selection to develop robust Machine learning (ML) models. Five different conventional ML classifiers were employed, such as RF, ET, KNN, LR, and Ensemble methods. Results: A total of 15 ML models were developed using 2D, FP, and Hybrid feature sets, with the ET model with hybrid features achieving the highest accuracy of 92% and an AUC of 0.97 on the independent test dataset. Conclusions: This study provides an effective method for screening AISMs, potentially impacting drug discovery and design.

Breakthroughs in AI and multi-omics for cancer drug discovery: A review

Cancer is one of the biggest medical challenges we face today. It is characterized by abnormal, uncontrolled growth of cells that can spread to different parts of the body. Cancer is extremely complex, with genetic variations and the ability to adapt and evolve. This means we must continuously pursue innovative approaches to developing new cancer drugs. While traditional drug discovery methods have led to important breakthroughs, they also have significant limitations that make it difficult to efficiently create new, cost-effective cancer therapies. Integrating computational tools into the cancer drug discovery process is a major step forward. By harnessing computing power, we can overcome some of the inherent barriers of traditional methods. This review examines the range of computational techniques now being used, such as molecular docking, QSAR models, virtual screening, and pharmacophore modeling. It looks at recent advances in areas like machine learning and molecular simulations. The review also discusses the current challenges with these technologies and envisions future directions, underscoring how transformative these computational tools can be for creating targeted, new cancer treatments.

Design and Optimization of Quinazoline Derivatives as Potent EGFR Inhibitors for Lung Cancer Treatment: A Comprehensive QSAR, ADMET, and Molecular Modeling Investigation

The epidermal growth factor receptor (EGFR) is part of a protein family that controls cell growth and development. Due to its importance, EGFR has been identified as a suitable target for creating novel drugs. For this research, we conducted a 2D-QSAR analysis on a set of 31 molecules derived from quinazoline, which exhibited inhibitory activity against human lung cancer. This investigation incorporated principal component analysis (PCA) and multiple linear regression (MLR), leading to the development of QSAR models with strong predictive capabilities (R 2 = 0.745, R 2_adj = 0.723, MSE = 0.061, R 2_test = 0.941, and Q 2_cv = 0.669). The reliability of these models was confirmed through internal, external, Y-randomization, and applicability domain validations. Leveraging the predictions from the QSAR model, we designed 18 new molecules based on the modifications at the N-3 and C-6 positions of the quinazoline ring, with electronegative substituents at these positions fostering optimal polar interactions and hydrophobic contacts within the ATP-binding site of EGFR, significantly enhancing the inhibitory activity against the lung cancer cell line. Subsequently, ADMET predictions were conducted for these 18 compounds, revealing outstanding ADMET profiles. Molecular docking analyses were performed to investigate the interactions between the newly designed molecules-Pred15, Pred17, Pred20, Pred21-and the EGFR protein, indicating high affinity of these proposed compounds to EGFR. Furthermore, molecular dynamics (MD) simulations were utilized to assess the stability and binding modes of compounds Pred17, Pred20, and Pred21, reinforcing their potential as novel inhibitors against human lung cancer. Overall, our findings suggest that these investigated compounds can serve as effective inhibitors, showcasing the utility of our analytical and design approach in the identification of promising therapeutic agents.

Improved QSAR models for PARP-1 inhibition using data balancing, interpretable machine learning, and matched molecular pair analysis

The poly (ADP-ribose) polymerase-1 (PARP-1) enzyme is an important target in the treatment of breast cancer. Currently, treatment options include the drugs Olaparib, Niraparib, Rucaparib, and Talazoparib; however, these drugs can cause severe side effects including hematological toxicity and cardiotoxicity. Although in silico models for the prediction of PARP-1 activity have been developed, the drawbacks of these models include low specificity, a narrow applicability domain, and a lack of interpretability. To address these issues, a comprehensive machine learning (ML)-based quantitative structure-activity relationship (QSAR) approach for the informed prediction of PARP-1 activity is presented. Classification models built using the Synthetic Minority Oversampling Technique (SMOTE) for data balancing gave robust and predictive models based on the K-nearest neighbor algorithm (accuracy 0.86, sensitivity 0.88, specificity 0.80). Regression models were built on structurally congeneric datasets, with the models for the phthalazinone class and fused cyclic compounds giving the best performance. In accordance with the Organization for Economic Cooperation and Development (OECD) guidelines, a mechanistic interpretation is proposed using the Shapley Additive Explanations (SHAP) to identify the important topological features to differentiate between PARP-1 actives and inactives. Moreover, an analysis of the PARP-1 dataset revealed the prevalence of activity cliffs, which possibly negatively impacts the model's predictive performance. Finally, a set of chemical transformation rules were extracted using the matched molecular pair analysis (MMPA) which provided mechanistic insights and can guide medicinal chemists in the design of novel PARP-1 inhibitors.

Support Vector Machine-Based Prediction Models for Drug Repurposing and Designing Novel Drugs for Colorectal Cancer

Colorectal cancer (CRC) has witnessed a concerning increase in incidence and poses a significant therapeutic challenge due to its poor prognosis. There is a pressing demand to identify novel drug therapies to combat CRC. In this study, we addressed this need by utilizing the pharmacological profiles of anticancer drugs from the Genomics of Drug Sensitivity in Cancer (GDSC) database and developed QSAR models using the Support Vector Machine (SVM) algorithm for prediction of alternative and promiscuous anticancer compounds for CRC treatment. Our QSAR models demonstrated their robustness by achieving a high correlation of determination (R2) after 10-fold cross-validation. For 12 CRC cell lines, R2 ranged from 0.609 to 0.827. The highest performance was achieved for SW1417 and GP5d cell lines with R2 values of 0.827 and 0.786, respectively. Further, we listed the most common chemical descriptors in the drug profiles of the CRC cell lines and we also further reported the correlation of these descriptors with drug activity. The KRFP314 fingerprint was the predominantly occurring descriptor, with the KRFPC314 fingerprint following closely in prevalence within the drug profiles of the CRC cell lines. Beyond predictive modeling, we also confirmed the applicability of our developed QSAR models via in silico methods by conducting descriptor-drug analyses and recapitulating drug-to-oncogene relationships. We also identified two potential anti-CRC FDA-approved drugs, viomycin and diamorphine, using QSAR models. To ensure the easy accessibility and utility of our research findings, we have incorporated these models into a user-friendly prediction Web server named "ColoRecPred", available at https://project.iith.ac.in/cgntlab/colorecpred. We anticipate that this Web server can be used for screening of chemical libraries to identify potential anti-CRC drugs.

VGAE-MCTS: A New Molecular Generative Model Combining the Variational Graph Auto-Encoder and Monte Carlo Tree Search

Molecular generation is crucial for advancing drug discovery, materials science, and chemical exploration. It expedites the search for new drug candidates, facilitates tailored material creation, and enhances our understanding of molecular diversity. By employing artificial intelligence techniques such as molecular generative models based on molecular graphs, researchers have tackled the challenge of identifying efficient molecules with desired properties. Here, we propose a new molecular generative model combining a graph-based deep neural network and a reinforcement learning technique. We evaluated the validity, novelty, and optimized physicochemical properties of the generated molecules. Importantly, the model explored uncharted regions of chemical space, allowing for the efficient discovery and design of new molecules. This innovative approach has considerable potential to revolutionize drug discovery, materials science, and chemical research for accelerating scientific innovation. By leveraging advanced techniques and exploring previously unexplored chemical spaces, this study offers promising prospects for the efficient discovery and design of new molecules in the field of drug development.

Optimizing drug discovery using multitasking models for quantitative structure-biological effect relationships: an update of the literature

INTRODUCTION

Drug discovery has provided modern societies with the means to fight against many diseases. In this sense, computational methods have been at the forefront, playing an important role in rationalizing the search for novel drugs. Yet, tackling phenomena such as the multi-genic nature of diseases and drug resistance are limitations of the current computational methods. Multi-tasking models for quantitative structure-biological effect relationships (mtk-QSBER) have emerged to overcome such limitations.

AREAS COVERED

The present review describes an update on the fundamentals and applications of the mtk-QSBER models as tools to accelerate multiple stages/substages of the drug discovery process.

EXPERT OPINION

Computational approaches are extremely important for the rationalization of the search for novel and efficacious therapeutic agents. However, they need to focus more on the multi-target drug discovery paradigm. In this sense, mtk-QSBER models are particularly suited for multi-target drug discovery, offering encouraging opportunities across multiple therapeutic areas and scientific disciplines associated with drug discovery.

Nano-QSAR modeling for predicting the cytotoxicity of metallic and metal oxide nanoparticles: A review

Given the rapid development of nanotechnology, it is crucial to understand the effects of nanoparticles on living organisms. However, it is laborious to perform toxicological tests on a case-by-case basis. Quantitative structure-activity relationship (QSAR) is an effective computational technique because it saves time, costs, and animal sacrifice. Therefore, this review presents general procedures for the construction and application of nano-QSAR models of metal-based and metal-oxide nanoparticles (MBNPs and MONPs). We also provide an overview of available databases and common algorithms. The molecular descriptors and their roles in the toxicological interpretation of MBNPs and MONPs are systematically reviewed and the future of nano-QSAR is discussed. Finally, we address the growing demand for novel nano-specific descriptors, new computational strategies to address the data shortage, in situ data for regulatory concerns, a better understanding of the physicochemical properties of NPs with bioactivity, and, most importantly, the design of nano-QSAR for real-life environmental predictions rather than laboratory simulations.

Artificial intelligence and machine-learning approaches in structure and ligand-based discovery of drugs affecting central nervous system

CNS disorders are indications with a very high unmet medical needs, relatively smaller number of available drugs, and a subpar satisfaction level among patients and caregiver. Discovery of CNS drugs is extremely expensive affair with its own unique challenges leading to extremely high attrition rates and low efficiency. With explosion of data in information age, there is hardly any aspect of life that has not been touched by data driven technologies such as artificial intelligence (AI) and machine learning (ML). Drug discovery is no exception, emergence of big data via genomic, proteomic, biological, and chemical technologies has driven pharmaceutical giants to collaborate with AI oriented companies to revolutionise drug discovery, with the goal of increasing the efficiency of the process. In recent years many examples of innovative applications of AI and ML techniques in CNS drug discovery has been reported. Research on therapeutics for diseases such as schizophrenia, Alzheimer's and Parkinsonism has been provided with a new direction and thrust from these developments. AI and ML has been applied to both ligand-based and structure-based drug discovery and design of CNS therapeutics. In this review, we have summarised the general aspects of AI and ML from the perspective of drug discovery followed by a comprehensive coverage of the recent developments in the applications of AI/ML techniques in CNS drug discovery.

In Silico Drug Repurposing for Anti-Inflammatory Therapy: Virtual Search for Dual Inhibitors of Caspase-1 and TNF-Alpha

Inflammation involves a complex biological response of the body tissues to damaging stimuli. When dysregulated, inflammation led by biomolecular mediators such as caspase-1 and tumor necrosis factor-alpha (TNF-alpha) can play a detrimental role in the progression of different medical conditions such as cancer, neurological disorders, autoimmune diseases, and cytokine storms caused by viral infections such as COVID-19. Computational approaches can accelerate the search for dual-target drugs able to simultaneously inhibit the aforementioned proteins, enabling the discovery of wide-spectrum anti-inflammatory agents. This work reports the first multicondition model based on quantitative structure–activity relationships and a multilayer perceptron neural network (mtc-QSAR-MLP) for the virtual screening of agency-regulated chemicals as versatile anti-inflammatory therapeutics. The mtc-QSAR-MLP model displayed accuracy higher than 88%, and was interpreted from a physicochemical and structural point of view. When using the mtc-QSAR-MLP model as a virtual screening tool, we could identify several agency-regulated chemicals as dual inhibitors of caspase-1 and TNF-alpha, and the experimental information later retrieved from the scientific literature converged with our computational results. This study supports the capabilities of our mtc-QSAR-MLP model in anti-inflammatory therapy with direct applications to current health issues such as the COVID-19 pandemic.

Artificial intelligence for assisting cancer diagnosis and treatment in the era of precision medicine

Over the past decade, artificial intelligence (AI) has contributed substantially to the resolution of various medical problems, including cancer. Deep learning (DL), a subfield of AI, is characterized by its ability to perform automated feature extraction and has great power in the assimilation and evaluation of large amounts of complicated data. On the basis of a large quantity of medical data and novel computational technologies, AI, especially DL, has been applied in various aspects of oncology research and has the potential to enhance cancer diagnosis and treatment. These applications range from early cancer detection, diagnosis, classification and grading, molecular characterization of tumors, prediction of patient outcomes and treatment responses, personalized treatment, automatic radiotherapy workflows, novel anti-cancer drug discovery, and clinical trials. In this review, we introduced the general principle of AI, summarized major areas of its application for cancer diagnosis and treatment, and discussed its future directions and remaining challenges. As the adoption of AI in clinical use is increasing, we anticipate the arrival of AI-powered cancer care.

Part-II- in silico drug design: application and success

In silico tools have indeed reframed the steps involved in traditional drug discovery and development process and the term in silico has become a familiar term in pharmaceutical sector like the terms in vitro and in vivo. The successful design of HIV protease inhibitors, Saquinavir, Indinavir and other important medicinal agents, initiated interest of researchers in structure based drug design approaches (SBDD). The interactions between biomolecules and a ligand, binding energy, free energy and stability of biomolecule-ligand complex can be envisioned and predicted by applying molecular docking studies. Protein-ligand, protein-protein, DNA-ligand interactions etc. aid in elucidating molecular level mechanisms of drug molecules. In the Ligand based drug design (LBDD) approaches, QSAR studies have tremendously contributed to the development of antimicrobial, anticancer, antimalarial agents. In the recent years, multiQSAR (mt-QSAR) approaches have been successfully employed for designing drugs against multifactorial diseases. Output of a research in several instances is rewarding when both SBDD and LBDD approaches are combined. Application of in silico studies for prediction of pharmacokinetics was once a real challenge but one can see unlimited number publications comprising tools, data bases which can accurately predict almost all the pharmacokinetic parameters. Absorption, distribution, metabolism, transporters, blood brain barrier permeability, hERG toxicity, P-gp affinity and several toxicological end points can be accurately predicted for a candidate molecule before its synthesis. In silico approaches are greatly encouraged a result of growing limitations and new legislations related to the animal use for research. The combined use of in vitro data and in silico tools will definitely decrease the use of animal testing in the future.In this chapter, in silico approaches and their applications are reviewed and discussed giving suitable examples.

Structure based designing of thiazolidinone-pyrimidine derivatives as ERK2 inhibitors: Synthesis and in vitro evaluation

Breast cancer has been associated with an overexpression of various molecular targets; accordingly, various target-specific chemotherapeutic agents have been developed. Inhibition of ERK2, a member of MAPK pathway, is an important target involved in the treatment of both oestrogen receptor-positive and triple-negative breast cancer. Thus, in continuation of our previous work on the ERK2 target, we here report novel inhibitors of this kinase. Out of three lead molecules reported in our previous study, we selected the thiazolidinone-pyrimidine scaffold for further development of small molecule inhibitors of ERK2. Analogues of the lead molecule were docked in the target kinase, followed by molecular dynamic simulations and MM-GBSA calculations. Analogues maintaining key interactions with amino acid residues in the ATP-binding domain of ERK2 were selected and duly synthesized. In vitro biochemical evaluation of these molecules against ERK2 kinase disclosed that two molecules possess significant kinase inhibitory potential with IC₅₀ values ≤ 0.5 µM.

The QSAR Paradigm in Fragment-Based Drug Discovery: From the Virtual Generation of Target Inhibitors to Multi-Scale Modeling

Fragment-Based Drug Design (FBDD) has established itself as a promising approach in modern drug discovery, accelerating and improving lead optimization, while playing a crucial role in diminishing the high attrition rates at all stages in the drug development process. On the other hand, FBDD has benefited from the application of computational methodologies, where the models derived from the Quantitative Structure-Activity Relationships (QSAR) have become consolidated tools. This mini-review focuses on the evolution and main applications of the QSAR paradigm in the context of FBDD in the last five years. This report places particular emphasis on the QSAR models derived from fragment-based topological approaches to extract physicochemical and/or structural information, allowing to design potentially novel mono- or multi-target inhibitors from relatively large and heterogeneous databases. Here, we also discuss the need to apply multi-scale modeling, to exemplify how different datasets based on target inhibition can be simultaneously integrated and predicted together with other relevant endpoints such as the biological activity against non-biomolecular targets, as well as in vitro and in vivo toxicity and pharmacokinetic properties. In this context, seminal papers are briefly analyzed. As huge amounts of data continue to accumulate in the domains of the chemical, biological and biomedical sciences, it has become clear that drug discovery must be viewed as a multi-scale optimization process. An ideal multi-scale approach should integrate diverse chemical and biological data and also serve as a knowledge generator, enabling the design of potentially optimal chemicals that may become therapeutic agents.

Pyrazole Substituted 9-Anilinoacridines as HER2 Inhibitors Targeting Breast Cancer - An In-Silico Approach

BACKGROUND

Breast cancer is one of the malignant tumours which mainly affect the female population. Total 20% of the cases of breast cancer are due to overexpression of Human epidermal growth factor receptor-2 (HER2), which is the dominant tyrosine kinase receptor. In general, 9-anilinoacridine derivatives play an important role as antitumor agents due to their DNA-intercalating properties.

OBJECTIVE

Some novel 9-anilinoacridines substituted with pyrazole moiety(1a-z) were designed, and their HER2enzyme (PDB id-3PP0) inhibition activity was evaluated by molecular docking studies using the Glide module of Schrodinger suite 2019-4.

METHODS

Glide module of the Schrodinger suite was used to perform docking studies, qikprop module was used for in-silico ADMET screening, and the Prime-MM-GBSA module was used for free binding energy calculations. Using GLIDE scoring functions, we can determine the binding affinity of ligands (1a-z) towards HER2.

RESULTS

The inhibitory activity of ligands against HER2 was mainly due to the strong hydrophobic and hydrogen bonding interactions. Almost all the compounds 1a-z have a good binding affinity with Glide scores in the range of -4.9 to -9.75 compared to the standard drugs CK0403(-4.105) and Tamoxifen (-3.78). From the results of in-silico ADMET properties, most of the compounds fall within the recommended values. MM-GBSA binding calculations of the most potent inhibitors are more favourable.

CONCLUSION

The results of in-silico studies provide strong evidence for the consideration of valuable ligands in pyrazole substituted 9-anilinoacridines as potential HER2 inhibitors, and the compounds, 1v,s,r,d, a,o with significant Glide scores may produce significant anti-breast cancer activity for further development.

Artificial intelligence and machine learning-aided drug discovery in central nervous system diseases: State-of-the-arts and future directions

Neurological disorders significantly outnumber diseases in other therapeutic areas. However, developing drugs for central nervous system (CNS) disorders remains the most challenging area in drug discovery, accompanied with the long timelines and high attrition rates. With the rapid growth of biomedical data enabled by advanced experimental technologies, artificial intelligence (AI) and machine learning (ML) have emerged as an indispensable tool to draw meaningful insights and improve decision making in drug discovery. Thanks to the advancements in AI and ML algorithms, now the AI/ML-driven solutions have an unprecedented potential to accelerate the process of CNS drug discovery with better success rate. In this review, we comprehensively summarize AI/ML-powered pharmaceutical discovery efforts and their implementations in the CNS area. After introducing the AI/ML models as well as the conceptualization and data preparation, we outline the applications of AI/ML technologies to several key procedures in drug discovery, including target identification, compound screening, hit/lead generation and optimization, drug response and synergy prediction, de novo drug design, and drug repurposing. We review the current state-of-the-art of AI/ML-guided CNS drug discovery, focusing on blood-brain barrier permeability prediction and implementation into therapeutic discovery for neurological diseases. Finally, we discuss the major challenges and limitations of current approaches and possible future directions that may provide resolutions to these difficulties.

AKT Inhibitors: The Road Ahead to Computational Modeling-Guided Discovery

AKT, is a serine/threonine protein kinase comprising three isoforms-namely: AKT1, AKT2 and AKT3, whose inhibitors have been recognized as promising therapeutic targets for various human disorders, especially cancer. In this work, we report a systematic evaluation of multi-target Quantitative Structure-Activity Relationship (mt-QSAR) models to probe AKT' inhibitory activity, based on different feature selection algorithms and machine learning tools. The best predictive linear and non-linear mt-QSAR models were found by the genetic algorithm-based linear discriminant analysis (GA-LDA) and gradient boosting (Xgboost) techniques, respectively, using a dataset containing 5523 inhibitors of the AKT isoforms assayed under various experimental conditions. The linear model highlighted the key structural attributes responsible for higher inhibitory activity whereas the non-linear model displayed an overall accuracy higher than 90%. Both these predictive models, generated through internal and external validation methods, were then used for screening the Asinex kinase inhibitor library to identify the most potential virtual hits as pan-AKT inhibitors. The virtual hits identified were then filtered by stepwise analyses based on reverse pharmacophore-mapping based prediction. Finally, results of molecular dynamics simulations were used to estimate the theoretical binding affinity of the selected virtual hits towards the three isoforms of enzyme AKT. Our computational findings thus provide important guidelines to facilitate the discovery of novel AKT inhibitors.

QSAR Modeling for Multi-Target Drug Discovery: Designing Simultaneous Inhibitors of Proteins in Diverse Pathogenic Parasites

Parasitic diseases remain as unresolved health issues worldwide. While for some parasites the treatments involve drug combinations with serious side effects, for others, chemical therapies are inefficient due to the emergence of drug resistance. This urges the search for novel antiparasitic agents able to act through multiple mechanisms of action. Here, we report the first multi-target model based on quantitative structure-activity relationships and a multilayer perceptron neural network (mt-QSAR-MLP) to virtually design and predict versatile inhibitors of proteins involved in the survival and/or infectivity of different pathogenic parasites. The mt-QSAR-MLP model exhibited high accuracy (>80%) in both training and test sets for the classification/prediction of protein inhibitors. Several fragments were directly extracted from the physicochemical and structural interpretations of the molecular descriptors in the mt-QSAR-MLP model. Such interpretations enabled the generation of four molecules that were predicted as multi-target inhibitors against at least three of the five parasitic proteins reported here with two of the molecules being predicted to inhibit all the proteins. Docking calculations converged with the mt-QSAR-MLP model regarding the multi-target profile of the designed molecules. The designed molecules exhibited drug-like properties, complying with Lipinski’s rule of five, as well as Ghose’s filter and Veber’s guidelines.

Synthesis, Pharmacological, and Biological Evaluation of 2-Furoyl-Based MIF-1 Peptidomimetics and the Development of a General-Purpose Model for Allosteric Modulators (ALLOPTML)

This work describes the synthesis and pharmacological evaluation of 2-furoyl-based Melanostatin (MIF-1) peptidomimetics as dopamine D₂ modulating agents. Eight novel peptidomimetics were tested for their ability to enhance the maximal effect of tritiated N-propylapomorphine ([³H]-NPA) at D₂ receptors (D₂R). In this series, 2-furoyl-l-leucylglycinamide (6a) produced a statistically significant increase in the maximal [³H]-NPA response at 10 pM (11 ± 1%), comparable to the effect of MIF-1 (18 ± 9%) at the same concentration. This result supports previous evidence that the replacement of proline residue by heteroaromatic scaffolds are tolerated at the allosteric binding site of MIF-1. Biological assays performed for peptidomimetic 6a using cortex neurons from 19-day-old Wistar-Kyoto rat embryos suggest that 6a displays no neurotoxicity up to 100 μM. Overall, the pharmacological and toxicological profile and the structural simplicity of 6a makes this peptidomimetic a potential lead compound for further development and optimization, paving the way for the development of novel modulating agents of D₂R suitable for the treatment of CNS-related diseases. Additionally, the pharmacological and biological data herein reported, along with >20 000 outcomes of preclinical assays, was used to seek a general model to predict the allosteric modulatory potential of molecular candidates for a myriad of target receptors, organisms, cell lines, and biological activity parameters based on perturbation theory (PT) ideas and machine learning (ML) techniques, abbreviated as ALLOPTML. By doing so, ALLOPTML shows high specificity Sp = 89.2/89.4%, sensitivity Sn = 71.3/72.2%, and accuracy Ac = 86.1%/86.4% in training/validation series, respectively. To the best of our knowledge, ALLOPTML is the first general-purpose chemoinformatic tool using a PTML-based model for the multioutput and multicondition prediction of allosteric compounds, which is expected to save both time and resources during the early drug discovery of allosteric modulators.

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.

Cell-based multi-target QSAR model for design of virtual versatile inhibitors of liver cancer cell lines

Liver cancers are one of the leading fatal diseases among malignant neoplasms. Current chemotherapeutic treatments used to fight these illnesses have become less efficient in terms of both efficacy and safety. Therefore, there is a great need of search for new anti-liver cancer agents and this can be accelerated by using computer-aided drug discovery approaches. In this work, we report the development of the first cell-based multi-target model based on quantitative structure-activity relationships (CBMT-QSAR) for the design and prediction of chemicals as anticancer agents against 17 liver cancer cell lines. While having a good quality and predictive power (accuracy higher than 80%) in the training and test sets, respectively, the CBMT-QSAR model was employed as a tool to directly extract suitable fragments from the physicochemical and structural interpretations of the molecular descriptors. Some of these desirable fragments were assembled, leading to the virtual design of eight molecules with drug-like properties, with six of them being predicted as versatile anticancer agents against the 17 liver cancer cell lines reported here.

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.

PTML Modeling for Alzheimer’s Disease: Design and Prediction of Virtual Multi-Target Inhibitors of GSK3B, HDAC1, and HDAC6

Background:

Alzheimer’s disease is characterized by a progressive pattern of cognitive and functional impairment, which ultimately leads to death. Computational approaches have played an important role in the context of drug discovery for anti-Alzheimer's therapies. However, most of the computational models reported to date have been focused on only one protein associated with Alzheimer's, while relying on small datasets of structurally related molecules.

Objective:

We introduce the first model combining perturbation theory and machine learning based on artificial neural networks (PTML-ANN) for simultaneous prediction and design of inhibitors of three Alzheimer’s disease-related proteins, namely glycogen synthase kinase 3 beta (GSK3B), histone deacetylase 1 (HDAC1), and histone deacetylase 6 (HDAC6).

Methods:

The PTML-ANN model was obtained from a dataset retrieved from ChEMBL, and it relied on a classification approach to predict chemicals as active or inactive.

Results:

The PTML-ANN model displayed sensitivity and specificity higher than 85% in both training and test sets. The physicochemical and structural interpretation of the molecular descriptors in the model permitted the direct extraction of fragments suggested to favorably contribute to enhancing the multitarget inhibitory activity. Based on this information, we assembled ten molecules from several fragments with positive contributions. Seven of these molecules were predicted as triple target inhibitors while the remaining three were predicted as dual-target inhibitors. The estimated physicochemical properties of the designed molecules complied with Lipinski’s rule of five and its variants.

Conclusion:

This work opens new horizons toward the design of multi-target inhibitors for anti- Alzheimer's therapies.

Machine Learning Algorithm Identifies an Antibiotic Vocabulary for Permeating Gram-Negative Bacteria

Drug discovery faces a crisis. The industry has used up the "obvious" space in which to find novel drugs for biomedical applications, and productivity is declining. One strategy to combat this is rational approaches to expand the search space without relying on chemical intuition, to avoid rediscovery of similar spaces. In this work, we present proof of concept of an approach to rationally identify a "chemical vocabulary" related to a specific drug activity of interest without employing known rules. We focus on the pressing concern of multidrug resistance in Pseudomonas aeruginosa by searching for submolecules that promote compound entry into this bacterium. By synergizing theory, computation, and experiment, we validate our approach, explain the molecular mechanism behind identified fragments promoting compound entry, and select candidate compounds from an external library that display good permeation ability.

A unified in silico model based on perturbation theory for assessing the genotoxicity of metal oxide nanoparticles

Nanomaterials (NMs) are an ever-increasing field of interest, due to their wide range of applications in science and technology. However, despite providing solutions to many societal problems and challenges, NMs are associated with adverse effects with potential severe damages towards biological species and their ecosystems. Particularly, it has been confirmed that NMs may induce serious genotoxic effects on various biological targets. Given the difficulties of experimental assays for estimating the genotoxicity of many NMs on diverse biological targets, development of alternative methodologies is crucial to establish their level of safety. In silico modelling approaches, such as Quantitative Structure-Toxicity Relationships (QSTR), are now considered a promising solution for such purpose. In this work, a perturbation theory machine learning (PTML) based QSTR approach is proposed for predicting the genotoxicity of metal oxide NMs under various experimental assay conditions. The application of such perturbation approach to 6084 NM-NM pair cases, set up from 78 unique NMs, afforded a final PTML-QSTR model with an accuracy better than 96% for both training and test sets. This model was then used to predict the genotoxicity of some NMs not included in the modelling dataset. The results for this independent data set were in excellent agreement with the experimental ones. Overall, that thus suggests that the derived PTML-QSTR model is a reliable in silico tool to rapidly and cost-efficiently assess the genotoxicity of metal oxide NMs. Finally, and most importantly, the model provides important insights regarding the mechanism of the genotoxicity triggered by these NMs.

A Review on Applications of Computational Methods in Drug Screening and Design

Drug development is one of the most significant processes in the pharmaceutical industry. Various computational methods have dramatically reduced the time and cost of drug discovery. In this review, we firstly discussed roles of multiscale biomolecular simulations in identifying drug binding sites on the target macromolecule and elucidating drug action mechanisms. Then, virtual screening methods (e.g., molecular docking, pharmacophore modeling, and QSAR) as well as structure- and ligand-based classical/de novo drug design were introduced and discussed. Last, we explored the development of machine learning methods and their applications in aforementioned computational methods to speed up the drug discovery process. Also, several application examples of combining various methods was discussed. A combination of different methods to jointly solve the tough problem at different scales and dimensions will be an inevitable trend in drug screening and design.

Interaction between piperine and genes associated with sciatica and its mechanism based on molecular docking technology and network pharmacology

Abstract

Piperine is the main active component of Piper longum L., which is also the main component of anti-sciatica Mongolian medicine Naru Sanwei pill. It has many pharmacological activities such as anti-inflammatory and immune regulation. This paper aims to preliminarily explore the potential mechanism of piperine in the treatment of sciatica through network pharmacology and molecular docking. TCMSP, ETCM database and literature mining were used to collect the active compounds of Piper longum L. Swiss TargetPrediction and SuperPred server were used to find the targets of compounds. At the same time, CTD database was used to collect the targets of sciatica. Then the above targets were compared and analyzed to select the targets of anti-sciatica in Piper longum L. The Go (gene ontology) annotation and KEGG pathway of the targets were enriched and analyzed by Metascape database platform. The molecular docking between the effective components and the targets was verified by Autodock. After that, the sciatica model of rats was established and treated with piperine. The expression level of inflammatory factors and proteins in the serum and tissues of rat sciatic nerve were detected by ELISA and Western blot. HE staining and immunohistochemistry were carried out on the sciatica tissues of rats. The results showed that Piper longum L. can regulate the development of sciatica and affect the expressions of PPARG and NF-kB1 through its active ingredient piperine, and there is endogenous interaction between PPARG and NF-kB1.

Graphic abstract

Molecular mechanism of inhibitor bindings to bromodomain-containing protein 9 explored based on molecular dynamics simulations and calculations of binding free energies

Recently, bromodomain-containing protein 9 (BRD9) has been a prospective therapeutic target for anticancer drug design. Molecular dynamics (MD) simulations combined with molecular mechanics generalized Born surface area (MM-GBSA) method were adopted to explore binding modes of three inhibitors (5SW, 5U2, and 5U6) to BRD9 and identify the hot spot of the inhibitor-BRD9 binding. The results indicate that the inhibitor 5SW has the strongest binding ability to BRD9 among the current three inhibitors. Furthermore, the rank of the binding free energies predicted by MM-GBSA approach agrees with that determined by the experimental values. In addition, inhibitor-residue interactions were computed by using residue-based free-energy decomposition method and the results suggest that residue His42 produces the CH-H interactions, residues Asn100, Ile53 and Val49 produce the CH-[Formula: see text] interactions with three inhibitors and Tyr106, Phe45 and Phe44 generate the π-π interactions with inhibitors. Notably, the residue Asn140 forms hydrogen bonding interactions with three inhibitors. This research is expected to provide useful molecular basis and dynamics information at atomic levels for the design of potent inhibitors inhibiting the activity of BRD9.

MCDCalc: Markov Chain Molecular Descriptors Calculator for Medicinal Chemistry

AIMS

Cheminformatics models are able to predict different outputs (activity, property, chemical reactivity) in single molecules or complex molecular systems (catalyzed organic synthesis, metabolic reactions, nanoparticles, etc.).

BACKGROUND

Cheminformatics models are able to predict different outputs (activity, property, chemical reactivity) in single molecules or complex molecular systems (catalyzed organic synthesis, metabolic reactions, nanoparticles, etc.).

OBJECTIVE

Cheminformatics prediction of complex catalytic enantioselective reactions is a major goal in organic synthesis research and chemical industry. Markov Chain Molecular Descriptors (MCDs) have been largely used to solve Cheminformatics problems. There are different types of Markov chain descriptors such as Markov-Shannon entropies (Shk), Markov Means (Mk), Markov Moments (πk), etc. However, there are other possible MCDs that have not been used before. In addition, the calculation of MCDs is done very often using specific software not always available for general users and there is not an R library public available for the calculation of MCDs. This fact, limits the availability of MCMDbased Cheminformatics procedures.

METHODS

We studied the enantiomeric excess ee(%)[Rcat] for 324 α-amidoalkylation reactions. These reactions have a complex mechanism depending on various factors. The model includes MCDs of the substrate, solvent, chiral catalyst, product along with values of time of reaction, temperature, load of catalyst, etc. We tested several Machine Learning regression algorithms. The Random Forest regression model has R2 > 0.90 in training and test. Secondly, the biological activity of 5644 compounds against colorectal cancer was studied.

RESULTS

We developed very interesting model able to predict with Specificity and Sensitivity 70-82% the cases of preclinical assays in both training and validation series.

CONCLUSION

The work shows the potential of the new tool for computational studies in organic and medicinal chemistry.

Construction of Quantitative Structure Activity Relationship (QSAR) Models to Predict Potency of Structurally Diversed Janus Kinase 2 Inhibitors

Janus kinase 2 (JAK2) inhibitors represent a promising therapeutic class of anticancer agents against many myeloproliferative disorders. Bioactivity data on pIC50 of 2229 JAK2 inhibitors were employed in the construction of quantitative structure-activity relationship (QSAR) models. The models were built from 100 data splits using decision tree (DT), support vector machine (SVM), deep neural network (DNN) and random forest (RF). The predictive power of RF models were assessed via 10-fold cross validation, which afforded excellent predictive performance with R2 and RMSE of 0.74 ± 0.05 and 0.63 ± 0.05, respectively. Moreover, test set has excellent performance of R2 (0.75 ± 0.03) and RMSE (0.62 ± 0.04). In addition, Y-scrambling was utilized to evaluate the possibility of chance correlation of the predictive model. A thorough analysis of the substructure fingerprint count was conducted to provide insights on the inhibitory properties of JAK2 inhibitors. Molecular cluster analysis revealed that pyrazine scaffolds have nanomolar potency against JAK2.

Phytochemicals in Anticancer Drug Development

BACKGROUND

In spite of major technological advances in conventional therapies, cancer continues to remain the leading cause of mortality worldwide. Phytochemicals are gradually emerging as a rich source of effective but safer agents against many life-threatening diseases.

METHODS

Various phytochemicals with reported anticancer activity have been simply categorized into major phytoconstituents- alkaloids, polyphenols, saponins, tannins and terpenoids.

RESULTS

The adverse effects associated with currently available anticancer medications may be overcome by using plant-derived compounds either alone or in combination. Exploration of plant kingdom may provide new leads for the accelerated development of new anticancer agents.

CONCLUSION

Although numerous potent synthetic drugs have been introduced for cancer chemotherapy, yet their serious toxicity concerns to normal cells apart from drug resistance have emerged as the major obstacles for their clinical utility over a prolonged duration of time. Current status and potential of phytochemicals and their derivatives in cancer therapy have been briefly reviewed in the present manuscript.

Multi-Target Chemometric Modelling, Fragment Analysis and Virtual Screening with ERK Inhibitors as Potential Anticancer Agents

Two isoforms of extracellular regulated kinase (ERK), namely ERK-1 and ERK-2, are associated with several cellular processes, the aberration of which leads to cancer. The ERK-1/2 inhibitors are thus considered as potential agents for cancer therapy. Multitarget quantitative structure–activity relationship (mt-QSAR) models based on the Box–Jenkins approach were developed with a dataset containing 6400 ERK inhibitors assayed under different experimental conditions. The first mt-QSAR linear model was built with linear discriminant analysis (LDA) and provided information regarding the structural requirements for better activity. This linear model was also utilised for a fragment analysis to estimate the contributions of ring fragments towards ERK inhibition. Then, the random forest (RF) technique was employed to produce highly predictive non-linear mt-QSAR models, which were used for screening the Asinex kinase library and identify the most potential virtual hits. The fragment analysis results justified the selection of the hits retrieved through such virtual screening. The latter were subsequently subjected to molecular docking and molecular dynamics simulations to understand their possible interactions with ERK enzymes. The present work, which utilises in-silico techniques such as multitarget chemometric modelling, fragment analysis, virtual screening, molecular docking and dynamics, may provide important guidelines to facilitate the discovery of novel ERK inhibitors.

Development of Multi-Target Chemometric Models for the Inhibition of Class I PI3K Enzyme Isoforms: A Case Study Using QSAR-Co Tool

The present work aims at establishing multi-target chemometric models using the recently launched quantitative structure–activity relationship (QSAR)-Co tool for predicting the activity of inhibitor compounds against different isoforms of phosphoinositide 3-kinase (PI3K) under various experimental conditions. The inhibitors of class I phosphoinositide 3-kinase (PI3K) isoforms have emerged as potential therapeutic agents for the treatment of various disorders, especially cancer. The cell-based enzyme inhibition assay results of PI3K inhibitors were curated from the CHEMBL database. Factors such as the nature and mutation of cell lines that may significantly alter the assay outcomes were considered as important experimental elements for mt-QSAR model development. The models, in turn, were developed using two machine learning techniques as implemented in QSAR-Co: linear discriminant analysis (LDA) and random forest (RF). Both techniques led to models with high accuracy (ca. 90%). Several molecular fragments were extracted from the current dataset, and their quantitative contributions to the inhibitory activity against all the proteins and experimental conditions under study were calculated. This case study also demonstrates the utility of QSAR-Co tool in solving multi-factorial and complex chemometric problems. Additionally, the combination of different in silico methods employed in this work can serve as a valuable guideline to speed up early discovery of PI3K inhibitors.

Concepts of Artificial Intelligence for Computer-Assisted Drug Discovery

Artificial intelligence (AI), and, in particular, deep learning as a subcategory of AI, provides opportunities for the discovery and development of innovative drugs. Various machine learning approaches have recently (re)emerged, some of which may be considered instances of domain-specific AI which have been successfully employed for drug discovery and design. This review provides a comprehensive portrayal of these machine learning techniques and of their applications in medicinal chemistry. After introducing the basic principles, alongside some application notes, of the various machine learning algorithms, the current state-of-the art of AI-assisted pharmaceutical discovery is discussed, including applications in structure- and ligand-based virtual screening, de novo drug design, physicochemical and pharmacokinetic property prediction, drug repurposing, and related aspects. Finally, several challenges and limitations of the current methods are summarized, with a view to potential future directions for AI-assisted drug discovery and design.

QSAR-Co: An Open Source Software for Developing Robust Multitasking or Multitarget Classification-Based QSAR Models

Quantitative structure-activity relationships (QSAR) modeling is a well-known computational technique with wide applications in fields such as drug design, toxicity predictions, nanomaterials, etc. However, QSAR researchers still face certain problems to develop robust classification-based QSAR models, especially while handling response data pertaining to diverse experimental and/or theoretical conditions. In the present work, we have developed an open source standalone software "QSAR-Co" (available to download at https://sites.google.com/view/qsar-co ) to setup classification-based QSAR models that allow mining the response data coming from multiple conditions. The software comprises two modules: (1) the Model development module and (2) the Screen/Predict module. This user-friendly software provides several functionalities required for developing a robust multitasking or multitarget classification-based QSAR model using linear discriminant analysis or random forest techniques, with appropriate validation, following the principles set by the Organisation for Economic Co-operation and Development (OECD) for applying QSAR models in regulatory assessments.

In-silico Design, ADMET Screening, MM-GBSA Binding Free Energy of Some Novel Isoxazole Substituted 9-Anilinoacridines as HER2 Inhibitors Targeting Breast Cancer

BACKGROUND

Human Epidermal development factor Receptor-2 (HER2) is a membrane tyrosine kinase which is overexpressed and gene amplified in human breast cancers. HER2 amplification and overexpression have been linked to important tumor cell proliferation and survival pathways for 20% of instances of breast cancer. 9-aminoacridines are significant DNA-intercalating agents because of their antiproliferative properties.

OBJECTIVE

Some novel isoxazole substituted 9-anilinoacridines(1a-z) were designed by in-silico technique for their HER2 inhibitory activity. Docking investigations of compounds 1a-z are performed against HER2 (PDB id-3PP0) by using Schrodinger suit 2016-2.

METHODS

Molecular docking study for the designed molecules 1a-z are performed by Glide module, in-silico ADMET screening by QikProp module and binding free energy by Prime-MMGBSA module of Schrodinger suit. The binding affinity of designed molecules 1a-z towards HER2 was chosen based on GLIDE score.

RESULTS

Many compounds showed good hydrophobic communications and hydrogen bonding associations to hinder HER2. The compounds 1a-z, aside from 1z have significant Glide scores in the scope of - 4.91 to - 10.59 when compared with the standard Ethacridine (- 4.23) and Tamoxifen (- 3.78). The in-silico ADMET properties are inside the suggested about drug likeness. MM-GBSA binding of the most intense inhibitor is positive.

CONCLUSION

The outcomes reveal that this study provides evidence for the consideration of isoxazole substituted 9-aminoacridine derivatives as potential HER2 inhibitors. The compounds, 1s,x,v,a,j,r with significant Glide scores may produce significant anti breast cancer activity and further in vitro and in vivo investigations may prove their therapeutic potential.

Combining Ensemble Learning with a Fragment-Based Topological Approach To Generate New Molecular Diversity in Drug Discovery: In Silico Design of Hsp90 Inhibitors

Machine learning methods have revolutionized modern science, providing fast and accurate solutions to multiple problems. However, they are commonly treated as "black boxes". Therefore, in important scientific fields such as medicinal chemistry and drug discovery, machine learning methods are restricted almost exclusively to the task of performing predictions of large and heterogeneous data sets of chemicals. The lack of interpretability prevents the full exploitation of the machine learning models as generators of new chemical knowledge. This work focuses on the development of an ensemble learning model for the prediction and design of potent dual heat shock protein 90 (Hsp90) inhibitors. The model displays accuracy higher than 80% in both training and test sets. To use the ensemble model as a generator of new chemical knowledge, three steps were followed. First, a physicochemical and/or structural interpretation was provided for each molecular descriptor present in the ensemble learning model. Second, the term "pseudolinear equation" was introduced within the context of machine learning to calculate the relative quantitative contributions of different molecular fragments to the inhibitory activity against the two Hsp90 isoforms studied here. Finally, by assembling the fragments with positive contributions, new molecules were designed, being predicted as potent Hsp90 inhibitors. According to Lipinski's rule of five, the designed molecules were found to exhibit potentially good oral bioavailability, a primordial property that chemicals must have to pass early stages in drug discovery. The present approach based on the combination of ensemble learning and fragment-based topological design holds great promise in drug discovery, and it can be adapted and applied to many different scientific disciplines.

PTML Combinatorial Model of ChEMBL Compounds Assays for Multiple Types of Cancer

Determining the target proteins of new anticancer compounds is a very important task in Medicinal Chemistry. In this sense, chemists carry out preclinical assays with a high number of combinations of experimental conditions (c _j). In fact, ChEMBL database contains outcomes of 65 534 different anticancer activity preclinical assays for 35 565 different chemical compounds (1.84 assays per compound). These assays cover different combinations of c _j formed from >70 different biological activity parameters ( c₀), >300 different drug targets ( c₁), >230 cell lines ( c₂), and 5 organisms of assay ( c₃) or organisms of the target ( c₄). It include a total of 45 833 assays in leukemia, 6227 assays in breast cancer, 2499 assays in ovarian cancer, 3499 in colon cancer, 3159 in lung cancer, 2750 in prostate cancer, 601 in melanoma, etc. This is a very complex data set with multiple Big Data features. This data is hard to be rationalized by researchers to extract useful relationships and predict new compounds. In this context, we propose to combine perturbation theory (PT) ideas and machine learning (ML) modeling to solve this combinatorial-like problem. In this work, we report a PTML (PT + ML) model for ChEMBL data set of preclinical assays of anticancer compounds. This is a simple linear model with only three variables. The model presented values of area under receiver operating curve = AUROC = 0.872, specificity = Sp(%) = 90.2, sensitivity = Sn(%) = 70.6, and overall accuracy = Ac(%) = 87.7 in training series. The model also have Sp(%) = 90.1, Sn(%) = 71.4, and Ac(%) = 87.8 in external validation series. The model use PT operators based on multicondition moving averages to capture all the complexity of the data set. We also compared the model with nonlinear artificial neural network (ANN) models obtaining similar results. This confirms the hypothesis of a linear relationship between the PT operators and the classification as anticancer compounds in different combinations of assay conditions. Last, we compared the model with other PTML models reported in the literature concluding that this is the only one PTML model able to predict activity against multiple types of cancer. This model is a simple but versatile tool for the prediction of the targets of anticancer compounds taking into consideration multiple combinations of experimental conditions in preclinical assays.

BET bromodomain inhibitors: fragment-based in silico design using multi-target QSAR models

Epigenetics has become a focus of interest in drug discovery. In this sense, bromodomain-containing proteins have emerged as potential epigenetic targets in cancer research and other therapeutic areas. Several computational approaches have been applied to the prediction of bromodomain inhibitors. Nevertheless, such approaches have several drawbacks such as the fact that they predict activity against only one bromodomain-containing protein, using structurally related compounds. Also, there are no reports focused on meaningfully analyzing the physicochemical/structural features that are necessary for the design of a bromodomain inhibitor. This work describes the development of two different multi-target models based on quantitative structure-activity relationships (mt-QSAR) for the prediction and in silico design of multi-target bromodomain inhibitors against the proteins BRD2, BRD3, and BRD4. The first model relied on linear discriminant analysis (LDA) while the second focused on artificial neural networks. Both models exhibited accuracies higher than 85% in the dataset. Several molecular fragments were extracted, and their contributions to the inhibitory activity against the three BET proteins were calculated by the LDA model. Six molecules were designed by assembling the fragments with positive contributions, and they were predicted as multi-target BET bromodomain inhibitors by the two mt-QSAR models. Molecular docking calculations converged with the predictions performed by the mt-QSAR models, suggesting that the designed molecules can exhibit potent activity against the three BET proteins. These molecules complied with the Lipinski's rule of five.

Computational Fragment-Based Drug Design: Current Trends, Strategies, and Applications

Fragment-based drug design (FBDD) has become an effective methodology for drug development for decades. Successful applications of this strategy brought both opportunities and challenges to the field of Pharmaceutical Science. Recent progress in the computational fragment-based drug design provide an additional approach for future research in a time- and labor-efficient manner. Combining multiple in silico methodologies, computational FBDD possesses flexibilities on fragment library selection, protein model generation, and fragments/compounds docking mode prediction. These characteristics provide computational FBDD superiority in designing novel and potential compounds for a certain target. The purpose of this review is to discuss the latest advances, ranging from commonly used strategies to novel concepts and technologies in computational fragment-based drug design. Particularly, in this review, specifications and advantages are compared between experimental and computational FBDD, and additionally, limitations and future prospective are discussed and emphasized.

Discovery of Novel KRAS-PDEδ Inhibitors by Fragment-Based Drug Design

Targeting KRAS-PDEδ protein-protein interactions with small molecules represents a promising opportunity for developing novel antitumor agents. However, current KRAS-PDEδ inhibitors are limited by poor cellular antitumor potency and the druggability of the target remains to be validated by new inhibitors. To tackle these challenges, herein, novel, highly potent KRAS-PDEδ inhibitors were identified by fragment-based drug design, providing promising lead compounds or chemical probes for investigating the biological functions and druggability of KRAS-PDEδ interaction.

Discovery of VEGFR2 inhibitors by integrating naïve Bayesian classification, molecular docking and drug screening approaches

The high morbidity and mortality of cancer make it one of the leading causes of global death, thus it is an urgent need to develop effective drugs for cancer therapy.

From Classification to Regression Multitasking QSAR Modeling Using a Novel Modular Neural Network: Simultaneous Prediction of Anticonvulsant Activity and Neurotoxicity of Succinimides

Succinimides, which contain a pharmacophore responsible for anticonvulsant activity, are frequently used antiepileptic drugs and the synthesis of their new derivatives with improved efficacy and tolerability presents an important task. Nowadays, multitarget/tasking methodologies focused on quantitative-structure activity relationships (mt-QSAR/mtk-QSAR) have an important role in the rational design of drugs since they enable simultaneous prediction of several standard measures of biological activities at diverse experimental conditions and against different biological targets. Relating to this very topic, the mt-QSAR/mtk-QSAR methodology can give only binary classification models, and as such, in this study a regression mtk-QSAR (rmtk-QSAR) model based on a novel modular neural network (MNN) has been proposed. The MNN uses standard classification mtk-QSAR models as input modules, while the regression is performed by the output module. The rmtk-QSAR model has been successfully developed for the simultaneous prediction of anticonvulsant activity and neurotoxicity of succinimides, with a satisfactory accuracy in testing (R² = 0.87). Thus, the proposed mtk-QSAR regression method can be regarded as a viable alternative to the standard QSAR methodology.

Exploring differential evolution for inverse QSAR analysis

Inverse quantitative structure-activity relationship (QSAR) modeling encompasses the generation of compound structures from values of descriptors corresponding to high activity predicted with a given QSAR model. Structure generation proceeds from descriptor coordinates optimized for activity prediction. Herein, we concentrate on the first phase of the inverse QSAR process and introduce a new methodology for coordinate optimization, termed differential evolution (DE), that originated from computer science and engineering. Using simulation and compound activity data, we demonstrate that DE in combination with support vector regression (SVR) yields effective and robust predictions of optimized coordinates satisfying model constraints and requirements. For different compound activity classes, optimized coordinates are obtained that exclusively map to regions of high activity in feature space, represent novel positions for structure generation, and are chemically meaningful.