Addressing docking pose selection with structure-based deep learning: Recent advances, challenges and opportunities.

Molecular docking is a widely used technique in drug discovery to predict the binding mode of a given ligand to its target. However, the identification of the near-native binding pose in docking experiments still represents a challenging task as the scoring functions currently employed by docking programs are parametrized to predict the binding affinity, and, therefore, they often fail to correctly identify the ligand native binding conformation. Selecting the correct binding mode is crucial to obtaining meaningful results and to conveniently optimizing new hit compounds. Deep learning (DL) algorithms have been an area of a growing interest in this sense for their capability to extract the relevant information directly from the protein-ligand structure. Our review aims to present the recent advances regarding the development of DL-based pose selection approaches, discussing limitations and possible future directions. Moreover, a comparison between the performances of some classical scoring functions and DL-based methods concerning their ability to select the correct binding mode is reported. In this regard, two novel DL-based pose selectors developed by us are presented.

"DompeKeys": a set of novel substructure-based descriptors for efficient chemical space mapping, development and structural interpretation of machine learning models, and indexing of large databases.

The conversion of chemical structures into computer-readable descriptors, able to capture key structural aspects, is of pivotal importance in the field of cheminformatics and computer-aided drug design. Molecular fingerprints represent a widely employed class of descriptors; however, their generation process is time-consuming for large databases and is susceptible to bias. Therefore, descriptors able to accurately detect predefined structural fragments and devoid of lengthy generation procedures would be highly desirable. To meet additional needs, such descriptors should also be interpretable by medicinal chemists, and suitable for indexing databases with trillions of compounds. To this end, we developed-as integral part of EXSCALATE, Dompé's end-to-end drug discovery platform-the DompeKeys (DK), a new substructure-based descriptor set, which encodes the chemical features that characterize compounds of pharmaceutical interest. DK represent an exhaustive collection of curated SMARTS strings, defining chemical features at different levels of complexity, from specific functional groups and structural patterns to simpler pharmacophoric points, corresponding to a network of hierarchically interconnected substructures. Because of their extended and hierarchical structure, DK can be used, with good performance, in different kinds of applications. In particular, we demonstrate how they are very well suited for effective mapping of chemical space, as well as substructure search and virtual screening. Notably, the incorporation of DK yields highly performing machine learning models for the prediction of both compounds' activity and metabolic reaction occurrence. The protocol to generate the DK is freely available at https://dompekeys.exscalate.eu and is fully integrated with the Molecular Anatomy protocol for the generation and analysis of hierarchically interconnected molecular scaffolds and frameworks, thus providing a comprehensive and flexible tool for drug design applications.

Active Learning Approach for Guiding Site-of-Metabolism Measurement and Annotation.

The ability to determine and predict metabolically labile atom positions in a molecule (also called "sites of metabolism" or "SoMs") is of high interest to the design and optimization of bioactive compounds, such as drugs, agrochemicals, and cosmetics. In recent years, several in silico models for SoM prediction have become available, many of which include a machine-learning component. The bottleneck in advancing these approaches is the coverage of distinct atom environments and rare and complex biotransformation events with high-quality experimental data. Pharmaceutical companies typically have measured metabolism data available for several hundred to several thousand compounds. However, even for metabolism experts, interpreting these data and assigning SoMs are challenging and time-consuming. Therefore, a significant proportion of the potential of the existing metabolism data, particularly in machine learning, remains dormant. Here, we report on the development and validation of an active learning approach that identifies the most informative atoms across molecular data sets for SoM annotation. The active learning approach, built on a highly efficient reimplementation of SoM predictor FAME 3, enables experts to prioritize their SoM experimental measurements and annotation efforts on the most rewarding atom environments. We show that this active learning approach yields competitive SoM predictors while requiring the annotation of only 20% of the atom positions required by FAME 3. The source code of the approach presented in this work is publicly available.

