Computational analyses on the structure-function relation in ion channels

Ion channels are protein molecules, embedded in the lipid bilayer of the cell membranes. They act as powerful sensing elements switching chemicalphysical stimuli into ion-fluxes. At a glance, ion channels are water-filled pores, which can open and close in response to different stimuli (gating), and one once open select the permeating ion species (selectivity). They play a crucial role in several physiological functions, like nerve transmission, muscular contraction, and secretion. Besides, ion channels can be used in technological applications for different purpose (sensing of organic molecules, DNA sequencing). As a result, there is remarkable interest in understanding the molecular determinants of the channel functioning. Nowadays, both the functional and the structural characteristics of ion channels can be experimentally solved. The purpose of this thesis was to investigate the structure-function relation in ion channels, by computational techniques. Most of the analyses focused on the mechanisms of ion conduction, and the numerical methodologies to compute the channel conductance. The standard techniques for atomistic simulation of complex molecular systems (Molecular Dynamics) cannot be routinely used to calculate ion fluxes in membrane channels, because of the high computational resources needed. The main step forward of the PhD research activity was the development of a computational algorithm for the calculation of ion fluxes in protein channels. The algorithm - based on the electrodiffusion theory - is computational inexpensive, and was used for an extensive analysis on the molecular determinants of the channel conductance. The first record of ion-fluxes through a single protein channel dates back to 1976, and since then measuring the single channel conductance has become a standard experimental procedure. Chapter 1 introduces ion channels, and the experimental techniques used to measure the channel currents. The abundance of functional data (channel currents) does not match with an equal abundance of structural data. The bacterial potassium channel KcsA was the first selective ion channels to be experimentally solved (1998), and after KcsA the structures of four different potassium channels were revealed. These experimental data inspired a new era in ion channel modeling. Once the atomic structures of channels are known, it is possible to define mathematical models based on physical descriptions of the molecular systems. These physically based models can provide an atomic description of ion channel functioning, and predict the effect of structural changes. Chapter 2 introduces the computation methods used throughout the thesis to model ion channels functioning at the atomic level. In Chapter 3 and Chapter 4 the ion conduction through potassium channels is analyzed, by an approach based on the Poisson-Nernst-Planck electrodiffusion theory. In the electrodiffusion theory ion conduction is modeled by the drift-diffusion equations, thus describing the ion distributions by continuum functions. The numerical solver of the Poisson- Nernst-Planck equations was tested in the KcsA potassium channel (Chapter 3), and then used to analyze how the atomic structure of the intracellular vestibule of potassium channels affects the conductance (Chapter 4). As a major result, a correlation between the channel conductance and the potassium concentration in the intracellular vestibule emerged. The atomic structure of the channel modulates the potassium concentration in the vestibule, thus its conductance. This mechanism explains the phenotype of the BK potassium channels, a sub-family of potassium channels with high single channel conductance. The functional role of the intracellular vestibule is also the subject of Chapter 5, where the affinity of the potassium channels hEag1 (involved in tumour-cell proliferation) and hErg (important in the cardiac cycle) for several pharmaceutical drugs was compared. Both experimental measurements and molecular modeling were used in order to identify differences in the blocking mechanism of the two channels, which could be exploited in the synthesis of selective blockers. The experimental data pointed out the different role of residue mutations in the blockage of hEag1 and hErg, and the molecular modeling provided a possible explanation based on different binding sites in the intracellular vestibule. Modeling ion channels at the molecular levels relates the functioning of a channel to its atomic structure (Chapters 3-5), and can also be useful to predict the structure of ion channels (Chapter 6-7). In Chapter 6 the structure of the KcsA potassium channel depleted from potassium ions is analyzed by molecular dynamics simulations. Recently, a surprisingly high osmotic permeability of the KcsA channel was experimentally measured. All the available crystallographic structure of KcsA refers to a channel occupied by potassium ions. To conduct water molecules potassium ions must be expelled from KcsA. The structure of the potassium-depleted KcsA channel and the mechanism of water permeation are still unknown, and have been investigated by numerical simulations. Molecular dynamics of KcsA identified a possible atomic structure of the potassium-depleted KcsA channel, and a mechanism for water permeation. The depletion from potassium ions is an extreme situation for potassium channels, unlikely in physiological conditions. However, the simulation of such an extreme condition could help to identify the structural conformations, so the functional states, accessible to potassium ion channels. The last chapter of the thesis deals with the atomic structure of the !- Hemolysin channel. !-Hemolysin is the major determinant of the Staphylococcus Aureus toxicity, and is also the prototype channel for a possible usage in technological applications. The atomic structure of !- Hemolysin was revealed by X-Ray crystallography, but several experimental evidences suggest the presence of an alternative atomic structure. This alternative structure was predicted, combining experimental measurements of single channel currents and numerical simulations. This thesis is organized in two parts, in the first part an overview on ion channels and on the numerical methods adopted throughout the thesis is provided, while the second part describes the research projects tackled in the course of the PhD programme. The aim of the research activity was to relate the functional characteristics of ion channels to their atomic structure. In presenting the different research projects, the role of numerical simulations to analyze the structure-function relation in ion channels is highlighted.

