published_or_final_version / Paediatrics and Adolescent Medicine / Doctoral / Doctor of Philosophy

Influenza B virus is a major human pathogen causing highly contagious respiratory disease. It accounts for approximately ~30% of influenza virus infection per year. The effector domain of the NS1 protein of influenza B virus (NS1B protein), encompassing the carboxy terminal two thirds of the protein, suppresses interferon-β (IFN-β) synthesis in virus-infected cells by unknown mechanism(s). The induced IFN-β mediates innate immunity. To elucidate the mechanism by which the NS1B effector domain suppresses the production of IFN-β, we identified cellular proteins that interact with the NS1B effector domain. Two approaches were used. The approach that succeeded employed the transfection into cells of plasmids expressing the NS1B effector domain containing two affinity tags. After double affinity purification, co-purified cellular proteins were identified by mass spectrometry. We identified Brd2 as a cellular protein that interacts with the NS1B protein. We established that Brd2 specifically binds to the NS1B effector domain in vitro, in vivo, and in virus-infected cells. Serial mutagenesis experiments showed the phenylalanine at position 171 (F171) of the NS1B protein is essential for Brd2 binding. To determine the function of the interaction of Brd2 with the NS1B protein, we generated a recombinant virus encoding an NS1B protein in which F at position 171 was replaced by an alanine. The F171A mutant virus was attenuated, and unlike the wild-type virus, induced the synthesis of IFN-β mRNA. IRF3, a key transcription factor for transcription of the IFN-β gene, was activated in mutant virusinfected cells, but not in wild-type virus-infected cells. Transfection assays implicated the activation of the TBK1 kinase as the step in IRF3 activation that is induced in mutant virus-infected cells. We interpreted these results as showing that Brd2 binding to the NS1B protein is required for suppressing IRF3 activation and IFN-β induction. Attempts at further confirmation by depletion of endogenous Brd2 using RNA interference were not successful because of inefficient knock-down efficiency and nonspecific IFN-β induction. A further complication is that another bromodomain protein, Brd4, interacts with the NS1B protein and could compensate for depletion of Brd2. / text

An increasing number of proteins with unknown function have their three-dimensional structure solved at high resolution. This situation, largely due to structural genomics initiatives, has been stimulating the development of automated structure-based function prediction methods. Knowledge of residues important for function – and more particularly – for binding can help automated prediction of function in different ways. The properties of a binding site such as its shape or amino acid composition can provide clues on the ligand that may bind to it. Also, having information on functionally important regions in similar proteins can refine the process of annotation transfer between homologues. Experimental results indicate that functional residues often have an unfavourable contribution to the stability of the folded state of a protein. This observation is the underlying principle of several computational methods for predicting the location of functional sites in protein structures. These methods search protein structures for destabilising residues, with the assumption that these are likely to be important for function. We have developed a method to detect clusters of destabilising residues which are in close spatial proximity within a protein structure. Individual residue contributions to protein stability are evaluated using detailed atomic models and an energy function based on fundamental physico-chemical principles. Our overall aim in this work was to evaluate the overlap between these clusters of destabilising residues and known binding sites in proteins. Unfortunately, reliable benchmark datasets of known binding sites in proteins are sorely lacking. Therefore, we have undertaken a comprehensive approach to define binding sites unambiguously from structural data. We have rigorously identified seven issues which should be considered when constructing datasets of binding sites to validate prediction methods, and we present the construction of two new datasets in which these problems are handled. In this regard, our work constitute a major improvement over previous studies in the field. Our first dataset consists of 70 proteins with binding sites for diverse types of ligands (e.g. nucleic acids, metal ions) and was constructed using all available data, including literature curation. The second dataset contains 192 proteins with binding sites for small ligands and polysaccharides, does not require literature curation, and can therefore be automatically updated. We have used our dataset of 70 proteins to evaluate the overlap between destabilising regions and binding sites (the second dataset of 192 proteins was not used for that evaluation as it constitutes a later improvement). The overlap is on average limited but significantly larger than random. The extent of the overlap varies with the type of bound ligand. Significant overlap is obtained for most polysaccharide- and small ligand-binding sites, whereas no overlap is observed for nucleic acid-binding sites. These differences are rationalised in terms of the geometry and energetics of the binding sites. Although destabilising regions, as detected in this work, can in general not be used to predict all types of binding sites in protein structures, they can provide useful information, particularly on the location of binding sites for polysaccharides and small ligands. In addition, our datasets of binding sites in proteins should help other researchers to derive and validate new function prediction methods. We also hope that the criteria which we use to define binding sites may be useful in setting future standards in other analyses.

