Reconstruction and analysis of the NF-kB pathway interactome: a systems biology approach

Background. One of the phenomena observed in human aging is the progressive increase of a systemic inflammatory state, a condition referred to as “inflammaging”, negatively correlated with longevity. A prominent mediator of inflammation is the transcription factor NF-kB, that acts as key transcriptional regulator of many genes coding for pro-inflammatory cytokines. Many different signaling pathways activated by very diverse stimuli converge on NF-kB, resulting in a regulatory network characterized by high complexity. NF-kB signaling has been proposed to be responsible of inflammaging. Scope of this analysis is to provide a wider, systemic picture of such intricate signaling and interaction network: the NF-kB pathway interactome. Methods. The study has been carried out following a workflow for gathering information from literature as well as from several pathway and protein interactions databases, and for integrating and analyzing existing data and the relative reconstructed representations by using the available computational tools. Strong manual intervention has been necessarily used to integrate data from multiple sources into mathematically analyzable networks. The reconstruction of the NF-kB interactome pursued with this approach provides a starting point for a general view of the architecture and for a deeper analysis and understanding of this complex regulatory system. Results. A “core” and a “wider” NF-kB pathway interactome, consisting of 140 and 3146 proteins respectively, were reconstructed and analyzed through a mathematical, graph-theoretical approach. Among other interesting features, the topological characterization of the interactomes shows that a relevant number of interacting proteins are in turn products of genes that are controlled and regulated in their expression exactly by NF-kB transcription factors. These “feedback loops”, not always well-known, deserve deeper investigation since they may have a role in tuning the response and the output consequent to NF-kB pathway initiation, in regulating the intensity of the response, or its homeostasis and balance in order to make the functioning of such critical system more robust and reliable. This integrated view allows to shed light on the functional structure and on some of the crucial nodes of thet NF-kB transcription factors interactome. Conclusion. Framing structure and dynamics of the NF-kB interactome into a wider, systemic picture would be a significant step toward a better understanding of how NF-kB globally regulates diverse gene programs and phenotypes. This study represents a step towards a more complete and integrated view of the NF-kB signaling system.

Preliminary studies for the establishment of a tsunami early detection algorithm to be used in the frame of a Tsunami Warning System

This thesis presents and discusses TEDA, an algorithm for the automatic detection in real-time of tsunamis and large amplitude waves on sea level records. TEDA has been developed in the frame of the Tsunami Research Team of the University of Bologna for coastal tide gauges and it has been calibrated and tested for the tide gauge station of Adak Island, in Alaska. A preliminary study to apply TEDA to offshore buoys in the Pacific Ocean is also presented.

Sonic Hedgehog pathway impairment in Neural Precursor Cells of the Ts65Dn mouse, an animal model of Down syndrome

Mental retardation in Down syndrome (DS) has been imputed to the decreased brain volume, which is evident starting from the early phases of development. Recent studies in a widely used mouse model of DS, the Ts65Dn mouse, have shown that neurogenesis is severely impaired during the early phases of brain development, suggesting that this defect may be a major determinant of brain hypotrophy and mental retardation in individuals with DS. Recently, it has been found that in the cerebellum of Ts65Dn mice there is a defective responsiveness to Sonic Hedgehog (Shh), a potent mitogen that controls cell division during brain development, suggesting that failure of Shh signaling may underlie the reduced proliferation potency in DS. Based on these premises, we sought to identify the molecular mechanisms underlying derangement of the Shh pathway in neural precursor cells (NPCs) from Ts65Dn mice. We found that the expression levels of the Shh receptor Patched1 (Ptch1) were increased compared to controls both at the RNA and protein level. Partial silencing of Ptch1 expression in trisomic NPCs restored cell proliferation, indicating that proliferation impairment was due to Ptch1 overexpression. We further found that the overexpression of Ptch1 in trisomic NPCs is related to increased levels of AICD, a transcription-promoting fragment of amyloid precursor protein (APP). Increased AICD binding to the Ptch1 promoter favored its acetylated status, thus enhancing Ptch1 expression. Taken together, these data provide novel evidence that Ptch1 over expression underlies derangement of the Shh pathway in trisomic NPCs, with consequent proliferation impairment. The demonstration that Ptch1 over expression in trisomic NPCs is due to an APP fragment provides a link between this trisomic gene and the defective neuronal production that characterizes the DS brain.

