Search

The management of patients with type 2 diabetes mellitus (T2DM) is shifting from cardio-centric to weight-centric or, even better, adipose-centric treatments. Considering the downsides of multidrug therapies and the relevance of dipeptidyl peptidase IV (DPP IV) and carbonic anhydrases (CAs II and V) in T2DM and in the weight loss, we report a new class of multitarget ligands targeting the mentioned enzymes. We started from the known α

Bacterial resistance to antibiotics has been long recognized as a priority to address for human health. Among all micro-organisms, the so-called multi-drug resistant (MDR) bacteria, which are resistant to most, if not all drugs in our current arsenal, are particularly worrisome. The World Health Organization has prioritized the ESKAPE ( Enterococcus faecium , Staphylococcus aureus , Klebsiella pneumoniae , Acinetobacter baumannii , Pseudomonas aeruginosa and Enterobacter species) pathogens, which include four Gram-negative bacterial species. In these bacteria, active extrusion of antimicrobial compounds out of the cell by means of ‘molecular guns’ known as efflux pumps is a main determinant of MDR phenotypes. The resistance-nodulation-cell division (RND) superfamily of efflux pumps connecting the inner and outer membrane in Gram-negative bacteria is crucial to the onset of MDR and virulence, as well as biofilm formation. Thus, understanding the molecular basis of the interaction of antibiotics and inhibitors with these pumps is key to the design of more effective therapeutics. With the aim to contribute to this challenge, and complement and inspire experimental research, in silico studies on RND efflux pumps have flourished in recent decades. Here, we review a selection of such investigations addressing the main determinants behind the polyspecificity of these pumps, the mechanisms of substrate recognition, transport and inhibition, as well as the relevance of their assembly for proper functioning, and the role of protein–lipid interactions. The journey will end with a perspective on the role of computer simulations in addressing the challenges posed by these beautifully complex machineries and in supporting the fight against the spread of MDR bacteria.

Putative binding modes (BMs) of quinolones to the bacterial efflux transporter MexB were identified. Multiple interaction patterns are possible, supporting the hypothesis that substrates oscillate between different BMs with similar affinity.

Zinc metalloproteins are ubiquitous, with protein zinc centers of structural and functional importance, involved in interactions with ligands and substrates and often of pharmacological interest. Biomolecular simulations are increasingly prominent in investigations of protein structure, dynamics, ligand interactions, and catalysis, but zinc poses a particular challenge, in part because of its versatile, flexible coordination. A computational workflow generating reliable models of ligand complexes of biological zinc centers would find broad application. Here, we evaluate the ability of alternative treatments, using (nonbonded) molecular mechanics (MM) and quantum mechanics/molecular mechanics (QM/MM) at semiempirical (DFTB3) and density functional theory (DFT) levels of theory, to describe the zinc centers of ligand complexes of six metalloenzyme systems differing in coordination geometries, zinc stoichiometries (mono- and dinuclear), and the nature of interacting groups (specifically the presence of zinc-sulfur interactions). MM molecular dynamics (MD) simulations can overfavor octahedral geometries, introducing additional water molecules to the zinc coordination shell, but this can be rectified by subsequent semiempirical (DFTB3) QM/MM MD simulations. B3LYP/MM geometry optimization further improved the accuracy of the description of coordination distances, with the overall effectiveness of the approach depending upon factors, including the presence of zinc-sulfur interactions that are less well described by semiempirical methods. We describe a workflow comprising QM/MM MD using DFTB3 followed by QM/MM geometry optimization using DFT (e.g., B3LYP) that well describes our set of zinc metalloenzyme complexes and is likely to be suitable for creating accurate models of zinc protein complexes when structural information is more limited.

TRPM8 is a calcium ion channel that is activated by multiple factors, such as temperature, voltage, pressure, and osmolality. It is a therapeutic target for anticancer drug development, and its modulators can be utilized for several pathological conditions. Here, we present a cryo-electron microscopy structure of a human TRPM8 channel in the closed state that was solved at 2.7 Å resolution. Our structure comprises the most complete model of the N-terminal pre-melastatin homology region. We also visualized several ligands that are bound by the protein and modeled how the human channel interacts with icilin. Analyses of pore helices in available TRPM structures showed that all these structures can be grouped into different closed, desensitized and open state conformations based on the register of the pore helix S6 which positions particular amino acid residues at the channel constriction. Our structure is the first for the human TRPM8 protein and it is among the most complete and the highest resolution structures of any TRPM8 channel available.

