Discovery of novel anti-inflammatory drug-like compounds by aligning in silico and in vivo screening: The nitroindazolinone chemotype

PTML Modeling for Pancreatic Cancer Research: In Silico Design of Simultaneous Multi-Protein and Multi-Cell Inhibitors

Pancreatic cancer (PANC) is a dangerous type of cancer that is a major cause of mortality worldwide and exhibits a remarkably poor prognosis. To date, discovering anti-PANC agents remains a very complex and expensive process. Computational approaches can accelerate the search for anti-PANC agents. We report for the first time two models that combined perturbation theory with machine learning via a multilayer perceptron network (PTML-MLP) to perform the virtual design and prediction of molecules that can simultaneously inhibit multiple PANC cell lines and PANC-related proteins, such as caspase-1, tumor necrosis factor-alpha (TNF-alpha), and the insulin-like growth factor 1 receptor (IGF1R). Both PTML-MLP models exhibited accuracies higher than 78%. Using the interpretation from one of the PTML-MLP models as a guideline, we extracted different molecular fragments desirable for the inhibition of the PANC cell lines and the aforementioned PANC-related proteins and then assembled some of those fragments to form three new molecules. The two PTML-MLP models predicted the designed molecules as potentially versatile anti-PANC agents through inhibition of the three PANC-related proteins and multiple PANC cell lines. Conclusions: This work opens new horizons for the application of the PTML modeling methodology to anticancer research.

A Comprehensive In Silico Exploration of Pharmacological Properties, Bioactivities, Molecular Docking, and Anticancer Potential of Vieloplain F from Xylopia vielana Targeting B-Raf Kinase

Compounds derived from plants have several anticancer properties. In the current study, one guaiane-type sesquiterpene dimer, vieloplain F, isolated from Xylopia vielana species, was tested against B-Raf kinase protein (PDB: 3OG7), a potent target for melanoma. A comprehensive in silico analysis was conducted in this research to understand the pharmacological properties of a compound encompassing absorption, distribution, metabolism, excretion, and toxicity (ADMET), bioactivity score predictions, and molecular docking. During ADMET estimations, the FDA-approved medicine vemurafenib was hepatotoxic, cytochrome-inhibiting, and non-cardiotoxic compared to the vieloplain F. The bioactivity scores of vieloplain F were active for nuclear receptor ligand and enzyme inhibitor. During molecular docking experiments, the compound vieloplain F has displayed a higher binding potential with −11.8 kcal/mol energy than control vemurafenib −10.2 kcal/mol. It was shown that intermolecular interaction with the B-Raf complex and the enzyme’s active gorge through hydrogen bonding and hydrophobic contacts was very accurate for the compound vieloplain F, which was then examined for MD simulations. In addition, simulations using MM-GBSA showed that vieloplain F had the greatest propensity to bind to active site residues. The vieloplain F has predominantly represented a more robust profile compared to control vemurafenib, and these results opened the road for vieloplain F for its utilization as a plausible anti-melanoma agent and anticancer drug in the next era.

Computational Exploration of Anti-Cancer Potential of Guaiane Dimers from Xylopia vielana by Targeting B-Raf Kinase Using Chemo-Informatics, Molecular Docking and MD Simulation Studies

BACKGROUND

Natural products from herbs are prolific to display robust anticancer activities.

OBJECTIVES

In the current study, B-Raf kinase protein (PDB: 3OG7), a potent target for melanoma, was tested against two guaiane-type sesquiterpene dimers, xylopin E-F, obtained from Xylopia vielana.

METHODS

In this work, a systematic in silico study using ADMET analysis, bioactivity score forecasts, molecular docking, and its simulations were conducted to understand compounds' pharmacological properties.

RESULTS

During ADMET predictions of both the compounds, Xylopin E-F has displayed a safer profile in hepatotoxicity, cytochrome inhibition, and only xylopin F displayed as non-cardiotoxic compared to FDA approved drug vemurafenib. Both the compounds were proceeded to molecular docking experiments using Autodock docking software and both the compounds Xylopin E-F have displayed higher binding potential with -11.5Kcal/mol energy compared to control vemurafenib -10.2 Kcal/mol. All the compounds were further evaluated for their MD simulations and their molecular interactions with the B-Raf kinase complex displayed precise interactions with the active gorge of the enzyme by hydrogen bonding.

CONCLUSIONS

Overall, xylopin F had a better profile relative to xylopin E and vemurafenib, and these findings indicated that this bio-molecule could be used as an anti-melanoma agent and as a possible anticancer drug in the future. Therefore, this is a systematic optimized in silico approach to creating an anticancer pathway for guaiane dimers against the backdrop of its potential for future drug development.