adaptive capabilities of the pi3k/akt/mtor pathway in acute myeloid leukemia revealed by the use of selective inhibitors

Because of its aberrant activation, the PI3K/AKT/mTOR signaling pathway represents a pharmacological target in blast cells from patients with acute myelogenous leukemia (AML). Using Reverse Phase Protein Microarrays (RPMA), we have analyzed 20 phosphorylated epitopes of the PI3K/Akt/mTor signal pathway of peripheral blood and bone marrow specimens of 84 patients with newly diagnosed AML. Fresh blast cells were grown for 2 h, 4 h or 20 h untreated or treated with a panel of phase I or phase II Akt allosteric inhibitors, either alone or in combination with the mTOR kinase inhibitor Torin1 or the broad RTK inhibitor Sunitinib. By unsupervised hierarchical clustering a strong phosphorylation/activity of most of the sampled members of the PI3K/Akt/mTOR pathway was observed in 70% of samples from AML patients. Remarkably, however, we observed that inhibition of Akt phosphorylation, as well as of its substrates, was transient, and recovered or even increased far above basal level after 20 h in 60% samples. We demonstrated that inhibition of Akt induces FOXO-dependent insulin receptor expression and IRS-1 activation, attenuating the effect of drug treatment by reactivation of PI3K/Akt. Consistent with this model we found that combined inhibition of Akt and RTKs is much more effective than either alone, revealing the adaptive capabilities of signaling networks in blast cells and highliting the limations of these drugs if used as monotherapy.

A parallel community detection algorithm

Complex networks analysis is a very popular topic in computer science. Unfortunately this networks, extracted from different contexts, are usually very large and the analysis may be very complicated: computation of metrics on these structures could be very complex. Among all metrics we analyse the extraction of subnetworks called communities: they are groups of nodes that probably play the same role within the whole structure. Communities extraction is an interesting operation in many different fields (biology, economics,...). In this work we present a parallel community detection algorithm that can operate on networks with huge number of nodes and edges. After an introduction to graph theory and high performance computing, we will explain our design strategies and our implementation. Then, we will show some performance evaluation made on a distributed memory architectures i.e. the supercomputer IBM-BlueGene/Q "Fermi" at the CINECA supercomputing center, Italy, and we will comment our results.

SMO inhibitor specifically targets the Hedgehog Pathway and reverts the drug-resistance of Leukemic Stem Cells

Abnormal Hedgehog signaling is associated with human malignancies. Smo, a key player of that signaling, is the most suitable target to inhibit this pathway. To this aim several molecules, antagonists of Smo, have been synthesized, and some of them have started the phase I in clinical trials. Our hospital participated to one of these studies which investigated the oral administration of a new selective inhibitor of Smo (SMOi). To evaluate ex vivo SMOi efficacy and to identify new potential clinical biomarkers of responsiveness, we separated bone marrow CD34+ cells from 5 acute myeloid leukemia (AML), 1 myelofibrosis (MF), 2 blastic phases chronic myeloid leukemia (CML) patients treated with SMOi by immunomagnetic separation, and we analysed their gene expression profile using Affimetrix HG-U133 Plus 2.0 platform. This analysis, showed differential expression after 28 days start of therapy (p-value ≤ 0.05) of 1,197 genes in CML patients and 589 genes in AML patients. This differential expression is related to Hedgehog pathway with a p-value = 0.003 in CML patients and with a p-value = 0.0002 in AML patients, suggesting that SMOi targets specifically this pathway. Among the genes differentially expressed we observed strong up-regulation of Gas1 and Kif27 genes, which may work as biomarkers of responsiveness of SMOi treatment in CML CD34+ cells whereas Hedgehog target genes (such as Smo, Gli1, Gli2, Gli3), Bcl2 and Abca2 were down-regulated, in both AML and CML CD34+ cells. It has been reported that Bcl-2 expression could be correlated with cancer therapy resistance and that Hedgehog signaling modulate ATP-binding (ABC) cassette transporters, whose expression has been correlated with chemoresistance. Moreover we confirmed that in vitro SMOi treatment targets Hedgehog pathway, down-regulate ABC transporters, Abcg2 and Abcb1 genes, and in combination with tyrosine kinase inhibitors (TKIs) could revert the chemoresistance mechanism in K562 TKIs-resistant cell line.