In the search for novel bitopic compounds targeting the dopamine D3 receptor (D3R), the N-(2,3-dichlorophenyl)piperazine nucleus (primary pharmacophore) has been linked to the 6,6- or 5,5-diphenyl-1,4-dioxane-2-carboxamide or the 1,4-benzodioxane-2-carboxamide scaffold (secondary pharmacophore) by an unsubstituted or 3-F-/3-OH-substituted butyl chain. This scaffold hybridization strategy led to the discovery of potent D3R-selective or multitarget ligands potentially useful for central nervous system disorders. In particular, the 6,6-diphenyl-1,4-dioxane derivative 3 showed a D3R-preferential profile, while an interesting multitarget behavior has been highlighted for the 5,5-diphenyl-1,4-dioxane and 1,4-benzodioxane derivatives 6 and 9, respectively, which displayed potent D2R antagonism, 5-HT1AR and D4R agonism, as well as potent D3R partial agonism. They also behaved as low-potency 5-HT2AR antagonists and 5-HT2CR partial agonists. Such a profile might be a promising starting point for the discovery of novel antipsychotic agents.

Cyclo(His-Pro) (CHP) is a cyclic dipeptide which is endowed with favorable pharmacokinetic properties combined with a variety of biological activities. CHP is found in a number of protein-rich foods and dietary supplements. While being stable at physiological pH, CHP can open yielding two symmetric dipeptides (His-Pro, Pro-His), the formation of which might be particularly relevant from dietary CHP due to the gastric acidic environment. The antioxidant and protective CHP properties were repeatedly reported although the non-enzymatic mechanisms were scantly investigated. The CHP detoxifying activity towards α,β unsaturated carbonyls was never investigated in detail, although its open dipeptides might be effective as already observed for histidine containing dipeptides. Hence, this study investigated the scavenging properties of TRH, CHP and its open derivatives towards 4-hydroxy-2-nonenal. The obtained results revealed that Pro-His possesses a marked activity and is more reactive than l-carnosine. As investigated by DFT calculations, the enhanced reactivity can be ascribed to the greater electrophilicity of the involved iminium intermediate. These findings emphasize that the primary amine (as seen in l-carnosine) can be replaced by secondary amines with beneficial effects on the quenching mechanisms. Serum stability of the tested peptides was also evaluated, showing that Pro-His is characterized by a greater stability than l-carnosine. Docking simulations suggested that its hydrolysis can be catalyzed by serum carnosinase. Altogether, the reported results evidence that the antioxidant CHP properties can be also due to the detoxifying activity of its open dipeptides, which might be thus responsible for the beneficial effects induced by CHP containing food.

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, e.g. 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 non-bonded 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, e.g. DFT/MM, for reliable structural predictions. This work demonstrates a reliable computational pipeline that can be applied to inhibitor design for MBLs and other zinc-metalloenzyme systems. This article is protected by copyright. All rights reserved.

In the search for novel bitopic compounds targeting the dopamine D3 receptor (D3R), the N-(2,3-dichlorophenyl)piperazine nucleus (primary pharmacophore) has been linked to the 6,6- or 5,5-diphenyl-1,4-dioxane-2-carboxamide or the 1,4-benzodioxane-2-carboxamide scaffold (secondary pharmacophore) by an unsubstituted or 3-F-/3-OH-substituted butyl chain. This scaffold hybridization strategy led to the discovery of potent D3R-selective or multitarget ligands potentially useful for central nervous system disorders. In particular, the 6,6-diphenyl-1,4-dioxane derivative 3 showed a D3R-preferential profile, while an interesting multitarget behavior has been highlighted for the 5,5-diphenyl-1,4-dioxane and 1,4-benzodioxane derivatives 6 and 9, respectively, which displayed potent D2R antagonism, 5-HT1AR and D4R agonism, as well as potent D3R partial agonism. They also behaved as low-potency 5-HT2AR antagonists and 5-HT2CR partial agonists. Such a profile might be a promising starting point for the discovery of novel antipsychotic agents.

Crop disease management often implies repeated application of fungicides. However, the increasing emergence of fungicide-resistant pathogens requires their rotation or combined use. Tank-mix combinations using fungicides with different modes of action are often hard to manage by farmers. An alternative and unexploited strategy are bifunctional fungicides, i.e. compounds resulting from conjugation of the pharmacophores of fungicides with different mechanisms of action. In this paper we describe a new approach to antifungal treatments based on the synthesis of dual agents, obtained by merging the strobilurin and succinate dehydrogenase inhibitor pharmacophores into a new entity. The compounds were tested against important fungal plant pathogens and showed good inhibition of Pyricularia oryzae and Sclerotinia sclerotiorum with activity comparable to commercial fungicides. The inhibition of the cytochrome bc1 and the succinate dehydrogenase enzyme activity confirmed that the new molecules are endowed with a dual mechanism of action. These results were further supported by molecular modelling which showed that selected compounds form stable complexes with both cytochrome b subunit and succinate dehydrogenase enzyme. This work can be considered an important first step towards the development of novel dual-action agents with optimized structure and improved interaction with the targets.