Domain swapping is a wide spread phenomenon which involves the association between two or more protein subunits such that intra-molecular interactions between domains in each subunit are replaced by equivalent inter-molecular interactions between the same domains in different subunits. This thesis is devoted to the analysis of the factors that drive proteins to undergo such association modes. The specific system analyzed is the monomer to swapped dimer formation of the B1 domain of the immunoglobulin G binding protein (GB1). The formation of this dimer was shown to be fostered by 4 amino acid substitutions (L5V, F30V, Y33F, A34F) (Byeon et al., 2003). In this work, computational protein design and molecular dynamics simulations, both with detailed atomic models, were used to gain insight into how these 4 mutations may promote the domain swapping reaction. The stability of the wt and quadruple mutant GB1 monomers was assessed using the software DESIGNER, a fully automatic procedure that selects amino acid sequences likely to stabilize a given backbone structure (Wernisch et al., 2000). Results suggest that 3 of the mutations (L5V, F30V, A34F) have a destabilizing effect. The first mutation (L5V) forms destabilizing interactions with surrounding residues, while the second (F30V) is engaged in unfavorable interactions with the protein backbone, consequently causing local strain. Although the A34F substitution itself is found to contribute favorably to the stability of the monomer, this is achieved only at the expense of forcing the wild type W43 into a highly strained conformation concomitant with the formation of unfavorable interactions with both W43 and V54. Finally, we also provide evidence that A34F mutation stabilizes the swapped dimer structure. Although we were unable to perform detailed protein design calculations on the dimer, due to the lower accuracy of the model, inspection of its 3D structure reveals that the 34F side chains pack against one another in the core of the swapped structure, thereby forming extensive non-native interactions that have no counterparts in the individual monomers. Their replacement by the much smaller Ala residue is suggested to be significantly destabilizing by creating a large internal cavity, a phenomenon, well known to be destabilizing in other proteins. Our analysis hence proposes that the A34F mutation plays a dual role, that of destabilizing the GB1 monomer structure while stabilizing the swapped dimer conformation. In addition to the above study, molecular dynamics simulations of the wild type and modeled quadruple mutant GB1 structures were carried out at room and elevated temperatures (450 K) in order to sample the conformational landscape of the protein near its native monomeric state, and to characterize the deformations that occur during early unfolding. This part of the study was aimed at investigating the influence of the amino acid sequence on the conformational properties of the GB1 monomer and the possible link between these properties and the swapping process. Analysis of the room temperature simulations indicates that the mutant GB1 monomer fluctuates more than its wild type counter part. In addition, we find that the C-terminal beta-hairpin is pushed away from the remainder of the structure, in agreement with the fact that this hairpin is the structural element that is exchanged upon domain swapping. The simulations at 450 K reveal that the mutant protein unfolds more readily than the wt, in agreement with its decreased stability. Also, among the regions that unfold early is the alpha-helix C-terminus, where 2 out of the 4 mutations reside. NMR experiments by our collaborators have shown this region to display increased flexibility in the monomeric state of the quadruple mutant. Our atomic scale investigation has thus provided insights into how sequence modifications can foster domain swapping of GB1. Our findings indicate that the role of the amino acid substitutions is to decrease the stability of individual monomers while at the same time increase the stability of the swapped dimer, through the formation of non-native interactions. Both roles cooperate to foster swapping.