Coordination and Control of Autonomous Mobile Robots Swarms by using Particle Swarm Optimization Algorithm and Consensus Theory

This thesis presents some different techniques designed to drive a swarm of robots in an a-priori unknown environment in order to move the group from a starting area to a final one avoiding obstacles. The presented techniques are based on two different theories used alone or in combination: Swarm Intelligence (SI) and Graph Theory. Both theories are based on the study of interactions between different entities (also called agents or units) in Multi- Agent Systems (MAS). The first one belongs to the Artificial Intelligence context and the second one to the Distributed Systems context. These theories, each one from its own point of view, exploit the emergent behaviour that comes from the interactive work of the entities, in order to achieve a common goal. The features of flexibility and adaptability of the swarm have been exploited with the aim to overcome and to minimize difficulties and problems that can affect one or more units of the group, having minimal impact to the whole group and to the common main target. Another aim of this work is to show the importance of the information shared between the units of the group, such as the communication topology, because it helps to maintain the environmental information, detected by each single agent, updated among the swarm. Swarm Intelligence has been applied to the presented technique, through the Particle Swarm Optimization algorithm (PSO), taking advantage of its features as a navigation system. The Graph Theory has been applied by exploiting Consensus and the application of the agreement protocol with the aim to maintain the units in a desired and controlled formation. This approach has been followed in order to conserve the power of PSO and to control part of its random behaviour with a distributed control algorithm like Consensus.

Parallelization of the algorithm WHAM with NVIDIA CUDA.

The aim of my thesis is to parallelize the Weighting Histogram Analysis Method (WHAM), which is a popular algorithm used to calculate the Free Energy of a molucular system in Molecular Dynamics simulations. WHAM works in post processing in cooperation with another algorithm called Umbrella Sampling. Umbrella Sampling has the purpose to add a biasing in the potential energy of the system in order to force the system to sample a specific region in the configurational space. Several N independent simulations are performed in order to sample all the region of interest. Subsequently, the WHAM algorithm is used to estimate the original system energy starting from the N atomic trajectories. The parallelization of WHAM has been performed through CUDA, a language that allows to work in GPUs of NVIDIA graphic cards, which have a parallel achitecture. The parallel implementation may sensibly speed up the WHAM execution compared to previous serial CPU imlementations. However, the WHAM CPU code presents some temporal criticalities to very high numbers of interactions. The algorithm has been written in C++ and executed in UNIX systems provided with NVIDIA graphic cards. The results were satisfying obtaining an increase of performances when the model was executed on graphics cards with compute capability greater. Nonetheless, the GPUs used to test the algorithm is quite old and not designated for scientific calculations. It is likely that a further performance increase will be obtained if the algorithm would be executed in clusters of GPU at high level of computational efficiency. The thesis is organized in the following way: I will first describe the mathematical formulation of Umbrella Sampling and WHAM algorithm with their apllications in the study of ionic channels and in Molecular Docking (Chapter 1); then, I will present the CUDA architectures used to implement the model (Chapter 2); and finally, the results obtained on model systems will be presented (Chapter 3).

A multi-objective genetic algorithm optimisation using variable speed pumps in water distribution systems