MetaSpot: A General Approach for Recognizing the Reactive Atoms Undergoing Metabolic Reactions Based on the MetaQSAR Database.

The prediction of drug metabolism is attracting great interest for the possibility of discarding molecules with unfavorable ADME/Tox profile at the early stage of the drug discovery process. In this context, artificial intelligence methods can generate highly performing predictive models if they are trained by accurate metabolic data. MetaQSAR-based datasets were collected to predict the sites of metabolism for most metabolic reactions. The models were based on a set of structural, physicochemical, and stereo-electronic descriptors and were generated by the random forest algorithm. For each considered biotransformation, two types of models were developed: the first type involved all non-reactive atoms and included atom types among the descriptors, while the second type involved only non-reactive centers having the same atom type(s) of the reactive atoms. All the models of the first type revealed very high performances; the models of the second type show on average worst performances while being almost always able to recognize the reactive centers; only conjugations with glucuronic acid are unsatisfactorily predicted by the models of the second type. Feature evaluation confirms the major role of lipophilicity, self-polarizability, and H-bonding for almost all considered reactions. The obtained results emphasize the possibility of recognizing the sites of metabolism by classification models trained on MetaQSAR database. The two types of models can be synergistically combined since the first models identify which atoms can undergo a given metabolic reactions, while the second models detect the truly reactive centers. The generated models are available as scripts for the VEGA program.

The VEGA web service: multipurpose online tools for molecular modelling and docking analyses.

The paper presents the VEGA Online web service, which includes a set of freely available tools deriving from the development of the VEGA suite of programs. In detail, the paper is focused on two tools: the VEGA Web Edition (WE) and the Score tool. The former is a versatile file format converter including relevant features for 2D/3D conversion, for surface mapping and for editing/preparing input files. The Score application allows rescoring docking poses and in particular includes the MLP Interactions Scores (MLPInS) for describing hydrophobic interactions. To the best of our knowledge, this web service is the only available resource by which one can calculate both the virtual log P of a given input molecule according to the MLP approach plus the corresponding MLP surface.

Ensemble of structure and ligand-based classification models for hERG liability profiling.

Drug-induced cardiotoxicity represents one of the most critical safety concerns in the early stages of drug development. The blockade of the human ether-à-go-go-related potassium channel (hERG) is the most frequent cause of cardiotoxicity, as it is associated to long QT syndrome which can lead to fatal arrhythmias. Therefore, assessing hERG liability of new drugs candidates is crucial to avoid undesired cardiotoxic effects. In this scenario, computational approaches have emerged as useful tools for the development of predictive models able to identify potential hERG blockers. In the last years, several efforts have been addressed to generate ligand-based (LB) models due to the lack of experimental structural information about hERG channel. However, these methods rely on the structural features of the molecules used to generate the model and often fail in correctly predicting new chemical scaffolds. Recently, the 3D structure of hERG channel has been experimentally solved enabling the use of structure-based (SB) strategies which may overcome the limitations of the LB approaches. In this study, we compared the performances achieved by both LB and SB classifiers for hERG-related cardiotoxicity developed by using Random Forest algorithm and employing a training set containing 12789 hERG binders. The SB models were trained on a set of scoring functions computed by docking and rescoring calculations, while the LB classifiers were built on a set of physicochemical descriptors and fingerprints. Furthermore, models combining the LB and SB features were developed as well. All the generated models were internally validated by ten-fold cross-validation on the TS and further verified on an external test set. The former revealed that the best performance was achieved by the LB model, while the model combining the LB and the SB attributes displayed the best results when applied on the external test set highlighting the usefulness of the integration of LB and SB features in correctly predicting unseen molecules. Overall, our predictive models showed satisfactory performances providing new useful tools to filter out potential cardiotoxic drug candidates in the early phase of drug discovery.

Obesity and Type 2 Diabetes: Adiposopathy as a Triggering Factor and Therapeutic Options.

