Evolutionary dynamics of Chlamydia trachomatis genome and identification of molecular patterns of hypothetical protein coding genes

The obligate intracellular bacterium Chlamydia trachomatis is a human pathogen of major public health significance. Strains can be classified into 15 main serovars (A to L3) that preferentially cause ocular infections (A-C), genital infections (D-K) or lymphogranuloma venereum (LGV) (L1-L3), but the molecular basis behind their distinct tropism, ecological success and pathogenicity is not welldefined. Most chlamydial research demands culture in eukaryotic cell lines, but it is not known if stains become laboratory adapted. By essentially using genomics and transcriptomics, we aimed to investigate the evolutionary patterns underlying the adaptation of C. trachomatis to the different human tissues, given emphasis to the identification of molecular patterns of genes encoding hypothetical proteins, and to understand the adaptive process behind the C. trachomatis in vivo to in vitro transition. Our results highlight a positive selection-driven evolution of C. trachomatis towards nichespecific adaptation, essentially targeting host-interacting proteins, namely effectors and inclusion membrane proteins, where some of them also displayed niche-specific expression patterns. We also identified potential "ocular-specific" pseudogenes, and pointed out the major gene targets of adaptive mutations associated with LGV infections. We further observed that the in vivo-derived genetic makeup of C. trachomatis is not significantly compromised by its long-term laboratory propagation. In opposition, its introduction in vitro has the potential to affect the phenotype, likely yielding virulence attenuation. In fact, we observed a "genital-specific" rampant inactivation of the virulence gene CT135, which may impact the interpretation of data derived from studies requiring culture. Globally, the findings presented in this Ph.D. thesis contribute for the understanding of C.trachomatis adaptive evolution and provides new insights into the biological role of C. trachomatishypothetical proteins. They also launch research questions for future functional studies aiming toclarify the determinants of tissue tropism, virulence or pathogenic dissimilarities among C. trachomatisstrains.

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-based assignment of the biochemical function of a hypothetical protein: A test case of structural genomics

Many small bacterial, archaebacterial, and eukaryotic genomes have been sequenced, and the larger eukaryotic genomes are predicted to be completely sequenced within the next decade. In all genomes sequenced to date, a large portion of these organisms’ predicted protein coding regions encode polypeptides of unknown biochemical, biophysical, and/or cellular functions. Three-dimensional structures of these proteins may suggest biochemical or biophysical functions. Here we report the crystal structure of one such protein, MJ0577, from a hyperthermophile, Methanococcus jannaschii, at 1.7-Å resolution. The structure contains a bound ATP, suggesting MJ0577 is an ATPase or an ATP-mediated molecular switch, which we confirm by biochemical experiments. Furthermore, the structure reveals different ATP binding motifs that are shared among many homologous hypothetical proteins in this family. This result indicates that structure-based assignment of molecular function is a viable approach for the large-scale biochemical assignment of proteins and for discovering new motifs, a basic premise of structural genomics.

trEST, trGEN and Hits: access to databases of predicted protein sequences.