Due to its practical importance and inherent complexity, the optimisation of distribution networks for supplying drinking water has been the subject of extensive study for the past 30 years. The optimization is governed by sizing the pipes in the water distribution network (WDN) and / or optimises specific parts of the network such as pumps, tanks etc. or try to analyse and optimise the reliability of a WDN. In this thesis, the author has analysed two different WDNs (Anytown City and Cabrera city networks), trying to solve and optimise a multi-objective optimisation problem (MOOP). The main two objectives in both cases were the minimisation of Energy Cost (€) or Energy consumption (kWh), along with the total Number of pump switches (TNps) during a day. For this purpose, a decision support system generator for Multi-objective optimisation used. Its name is GANetXL and has been developed by the Center of Water System in the University of Exeter. GANetXL, works by calling the EPANET hydraulic solver, each time a hydraulic analysis has been fulfilled. The main algorithm used, was a second-generation algorithm for multi-objective optimisation called NSGA_II that gave us the Pareto fronts of each configuration. The first experiment that has been carried out was the network of Anytown city. It is a big network with a pump station of four fixed speed parallel pumps that are boosting the water dynamics. The main intervention was to change these pumps to new Variable speed driven pumps (VSDPs), by installing inverters capable to diverse their velocity during the day. Hence, it’s been achieved great Energy and cost savings along with minimisation in the number of pump switches. The results of the research are thoroughly illustrated in chapter 7, with comments and a variety of graphs and different configurations. The second experiment was about the network of Cabrera city. The smaller WDN had a unique FS pump in the system. The problem was the same as far as the optimisation process was concerned, thus, the minimisation of the energy consumption and in parallel the minimisation of TNps. The same optimisation tool has been used (GANetXL).The main scope was to carry out several and different experiments regarding a vast variety of configurations, using different pump (but this time keeping the FS mode), different tank levels, different pipe diameters and different emitters coefficient. All these different modes came up with a large number of results that were compared in the chapter 8. Concluding, it should be said that the optimisation of WDNs is a very interested field that has a vast space of options to deal with. This includes a large number of algorithms to choose from, different techniques and configurations to be made and different support system generators. The researcher has to be ready to “roam” between these choices, till a satisfactory result will convince him/her that has reached a good optimisation point.

Genus Bifidobacterium: taxonomy studies and gene expression analysis on folate pathway

Folates (vitamin B9) are essential water soluble vitamins, whose deficiency in humans may contribute to the onset of several diseases, such as anaemia, cancer, cardiovascular diseases, neurological problems as well as defects in embryonic development. Human and other mammals are unable to synthesize ex novo folate obtaining it from exogenous sources, via intestinal absorption. Recently the gut microbiota has been identified as an important source of folates and the selection and use of folate producing microorganisms represents an innovative strategy to increase human folate levels. The aim of this thesis was to gain a fundamental understanding of folate metabolism in Bifidobacterium adolescentis. The work was subdivided in three main phases, also aimed to solve different problems encountered working with Bifidobacterium strains. First, a new identification method (based on PCR-RFLP of hsp60 gene) was specifically developed to identify Bifidobacterium strains. Secondly, Bifidobacterium adolescentis biodiversity was explored in order to recognize representing strains of this species to be screened for their folate production ability. Results showed that this species is characterized by a wide variability and support the idea that a possible new taxonomic re-organization would be required. Finally B. adolescentis folate metabolism was studied using a double approach. A quantitative analysis of folate content was complemented by the examination of expression levels of genes involved in folate related pathways. For the normalization process, required to increase the robustness of the qRT-PCR analysis, an appropriate set of reference genes was tested using two different algorithms. Results demonstrate that B.adolescentis strains may represent an endogenous source of natural folate and they could be used to fortify fermented dairy products. This bio-fortification strategy presents many advantages for the consumer, providing native folate forms more bio-available, and not implicated in the discussed controversy concerning the safety of high intake of synthetic folic acid.

An optical sensor for online hematocrit measurement: characterization and fitting algorithm development