In the past decade, automatisation has led to an immense increase of data in biology. Next generation sequencing techniques will produce a vast amount of sequences across all species in the coming years. In many cases, identifying the function and biological role of a protein from its sequence can be a complicated and time-intensive task. The identification of a protein's interaction partners is a tremendous help for understanding the biological context in which it is involved. In order to fully characterise a protein-protein interaction (PPIs), it is necessary to know the three-dimensional structure of the interacting partners. Despite optimisation efforts from projects such as the Protein Structure Initivative, determining the structure of a protein through crystallography remains a time- and cost-intensive procedure. The primary aim of the research described in this dissertation was to produce a World Wide Web resource that facilitates visual exploration and validation (or questioning) of data derived from functional genomics experiments, by building upon existing structural information about direct physical PPIs. Secondary aims were (i) to demonstrate the utility of the new resource, and (ii) its application in biological research. We created a database that emphasises specifically the intersection between the PPIs-results emerging from the structural biology and functional genomics communities. The BISC database holds BInary SubComplexes and Modellable Interactions in current functional genomics databases (BICS-MI). It is publicly available at hyyp://bisc.cse.ucsc.edu. BISC is divided in three sections that deliver three types of information of interest to users seeking to investigate or browse PPIs. The template section (BISCHom and BISCHet) is devoted to those PPIs that are characterised in structural detail, i.e. binary SCs extracted from experimentally determined three-dimensional structures. BISCHom and BISCHet contain the homodimeric (13,583 records) and heterodimeric (5612 records) portions of these, respectively. Besides interactive, embedded Jmol displays emphasising the interface, standard information and links are provided, e.g. sequence information and SPOP classification for both partners, interface size and energy scores (PISA). An automated launch of the MolSurfer program enables the user to investigate electrostatic and hydrophobic correlation between the partners, at the inter-molecular interface. The modellable interactions section (BISC0MI) identifies potentially modellable interactions in three major functional genomics interaction databases (BioGRID), IntAct, HPRD). To create BISC-MI all PPIs that are amenable to automated homology modelling based on conservative similarity cut-offs and whose partner protein sequences have recrods in the UniProt database, have been extracted. The modellable interaction services (BISC-MI Services) section offers, upon user request, modelled SC-structures for any PPIs in BISC-MI. This is enabled through an untomated template-based (homology) modelling protocol using the popular MODELLER program. First, a multiple sequence alignment (MSA) is generated using MUSCLE, between the target and homologous proteins collected from UniProt (only reviewed proteins from organisms whose genome has been completely sequenced are included to find putative orthologs). Then a sequence-to-profile alignment is generated to integrate the template structure in the MSA. All models are produced upon user request to ensure that the most recent sequence data for the MSAs are used. Models generated through this protocol are expected to be more accurate generally than models offered by other automated resources that rely on pairwise alignments, e.g. ModBase. Two small studies were carried out to demonstrate the usability and utility of BISC in biological research. (1) Interaction data in functional genomics databases often suffers from insufficient experimental and reporting standards. For example, multiple protein complexes are typically recorded as an inferred set of binary interactions. Using the 20S core particle of the yeast proteasome as an example, we demonstrate how the BISC Web resource can be used as a starting point for further investigation of such inferred interactions. (2) Malaria, a mosquito-borne disease, affects 3500-500 million people worldwide. Still very little is known about the malarial parasites' genes and their protein functions. For Plasmodium falciparum, the most lethal among the malaria parasites, only one experimentally derived medium scale PPIs set is available. The validity of this set has been doubted in the the malarial biologist community. We modelled and investigated eleven binary interactions from this set using the BISC modelling pipeline. Alongside we compared the BISC models of the individual partners to those obtained from ModBase.

