Machine-learning methods for structure prediction of β-barrel membrane proteins

Different types of proteins exist with diverse functions that are essential for living organisms. An important class of proteins is represented by transmembrane proteins which are specifically designed to be inserted into biological membranes and devised to perform very important functions in the cell such as cell communication and active transport across the membrane. Transmembrane β-barrels (TMBBs) are a sub-class of membrane proteins largely under-represented in structure databases because of the extreme difficulty in experimental structure determination. For this reason, computational tools that are able to predict the structure of TMBBs are needed. In this thesis, two computational problems related to TMBBs were addressed: the detection of TMBBs in large datasets of proteins and the prediction of the topology of TMBB proteins. Firstly, a method for TMBB detection was presented based on a novel neural network framework for variable-length sequence classification. The proposed approach was validated on a non-redundant dataset of proteins. Furthermore, we carried-out genome-wide detection using the entire Escherichia coli proteome. In both experiments, the method significantly outperformed other existing state-of-the-art approaches, reaching very high PPV (92%) and MCC (0.82). Secondly, a method was also introduced for TMBB topology prediction. The proposed approach is based on grammatical modelling and probabilistic discriminative models for sequence data labeling. The method was evaluated using a newly generated dataset of 38 TMBB proteins obtained from high-resolution data in the PDB. Results have shown that the model is able to correctly predict topologies of 25 out of 38 protein chains in the dataset. When tested on previously released datasets, the performances of the proposed approach were measured as comparable or superior to the current state-of-the-art of TMBB topology prediction.

From "wet biology" to statistical analysis of structural features with bioinformatics tools

Many new Escherichia coli outer membrane proteins have recently been identified by proteomics techniques. However, poorly expressed proteins and proteins expressed only under certain conditions may escape detection when wild-type cells are grown under standard conditions. Here, we have taken a complementary approach where candidate outer membrane proteins have been identified by bioinformatics prediction, cloned and overexpressed, and finally localized by cell fractionation experiments. Out of eight predicted outer membrane proteins, we have confirmed the outer membrane localization for five—YftM, YaiO, YfaZ, CsgF, and YliI—and also provide preliminary data indicating that a sixth—YfaL—may be an outer membrane autotransporter.

Interrogation of human monoclonal antibodies induced by 4C-MenB to identify protective antigens contained in the OMV component

Neisseria meningitidis is a gram negative human obligated pathogen, mostly found as a commensal in the oropharyngeal mucosa of healthy individuals. It can invade this epithelium determining rare but devastating and fast progressing outcomes, such as meningococcal meningitidis and septicemia, leading to death (about 135000 per year worldwide). Conjugated vaccines for serogroups A, C, W135, X and Y were developed, while for N. meningitidis serogroup B (MenB) the vaccines were based on Outern Membrane Vesicles (OMV). One of them is the 4C-MenB (Bexsero). The antigens included in this vaccine’s formulation are, in addition to the OMV from New Zeland epidemic strain 98/254, three recombinant proteins: NadA, NHBA and fHbp. While the role of these recombinant components was deeply characterized, the vesicular contribution in 4C-MenB elicited protection is mediated mainly by porin A and other unidentified antigens. To unravel the relative contribution of these different antigens in eliciting protective antibody responses, we isolated human monoclonal antibodies (mAbs) from single-cell sorted plasmablasts of 3 adult vaccinees peripheral blood. mAbs have been screened for binding to 4C-MenB components by Luminex bead-based assay. OMV-specific mAbs were purified and tested for functionality by serum bactericidal assay (SBA) on 18 different MenB strains and characterized in a protein microarray containing a panel of prioritized meningococcal proteins. The bactericidal mAbs identified to recognize the outer membrane proteins PorA and PorB, stating the importance of PorB in cross-strain protection. In addition, RmpM, BamE, Hyp1065 and ComL were found as immunogenic components of the 4C-MenB vaccine.

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.

Use of E. coli OMVs as delivery system to elicit neutralizing antibodies against Chlamydia infectivity: the HtrA protein example

