Search

Fragment-Based Drug Discovery (FBDD) approaches have gained popularity not only in industry but also in academic research institutes. However, the computational prediction of the binding mode adopted by fragment-like molecules within a protein binding site is still a very challenging task. One of the most crucial aspects of fragment binding is related to the large amounts of bound waters in the targeted binding pocket. The binding affinity of fragments may not be sufficient to displace the bound water molecules. In the present work, we confirmed the importance of the bound water molecules in the correct prediction of the fragment binding mode. Moreover, we investigate whether the use of methods based on explicit solvent molecular dynamics simulations can improve the accuracy of fragment posing. The protein chosen for this study is HSP-90.

Fragment-based drug discovery (FBDD) is an innovative approach, progressively more applied in the academic and industrial context, to enhance hit identification for previously considered undruggable biological targets. In particular, FBDD discovers low-molecular-weight (LMW) ligands (<300Da) able to bind to therapeutically relevant macromolecules in an affinity range from the micromolar (μM) to millimolar (mM). X-ray crystallography (XRC) and nuclear magnetic resonance (NMR) spectroscopy are commonly the methods of choice to obtain 3D information about the bound ligand-protein complex, but this can occasionally be problematic, mainly for early, low-affinity fragments. The recent development of computational fragment-based approaches provides a further strategy for improving the identification of fragment hits. In this review, we summarize the state of the art of molecular dynamics simulations approaches used in FBDD, and discuss limitations and future perspectives for these approaches., (Copyright © 2020 The Author(s). Published by Elsevier Ltd.. All rights reserved.)

While a plethora of different protein-ligand docking protocols have been developed over the past twenty years, their performances greatly depend on the provided input protein-ligand pair. In this study, we developed a machine-learning model that uses a combination of convolutional and fully connected neural networks for the task of predicting the performance of several popular docking protocols given a protein structure and a small compound. We also rigorously evaluated the performance of our model using a widely available database of protein-ligand complexes and different types of data splits. We further open-source all code related to this study so that potential users can make informed selections on which protocol is best suited for their particular protein-ligand pair.

Adenosine receptors (ARs), like many otherGprotein-coupledreceptors (GPCRs), are targets of primary interest indrug design. However, one of the main limits for the development of drugs for this class of GPCRs is the complex selectivity profile usually displayed by ligands. Numerous efforts have been madefor clarifying the selectivity of ARs, leading to the development of many ligand-based models. The structure of the AR subtype A 1 (A 1 AR) has been recently solved,providing important structural insights. In the present work, we rationalized the selectivity profile of two selective A 1 AR and A 2A AR antagonists, investigating their recognition trajectories obtained by Supervised Molecular Dynamics from an unbound state and monitoring the role of the water molecules in the binding site., Competing Interests: The authors declare no conflict of interest.

Recent studies have highlighted the key role of Casein kinase 1 δ (CK1δ) in the development of several neurodegenerative pathologies, such as Alzheimer's disease (AD), Parkinson's disease (PD), and amyotrophic lateral sclerosis (ALS). So far, CK1δ inhibitors are noncovalent ATP competitive ligands and no drugs are currently available for this molecular target, hence the interest in developing new CK1δ inhibitors. The study aims to identify new inhibitors able to bind the enzyme; by a dual approach in silico/in vitro, the virtual screening has been performed on an in-house chemical library, which was previously designed and synthesized for other targets. The work can, therefore, be seen in the scaffold repurposing logic. The proposed strategy has led to the identification of two hits, having a novel scaffold in the landscape of CK1δ inhibitors and with an activity in the micromolar range., Competing Interests: The authors declare no competing financial interest., (Copyright © 2020 American Chemical Society.)

Although proteins have represented the molecular target of choice in the development of new drug candidates, the pharmaceutical importance of ribonucleic acids has gradually been growing. The increasing availability of structural information has brought to light the existence of peculiar three-dimensional RNA arrangements, which can, contrary to initial expectations, be recognized and selectively modulated through small chemical entities or peptides. The application of classical computational methodologies, such as molecular docking, for the rational development of RNA-binding candidates is, however, complicated by the peculiarities characterizing these macromolecules, such as the marked conformational flexibility, the singular charges distribution, and the relevant role of solvent molecules. In this work, we have thus validated and extended the applicability domain of SuMD, an all-atoms molecular dynamics protocol that allows to accelerate the sampling of molecular recognition events on a nanosecond timescale, to ribonucleotide targets of pharmaceutical interest. In particular, we have proven the methodological ability by reproducing the binding mode of viral or prokaryotic ribonucleic complexes, as well as that of artificially engineered aptamers, with an impressive degree of accuracy., (Copyright © 2020 Bissaro, Sturlese and Moro.)

Mitochondrial Ca 2+ uptake depends on the mitochondrial calcium uniporter (MCU) complex, a highly selective channel of the inner mitochondrial membrane (IMM). Here, we screen a library of 44,000 non-proprietary compounds for their ability to modulate mitochondrial Ca 2+ uptake. Two of them, named MCU-i4 and MCU-i11, are confirmed to reliably decrease mitochondrial Ca 2+ influx. Docking simulations reveal that these molecules directly bind a specific cleft in MICU1, a key element of the MCU complex that controls channel gating. Accordingly, in MICU1-silenced or deleted cells, the inhibitory effect of the two compounds is lost. Moreover, MCU-i4 and MCU-i11 fail to inhibit mitochondrial Ca 2+ uptake in cells expressing a MICU1 mutated in the critical amino acids that forge the predicted binding cleft. Finally, these compounds are tested ex vivo, revealing a primary role for mitochondrial Ca 2+ uptake in muscle growth. Overall, MCU-i4 and MCU-i11 represent leading molecules for the development of MICU1-targeting drugs., Competing Interests: Declaration of Interests B.J., C.B., V.T.-E., D.H., B.O., S.S., G.M., and M.F. are employees of Novartis Pharma AG, and some are also shareholders of Novartis., (Copyright © 2020 The Author(s). Published by Elsevier Inc. All rights reserved.)

