Procedures of quality control and data analysis of multi-site ground-based observations for the absolute flux calibration of Gaia

The Gaia space mission is a major project for the European astronomical community. As challenging as it is, the processing and analysis of the huge data-flow incoming from Gaia is the subject of thorough study and preparatory work by the DPAC (Data Processing and Analysis Consortium), in charge of all aspects of the Gaia data reduction. This PhD Thesis was carried out in the framework of the DPAC, within the team based in Bologna. The task of the Bologna team is to define the calibration model and to build a grid of spectro-photometric standard stars (SPSS) suitable for the absolute flux calibration of the Gaia G-band photometry and the BP/RP spectrophotometry. Such a flux calibration can be performed by repeatedly observing each SPSS during the life-time of the Gaia mission and by comparing the observed Gaia spectra to the spectra obtained by our ground-based observations. Due to both the different observing sites involved and the huge amount of frames expected (≃100000), it is essential to maintain the maximum homogeneity in data quality, acquisition and treatment, and a particular care has to be used to test the capabilities of each telescope/instrument combination (through the “instrument familiarization plan”), to devise methods to keep under control, and eventually to correct for, the typical instrumental effects that can affect the high precision required for the Gaia SPSS grid (a few % with respect to Vega). I contributed to the ground-based survey of Gaia SPSS in many respects: with the observations, the instrument familiarization plan, the data reduction and analysis activities (both photometry and spectroscopy), and to the maintenance of the data archives. However, the field I was personally responsible for was photometry and in particular relative photometry for the production of short-term light curves. In this context I defined and tested a semi-automated pipeline which allows for the pre-reduction of imaging SPSS data and the production of aperture photometry catalogues ready to be used for further analysis. A series of semi-automated quality control criteria are included in the pipeline at various levels, from pre-reduction, to aperture photometry, to light curves production and analysis.

Influence of nitrogen and soil physical characteristics on belowground carbon flux dynamics of woody plants

At ecosystem level soil respiration (Rs) represents the largest carbon (C) flux after gross primary productivity, being mainly generated by root respiration (autotrophic respiration, Ra) and soil microbial respiration (heterotrophic respiration, Rh). In the case of terrestrial ecosystems, soils contain the largest C-pool, storing twice the amount of C contained in plant biomass. Soil organic matter (SOM), representing the main C storage in soil, is decomposed by soil microbial community. This process produces CO2 which is mainly released as Rh. It is thus relevant to understand how microbial activity is influenced by environmental factors like soil temperature, soil moisture and nutrient availability, since part of the CO2 produced by Rh, directly increases atmospheric CO2 concentration and therefore affects the phenomenon of climate change. Among terrestrial ecosystems, agricultural fields have traditionally been considered as sources of atmospheric CO2. In agricultural ecosystems, in particular apple orchards, I identified the role of root density, soil temperature, soil moisture and nitrogen (N) availability on Rs and on its two components, Ra and Rh. To do so I applied different techniques to separate Rs in its two components, the ”regression technique” and the “trenching technique”. I also studied the response of Ra to different levels of N availability, distributed either in a uniform or localized way, in the case of Populus tremuloides trees. The results showed that Rs is mainly driven by soil temperature, to which it is positively correlated, that high levels of soil moisture have inhibiting effects, and that N has a negligible influence on total Rs, as well as on Ra. Further I found a negative response of Rh to high N availability, suggesting that microbial decomposition processes in the soil are inhibited by the presence of N. The contribution of Ra to Rs was of 37% on average.

Characterization of the Staphylococcus aureus bone sialoprotein-binding protein SdrE and the serine protease EpiP