The obligate intracellular pathogen Chlamydia trachomatis is a gram negative bacterium which infects epithelial cells of the reproductive tract. C. trachomatis is the leading cause of bacterial sexually transmitted disease worldwide and a vaccine against this pathogen is highly needed. Many evidences suggest that both antigen specific-Th1 cells and antibodies may be important to provide protection against Chlamydia infection. In a previous study we have identified eight new Chlamydia antigens inducing CD4-Th1 and/or antibody responses that, when combined properly, can protect mice from Chlamydia infection. However, all selected recombinant antigens, upon immunization in mice, elicited antibodies not able to neutralize Chlamydia infectivity in vitro. With the aim to improve the quality of the immune response by inducing effective neutralizing antibodies, we used a novel delivery system based on the unique capacity of E. coli Outer Membrane Vesicles (OMV) to present membrane proteins in their natural composition and conformation. We have expressed Chlamydia antigens, previously identified as vaccine candidates, in the OMV system. Among all OMV preparations, the one expressing HtrA Chlamydia antigen (OMV-HtrA), showed to be the best in terms of yield and quantity of expressed protein, was used to produce mice immune sera to be tested in neutralization assay in vitro. We observed that OMV-HtrA elicited specific antibodies able to neutralize efficiently Chlamydia infection in vitro, indicating that the presentation of the antigens in their natural conformation is crucial to induce an effective immune response. This is one of the first examples in which antibodies directed against a new Chlamydia antigen, other than MOMP (the only so far known antigen inducing neutralizing antibodies), are able to block the Chlamydia infectivity in vitro. Finally, by performing an epitope mapping study, we investigated the specificity of the antibody response induced by the recombinant HtrA and by OMV-HtrA. In particular, we identified some linear epitopes exclusively recognized by antibodies raised with the OMV-HtrA system, detecting in this manner the antigen regions likely responsible of the neutralizing effect.

Molecular variability of Meningococcal antigens in carriage and disease strains and in other Neisseria species

This PhD thesis is focused on the study of the molecular variability of some specific proteins, part of the outer membrane of the pathogen Neisseria meningitidis, and described as protective antigens and important virulence factors. These antigens have been employed as components of the vaccine developed by Novartis Vaccines against N. meningitidis of serogroup B, and their variability in the meningococcal population is a key aspect when the effect of the vaccine is evaluated. The PhD project has led to complete three major studies described in three different manuscritps, of which two have been published and the third is in preparation. The thesis is structured in three main chapters, each of them dedicated to the three studies. The first, described in Chapter 1, is specifically dedicated to the analysis of the molecular conservation of meningococcal antigens in the genomes of all species classified in the genus Neisseria (Conservation of Meningococcal Antigens in the Genus Neisseria. A. Muzzi et al.. 2013. mBio 4 (3)). The second study, described in Chapter 2, focuses on the analysis of the presence and conservation of the antigens in a panel of bacterial isolates obtained from cases of the disease and from healthy individuals, and collected in the same year and in the same geographical area (Conservation of fHbp, NadA, and NHBA in carrier and pathogenic isolates of Neisseria meningitidis collected in the Czech Republic in 1993. A. Muzzi et al.. Manuscript in preparation). Finally, Chapter 3 describes the molecular features of the antigens in a panel of bacterial isolates collected over a period of 50 years, and representatives of the epidemiological history of meningococcal disease in the Netherlands (An Analysis of the Sequence Variability of Meningococcal fHbp, NadA and NHBA over a 50-Year Period in the Netherlands. S. Bambini et al.. 2013. PloS one e65043).

Dissection of the protective response elicited by the DOMV component of the 4CMenB vaccine

Neisseria meningitidis serogroup B is the major etiological agent of meningitis and life-threatening sepsis, against which two vaccines are licensed. The 4CMenB vaccine is composed of three major protein antigens (fHbp, NHBA and NadA) and detergent-extracted outer membrane vesicles (DOMV) from the NZ98/254 strain. DOMV are safe, immunogenic and able to raise bactericidal antibodies, mainly attributed to the immunodominant PorA protein. Nevertheless, DOMV offer a complex reservoir of potentially immunogenic proteins, whose relative contribution in protection is still poorly characterized. By testing antisera from vaccinated infants in serum bactericidal assay, we observed that the addition of DOMV in the vaccine formulation enhanced breadth of coverage compared to recombinant proteins alone against a panel of 11 meningococcal strains mismatched for the vaccine antigens. To unravel the DOMV components involved in such protection, 30 DOMV antigens were cloned and expressed in Escherichia coli as recombinant proteins and/or in vesicles to maintain their native conformation. Samples obtained were both included in tailor-made protein-microarrays to immunoprofile the antibody repertoire raised by DOMV-containing formulations and were individually used for mouse immunization studies to assess their ability to induce bactericidal antibodies. The protein-array immunosignature of mouse DOMV/4CMenB antisera unveiled a subset of 8 DOMV-reactive proteins potentially responsible for the additional protective responses. The antisera derived from mouse immunizations showed high levels of antibodies and recognized the corresponding antigen across different meningococcal strains. Among the protein-array reactive antigens, OpcA, NspA and PorB induced antibodies able to kill 10 of the 11 genetically diverse meningococcal strains and the specificity of the protective role of OpcA and PorB was also confirmed in 4CMenB infant vaccinee sera. In conclusion, we identified additional PorA-independent antigens within DOMV involved in broadening the coverage of 4CMenB, thus supporting the key role played by vesicles in this multivalent formulation.

Machine learning tools for protein annotation: the cases of transmembrane β-barrel and myristoylated proteins