The posterior parietal cortex: a bridge between vision and action

The present work takes into account three posterior parietal areas, V6, V6A, and PEc, all operating on different subsets of signals (visual, somatic, motor). The work focuses on the study of their functional properties, to better understand their respective contribution in the neuronal circuits that make possible the interactions between subject and external environment. In the caudalmost pole of parietal lobe there is area V6. Functional data suggest that this area is related to the encoding of both objects motion and ego-motion. However, the sensitivity of V6 neurons to optic flow stimulations has been tested only in human fMRI experiments. Here we addressed this issue by applying on monkey the same experimental protocol used in human studies. The visual stimulation obtained with the Flow Fields stimulus was the most effective and powerful to activate area V6 in monkey, further strengthening this homology between the two primates. The neighboring areas, V6A and PEc, show different cytoarchitecture and connectivity profiles, but are both involved in the control of reaches. We studied the sensory responses present in these areas, and directly compared these.. We also studied the motor related discharges of PEc neurons during reaching movements in 3D space comparing also the direction and depth tuning of PEc cells with those of V6A. The results show that area PEc and V6A share several functional properties. Area PEc, unlike V6A, contains a richer and more complex somatosensory input, and a poorer, although complex visual one. Differences emerged also comparing the motor-related properties for reaches in depth: the incidence of depth modulations in PEc and the temporal pattern of modulation for depth and direction allow to delineate a trend among the two parietal visuomotor areas.

Identification of novel genetic alterations in pediatric cytogenetically normal acute myeloid leukemia by next-generation sequencing

Pediatric acute myeloid leukemia (AML) is a molecularly heterogeneous disease that arises from genetic alterations in pathways that regulate self-renewal and myeloid differentiation. While the majority of patients carry recurrent chromosomal translocations, almost 20% of childhood AML do not show any recognizable cytogenetic alteration and are defined as cytogenetically normal (CN)-AML. CN-AML patients have always showed a great variability in response to therapy and overall outcome, underlining the presence of unknown genetic changes, not detectable by conventional analyses, but relevant for pathogenesis, and outcome of AML. The development of novel genome-wide techniques such as next-generation sequencing, have tremendously improved our ability to interrogate the cancer genome. Based on this background, the aim of this research study was to investigate the mutational landscape of pediatric CN-AML patients negative for all the currently known somatic mutations reported in AML through whole-transcriptome sequencing (RNA-seq). RNA-seq performed on diagnostic leukemic blasts from 19 pediatric CN-AML cases revealed a considerable incidence of cryptic chromosomal rearrangements, with the identification of 21 putative fusion genes. Several of the fusion genes that were identified in this study are recurrent and might have a prognostic and/or therapeutic relevance. A paradigm of that is the CBFA2T3-GLIS2 fusion, which has been demonstrated to be a common alteration in pediatric CN-AML, predicting poor outcome. Important findings have been also obtained in the identification of novel therapeutic targets. On one side, the identification of NUP98-JARID1A fusion suggests the use of disulfiram; on the other, here we describe alteration-activating tyrosine kinases, providing functional data supporting the use of tyrosine kinase inhibitors to specifically inhibit leukemia cells. This study provides new insights in the knowledge of genetic alterations underlying pediatric AML, defines novel prognostic markers and putative therapeutic targets, and prospectively ensures a correct risk stratification and risk-adapted therapy also for the “all-neg” AML subgroup.

Transcriptional regulation of human mu-opioid receptor gene: functional characterization of activating and inhibitory transcription factors