High throughput genome (HTG) and expressed sequence tag (EST) sequences are currently the most abundant nucleotide sequence classes in the public database. The large volume, high degree of fragmentation and lack of gene structure annotations prevent efficient and effective searches of HTG and EST data for protein sequence homologies by standard search methods. Here, we briefly describe three newly developed resources that should make discovery of interesting genes in these sequence classes easier in the future, especially to biologists not having access to a powerful local bioinformatics environment. trEST and trGEN are regularly regenerated databases of hypothetical protein sequences predicted from EST and HTG sequences, respectively. Hits is a web-based data retrieval and analysis system providing access to precomputed matches between protein sequences (including sequences from trEST and trGEN) and patterns and profiles from Prosite and Pfam. The three resources can be accessed via the Hits home page (http://hits. isb-sib.ch).

Vaccine potentials of an intrinsically unstructured fragment derived from the blood stage-associated Plasmodium falciparum protein PFF0165c.

We have identified new malaria vaccine candidates through the combination of bioinformatics prediction of stable protein domains in the Plasmodium falciparum genome, chemical synthesis of polypeptides, in vitro biological functional assays, and association of an antigen-specific antibody response with protection against clinical malaria. Within the predicted open reading frame of P. falciparum hypothetical protein PFF0165c, several segments with low hydrophobic amino acid content, which are likely to be intrinsically unstructured, were identified. The synthetic peptide corresponding to one such segment (P27A) was well recognized by sera and peripheral blood mononuclear cells of adults living in different regions where malaria is endemic. High antibody titers were induced in different strains of mice and in rabbits immunized with the polypeptide formulated with different adjuvants. These antibodies recognized native epitopes in P. falciparum-infected erythrocytes, formed distinct bands in Western blots, and were inhibitory in an in vitro antibody-dependent cellular inhibition parasite-growth assay. The immunological properties of P27A, together with its low polymorphism and association with clinical protection from malaria in humans, warrant its further development as a malaria vaccine candidate.

trome, trEST and trGEN: databases of predicted protein sequences.

We previously introduced two new protein databases (trEST and trGEN) of hypothetical protein sequences predicted from EST and HTG sequences, respectively. Here, we present the updates made on these two databases plus a new database (trome), which uses alignments of EST data to HTG or full genomes to generate virtual transcripts and coding sequences. This new database is of higher quality and since it contains the information in a much denser format it is of much smaller size. These new databases are in a Swiss-Prot-like format and are updated on a weekly basis (trEST and trGEN) or every 3 months (trome). They can be downloaded by anonymous ftp from ftp://ftp.isrec.isb-sib.ch/pub/databases.

Assignment of Putative Functions to Membrane "Hypothetical Proteins" from the Trypanosoma cruzi Genome

Protozoan parasites cause thousands of deaths each year in developing countries. The genome projects of these parasites opened a new era in the identification of therapeutic targets. However, the putative function could be predicted for fewer than half of the protein-coding genes. In this work, all Trypanosoma cruzi proteins containing predicted transmembrane spans were processed through an automated computational routine and further analyzed in order to assign the most probable function. The analysis consisted of dissecting the whole predicted protein in different regions. More than 5,000 sequences were processed, and the predicted biological functions were grouped into 19 categories according to the hits obtained after analysis. One focus of interest, due to the scarce information available on trypanosomatids, is the proteins involved in signal-transduction processes. In the present work, we identified 54 proteins belonging to this group, which were individually analyzed. The results show that by means of a simple pipeline it was possible to attribute probable functions to sequences annotated as coding for "hypothetical proteins.'' Also, we successfully identified the majority of candidates participating in the signal-transduction pathways in T. cruzi.

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.

trEST, trGEN and Hits: access to databases of predicted protein sequences

High throughput genome (HTG) and expressed sequence tag (EST) sequences are currently the most abundant nucleotide sequence classes in the public database. The large volume, high degree of fragmentation and lack of gene structure annotations prevent efficient and effective searches of HTG and EST data for protein sequence homologies by standard search methods. Here, we briefly describe three newly developed resources that should make discovery of interesting genes in these sequence classes easier in the future, especially to biologists not having access to a powerful local bioinformatics environment. trEST and trGEN are regularly regenerated databases of hypothetical protein sequences predicted from EST and HTG sequences, respectively. Hits is a web-based data retrieval and analysis system providing access to precomputed matches between protein sequences (including sequences from trEST and trGEN) and patterns and profiles from Prosite and Pfam. The three resources can be accessed via the Hits home page (http://hits.isb-sib.ch).

Cytokine-related genes identified from the RIKEN full-length mouse cDNA data set

To identify novel cytokine-related genes, we searched the set of 60,770 annotated RIKEN mouse cDNA clones (FANTOM2 clones), using keywords such as cytokine itself or cytokine names (such as interferon, interleukin, epidermal growth factor, fibroblast growth factor, and transforming growth factor). This search produced 108 known cytokines and cytokine-related products such as cytokine receptors, cytokine-associated genes, or their products (enhancers, accessory proteins, cytokine-induced genes). We found 15 clusters of FANTOM2 clones that are candidates for novel cytokine-related genes. These encoded products with strong sequence similarity to guanylate-binding protein (GBP-5), interleukin-1 receptor-associated kinase 2 (IRAK-2), interleukin 20 receptor alpha isoform 3, a member of the interferon-inducible proteins of the Ifi 200 cluster, four members of the membrane-associated family 1-8 of interferon-inducible proteins, one p27-like protein, and a hypothetical protein containing a Toll/Interleukin receptor domain. All four clones representing novel candidates of gene products from the family contain a novel highly conserved cross-species domain. Clones similar to growth factor-related products included transforming growth factor beta-inducible early growth response protein 2 (TIEG-2), TGFbeta-induced factor 2, integrin beta-like 1, latent TGF-binding protein 4S, and FGF receptor 4B. We performed a detailed sequence analysis of the candidate novel genes to elucidate their likely functional properties.

Karyotype variability in KP1(+) and KP1(-) strains of Trypanosoma rangeli isolated in Brazil and Colombia

In the present study, the molecular karyotypes of 12 KP1(+) and KP1(-) Trypanosoma rangeli strains were determined and 10 different molecular markers were hybridized to the chromosomes of the parasite, including seven obtained from T. rangeli [ubiquitin hydrolase (UH), a predicted serine/threonine protein kinase (STK), hexose transporter, hypothetical protein, three anonymous sequences] and three from Trypanosoma cruzi [ubiquitin-conjugating enzyme E2 (UBE2), ribosomal RNA methyltransferase (rRNAmtr), proteasome non-ATPase regulatory subunit 6 (PSMD6)]. Despite intraspecific variation, analysis of the karyotype profiles permitted the division of the T rangeli strains into two groups coinciding with the KP1(+) and KP1(-) genotypes. Southern blot hybridization showed that, except for the hexose transporter probe, all other probes produced distinct patterns able to differentiate the KP1(+) and KP1(-) genotypes. The UH, STK and An-1A04 probes exclusively hybridized to the chromosomes of KP1(+) strains and can be used as markers of this group. In addition, the UBE2, rRNAmtr and PSMD6 markers, which are present in a conserved region in all trypanosomatid species sequenced so far, co-hybridized to the same T. rangeli chromosomal bands, suggesting the occurrence of gene synteny in these species. The finding of distinct molecular karyotypes in KP1(+) and KP1 (-) strains of T rangeli is noteworthy and might be used as a new approach to the study of genetic variability in this parasite. Together with the Southern blot hybridization results, these findings demonstrate that differences at the kDNA level might be associated with variations in nuclear DNA. (c) 2009 Elsevier BY. All rights reserved.

Identification of Schistosoma mansoni candidate antigens for diagnosis of schistosomiasis

The development of a more sensitive diagnostic test for schistosomiasis is needed to overcome the limitations of the use of stool examination in low endemic areas. Using parasite antigens in enzyme linked immunosorbent assay is a promising strategy, however a more rational selection of parasite antigens is necessary. In this study we performed in silico analysis of the Schistosoma mansoni genome, using SchistoDB database and bioinformatic tools for screening immunogenic antigens. Based on evidence of expression in all parasite life stage within the definitive host, extracellular or plasmatic membrane localization, low similarity to human and other helminthic proteins and presence of predicted B cell epitopes, six candidates were selected: a glycosylphosphatidylinositol-anchored 200 kDa protein, two putative cytochrome oxidase subunits, two expressed proteins and one hypothetical protein. The recognition in unidimensional and bidimensional Western blot of protein with similar molecular weight and isoelectric point to the selected antigens by sera from S. mansoni infected mice indicate a good correlation between these two approaches in selecting immunogenic proteins.

Longitudinal analyses of immune responses to Plasmodium falciparum derived peptides corresponding to novel blood stage antigens in coastal Kenya

We have recently described 95 predicted alpha-helical coiled-coil peptides derived from putative Plasmodium falciparum erythrocytic stage proteins. Seventy peptides recognized with the highest level of prevalence by sera from three endemic areas were selected for further studies. In this study, we sequentially examined antibody responses to these synthetic peptides in two cohorts of children at risk of clinical malaria in Kilifi district in coastal Kenya, in order to characterize the level of peptide recognition by age, and the role of anti-peptide antibodies in protection from clinical malaria. Antibody levels from 268 children in the first cohort (Chonyi) were assayed against 70 peptides. Thirty-nine peptides were selected for further study in a second cohort (Junju). The rationale for the second cohort was to confirm those peptides identified as protective in the first cohort. The Junju cohort comprised of children aged 1-6 years old (inclusive). Children were actively followed up to identify episodes of febrile malaria in both cohorts. Of the 70 peptides examined, 32 showed significantly (p<0.05) increased antibody recognition in older children and 40 showed significantly increased antibody recognition in parasitaemic children. Ten peptides were associated with a significantly reduced odds ratio (OR) for an episode of clinical malaria in the first cohort of children and two of these peptides (LR146 and AS202.11) were associated with a significantly reduced OR in both cohorts. LR146 is derived from hypothetical protein PFB0145c in PlasmoDB. Previous work has identified this protein as a target of antibodies effective in antibody dependent cellular inhibition (ADCI). The current study substantiates further the potential of protein PFB0145c and also identifies protein PF11_0424 as another likely target of protective antibodies against P. falciparum malaria

Estudi d'un receptor quinasa atípic (Mark) i de les proteïnes que interaccionen amb el seu domini intracel·lular. Transducció del senyal i desenvolupament en blat de moro (Zea mays)

Aquesta tesi doctoral s'engloba dins d'un projecte general d'estudi de gens implicats en l'embriogènesi del blat de moro. L'embriogènesi del blat de moro, i en general la de totes les plantes superiors, es dóna en tres etapes: una primera etapa on es diferencien tots els diversos teixits que formaran l'embrió, una segona etapa on l'embrió acumula productes de reserva i un tercer període, la dormància, que finalitza quan les condicions ambientals són les idònies per a la germinació. En el laboratori estàvem interessats, concretament, en l'estudi de gens implicats en la primera etapa morfogenètica, on els diferents teixits i estructures embrionàries queden definides. Per tal d'estudiar gens que s'expressaven en aquest període, una de les estratègies que es va realitzar fou un crivellat diferencial entre teixit embrionari i teixit de planta adulta. D'entre els diferents clons obtinguts, un corresponia a un clon parcial que presentava similitud amb receptors quinasa i que fou objecte d'estudi. A partir d'aquest clon es va obtenir el clon complet i es va anomenar MARK (per Maize Atypical Receptor Kinase). MARK presenta una estructura típica d'un receptor quinasa amb un domini extracel.lular, que conté 6 còpies imperfectes de LRR (Leucine- Rich Repeats), un únic domini transmembrana i un domini quinasa intracel.lular. El domini quinasa de MARK presenta, però, algunes variacions en els residus aminoacídics que es consideren claus per a la funció catalítica dels dominis quinasa. En concret cinc dels aminoàcids considerats essencials per a la fosforilació es troben substituits en el domini quinasa de MARK (DK-MARK). Els experiments de fosforilació in vitro que es van realitzar al laboratori, van mostrar com MARK era incapaç de fosforilar in vitro. Aquesta característica no és, però, exclusiva de MARK. Una búsqueda en les bases de dades ens van permetre identificar altres seqüències que també presentaven els mateixos o altres canvis en aquestes posicions aminoacídiques. En les bases de dades de plantes es van identificar un conjunt de seqüències genòmiques o ESTs amb aquestes característiques i només una d'elles, la proteïna TMKL1 d'Arabidopsis, ha sigut descrita com un receptor quinasa incapaç de fosforilar in vitro. Respecte a la búsqueda de receptors similars a MARK en les bases de dades d'animals, es van identificar també un conjunt de proteïnes que, en alguns casos, s'ha descrit que no tenen activitat quinasa in vivo. Per exemple, un dels casos més ben estudiats és el del receptor erbB3 que forma part de la família de receptors del EGF (Epidermal Growth Factor). Aquesta família de receptors està formada per 4 receptors: erbB1, erbB2, erbB3 i erbB4, dels quals només l'erbB3 no presenta activitat catalítica. S'ha descrit que erbB3 és capaç, tot i no fosforilar in vivo, de participar activament en la transducció del senyal formant heterodímers amb els altres membres de la família. Així, erbB3 és fosforilat pel seu partner i pot iniciar la cascada de transducció del senyal. La participació d'erbB3 en la transducció del senyal és essencial ja que embrions de ratolí knock-out pel gen erbB3 són inviables. Així doncs, el fet que receptors quinasa catalíticament inactius participin en les cascades de transducció del senyal, suggereix l'existència de nous mecanismes d'acció per a la transducció del senyal. Per tant, l'objectiu d'aquest treball fou l'estudi del mecanisme d'acció de MARK mitjançant la caracterització les proteïnes capaces d'interaccionar amb el seu domini quinasa. Per tal d'assolir aquest objectiu, es va realitzar un crivellat de doble-híbrid amb una llibreria de cDNA d'embrions de blat de moro de 7 DAP. D'aquest crivellat es va obtenir un conjunt de possibles clons positius que foren seqüenciats i entre els quals es van escollir per un estudi més detallat aquells que s'havien obtingut més vegades com a clons independents. Aquests clons codificaven per: una SAMDC (S-Adenosil Descarboxilasa), una eIF5 (Eukaryotic translation initiation), una hypothetical protein, una unknown protein, una gamma-adaptina i una MAP4K. Amb aquests 6 clons es van fer estudis in vitro i in vivo per tal de confirmar al seva interacció amb DK-MARK. Els estudis in vivo es van realitzar amb la soca de llevat AH109, una soca més astringent que la utilitzada en el crivellat, ja que presenta tres gens marcadors: Histidina, Adenina i Lacz. Els resultats obtinguts van mostrar que els clons codificants per SAMDC i eIF5 no van créixer en un medi selectiu per His i Ade i, per tant o es tracta de falsos positius del sistema o la seva interacció amb DK-MARK és dèbil. D'altra banda, la resta dels clons analitzats (proteïna hipotètica, una proteïna de funció desconeguda, la gamma-adaptina i una MAP4K) van créixer en medis en absència de Histidina i Adenina. Els assatjos de b-galactosidasa van ser tots positius a excepció de la proteïna hipotètica suggerint que potser aquesta interacció sigui més feble. D'altra banda també es van realitzar estudis in vitro amb la tècnica del pull-down. Els resultats obtinguts amb aquesta tècnica van recolzar els obtinguts en cèl.lules de llevat, ja que tots els clons analitzats a excepció dels codificants per SAMDC i eIF5 van donar un resultat d'interacció amb KD-MARK in vitro positiu. Davant aquests resultats ens vam centrar en l'estudi de la proteïna similar a MAP4K, doncs algunes proteïnes de la seva família s'han relacionat amb receptors de membrana. Els clons que es va obtenir del crivellat codificaven per una proteïna similar amb el domini C-terminal a les proteïnes BnMAP4Ka1 i a2 de Brassica napus. Aquestes proteïnes presenten una forta similitud de seqüència amb proteïnes de la família GCK/SPS1 que formen part d'un grup particular de MAPK relacionades amb la proteïna Ste20 (sterile 20 protein) de llevat. Ste20p activa la MAP3K de llevat Ste11 directament per fosforilació, transduint d'aquesta manera el senyal del receptor de feromones de creuament de les cèl.lules de llevat i es pot, doncs, considerar com una proteïna del tipus MAP4K (mitogen-activated protein kinase kinase kinase kinase). En els darrers anys, s'han identificat un gran nombre de proteïnes similars a Ste20: fins a una trentena en mamífers, en Drosophila, en Caenorhabditis elegans i en altres organismes. Segons la seva estructura aminoacídica, la família Ste20 s'ha classificat en dues subfamílies: les proteïnes STE20/PAK (p21-activated kinases) i la subfamília GCK/SPS1 (germinal center kinases). Les dues subfamílies estan formades per proteïnes que contenen un domini quinasa i un domini regulador, però, mentre que les proteïnes PAK presenten el domini quinasa en la part C-terminal, les GCKs el presenten en la regió N terminal. Les proteïnes GCK presenten una elevada diversitat estructural en el domini regulador permetent la seva classificació en 6 subfamílies. Mitjançant la tècnica del RACE es va obtenir el clon de cDNA complet que es va anomenar MIK (MARK Interacting Kinase). Amb la tècnica del Southern blot es va poder determinar que el gen MIK és un gen de còpia única en el genoma de blat de moro. Per tal d'analitzar la possible interacció entre DK-MARK i MIK, es va estudiar tant el patró d'expressió d'ambdós gens com el seu patró d'acumulació d'ambdues proteïnes durant l'embriogènesi del blat de moro. El patró d'expressió, analitzat per Northen blot va mostrar uns patrons coincidents al llarg de l'embriogènesi des del seu inici fins als 20 DAP amb una acumulació màxima de mRNA en embrions de 15 DAP. D'altra banda per tal d'estudiar el patró d'acumulació de la proteïna MIK així com per comparar-lo amb el de MARK, es van realitzar estudis de Westerns blot. Els resultats també van mostrar una coincidència en el temps de l'acumulació de les proteïnes MARK i MIK durant l'embriogènesi de blat de moro amb una major acumulació en embrions de 15 i 20 DAP. Es van dur a terme també estudis d'immunolocalitzacions sobre embrions de blat de moro de 15 DAP per tal d'estudiar en quins teixits s'acumulaven ambdues proteïnes. Les immunolocalitzacions van mostrar una major acumulació tant de MARK com de MIK en les zones meristemàtiques i en el teixit vascular sobretot del coleòptil on s'aprecia una forta co-localització de MARK i MIK. Totes aquestes dades són compatibles, doncs, amb una possible interacció de les proteïnes MARK i MIK, tot i que no la demostren. Per tal de demostrar la interacció es van realitzar experiments d'immunoprecipitació in vivo a partir d'extractes d'embrions. Malauradament, els resultats no són clars i en aquests moments en el laboratori s'estan posant a punt aquests experiments. També es van realitzar estudis comparatius de seqüència amb diferents proteïnes de la família GCK, mostrant una major similitud amb les proteïnes de la subfamília GCK-III. La subfamília GCK-III ha estat molt poc estudiada i en formen part un conjunt de proteïnes amb funcions molt diverses des de l'apoptosi, la citoquinesi o l'anòxia cel.lular. Per tant, la similitud de seqüència possiblement fa referència a una conservació en el mecanisme d'acció més que no pas a una conservació funcional. La possible interacció de MARK amb el domini C-terminal de MIK (el domini regulador) podria activar aquesta última iniciant una cascada de transducció del senyal en un model en el que una proteïna del tipus GCK-III faria de lligam directa entre un receptor de membrana i una cascada de senyalització intracel.lular. Aquest tipus de lligam entre un recepctor de membrana i mòduls intracel.lulars de senyalització s'ha descrit per a altres proteïnes GCK, si bé no directament sinó a través de proteïnes adaptadores. D'altra banda, la interacció directa de MARK, un receptor quinasa atípic que no té activitat catalítica, amb MIK suggereix un mecanisme on receptors atípics podrien interaccionar en la transducció del senyal activant la via de les MAPK.

Predicting metal-binding site residues in low-resolution structural models

The accurate prediction of the biochemical function of a protein is becoming increasingly important, given the unprecedented growth of both structural and sequence databanks. Consequently, computational methods are required to analyse such data in an automated manner to ensure genomes are annotated accurately. Protein structure prediction methods, for example, are capable of generating approximate structural models on a genome-wide scale. However, the detection of functionally important regions in such crude models, as well as structural genomics targets, remains an extremely important problem. The method described in the current study, MetSite, represents a fully automatic approach for the detection of metal-binding residue clusters applicable to protein models of moderate quality. The method involves using sequence profile information in combination with approximate structural data. Several neural network classifiers are shown to be able to distinguish metal sites from non-sites with a mean accuracy of 94.5%. The method was demonstrated to identify metal-binding sites correctly in LiveBench targets where no obvious metal-binding sequence motifs were detectable using InterPro. Accurate detection of metal sites was shown to be feasible for low-resolution predicted structures generated using mGenTHREADER where no side-chain information was available. High-scoring predictions were observed for a recently solved hypothetical protein from Haemophilus influenzae, indicating a putative metal-binding site.