OPTIMIZATION AND APPLICATION OF COMPUTATIONAL METHODS FOR THE DESIGN OF PROTEIN-PROTEIN INTERACTIONS MODULATORS

In the wide field of PPIs, this PhD project has been focused on the optimization and application of computational methods for the design of PPIs modulators, with a particular interest toward peptide modulators targeting PPIs involving helical motifs. In this contest, the first part of the project has been aimed to define the rationales behind the helical secondary structure stabilization and the helical screw sense selectivity exerted by chiral C\u3b1-tetrasubstituted amino acids (cCTAAs) through REMD simulations and QTAIM analyses, and the mechanisms responsible of the helical screw sense inversion through PNEB simulations. In detail, it has been found that the helical motif is stabilized by two complementary mechanisms: the first depends on the steric hindrance exerted by the cCTAA in an area parallel to the peptide helix axis and downstream of the cCTAA itself, whereas the second consists in the strengthening of the helical H-bond network thanks to peculiar C-H\ub7\ub7\ub7O=C interactions. Analogously, P-helical screw sense selectivity is ascribable to the cCTAA steric hindrance parallel to the peptide helix axis, without particular preferences for the region downstream and upstream of the cCTAA, together with quite strong noncovalent interactions, consisting of classical N \u2013 H\ub7\ub7\ub7O=C H-bonds and weak C \u2013 H\ub7\ub7\ub7O=C interactions. Furthermore, PNEB simulations performed on achiral peptides of different lengths suggest that the helical screw sense inversion requires the formation of \u3b3-turns, although a preferential screw sense inversion direction was not found. Therefore, the knowledge gained from these studies could be helpful in designing stable helical peptides, having a preferential screw sense and that can be in principle activated in situ by inducing a conformational switch from P to M helix or vice versa. Conversely, the second part of the project has been focused on the optimization of an MMGBSA based method, called Nwat-MMGBSA, aimed to improve the correlation between predicted binding energies of PPI complexes and experimental data. This approach, consisting in the inclusion, as part of the receptor, of hydration shells around the ligand during the MMGBSA calculations, was initially tested on classical receptor-ligand complexes and, then, automatized, optimized and tested on PPI complexes. This approach turned out to be good for the evaluation of PPI modulators activities, from different points of view. First of all, when water played a significant role in mediating protein-ligand interactions, the application of Nwat-MMGBSA improved the correlation between predicted and experimental data. On the other hand, if the solvent does not explicitly participate to the interaction, it did not give detrimental results compared to those obtained with the standard approach. In addition, the protocol proved to be robust and reproducible, giving equivalent results by using different setups. Furthermore, although an optimal number of water molecules to include in the hydration shell could not be found, in the case of PPI interactions inhibited by small molecules the inclusion of 50 \u2013 60 water molecules appears to be a good choice. A non-negligible advantage of this approach is represented by the possibility to automatize it, making it applicable for drug design/discovery purposes. Therefore, although further evaluations are needed, most of all on larger datasets, the knowledge coming from the combination of both parts of the project can be exploited for the design of stable non-natural peptides targeting PPIs

Mechanism of stabilization of helix secondary structure by constrained Cα-tetrasubstituted alpha-amino acids

The theoretical basis behind the ability of constrained C\u3b1-tetrasubstituted amino acids (CTAAs) to induce stable helical conformations has been studied through Replica Exchange Molecular Dynamics Potential of Mean Force Quantum Theory of Atoms In Molecules calculations on Ac-l-Ala-CTAA-l-Ala-Aib-l-Ala-NHMe peptide models. We found that the origin of helix stabilization by CTAAs can be ascribed to at least two complementary mechanisms limiting the backbone conformational freedom: steric hindrance predominantly in the (+x,+y,\u2013z) sector of a right-handed 3D Cartesian space, where the z axis coincides with the helical axis and the C\u3b1 of the CTAA lies on the +y axis (0,+y,0), and the establishment of additional and relatively strong C\u2013H\ub7\ub7\ub7O interactions involving the CTAA

Class II Phosphoinositide 3-Kinases Contribute to Endothelial Cells Morphogenesis

PMCID: PMC3539993This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited

Advances in the treatment of explicit water molecules in docking and binding free energy calculations

