From protein structure to function with bioinformatics

It has long been known that amino acids are the building blocks for proteins and govern their folding into specific three-dimensional structures. However, the details of this process are still unknown and represent one of the main problems in structural bioinformatics, which is a highly active research area with the focus on the prediction of three-dimensional structure and its relationship to protein function. The protein structure prediction procedure encompasses several different steps from searches and analyses of sequences and structures, through sequence alignment to the creation of the structural model. Careful evaluation and analysis ultimately results in a hypothetical structure, which can be used to study biological phenomena in, for example, research at the molecular level, biotechnology and especially in drug discovery and development. In this thesis, the structures of five proteins were modeled with templatebased methods, which use proteins with known structures (templates) to model related or structurally similar proteins. The resulting models were an important asset for the interpretation and explanation of biological phenomena, such as amino acids and interaction networks that are essential for the function and/or ligand specificity of the studied proteins. The five proteins represent different case studies with their own challenges like varying template availability, which resulted in a different structure prediction process. This thesis presents the techniques and considerations, which should be taken into account in the modeling procedure to overcome limitations and produce a hypothetical and reliable three-dimensional structure. As each project shows, the reliability is highly dependent on the extensive incorporation of experimental data or known literature and, although experimental verification of in silico results is always desirable to increase the reliability, the presented projects show that also the experimental studies can greatly benefit from structural models. With the help of in silico studies, the experiments can be targeted and precisely designed, thereby saving both money and time. As the programs used in structural bioinformatics are constantly improved and the range of templates increases through structural genomics efforts, the mutual benefits between in silico and experimental studies become even more prominent. Hence, reliable models for protein three-dimensional structures achieved through careful planning and thoughtful executions are, and will continue to be, valuable and indispensable sources for structural information to be combined with functional data.

Probing MST1 Kinase Sequence and Structure for Rational Drug Discovery (Targeting diabetes by blocking protein-protein interactions)

Apoptotic beta cell death is an underlying cause majorly for type I and to a lesser extent for type II diabetes. Recently, MST1 kinase was identified as a key apoptotic agent in diabetic condition. In this study, I have examined MST1 and closely related kinases namely, MST2, MST3 and MST4, aiming to tackle diabetes by exploring ways to selectively block MST1 kinase activity. The first investigation was directed towards evaluating possibilities of selectively blocking the ATP binding site of MST1 kinase that is essential for the activity of the enzymes. Structure and sequence analyses of this site however revealed a near absolute conservation between the MSTs and very few changes with other kinases. The observed residue variations also displayed similar physicochemical properties making it hard for selective inhibition of the enzyme. Second, possibilities for allosteric inhibition of the enzyme were evaluated. Analysis of the recognized allosteric site also posed the same problem as the MSTs shared almost all of the same residues. The third analysis was made on the SARAH domain, which is required for the dimerization and activation of MST1 and MST2 kinases. MST3 and MST4 lack this domain, hence selectivity against these two kinases can be achieved. Other proteins with SARAH domains such as the RASSF proteins were also examined. Their interaction with the MST1 SARAH domain were evaluated to mimic their binding pattern and design a peptide inhibitor that interferes with MST1 SARAH dimerization. In molecular simulations the RASSF5 SARAH domain was shown to strongly interact with the MST1 SARAH domain and possibly preventing MST1 SARAH dimerization. Based on this, the peptidic inhibitor was suggested to be based on the sequence of RASSF5 SARAH domain. Since the MST2 kinase also interacts with RASSF5 SARAH domain, absolute selectivity might not be achieved.

The effect of age on the structure and compostition of the cell walls of Choanephora cucurbitarum