A small number of fluorinated 7-phenyl-pyrroloquinolinone (7-PPyQ) derivatives was synthesized in an attempt to improve the metabolic stability of 3N-ethyl-7-PPyQ and 3N-benzoyl-7-PPyQ. The possible impacts of the fluorine-hydrogen isosterism on both biological activity and metabolic stability were evaluated. Introduction of a fluorine atom in the 2 or 3 position of the 7-phenyl ring yielded the 7-PPyQ derivatives 12, 13 and 15, which showed potent cytotoxicity (low micromolar and sub-nanomolar GI 50 s) both in human leukemic and solid tumor cell lines. None of them induced significant cell death in quiescent and proliferating human lymphocytes. Moreover, 12, 13 and 15 exhibited remarkable cytotoxic activity in the multidrug-resistant cell line CEM Vbl100 , suggesting that they are not substrates for P-glycoprotein. All compounds inhibited tubulin assembly and the binding of [ 3 H]colchicine to tubulin, with the best activity occurring with compound 15. Mechanistic studies carried out on compound 12 indicated that it caused (a) a strong G2/M arrest; (b) apoptosis in a time- and concentration-dependent manner; (c) a significant production of ROS (in good agreement with the observed mitochondrial depolarization); (d) caspase-3 and poly (ADP-ribose) polymerase activation; and (e) a decrease in the expression of anti-apoptotic proteins. In vivo experiments in a murine syngeneic tumor model demonstrated that compounds 12 and 15 significantly reduced tumor mass at doses four times lower than that required for the reference compound combretastatin A-4 phosphate. Neither monofluorination of the 7-phenyl ring of 3N-ethyl-7-PPyQ nor replacement of the benzoyl function of 3N-benzoyl-7-PPyQ with a 2-fluorobenzoyl moiety led to any improvement in the metabolic stability., (Copyright © 2019 Elsevier Masson SAS. All rights reserved.)

The use and misuse of antibiotics has resulted in critical conditions for drug-resistant bacteria emergency, accelerating the development of antimicrobial resistance (AMR). In this context, the co-administration of an antibiotic with a compound able to restore sufficient antibacterial activity may be a successful strategy. In particular, the identification of efflux pump inhibitors (EPIs) holds promise for new antibiotic resistance breakers (ARBs). Indeed, bacterial efflux pumps have a key role in AMR development; for instance, NorA efflux pump contributes to Staphylococcus aureus ( S. aureus ) resistance against fluoroquinolone antibiotics (e.g., ciprofloxacin) by promoting their active extrusion from the cells. Even though NorA efflux pump is known to be a potential target for EPIs development, the absence of structural information about this protein and the little knowledge available on its mechanism of action have strongly hampered rational drug discovery efforts in this area. In the present work, we investigated at the molecular level the substrate recognition pathway of NorA through a Supervised Molecular Dynamics (SuMD) approach, using a NorA homology model. Specific amino acids were identified as playing a key role in the efflux pump-mediated extrusion of its substrate, paving the way for a deeper understanding of both the mechanisms of action and the inhibition of such efflux pumps.

One of the most intriguing findings highlighted from G protein-coupled receptor (GPCR) crystallography is the presence, in many members of class A, of a partially hydrated sodium ion in the middle of the seven transmembrane helices (7TM) bundle. In particular, the human adenosine A 2A receptor (A 2A AR) is the first GPCR in which a monovalent sodium ion was crystallized in a distal site from the canonical orthosteric one, corroborating, from a structural point of view, its role as a negative allosteric modulator. However, the molecular mechanism by which the sodium ion influences the recognition of the A 2A AR agonists is not yet fully understood. In this study, the supervised molecular dynamics (SuMD) technique was exploited to analyse the sodium ion recognition mechanism and how its presence influences the binding of the endogenous agonist adenosine. Due to a higher degree of flexibility of the receptor extracellular (EC) vestibule, we propose the sodium-bound A 2A AR as less efficient in stabilizing the adenosine during the different steps of binding.

The number of entries in the Protein Data Bank (PDB) has doubled in the last decade, and it has increased tenfold in the last twenty years. The availability of an ever-growing number of structures is having a huge impact on the Structure-Based Drug Discovery (SBDD), allowing investigation of new targets and giving the possibility to have multiple structures of the same macromolecule in a complex with different ligands. Such a large resource often implies the choice of the most suitable complex for molecular docking calculation, and this task is complicated by the plethora of possible posing and scoring function algorithms available, which may influence the quality of the outcomes. Here, we report a large benchmark performed on the PDBbind database containing more than four thousand entries and seventeen popular docking protocols. We found that, even in protein families wherein docking protocols generally showed acceptable results, certain ligand-protein complexes are poorly reproduced in the self-docking procedure. Such a trend in certain protein families is more pronounced, and this underlines the importance in identification of a suitable protein-ligand conformation coupled to a well-performing docking protocol.

Exploring at the molecular level, all possible ligand-protein approaching pathways and, consequently, identifying the energetically favorable binding sites is considered crucial to depict a clear picture of the whole scenario of ligand-protein binding. In fact, a ligand can recognize a protein in multiple binding sites, adopting multiple conformations in every single binding site and inducing protein modifications upon binding. In the present work, we would like to present how it is possible to couple a supervised molecular dynamics (SuMD) approach to explore, from an unbound state, the most energetically favorable recognition pathways of the ligand to its protein, with an enthalpic and entropic characterization of the most stable ligand-protein bound states, using the protein kinase CK2α as a prototype study. We identified two accessory binding pockets surrounding the ATP-binding site having a strong enthalpic contribution but a different configurational entropy contribution, suggesting that they play a different role., Competing Interests: The authors declare no competing financial interest.

Riluzole, approved by the US Food and Drug Administration (FDA) in 1995, is the most widespread oral treatment for the fatal neurodegenerative disorder amyotrophic lateral sclerosis (ALS). The drug, whose mechanism of action is still obscure, mitigates progression of the illness, but unfortunately with only limited improvements. Herein we report the first demonstration, using a combination of computational and in vitro studies, that riluzole is an ATP-competitive inhibitor of the protein kinase CK1 isoform δ, with an IC 50 value of 16.1 μm. This allows us to rewrite its possible molecular mechanism of action in the treatment of ALS. The inhibition of CK1δ catalytic activity indeed links the two main pathological hallmarks of ALS: transactive response DNA-binding protein of 43 kDa (TDP-43) proteinopathy and glutamate excitotoxicity, exacerbated by the loss of expression of glial excitatory amino acid transporter-2 (EAAT2)., (© 2018 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim.)

