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.

Studio del pathway della Displasia Ectodermica Anidrotica (EDA)

Anhidrotic Ectodermal Dysplasia (EDA), is the most frequent form among Ectodermal Dysplasias, hereditary genetic disorders causing ectodermal appendages defective development. Indeed, EDA is characterized by defective formation of hair follicles, sweat glands and teeth both in human patients and animals. EDA, the gene mutated in Anhidrotic Ectodermal Dysplasia, encodes Ectodysplasin, a TNF family member that activates NF-kB mediated transcription. This disease can occur with mutations in other EDA-NF-kB pathway members, as EDA receptor, EDAR and its adapter, EDARADD. Moreover, mutations in TRAF6, NEMO, IKB and NF-kBs genes are responsible for Immunodeficiency associated EDA (EDA-ID). Several molecules, as SHH, WNT/DKK, BMP and LTβ, have already been reported to be EDA pathway regulators or effectors although the knowledge of the full spectrum of EDA targets remains incomplete. During the first part of the research project a gene expression analysis was performed in primary keratinocytes from Wild-type and Tabby (EDA model mouse) mice to identify novel EDA target genes. Earlier expression profiling at various developmental time points in Tabby and Wild-type mouse skin reported genes differentially expressed in the two samples and, to increase the resolution to find genes whose expression may be restricted to epidermal cells, the study was extended to primary keratinocyte cultures established from E19 Wild-type and Tabby skin. Using microarrays bearing 44,000 gene probes, we found 385 “preliminary candidate” genes whose expression was significantly affected by Eda defect. By comparing expression profiles to those from Eda-A1 (where Eda-A1 is highly expressed) transgenic skin, we restricted the list to 38 “candidate EDA targets”, 14 of which were already known to be expressed in hair follicles or epidermis. This work confirmed expression changes for 3 selected genes, Tbx1, Bmp7, and Jag1, both in primary keratinocytes and in Wild-type and Tabby whole skin, by Q-PCR and Western blotting analyses. Thus, this study detected novel candidate pathways downstream of EDA. In the second part of the research project, plasmid constructs were produced and analyzed to create a transgenic mouse model for Immunodeficiency associated EDA disease (XL-EDA-ID). In particular, plasmids containing mouse Wild-type and mutated Nemo cDNA under K-17 epidermis-specific promoter control and a Flag tag, were prepared, on the way to confine transgene expression to mice epidermis and to determine EDA phenotype without immunodeficiency for a comparison to Tabby model phenotype. EDA-ID mutations reported in patients and selected for this study are: C417R (C409R in mouse), causing Zinc Finger protein domain destabilization and A288G (A282G in mouse) affecting oligomerization of the protein. Moreover, the ex-novo mutation, ZnF, C-terminal Zinc Finger domain deletion, was tested. Thus, the constructs were analyzed by transient transfection, Western blotting and luciferase assays techniques, detecting Nemo Wild-type and mutant protein products and residue NF-kB activity in presence of mutants, after TNF stimulation. In particular, MEF_Nemo-/- cell line was used to monitor NF-kB activity without endogenous Nemo gene. Results show reduced NF-kB activity in presence of mutated Nemo forms compared to Wild-type: 81% for A282G (A288G in human); 24% for C409R (C417R in human); 15% for ZnF. C409R mutation (C417R in human), reported in 6 EDA-ID human patients, was selected to prepare transgenic model mouse. Mice (white, FVP) born following K17-promoter-Flag-Nemo_C409R plasmid region pronuclear injection, were analyzed for the transgene presence in the genotype and a preliminar examination of their phenotype was performed. In particular, one mouse showed considerable coat defects if compared to Wild-type mice. This preliminar analysis suggests a possible influence of Nemo mutant over-expression in epidermis without immunodeficiency. Still, more microscopic studies to analyze hair subtypes, Guard, Awl and Zigzag (usually alterated inTabby mouse model), Immunohistochemistry experiments to detect epidermis restricted Nemo expression and sweat glands analysis, will follow. This and other transgene positive mice will be crossed with black mice C57BL6 to obtain at least two indipendent agouti lines to analyze. Theses mice will be used in EDA target genes detection through microarrays. Following, plasmid constructs containing other Nemo mutant forms (A282G and ZnF) might be studied by the same experimental approaches to prepare more transgenic model mice to compare to Nemo_C409R and Tabby mouse models.

Studio molecolare dei meccanismi dell'autoincompatibilità gametofitica in pero europeo (Pyrus communis)