Adiponectin, a hormone that homo-oligomerizes into trimer, hexamer, or higher molecular weight (HMW) species, is involved in maintaining insulin sensitivity in muscle and liver. Interestingly, its functions appear to be oligomer-specific. Recent data suggest that HMW levels are decreased in obesity and insulin resistance, although, the cause for this decrease is not known. Impaired assembly to the octadecamer represents one possible reason for decreased HMW adiponectin in insulin resistance and type 2 diabetes, but mechanisms by which HMW adiponectin assembles are unknown. This dissertation discusses the progress that we have made regarding formation of HMW adiponectin in vitro.I found that disulfide bonds are important in the assembly process to octadecameric adiponectin, but are not required for stability of the octadecamer itself. We showed that hydrogen peroxide accelerated oligomerization to the octadecamer through formation of disulfide bonds, while alkylation of the cysteines led to inhibition of both oligomerization and disulfide bond formation. Using comparative native/denaturing polyacrylamide gel electrophoresis (PAGE), dynamic light scattering, and tandem mass spectrometry, we demonstrated that octadecamer is stable in the absence of disulfide bonds by using multiple biochemical and biophysical assays. In addition, oxidized adiponectin oligomerizes to octadecamer far slower than reduced adiponectin. To further evaluate the role of disulfide bonds in the formation to octadecamer, we analyzed the role of reduction potential on adiponectin oligomerization. We observed that under immediate oxidizing conditions, hexamers and trimers form. Oxidized hexamer can form HMW adiponectin through disulfide bond rearrangement using beta-mercaptoethanol (βME) or increasing the total concentration of glutathione under oxidizing conditions. To further understand the role of disulfide bonds, we showed that zinc increased the oligomerization to octadecamer. This effect was associated with decreased initial disulfide bonding during the assembly to the octadecamer. In summary, these data suggest the rate of disulfide bond formation and the ability to undergo disulfide bond isomerization are important in the oligomerization process of HMW adiponectin.

Eukaryotic protein kinases are one of the most important classes of human proteins, and a great deal of research has focused on the development of small molecule inhibitors as biological probes for the determination of their cellular function or as therapeutics for the treatment of disease, such as cancer. The need for new selective inhibitors and a better understanding of the selectivities of existing small molecules is readily apparent. Towards the goal of better understanding protein kinases and the molecules that inhibit them, I have developed a split-protein-based approach for the investigation of these kinase-small molecule interactions. Employing split-firefly luciferase as a reporter domain, we engineered a three-hybrid system capable of determining kinase inhibition through competitive interactions between an active site-directed ligand and a small molecule of interest. This method measures luciferase activity as a function of ligand binding, as opposed to the more traditional assays which quantify kinase activity directly, and alleviates the laborious process of protein purification. The model kinase PKA and the promiscuous ligand staurosporine were used in an initial test case to successfully validate the general design principles of our assay. The modular nature inherent to the assay's design enabled us to adapt it to roughly 300 additional protein kinases and two different ligands. We were able to establish a protocol for rapidly ascertaining the inhibition of a kinase by a library of 80 commercially available kinase inhibitors in a 96-well, high-throughput format. This protocol was then systematically applied to the AGC group of kinases to observe patterns of inhibition across similarly related kinases. We have further shown how these results might be correlated with the sequence identity between kinases to better anticipate inhibitor promiscuity. Finally, we were able to illustrate how a kinase-centric approach could be applied to correlate alterations to the kinase domain with changes in luminescence. This has use for the interrogation of different modes of inhibition as well as in identifying the specific determinants of inhibitor binding. In total, these efforts represent the optimization of a new, general platform for determining kinase inhibitor selectivity across the kinome, and it could potentially be applied universally to the interrogation of protein-ligand interactions.

Non-competitive antagonists of the glutamate N-methyl-D-aspartate receptor (NMDAR) induce a broad range of schizophrenia-like symptoms in humans. Consequently hypothesis has emerged suggesting that glutamate or NMDAR hypofunction may occur in schizophrenia. The NMDAR is localised at dendritic spines of neurons and is embedded in a multi-protein complex called the post-synaptic density (PSD). The biochemical composition of the postsynaptic membrane and the structure of dendritic spines are continuously modulated by glutamatergic synaptic activity. The activity-dependent interaction between glutamate receptors and proteins of the PSD stimulate intracellular signalling pathways underlying learning and memory processes. These may be disturbed in schizophrenia. In the present study we hypothesised that molecules of the PSD may be disturbed in expression in the premotor cortex of patients with schizophrenia. Postmortem premotor cortex from patients with schizophrenia, major depressive disorder, bipolar disorder and healthy controls were processed for PSD extraction and purification. Protein expression of the PSD fraction was assessed using co-immunoprecipitation (co-IP) and Western blotting (WB) methods. The expression of NMDAR subunit NR2A, PSD-95, Ca2+/calmodulin-dependent protein kinase II subunit β (CaMKIIβ) and truncated isoform of the tropomyosin receptor kinase type B (TrkB-T1) were significantly reduced in schizophrenia. A significant decrease in the expression of NR2A was also observed in patients with major depressive disorder relative to controls. A decrease in the abundance of key PSD proteins in schizophrenia provides strong evidence that PSD function and possibly synaptic plasticity may be disturbed in the premotor cortex in the disease. There may also be more subtle disturbances in PSD function in major depressive disorder.