Computational Drug Repurposing for Antituberculosis Therapy: Discovery of Multi-Strain Inhibitors

Tuberculosis remains the most afflicting infectious disease known by humankind, with one quarter of the population estimated to have it in the latent state. Discovering antituberculosis drugs is a challenging, complex, expensive, and time-consuming task. To overcome the substantial costs and accelerate drug discovery and development, drug repurposing has emerged as an attractive alternative to find new applications for “old” drugs and where computational approaches play an essential role by filtering the chemical space. This work reports the first multi-condition model based on quantitative structure–activity relationships and an ensemble of neural networks (mtc-QSAR-EL) for the virtual screening of potential antituberculosis agents able to act as multi-strain inhibitors. The mtc-QSAR-EL model exhibited an accuracy higher than 85%. A physicochemical and fragment-based structural interpretation of this model was provided, and a large dataset of agency-regulated chemicals was virtually screened, with the mtc-QSAR-EL model identifying already proven antituberculosis drugs while proposing chemicals with great potential to be experimentally repurposed as antituberculosis (multi-strain inhibitors) agents. Some of the most promising molecules identified by the mtc-QSAR-EL model as antituberculosis agents were also confirmed by another computational approach, supporting the capabilities of the mtc-QSAR-EL model as an efficient tool for computational drug repurposing.

A novel indazole derivative, compound Cyy-272, attenuates LPS-induced acute lung injury by inhibiting JNK phosphorylation

Acute lung injury (ALI) is a diffuse lung dysfunction disease characterized by high prevalence and high mortality. Thus far, no effective pharmacological treatment has been made for ALI in clinics. Inflammation is critical to the development of ALI. Therefore, anti-inflammation may be a potential therapy strategy for ALI. Indazole-containing derivatives, representing one of the most important heterocycles in drug molecules, are endowed with a broad range of biological properties, such as anti-cancer and anti-inflammation. In the current study, we investigated the biological effects of Cyy-272, a newly synthesized indazole compound, on LPS-induced ALI both in vivo and in vitro. Results show that Cyy-272 can inhibit the release of inflammatory cytokines in LPS-stimulated macrophage and alleviate LPS induced ALI. Further experiment revealed that Cyy-272 exhibit anti-inflammation activity by inhibiting JNK phosphorylation. Overall, our studies show that an indazole derivative, Cyy-272, is effective in suppressing LPS-induced JNK activation and inflammatory signaling.

Zebrafish Models for the Safety and Therapeutic Testing of Nanoparticles with a Focus on Macrophages

New nanoparticles and biomaterials are increasingly being used in biomedical research for drug delivery, diagnostic applications, or vaccines, and they are also present in numerous commercial products, in the environment and workplaces. Thus, the evaluation of the safety and possible therapeutic application of these nanomaterials has become of foremost importance for the proper progress of nanotechnology. Due to economical and ethical issues, in vitro and in vivo methods are encouraged for the testing of new compounds and/or nanoparticles, however in vivo models are still needed. In this scenario, zebrafish (Danio rerio) has demonstrated potential for toxicological and pharmacological screenings. Zebrafish presents an innate immune system, from early developmental stages, with conserved macrophage phenotypes and functions with respect to humans. This fact, combined with the transparency of zebrafish, the availability of models with fluorescently labelled macrophages, as well as a broad variety of disease models offers great possibilities for the testing of new nanoparticles. Thus, with a particular focus on macrophage-nanoparticle interaction in vivo, here, we review the studies using zebrafish for toxicological and biodistribution testing of nanoparticles, and also the possibilities for their preclinical evaluation in various diseases, including cancer and autoimmune, neuroinflammatory, and infectious diseases.

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.

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.

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.

PTML Model of ChEMBL Compounds Assays for Vitamin Derivatives