Glutaredoxins (Grx) are small proteins conserved throughout all the kingdoms of life that are engaged in a wide variety of biological processes and share a common thioredoxin-fold. Among them, class II Grx are redox-inactive proteins involved in iron-sulfur (FeS) metabolism. They contain a single thiol group in their active site and use low molecular mass thiols such as glutathione as ligand for binding FeS-clusters. In this study, we investigated molecular aspects of 1CGrx1 from the pathogenic parasite Trypanosoma brucei brucei, a mitochondrial class II Grx that fulfills an indispensable role in vivo. Mitochondrial 1CGrx1 from trypanosomes differs from orthologues in several features including the presence of a parasite-specific N-terminal extension (NTE) whose role has yet to be elucidated. Previously we have solved the structure of a truncated form of 1CGrx1 containing only the conserved glutaredoxin domain but lacking the NTE. Our aim here is to investigate the effect of the NTE on the conformation of the protein. We therefore solved the NMR structure of the full-length protein, which reveals subtle but significant differences with the structure of the NTE-less form. By means of different experimental approaches, the NTE proved to be intrinsically disordered and not involved in the non-redox dependent protein dimerization, as previously suggested. Interestingly, the portion comprising residues 65-76 of the NTE modulates the conformational dynamics of the glutathione-binding pocket, which may play a role in iron-sulfur cluster assembly and delivery. Furthermore, we disclosed that the class II-strictly conserved loop that precedes the active site is critical for stabilizing the protein structure. So far, this represents the first communication of a Grx containing an intrinsically disordered region that defines a new protein subgroup within class II Grx.

The complex formation between the proteins apoptosis-inducing factor (AIF) and cyclophilin A (CypA) following oxidative stress in neuronal cells has been suggested as a main target for reverting ischemia-stroke damage. Recently, a peptide encompassing AIF residues 370-394 has been developed to target the AIF-binding site on CypA, to prevent the association between the two proteins and suppress glutamate-induced cell death in neuronal cells. Using a combined approach based on NMR spectroscopy, synthesis and in vitro testing of all Ala-scan mutants of the peptide and molecular docking/molecular dynamics, we have generated a detailed model of the AIF (370-394)/CypA complex. The model suggests us that the central region of the peptide spanning residues V374-K384 mostly interacts with the protein and that for efficient complex inhibition and preservation of CypA activity, it is bent around amino acids F46-G75 of the protein. The model is consistent with experimental data also from previous works and supports the concept that the peptide does not interfere with other CypA activities unrelated to AIF activation; therefore, it may serve as an ideal template for generating future non-peptidic antagonists., (© 2018 The Author(s). Published by Portland Press Limited on behalf of the Biochemical Society.)

Unquestionably, water appears to be an active player in the noncovalent protein-ligand binding process, as it can either bridge interactions between protein and ligand or can be replaced by the bound ligand. Accordingly, in the last decade, alternative computational methodologies have been sought with the aim of predicting the position and thermodynamic profile of water molecules (i.e., hydration sites) in the binding site using either the ligand-bound or ligand-free protein conformation. Herein, we present an alternative approach, named AquaMMapS, that provides a three-dimensional sampling of putative hydration sites. Interestingly, AquaMMapS can post-inspect molecular dynamics (MD) trajectories obtained from different MD engines using indifferently crystallographic or docking-driven structures as a starting point. Moreover, AquaMMapS is naturally integrated into supervised molecular dynamics (SuMD) simulations, presenting the possibility to inspect hydration sites during the ligand-protein association process. Finally, a penalty scoring method, named AquaMMapScoring(AMS), was developed to evaluate the number and nature of the water molecules displaced by a ligand approaching its binding site during the binding event, guiding a medicinal chemist to explore the most suitable regions of a ligand that can be decorated either with or without interfering with the interaction networks mediated by water molecules with specific recognition regions of the protein., (© 2018 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim.)

Significance: Major pathogenic enterobacteria and protozoan parasites from the phylum Euglenozoa, such as trypanosomatids, are endowed with glutathione (GSH)-spermidine (Sp) derivatives that play important roles in signaling and metal and thiol-redox homeostasis. For some Euglenozoa lineages, the GSH-Sp conjugates represent the main redox cosubstrates around which entire new redox systems have evolved. Several proteins underwent molecular adaptations to synthesize and utilize the new polyamine-based thiols. Recent Advances: The genomes of closely related organisms have recently been sequenced, which allows mining and analysis of gene sequences that belong to these peculiar redox systems. Similarly, the three-dimensional structures of several of these proteins have been solved, which allows for comparison with their counterparts in classical redox systems that rely on GSH/glutaredoxin and thioredoxin., Critical Issues: The evolutionary and structural aspects related to the emergence and use of GSH-Sp conjugates in Euglenozoa are reviewed focusing on unique structural specializations that proteins developed to use N 1 ,N 8 -bisglutathionylspermidine (trypanothione) as redox cosubstrate. An updated overview on the biochemical and biological significance of the major enzymatic activities is also provided., Future Directions: A thiol-redox system strictly dependent on trypanothione is a feature unique to trypanosomatids. The physicochemical properties of the polyamine-GSH conjugates were a major driving force for structural adaptation of proteins that use these thiols as ligand and redox cofactor. In fact, the structural differences of indispensable components of this system can be exploited toward selective drug development. Future research should clarify whether additional cellular processes are regulated by the trypanothione system. Antioxid. Redox Signal. 28, 463-486.

Supervised MD (SuMD) is a computational method that enables the exploration of ligand-receptor recognition pathway in a reduced timescale. The performance speedup is due to the incorporation of a tabu-like supervision algorithm on the ligand-receptor approaching distance into a classic molecular dynamics (MD) simulation. SuMD enables the investigation of ligand-receptor binding events independently from the starting position, chemical structure of the ligand (small molecules or peptides), and also from its receptor-binding affinity. The application of SuMD highlights an appreciable capability of the technique to reproduce the crystallographic structures of several ligand-protein complexes and can provide high-quality protein-ligand models of for which yet experimental confirmation of binding mode is not available.

