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.

Identificazione delle Lamine di tipo A come nuovi substrati della protein chinasi Akt/PKB

Akt (also called PKB) is a 63 kDa serine/threonine kinase involved in promotion of cell survival, proliferation a nd metabolic responses downstream the phosphoinositide-3-kinase (PI 3-kinase) signaling pathway. In resting cells, Akt is a predominantly cytosolic enzyme; however generation of PI 3-kinase lipid products recruits Akt to the plasma membrane, resulting in a conformational change which confers full enzymatic activity through the phosphorylation of the membrane-bound protein at two residues, Thr308, and Ser473. Activated Akt redistributes to cytoplasm and nucleus, where phosphorylation of specific substrates occurs. Both the presence and the activity of Akt in the nucleus have been described. An interesting mechanism that mediates nuclear translocation of Akt has been described in human mature T-cell leukemia: the product of TCL1 gene, Tcl1, interacts with the PH domain of phosphorylated Akt, thus driving Akt to the nucleus. In this context, Tcl1 may act as a direct transporter of Akt or may contribute to the formation of a complex that promotes the transport of active Akt to the nucleus, where it can phosphorylate nuclear substrates. A well described nuclear substrate if Foxo. IGF-1 triggers phosphorylation of Foxo by Akt inside the nucleus, where phospho-Foxo associates to 14.3.3 proteins that, in turn, promote its export to the cytoplasm where it is sequestered. Remarkably, Foxo phosphorylation by Akt has been shown to be a crucial event in Akt-dependent myogenesis. However, most Akt nuclear substrates have so far remained elusive, as well as nuclear Akt functions. This lack of information prompted us to undertake a search of substrates of Akt in the nucleus, by the combined use of 2D-separation/mass spectrometry and anti-Akt-phosphosubstrate antibody. This study presents evidence of A-type lamins as novel nuclear substrates of Akt. Lamins are type V intermediate filaments proteins found in the nucleus of higher eukaryotes where, together with lamin-binding proteins, they form the lamina at the nuclear envelope, providing mechanical stability for the nuclear membrane. By coimmunoprecipitation, it is demonstrated here that endogenous lamin A and Akt interact, and that A-type lamins are phosphorylated by Akt both in vitro and in vivo. Moreover, by phosphoaminoacid analysis and mutagenesis, it is further demonstrated that Akt phosphorylates lamin A at Ser404, and, more importantly, that while lamin A/C phosphorylation is stable throughout the cell cycle, phosphorylation of the precursor prelamin A becomes detectable as cells enter the G2 phase, picking at G2/M. This study also shows that lamin phosphorylation by Akt creates a binding site for 14.3.3 adaptors which, in turn, promote prelamin A degradation. While this mechanism is in agreement with a general role of Akt in the regulation of a subset of its substrates, opposite to what has been described, degradation is not mediated through a ubiquitination and proteasomal mechanism but through a lysosomal pathway, as indicated by the reverting action of the lysosomal inhibitor cloroquine. Phosphorylation is a key event in the mitotic breakdown of the nuclear lamina. However, the kinases and the precise sites of phosphorylation are scarcely known. Therefore, these results represent an important breakthrough in this very significant but understudied area. The phosphorylation of the precursor protein prelamin A and its subsequent degradation at G2/M, when both the nuclear envelop and the nuclear lamina disassemble, can be view as part of a mechanism to dispose off the precursor that is not needed in this precise context. The recently reported finding that patients affected by Emery-Dreifuss muscular dystrophy carry a mutation at Arg 401, in the Akt phosphorylation motif, open new perspective that warrant further investigation in this very important field.

Computational methods for the analysis of protein structure and function