The effect of age on the structure and composition of isolated and purified cell walls from cultures of Choanephora cucurbitarum was investigated by microchemical analyses, visible and infrared spectrophotometry, x-ray diffractometry and electron microscopy. Qualitative evaluation revealed the presence of lipids, proteins, neutral sugars, strong alkali soluble sugars, chitin, chitosan and uronic acids in the cell walls of both the 1 and 7 day old cultures. As the mycelium aged, there was a slight but statistically significant increase in the protein content, and a pronounced rise in the chitin and neutral sugar constituents of the cell walls. Conversely, the decrease in the chitosan content during this period had the net effect of altering the chitin: chitosan ratio from near unity in the younger cultures, to a 2:1 ratio in the 7 day old cell wall samples. Glutaraldehyde-osmium fixed thin sections of the 1 day old vegetative hyphae of £. curbitarum revealed the presence of a monolayered cell wall, which upon aging became bilayered. Replicas of acid hydrolysed cell walls demonstrated that both the 1 and 7 day old samples possessed an outer layer which was composed of finely granular amorphous material and randomly distributed microfibrils. The deposition of an inner secondary layer composed of parallel oriented microfibrils in the older hypha was correlated with an increase in the chitin content in the cell wall. The significance of these results with respect to the intimate relationship between composition and structure is discussed.

The structure of arabidopsis NPR1 : its function as a salicylic acid receptor and a metal-binding protein

The Arabidopsis NPRI protein regulates systemic acquired resistance dependent on salicylic acid. Analyses by plant two-hybrid analysis in vivo and pull-down assays in vitro showed that the BTB/POZ domain of NPRI at the N-terminus serves as an autoinhibitory domain to negate the function of the transactivation domain at the C-terminus through direct binding of these two domains. I t was also shown that the binding of the BTB/POZ domain to the C-terminus of NPRI was abolished by SA treatment, suggesting that SA could interfere directly with this binding. By gel filtration, it was demonstrated that SA affects the conformation of full-length NPRl , confirming the role of NPRI as an SA receptor. Gel filtration analysis also indicated that NPRI could be converted from an oligomer to a dimer with SA treatment. Furthermore, one N-terminal deletion ~513 has been shown to act as a metal-binding protein and its two Cys-521 and Cys-529 are important for binding to Ni 2 + by pull-down assays.

Modeling Metal Protein Complexes from Experimental Extended X-ray Absorption Fine Structure using Computational Intelligence

Experimental Extended X-ray Absorption Fine Structure (EXAFS) spectra carry information about the chemical structure of metal protein complexes. However, pre- dicting the structure of such complexes from EXAFS spectra is not a simple task. Currently methods such as Monte Carlo optimization or simulated annealing are used in structure refinement of EXAFS. These methods have proven somewhat successful in structure refinement but have not been successful in finding the global minima. Multiple population based algorithms, including a genetic algorithm, a restarting ge- netic algorithm, differential evolution, and particle swarm optimization, are studied for their effectiveness in structure refinement of EXAFS. The oxygen-evolving com- plex in S1 is used as a benchmark for comparing the algorithms. These algorithms were successful in finding new atomic structures that produced improved calculated EXAFS spectra over atomic structures previously found.

Structure-fonction de MARCH1, une E3 ubiquitine ligase régulant la présentation antigénique par le CMH II