Molecular docking is a powerful tool in the field of computer-aided molecular design. In particular, it is the technique of choice for the prediction of a ligand pose within its target binding site. A multitude of docking methods is available nowadays, whose performance may vary depending on the data set. Therefore, some non-trivial choices should be made before starting a docking simulation. In the same framework, the selection of the target structure to use could be challenging, since the number of available experimental structures is increasing. Both issues have been explored within this work. The pose prediction of a pool of 36 compounds provided by D3R Grand Challenge 2 organizers was preceded by a pipeline to choose the best protein/docking-method couple for each blind ligand. An integrated benchmark approach including ligand shape comparison and cross-docking evaluations was implemented inside our DockBench software. The results are encouraging and show that bringing attention to the choice of the docking simulation fundamental components improves the results of the binding mode predictions.

Earlier studies had confirmed that the 7-phenylpyrroloquinolinone (7-PPyQ) nucleus was an important scaffold for new chemotherapeutic drugs targeting microtubules. For wide-ranging SARs, a series of derivatives were synthesized through a robust procedure. For comparison with the reference 3-ethyl-7-PPyQ 31, the angular geometry and substituents at the 3 and 7 positions were varied to explore interactions inside the colchicine site of tubulin. Of the new compounds synthesized, potent cytotoxicity (low and sub-nanomolar GI 50 values) was observed with 21 and 24, both more potent than 31, in both leukemic and solid tumor cell lines. Neither compound 21 nor 24 induced significant cell death in normal human lymphocytes, suggesting that the compounds may be selectively active against cancer cells. In particular, 24 was a potent inducer of apoptosis in the A549 and HeLa cell lines. With both compounds, induction of apoptosis was associated with dissipation of the mitochondrial transmembrane potential and production of reactive oxygen species, indicating that cells treated with the compounds followed the intrinsic pathway of apoptosis. Moreover, immunoblot analysis revealed that compound 24 even at 50 nM reduced the expression of anti-apoptotic proteins such as Bcl-2 and Mcl-1. Finally, molecular docking studies of the newly synthesized compounds demonstrate that active pyrroloquinolinone derivatives strongly bind in the colchicine site of β-tubulin., (Copyright © 2017 Elsevier Masson SAS. All rights reserved.)

New 7-amino-2-phenylpyrazolo[4,3-d]pyrimidine derivatives, substituted at the 5-position with aryl(alkyl)amino- and 4-substituted-piperazin-1-yl- moieties, were synthesized with the aim of targeting human (h) adenosine A 1 and/or A 2A receptor subtypes. On the whole, the novel derivatives 1-24 shared scarce or no affinities for the off-target hA 2B and hA 3 ARs. The 5-(4-hydroxyphenethylamino)- derivative 12 showed both good affinity (K i = 150 nM) and the best selectivity for the hA 2A AR while the 5-benzylamino-substituted 5 displayed the best combined hA 2A (K i = 123 nM) and A 1 AR affinity (K i = 25 nM). The 5-phenethylamino moiety (compound 6) achieved nanomolar affinity (K i = 11 nM) and good selectivity for the hA 1 AR. The 5-(N 4 -substituted-piperazin-1-yl) derivatives 15-24 bind the hA 1 AR subtype with affinities falling in the high nanomolar range. A structure-based molecular modeling study was conducted to rationalize the experimental binding data from a molecular point of view using both molecular docking studies and Interaction Energy Fingerprints (IEFs) analysis.[Formula: see text].

The serine-threonine checkpoint kinase 1 (Chk1) plays a critical role in the cell cycle arrest in response to DNA damage. In the last decade, Chk1 inhibitors have emerged as a novel therapeutic strategy to potentiate the anti-tumour efficacy of cytotoxic chemotherapeutic agents. In the search for new Chk1 inhibitors, a congeneric series of 2-aryl-2 H-pyrazolo[4,3-c]quinolin-3-one (PQ) was evaluated by in-vitro and in-silico approaches for the first time. A total of 30 PQ structures were synthesised in good to excellent yields using conventional or microwave heating, highlighting that 14 of them are new chemical entities. Noteworthy, in this preliminary study two compounds 4e 2 and 4h 2 have shown a modest but significant reduction in the basal activity of the Chk1 kinase. Starting from these preliminary results, we have designed the second generation of analogous in this class and further studies are in progress in our laboratories.

A novel series of variously substituted N-[3-(9H-carbazol-9-yl)-2-hydroxypropyl]-arylsulfonamides has been synthesized and assayed for β-Secretase (BACE1) inhibitory activity. BACE1 is a widely recognized drug target for the prevention and treatment of Alzheimer's Disease (AD). The introduction of benzyl substituents on the nitrogen atom of the arylsulfonamide moiety has so far led to the best results, with three derivatives showing IC 50 values ranging from 1.6 to 1.9 μM. Therefore, a significant improvement over the previously reported series of N-carboxamides (displaying IC 50 's ≥ 2.5 μM) has been achieved, thus suggesting an active role of the sulfonamido-portion in the inhibition process. Preliminary molecular modeling studies have been carried out to rationalize the observed structure-activity relationships., (Copyright © 2017 Elsevier Ltd. All rights reserved.)

Many cysteine-stabilized antimicrobial peptides from a variety of living organisms could be good candidates for the development of anti-infective agents. In the absence of experimentally obtained structural data, peptide modeling is an essential tool for understanding structure-activity relationships and for optimizing the bioactive moieties. Focusing on cysteine-rich peptide structures, we reproduced the case of structure predictions in the so-called midnight zone. We developed our protocol on a training set derived by clustering the available cysteine-stabilized αβ (CSαβ) structures in nine different representative families and tested it on peptides randomly selected from each family. Starting from draft models, we tested a structure-based disulfide predictor and we used cysteine distances as constraints during molecular dynamics. Finally, we proposed an analysis for final structure selection. Accordingly, we obtained a mean root mean square deviation improvement of 21% for the test set. Our findings demonstrate that it is possible to predict the network of disulfide bridges in cysteine-stabilized peptides and to use this result to improve the accuracy of structural predictions. Finally, we applied the methods to predict the structure of royalisin, a cysteine-rich peptide with unknown structure.