The organization of the nervous and immune systems is characterized by obvious differences and striking parallels. Both systems need to relay information across very short and very long distances. The nervous system communicates over both long and short ranges primarily by means of more or less hardwired intercellular connections, consisting of axons, dendrites, and synapses. Longrange communication in the immune system occurs mainly via the ordered and guided migration of immune cells and systemically acting soluble factors such as antibodies, cytokines, and chemokines. Its short-range communication either is mediated by locally acting soluble factors or transpires during direct cell–cell contact across specialized areas called “immunological synapses” (Kirschensteiner et al., 2003). These parallels in intercellular communication are complemented by a complex array of factors that induce cell growth and differentiation: these factors in the immune system are called cytokines; in the nervous system, they are called neurotrophic factors. Neither the cytokines nor the neurotrophic factors appear to be completely exclusive to either system (Neumann et al., 2002). In particular, mounting evidence indicates that some of the most potent members of the neurotrophin family, for example, nerve growth factor (NGF) and brainderived neurotrophic factor (BDNF), act on or are produced by immune cells (Kerschensteiner et al., 1999) There are, however, other neurotrophic factors, for example the insulin-like growth factor-1 (IGF-1), that can behave similarly (Kermer et al., 2000). These factors may allow the two systems to “cross-talk” and eventually may provide a molecular explanation for the reports that inflammation after central nervous system (CNS) injury has beneficial effects (Moalem et al., 1999). In order to shed some more light on such a cross-talk, therefore, transcription factors modulating mu-opioid receptor (MOPr) expression in neurons and immune cells are here investigated. More precisely, I focused my attention on IGF-I modulation of MOPr in neurons and T-cell receptor induction of MOPr expression in T-lymphocytes. Three different opioid receptors [mu (MOPr), delta (DOPr), and kappa (KOPr)] belonging to the G-protein coupled receptor super-family have been cloned. They are activated by structurallyrelated exogenous opioids or endogenous opioid peptides, and contribute to the regulation of several functions including pain transmission, respiration, cardiac and gastrointestinal functions, and immune response (Zollner and Stein 2007). MOPr is expressed mainly in the central nervous system where it regulates morphine-induced analgesia, tolerance and dependence (Mayer and Hollt 2006). Recently, induction of MOPr expression in different immune cells induced by cytokines has been reported (Kraus et al., 2001; Kraus et al., 2003). The human mu-opioid receptor gene (OPRM1) promoter is of the TATA-less type and has clusters of potential binding sites for different transcription factors (Law et al. 2004). Several studies, primarily focused on the upstream region of the OPRM1 promoter, have investigated transcriptional regulation of MOPr expression. Presently, however, it is still not completely clear how positive and negative transcription regulators cooperatively coordinate cellor tissue-specific transcription of the OPRM1 gene, and how specific growth factors influence its expression. IGF-I and its receptors are widely distributed throughout the nervous system during development, and their involvement in neurogenesis has been extensively investigated (Arsenijevic et al. 1998; van Golen and Feldman 2000). As previously mentioned, such neurotrophic factors can be also produced and/or act on immune cells (Kerschenseteiner et al., 2003). Most of the physiologic effects of IGF-I are mediated by the type I IGF surface receptor which, after ligand binding-induced autophosphorylation, associates with specific adaptor proteins and activates different second messengers (Bondy and Cheng 2004). These include: phosphatidylinositol 3-kinase, mitogen-activated protein kinase (Vincent and Feldman 2002; Di Toro et al. 2005) and members of the Janus kinase (JAK)/STAT3 signalling pathway (Zong et al. 2000; Yadav et al. 2005). REST plays a complex role in neuronal cells by differentially repressing target gene expression (Lunyak et al. 2004; Coulson 2005; Ballas and Mandel 2005). REST expression decreases during neurogenesis, but has been detected in the adult rat brain (Palm et al. 1998) and is up-regulated in response to global ischemia (Calderone et al. 2003) and induction of epilepsy (Spencer et al. 2006). Thus, the REST concentration seems to influence its function and the expression of neuronal genes, and may have different effects in embryonic and differentiated neurons (Su et al. 2004; Sun et al. 2005). In a previous study, REST was elevated during the early stages of neural induction by IGF-I in neuroblastoma cells. REST may contribute to the down-regulation of genes not yet required by the differentiation program, but its expression decreases after five days of treatment to allow for the acquisition of neural phenotypes. Di Toro et al. proposed a model in which the extent of neurite outgrowth in differentiating neuroblastoma cells was affected by the disappearance of REST (Di Toro et al. 2005). The human mu-opioid receptor gene (OPRM1) promoter contains a DNA sequence binding the repressor element 1 silencing transcription factor (REST) that is implicated in transcriptional repression. Therefore, in the fist part of this thesis, I investigated whether insulin-like growth factor I (IGF-I), which affects various aspects of neuronal induction and maturation, regulates OPRM1 transcription in neuronal cells in the context of the potential influence of REST. A series of OPRM1-luciferase promoter/reporter constructs were transfected into two neuronal cell models, neuroblastoma-derived SH-SY5Y cells and PC12 cells. In the former, endogenous levels of human mu-opioid receptor (hMOPr) mRNA were evaluated by real-time PCR. IGF-I upregulated OPRM1 transcription in: PC12 cells lacking REST, in SH-SY5Y cells transfected with constructs deficient in the REST DNA binding element, or when REST was down-regulated in retinoic acid-differentiated cells. IGF-I activates the signal transducer and activator of transcription-3 (STAT3) signaling pathway and this transcription factor, binding to the STAT1/3 DNA element located in the promoter, increases OPRM1 transcription. T-cell receptor (TCR) recognizes peptide antigens displayed in the context of the major histocompatibility complex (MHC) and gives rise to a potent as well as branched intracellular signalling that convert naïve T-cells in mature effectors, thus significantly contributing to the genesis of a specific immune response. In the second part of my work I exposed wild type Jurkat CD4+ T-cells to a mixture of CD3 and CD28 antigens in order to fully activate TCR and study whether its signalling influence OPRM1 expression. Results were that TCR engagement determined a significant induction of OPRM1 expression through the activation of transcription factors AP-1, NF-kB and NFAT. Eventually, I investigated MOPr turnover once it has been expressed on T-cells outer membrane. It turned out that DAMGO induced MOPr internalisation and recycling, whereas morphine did not. Overall, from the data collected in this thesis we can conclude that that a reduction in REST is a critical switch enabling IGF-I to up-regulate human MOPr, helping these findings clarify how human MOPr expression is regulated in neuronal cells, and that TCR engagement up-regulates OPRM1 transcription in T-cells. My results that neurotrophic factors a and TCR engagement, as well as it is reported for cytokines, seem to up-regulate OPRM1 in both neurons and immune cells suggest an important role for MOPr as a molecular bridge between neurons and immune cells; therefore, MOPr could play a key role in the cross-talk between immune system and nervous system and in particular in the balance between pro-inflammatory and pro-nociceptive stimuli and analgesic and neuroprotective effects.