In an attempt to develop a Staphylococcus aureus vaccine, we have applied reverse vaccinology approach, mainly based on in silico screening and proteomics. By using this approach SdrE, a protein belonging to serine-aspartate repeat protein family was identified as potential vaccine antigen against S. aureus. We have investigated the biochemical properties as well as the vaccine potential of SdrE and its highly conserved CnaBE3 domain. We found the protein SdrE to be resistant to trypsin. Further analysis of the resistant fragment revealed that it comprises a CnaBE3 domain, which also showed partial trypsin resistant behavior. Furthermore, intact mass spectrometry of rCnaBE3 suggested the possible presence of isopeptide bond or some other post-translational modification in the protein.However, this observation needs further investigation. Differential Scanning Fluorimetry study reveals that calcium play role in protein folding and provides stability to SdrE. At the end we have demonstrated that SdrE is immunogenic against clinical strain of S. aureus in murine abscess model. In the second part, I characterized a protein, annotated as epidermin leader peptide processing serine protease (EpiP), as a novel S. aureus vaccine candidate. The crystal structure of the rEpiP was solved at 2.05 Å resolution by x-ray crystallography . The structure showed that rEpiP was cleaved somewhere between residues 95 and 100 and cleavage occurs through an autocatalytic intra-molecular mechanism. In addition, the protein expressed by S. aureus cells also appeared to undergo a similar processing event. To determine if the protein acts as a serine protease, we mutated the catalytic serine 393 residue to alanine, generating rEpiP-S393A and solved its crystal structure at a resolution of 1.95 Å. rEpiP-S393A was impaired in its protease activity, as expected. Protective efficacy of rEpiP and the non-cleaving mutant protein was comparable, implying that the two forms are interchangeable for vaccination purposes.

Analysis of ion channel stochastic signals for biosensing

In biological world, life of cells is guaranteed by their ability to sense and to respond to a large variety of internal and external stimuli. In particular, excitable cells, like muscle or nerve cells, produce quick depolarizations in response to electrical, mechanical or chemical stimuli: this means that they can change their internal potential through a quick exchange of ions between cytoplasm and the external environment. This can be done thanks to the presence of ion channels, proteins that span the lipid bilayer and act like switches, allowing ionic current to flow opening and shutting in a stochastic way. For a particular class of ion channels, ligand-gated ion channels, the gating processes is strongly influenced by binding between receptive sites located on the channel surface and specific target molecules. These channels, inserted in biomimetic membranes and in presence of a proper electronic system for acquiring and elaborating the electrical signal, could give us the possibility of detecting and quantifying concentrations of specific molecules in complex mixtures from ionic currents across the membrane; in this thesis work, this possibility is investigated. In particular, it reports a description of experiments focused on the creation and the characterization of artificial lipid membranes, the reconstitution of ion channels and the analysis of their electrical and statistical properties. Moreover, after a chapter about the basis of the modelling of the kinetic behaviour of ligand gated ion channels, a possible approach for the estimation of the target molecule concentration, based on a statistical analysis of the ion channel open probability, is proposed. The fifth chapter contains a description of the kinetic characterisation of a ligand gated ion channel: the homomeric α2 isoform of the glycine receptor. It involved both experimental acquisitions and signal analysis. The last chapter represents the conclusions of this thesis, with some remark on the effective performance that may be achieved using ligand gated ion channels as sensing elements.

Identification of Drosophila heart-specific Cis-Regulatory Modules under Hox control