The vast majority of known proteins have not yet been experimentally characterized and little is known about their function. The design and implementation of computational tools can provide insight into the function of proteins based on their sequence, their structure, their evolutionary history and their association with other proteins. Knowledge of the three-dimensional (3D) structure of a protein can lead to a deep understanding of its mode of action and interaction, but currently the structures of <1% of sequences have been experimentally solved. For this reason, it became urgent to develop new methods that are able to computationally extract relevant information from protein sequence and structure. The starting point of my work has been the study of the properties of contacts between protein residues, since they constrain protein folding and characterize different protein structures. Prediction of residue contacts in proteins is an interesting problem whose solution may be useful in protein folding recognition and de novo design. The prediction of these contacts requires the study of the protein inter-residue distances related to the specific type of amino acid pair that are encoded in the so-called contact map. An interesting new way of analyzing those structures came out when network studies were introduced, with pivotal papers demonstrating that protein contact networks also exhibit small-world behavior. In order to highlight constraints for the prediction of protein contact maps and for applications in the field of protein structure prediction and/or reconstruction from experimentally determined contact maps, I studied to which extent the characteristic path length and clustering coefficient of the protein contacts network are values that reveal characteristic features of protein contact maps. Provided that residue contacts are known for a protein sequence, the major features of its 3D structure could be deduced by combining this knowledge with correctly predicted motifs of secondary structure. In the second part of my work I focused on a particular protein structural motif, the coiled-coil, known to mediate a variety of fundamental biological interactions. Coiled-coils are found in a variety of structural forms and in a wide range of proteins including, for example, small units such as leucine zippers that drive the dimerization of many transcription factors or more complex structures such as the family of viral proteins responsible for virus-host membrane fusion. The coiled-coil structural motif is estimated to account for 5-10% of the protein sequences in the various genomes. Given their biological importance, in my work I introduced a Hidden Markov Model (HMM) that exploits the evolutionary information derived from multiple sequence alignments, to predict coiled-coil regions and to discriminate coiled-coil sequences. The results indicate that the new HMM outperforms all the existing programs and can be adopted for the coiled-coil prediction and for large-scale genome annotation. Genome annotation is a key issue in modern computational biology, being the starting point towards the understanding of the complex processes involved in biological networks. The rapid growth in the number of protein sequences and structures available poses new fundamental problems that still deserve an interpretation. Nevertheless, these data are at the basis of the design of new strategies for tackling problems such as the prediction of protein structure and function. Experimental determination of the functions of all these proteins would be a hugely time-consuming and costly task and, in most instances, has not been carried out. As an example, currently, approximately only 20% of annotated proteins in the Homo sapiens genome have been experimentally characterized. A commonly adopted procedure for annotating protein sequences relies on the "inheritance through homology" based on the notion that similar sequences share similar functions and structures. This procedure consists in the assignment of sequences to a specific group of functionally related sequences which had been grouped through clustering techniques. The clustering procedure is based on suitable similarity rules, since predicting protein structure and function from sequence largely depends on the value of sequence identity. However, additional levels of complexity are due to multi-domain proteins, to proteins that share common domains but that do not necessarily share the same function, to the finding that different combinations of shared domains can lead to different biological roles. In the last part of this study I developed and validate a system that contributes to sequence annotation by taking advantage of a validated transfer through inheritance procedure of the molecular functions and of the structural templates. After a cross-genome comparison with the BLAST program, clusters were built on the basis of two stringent constraints on sequence identity and coverage of the alignment. The adopted measure explicity answers to the problem of multi-domain proteins annotation and allows a fine grain division of the whole set of proteomes used, that ensures cluster homogeneity in terms of sequence length. A high level of coverage of structure templates on the length of protein sequences within clusters ensures that multi-domain proteins when present can be templates for sequences of similar length. This annotation procedure includes the possibility of reliably transferring statistically validated functions and structures to sequences considering information available in the present data bases of molecular functions and structures.

B-type cytochromes of the DOMON domain superfamily (Novel redox elements of plant plasma membrane)