Morphological and functional study of cerebral aging in rat

Aging is a physiological process characterized by a progressive decline of the “cellular homeostatic reserve”, refereed as the capability to respond suitably to exogenous and endogenous stressful stimuli. Due to their high energetic requests and post-mitotic nature, neurons are peculiarly susceptible to this phenomenon. However, the aged brain maintains a certain level of adaptive capacities and if properly stimulated may warrant a considerable functional recovery. Aim of the present research was to verify the plastic potentialities of the aging brain of rats subjected to two kind of exogenous stimuli: A) the replacement of the standard diet with a ketogenic regimen (the change forces the brain to use ketone bodies (KB) in alternative to glucose to satisfy the energetic needs) and B) a behavioural task able to induce the formation of inhibitory avoidance memory. A) Fifteen male Wistar rats of 19 months of age were divided into three groups (average body weight pair-matched), and fed for 8 weeks with different dietary regimens: i) diet containing 10% medium chain triglycerides (MCT); ii) diet containing 20% MCT; iii) standard commercial chow. Five young (5 months of age) and five old (26-27 months of age) animals fed with the standard diet were used as further controls. The following morphological parameters reflecting synaptic plasticity were evaluated in the stratum moleculare of the hippocampal CA1 region (SM CA1), in the outer molecular layer of the hippocampal dentate gyrus (OML DG), and in the granule cell layer of the cerebellar cortex (GCL-CCx): average area (S), numeric density (Nvs), and surface density (Sv) of synapses, and average volume (V), numeric density (Nvm), and volume density (Vv) of synaptic mitochondria. Moreover, succinic dehydrogenase (SDH) activity was cytochemically determined in Purkinje cells (PC) and V, Nvm, Vv, and cytochemical precipitate area/mitochondrial area (R) of SDH-positive mitochondria were evaluated. In SM CA1, MCT-KDs induced the early appearance of the morphological patterns typical of old animals: higher S and V, and lower Nvs and Nvm. On the contrary, in OML DG, Sv and Vv of MCT-KDs-fed rats were higher (as a result of higher Nvs and Nvm) vs. controls; these modifications are known to improve synaptic function and metabolic supply. The opposite effects of MCT-KDs might reflect the different susceptibility of these brain regions to the aging processes: OML DG is less vulnerable than SM CA1, and the reactivation of ketone bodies uptake and catabolism might occur more efficiently in this region, allowing the exploitation of their peculiar metabolic properties. In GCL-CCx, the results described a new scenario in comparison to that found in the hippocampal formation: 10%MCT-KD induced the early appearance of senescent patterns (decreased Nvs and Nvm; increased V), whereas 20%MCT-KD caused no changes. Since GCL-CCx is more vulnerable to age than DG, and less than CA1, these data further support the hypothesis that MCT-KDs effects in the aging brain critically depend on neuronal vulnerability to age, besides MCT percentage. Regarding PC, it was decided to evaluate only the metabolic effect of the dietetic regimen (20%MCT-KD) characterized by less side effects. KD counteracted age-related decrease in numeric density of SDH-positive mitochondria, and enhanced their energetic efficiency (R was significantly higher in MCT-KD-fed rats vs. all the controls). Since it is well known that Purkinje and dentate gyrus cells are less vulnerable to aging than CA1 neurons, these results corroborate our previous hypothesis. In conclusion, the A) experimental line provides the first evidence that morphological and functional parameters reflecting synaptic plasticity and mitochondrial metabolic competence may be modulated by MCT-KDs in the pre-senescent central nervous system, and that the effects may be heterogeneous in different brain regions. MCT-KDs seem to supply high energy metabolic intermediates and to be beneficial (“anti-aging”) for those neurons that maintain the capability to exploit them. This implies risks but also promising potentialities for the therapeutic use of these diets during aging B) Morphological parameters of synapses and synaptic mitochondria in SM CA1 were investigated in old (26-27 month-old) female Wistar rats following a single trial inhibitory avoidance task. In this memory protocol animals learn to avoid a dark compartment in which they received a mild, inescapable foot-shock. Rats were tested 3 and 6 or 9 hours after the training, divided into good and bad responders according to their performance (retention times above or below 100 s, respectively) and immediately sacrificed. Nvs, S, Sv, Nvm, V, and Vv were evaluated. In the good responder group, the numeric density of synapses and mitochondria was significantly higher and the average mitochondrial volume was significantly smaller 9 hours vs. 6 hours after the training. No significant differences were observed among bad responders. Thus, better performances in passive avoidance memory task are correlated with more efficient plastic remodeling of synaptic contacts and mitochondria in hippocampal CA1. These findings indicate that maintenance of synaptic plastic reactivity during aging is a critical requirement for preserving long-term memory consolidation.

