Thermolysin is a bacterial proteolytic enzyme, considered by many authors as a pharmacological and biological model of other mammalian enzymes, with similar structural characteristics, such as angiotensin converting enzyme and neutral endopeptidase. Inhibitors of these enzymes are considered therapeutic targets for common diseases, such as hypertension and heart failure. In this report, a mathematical model of Multiple Linear Regression, for ordinary least squares, and genetic algorithm, for selection of variables, are developed and implemented in QSARINS software, with appropriate parameters for its fitting. The model is extensively validated according to OECD standards, so that its robustness, stability, low correlation of descriptors and good predictive power are proven. In addition, it is found that the model fit is not the product of a random correlation. Two possible outliers are identified in the model application domain but, in a molecular docking study, they show good activity, so we decide to keep both in our database. Finally, 141 and 69 compounds (2.5⩽pKi<3.5 and pKi⩾3.5, respectively) are identified as potential Thermolysin inhibitors, concluding that the proposed computational tools are an efficient method for the identification of new drugs that could inhibit this enzyme

NADPH oxidase (NOX2) is responsible for reactive oxygen species (ROS) production in neutrophils and has been recognized as a key mediator in inflammatory and cardiovascular pathologies. Nevertheless, there is a lack of specific NOX2 pharmacological inhibitors. In medicinal chemistry, heterocyclic compounds are essential scaffolds for drug design, and among them, indole is a very versatile pharmacophore. We tested the hypothesis that indole heteroaryl-acrylonitrile derivatives may serve as NOX2 inhibitors by evaluating the capacity of 19 of these molecules to inhibit NOX2-derived ROS production in human neutrophils (HL-60 cells). Of these compounds, C6 and C14 exhibited concentration-dependent inhibition of NOX2 (IC₅₀~1 µM). These molecules also reduced NOX2-derived oxidative stress in cardiomyocytes and prevented cardiac damage induced by ischemia-reperfusion. Compound C6 significantly reduced the membrane translocation of p47^phox, a cytosolic subunit that is required for NOX2 activation. Molecular docking analyses of the binding modes of these molecules with p47^phox indicated that C6 and C14 interact with specific residues in the inner part of the groove of p47^phox, the binding cavity for p22^phox. This combination of methods showed that novel indole heteroaryl acrylonitriles represent interesting lead compounds for developing specific and potent NOX2 inhibitors.

In this work we report the first attempt to study the effect of activity cliffs over the generalization ability of machine learning (ML) based QSAR classifiers, using as study case a previously reported diverse and noisy dataset focused on drug induced liver injury (DILI) and more than 40 ML classification algorithms. Here, the hypothesis of structure-activity relationship (SAR) continuity restoration by activity cliffs removal is tested as a potential solution to overcome such limitation. Previously, a parallelism was established between activity cliffs generators (ACGs) and instances that should be misclassified (ISMs), a related concept from the field of machine learning. Based on this concept we comparatively studied the classification performance of multiple machine learning classifiers as well as the consensus classifier derived from predictive classifiers obtained from training sets including or excluding ACGs. The influence of the removal of ACGs from the training set over the virtual screening performance was also studied for the respective consensus classifiers algorithms. In general terms, the removal of the ACGs from the training process slightly decreased the overall accuracy of the ML classifiers and multi-classifiers, improving their sensitivity (the weakest feature of ML classifiers trained with ACGs) but decreasing their specificity. Although these results do not support a positive effect of the removal of ACGs over the classification performance of ML classifiers, the “balancing effect” of ACG removal demonstrated to positively influence the virtual screening performance of multi-classifiers based on valid base ML classifiers. Specially, the early recognition ability was significantly favored after ACGs removal. The results presented and discussed in this work represent the first step towards the application of a remedial solution to the activity cliffs problem in QSAR studies.

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

In this study, a methodology is proposed, combining ligand- and structure-based virtual screening tools, for the identification of phosphorus-containing compounds as inhibitors of zinc metalloproteases. First, we use Dragon molecular descriptors to develop a Linear Discriminant Analysis classification model, which is widely validated according to the OECD principles. This model is simple, robust, stable and has good discriminating power. Furthermore, it has a defined applicability domain and it is used for virtual screening of the DrugBank database. Second, docking experiments are carried out on the identified compounds that showed good binding energies to the enzyme thermolysin. Considering the potential toxicity of phosphorus-containing compounds, their toxicological profile is evaluated according to Protox II. Of the five molecules evaluated, two show carcinogenic and mutagenic potential at small LD50, not recommended as drugs, while three of them are classified as non-toxic, and could constitute a starting point for the development of new vasoactive metalloprotease inhibitor drugs. According to molecular dynamics simulation, two of them show stable interactions with the active site maintaining coordination with the metal. A high agreement is evident between QSAR, docking and molecular dynamics results, demonstrating the potentialities of the combination of these tools.