The aim of the present study is understanding the properties of a new group of redox proteins having in common a DOMON-type domain with characteristics of cytochromes b. The superfamily of proteins containing a DOMON of this type includes a few protein families. With the aim of better characterizing this new protein family, the present work addresses both a CyDOM protein (a cytochrome b561) and a protein only comprised of DOMON(AIR12), both of plant origin. Apoplastic ascorbate can be regenerated from monodehydroascorbate by a trans-plasma membrane redox system which uses cytosolic ascorbate as a reductant and comprises a high potential cytochrome b. We identified the major plasma membrane (PM) ascorbate-reducible b-type cytochrome of bean (Phaseolus vulgaris) and soybean (Glycine max) hypocotyls as orthologs of Arabidopsis auxin-responsive gene air12. The protein, which is glycosylated and glycosylphosphatidylinositol-anchored to the external side of the PM in vivo, was expressed in Pichia pastoris in a recombinant form, lacking the glycosylphosphatidylinositol-modification signal, and purified from the culture medium. Recombinant AIR12 is a soluble protein predicted to fold into a β-sandwich domain and belonging to the DOMON superfamily. It is shown to be a b-type cytochrome with a symmetrical α-band at 561 nm, to be fully reduced by ascorbate and fully oxidized by monodehydroascorbate. Redox potentiometry suggests that AIR12 binds two high-potential hemes (Em,7 +135 and +236 mV). Phylogenetic analyses reveal that the auxin-responsive genes AIR12 constitute a new family of plasma membrane b-type cytochromes specific to flowering plants. Although AIR12 is one of the few redox proteins of the PM characterized to date, the role of AIR12 in trans-PM electron transfer would imply interaction with other partners which are still to be identified. Another part of the present project was aimed at understanding of a soybean protein comprised of a DOMON fused with a well-defined b561 cytochrome domain (CyDOM). Various bioinformatic approaches show this protein to be composed of an N-terminal DOMON followed by b561 domain. The latter contains five transmembrane helices featuring highly conserved histidines, which might bind haem groups. The CyDOM has been cloned and expressed in the yeast Pichia pastoris, and spectroscopic analyses have been accomplished on solubilized yeast membranes. CyDOM clearly reveal the properties of b-type cytochrome. The results highlight the fact that CyDOM is clearly able to lead an electron flux through the plasmamembrane. Voltage clamp experiments demonstrate that Xenopus laevis oocytes transformed with CyDOM of soybean exhibit negative electrical currents in presence of an external electron acceptor. Analogous investigations were carried out with SDR2, a CyDOM of Drosophila melanogaster which shows an electron transport capacity even higher than plant CyDOM. As quoted above, these data reinforce those obtained in plant CyDOM on the one hand, and on the other hand allow to attribute to SDR2-like proteins the properties assigned to CyDOM. Was expressed in Regenerated tobacco roots, transiently transformed with infected a with chimeral construct GFP: CyDOM (by A. rhizogenes infection) reveals a plasmamembrane localization of CyDOM both in epidermal cells of the elongation zone of roots and in root hairs. In conclusion. Although the data presented here await to be expanded and in part clarified, it is safe to say they open a new perspective about the role of this group of proteins. The biological relevance of the functional and physiological implications of DOMON redox domains seems noteworthy, and it can but increase with future advances in research. Beyond the very finding, however interesting in itself, of DOMON domains as extracellular cytochromes, the present study testifies to the fact that cytochrome proteins containing DOMON domains of the type of “CyDOM” can transfer electrons through membranes and may represent the most important redox component of the plasmamembrane as yet discovered.

Studies of OPA1 pathogenic mechanisms in Dominant Optic Atrophy and novel protein function in mitochondrial DNA stability