Supramolecular functional assemblies with photo or chemical addressability

Supramolecular self-assembly represents a key technology for the spontaneous construction of nanoarchitectures and for the fabrication of materials with enhanced physical and chemical properties. In addition, a significant asset of supramolecular self-assemblies rests on their reversible formation, thanks to the kinetic lability of their non-covalent interactions. This dynamic nature can be exploited for the development of “self-healing” and “smart” materials towards the tuning of their functional properties upon various external factors. One particular intriguing objective in the field is to reach a high level of control over the shape and size of the supramolecular architectures, in order to produce well-defined functional nanostructures by rational design. In this direction, many investigations have been pursued toward the construction of self-assembled objects from numerous low-molecular weight scaffolds, for instance by exploiting multiple directional hydrogen-bonding interactions. In particular, nucleobases have been used as supramolecular synthons as a result of their efficiency to code for non-covalent interaction motifs. Among nucleobases, guanine represents the most versatile one, because of its different H-bond donor and acceptor sites which display self-complementary patterns of interactions. Interestingly, and depending on the environmental conditions, guanosine derivatives can form various types of structures. Most of the supramolecular architectures reported in this Thesis from guanosine derivatives require the presence of a cation which stabilizes, via dipole-ion interactions, the macrocyclic G-quartet that can, in turn, stack in columnar G-quadruplex arrangements. In addition, in absence of cations, guanosine can polymerize via hydrogen bonding to give a variety of supramolecular networks including linear ribbons. This complex supramolecular behavior confers to the guanine-guanine interactions their upper interest among all the homonucleobases studied. They have been subjected to intense investigations in various areas ranging from structural biology and medicinal chemistry – guanine-rich sequences are abundant in telomeric ends of chromosomes and promoter regions of DNA, and are capable of forming G-quartet based structures– to material science and nanotechnology. This Thesis, organized into five Chapters, describes mainly some recent advances in the form and function provided by self-assembly of guanine based systems. More generally, Chapter 4 will focus on the construction of supramolecular self-assemblies whose self-assembling process and self-assembled architectures can be controlled by light as external stimulus. Chapter 1 will describe some of the many recent studies of G-quartets in the general area of nanoscience. Natural G- quadruplexes can be useful motifs to build new structures and biomaterials such as self-assembled nanomachines, biosensors, therapeutic aptamer and catalysts. In Chapters 2-4 it is pointed out the core concept held in this PhD Thesis, i.e. the supramolecular organization of lipophilic guanosine derivatives with photo or chemical addressability. Chapter 2 will mainly focus on the use of cation-templated guanosine derivatives as a potential scaffold for designing functional materials with tailored physical properties, showing a new way to control the bottom-up realization of well-defined nanoarchitectures. In section 2.6.7, the self-assembly properties of compound 28a may be considered an example of open-shell moieties ordered by a supramolecular guanosine architecture showing a new (magnetic) property. Chapter 3 will report on ribbon-like structures, supramolecular architectures formed by guanosine derivatives that may be of interest for the fabrication of molecular nanowires within the framework of future molecular electronic applications. In section 3.4 we investigate the supramolecular polymerizations of derivatives dG 1 and G 30 by light scattering technique and TEM experiments. The obtained data reveal the presence of several levels of organization due to the hierarchical self-assembly of the guanosine units in ribbons that in turn aggregate in fibrillar or lamellar soft structures. The elucidation of these structures furnishes an explanation to the physical behaviour of guanosine units which display organogelator properties. Chapter 4 will describe photoresponsive self-assembling systems. Numerous research examples have demonstrated that the use of photochromic molecules in supramolecular self-assemblies is the most reasonable method to noninvasively manipulate their degree of aggregation and supramolecular architectures. In section 4.4 we report on the photocontrolled self-assembly of modified guanosine nucleobase E-42: by the introduction of a photoactive moiety at C8 it is possible to operate a photocontrol over the self-assembly of the molecule, where the existence of G-quartets can be alternately switched on and off. In section 4.5 we focus on the use of cyclodextrins as photoresponsive host-guest assemblies: αCD–azobenzene conjugates 47-48 (section 4.5.3) are synthesized in order to obtain a photoresponsive system exhibiting a fine photocontrollable degree of aggregation and self-assembled architecture. Finally, Chapter 5 contains the experimental protocols used for the research described in Chapters 2-4.

Functional Genomics and Cell Biology of the Dolphin (Tursiops runcatus): Establishment of Novel Molecular Tools to Study Marine Mammals in Changing Environments