Lo scopo del presente lavoro di tesi riguarda la caratterizzazione di un sensore ottico per la lettura di ematocrito e lo sviluppo dell’algoritmo di calibrazione del dispositivo. In altre parole, utilizzando dati ottenuti da una sessione di calibrazione opportunamente pianificata, l’algoritmo sviluppato ha lo scopo di restituire la curva di interpolazione dei dati che caratterizza il trasduttore. I passi principali del lavoro di tesi svolto sono sintetizzati nei punti seguenti: 1) Pianificazione della sessione di calibrazione necessaria per la raccolta dati e conseguente costruzione di un modello black box.  Output: dato proveniente dal sensore ottico (lettura espressa in mV)  Input: valore di ematocrito espresso in punti percentuali ( questa grandezza rappresenta il valore vero di volume ematico ed è stata ottenuta con un dispositivo di centrifugazione sanguigna) 2) Sviluppo dell’algoritmo L’algoritmo sviluppato e utilizzato offline ha lo scopo di restituire la curva di regressione dei dati. Macroscopicamente, il codice possiamo distinguerlo in due parti principali: 1- Acquisizione dei dati provenienti da sensore e stato di funzionamento della pompa bifasica 2- Normalizzazione dei dati ottenuti rispetto al valore di riferimento del sensore e implementazione dell’algoritmo di regressione. Lo step di normalizzazione dei dati è uno strumento statistico fondamentale per poter mettere a confronto grandezze non uniformi tra loro. Studi presenti, dimostrano inoltre un mutazione morfologica del globulo rosso in risposta a sollecitazioni meccaniche. Un ulteriore aspetto trattato nel presente lavoro, riguarda la velocità del flusso sanguigno determinato dalla pompa e come tale grandezza sia in grado di influenzare la lettura di ematocrito.

The constitutive activation of the DNA damage response pathway is a novel therapeutic target in aggressive B-cell lymphoma

The recent finding that MYC-driven cancers are sensitive to inhibition of the DNA damage response (DDR) pathway, prompted us to investigate the role of DDR pathway as therapeutic target in diffuse large B-cell lymphoma (DLBCL), which frequently overexpresses the MYC oncogene. In a preliminary immunohistochemical study conducted on 99 consecutive DLBCL patients, we found that about half of DLBCLs showed constitutive expression of the phosphorylated forms of checkpoint kinases (CHK) and CDC25c, markers of DDR activation, and of phosphorylated histone H2AX (γH2AX), marker of DNA damage and genomic instability. Constitutive γH2AX expression correlated with c-MYC levels and DDR activation, and defined a subset of tumors characterised by poor outcome. Next, we used the CHK inhibitor PF-0477736 as a tool to investigate whether the inhibition of the DDR pathway might represent a novel therapeutic approach in DLBCL. Submicromolar concentrations of PF-0477736 hindered proliferation in DLBCL cell lines with activated DDR pathway. These results were fully recapitulated with a different CHK inhibitor (AZD-7762). Inhibition of checkpoint kinases induced rapid DNA damage accumulation and apoptosis in DLBCL cell lines and primary cells. These data suggest that pharmacologic inhibition of DDR through targeting of CHK kinases may represent a novel therapeutic strategy in DLBCL. The second part of this work is the clinical, molecular and functional description of a paradigmatic case of primary refractory Burkitt lymphoma characterized by spatial intratumor heterogeneity for the TP53 mutational status, high expression levels of genomic instability and DDR activation markers, primary resistance to chemotherapy and exquisite sensitivity to DDR inhibitors.

Study of PI3K/Akt signaling pathway as potential molecular target for T-cell acute lymphoblastic leukemia (T-ALL) treatment: pan-inhibition of PI3K catalitic isoforms as better therapeutic approach