Determining the biological activity of vitamin derivatives is needed given that organic synthesis of analogs of vitamins is an active field of interest for medicinal chemistry, pharmaceuticals, and food additives. Accordingly, scientists from different disciplines perform preclinical assays (n_ij) with a considerable combination of assay conditions (c_j). Indeed, the ChEMBL platform contains a database that includes results from 36 220 different biological activity bioassays of 21 240 different vitamins and vitamin derivatives. These assays present are heterogeneous in terms of assay combinations of c_j. They are focused on >500 different biological activity parameters (c₀), >340 different targets (c₁), >6200 types of cell (c₂), >120 organisms of assay (c₃), and >60 assay strains (c₄). It includes a total of >1850 niacin assays, >1580 tretinoin assays, >1580 retinol assays, 857 ascorbic acid assays, etc. Given the complexity of this combinatorial data in terms of being assimilated by researchers, we propose to build a model by combining perturbation theory (PT) and machine learning (ML). Through this study, we propose a PTML (PT + ML) combinatorial model for ChEMBL results on biological activity of vitamins and vitamins derivatives. The linear discriminant analysis (LDA) model presented the following results for training subset a: specificity (%) = 90.38, sensitivity (%) = 87.51, and accuracy (%) = 89.89. The model showed the following results for the external validation subset: specificity (%) = 90.58, sensitivity (%) = 87.72, and accuracy (%) = 90.09. Different types of linear and nonlinear PTML models, such as logistic regression (LR), classification tree (CT), näive Bayes (NB), and random Forest (RF), were applied to contrast the capacity of prediction. The PTML-LDA model predicts with more accuracy by applying combinatorial descriptors. In addition, a PCA experiment with chemical structure descriptors allowed us to characterize the high structural diversity of the chemical space studied. In any case, PTML models using chemical structure descriptors do not improve the performance of the PTML-LDA model based on ALOGP and PSA. We can conclude that the three variable PTML-LDA model is a simplified and adaptable tool for the prediction, for different experiment combinations, the biological activity of derivative vitamins.

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.

Ensemble-Based Modeling of Chemical Compounds with Antimalarial Activity

Background:

Malaria or Paludism is a tropical disease caused by parasites of the Plasmodium genre and transmitted to humans through the bite of infected mosquitos of the Anopheles genre. This pathology is considered one of the first causes of death in tropical countries and, despite several existing therapies, they have a high toxicity. Computational methods based on Quantitative Structure- Activity Relationship studies have been widely used in drug design work flows.

Objective:

The main goal of the current research is to develop computational models for the identification of antimalarial hit compounds.

Materials and Methods:

For this, a data set suitable for the modeling of the antimalarial activity of chemical compounds was compiled from the literature and subjected to a thorough curation process. In addition, the performance of a diverse set of ensemble-based classification methodologies was evaluated and one of these ensembles was selected as the most suitable for the identification of antimalarial hits based on its virtual screening performance. Data curation was conducted to minimize noise. Among the explored ensemble-based methods, the one combining Genetic Algorithms for the selection of the base classifiers and Majority Vote for their aggregation showed the best performance.

Results:

Our results also show that ensemble modeling is an effective strategy for the QSAR modeling of highly heterogeneous datasets in the discovery of potential antimalarial compounds.

Conclusion:

It was determined that the best performing ensembles were those that use Genetic Algorithms as a method of selection of base models and Majority Vote as the aggregation method.

Modeling Antibacterial Activity with Machine Learning and Fusion of Chemical Structure Information with Microorganism Metabolic Networks

Predicting the activity of new chemical compounds over pathogenic microorganisms with different metabolic reaction networks (MRN _s) is an important goal due to the different susceptibility to antibiotics. The ChEMBL database contains >160 000 outcomes of preclinical assays of antimicrobial activity for 55 931 compounds with >365 parameters of activity (MIC, IC₅₀, etc.) and >90 bacteria strains of >25 bacterial species. In addition, the Leong and Barabàsi data set includes >40 MRNs of microorganisms. However, there are no models able to predict antibacterial activity for multiple assays considering both drug and MRN structures at the same time. In this work, we combined perturbation theory, machine learning, and information fusion techniques to develop the first PTMLIF model. The best linear model found presented values of specificity = 90.31/90.40 and sensitivity = 88.14/88.07 in training/validation series. We carried out a comparison to nonlinear artificial neural network (ANN) techniques and previous models from the literature. Next, we illustrated the practical use of the model with an experimental case of study. We reported for the first time the isolation and characterization of terpenes from the plant Cissus incisa. The antibacterial activity of the terpenes was experimentally determined. The more active compounds were phytol and α-amyrin, with MIC = 100 μg/mL for Vancomycin-resistant Enterococcus faecium and Acinetobacter baumannii resistant to carbapenems. These compounds are already known from other sources. However, they have been isolated and evaluated for the first time here against several strains of multidrug-resistant bacteria including World Health Organization (WHO) priority pathogens. Last, we used the model to predict the activity of these compounds versus other microorganisms with different MRNs in order to find other potential targets.

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.

Discrete Fourier Transform-Based Multivariate Image Analysis: Application to Modeling of Aromatase Inhibitory Activity