Obesity and type 2 diabetes (T2DM) are major public health concerns associated with serious morbidity and increased mortality. Both obesity and T2DM are strongly associated with adiposopathy, a term that describes the pathophysiological changes of the adipose tissue. In this review, we have highlighted adipose tissue dysfunction as a major factor in the etiology of these conditions since it promotes chronic inflammation, dysregulated glucose homeostasis, and impaired adipogenesis, leading to the accumulation of ectopic fat and insulin resistance. This dysfunctional state can be effectively ameliorated by the loss of at least 15% of body weight, that is correlated with better glycemic control, decreased likelihood of cardiometabolic disease, and an improvement in overall quality of life. Weight loss can be achieved through lifestyle modifications (healthy diet, regular physical activity) and pharmacotherapy. In this review, we summarized different effective management strategies to address weight loss, such as bariatric surgery and several classes of drugs, namely metformin, GLP-1 receptor agonists, amylin analogs, and SGLT2 inhibitors. These drugs act by targeting various mechanisms involved in the pathophysiology of obesity and T2DM, and they have been shown to induce significant weight loss and improve glycemic control in obese individuals with T2DM.

Target Prediction by Multiple Virtual Screenings: Analyzing the SARS-CoV-2 Phenotypic Screening by the Docking Simulations Submitted to the MEDIATE Initiative.

Phenotypic screenings are usually combined with deconvolution techniques to characterize the mechanism of action for the retrieved hits. These studies can be supported by various computational analyses, although docking simulations are rarely employed. The present study aims to assess if multiple docking calculations can prove successful in target prediction. In detail, the docking simulations submitted to the MEDIATE initiative are utilized to predict the viral targets involved in the hits retrieved by a recently published cytopathic screening. Multiple docking results are combined by the EFO approach to develop target-specific consensus models. The combination of multiple docking simulations enhances the performances of the developed consensus models (average increases in EF1% value of 40% and 25% when combining three and two docking runs, respectively). These models are able to propose reliable targets for about half of the retrieved hits (31 out of 59). Thus, the study emphasizes that docking simulations might be effective in target identification and provide a convincing validation for the collaborative strategies that inspire the MEDIATE initiative. Disappointingly, cross-target and cross-program correlations suggest that common scoring functions are not specific enough for the simulated target.

Computational Insights into the Sequence-Activity Relationships of the NGF(1-14) Peptide by Molecular Dynamics Simulations.

The Nerve Growth Factor (NGF) belongs to the neurothrophins protein family involved in the survival of neurons in the nervous system. The interaction of NGF with its high-affinity receptor TrkA mediates different cellular pathways related to Alzheimer's disease, pain, ocular dysfunction, and cancer. Therefore, targeting NGF-TrkA interaction represents a valuable strategy for the development of new therapeutic agents. In recent years, experimental studies have revealed that peptides belonging to the N-terminal domain of NGF are able to partly mimic the biological activity of the whole protein paving the way towards the development of small peptides that can selectively target specific signaling pathways. Hence, understanding the molecular basis of the interaction between the N-terminal segment of NGF and TrkA is fundamental for the rational design of new peptides mimicking the NGF N-terminal domain. In this study, molecular dynamics simulation, binding free energy calculations and per-residue energy decomposition analysis were combined in order to explore the molecular recognition pattern between the experimentally active NGF(1-14) peptide and TrkA. The results highlighted the importance of His4, Arg9 and Glu11 as crucial residues for the stabilization of NGF(1-14)-TrkA interaction, thus suggesting useful insights for the structure-based design of new therapeutic peptides able to modulate NGF-TrkA interaction.

Extensive Sampling of Molecular Dynamics Simulations to Identify Reliable Protein Structures for Optimized Virtual Screening Studies: The Case of the hTRPM8 Channel.