Self-incompatibility (SI) systems have evolved in many flowering plants to prevent self-fertilization and thus promote outbreeding. Pear and apple, as many of the species belonging to the Rosaceae, exhibit RNase-mediated gametophytic self-incompatibility, a widespread system carried also by the Solanaceae and Plantaginaceae. Pear orchards must for this reason contain at least two different cultivars that pollenize each other; to guarantee an efficient cross-pollination, they should have overlapping flowering periods and must be genetically compatible. This compatibility is determined by the S-locus, containing at least two genes encoding for a female (pistil) and a male (pollen) determinant. The female determinant in the Rosaceae, Solanaceae and Plantaginaceae system is a stylar glycoprotein with ribonuclease activity (S-RNase), that acts as a specific cytotoxin in incompatible pollen tubes degrading cellular RNAs. Since its identification, the S-RNase gene has been intensively studied and the sequences of a large number of alleles are available in online databases. On the contrary, the male determinant has been only recently identified as a pollen-expressed protein containing a F-box motif, called S-Locus F-box (abbreviated SLF or SFB). Since F-box proteins are best known for their participation to the SCF (Skp1 - Cullin - F-box) E3 ubiquitine ligase enzymatic complex, that is involved in protein degradation through the 26S proteasome pathway, the male determinant is supposed to act mediating the ubiquitination of the S-RNases, targeting them for the degradation in compatible pollen tubes. Attempts to clone SLF/SFB genes in the Pyrinae produced no results until very recently; in apple, the use of genomic libraries allowed the detection of two F-box genes linked to each S haplotype, called SFBB (S-locus F-Box Brothers). In Japanese pear, three SFBB genes linked to each haplotype were cloned from pollen cDNA. The SFBB genes exhibit S haplotype-specific sequence divergence and pollen-specific expression; their multiplicity is a feature whose interpretation is unclear: it has been hypothesized that all of them participate in the S-specific interaction with the RNase, but it is also possible that only one of them is involved in this function. Moreover, even if the S locus male and female determinants are the only responsible for the specificity of the pollen-pistil recognition, many other factors are supposed to play a role in GSI; these are not linked to the S locus and act in a S-haplotype independent manner. They can have a function in regulating the expression of S determinants (group 1 factors), modulating their activity (group 2) or acting downstream, in the accomplishment of the reaction of acceptance or rejection of the pollen tube (group 3). This study was aimed to the elucidation of the molecular mechanism of GSI in European pear (Pyrus communis) as well as in the other Pyrinae; it was divided in two parts, the first focusing on the characterization of male determinants, and the second on factors external to the S locus. The research of S locus F-box genes was primarily aimed to the identification of such genes in European pear, for which sequence data are still not available; moreover, it allowed also to investigate about the S locus structure in the Pyrinae. The analysis was carried out on a pool of varieties of the three species Pyrus communis (European pear), Pyrus pyrifolia (Japanese pear), and Malus × domestica (apple); varieties carrying S haplotypes whose RNases are highly similar were chosen, in order to check whether or not the same level of similarity is maintained also between the male determinants. A total of 82 sequences was obtained, 47 of which represent the first S-locus F-box genes sequenced from European pear. The sequence data strongly support the hypothesis that the S locus structure is conserved among the three species, and presumably among all the Pyrinae; at least five genes have homologs in the analysed S haplotypes, but the number of F-box genes surrounding the S-RNase could be even greater. The high level of sequence divergence and the similarity between alleles linked to highly conserved RNases, suggest a shared ancestral polymorphism also for the F-box genes. The F-box genes identified in European pear were mapped on a segregating population of 91 individuals from the cross 'Abbé Fétel' × 'Max Red Bartlett'. All the genes were placed on the linkage group 17, where the S locus has been placed both in pear and apple maps, and resulted strongly associated to the S-RNase gene. The linkage with the RNase was perfect for some of the F-box genes, while for others very rare single recombination events were identified. The second part of this study was focused on the research of other genes involved in the SI response in pear; it was aimed on one side to the identification of genes differentially expressed in compatible and incompatible crosses, and on the other to the cloning and characterization of the transglutaminase (TGase) gene, whose role may be crucial in pollen rejection. For the identification of differentially expressed genes, controlled pollinations were carried out in four combinations (self pollination, incompatible, half-compatible and fully compatible cross-pollination); expression profiles were compared through cDNA-AFLP. 28 fragments displaying an expression pattern related to compatibility or incompatibility were identified, cloned and sequenced; the sequence analysis allowed to assign a putative annotation to a part of them. The identified genes are involved in very different cellular processes or in defense mechanisms, suggesting a very complex change in gene expression following the pollen/pistil recognition. The pool of genes identified with this technique offers a good basis for further study toward a better understanding of how the SI response is carried out. Among the factors involved in SI response, moreover, an important role may be played by transglutaminase (TGase), an enzyme involved both in post-translational protein modification and in protein cross-linking. The TGase activity detected in pear styles was significantly higher when pollinated in incompatible combinations than in compatible ones, suggesting a role of this enzyme in the abnormal cytoskeletal reorganization observed during pollen rejection reaction. The aim of this part of the work was thus to identify and clone the pear TGase gene; the PCR amplification of fragments of this gene was achieved using primers realized on the alignment between the Arabidopsis TGase gene sequence and several apple EST fragments; the full-length coding sequence of the pear TGase gene was then cloned from cDNA, and provided a precious tool for further study of the in vitro and in vivo action of this enzyme.

Computational investigation of the Plasmodium falciparum fatty acid biosynthetic pathway toward the discovery of novel antimalarials

The structural peculiarities of a protein are related to its biological function. In the fatty acid elongation cycle, one small carrier protein shuttles and delivers the acyl intermediates from one enzyme to the other. The carrier has to recognize several enzymatic counterparts, specifically interact with each of them, and finally transiently deliver the carried substrate to the active site. Carry out such a complex game requires the players to be flexible and efficiently adapt their structure to the interacting protein or substrate. In a drug discovery effort, the structure-function relationships of a target system should be taken into account to optimistically interfere with its biological function. In this doctoral work, the essential role of structural plasticity in key steps of fatty acid biosynthesis in Plasmodium falciparum is investigated by means of molecular simulations. The key steps considered include the delivery of acyl substrates and the structural rearrangements of catalytic pockets upon ligand binding. The ground-level bases for carrier/enzyme recognition and interaction are also put forward. The structural features of the target have driven the selection of proper drug discovery tools, which captured the dynamics of biological processes and could allow the rational design of novel inhibitors. The model may be perspectively used for the identification of novel pathway-based antimalarial compounds.

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.

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.

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.

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.