The mitochondrion is an essential cytoplasmic organelle that provides most of the energy necessary for eukaryotic cell physiology. Mitochondrial structure and functions are maintained by proteins of both mitochondrial and nuclear origin. These organelles are organized in an extended network that dynamically fuses and divides. Mitochondrial morphology results from the equilibrium between fusion and fission processes, controlled by a family of “mitochondria-shaping” proteins. It is becoming clear that defects in mitochondrial dynamics can impair mitochondrial respiration, morphology and motility, leading to apoptotic cell death in vitro and more or less severe neurodegenerative disorders in vivo in humans. Mutations in OPA1, a nuclear encoded mitochondrial protein, cause autosomal Dominant Optic Atrophy (DOA), a heterogeneous blinding disease characterized by retinal ganglion cell degeneration leading to optic neuropathy (Delettre et al., 2000; Alexander et al., 2000). OPA1 is a mitochondrial dynamin-related guanosine triphosphatase (GTPase) protein involved in mitochondrial network dynamics, cytochrome c storage and apoptosis. This protein is anchored or associated on the inner mitochondrial membrane facing the intermembrane space. Eight OPA1 isoforms resulting from alternative splicing combinations of exon 4, 4b and 5b have been described (Delettre et al., 2001). These variants greatly vary among diverse organs and the presence of specific isoforms has been associated with various mitochondrial functions. The different spliced exons encode domains included in the amino-terminal region and contribute to determine OPA1 functions (Olichon et al., 2006). It has been shown that exon 4, that is conserved throughout evolution, confers functions to OPA1 involved in maintenance of the mitochondrial membrane potential and in the fusion of the network. Conversely, exon 4b and exon 5b, which are vertebrate specific, are involved in regulation of cytochrome c release from mitochondria, and activation of apoptosis, a process restricted to vertebrates (Olichon et al., 2007). While Mgm1p has been identified thanks to its role in mtDNA maintenance, it is only recently that OPA1 has been linked to mtDNA stability. Missense mutations in OPA1 cause accumulation of multiple deletions in skeletal muscle. The syndrome associated to these mutations (DOA-1 plus) is complex, consisting of a combination of dominant optic atrophy, progressive external ophtalmoplegia, peripheral neuropathy, ataxia and deafness (Amati- Bonneau et al., 2008; Hudson et al., 2008). OPA1 is the fifth gene associated with mtDNA “breakage syndrome” together with ANT1, PolG1-2 and TYMP (Spinazzola et al., 2009). In this thesis we show for the first time that specific OPA1 isoforms associated to exon 4b are important for mtDNA stability, by anchoring the nucleoids to the inner mitochondrial membrane. Our results clearly demonstrate that OPA1 isoforms including exon 4b are intimately associated to the maintenance of the mitochondrial genome, as their silencing leads to mtDNA depletion. The mechanism leading to mtDNA loss is associated with replication inhibition in cells where exon 4b containing isoforms were down-regulated. Furthermore silencing of exon 4b associated isoforms is responsible for alteration in mtDNA-nucleoids distribution in the mitochondrial network. In this study it was evidenced that OPA1 exon 4b isoform is cleaved to provide a 10kd peptide embedded in the inner membrane by a second transmembrane domain, that seems to be crucial for mitochondrial genome maintenance and does correspond to the second transmembrane domain of the yeasts orthologue encoded by MGM1 or Msp1, which is also mandatory for this process (Diot et al., 2009; Herlan et al., 2003). Furthermore in this thesis we show that the NT-OPA1-exon 4b peptide co-immuno-precipitates with mtDNA and specifically interacts with two major components of the mitochondrial nucleoids: the polymerase gamma and Tfam. Thus, from these experiments the conclusion is that NT-OPA1- exon 4b peptide contributes to the nucleoid anchoring in the inner mitochondrial membrane, a process that is required for the initiation of mtDNA replication and for the distribution of nucleoids along the network. These data provide new crucial insights in understanding the mechanism involved in maintenance of mtDNA integrity, because they clearly demonstrate that, besides genes implicated in mtDNA replications (i.e. polymerase gamma, Tfam, twinkle and genes involved in the nucleotide pool metabolism), OPA1 and mitochondrial membrane dynamics play also an important role. Noticeably, the effect on mtDNA is different depending on the specific OPA1 isoforms down-regulated, suggesting the involvement of two different combined mechanisms. Over two hundred OPA1 mutations, spread throughout the coding region of the gene, have been described to date, including substitutions, deletions or insertions. Some mutations are predicted to generate a truncated protein inducing haploinsufficiency, whereas the missense nucleotide substitutions result in aminoacidic changes which affect conserved positions of the OPA1 protein. So far, the functional consequences of OPA1 mutations in cells from DOA patients are poorly understood. Phosphorus MR spectroscopy in patients with the c.2708delTTAG deletion revealed a defect in oxidative phosphorylation in muscles (Lodi et al., 2004). An energetic impairment has been also show in fibroblasts with the severe OPA1 R445H mutation (Amati-Bonneau et al., 2005). It has been previously reported by our group that OPA1 mutations leading to haploinsufficiency are associated in fibroblasts to an oxidative phosphorylation dysfunction, mainly involving the respiratory complex I (Zanna et al., 2008). In this study we have evaluated the energetic efficiency of a panel of skin fibroblasts derived from DOA patients, five fibroblast cell lines with OPA1 mutations causing haploinsufficiency (DOA-H) and two cell lines bearing mis-sense aminoacidic substitutions (DOA-AA), and compared with control fibroblasts. Although both types of DOA fibroblasts maintained a similar ATP content when incubated in a glucose-free medium, i.e. when forced to utilize the oxidative phosphorylation only to produce ATP, the mitochondrial ATP synthesis through complex I, measured in digitonin-permeabilized cells, was significantly reduced in cells with OPA1 haploinsufficiency only, whereas it was similar to controls in cells with the missense substitutions. Furthermore, evaluation of the mitochondrial membrane potential (DYm) in the two fibroblast lines DOA-AA and in two DOA-H fibroblasts, namely those bearing the c.2819-2A>C mutation and the c.2708delTTAG microdeletion, revealed an anomalous depolarizing response to oligomycin in DOA-H cell lines only. This finding clearly supports the hypothesis that these mutations cause a significant alteration in the respiratory chain function, which can be unmasked only when the operation of the ATP synthase is prevented. Noticeably, oligomycin-induced depolarization in these cells was almost completely prevented by preincubation with cyclosporin A, a well known inhibitor of the permeability transition pore (PTP). This results is very important because it suggests for the first time that the voltage threshold for PTP opening is altered in DOA-H fibroblasts. Although this issue has not yet been addressed in the present study, several are the mechanisms that have been proposed to lead to PTP deregulation, including in particular increased reactive oxygen species production and alteration of Ca2+ homeostasis, whose role in DOA fibroblasts PTP opening is currently under investigation. Identification of the mechanisms leading to altered threshold for PTP regulation will help our understanding of the pathophysiology of DOA, but also provide a strategy for therapeutic intervention.