(1) Background: Virtual screening campaigns require target structures in which the pockets are properly arranged for binding. Without these, MD simulations can be used to relax the available target structures, optimizing the fine architecture of their binding sites. Among the generated frames, the best structures can be selected based on available experimental data. Without experimental templates, the MD trajectories can be filtered by energy-based criteria or sampled by systematic analyses. (2) Methods: A blind and methodical analysis was performed on the already reported MD run of the hTRPM8 tetrameric structures; a total of 50 frames underwent docking simulations by using a set of 1000 ligands including 20 known hTRPM8 modulators. Docking runs were performed by LiGen program and involved the frames as they are and after optimization by SCRWL4.0. For each frame, all four monomers were considered. Predictive models were developed by the EFO algorithm based on the sole primary LiGen scores. (3) Results: On average, the MD simulation progressively enhances the performance of the extracted frames, and the optimized structures perform better than the non-optimized frames (EF1% mean: 21.38 vs. 23.29). There is an overall correlation between performances and volumes of the explored pockets and the combination of the best performing frames allows to develop highly performing consensus models (EF1% = 49.83). (4) Conclusions: The systematic sampling of the entire MD run provides performances roughly comparable with those previously reached by using rationally selected frames. The proposed strategy appears to be helpful when the lack of experimental data does not allow an easy selection of the optimal structures for docking simulations. Overall, the reported docking results confirm the relevance of simulating all the monomers of an oligomer structure and emphasize the efficacy of the SCRWL4.0 method to optimize the protein structures for docking calculations.

A multiscale approach to predict the binding mode of metallo beta-lactamase inhibitors.

Antibiotic resistance is a major threat to global public health. ß-lactamases, which catalyze breakdown of ß-lactam antibiotics, are a principal cause. Metallo ß-lactamases (MBLs) represent a particular challenge because they hydrolyze almost all ß-lactams and to date no MBL inhibitor has been approved for clinical use. Molecular simulations can aid drug discovery, for example, predicting inhibitor complexes, but empirical molecular mechanics (MM) methods often perform poorly for metalloproteins. Here we present a multiscale approach to model thiol inhibitor binding to IMP-1, a clinically important MBL containing two catalytic zinc ions, and predict the binding mode of a 2-mercaptomethyl thiazolidine (MMTZ) inhibitor. Inhibitors were first docked into the IMP-1 active site, testing different docking programs and scoring functions on multiple crystal structures. Complexes were then subjected to molecular dynamics (MD) simulations and subsequently refined through QM/MM optimization with a density functional theory (DFT) method, B3LYP/6-31G(d), increasing the accuracy of the method with successive steps. This workflow was tested on two IMP-1:MMTZ complexes, for which it reproduced crystallographically observed binding, and applied to predict the binding mode of a third MMTZ inhibitor for which a complex structure was crystallographically intractable. We also tested a 12-6-4 nonbonded interaction model in MD simulations and optimization with a SCC-DFTB QM/MM approach. The results show the limitations of empirical models for treating these systems and indicate the need for higher level calculations, for example, DFT/MM, for reliable structural predictions. This study demonstrates a reliable computational pipeline that can be applied to inhibitor design for MBLs and other zinc-metalloenzyme systems.

Erratum: Altomare et al. In-Depth AGE and ALE Profiling of Human Albumin in Heart Failure: Ex Vivo Studies. Antioxidants 2021, 10, 358.

The authors of this paper [...].

MetaClass, a Comprehensive Classification System for Predicting the Occurrence of Metabolic Reactions Based on the MetaQSAR Database.