Hybrid molecules targeting simultaneously DNA polymerase α (POLA1) and histone deacetylases (HDACs) were designed and synthesized to exploit a potential synergy of action. Among a library of screened molecules, MIR002 and GEM144 showed antiproliferative activity at nanomolar concentrations on a panel of human solid and haematological cancer cell lines. In vitro functional assays confirmed that these molecules inhibited POLA1 primer extension activity, as well as HDAC11. Molecular docking studies also supported these findings. Mechanistically, MIR002 and GEM144 induced acetylation of p53, activation of p21, G1/S cell cycle arrest, and apoptosis. Oral administration of these inhibitors confirmed their antitumor activity in in vivo models. In human non-small cancer cell (H460) xenografted in nude mice MIR002 at 50 mg/kg, Bid (qd × 5 × 3w) inhibited tumor growth (TGI = 61%). More interestingly, in POLA1 inhibitor resistant cells (H460-R9A), the in vivo combination of MIR002 with cisplatin showed an additive antitumor effect with complete disappearance of tumor masses in two animals at the end of the treatment. Moreover, in two human orthotopic malignant pleural mesothelioma xenografts (MM473 and MM487), oral treatments with MIR002 and GEM144 confirmed their significant antitumor activity (TGI = 72–77%). Consistently with recent results that have shown an inverse correlation between POLA1 expression and type I interferon levels, MIR002 significantly upregulated interferon-α in immunocompetent mice.

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.

Atypical retinoids (AR) or retinoid-related molecules (RRMs) represent a promising class of antitumor compounds. Among AR, E-3-(3'-adamantan-1-yl-4'-hydroxybiphenyl-4-yl)acrylic acid (adarotene), has been extensively investigated. In the present work we report the results of our efforts to develop new adarotene-related atypical retinoids endowed also with POLA1 inhibitory activity. The effects of the synthesized compounds on cell growth were determined on a panel of human and hematological cancer cell lines. The most promising compounds showed antitumor activity against several tumor histotypes and increased cytotoxic activity against an adarotene-resistant cell line, compared to the parent molecule. The antitumor activity of a selected compound was evaluated on HT-29 human colon carcinoma and human mesothelioma (MM487) xenografts. Particularly significant was the in vivo activity of the compound as a single agent compared to adarotene and cisplatin, against pleural mesothelioma MM487. No reduction of mice body weight was observed, thus suggesting a higher tolerability with respect to the parent compound adarotene.

Background: There is an increasing interest in TRPM8 ligands of medicinal interest, the rational design of which can be nowadays supported by structure-based in silico studies based on the recently resolved TRPM8 structures. Methods: The study involves the generation of a reliable hTRPM8 homology model, the reliability of which was assessed by a 1.0 &mu, s MD simulation which was also used to generate multiple receptor conformations for the following structure-based virtual screening (VS) campaigns, docking simulations utilized different programs and involved all monomers of the selected frames, the so computed docking scores were combined by consensus approaches based on the EFO algorithm. Results: The obtained models revealed very satisfactory performances, LiGen&trade, provided the best results among the tested docking programs, the combination of docking results from the four monomers elicited a markedly beneficial effect on the computed consensus models. Conclusions: The generated hTRPM8 model appears to be amenable for successful structure-based VS studies, cross-talk modulating effects between interacting monomers on the binding sites can be accounted for by combining docking simulations as performed on all the monomers, this strategy can have general applicability for docking simulations involving quaternary protein structures with multiple identical binding pockets.

Protein carbonylation represents one of the most important oxidative-based modifications involving nucleophilic amino acids and affecting protein folding and function. Protein carbonylation is induced by electrophilic carbonyl species and is an highly selective process since few nucleophilic residues are carbonylated within each protein. While considering the great interest for protein carbonylation, few studies investigated the factors which render a nucleophilic residue susceptible to carbonylation. Hence, the present study is aimed to delve into the factors which modulate the reactivity of cysteine, histidine and lysine residues towards α,β unsaturated carbonyls by a retrospective analysis of the available studies which identified the adducted residues for proteins, the structure of which was resolved. Such an analysis involved different parameters including exposure, nucleophilicity, surrounding residues and capacity to attract carbonyl species (as derived by docking simulations). The obtained results allowed a meaningful clustering of the analyzed proteins suggesting that on average carbonylation selectivity increases with protein size. The comparison between adducted and unreactive residues revealed differences in all monitored parameters which are markedly more pronounced for cysteines compared to lysines and histidines. Overall, these results suggest that cysteine's carbonylation is a finely (and reasonably purposely) modulated process, while the carbonylation of lysines and histidines seems to be a fairly random event in which limited differences influence their reactivity.