Class I phosphatidylinositol 3-kinases (PI3Ks) are heterodimeric lipid kinases consisting of a regulatory subunit and one of four catalytic subunits (p110α, p110β, p110γ or p110δ). p110γ/p110δ PI3Ks are highly enriched in leukocytes. In general, PI3Ks regulate a variety of cellular processes including cell proliferation, survival and metabolism, by generating the second messenger phosphatidylinositol-3,4,5-trisphosphate (PtdIns(3,4,5)P3). Their activity is tightly regulated by the phosphatase and tensin homolog (PTEN) lipid phosphatase. PI3Ks are widely implicated in human cancers, and in particular are upregulated in T-cell acute lymphoblastic leukemia (T-ALL), mainly due to loss of PTEN function. These observations lend compelling weight to the application of PI3K inhibitors in the therapy of T-ALL. At present different compounds which target single or multiple PI3K isoforms have entered clinical trials. In the present research, it has been analyzed the therapeutic potential of the pan-PI3K inhibitor BKM120, an orally bioavailable 2,6-dimorpholino pyrimidine derivative, which has entered clinical trials for solid tumors, on both T-ALL cell lines and patient samples. BKM120 treatment resulted in cell cycle arrest and apoptosis, being cytotoxic to a panel of T-ALL cell lines and patient T-lymphoblasts. Remarkably, BKM120 synergized with chemotherapeutic agents currently used for treating T-ALL patients. BKM120 efficacy was confirmed in in vivo studies to a subcutaneous xenotransplant model of human T-ALL. Because it is still unclear which agents among isoform-specific or pan inhibitors can achieve the greater efficacy, further analyses have been conducted to investigate the effects of PI3K inhibition, in order to elucidate the mechanisms responsible for the proliferative impairment of T-ALL. Overall, these results indicated that BKM120 may be an efficient treatment for T-ALLs that have aberrant up-regulation of the PI3K signaling pathway and strongly support clinical application of pan-class I PI3K rather than single-isoform inhibitors in T-ALL treatment.

Image quality and dose evaluation of filtered back projection versus iterative reconstruction algorithm in multislice computed tomography.

Il presente lavoro di tesi è stato svolto presso il servizio di Fisica Sanitaria del Policlinico Sant'Orsola-Malpighi di Bologna. Lo studio si è concentrato sul confronto tra le tecniche di ricostruzione standard (Filtered Back Projection, FBP) e quelle iterative in Tomografia Computerizzata. Il lavoro è stato diviso in due parti: nella prima è stata analizzata la qualità delle immagini acquisite con una CT multislice (iCT 128, sistema Philips) utilizzando sia l'algoritmo FBP sia quello iterativo (nel nostro caso iDose4). Per valutare la qualità delle immagini sono stati analizzati i seguenti parametri: il Noise Power Spectrum (NPS), la Modulation Transfer Function (MTF) e il rapporto contrasto-rumore (CNR). Le prime due grandezze sono state studiate effettuando misure su un fantoccio fornito dalla ditta costruttrice, che simulava la parte body e la parte head, con due cilindri di 32 e 20 cm rispettivamente. Le misure confermano la riduzione del rumore ma in maniera differente per i diversi filtri di convoluzione utilizzati. Lo studio dell'MTF invece ha rivelato che l'utilizzo delle tecniche standard e iterative non cambia la risoluzione spaziale; infatti gli andamenti ottenuti sono perfettamente identici (a parte le differenze intrinseche nei filtri di convoluzione), a differenza di quanto dichiarato dalla ditta. Per l'analisi del CNR sono stati utilizzati due fantocci; il primo, chiamato Catphan 600 è il fantoccio utilizzato per caratterizzare i sistemi CT. Il secondo, chiamato Cirs 061 ha al suo interno degli inserti che simulano la presenza di lesioni con densità tipiche del distretto addominale. Lo studio effettuato ha evidenziato che, per entrambi i fantocci, il rapporto contrasto-rumore aumenta se si utilizza la tecnica di ricostruzione iterativa. La seconda parte del lavoro di tesi è stata quella di effettuare una valutazione della riduzione della dose prendendo in considerazione diversi protocolli utilizzati nella pratica clinica, si sono analizzati un alto numero di esami e si sono calcolati i valori medi di CTDI e DLP su un campione di esame con FBP e con iDose4. I risultati mostrano che i valori ricavati con l'utilizzo dell'algoritmo iterativo sono al di sotto dei valori DLR nazionali di riferimento e di quelli che non usano i sistemi iterativi.