Les molécules classiques du CMH de classe II sont responsables de la présentation de peptides exogènes par les cellules présentatrices d’antigène aux lymphocytes T CD4+. Cette présentation antigénique est essentielle à l’établissement d’une réponse immunitaire adaptative. Cependant, la reconnaissance d’auto-antigènes ainsi que l’élimination des cellules du Soi sont des problèmes à l’origine de nombreuses maladies auto-immunes. Notamment, le diabète et la sclérose en plaque. D’éventuels traitements de ces maladies pourraient impliquer la manipulation de la présentation antigénique chez les cellules dont la reconnaissance et l’élimination engendrent ces maladies. Il est donc primordial d’approfondir nos connaissances en ce qui concerne les mécanismes de régulation de la présentation antigénique. La présentation antigénique est régulée tant au niveau transcriptionnel que post-traductionnel. Au niveau post-traductionnel, diverses cytokines affectent le processus. Parmi celles-ci, l’IL-10, une cytokine anti-inflammatoire, cause une rétention intracellulaire des molécules du CMH II. Son mécanisme d’action consiste en l’ubiquitination de la queue cytoplasmique de la chaîne bêta des molécules de CMH II. Cette modification protéique est effectuée par MARCH1, une E3 ubiquitine ligase dont l’expression est restreinte aux organes lymphoïdes secondaires. Jusqu’à tout récemment, il y avait très peu de connaissance concernant la structure et les cibles de MARCH1. Considérant son impact majeur sur la présentation antigénique, nous nous sommes intéressé à la structure-fonction de cette molécule afin de mieux caractériser sa régulation ainsi que les diverses conditions nécessaires à son fonctionnement. Dans un premier article, nous avons étudié la régulation de l’expression de MARCH1 au niveau protéique. Nos résultats ont révélé l’autorégulation de la molécule par formation de dimères et son autoubiquitination. Nous avons également démontré l’importance des domaines transmembranaires de MARCH1 dans la formation de dimères et l’interaction avec le CMH II. Dans un second article, nous avons investigué l’importance de la localisation de MARCH1 pour sa fonction. Les résultats obtenus montrent la fonctionnalité des motifs de localisation de la portion C-terminale de MARCH1 ainsi que la présence d’autres éléments de localisation dans la portion N-terminale de la protéine. Les nombreux mutants utilisés pour ce projet nous ont permis d’identifier un motif ‘‘VQNC’’, situé dans la portion cytoplasmique C-terminale de MARCH1, dont la valine est requise au fonctionnement optimal de la molécule. En effet, la mutation de la valine engendre une diminution de la fonction de la molécule et des expériences de BRET ont démontré une modification de l’orientation spatiale des queues cytoplasmiques. De plus, une recherche d’homologie de séquence a révélé la présence de ce même motif dans d’autres ubiquitines ligases, dont Parkin. Parkin est fortement exprimée dans le cerveau et agirait, entre autre, sur la dégradation des agrégats protéiques. La dysfonction de Parkin cause l’accumulation de ces agrégats, nommés corps de Lewy, qui entraînent des déficiences au niveau du fonctionnement neural observé chez les patients atteints de la maladie de Parkinson. La valine comprise dans le motif ‘’VQNC’’ a d’ailleurs été identifiée comme étant mutée au sein d’une famille où cette maladie est génétiquement transmise. Nous croyons que l’importance de ce motif ne se restreint pas à MARCH1, mais serait généralisée à d’autres E3 ligases. Ce projet de recherche a permis de caractériser des mécanismes de régulation de MARCH1 ainsi que de découvrir divers éléments structuraux requis à sa fonction. Nos travaux ont permis de mieux comprendre les mécanismes de contrôle de la présentation antigénique par les molécules de CMH II.

Nuclear magnetic resonance structure of the nucleic acid-binding domain of severe acute respiratory syndrome coronavirus nonstructural protein 3

The nuclear magnetic resonance (NMR) structure of a globular domain of residues 1071 to 1178 within the previously annotated nucleic acid-binding region (NAB) of severe acute respiratory syndrome coronavirus nonstructural protein 3 (nsp3) has been determined, and N- and C-terminally adjoining polypeptide segments of 37 and 25 residues, respectively, have been shown to form flexibly extended linkers to the preceding globular domain and to the following, as yet uncharacterized domain. This extension of the structural coverage of nsp3 was obtained from NMR studies with an nsp3 construct comprising residues 1066 to 1181 [ nsp3(1066-1181)] and the constructs nsp3(1066-1203) and nsp3(1035-1181). A search of the protein structure database indicates that the globular domain of the NAB represents a new fold, with a parallel four-strand beta-sheet holding two alpha-helices of three and four turns that are oriented antiparallel to the beta-strands. Two antiparallel two-strand beta-sheets and two 3(10)-helices are anchored against the surface of this barrel-like molecular core. Chemical shift changes upon the addition of single-stranded RNAs (ssRNAs) identified a group of residues that form a positively charged patch on the protein surface as the binding site responsible for the previously reported affinity for nucleic acids. This binding site is similar to the ssRNA-binding site of the sterile alpha motif domain of the Saccharomyces cerevisiae Vts1p protein, although the two proteins do not share a common globular fold.

Nuclear magnetic resonance structure of the N-terminal domain of nonstructural protein 3 from the severe acute respiratory syndrome coronavirus