Peptides have gained increased interest as therapeutic agents during recent years. The high specificity and relatively low toxicity of peptide drugs derive from their extremely tight binding to their targets. Indeed, understanding the molecular mechanism of protein-peptide recognition has important implications in the fields of biology, medicine, and pharmaceutical sciences. Even if crystallography and nuclear magnetic resonance are offering valuable atomic insights into the assembling of the protein-peptide complexes, the mechanism of their recognition and binding events remains largely unclear. In this work we report, for the first time, the use of a supervised molecular dynamics approach to explore the possible protein-peptide binding pathways within a timescale reduced up to three orders of magnitude compared with classical molecular dynamics. The better and faster understating of the protein-peptide recognition pathways could be very beneficial in enlarging the applicability of peptide-based drug design approaches in several biotechnological and pharmaceutical fields., (Copyright © 2017 Elsevier Ltd. All rights reserved.)

A small library of 7-pyrrolo[3,2-f]quinolinones was obtained by introducing benzoyl, sulfonyl and carbamoyl side chains at the 3-N position, and their cytotoxicity against a panel of leukemic and solid tumor cell lines was evaluated. Most of them showed high antiproliferative activity with GI 50 s ranging from micro-to sub-nanomolar values, and these values correlated well with the inhibitory activities of the compounds against tubulin polymerization. Based on a recently proposed colchicine bind site inhibitors (CBSIs) pharmacophore, the interactions of the novel 7-PPyQs at the colchicine domain were rationalized. The most active compounds (4a and 4b) did not induce significant cell death in normal human lymphocytes, suggesting that the compounds may be selective against cancer cells. In particular, 4a was a potent inducer of apoptosis in both the HeLa and Jurkat cell lines. On the other hand, the sulfonyl derivative 4b exhibited a lower potency in comparison with 4a. With both compounds, induction of apoptosis was associated with dissipation of the mitochondrial transmembrane potential and production of reactive oxygen species, suggesting that cells treated with the compounds followed the intrinsic pathway of apoptosis., (Copyright © 2016 Elsevier Masson SAS. All rights reserved.)

In this review, we present a survey of the recent advances carried out by our research groups in the field of ligand-GPCRs recognition process simulations recently implemented at the Molecular Modeling Section (MMS) of the University of Padova. We briefly describe a platform of tools we have tuned to aid the identification of novel GPCRs binders and the better understanding of their binding mechanisms, based on two extensively used computational techniques such as molecular docking and MD simulations. The developed methodologies encompass: (i) the selection of suitable protocols for docking studies, (ii) the exploration of the dynamical evolution of ligand-protein interaction networks, (iii) the detailed investigation of the role of water molecules upon ligand binding, and (iv) a glance at the way the ligand might go through prior reaching the binding site., (© 2016 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim.)

Structure-based drug design (SBDD) has matured within the last two decades as a valuable tool for the optimization of low molecular weight lead compounds to highly potent drugs. The key step in SBDD requires knowledge of the three-dimensional structure of the target-ligand complex, which is usually determined by X-ray crystallography. In the absence of structural information for the complex, SBDD relies on the generation of plausible molecular docking models. However, molecular docking protocols suffer from inaccuracies in the description of the interaction energies between the ligand and the target molecule, and often fail in the prediction of the correct binding mode. In this context, the appropriate selection of the most accurate docking protocol is absolutely relevant for the final molecular docking result, even if addressing this point is absolutely not a trivial task. D3R Grand Challenge 2015 has represented a precious opportunity to test the performance of DockBench, an integrate informatics platform to automatically compare RMDS-based molecular docking performances of different docking/scoring methods. The overall performance resulted in the blind prediction are encouraging in particular for the pose prediction task, in which several complex were predicted with a sufficient accuracy for medicinal chemistry purposes.

In this work, a sensitive and convenient protease-based fluorimetric high-throughput screening (HTS) assay for determining peptidyl-prolyl cis-trans isomerase activity was developed. The assay was based on a new intramolecularly quenched substrate, whose fluorescence and structural properties were examined together with kinetic constants and the effects of solvents on its isomerization process. Pilot screens performed using the Library of Pharmacologically Active Compounds (LOPAC) and cyclophilin A (CypA), as isomerase model enzyme, indicated that the assay was robust for HTS, and that comparable results were obtained with a CypA inhibitor tested both manually and automatically. Moreover, a new compound that inhibits CypA activity with an IC50 in the low micromolar range was identified. Molecular docking studies revealed that the molecule shows a notable shape complementarity with the catalytic pocket confirming the experimental observations. Due to its simplicity and precision in the determination of extent of inhibition and reaction rates required for kinetic analysis, this assay offers many advantages over other commonly used assays., (© 2016 Society for Laboratory Automation and Screening.)

Molecular recognition is a crucial issue when aiming to interpret the mechanism of known active substances as well as to develop novel active candidates. Unfortunately, simulating the binding process is still a challenging task because it requires classical MD experiments in a long microsecond time scale that are affordable only with a high-level computational capacity. In order to overcome this limiting factor, we have recently implemented an alternative MD approach, named supervised molecular dynamics (SuMD), and successfully applied it to G protein-coupled receptors (GPCRs). SuMD enables the investigation of ligand-receptor binding events independently from the starting position, chemical structure of the ligand, and also from its receptor binding affinity. In this article, we present an extension of the SuMD application domain including different types of proteins in comparison with GPCRs. In particular, we have deeply analyzed the ligand-protein recognition pathways of six different case studies that we grouped into two different classes: globular and membrane proteins. Moreover, we introduce the SuMD-Analyzer tool that we have specifically implemented to help the user in the analysis of the SuMD trajectories. Finally, we emphasize the limit of the SuMD applicability domain as well as its strengths in analyzing the complexity of ligand-protein recognition pathways.