BACKGROUND: The inclusion of direct effects mediated by water during the ligand-receptor recognition is a hot-topic of modern computational chemistry applied to drug discovery and development. Docking or virtual screening with explicit hydration is still debatable, despite the successful cases that have been presented in the last years. Indeed, how to select the water molecules that will be included in the docking process or how the included waters should be treated remain open questions. OBJECTIVE: In this review, we will discuss some of the most recent methods that can be used in computational drug discovery and drug development when the effect of a single water, or of a small network of interacting waters, needs to be explicitly considered. RESULTS: Here, we analyse software to aid the selection, or to predict the position, of water molecules that are going to be explicitly considered in later docking studies. We also present software and protocols able to efficiently treat flexible water molecules during docking, including examples of applications. Finally, we discuss methods based on molecular dynamics simulations that can be used to integrate docking studies or to reliably and efficiently compute binding energies of ligands in presence of interfacial or bridging water molecules. CONCLUSIONS: Software applications aiding the design of new drugs that exploit water molecules, either as displaceable residues or as bridges to the receptor, are constantly being developed. Although further validation is needed, workflows that explicitly consider water will probably become a standard for computational drug discovery soon

An Efficient Implementation of the Nwat-MMGBSA Method to Rescore Docking Results in Medium-Throughput Virtual Screenings

Nwat-MMGBSA is a variant of MM-PB/GBSA based on the inclusion of a number of explicit water molecules that are the closest to the ligand in each frame of a molecular dynamics trajectory. This method demonstrated improved correlations between calculated and experimental binding energies in both protein-protein interactions and ligand-receptor complexes, in comparison to the standard MM-GBSA. A protocol optimization, aimed to maximize efficacy and efficiency, is discussed here considering penicillopepsin, HIV1-protease, and BCL-XL as test cases. Calculations were performed in triplicates on both classic HPC environments and on standard workstations equipped by a GPU card, evidencing no statistical differences in the results. No relevant differences in correlation to experiments were also observed when performing Nwat-MMGBSA calculations on 4 or 1 ns long trajectories. A fully automatic workflow for structure-based virtual screening, performing from library set-up to docking and Nwat-MMGBSA rescoring, has then been developed. The protocol has been tested against no rescoring or standard MM-GBSA rescoring within a retrospective virtual screening of inhibitors of AmpC \u3b2-lactamase and of the Rac1-Tiam1 protein-protein interaction. In both cases, Nwat-MMGBSA rescoring provided a statistically significant increase in the ROC AUCs of between 20 and 30%, compared to docking scoring or to standard MM-GBSA rescoring

Inositol 1,3,4,5,6-pentakisphosphate 2-kinase is a distant IPK member with a singular inositide binding site for axial 2-OH recognition

Inositol phosphates (InsPs) are signaling molecules with multiple roles in cells. In particular Graphic (InsP6) is involved in mRNA export and editing or chromatin remodeling among other events. InsP6 accumulates as mixed salts (phytate) in storage tissues of plants and plays a key role in their physiology. Human diets that are exclusively grain-based provide an excess of InsP6 that, through chelation of metal ions, may have a detrimental effect on human health. Ins(1,3,4,5,6)P5 2-kinase (InsP5 2-kinase or Ipk1) catalyses the synthesis of InsP6 from InsP5 and ATP, and is the only enzyme that transfers a phosphate group to the axial 2-OH of the myo-inositide. We present the first structure for an InsP5 2-kinase in complex with both substrates and products. This enzyme presents a singular structural region for inositide binding that encompasses almost half of the protein. The key residues in substrate binding are identified, with Asp368 being responsible for recognition of the axial 2-OH. This study sheds light on the unique molecular mechanism for the synthesis of the precursor of inositol pyrophosphates

Novel roles for class II Phosphoinositide 3-Kinase C2 beta in signalling pathways involved in prostate cancer cell invasion

Phosphoinositide 3-kinases (PI3Ks) regulate several cellular functions such as proliferation, growth, survival and migration. The eight PI3K isoforms are grouped into three classes and the three enzymes belonging to the class II subfamily (PI3K-C2a, ß and ?) are the least investigated amongst all PI3Ks. Interest on these isoforms has been recently fuelled by the identification of specific physiological roles for class II PI3Ks and by accumulating evidence indicating their involvement in human diseases. While it is now established that these isoforms can regulate distinct cellular functions compared to other PI3Ks, there is still a limited understanding of the signalling pathways that can be specifically regulated by class II PI3Ks. Here we show that PI3K-C2ß regulates mitogen-activated protein kinase kinase (MEK1/2) and extracellular signal-regulated kinase (ERK1/2) activation in prostate cancer (PCa) cells. We further demonstrate that MEK/ERK and PI3K-C2ß are required for PCa cell invasion but not proliferation. In addition we show that PI3K-C2ß but not MEK/ERK regulates PCa cell migration as well as expression of the transcription factor Slug. These data identify novel signalling pathways specifically regulated by PI3K-C2ß and they further identify this enzyme as a key regulator of PCa cell migration and invasion

