Search

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.

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.

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.

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.

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.

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.