We recently generalized the formerly alignment-dependent multivariate image analysis applied to quantitative structure-activity relationships (MIA-QSAR) method through the application of the discrete Fourier transform (DFT), allowing for its application to noncongruent and structurally diverse chemical compound data sets. Here we report the first practical application of this method in the screening of molecular entities of therapeutic interest, with human aromatase inhibitory activity as the case study. We developed an ensemble classification model based on the two-dimensional (2D) DFT MIA-QSAR descriptors, with which we screened the NCI Diversity Set V (1593 compounds) and obtained 34 chemical compounds with possible aromatase inhibitory activity. These compounds were docked into the aromatase active site, and the 10 most promising compounds were selected for in vitro experimental validation. Of these compounds, 7419 (nonsteroidal) and 89 201 (steroidal) demonstrated satisfactory antiproliferative and aromatase inhibitory activities. The obtained results suggest that the 2D-DFT MIA-QSAR method may be useful in ligand-based virtual screening of new molecular entities of therapeutic utility.

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

Breast cancer is the most frequent cancer reported in women, being responsible for hundreds of thousands of deaths. Chemotherapy has proven to be effective against this malignant neoplasm depending on different biological factors such as the histopathology, grade, and stage, among others. However, breast cancer cells have become resistant to current chemotherapeutic regimens, urging the discovery of new anti-breast cancer drugs. Computational approaches have the potential to offer promising alternatives to accelerate the search for potent and versatile anti-breast cancer agents. In the present work, we introduce the first multitasking (mtk) computational model devoted to the in silico fragment-based design of new molecules with high inhibitory activity against 19 different proteins involved in breast cancer. The mtk-computational model was created from a dataset formed by 24,285 cases, and it exhibited accuracy around 93% in both training and prediction (test) sets. Several molecular fragments were extracted from the molecules present in the dataset, and their quantitative contributions to the inhibitory activities against all the proteins under study were calculated. The combined use of the fragment contributions and the physicochemical interpretations of the different molecular descriptors in the mtk-computational model allowed the design of eight new molecular entities not reported in our dataset. These molecules were predicted as potent multi-target inhibitors against all the proteins, and they exhibited a desirable druglikeness according to the Lipinski's rule of five and its variants.

Molecular topology applied to the discovery of 1-benzyl-2-(3-fluorophenyl)-4-hydroxy-3-(3-phenylpropanoyl)-2H-pyrrole-5-one as a non-ligand-binding-pocket antiandrogen

We report the discovery of 1-benzyl-2-(3-fluorophenyl)-4-hydroxy-3-(3-phenylpropanoyl)-2H-pyrrole-5-one as a novel non-ligand binding pocket (non-LBP) antagonist of the androgen receptor (AR) through the application of molecular topology techniques. This compound, validated through time-resolved fluorescence resonance energy transfer and fluorescence polarization biological assays, provides the basis for lead optimization and structure-activity relationship analysis of a new series of non-LBP AR antagonists. Induced-fit docking and molecular dynamics studies have been performed to establish a consistent hypothesis for the interaction of the new active molecule on the AR surface.

Optimization and pharmacological validation of a leukocyte migration assay in zebrafish larvae for the rapid in vivo bioactivity analysis of anti-inflammatory secondary metabolites

Over the past decade, zebrafish (Danio rerio) have emerged as an attractive model for in vivo drug discovery. In this study, we explore the suitability of zebrafish larvae to rapidly evaluate the anti-inflammatory activity of natural products (NPs) and medicinal plants used in traditional medicine for the treatment of inflammatory disorders. First, we optimized a zebrafish assay for leukocyte migration. Inflammation was induced in four days post-fertilization (dpf) zebrafish larvae by tail transection and co-incubation with bacterial lipopolysaccharides (LPS), resulting in a robust recruitment of leukocytes to the zone of injury. Migrating zebrafish leukocytes were detected in situ by myeloperoxidase (MPO) staining, and anti-inflammatory activity was semi-quantitatively scored using a standardized scale of relative leukocyte migration (RLM). Pharmacological validation of this optimized assay was performed with a panel of anti-inflammatory drugs, demonstrating a concentration-responsive inhibition of leukocyte migration for both steroidal and non-steroidal anti-inflammatory drugs (SAIDs and NSAIDs). Subsequently, we evaluated the bioactivity of structurally diverse NPs with well-documented anti-inflammatory properties. Finally, we further used this zebrafish-based assay to quantify the anti-inflammatory activity in the aqueous and methanolic extracts of several medicinal plants. Our results indicate the suitability of this LPS-enhanced leukocyte migration assay in zebrafish larvae as a front-line screening platform in NP discovery, including for the bioassay-guided isolation of anti-inflammatory secondary metabolites from complex NP extracts.

Fragment-based optimization of small molecule CXCL12 inhibitors for antagonizing the CXCL12/CXCR4 interaction

