Computational strategies to include protein flexibility in Ligand Docking and Virtual Screening

The dynamic character of proteins strongly influences biomolecular recognition mechanisms. With the development of the main models of ligand recognition (lock-and-key, induced fit, conformational selection theories), the role of protein plasticity has become increasingly relevant. In particular, major structural changes concerning large deviations of protein backbones, and slight movements such as side chain rotations are now carefully considered in drug discovery and development. It is of great interest to identify multiple protein conformations as preliminary step in a screening campaign. Protein flexibility has been widely investigated, in terms of both local and global motions, in two diverse biological systems. On one side, Replica Exchange Molecular Dynamics has been exploited as enhanced sampling method to collect multiple conformations of Lactate Dehydrogenase A (LDHA), an emerging anticancer target. The aim of this project was the development of an Ensemble-based Virtual Screening protocol, in order to find novel potent inhibitors. On the other side, a preliminary study concerning the local flexibility of Opioid Receptors has been carried out through ALiBERO approach, an iterative method based on Elastic Network-Normal Mode Analysis and Monte Carlo sampling. Comparison of the Virtual Screening performances by using single or multiple conformations confirmed that the inclusion of protein flexibility in screening protocols has a positive effect on the probability to early recognize novel or known active compounds.

Ruolo dell'Enantiomero (R)-9-HSA nel controllo della proliferazione in una linea di Adenocarcinoma del Colon umano

9-hydroxystearic acid (9-HSA) is an endogenous lipoperoxidation product and its administration to HT29, a colon adenocarcinoma cell line, induced a proliferative arrest in G0/G1 phase mediated by a direct activation of the p21WAF1 gene, bypassing p53. We have previously shown that 9-HSA controls cell growth and differentiation by inhibiting histone deacetylase 1 (HDAC1) activity, showing interesting features as a new anticancer drug. The interaction of 9-HSA with the catalytic site of the 3D model has been tested with a docking procedure: noticeably, when interacting with the site, the (R)-9-enantiomer is more stable than the (S) one. Thus, in this study, (R)- and (S)-9-HSA were synthesized and their biological activity tested in HT29 cells. At the concentration of 50 M (R)-9-HSA showed a stronger antiproliferative effect than the (S) isomer, as indicated by the growth arrest in G0/G1. The inhibitory effect of (S)-9-HSA on HDAC1, HDAC2 and HDAC3 activity was less effective than that of the (R)-9-HSA in vitro, and the inhibitory activity of both the (R)- and the (S)-9-HSA isomer, was higher on HDAC1 compared to HDAC2 and HDAC3, thus demonstrating the stereospecific and selective interaction of 9-HSA with HDAC1. In addition, histone hyperacetylation caused by 9-HSA treatment was examined by an innovative HPLC/ESI/MS method. Analysis on histones isolated from control and treated HT29 confirmed the higher potency of (R)-9-HSA compared to (S)-9-HSA, severely affecting H2A-2 and H4 acetylation. On the other side, it seemed of interest to determine whether the G0/G1 arrest of HT29 cell proliferation could be bypassed by the stimulation with the growth factor EGF. Our results showed that 9-HSA-treated cells were not only prevented from proliferating, but also showed a decreased [3H]thymidine incorporation after EGF stimulation. In this condition, HT29 cells expressed very low levels of cyclin D1, that didn’t colocalize with HDAC1. These results suggested that the cyclin D1/HDAC1 complex is required for proliferation. Furthermore, in the effort of understanding the possible mechanisms of this effect, we have analyzed the degree of internalization of the EGF/EGFR complex and its interactions with HDAC1. EGF/EGFR/HDAC1 complex quantitatively increases in 9-HSA-treated cells but not in serum starved cells after EGF stimulation. Our data suggested that 9-HSA interaction with the catalytic site of the HDAC1 disrupts the HDAC1/cyclin D1 complex and favors EGF/EGFR recruitment by HDAC1, thus enhancing 9-HSA antiproliferative effects. In conclusion 9-HSA is a promising HDAC inhibitor with high selectivity and specificity, capable of inducing cell cycle arrest and histone hyperacetylation, but also able to modulate HDAC1 protein interaction. All these aspects may contribute to the potency of this new antitumor agent.

Molecular interactions: metal ions and protein chaperones in the urease system from Helicobacter pylori

Although nickel is a toxic metal for living organisms in its soluble form, its importance in many biological processes recently emerged. In this view, the investigation of the nickel-dependent enzymes urease and [NiFe]-hydrogenase, especially the mechanism of nickel insertion into their active sites, represent two intriguing case studies to understand other analogous systems and therefore to lead to a comprehension of the nickel trafficking inside the cell. Moreover, these two enzymes have been demonstrated to ensure survival and colonization of the human pathogen H. pylori, the only known microorganism able to proliferate in the gastric niche. The right nickel delivering into the urease active site requires the presence of at least four accessory proteins, UreD, UreE, UreF and UreG. Similarly, analogous process is principally mediated by HypA and HypB proteins in the [NiFe]-hydrogenase system. Indeed, HpHypA and HpHypB also have been proposed to act in the activation of the urease enzyme from H. pylori, probably mobilizing nickel ions from HpHypA to the HpUreE-HpUreG complex. A complete comprehension of the interaction mechanism between the accessory proteins and the crosstalk between urease and hydrogenase accessory systems requires the determination of the role of each protein chaperone that strictly depends on their structural and biochemical properties. The availability of HpUreE, HpUreG and HpHypA proteins in a pure form is a pre-requisite to perform all the subsequent protein characterizations, thus their purification was the first aim of this work. Subsequently, the structural and biochemical properties of HpUreE were investigated using multi-angle and quasi-elastic light scattering, as well as NMR and circular dichroism spectroscopy. The thermodynamic parameters of Ni2+ and Zn2+ binding to HpUreE were principally established using isothermal titration calorimetry and the importance of key histidine residues in the process of binding metal ions was studied using site-directed mutagenesis. The molecular details of the HpUreE-HpUreG and HpUreE-HpHypA protein-protein assemblies were also elucidated. The interaction between HpUreE and HpUreG was investigated using ITC and NMR spectroscopy, and the influence of Ni2+ and Zn2+ metal ions on the stabilization of this association was established using native gel electrophoresis, light scattering and thermal denaturation scanning followed by CD spectroscopy. Preliminary HpUreE-HpHypA interaction studies were conducted using ITC. Finally, the possible structural architectures of the two protein-protein assemblies were rationalized using homology modeling and docking computational approaches. All the obtained data were interpreted in order to achieve a more exhaustive picture of the urease activation process, and the correlation with the accessory system of the hydrogenase enzyme, considering the specific role and activity of the involved protein players. A possible function for Zn2+ in the chaperone network involved in Ni2+ trafficking and urease activation is also envisaged.