Trypanosomatids are parasites responsible for several tropical and subtropical diseases, such as Chaga's disease, sleeping sickness and Leishmaniasis. In contrast to the mammalian host, the thiol-redox metabolism of these pathogens depends on trypanothione [bis-glutathionylspermidine, T(SH)2] instead of glutathione (GSH) providing a set of lineage-specific proteins as drug target candidates. Glutaredoxins (Grx) are ubiquitous small thiol-disulfide oxidoreductases that belong to the thioredoxin-fold family. They play a central role in redox homeostasis and iron sulfur-cluster biogenesis. Each species, including trypanosomes, possesses its own set of isoforms distributed in different subcellular compartments. The genome of trypanosomatids encodes for two class I (dithiolic) Grxs named 2-C-Grx1 and 2-C-Grx2. Both proteins were shown to efficiently reduce different disulfides at the expenses of T(SH)2 using a mechanism that involves the two cysteines in the active site. Moreover, the cytosolic Trypanosoma brucei 2-C-Grx1 but not the mitochondrial 2-C-Grx2 was able to coordinate an iron-sulfur cluster with T(SH)2 or GSH as ligand. As a first step to unravel the structural basis for the specificity observed in the trypanosomal glutaredoxins, we present here the NMR resonance assignment of 2-C-Grx1 from the parasite T. brucei brucei.

In silico three-dimensional (3D) molecular modeling tools based upon the receptor/enzyme-ligand docking simulation in protein crystal structures and/or homology modeling of receptors have been reliably used in pharmacological research and development for decades. Molecular docking methodologies are helpful for revealing facets of activation and inactivation, thus improving mechanistic understanding and predicting molecular ligand binding activity, and they can have a high level of accuracy, and have also been explored and applied in chemical risk assessment. This computational approach is, however, only applicable for chemical hazard identification situations where the specific target receptor for a given chemical is known and the crystal structure/homology model of the receptor is available.

A new series of pyrazolo[4,3-e][1,2,4]triazolo[1,5-c]pyrimidine (PTP) derivatives has been developed in order to explore their affinity and selectivity profile at the four adenosine receptor subtypes. In particular, the PTP scaffold was conjugated at the C2 position with the 1-(3-trifluoromethyl-benzyl)-1H-pyrazole, a group believed to confer potency and selectivity toward the human (h) A2B adenosine receptor (AR) to the xanthine ligand 8-(1-(3-(trifluoromethyl)benzyl)-1H-pyrazol-4-yl)-1,3-dimethyl-1H-purine-2,6(3H,7H)-dione (CVT 6975). Interestingly, the synthesized compounds turned out to be inactive at the hA2B AR but they displayed affinity at the hA3 AR in the nanomolar range. The best compound of the series (6) shows both high affinity (hA3 AR Ki = 11 nM) and selectivity (A1/A3 and A2A/A3 > 9090; A2B/A3 > 909) at the hA3 AR. To better rationalize these results, a molecular docking study on the four AR subtypes was performed for all the synthesized compounds. In addition, CTV 6975 and two close analogues have been subjected to the same molecular docking protocol to investigate the role of the 1-(3-trifluoromethyl-benzyl)-1H-pyrazole on the binding at the four ARs.

Nuclear magnetic resonance (NMR) spectroscopy and molecular docking are regularly being employed as helpful tools of drug discovery research. Molecular docking is an extremely rapid method to evaluate possible binders from a large chemical library in a fast and cheap manner. NMR techniques can directly detect a protein-ligand interaction, can determine the corresponding association constant, and can consistently identify the ligand binding cavity. Consequently, molecular docking and NMR techniques are naturally complementary techniques where the combination of the two has the potential to improve the overall efficiency of drug discovery process. In this review, we would like to summarize the state of the art of docking methods which have been recently bridged to NMR experiments to identify novel and effective therapeutic drug candidates., (© 2015 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim.)

Virtual screening (VS) is a computational methodology that streamlines the drug discovery process by reducing costs and required resources through the in silico identification of potential drug candidates. Structure-based VS (SBVS) exploits knowledge about the three-dimensional (3D) structure of protein targets and uses the docking methodology as search engine for novel hits. The success of a SBVS campaign strongly depends upon the accuracy of the docking protocol used to select the candidates from large chemical libraries. The identification of suitable protocols is therefore a crucial step in the setup of SBVS experiments. Carrying out extensive benchmark studies, however, is usually a tangled task that requires users' proficiency in handling different file formats and philosophies at the basis of the plethora of existing software packages. We present here DockBench 1.0, a platform available free of charge that eases the pipeline by automating the entire procedure, from docking benchmark to VS setups. In its current implementation, DockBench 1.0 handles seven docking software packages and offers the possibility to test up to seventeen different protocols. The main features of our platform are presented here and the results of the benchmark study of human Checkpoint kinase 1 (hChk1) are discussed as validation test.

Glutaredoxins (Grx) are small proteins, conserved throughout all the kingdoms of life, which are engaged in a wide variety of biological processes. According to the number of cysteines in their active site, Grx are classified as dithiolic or monothiolic (1-C-Grx). In most organisms, 1-C-Grx are implicated in iron-sulfur cluster (FeS) metabolism and utilize glutathione as cofactor. Trypanosomatids are parasitic protozoa of the order Kinetoplastida, which cause severe diseases in humans and domestic animals. These parasites exploit a unique thiol-dependent redox system based on bis(glutathionyl)spermidine (trypanothione) rather than on glutathione. Mitochondrial 1-C-Grx1 from trypanosomes differs from orthologues in several features including the use of trypanothione as ligand for FeS binding and the presence of a parasite-specific N-terminal extension. We have recently shown that 1-C-Grx1 from Trypanosoma brucei is indispensable for parasite survival in mouse, making this protein a potential drug target candidate against trypanosomiasis. However, structural information for the full-length form of 1-C-Grx1 is still lacking. Here, we report the NMR resonance assignment of the mature form of Tb1-C-Grx1 including an N-terminal tail, paving the way to disclose the role of this intrinsically disordered region in the protein function.