Herpesviruses are ubiquitous in both humans and animals and can cause life-threatening disease. The discovery of murine gammaherpesvirus 68 (MHV-68), which has many similarities in genome and pathogenesis as the human pathogens Epstein-Barr virus and Kaposi’s sarcoma-associated herpesvirus, provides a model for further investigation of the pathogenesis of gammaherpesviruses. The M4 gene was found to be at the left end region of MHV-68 genome. The presence of the M4 protein is required during the early establishment of MHV-68 latency. However, the function of M4 protein remains unclear. The aim of this project was to investigate the function of the M4 protein in vitro and during infection. By using an ELISA, the recombinant M4 protein was shown to bind several Cxc-chemokines and stop the interaction between Cxcl4 and GAGs. The role of M4 protein during MHV-68 lytic infection and in the early establishment of latency was studied by comparing the pathogenesis of virus which does not express M4 (M4stop) and wild type virus (WT). Compared to WT infection, this study found that M4stop was decreased in the lungs at day 8 post infection (p.i.). At the same time point, the viral loads were higher in M4stop infected spleens, which was accompanied by increased expression of the CD4+ T cell activation marker PD-1 and the macrophage activation marker CD69. However, at day 14 p.i., the M4stop infected spleens had lower viral loads, and the expression of CD69 was decreased on CD4+, CD8+ T cells, B cells and macrophages. Furthermore, gene expression PCR arrays were used to investigate how cellular activation and inflammation were transcriptionally regulated. It has been found that the transcription of several genes, which are involved in germinal centre development, was lower in the spleens of WT infected mice at day 12 and 14 p.i. compared to day 10 p.i. of WT infection, as well as day 12 and 14 p.i. of M4stop infection. In addition, the percentage of germinal centre B cells was found to be higher in spleens infected with M4stop at day 10 p.i.. However, there was no difference in percentages of TFH and plasma cells in the spleens. Finally, in order to understand the role of IFN-γ in control of infection in M4stop infected mice, IFN-γR-/- mice were infected with M4stop and WT. Although there were differences in pathogenesis between WT and M4Stop virus infected IFN-γR-/- mice, there was no clear evidence that M4 function is involved in inhibiting IFN-γ pathways. In this study, we found M4 can disturb the interaction of chemokine and GAGs and might delay virus trafficking to the spleen, which could lead to a reduction of cellular activation. M4 may also impair the development of germinal centres at the beginning of latent infection in the spleens.

Protein-protein interactions (PPIs) are key drug targets and recent breakthroughs in this area are providing insight into the types of molecules needed to selectively and potently inhibit a target traditionally seen as untractable. The rules that have been used to design classic substratecompetitive drugs (for example Lipinski's rule of five) may not apply in this new field in the same way. Here I present work performed in three systems that are well-validated drug targets for oncogenesis: the CDK2/cyclin A complex, the PLK1 Polobox domain and MDM2. In each case the site of the protein-protein interaction is defined and understood and the rationale for pharmaceutical intervention is clear. I use these as a model system to evaluate the characteristics of drugs that target protein-protein interaction sites and present work on the development of inhibitors as potential leads for subsequent drug development. In Chapter 1 I introduce the problems, challenges and rewards of PPI drug development; in Chapter 2 I present co-crystal structures of MDM2 with isoindolinone inhibitors; in Chapter 3 I detail attempts to co-crystallise the Plk1 Polobox with inhibitors and screen potential inhibitors; in Chapter 4 I present the results of screening to identify inhibitors of Cyclin A recruitment; and in Chapter 5 I discuss other strategies for inhibition of the CDK2/cyclin A complex, including results with a covalent inhibitor. Through these projects I have been able to demonstrate the wide applicability of the PPI inhibition approach, identify key features of drugs able to inhibit PPIs and contribute to drug design in each system.