The chemokine CXCL12 and its G protein-coupled receptor (GPCR) CXCR4 are high-priority clinical targets because of their involvement in metastatic cancers (also implicated in autoimmune disease and cardiovascular disease). Because chemokines interact with two distinct sites to bind and activate their receptors, both the GPCRs and chemokines are potential targets for small molecule inhibition. A number of chemokines have been validated as targets for drug development, but virtually all drug discovery efforts focus on the GPCRs. However, all CXCR4 receptor antagonists with the exception of MSX-122 have failed in clinical trials due to unmanageable toxicities, emphasizing the need for alternative strategies to interfere with CXCL12/CXCR4-guided metastatic homing. Although targeting the relatively featureless surface of CXCL12 was presumed to be challenging, focusing efforts at the sulfotyrosine (sY) binding pockets proved successful for procuring initial hits. Using a hybrid structure-based in silico/NMR screening strategy, we recently identified a ligand that occludes the receptor recognition site. From this initial hit, we designed a small fragment library containing only nine tetrazole derivatives using a fragment-based and bioisostere approach to target the sY binding sites of CXCL12. Compound binding modes and affinities were studied by 2D NMR spectroscopy, X-ray crystallography, molecular docking and cell-based functional assays. Our results demonstrate that the sY binding sites are conducive to the development of high affinity inhibitors with better ligand efficiency (LE) than typical protein-protein interaction inhibitors (LE ≤ 0.24). Our novel tetrazole-based fragment 18 was identified to bind the sY21 site with a K(d) of 24 μM (LE = 0.30). Optimization of 18 yielded compound 25 which specifically inhibits CXCL12-induced migration with an improvement in potency over the initial hit 9. The fragment from this library that exhibited the highest affinity and ligand efficiency (11: K(d) = 13 μM, LE = 0.33) may serve as a starting point for development of inhibitors targeting the sY12 site.

Integration of Microfractionation, qNMR and zebrafish screening for the in vivo bioassay-guided isolation and quantitative bioactivity analysis of natural products

Natural products (NPs) are an attractive source of chemical diversity for small-molecule drug discovery. Several challenges nevertheless persist with respect to NP discovery, including the time and effort required for bioassay-guided isolation of bioactive NPs, and the limited biomedical relevance to date of in vitro bioassays used in this context. With regard to bioassays, zebrafish have recently emerged as an effective model system for chemical biology, allowing in vivo high-content screens that are compatible with microgram amounts of compound. For the deconvolution of the complex extracts into their individual constituents, recent progress has been achieved on several fronts as analytical techniques now enable the rapid microfractionation of extracts, and microflow NMR methods have developed to the point of allowing the identification of microgram amounts of NPs. Here we combine advanced analytical methods with high-content screening in zebrafish to create an integrated platform for microgram-scale, in vivo NP discovery. We use this platform for the bioassay-guided fractionation of an East African medicinal plant, Rhynchosia viscosa, resulting in the identification of both known and novel isoflavone derivatives with anti-angiogenic and anti-inflammatory activity. Quantitative microflow NMR is used both to determine the structure of bioactive compounds and to quantify them for direct dose-response experiments at the microgram scale. The key advantages of this approach are (1) the microgram scale at which both biological and analytical experiments can be performed, (2) the speed and the rationality of the bioassay-guided fractionation – generic for NP extracts of diverse origin – that requires only limited sample-specific optimization and (3) the use of microflow NMR for quantification, enabling the identification and dose-response experiments with only tens of micrograms of each compound. This study demonstrates that a complete in vivo bioassay-guided fractionation can be performed with only 20 mg of NP extract within a few days.

Modeling natural anti-inflammatory compounds by molecular topology

One of the main pharmacological problems today in the treatment of chronic inflammation diseases consists of the fact that anti-inflammatory drugs usually exhibit side effects. The natural products offer a great hope in the identification of bioactive lead compounds and their development into drugs for treating inflammatory diseases. Computer-aided drug design has proved to be a very useful tool for discovering new drugs and, specifically, Molecular Topology has become a good technique for such a goal. A topological-mathematical model, obtained by linear discriminant analysis, has been developed for the search of new anti-inflammatory natural compounds. An external validation obtained with the remaining compounds (those not used in building up the model), has been carried out. Finally, a virtual screening on natural products was performed and 74 compounds showed actual anti-inflammatory activity. From them, 54 had been previously described as anti-inflammatory in the literature. This can be seen as a plus in the model validation and as a reinforcement of the role of Molecular Topology as an efficient tool for the discovery of new anti-inflammatory natural compounds.