(1) Background: Machine learning algorithms are finding fruitful applications in predicting the ADME profile of new molecules, with a particular focus on metabolism predictions. However, the development of comprehensive metabolism predictors is hampered by the lack of highly accurate metabolic resources. Hence, we recently proposed a manually curated metabolic database (MetaQSAR), the level of accuracy of which is well suited to the development of predictive models. (2) Methods: MetaQSAR was used to extract datasets to predict the metabolic reactions subdivided into major classes, classes and subclasses. The collected datasets comprised a total of 3788 first-generation metabolic reactions. Predictive models were developed by using standard random forest algorithms and sets of physicochemical, stereo-electronic and constitutional descriptors. (3) Results: The developed models showed satisfactory performance, especially for hydrolyses and conjugations, while redox reactions were predicted with greater difficulty, which was reasonable as they depend on many complex features that are not properly encoded by the included descriptors. (4) Conclusions: The generated models allowed a precise comparison of the propensity of each metabolic reaction to be predicted and the factors affecting their predictability were discussed in detail. Overall, the study led to the development of a freely downloadable global predictor, MetaClass, which correctly predicts 80% of the reported reactions, as assessed by an explorative validation analysis on an external dataset, with an overall MCC = 0.44.

MetaTREE, a Novel Database Focused on Metabolic Trees, Predicts an Important Detoxification Mechanism: The Glutathione Conjugation.

(1) Background: Data accuracy plays a key role in determining the model performances and the field of metabolism prediction suffers from the lack of truly reliable data. To enhance the accuracy of metabolic data, we recently proposed a manually curated database collected by a meta-analysis of the specialized literature (MetaQSAR). Here we aim to further increase data accuracy by focusing on publications reporting exhaustive metabolic trees. This selection should indeed reduce the number of false negative data. (2) Methods: A new metabolic database (MetaTREE) was thus collected and utilized to extract a dataset for metabolic data concerning glutathione conjugation (MT-dataset). After proper pre-processing, this dataset, along with the corresponding dataset extracted from MetaQSAR (MQ-dataset), was utilized to develop binary classification models using a random forest algorithm. (3) Results: The comparison of the models generated by the two collected datasets reveals the better performances reached by the MT-dataset (MCC raised from 0.63 to 0.67, sensitivity from 0.56 to 0.58). The analysis of the applicability domain also confirms that the model based on the MT-dataset shows a more robust predictive power with a larger applicability domain. (4) Conclusions: These results confirm that focusing on metabolic trees represents a convenient approach to increase data accuracy by reducing the false negative cases. The encouraging performances shown by the models developed by the MT-dataset invites to use of MetaTREE for predictive studies in the field of xenobiotic metabolism.

In-Depth AGE and ALE Profiling of Human Albumin in Heart Failure: Ex Vivo Studies.

Advanced glycation end-products (AGEs) and advanced lipoxidation end-products (ALEs), particularly carboxymethyl-lysine (CML), have been largely proposed as factors involved in the establishment and progression of heart failure (HF). Despite this evidence, the current literature lacks the comprehensive identification and characterization of the plasma AGEs/ALEs involved in HF (untargeted approach). This work provides the first ex vivo high-resolution mass spectrometry (HR-MS) profiling of AGEs/ALEs occurring in human serum albumin (HSA), the most abundant protein in plasma, characterized by several nucleophilic sites and thus representing the main protein substrate for AGE/ALE formation. A set of AGE/ALE adducts in pooled HF-HSA samples was defined, and a semi-quantitative analysis was carried out in order to finally select those presenting in increased amounts in the HF samples with respect to the control condition. These adducts were statistically confirmed by monitoring their content in individual HF samples by applying a targeted approach. Selected AGEs/ALEs proved to be mostly CML derivatives on Lys residues (i.e., CML-Lys12, CML-Lys378, CML-Lys402), and one deoxy-fructosyl derivative on the Lys 389 (DFK-Lys 389). The nature of CML adducts was finally confirmed using immunological methods and in vitro production of such adducts further confirmed by mass spectrometry.

Combining Different Docking Engines and Consensus Strategies to Design and Validate Optimized Virtual Screening Protocols for the SARS-CoV-2 3CL Protease.