The effect of osmolytes on protein fold stability at the single-molecule level

By pulling and releasing the tension on protein homomers with the Atomic Force Miscroscope (AFM) at different pulling speeds, dwell times and dwell distances, the observed force-response of the protein can be fitted with suitable theoretical models. In this respect we developed mathematical procedures and open-source computer codes for driving such experiments and fitting Bell’s model to experimental protein unfolding forces and protein folding frequencies. We applied the above techniques to the study of proteins GB1 (the B1 IgG-binding domain of protein G from Streptococcus) and I27 (a module of human cardiac titin) in aqueous solutions of protecting osmolytes such as dimethyl sulfoxide (DMSO), glycerol and trimethylamine N-oxide (TMAO). In order to get a molecular understanding of the experimental results we developed an Ising-like model for proteins that incorporates the osmophobic nature of their backbone. The model benefits from analytical thermodynamics and kinetics amenable to Monte-Carlo simulation. The prevailing view used to be that small protecting osmolytes bridge the separating beta-strands of proteins with mechanical resistance, presumably shifting the transition state to significantly higher distances that correlate with the molecular size of the osmolyte molecules. Our experiments showed instead that protecting osmolytes slow down protein unfolding and speed-up protein folding at physiological pH without shifting the protein transition state on the mechanical reaction coordinate. Together with the theoretical results of the Ising-model, our results lend support to the osmophobic theory according to which osmolyte stabilisation is a result of the preferential exclusion of the osmolyte molecules from the protein backbone. The results obtained during this thesis work have markedly improved our understanding of the strategy selected by Nature to strengthen protein stability in hostile environments, shifting the focus from hypothetical protein-osmolyte interactions to the more general mechanism based on the osmophobicity of the protein backbone.

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.

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.

OPA1 isoforms and protein domains in the rescue of mitochondrial dysfunctions

Mutations in OPA1 gene have been identified in the majority of patients with Dominant Optic Atrophy (DOA), a blinding disease, and the syndromic form DOA-plus. OPA1 protein is a mitochondrial GTPase involved in various mitochondrial functions, present in humans in eight isoforms, resulting from alternative splicing and proteolytic processing. In this study we have investigated the specific role of each isoform through expression in OPA-/- MEFs, by evaluating their ability to improve the defective mitochondrial phenotypes. All isoforms were able to rescue the energetic efficiency, mitochondrial DNA (mtDNA) content and cristae integrity, but only the presence of both long and short forms could recover the mitochondrial morphology. In order to identify the OPA1 protein domains crucial for its functions, we selected and modified the isoform 1, shown to be one of the most efficient in preserving mitochondrial phenotype, to express three specific OPA1 variants, namely: one with a different N-terminus portion, one unable to generate short form owing to deletion of S1 cleavage site and one with a defective GTPase domain. We demonstrated that the simultaneous presence of the N- and C-terminus of OPA1 was essential for the mtDNA maintenance; a cleavable isoform generating s-forms was necessary to completely rescue the energetic competence and the presence of the C-terminus was sufficient to partially recover the cristae ultrastructure. Lastly, several pathogenic OPA1 mutations were inserted in MEF clones and the biochemical features investigated, to correlate the defective phenotypes with the clinical severity of patients. Our results clearly indicate that this cell model reflects very well the clinical characteristics of the patients, and therefore can be proposed as an useful tool to shed light on the pathomechanism underlying DOA.