This paper describes the structure determination of nsp3a, the N-terminal domain of the severe acute respiratory syndrome coronavirus (SARS-CoV) nonstructural protein 3. nsp3a exhibits a ubiquitin-like globular fold of residues 1 to 112 and a flexibly extended glutamic acid-rich domain of residues 113 to 183. In addition to the four beta-strands and two alpha-helices that are common to ubiquitin-like folds, the globular domain of nsp3a contains two short helices representing a feature that has not previously been observed in these proteins. Nuclear magnetic resonance chemical shift perturbations showed that these unique structural elements are involved in interactions with single-stranded RNA. Structural similarities with proteins involved in various cell-signaling pathways indicate possible roles of nsp3a in viral infection and persistence.

The structure of echovirus type 12 bound to a two-domain fragment of its cellular attachment protein decay-accelerating factor (CD 55)

Echovirus type 12 (EV12), an enterovirus of the Picornaviridae family, uses the complement regulator, decay-accelerating factor (DAF, CD55) as a cellular receptor. We have calculated a three-dimensional reconstruction of EV12 bound to a fragment of DAF, consisting of short consensus repeat domains 3 and 4, from cryo-negative stain electron microscopy data (EMD #1057). This shows that, as for an earlier reconstruction of the related echovirus type 7 bound to DAF, attachment is not within the viral canyon but occurs close to the two-fold symmetry axes. Despite this general similarity, our reconstruction reveals a receptor interaction that is quite different from that observed for EV7. Fitting of the crystallographic co-ordinates for DAF34 and EV11 into the reconstruction shows a close agreement between the crystal structure of the receptor fragment and the density for the virus-bound receptor, allowing unambiguous positioning of the receptor with respect to the virion (PDB #1UPN). Our finding that the mode of virus-receptor interaction in EV12 is distinct from that seen for EV7 raises interesting questions regarding the evolution and biological significance of the DAF-binding phenotype in these viruses.

Structure-based energetics of protein interfaces guides foot-and-mouth disease virus vaccine design

Virus capsids are primed for disassembly, yet capsid integrity is key to generating a protective immune response. Foot-and-mouth disease virus (FMDV) capsids comprise identical pentameric protein subunits held together by tenuous noncovalent interactions and are often unstable. Chemically inactivated or recombinant empty capsids, which could form the basis of future vaccines, are even less stable than live virus. Here we devised a computational method to assess the relative stability of protein-protein interfaces and used it to design improved candidate vaccines for two poorly stable, but globally important, serotypes of FMDV: O and SAT2. We used a restrained molecular dynamics strategy to rank mutations predicted to strengthen the pentamer interfaces and applied the results to produce stabilized capsids. Structural analyses and stability assays confirmed the predictions, and vaccinated animals generated improved neutralizing-antibody responses to stabilized particles compared to parental viruses and wild-type capsids.

Activation of c-Jun N-terminal kinases and p38-mitogen-activated protein kinases in human heart failure secondary to ischaemic heart disease.

Three well-characterized mitogen-activated protein kinase (MAPK) subfamilies are expressed in rodent and rabbit hearts, and are activated by pathophysiological stimuli. We have determined and compared the expression and activation of these MAPKs in donor and failing human hearts. The amount and activation of MAPKs was assessed in samples from the left ventricles of 4 unused donor hearts and 12 explanted hearts from patients with heart failure secondary to ischaemic heart disease. Total MAPKs or dually phosphorylated (activated) MAPKs were detected by Western blotting and MAPK activities were measured by in gel kinase assays. As in rat heart, c-Jun N-terminal kinases (JNKs) were detected in human hearts as bands corresponding to 46 and 54 kDa; p38-MAPK(s) was detected as a band corresponding to approximately 40 kDa, and extracellularly regulated kinases, ERK1 and ERK2, were detected as 44- and 42-kDa bands respectively. The total amounts of 54 kDa JNK, p38-MAPK and ERK2 were similar in all samples, although 46-kDa JNK was reduced in the failing hearts. However, the mean activities of JNKs and p38-MAPK(s) were significantly higher in failing heart samples than in those from donor hearts (P<0.05). There was no significant difference in phosphorylated (activated) ERKs between the two groups. In conclusion, JNKs, p38-MAPK(s) and ERKs are expressed in the human heart and the activities of JNKs and p38-MAPK(s) were increased in heart failure secondary to ischaemic heart disease. These data indicate that JNKs and p38-MAPKs may be important in human cardiac pathology.

Crystallographic structure and substrate-binding interactions of the molybdate-binding protein of the phytopathogen Xanthomonas axonopodis pv. citri

In Xanthomonas axonopodis pv. citri (Xac or X citri), the modA gene codes for a periplasmic protein (ModA) that is capable of binding molybdate and tungstate as part of the ABC-type transporter required for the uptake of micronutrients. In this study, we report the crystallographic structure of the Xac ModA protein with bound molybdate. The Xac ModA structure is similar to orthologs with known three-dimensional structures and consists of two nearly symmetrical domains separated by a hinge region where the oxyanion-binding site lies. Phylogenetic analysis of different ModA orthologs based on sequence alignments revealed three groups of molybdate-binding proteins: bacterial phytopathogens, enterobacteria and soil bacteria. Even though the ModA orthologs are segregated into different groups, the ligand-binding hydrogen bonds are mostly conserved, except for Archaeglobus fulgidus ModA. A detailed discussion of hydrophobic interactions in the active site is presented and two new residues, Ala(38) and Ser(151), are shown to be part of the ligand-binding pocket. (c) 2007 Elsevier B.V All rights reserved.

The crystal structure of the leptospiral hypothetical protein LIC12922 reveals homology with the periplasmic chaperone SurA

Leptospirosis is a world spread zoonosis caused by members of the genus Leptospira. Although leptospires were identified as the causal agent of leptospirosis almost 100 years ago, little is known about their biology, which hinders the development of new treatment and prevention strategies. One of the several aspects of the leptospiral biology not yet elucidated is the process by which outer membrane proteins (OMPs) traverse the periplasm and are inserted into the outer membrane. The crystal structure determination of the conserved hypothetical protein LIC12922 from Leptospira interrogans revealed a two domain protein homologous to the Escherichia coli periplasmic chaperone SurA. The LIC12922 NC-domain is structurally related to the chaperone modules of E. coli SurA and trigger factor, whereas the parvulin domain is devoid of peptidyl prolyl cis-trans isomerase activity. Phylogenetic analyses suggest a relationship between LIC12922 and the chaperones PrsA, PpiD and SurA. Based on our structural and evolutionary analyses, we postulate that LIC12922 is a periplasmic chaperone involved in OMPs biogenesis in Leptospira spp. Since LIC12922 homologs were identified in all spirochetal genomes sequenced to date, this assumption may have implications for the OMPs biogenesis studies not only in leptospires but in the entire Phylum Spirochaetes. (C) 2010 Elsevier Inc. All rights reserved.

Structure, processing and midgut secretion of putative peritrophic membrane ancillary protein (PMAP) from Tenebrio molitor larvae

A cDNA coding for a Tenebrio molitor midgut protein named peritrophic membrane ancillary protein (PMAP) was cloned and sequenced. The complete cDNA codes for a protein of 595 amino acids with six insect-allergen-related-repeats that may be grouped in A (predicted globular)- and B (predicted nonglobular)-types forming an ABABAB structure. The PMAP-cDNA was expressed in Pichia pastoris and the recombinant protein (64 kDa) was purified to homogeneity and used to raise antibodies in rabbits. The specific antibody detected PMAP peptides (22 kDa) in the anterior and middle midgut tissue, luminal contents, peritrophic membrane and feces. These peptides derive from PMAP, as supported by mass spectrometry, and resemble those formed by the in vitro action of trypsin on recombinant PMAP. Both in vitro and in vivo PMAP processing seem to occur by attack of trypsin to susceptible bonds in the coils predicted to link AB pairs, thus releasing the putative functional AB structures. The AB-domain structure of PMAP is found in homologous proteins from several insect orders, except lepidopterans that have the apparently derived protein known as nitrile-specifier protein. Immunocytolocalization shows that PMAP is secreted by exocytosis and becomes entrapped in the glycocalyx, before being released into midgut contents. Circumstantial evidence suggests that PMAP-like proteins have a role in peritrophic membrane type 2 formation. (C) 2007 Elsevier Ltd. All rights reserved.