The 3CL-Protease appears to be a very promising medicinal target to develop anti-SARS-CoV-2 agents. The availability of resolved structures allows structure-based computational approaches to be carried out even though the lack of known inhibitors prevents a proper validation of the performed simulations. The innovative idea of the study is to exploit known inhibitors of SARS-CoV 3CL-Pro as a training set to perform and validate multiple virtual screening campaigns. Docking simulations using four different programs (Fred, Glide, LiGen, and PLANTS) were performed investigating the role of both multiple binding modes (by binding space) and multiple isomers/states (by developing the corresponding isomeric space). The computed docking scores were used to develop consensus models, which allow an in-depth comparison of the resulting performances. On average, the reached performances revealed the different sensitivity to isomeric differences and multiple binding modes between the four docking engines. In detail, Glide and LiGen are the tools that best benefit from isomeric and binding space, respectively, while Fred is the most insensitive program. The obtained results emphasize the fruitful role of combining various docking tools to optimize the predictive performances. Taken together, the performed simulations allowed the rational development of highly performing virtual screening workflows, which could be further optimized by considering different 3CL-Pro structures and, more importantly, by including true SARS-CoV-2 3CL-Pro inhibitors (as learning set) when available.

The VEGA suite of programs: an versatile platform for cheminformatics and drug design projects.

The purpose of the article is to offer an overview of the latest release of the VEGA suite of programs. This software has been constantly developed and freely released during the last 20 years and has now reached a significant diffusion and technology level as confirmed by the about 22 500 registered users. While being primarily developed for drug design studies, the VEGA package includes cheminformatics and modeling features, which can be fruitfully utilized in various contexts of the computational chemistry. To offer a glimpse of the remarkable potentials of the software, some examples of the implemented features in the cheminformatics field and for structure-based studies are discussed. Finally, the flexible architecture of the VEGA program which can be expanded and customized by plug-in technology or scripting languages will be described focusing attention on the HyperDrive library including highly optimized functions. AVAILABILITY AND IMPLEMENTATION: The VEGA suite of programs and the source code of the VEGA command-line version are available free of charge for non-profit organizations at http://www.vegazz.net.

Health- and skill-related physical fitness profile of Brazilian children and adolescents: a systematic review / Perfil da aptidão física relacionada à saúde e ao desempenho motor de crianças e adolescentes brasileiros: uma revisão sistemática

Since 1994, the Projeto Esporte Brasil (PROESP-Br) battery tests has been used to evaluate health- and skill-related physical fitness among aged 6-17 Brazilian schoolchildren. The aim of this study was to delineate the Brazilian children and youth's physical fitness profile from a systematic review over studies that used the PROESP-Br proposal. The search was carried at PubMed, ScienceDirect, Lilacs, SciELO and Google Scholar. Original studies published between 1994 and 2017 about physical fitness (health and/or motor performance) with schoolchildren (children and/or adolescents) that used the PROESP-Br battery test were included. A total of 13.582 participants were evaluated to health-related fitness and 276 to skill-related fitness from 18 included studies. The methodological quality was evaluated using the Newcastle-Ottawa quality assessment scale adapted version. The results show that 27-30% of youngsters are at health "risk zone" for Body Mass Index (BMI), 70% for cardiorespiratory fitness (CRF), 50% and 65% for flexibility (FLEX) and muscular strength (MST), respectively. The data concerning skill-related fitness were inconsistent. In summary, the results suggest that Brazilian children and adolescents have low cardiovascular health level (BMI/CRF), mainly regarding CRF, and low muscle health level (FLEX/MST). We emphasize that the lack of studies regarding skill-related fitness, make it impossible to describe the profile of the components of this construct