The dolphin (Tursiops truncatus) is a mammal that is adapted to life in a totally aquatic environment. Despite the popularity and even iconic status of the dolphin, our knowledge of its physiology, its unique adaptations and the effects on it of environmental stressors are limited. One approach to improve this limited understanding is the implementation of established cellular and molecular methods to provide sensitive and insightful information for dolphin biology. We initiated our studies with the analysis of wild dolphin peripheral blood leukocytes, which have the potential to be informative of the animal’s global immune status. Transcriptomic profiles from almost 200 individual samples were analyzed using a newly developed species-specific microarray to assess its value as a prognostic and diagnostic tool. Functional genomics analyses were informative of stress-induced gene expression profiles and also of geographical location specific transcriptomic signatures, determined by the interaction of genetic, disease and environmental factors. We have developed quantitative metrics to unambiguously characterize the phenotypic properties of dolphin cells in culture. These quantitative metrics can provide identifiable characteristics and baseline data which will enable identification of changes in the cells due to time in culture. We have also developed a novel protocol to isolate primary cultures from cryopreserved tissue of stranded marine mammals, establishing a tissue (and cell) biorepository, a new approach that can provide a solution to the limited availability of samples. The work presented represents the development and application of tools for the study of the biology, health and physiology of the dolphin, and establishes their relevance for future studies of the impact on the dolphin of environmental infection and stress.

Functional characterization of Streptococcus pyogenes pili

Group A Streptococcus is a Gram-positive human pathogen able to colonize both upper respiratory tract and skin. GAS is responsible for several acute diseases and autoimmune sequelae that account for half a million deaths worldwide every year (Cunningham et al., 2000). As other bacteria, GAS infections requires the capacity of the pathogen to adhere to host tissues and to form cell aggregates. The ability to persist in distinct host niches like the throat and the skin and to trigger infections is associated with the expression of different GAS virulence factors. GAS pili has been described as important virulence factors encoded by different FCT-operon regions. Based on this information, we decided to study the possible effect of environmental conditions that could regulate the pili expression. In this study we reported the influence of pH environment variations in biofilm formation for strains pertaining to a panel of different GAS FCT-types. The biofilm formation was promoted, excepted in the FCT-1 strains, by a changing in pH from physiological to acidic condition of growth in in vitro biofilm assay. By analyzing the possible association between biofilm formation and pH dependence, we have found that in FCT-2 and FCT-3 strains, the biofilm is promoted by pH reduction leading to an increase of pili expression. These data confirmed a direct link between pH dependent pilus expression and biofilm formation in GAS. As pili are a multi component structure we decided to investigate the functional role of one of its subunits, the AP-1 protein. AP-1 is highly conserved through the different FCT-types and suggests a possible essential role for the pili function. We focused our attention on the AP-1 protein encoded by the FCT-1 strains (M6). In particular this AP-1 protein contains the von Willebrand Factor A (VWFA) domain, which share an homology with the human VWFA domain that has been reported to be involved in adhesion process. We have demonstrated that the AP-1 protein binds to human epithelial cells by its VWFA domain, whereas the biofilm formation is mediated by the N-terminal region of AP-1 protein. Moreover, analyzing the importance of AP-1 in in vivo experiments we found a major capacity of tissue dissemination for the wild-type strain compared to the isogenic AP-1 deletion mutant. Pili have been also reported as potential vaccine candidates against Gram positive bacteria. For these reason we decided to investigate the relationship between cross reaction of sera raised against different GAS and GBS pilin subunits and the presence of a conserved Cna_B domain, in different pilin components. Our idea was to investigate if, using pilus conserved domains, a broad coverage vaccine against streptococcal infection could be possible.

Epidemiological and functional characterization of Streptococcus pneumoniae pili