Prestin (SLC26A5) is a member of the SLC26/SulP anion transporter family. Its unique quasi-piezoelectric mechanical activity generates fast cellular motility of cochlear outer hair cells, a key process underlying active amplification in the mammalian ear. Despite its established physiological role, it is essentially unknown how prestin can generate mechanical force, since structural information on SLC26/SulP proteins is lacking. Here we derive a structural model of prestin and related transporters by combining homology modelling, MD simulations and cysteine accessibility scanning. Prestin's transmembrane core region is organized in a 7+7 inverted repeat architecture. The model suggests a central cavity as the substrate-binding site located midway of the anion permeation pathway, which is supported by experimental solute accessibility and mutational analysis. Anion binding to this site also controls the electromotile activity of prestin. The combined structural and functional data provide a framework for understanding electromotility and anion transport by SLC26 transporters.

Aims: Monothiol glutaredoxins (1-C-Grxs) are small proteins linked to the cellular iron and redox metabolism. Trypanosoma brucei brucei, model organism for human African trypanosomiasis, expresses three 1-C-Grxs. 1-C-Grx1 is a highly abundant mitochondrial protein capable to bind an iron-sulfur cluster (ISC) in vitro using glutathione (GSH) as cofactor. We here report on the functional and structural analysis of 1-C-Grx1 in relation to its ISC-binding properties., Results: An N-terminal extension unique to 1-C-Grx1 from trypanosomatids affects the oligomeric structure and the ISC-binding capacity of the protein. The active-site Cys104 is essential for ISC binding, and the parasite-specific glutathionylspermidine and trypanothione can replace GSH as the ligands of the ISC. Interestingly, trypanothione forms stable protein-free ISC species that in vitro are incorporated into the dithiol T. brucei 2-C-Grx1, but not 1-C-Grx1. Overexpression of the C104S mutant of 1-C-Grx1 impairs disease progression in a mouse model. The structure of the Grx-domain of 1-C-Grx1 was solved by nuclear magnetic resonance spectroscopy. Despite the fact that several residues--which in other 1-C-Grxs are involved in the noncovalent binding of GSH--are conserved, different physicochemical approaches did not reveal any specific interaction between 1-C-Grx1 and free thiol ligands., Innovation: Parasite Grxs are able to coordinate an ISC formed with trypanothione, suggesting a new mechanism of ISC binding and a novel function for the parasite-specific dithiol. The first 3D structure and in vivo relevance of a 1-C-Grx from a pathogenic protozoan are reported., Conclusion: T. brucei 1-C-Grx1 is indispensable for mammalian parasitism and utilizes a new mechanism for ISC binding.

The physiological role of DJ-1, a protein involved in familial Parkinson disease is still controversial. One of the hypotheses proposed indicates a sensor role for oxidative stress, through oxidation of a conserved cysteine residue (Cys-106). The association of DJ-1 mutations with Parkinson disease suggests a loss of function, specific to dopaminergic neurons. Under oxidative conditions, highly reactive dopamine quinones (DAQs) can be produced, which can modify cysteine residues. In cellular models, DJ-1 was found covalently modified by dopamine. We analyzed the structural modifications induced on human DJ-1 by DAQs in vitro. We described the structural perturbations induced by DAQ adduct formation on each of the three cysteine residues of DJ-1 using specific mutants. Cys-53 is the most reactive residue and forms a covalent dimer also in SH-SY5Y DJ-1-transfected cells, but modification of Cys-106 induces the most severe structural perturbations; Cys-46 is not reactive. The relevance of these covalent modifications to the several functions ascribed to DJ-1 is discussed in the context of the cell response to a dopamine-derived oxidative insult.

Recent characterization of Mcl-1 as the primary anti-apoptotic Bcl-2 family member expressed in solid tumors, coupled with its ability to enable therapeutic resistance, has provided the impetus for further study into how Mcl-1 is involved in apoptosis signaling. Here, we employ Sabutoclax, a potent and effective Mcl-1 antagonist, as a competing agent to screen a randomized 12-residue phage display library for peptides that bind strongly to the Bcl-2 homology 3 (BH3) binding groove of Mcl-1. Although the screen identified a number of α-helical peptides with canonical BH3 domain sequences, it also isolated a pair of unique peptide sequences. These sequences exhibit a reverse organization of conserved hydrophobic and acidic residues when compared with canonical BH3 sequences, and we therefore refer to them as reverse BH3 (rBH3) peptides. Furthermore, studies of the rBH3 peptides using NMR spectroscopy, fluorescence polarization displacement assays, and alanine scanning data all suggest that they bind to the BH3 binding groove of Mcl-1 selectively over Bcl-x(L). A search for proteins containing the rBH3 motif has identified a number of interesting Mcl-1 protein partners, some of which have previously been associated with apoptosis regulation involving Mcl-1. These findings provide insights into the development of more specific Mcl-1 antagonists and open the way to the identification of a previously unknown family of apoptosis-regulating and Mcl-1 interacting proteins.

The N-terminal region of the parathyroid hormone (PTH) is sufficient to activate the G-protein-coupled PTH receptor 1 (PTHR1). The shortest PTH analogue displaying nanomolar potency is the undecapeptide H-Aib-Val-Aib-Glu-Ile-Gln-Leu-Nle-His-Gln-Har-NH(2) that contains two helix-stabilizing residues (Aib(1,3)). To increase the helical character and proteolytic stability of this linear peptide, we replaced Gln(6,10) with (a) Lys(6) and Glu(10) to introduce a lactam bridge and (b) Ser(6,10) to form a diester bridge upon cross-linking with adipic acid. These cyclopeptides were, respectively, 468-fold less and 12-fold more potent agonists than the linear analogue. Despite their different potencies, all three analogues adopted similar α-helix structures, as shown by NMR and molecular dynamics studies. However, the diester bridge could better mimic the orientation and chemical properties of the side chains of Gln(6) and Gln(10) in the linear PTH analogue than the lactam moiety. This is apparently important for efficient receptor activation and provides further insights into the receptor-bound ligand conformation.

The N-terminal 1-34 fragment of parathyroid hormone (PTH) is fully active in vitro and in vivo and it can reproduce all biological responses characteristic of the native intact PTH. Recently, analogues of PTH(1-11) fragments with helicity-enhancing substitutions have been demonstrated to yield potent analogues of PTH(1-34). The work describes the synthesis, biological activity and structure of analogues of the best modified PTH sequence H-Aib-Val-Aib-Glu-Ile-Gln-Leu-Nle-His-Gln-Har-NH2 (I). In particular, the effect of the Ala/Aib substitution at positions 1 and 3 as well as of the replacement of Nle in position 8 with D-Nle, L-(αMe)-Nle and D-(αMe)-Nle was studied. The resulting peptides were characterized structurally by CD spectroscopy, solution NMR and MD, and in vitro for activity with respect to the cognate receptor, parathyroid hormone receptor.