A novel inhibitor of the PI3K/Akt pathway based on the structure of inositol 1,3,4,5,6-pentakisphosphate

Background: Owing to its role in cancer, the phosphoinositide 3-kinase (PI3K)/Akt pathway is an attractive target for therapeutic intervention. We previously reported that the inhibition of Akt by inositol 1,3,4,5,6- pentakisphosphate (InsP5) results in anti-tumour properties. To further develop this compound we modified its structure to obtain more potent inhibitors of the PI3K/Akt pathway.Methods: Cell proliferation/survival was determined by cell counting, sulphorhodamine or acridine orange/ethidium bromide assay; Akt activation was determined by western blot analysis. In vivo effect of compounds was tested on PC3 xenografts, whereas in vitro activity on kinases was determined by SelectScreen Kinase Profiling Service.Results: The derivative 2-O-benzyl-myo-inositol 1,3,4,5,6-pentakisphosphate (2-O-Bn-InsP5) is active towards cancer types resistant to InsP5 in vitro and in vivo. 2-O-Bn-InsP5 possesses higher pro-apoptotic activity than InsP 5 in sensitive cells and enhances the effect of anti-cancer compounds. 2-O-Bn-InsP5 specifically inhibits 3-phosphoinositide- dependent protein kinase 1 (PDK1) in vitro (IC 50 in the low nanomolar range) and the PDK1-dependent phosphorylation of Akt in cell lines and excised tumours. It is interesting to note that 2-O-Bn-InsP5 also inhibits the mammalian target of rapamycin (mTOR) in vitro.Conclusions: InsP5 and 2-O-Bn-InsP5 may represent lead compounds to develop novel inhibitors of the PI3K/Akt pathway (including potential dual PDK1/mTOR inhibitors) and novel potential anti-cancer drugs

Vascular Endothelial Growth Factor Receptor-3 Directly Interacts with Phosphatidylinositol 3-Kinase to Regulate Lymphangiogenesis

Background Dysfunctional lymphatic vessel formation has been implicated in a number of pathological conditions including cancer metastasis, lymphedema, and impaired wound healing. The vascular endothelial growth factor (VEGF) family is a major regulator of lymphatic endothelial cell (LEC) function and lymphangiogenesis. Indeed, dissemination of malignant cells into the regional lymph nodes, a common occurrence in many cancers, is stimulated by VEGF family members. This effect is generally considered to be mediated via VEGFR-2 and VEGFR-3. However, the role of specific receptors and their downstream signaling pathways is not well understood. Methods and Results Here we delineate the VEGF-C/VEGF receptor (VEGFR)-3 signaling pathway in LECs and show that VEGF-C induces activation of PI3K/Akt and MEK/Erk. Furthermore, activation of PI3K/Akt by VEGF-C/VEGFR-3 resulted in phosphorylation of P70S6K, eNOS, PLCc1, and Erk1/2. Importantly, a direct interaction between PI3K and VEGFR-3 in LECs was demonstrated both in vitro and in clinical cancer specimens. This interaction was strongly associated with the presence of lymph node metastases in primary small cell carcinoma of the lung in clinical specimens. Blocking PI3K activity abolished VEGF-C-stimulated LEC tube formation and migration. Conclusions Our findings demonstrate that specific VEGFR-3 signaling pathways are activated in LECs by VEGF-C. The importance of PI3K in VEGF-C/VEGFR-3-mediated lymphangiogenesis provides a potential therapeutic target for the inhibition of lymphatic metastasis

Two PI 3-Kinases and One PI 3-Phosphatase Together Establish the Cyclic Waves of Phagosomal PtdIns(3)P Critical for the Degradation of Apoptotic Cells

Cyclic oscillations in the level of phosphatidylinositol 3-phosphate in phagosomes, regulated by two phosphoinositide kinases and one phosphatase, are critical for phagosome maturation and degradation of apoptotic cells