Biology is now a “Big Data Science” thanks to technological advancements allowing the characterization of the whole macromolecular content of a cell or a collection of cells. This opens interesting perspectives, but only a small portion of this data may be experimentally characterized. From this derives the demand of accurate and efficient computational tools for automatic annotation of biological molecules. This is even more true when dealing with membrane proteins, on which my research project is focused leading to the development of two machine learning-based methods: BetAware-Deep and SVMyr. BetAware-Deep is a tool for the detection and topology prediction of transmembrane beta-barrel proteins found in Gram-negative bacteria. These proteins are involved in many biological processes and primary candidates as drug targets. BetAware-Deep exploits the combination of a deep learning framework (bidirectional long short-term memory) and a probabilistic graphical model (grammatical-restrained hidden conditional random field). Moreover, it introduced a modified formulation of the hydrophobic moment, designed to include the evolutionary information. BetAware-Deep outperformed all the available methods in topology prediction and reported high scores in the detection task. Glycine myristoylation in Eukaryotes is the binding of a myristic acid on an N-terminal glycine. SVMyr is a fast method based on support vector machines designed to predict this modification in dataset of proteomic scale. It uses as input octapeptides and exploits computational scores derived from experimental examples and mean physicochemical features. SVMyr outperformed all the available methods for co-translational myristoylation prediction. In addition, it allows (as a unique feature) the prediction of post-translational myristoylation. Both the tools here described are designed having in mind best practices for the development of machine learning-based tools outlined by the bioinformatics community. Moreover, they are made available via user-friendly web servers. All this make them valuable tools for filling the gap between sequential and annotated data.

Molecular basis oh herpes simplex virus entry into the cell and retargeting of the viral tropism for the design of oncolytic herpesviruses

Herpes simplex virus 1 (HSV-1) infects oral epitelial cells, then spreads to the nerve endings and estabilishes latency in sensory ganglia, from where it may, or may not reactivate. Diseases caused by virus reactivation include mild diseases such as muco-cutaneous lesions, and more severe, and even life-threatening encephalitis, or systemic infections affecting diverse organs. Herpes simplex virus represents the most comprehensive example of virus receptor interaction in Herpesviridae family, and the prototype virus encoding multipartite entry genes. In fact, it encodes 11-12 glycoproteins and a number of additional membrane proteins: five of these proteins play key roles in virus entry into subsceptible cells. Thus, glycoprotein B (gB) and glycoprotein C (gC) interact with heparan sulfate proteoglycan to enable initial attachment to cell surfaces. In the next step, in the entry cascade, gD binds a specific surface receptor such as nectin1 or HVEM. The interaction of glycoprotein D with the receptor alters the conformation of gD to enable the activation of gB, glycoprotein H, and glycoprotein L, a trio of glycoproteins that execute the fusion of the viral envelope with the plasma membrane. In this thesis, I described two distinct projects: I. The retargeting of viral tropism for the design of oncolytic Herpesviruses: • capable of infecting cells through the human epitelial growth factor receptor 2 (HER2), overexpressed in highly malignant mammary and ovarian tumors and correlates with a poor prognosis; • detargeted from its natural receptors, HVEM and nectin1. To this end, we inserted a ligand to HER2 in gD. Because HER2 has no natural ligand, the selected ligand was a single chain antibody (scFv) derived from MAb4D5 (monoclonal antibody to HER2), herein designated scHER2. All recombinant viruses were targeted to HER2 receptor, but only two viruses (R-LM113 and R-LM249) were completely detargeted from HVEM and nectin1. To engineer R-LM113, we removed a large portion at the N-terminus of gD (from aa 6 to aa 38) and inserted scHER2 sequence plus 9-aa serine-glycine flexible linker at position 39. On the other hand, to engineer R-LM249, we replaced the Ig-folded core of gD (from aa 61 to aa 218) with scHER2 flanked by Ser-Gly linkers. In summary, these results provide evidence that: i. gD can tolerate an insert almost as big as gD itself; ii. the Ig-like domain of gD can be removed; iii. the large portion at the N-terminus of gD (from aa 6 to aa 38) can be removed without loss of key function; iv. R-LM113 and R-LM249 recombinants are ready to be assayed in animal models of mammary and ovary tumour. This finding and the avaibility of a large number of scFv greatly increase the collection of potential receptors to which HSV can be redirected. II. The production and purification of recombinant truncated form of the heterodimer gHgL. We cloned a stable insect cell line expressing a soluble form of gH in complex with gL under the control of a metalloprotein inducible promoter and purified the heterodimer by means of ONE-STrEP-tag system by IBA. With respect to biological function, the purified heterodimer is capable: • of reacting to antibodies that recognize conformation dependent epitopes and neutralize virion infectivity; • of binding a variety cells at cell surface. No doubt, the availability of biological active purified gHgL heterodimer, in sufficient quantities, will speed up the efforts to solve its crystal structure and makes it feasible to identify more clearly whether gHgL has a cellular partner, and what is the role of this interaction on virus entry.