Cardiac morphogenesis is a complex process governed by evolutionarily conserved transcription factors and signaling molecules. The Drosophila cardiac tube is linear, made of 52 pairs of cardiomyocytes (CMs), which express specific transcription factor genes that have human homologues implicated in Congenital Heart Diseases (CHDs) (NKX2-5, GATA4 and TBX5). The Drosophila cardiac tube is linear and composed of a rostral portion named aorta and a caudal one called heart, distinguished by morphological and functional differences controlled by Hox genes, key regulators of axial patterning. Overexpression and inactivation of the Hox gene abdominal-A (abd-A), which is expressed exclusively in the heart, revealed that abd-A controls heart identity. The aim of our work is to isolate the heart-specific cisregulatory sequences of abd-A direct target genes, the realizator genes granting heart identity. In each segment of the heart, four pairs of cardiomyocytes (CMs) express tinman (tin), homologous to NKX2-5, and acquire strong contractile and automatic rhythmic activities. By tyramide amplified FISH, we found that seven genes, encoding ion channels, pumps or transporters, are specifically expressed in the Tin-CMs of the heart. We initially used online available tools to identify their heart-specific cisregutatory modules by looking for Conserved Non-coding Sequences containing clusters of binding sites for various cardiac transcription factors, including Hox proteins. Based on these data we generated several reporter gene constructs and transgenic embryos, but none of them showed reporter gene expression in the heart. In order to identify additional abd-A target genes, we performed microarray experiments comparing the transcriptomes of aorta versus heart and identified 144 genes overexpressed in the heart. In order to find the heart-specific cis-regulatory regions of these target genes we developed a new bioinformatic approach where prediction is based on pattern matching and ordered statistics. We first retrieved Conserved Noncoding Sequences from the alignment between the D.melanogaster and D.pseudobscura genomes. We scored for combinations of conserved occurrences of ABD-A, ABD-B, TIN, PNR, dMEF2, MADS box, T-box and E-box sites and we ranked these results based on two independent strategies. On one hand we ranked the putative cis-regulatory sequences according to best scored ABD-A biding sites, on the other hand we scored according to conservation of binding sites. We integrated and ranked again the two lists obtained independently to produce a final rank. We generated nGFP reporter construct flies for in vivo validation. We identified three 1kblong heart-specific enhancers. By in vivo and in vitro experiments we are determining whether they are direct abd-A targets, demonstrating the role of a Hox gene in the realization of heart identity. The identified abd-A direct target genes may be targets also of the NKX2-5, GATA4 and/or TBX5 homologues tin, pannier and Doc genes, respectively. The identification of sequences coregulated by a Hox protein and the homologues of transcription factors causing CHDs, will provide a mean to test whether these factors function as Hox cofactors granting cardiac specificity to Hox proteins, increasing our knowledge on the molecular mechanisms underlying CHDs. Finally, it may be investigated whether these Hox targets are involved in CHDs.

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.

Synthesis and application of new ion-tagged ligands and organocatalysts

We report the synthesis and application of some ion-tagged catalysts in organometallic catalysis and organocatalysis. With the installation of an ionic group on the backbone of a known catalyst, two main effects are generally obtained. i) a modification of the solubility of the catalyst: if judicious choice of the ion pair is made, the ion-tag can confer to the catalyst a solubility profile suitable for catalyst recycling. ii) the ionic group can play a non-innocent role in the process considered: if stabilizing interaction between the ionic group and the developing charges in the transition state are established, the reaction can speed up. We describe the use of ion-tagged diphenylprolinol as Zn ligand. The chiral ligand grafted onto an ionic liquid (IL) was recycled 10 times with no loss of reactivity and selectivity, when it was employed in the first example of enantioselective addition of ZnEt2 to aldehydes in ILs. An ammonium-tagged phosphine displayed the capability to stabilize Pd catalysts for the Suzuki reaction in ILs. The ionic phase was recycled 6 times with no detectable loss of activity and very low Pd leaching in the organic phase. This catalytic system was also employed for the functionalization of the challenging substrate 5,11-dibromotetracene. In the field of organocatalysis, we prepared two ion-tagged derivatives of the McMillan imidazolidinone. The results of the asymmetric Diels-Alder reaction between trans-cinnamaldehyde and cyclopentadiene exhibited great dependence on the position and nature of the ionic group. Finally, when O-TMS-diphenylprolinol was tagged with an imidazolium ion, exploiting a silyl ether linker, an efficient catalyst for the asymmetric addition of aldehydes to nitroolefins was achieved. The catalyst displayed enhanced reactivity and the same high level of selectivity of the untagged parent catalyst and it could be employed in a wide range of reaction conditions, included use of water as solvent.