Streptococcus pneumoniae is an important life threatening human pathogen causing agent of invasive diseases such as otitis media, pneumonia, sepsis and meningitis, but is also a common inhabitant of the respiratory tract of children and healthy adults. Likewise most streptococci, S. pneumoniae decorates its surface with adhesive pili, composed of covalently linked subunits and involved in the attachment to epithelial cells and virulence. The pneumococcal pili are encoded by two genomic regions, pilus islet 1 (PI-1), and pilus islet-2 (PI-2), which are present in about 30% and 16% of the pneumococcal strains, respectively. PI-1 exists in three clonally related variants, whereas PI-2 is highly conserved. The presence of the islets does not correlate with the serotype of the strains, but with the genotype (as determined by Multi Locus Sequence Typing). The prevalence of PI-1 and PI-2 positive strains is similar in isolates from invasive disease and carriage. To better dissect a possible association between PIs presence and disease we evaluated the distribution of the two PIs in a panel of 113 acute otitis media (AOM) clinical isolates from Israel. PI-1 was present in 30.1% (N=34) of the isolates tested, and PI-2 in 7% (N=8). We found that 50% of the PI-1 positive isolates belonged to the international clones Spain9V-3 (ST156) and Taiwan19F-14 (ST236), and that PI-2 was not present in the absence of Pl-1. In conclusion, there was no correlation between PIs presence and AOM, and, in general, the observed differences in PIs prevalence are strictly dependent upon regional differences in the distribution of the clones. Finally, in the AOM collection the prevalence of PI-1 was higher among antibiotic resistant isolates, confirming previous indications obtained by the in silico analysis of the MLST database collection. Since the pilus-1 subunits were shown to confer protection in mouse models of infection both in active and passive immunization studies, and were regarded as potential candidates for a new generation of protein-based vaccines, the functional characterization was mainly focused on S. pneumoniae pilus -1 components. The pneumococcal pilus-1 is composed of three subunits, RrgA, RrgB and RrgC, each stabilized by intra-molecular isopeptide bonds and covalently polymerized by means of inter-molecular isopeptide bonds to form an extended fibre. The pilus shaft is a multimeric structure mainly composed by the RrgB backbone subunit. The minor ancillary proteins are located at the tip and at the base of the pilus, where they have been proposed to act as the major adhesin (RrgA) and as the pilus anchor (RrgC), respectively. RrgA is protective in in vivo mouse models, and exists in two variants (clades I and II). Mapping of the sequence variability onto the RrgA structure predicted from X-ray data showed that the diversity was restricted to the “head” of the protein, which contains the putative binding domains, whereas the elongated “stalk” was mostly conserved. To investigate whether this variability could influence the adhesive capacity of RrgA and to map the regions important for binding, two full-length protein variants and three recombinant RrgA portions were tested for adhesion to lung epithelial cells and to purified extracellular matrix (ECM) components. The two RrgA variants displayed similar binding abilities, whereas none of the recombinant fragments adhered at levels comparable to those of the full-length protein, suggesting that proper folding and structural arrangement are crucial to retain protein functionality. Furthermore, the two RrgA variants were shown to be cross-reactive in vitro and cross-protective in vivo in a murine model of passive immunization. Taken together, these data indicate that the region implicated in adhesion and the functional epitopes responsible for the protective ability of RrgA may be conserved and that the considerable level of variation found within the “head” domain of RrgA may have been generated by immunologic pressure without impairing the functional integrity of the pilus.

Trascriptome analysis and functional genomics of Neisseria meningitidis in human blood

Neisseria meningitidis (Nm) is the major cause of septicemia and meningococcal meningitis. During the course of infection, it must adapt to different host environments as a crucial factor for survival. Despite the severity of meningococcal sepsis, little is known about how Nm adapts to permit survival and growth in human blood. A previous time-course transcriptome analysis, using an ex vivo model of human whole blood infection, showed that Nm alters the expression of nearly 30% of ORFs of the genome: major dynamic changes were observed in the expression of transcriptional regulators, transport and binding proteins, energy metabolism, and surface-exposed virulence factors. Starting from these data, mutagenesis studies of a subset of up-regulated genes were performed and the mutants were tested for the ability to survive in human whole blood; Nm mutant strains lacking the genes encoding NMB1483, NalP, Mip, NspA, Fur, TbpB, and LctP were sensitive to killing by human blood. Then, the analysis was extended to the whole Nm transcriptome in human blood, using a customized 60-mer oligonucleotide tiling microarray. The application of specifically developed software combined with this new tiling array allowed the identification of different types of regulated transcripts: small intergenic RNAs, antisense RNAs, 5’ and 3’ untranslated regions and operons. The expression of these RNA molecules was confirmed by 5’-3’RACE protocol and specific RT-PCR. Here we describe the complete transcriptome of Nm during incubation in human blood; we were able to identify new proteins important for survival in human blood and also to identify additional roles of previously known virulence factors in aiding survival in blood. In addition the tiling array analysis demonstrated that Nm expresses a set of new transcripts, not previously identified, and suggests the presence of a circuit of regulatory RNA elements used by Nm to adapt to proliferate in human blood.

Lesion detection and classification in breast cancer: evaluation of approaches based on morphological features, tracer kinetic modelling and semi-quantitative parameters in MR functional imaging (DCE-MRI)

The diagnosis, grading and classification of tumours has benefited considerably from the development of DCE-MRI which is now essential to the adequate clinical management of many tumour types due to its capability in detecting active angiogenesis. Several strategies have been proposed for DCE-MRI evaluation. Visual inspection of contrast agent concentration curves vs time is a very simple yet operator dependent procedure, therefore more objective approaches have been developed in order to facilitate comparison between studies. In so called model free approaches, descriptive or heuristic information extracted from time series raw data have been used for tissue classification. The main issue concerning these schemes is that they have not a direct interpretation in terms of physiological properties of the tissues. On the other hand, model based investigations typically involve compartmental tracer kinetic modelling and pixel-by-pixel estimation of kinetic parameters via non-linear regression applied on region of interests opportunely selected by the physician. This approach has the advantage to provide parameters directly related to the pathophysiological properties of the tissue such as vessel permeability, local regional blood flow, extraction fraction, concentration gradient between plasma and extravascular-extracellular space. Anyway, nonlinear modelling is computational demanding and the accuracy of the estimates can be affected by the signal-to-noise ratio and by the initial solutions. The principal aim of this thesis is investigate the use of semi-quantitative and quantitative parameters for segmentation and classification of breast lesion. The objectives can be subdivided as follow: describe the principal techniques to evaluate time intensity curve in DCE-MRI with focus on kinetic model proposed in literature; to evaluate the influence in parametrization choice for a classic bi-compartmental kinetic models; to evaluate the performance of a method for simultaneous tracer kinetic modelling and pixel classification; to evaluate performance of machine learning techniques training for segmentation and classification of breast lesion.