O Projeto Esporte Brasil (PROESP-Br) propõe, desde 1994, uma bateria de medidas e testes para avaliação de escolares entre seis e 17 anos com o objetivo de delinear o perfil de crianças e jovens brasileiros no que se refere a aptidão física relacionada à saúde e ao desempenho motor. O objetivo deste estudo foi delinear o perfil da aptidão física de crianças e jovens brasileiros a partir de uma revisão sistemática da literatura sobre artigos que utilizaram da proposta do PROESP-Br. A busca foi realizada na PubMed, ScienceDirect, Lilacs, SciELO e Google Acadêmico. Estudos originais publicados entre 1994 e 2017 acerca da aptidão física (saúde e/ou desempenho motor) de escolares (crianças e/ou adolescentes) que utilizaram a bateria de testes do PROESP-Br foram incluídos. Um total de 13.582 sujeitos foram avaliados quanto a saúde e 276 quanto ao desempenho motor nos 18 estudos incluídos. A qualidade metodológica foi avaliada uma versão adaptada da Newcastle-Ottawa quality assessment scale. Os resultados evidenciam que 27-30% dos jovens estão na "zona de risco" à saúde para o Índice de Massa Corporal (IMC), 70% para a aptidão cardiorrespiratória (ApC) e 50 e 65% para flexibilidade e força muscular localizada (FML), respectivamente. Os dados sobre o desempenho motor são inconsistentes nesta revisão de literatura. Em síntese, os resultados indicam baixos níveis de saúde cardiovascular (IMC/ApC), principalmente quanto à ApC, assim como baixos níveis de saúde musculoesquelética (flexibilidade/FML) dos jovens. Ressalta-se a escassez de estudos quanto ao desempenho motor impossibilitando delinear o perfil dos componentes deste construto

Repositioning Dequalinium as Potent Muscarinic Allosteric Ligand by Combining Virtual Screening Campaigns and Experimental Binding Assays.

Structure-based virtual screening is a truly productive repurposing approach provided that reliable target structures are available. Recent progresses in the structural resolution of the G-Protein Coupled Receptors (GPCRs) render these targets amenable for structure-based repurposing studies. Hence, the present study describes structure-based virtual screening campaigns with a view to repurposing known drugs as potential allosteric (and/or orthosteric) ligands for the hM2 muscarinic subtype which was indeed resolved in complex with an allosteric modulator thus allowing a precise identification of this binding cavity. First, a docking protocol was developed and optimized based on binding space concept and enrichment factor optimization algorithm (EFO) consensus approach by using a purposely collected database including known allosteric modulators. The so-developed consensus models were then utilized to virtually screen the DrugBank database. Based on the computational results, six promising molecules were selected and experimentally tested and four of them revealed interesting affinity data; in particular, dequalinium showed a very impressive allosteric modulation for hM2. Based on these results, a second campaign was focused on bis-cationic derivatives and allowed the identification of other two relevant hM2 ligands. Overall, the study enhances the understanding of the factors governing the hM2 allosteric modulation emphasizing the key role of ligand flexibility as well as of arrangement and delocalization of the positively charged moieties.

Towards the Inhibition of Protein-Protein Interactions (PPIs) in STAT3: Insights into a New Class of Benzothiadiazole Derivatives.

Signal transducer and activator of transcription 3 (STAT3) is a validated anticancer target due to the relationship between its constitutive activation and malignant tumors. Through a virtual screening approach on the STAT3-SH2 domain, 5,6-dimethyl-1H,3H-2,1,3-benzothiadiazole-2,2-dioxide (1) was identified as a potential STAT3 inhibitor. Some benzothiadiazole derivatives were synthesized by employing a versatile methodology, and they were tested by an AlphaScreen-based assay. Among them, benzosulfamide 1 showed a significant activity with an IC50 = 15.8 ± 0.6 µM as a direct STAT3 inhibitor. Notably, we discovered that compound 1 was also able to interact with cysteine residues located around the SH2 domain. By applying mass spectrometry, liquid chromatography, NMR, and UV spectroscopy, an in-depth investigation was carried out, shedding light on its intriguing and unexpected mechanism of interaction.