MMsINC (http://mms.dsfarm.unipd.it/MMsINC/search) is a database of non-redundant, richly annotated and biomedically relevant chemical structures. A primary goal of MMsINC is to guarantee the highest quality and the uniqueness of each entry. MMsINC then adds value to these entries by including the analysis of crucial chemical properties, such as ionization and tautomerization processes, and the in silico prediction of 24 important molecular properties in the biochemical profile of each structure. MMsINC is consequently a natural input for different chemoinformatics and virtual screening applications. In addition, MMsINC supports various types of queries, including substructure queries and the novel 'molecular scissoring' query. MMsINC is interfaced with other primary data collectors, such as PubChem, Protein Data Bank (PDB), the Food and Drug Administration database of approved drugs and ZINC.

REST/NRSF is a multifunctional transcription factor that represses or silences many neuron-specific genes in both neural and non-neural cells by recruitment to its cognate RE1/NRSE regulatory sites. An increase in RE1/NRSE genomic binding is found in Huntington's disease (HD), resulting in the repression of REST/NRSF regulated gene transcription, among which BDNF, thus representing one of the possible detrimental effectors in HD. Three 2-aminothiazole derivatives were recently identified as potent modulators of the RE1/NRSE silencing activity through a cell-based gene reporter assay. In this study, the structure-activity relationships (SAR) of a library of commercially available 2-aminoisothiazoles diversely substituted at the amino group or at position 4 has been evaluated. A quantitative structure-activity relationship analysis performed using the Phase strategy yielded highly predictive 3D-QSAR pharmacophore model for in silico drug screening.

Peritonitis complicating diverticular disease may be treated by sigmoid resection (with or without primary anastomosis) or by a conservative surgical approach, either laparoscopically or by open surgery. The choice depends on the severity of the peritonitis (Hinchey), the patient's conditions (ASA) and the surgeon's experience. Sigmoid resection with primary anastomosis has a lower morbidity and mortality vs Hartmann's procedure. After the introduction of laparoscopy in colorectal surgery, exploratory laparoscopy combined with drainage has been proposed to treat acute episodes, followed by laparoscopic resection. Since 1982, over 1000 patients have been operated on for colorectal disease: 119 for complicated diverticulitis, 55 of which complicated by peritonitis. In the latter, we performed conservative surgery (25 patients) and resection (30 patients) laparoscopically or by open surgery. Our results show a higher morbidity and mortality for the Hartmann procedure vs sigmoid resection with primary anastomosis and a lower specific morbidity in patients undergoing laparoscopic exploration and drainage. Moreover, there was a low percentage (52%) of re-canalisations with the Hartmann procedure, with a morbidity of 32% associated with this procedure. In conclusion, we believe that a conservative laparoscopic surgical approach may be advocated in selected cases (Hinchey II and III without clear perforation), followed by laparoscopic sigmoidectomy, resection with primary anastomosis in Hinchey I or in cases of evident perforation with purulent or faecal peritonitis (possibly combined with a stoma), reserving the Hartmann procedure for compromised patients.

The authors present their initial experience with minimally invasive surgery of the oesophagus. Two emblematic cases of benign and malignant oesophageal pathology, treated by laparoscopic transhiatal oesophagectomy, are reported. The surgical procedure is described in detail and compared with other techniques reported in the literature. Its advantages in terms of postoperative pain and morbidity are stressed. The role of the procedure in neoplastic diseases of the oesophagus is still debated, above all with regard to the accuracy of the mediastinal node dissection.

The authors present two cases of inflammatory abdominal aortic aneurysms complicated by an aorto-caval fistula in one case and by an aorto-duodenal fistula in the other. The aetiology of such aneurysms is still debatable, the histological features are typical, and the preoperative clinical and instrumental findings do not always allow them to be differentiated from atherosclerotic ones. The complications may be rapidly fatal or asymptomatic, as in the case of occult aorto-caval fistulas. The surgical approach is technically more difficult because of inflammation and adhesions between the aneurysmal sac and contiguous structures which make dissection and aorto-iliac clamping more demanding. In the case of aorto-caval fistulas, clamping of the vein should be avoided and manual compression preferred to assure haemostasis; with aorto-enteric fistulas, it is necessary to reduce the contamination of the operative field. Ureterolysis is not always necessary in patients with hydronephrosis. Endovascular devices may be an alternative also for a number of complications.

The authors present two cases of toxic megacolon. A proctocolectomy with ileostomy was performed in both cases. One of the two cases was detected late; the patient underwent surgery in desperate conditions and died on postoperative day one. The authors regard a prompt diagnosis as fundamental for the correct timing of the surgical approach in patients with toxic megacolon and severe acute colitis. Clinical and laboratory findings and, above all, plain films of the abdomen should be evaluated very carefully. The authors suggest total proctocolectomy with section-suture of the rectal stump at levator level and ileostomy, which makes it possible to avoid postoperative rectal haemorrhages and recurrences.

Limb threatening ischemia is an acute step in the chronic course of peripheral arterial obstructive disease that requires some form of intervention. The objective of this study is to prove that reconstructive surgery as well as non reconstructive approaches are associated with positive results. Ours is a retrospective analysis of a ten year experience in the treatment of limb threatening ischemia. In the period 1983-93, 139 patients under-went 164 procedures. 67% of patients were diabetics. Early in the observation period the therapeutic strategy was non reconstructive, the procedure of choice was sympathectomy. Later vascular reconstructions have been recognized as the procedures of choice. In the cases not amenable to reconstructive procedures according to our group criteria (absence of a tibial vessel in continuity with a patent pedal arch), we have employed procedures such as prostanoid infusion; thrombolysis and epidural spinal cord stimulation. Reconstructive procedures have been associated with a decrease in the number of major amputations. Alternative procedures, employed in patients not amenable to reconstruction have proven worthwhile in terms of limb salvage even if this trend is limited to a short period of observation.