A clustering method for robust and reliable large scale functional and structural protein sequence annotation

Bioinformatics, in the last few decades, has played a fundamental role to give sense to the huge amount of data produced. Obtained the complete sequence of a genome, the major problem of knowing as much as possible of its coding regions, is crucial. Protein sequence annotation is challenging and, due to the size of the problem, only computational approaches can provide a feasible solution. As it has been recently pointed out by the Critical Assessment of Function Annotations (CAFA), most accurate methods are those based on the transfer-by-homology approach and the most incisive contribution is given by cross-genome comparisons. In the present thesis it is described a non-hierarchical sequence clustering method for protein automatic large-scale annotation, called “The Bologna Annotation Resource Plus” (BAR+). The method is based on an all-against-all alignment of more than 13 millions protein sequences characterized by a very stringent metric. BAR+ can safely transfer functional features (Gene Ontology and Pfam terms) inside clusters by means of a statistical validation, even in the case of multi-domain proteins. Within BAR+ clusters it is also possible to transfer the three dimensional structure (when a template is available). This is possible by the way of cluster-specific HMM profiles that can be used to calculate reliable template-to-target alignments even in the case of distantly related proteins (sequence identity < 30%). Other BAR+ based applications have been developed during my doctorate including the prediction of Magnesium binding sites in human proteins, the ABC transporters superfamily classification and the functional prediction (GO terms) of the CAFA targets. Remarkably, in the CAFA assessment, BAR+ placed among the ten most accurate methods. At present, as a web server for the functional and structural protein sequence annotation, BAR+ is freely available at http://bar.biocomp.unibo.it/bar2.0.

Functional and Genetic characterization of new genomic islands from an E. coli strain associated with neonatal meningitis.

Enterobacteriaceae genomes evolve through mutations, rearrangements and horizontal gene transfer (HGT). The latter evolutionary pathway works through the acquisition DNA (GEI) modules of foreign origin that enhances fitness of the host to a given environment. The genome of E. coli IHE3034, a strain isolated from a case of neonatal meningitis, has recently been sequenced and its subsequent sequence analysis has predicted 18 possible GEIs, of which: 8 have not been previously described, 5 fully meet the pathogenic island definition and at least 10 that seem to be of prophagic origin. In order to study the GEI distribution of our reference strain, we screened for the presence 18 GEIs a panel of 132 strains, representative of E. coli diversity. Also, using an inverse nested PCR approach we identified 9 GEI that can form an extrachromosomal circular intermediate (CI) and their respective attachment sites (att). Further, we set up a qPCR approach that allowed us to determine the excision rates of 5 genomic islands in different growth conditions. Four islands, specific for strains appertaining to the sequence type complex 95 (STC95), have been deleted in order to assess their function in a Dictyostelium discoideum grazing assays. Overall, the distribution data presented here indicate that 16 IHE3034 GEIs are more associated to the STC95 strains. Also the functional and genetic characterization has uncovered that GEI 13, 17 and 19 are involved in the resistance to phagocitation by Dictyostelium d thus suggesting a possible role in the adaptation of the pathogen during certain stages of infection.

Identification and functional characterization of human cytomegalovirus Fc binding proteins

ABSTRACT Human cytomegalovirus (HCMV) employs many different mechanisms to escape and subvert the host immune system surveillance. Among these different mechanisms the role of human IgG Fc receptors (FcγR) in HCMV pathogenesis is still unclear. In mammalians, FcγRs are expressed on the surface of all haematopoietic cells and have a multifaceted role in regulating the activity of antibodies to generate a well-balanced immune response. Viral proteins with Fcγ binding ability are highly diffuse among herpesviruses. They interfere with the host receptors functions in order to counteract immune system recognition. So far, two human HCMV Fcγ binding proteins have been described: UL119 and RL11. This work was aimed to the identification and characterization of HCMV Fcγ binding proteins. The study is divided in two parts: first the characterization of UL119 and RL11; second the identification and characterization of novel HCMV Fcγ binding proteins. Regarding the first part, we demonstrated that both UL119 and RL11 internalize Fcγ fragments from transfected cells surface through a clathrin dependent pathway. In infected cells both proteins were found in the viral assembly complex and on virions surface as envelope associated glycoproteins. Moreover, internalized Fcγ in infected cells do not undergo lysosomal degradation but rather traffic in early endosomes up to the viral assembly complex. Regarding the second part, we were able to identify two novels Fcγ binding protein coded by CMV: RL12 and RL13. The latter was also further characterized as recombinant protein in terms of cellular localization, Fc binding site and IgG internalization ability. Finally binding specificity of both RL12 and RL13 seems to be confined to human IgG1 and IgG2. Taken together, these data show that HCMV codes for up to 4 FcγR and that they could have a double role both on virus and on infected cells.