The study proposes a novel consensus strategy based on linear combinations of different docking scores to be used in the evaluation of virtual screening campaigns. The consensus models are generated by applying the recently proposed Enrichment Factor Optimization (EFO) method, which develops the linear equations by exhaustively combining the available docking scores and by optimizing the resulting enrichment factors. The performances of such a consensus strategy were evaluated by simulating the entire Directory of Useful Decoys (DUD datasets). In detail, the poses were initially generated by the PLANTS docking program and then rescored by ReScore+ with and without the minimization of the complexes. The so calculated scores were then used to generate the mentioned consensus models including two or three different scoring functions. The reliability of the generated models was assessed by a per target validation as performed by default by the EFO approach. The encouraging performances of the here proposed consensus strategy are emphasized by the average increase of the 17% in the Top 1% enrichment factor (EF) values when comparing the single best score with the linear combination of three scores. Specifically, kinases offer a truly convincing demonstration of the efficacy of the here proposed consensus strategy since their Top 1% EF average ranges from 6.4 when using the single best performing primary score to 23.5 when linearly combining scoring functions. The beneficial effects of this consensus approach are clearly noticeable even when considering the entire DUD datasets as evidenced by the area under the curve (AUC) averages revealing a 14% increase when combining three scores. The reached AUC values compare very well with those reported in literature by an extended set of recent benchmarking studies and the three-variable models afford the highest AUC average.

[Image: see text] Alzheimer’s disease (AD) represents a global problem, with an estimation of the majority of dementia patients in low- and middle-income countries by 2050. Thus, the development of sustainable drugs has attracted much attention in recent years. In light of this, taking inspiration from the HDAC inhibitor vorinostat (1), we develop the first HDAC inhibitors derived from cashew nut shell liquid (CNSL), an inexpensive agro-food waste material. CNSL derivatives 8 and 9 display a HDAC inhibitory profile similar to 1, together with a more promising safety for 9 compared to 1. Moreover, both compounds and particularly 9 were able to effectively modulate glial cell-induced inflammation and to revert the pro-inflammatory phenotype. All these results demonstrate that the use of inexpensive food waste materials could be successfully applied for the development of accessible and sustainable drug candidates for the treatment of AD.

Alzheimer's disease is likely to be caused by copathogenic factors including aggregation of Aβ peptides into oligomers and fibrils, neuroinflammation, and oxidative stress. To date, no effective treatments are available, and because of the multifactorial nature of the disease, it emerges the need to act on different and simultaneous fronts. Despite the multiple biological activities ascribed to curcumin as neuroprotector, its poor bioavailability and toxicity limit the success in clinical outcomes. To tackle Alzheimer's disease on these aspects, the curcumin template was suitably modified and a small set of analogues was attained. In particular, derivative 1 turned out to be less toxic than curcumin. As evidenced by capillary electrophoresis and transmission electron microscopy studies, 1 proved to inhibit the formation of large toxic Aβ oligomers, by shifting the equilibrium toward smaller nontoxic assemblies and to limit the formation of insoluble fibrils. These findings were supported by molecular docking and steered molecular dynamics simulations which confirmed the superior capacity of 1 to bind Aβ structures of different complexity. Remarkably, 1 also showed in vitro anti-inflammatory and antioxidant properties. In summary, the curcumin-based analogue 1 emerged as multipotent compound worthy to be further investigated and exploited in the Alzheimer's disease multitarget context.

Carnosine (β-alanyl-L-histidine) is a natural peptide that have been described as a potential pharmacological agent owing to some positive outcomes from several pharmacological tests in animal models of human diseases. However, carnosine has limited activity in humans since the peptide upon absorption is rapidly hydrolyzed in the serum by the enzyme carnosinase (i.e. CN1; E.C. 3.4.13.20). Over the years the main approaches aimed at limiting carnosine hydrolysis have been focused on obtaining CN1-stable derivatives with an increased bioavailability and unmodified or enhanced activity. Only recently the hypothesis of co-administration of carnosine and selective inhibitors of CN1 have been proposed. Such an approach requires reliable methods for screening the effect on carnosine hydrolysis rate operated by CN1 in a throughput scale allowing to test from few compounds up to whole compound libraries. The only assay with such features available in literature relies on ortho-phtalaldehyde (OPA) derivatization of the hydrolysis product (i.e. histidine), followed by a fluorimetric read. Herein, we propose an alternative method based on a direct measurement of the residual substrate by liquid chromatography-mass spectrometry (LC–MS). The assay demonstrated to be reliable since gave results comparable to literature data concerning the hydrolysis rate of carnosine as determined into human serum. Moreover, the method was quite flexible and easily adaptable to a substrate change, as demonstrated by the measurement of the hydrolysis rate of all the natural analogs of carnosine. In this context the data collected for anserine suggest that our method looked more reliable and substrate change can undergo an underestimation of hydrolytic activity in OPA -based assays.