Stereoselective synthesis of heterocycles as intermediates of biologically active compounds

The aim of this thesis was to investigate the synthesis of enantiomerically enriched heterocycles and dehydro-β-amino acid derivatives which can be used as scaffolds or intermediates of biologically active compounds, in particular as novel αvβ3 and α5β1 integrin ligands. The starting materials of all the compounds here synthesized are alkylideneacetoacetates. Alkylidene derivates are very usefull compounds, they are usually used as unsaturated electrophiles and they have the advantage of introducing different kind of functionality that may be further elaborated. In chapter 1, regio- and stereoselective allylic amination of pure carbonates is presented. The reaction proceeds via uncatalyzed or palladium-catalyzed conditions and affords enantiopure dehydro-β-amino esters that are useful precursor of biologically active compounds. Chapter 2 illustrates the synthesis of substituted isoxazolidines and isoxazolines via Michael addition followed by intramolecular hemiketalisation. The investigation on the effect of the Lewis acid catalysis on the regioselectivity of the addition it also reported. Isoxazolidines and isoxazolines are interesting heterocyclic compounds that may be regarded as unusual constrained -amino acids or as furanose mimetics. The synthesis of unusual cyclic amino acids precursors, that may be envisaged as proline analogues, as scaffolds for the design of bioactive peptidomimetics is presented in chapter 3. The synthesis of 2-substituted-3,4-dehydropyrrole derivatives starting from allylic carbonates via a two step allylic amination/ring closing metathesis (RCM) protocol is carried out. The reaction was optimized by testing different Grubbs’ catalysts and carbamate nitrogen protecting groups. Moreover, in view of a future application of these dehydro-β-amino acids as central core of peptidomimetics , the malonate chain was also used to protect nitrogen prior to RCM. Finally, chapter 4 presents the synthesis of two novel different classes of integrin antagonists, one derived from dehydro-β-amino acid prepared as described in chapter 1 and the other one has isoxazolidines synthesized in chapter 2 as rigid constrained core. Since that these compounds are promising RGD mimetics for αvβ3 and α5β1 integrins, they have been submitted to biological assay. and to interpret on a molecular basis their different affinities for the αvβ3 receptor, docking studies were performed using Glide program.

Synthesis of Modified Amino Acids and Insertion in Peptides and Mimetics. Structural Aspects and Impact on Biological Activity.

I studied the effects exerted by the modifications on structures and biological activities of the compounds so obtained. I prepared peptide analogues containing unusual amino acids such as halogenated, alkylated (S)- or (R)-tryptophans, useful for the synthesis of mimetics of the endogenous opioid peptide endomorphin-1, or 2-oxo-1,3-oxazolidine-4-carboxylic acids, utilized as pseudo-prolines having a clear all-trans configuration of the preceding peptide bond. The latter gave access to a series of constrained peptidomimetics with potential interest in medicinal chemistry and in the field of the foldamers. In particular, I have dedicated much efforts to the preparation of cyclopentapeptides containing D-configured, alfa-, or beta-aminoacids, and also of cyclotetrapeptides including the retro-inverso modification. The conformational analyses confirmed that these cyclic compounds can be utilized as rigid scaffolds mimicking gamma- or beta-turns, allowing to generate new molecular and 3D diversity. Much work has been dedicated to the structural analysis in solution and in the receptor-bound state, fundamental for giving a rationale to the experimentally determined bioactivity, as well as for predicting the activity of virtual compounds (in silico pre-screen). The conformational analyses in solution has been done mostly by NMR (2D gCosy, Roesy, VT, molecular dynamics, etc.). A special section is dedicated to the prediction of plausible poses of the ligands when bound to the receptors by Molecular Docking. This computational method proved to be a powerful tool for the investigation of ligand-receptor interactions, and for the design of selective agonists and antagonists. Another practical use of cyclic peptidomimetics was the synthesis and biological evaluation of cyclic analogues of endomorphin-1 lacking in a protonable amino group. The studies revealed that a inverse type II beta-turn on D-Trp-Phe constituted the bioactive conformation.

Computational Methods in Biophysics and Medicinal Chemistry: Applications and Challenges

In this thesis I described the theory and application of several computational methods in solving medicinal chemistry and biophysical tasks. I pointed out to the valuable information which could be achieved by means of computer simulations and to the possibility to predict the outcome of traditional experiments. Nowadays, computer represents an invaluable tool for chemists. In particular, the main topics of my research consisted in the development of an automated docking protocol for the voltage-gated hERG potassium channel blockers, and the investigation of the catalytic mechanism of the human peptidyl-prolyl cis-trans isomerase Pin1.