The Evolution of Function in the Rab family of Small GTPases

Dissertation presented to obtain the PhD degree in Computational Biology.

The InterPro database, an integrated documentation resource for protein families, domains and functional sites.

Signature databases are vital tools for identifying distant relationships in novel sequences and hence for inferring protein function. InterPro is an integrated documentation resource for protein families, domains and functional sites, which amalgamates the efforts of the PROSITE, PRINTS, Pfam and ProDom database projects. Each InterPro entry includes a functional description, annotation, literature references and links back to the relevant member database(s). Release 2.0 of InterPro (October 2000) contains over 3000 entries, representing families, domains, repeats and sites of post-translational modification encoded by a total of 6804 different regular expressions, profiles, fingerprints and Hidden Markov Models. Each InterPro entry lists all the matches against SWISS-PROT and TrEMBL (more than 1,000,000 hits from 462,500 proteins in SWISS-PROT and TrEMBL). The database is accessible for text- and sequence-based searches at http://www.ebi.ac.uk/interpro/. Questions can be emailed to interhelp@ebi.ac.uk.

Evidence for an episodic model of protein sequence evolution.

The evolution of protein function appears to involve alternating periods of conservative evolution and of relatively rapid change. Evidence for such episodic evolution, consistent with some theoretical expectations, comes from the application of increasingly sophisticated models of evolution to large sequence datasets. We present here some of the recent methods to detect functional shifts, using amino acid or codon models. Both provide evidence for punctual shifts in patterns of amino acid conservation, including the fixation of key changes by positive selection. Although a link to gene duplication, a presumed source of functional changes, has been difficult to establish, this episodic model appears to apply to a wide variety of proteins and organisms.

Hepatitis C virus proteins: from structure to function.

Great progress has been made over the past years in elucidating the structure and function of the hepatitis C virus (HCV) proteins, most of which are now actively being pursued as antiviral targets. The structural proteins, which form the viral particle, include the core protein and the envelope glycoproteins E1 and E2. The nonstructural proteins include the p7 viroporin, the NS2 protease, the NS3-4A complex harboring protease and NTPase/RNA helicase activities, the NS4B and NS5A proteins, and the NS5B RNA-dependent RNA polymerase. NS4B is a master organizer of replication complex formation while NS5A is a zinc-containing phosphoprotein involved in the regulation of HCV RNA replication versus particle production. Core to NS2 make up the assembly module while NS3 to NS5B represent the replication module (replicase). However, HCV proteins exert multiple functions during the viral life cycle, and these may be governed by different structural conformations and/or interactions with viral and/or cellular partners. Remarkably, each viral protein is anchored to intracellular membranes via specific determinants that are essential to protein function in the cell. This review summarizes current knowledge of the structure and function of the HCV proteins and highlights recent advances in the field.

Resolving the ortholog conjecture: orthologs tend to be weakly, but significantly, more similar in function than paralogs.

The function of most proteins is not determined experimentally, but is extrapolated from homologs. According to the "ortholog conjecture", or standard model of phylogenomics, protein function changes rapidly after duplication, leading to paralogs with different functions, while orthologs retain the ancestral function. We report here that a comparison of experimentally supported functional annotations among homologs from 13 genomes mostly supports this model. We show that to analyze GO annotation effectively, several confounding factors need to be controlled: authorship bias, variation of GO term frequency among species, variation of background similarity among species pairs, and propagated annotation bias. After controlling for these biases, we observe that orthologs have generally more similar functional annotations than paralogs. This is especially strong for sub-cellular localization. We observe only a weak decrease in functional similarity with increasing sequence divergence. These findings hold over a large diversity of species; notably orthologs from model organisms such as E. coli, yeast or mouse have conserved function with human proteins.

Proteomic analysis of membrane rafts of melanoma cells identifies protein patterns characteristic of the tumor progression stage.

The molecular mechanisms controlling the progression of melanoma from a localized tumor to an invasive and metastatic disease are poorly understood. In the attempt to start defining a functional protein profile of melanoma progression, we have analyzed by LC-MS/MS the proteins associated with detergent resistant membranes (DRMs), which are enriched in cholesterol/sphingolipids-containing membrane rafts, of melanoma cell lines derived from tumors at different stages of progression. Since membrane rafts are involved in several biological processes, including signal transduction and protein trafficking, we hypothesized that the association of proteins with rafts can be regulated during melanoma development and affect protein function and disease progression. We have identified a total of 177 proteins in the DRMs of the cell lines examined. Among these, we have found groups of proteins preferentially associated with DRMs of either less malignant radial growth phase/vertical growth phase (VGP) cells, or aggressive VGP and metastatic cells suggesting that melanoma cells with different degrees of malignancy have different DRM profiles. Moreover, some proteins were found in DRMs of only some cell lines despite being expressed at similar levels in all the cell lines examined, suggesting the existence of mechanisms controlling their association with DRMs. We expect that understanding the mechanisms regulating DRM targeting and the activity of the proteins differentially associated with DRMs in relation to cell malignancy will help identify new molecular determinants of melanoma progression.

Structure and function of the cystic fibrosis transmembrane conductance regulator

Cystic fibrosis (CF) is a lethal autosomal recessive genetic disease caused by mutations in the CF transmembrane conductance regulator (CFTR). Mutations in the CFTR gene may result in a defective processing of its protein and alter the function and regulation of this channel. Mutations are associated with different symptoms, including pancreatic insufficiency, bile duct obstruction, infertility in males, high sweat Cl-, intestinal obstruction, nasal polyp formation, chronic sinusitis, mucus dehydration, and chronic Pseudomonas aeruginosa and Staphylococcus aureus lung infection, responsible for 90% of the mortality of CF patients. The gene responsible for the cellular defect in CF was cloned in 1989 and its protein product CFTR is activated by an increase of intracellular cAMP. The CFTR contains two membrane domains, each with six transmembrane domain segments, two nucleotide-binding domains (NBDs), and a cytoplasmic domain. In this review we discuss the studies that have correlated the role of each CFTR domain in the protein function as a chloride channel and as a regulator of the outwardly rectifying Cl- channels (ORCCs).

New insights into small molecules inhibitors and protein-protein interactions of VirB8 : a critical conserved component of the type IV secretion system.

Les systèmes bactériens de sécrétion de type IV (T4SS) sont constitués d’un ensemble de 8 à 12 protéines conservées. Ces dernières sont utilisées lors de la translocation de protéines, la translocation de complexes ADN-protéines mais aussi pour le transport de ces derniers au travers de la membrane cellulaire. Les T4SS, en tant que facteurs de virulence pour beaucoup de pathogènes comme Brucella suis, sont donc d’excellents modèles cibles pour le développement de médicaments d’antivirulence. Ces médicaments, en privant le pathogène de son facteur essentiel de virulence : le T4SS, constituent une alternative ou encore une amélioration des traitements antibiotiques utilisés actuellement. VirB8, un facteur d’assemblage conservé dans le T4SS, forme des dimères qui sont importants pour la fonction des T4SS dans ces pathogènes. De par ses interactions multiples, VirB8 est un excellent modèle pour l’analyse des facteurs d’assemblage mais aussi en tant que cible de médicaments qui empêcheraient son interaction avec d’autres protéines et qui, in fine, désarmeraient les bactéries en les privant de leur fonctions essentielles de virulence. À ce jour, nous savons qu’il existe un équilibre monomère-dimère et un processus d’homodimerization de VirB8 dont l’importance est vitale pour la fonctionnement biologique des T4SSs. En se basant sur des essais quantitatifs d’interaction, nous avons identifié (i) des sites potentiels d’interaction avec d’autres protéines VirB du T4SS mais aussi (ii) isolé des petites molécules inhibitrices afin de tester la fonction protéique de VirB8. Afin de déterminer les acides aminés importants pour l’hétérodimérization de VirB8 avec VirB10, nous avons effectué des expériences de mutagenèse aléatoire, de phage display et d’arrimage moléculaire in silico. Ces expériences ont démontré l’importance de trois acides aminés localisés sur le feuillet β : R160, S162, T164 et I165. Ces derniers seraient importants pour l’association de VirB8 avec VirB10 étant donné que leur mutagenèse entraine une diminution de la formation du complexe VirB8-VirB10. L’objectif actuel de notre projet de recherche est de pouvoir mieux comprendre mais aussi d’évaluer le rôle de VirB8 dans l’assemblage du T4SS. Grace à un méthode de criblage adaptée à partir de la structure de VirB8, nous avons pu identifié une petite molécule inhibitrice BAR-068, qui aurait un rôle prometteur dans l’inhibition du T4SS. Nous avons utilisé la spectroscopie par fluorescence, l’essai à deux hybrides, le cross-linking et la cristallographie afin de déterminer le mécanisme d'interaction existant entre VirB8 et BAR-068. Ces travaux pourraient permettre de nombreuses avancées, notamment en termes de compréhension des mécanismes d’inhibition du T4SS. Notre objectif ultime est de pouvoir caractériser la séquence d’évènements essentiels à l’assemblage et au fonctionnement du T4SS. De manière globale, notre projet de recherche permettrait de révéler les grands principes d’assemblage des protéines membranaires, les processus de sécrétion de protéines chez les bactéries mais aussi de proposer une nouvelle stratégie lors du développement de drogues antimicrobiennes.

The FunFOLD2 server for the prediction of protein-ligand interactions

The FunFOLD2 server is a new independent server that integrates our novel protein–ligand binding site and quality assessment protocols for the prediction of protein function (FN) from sequence via structure. Our guiding principles were, first, to provide a simple unified resource to make our function prediction software easily accessible to all via a simple web interface and, second, to produce integrated output for predictions that can be easily interpreted. The server provides a clean web interface so that results can be viewed on a single page and interpreted by non-experts at a glance. The output for the prediction is an image of the top predicted tertiary structure annotated to indicate putative ligand-binding site residues. The results page also includes a list of the most likely binding site residues and the types of predicted ligands and their frequencies in similar structures. The protein–ligand interactions can also be interactively visualized in 3D using the Jmol plug-in. The raw machine readable data are provided for developers, which comply with the Critical Assessment of Techniques for Protein Structure Prediction data standards for FN predictions. The FunFOLD2 webserver is freely available to all at the following web site: http://www.reading.ac.uk/bioinf/FunFOLD/FunFOLD_form_2_0.html.

Genetic, structural, and molecular insights into the function of ras of complex proteins domains

Ras of complex proteins (ROC) domains were identified in 2003 as GTP binding modules in large multidomain proteins from Dictyostelium discoideum. Research into the function of these domains exploded with their identification in a number of proteins linked to human disease, including leucine-rich repeat kinase 2 (LRRK2) and death-associated protein kinase 1 (DAPK1) in Parkinson’s disease and cancer, respectively. This surge in research has resulted in a growing body of data revealing the role that ROC domains play in regulating protein function and signaling pathways. In this review, recent advances in the structural informa- tion available for proteins containing ROC domains, along with insights into enzymatic function and the integration of ROC domains as molecular switches in a cellular and organismal context, are explored.

Regulation of gene and protein expression in cardiac myocyte hypertrophy and apoptosis

Considerable efforts have been expended in elucidating the inter-cellular and intra-cellular signaling pathways which elicit cardiac myocyte hypertrophy or apoptosis, and in identifying the changes which are associated with the end-stage of the response. The challenge now is to link the two. Although some of the signaling effects will be the acute modulation of existing protein function, long-term effects which bring about and maintain the hypertrophic state or which culminate in cell death are mediated at the level of gene and protein expression. With the advances in micro-array technology and genome sequencing, it is now possible to obtain a picture of the global gene expression profile in myocytes or in whole heart which dictates the proteins which could be made. This is not the final picture since additional regulation at the level of translation modulates the relative proportions of each protein that can be made from the transcriptome. Even here, further regulation of protein stability and turnover means that ultimately it is still necessary to examine the proteome to determine what may cause the functional changes in a cell. Thus, in order to gain a full picture of events which regulate the response and gain some insight into possible points of intervention for therapy, it is necessary to examine gene expression, mRNA translation and protein expression in concert.

In silico identification and characterization of protein-ligand binding sites

Protein–ligand binding site prediction methods aim to predict, from amino acid sequence, protein–ligand interactions, putative ligands, and ligand binding site residues using either sequence information, structural information, or a combination of both. In silico characterization of protein–ligand interactions has become extremely important to help determine a protein’s functionality, as in vivo-based functional elucidation is unable to keep pace with the current growth of sequence databases. Additionally, in vitro biochemical functional elucidation is time-consuming, costly, and may not be feasible for large-scale analysis, such as drug discovery. Thus, in silico prediction of protein–ligand interactions must be utilized to aid in functional elucidation. Here, we briefly discuss protein function prediction, prediction of protein–ligand interactions, the Critical Assessment of Techniques for Protein Structure Prediction (CASP) and the Continuous Automated EvaluatiOn (CAMEO) competitions, along with their role in shaping the field. We also discuss, in detail, our cutting-edge web-server method, FunFOLD for the structurally informed prediction of protein–ligand interactions. Furthermore, we provide a step-by-step guide on using the FunFOLD web server and FunFOLD3 downloadable application, along with some real world examples, where the FunFOLD methods have been used to aid functional elucidation.

Mutational analysis of Arabidopsis thaliana plant uncoupling mitochondrial protein

In this study, point mutations were introduced in plant uncoupling mitochondrial protein AtUCP1, a typical member of the plant uncoupling protein (UCP) gene subfamily, in amino acid residues Lys147, Arg155 and Tyr269, located inside the so-called UCP-signatures, and in two more residues, Cys28 and His83, specific for plant UCPs. The effects of amino acid replacements on AtUCP1 biochemical properties were examined using reconstituted proteoliposomes. Residue Arg155 appears to be crucial for AtUCP1 affinity to linoleic acid (LA) whereas His83 plays an important role in AtUCP1 transport activity. Residues Cys28, Lys147, and also Tyr269 are probably essential for correct protein function, as their substitutions affected either the AtUCP1 affinity to LA and its transport activity, or sensitivity to inhibitors (purine nucleotides). Interestingly, Cys28 substitution reduced ATP inhibitory effect on AtUCP1, while Tyr269Phe mutant exhibited 2.8-fold increase in sensitivity to ATP, in accordance with the reverse mutation Phe267Tyr of mammalian UCP1. (C) 2007 Elsevier B.V. All fights reserved.

COX24 codes for a mitochondrial protein required for processing of the COX1 transcript

In most strains of Saccharomyces cerevisiae the mitochondrial gene COX1, for subunit 1 of cytochrome oxidase, contains multiple exons and introns. Processing of COX1 primary transcript requires accessory proteins factors, some of which are encoded by nuclear genes and others by reading frames residing in some of the introns of the COX1 and COB genes. Here we show that the low molecular weight protein product of open reading frame YLR204W, for which we propose the name COX24, is also involved in processing of COX1 RNA intermediates. The growth defect of cox24 mutants is partially rescued in strains harboring mitochondrial DNA lacking introns. Northern blot analyses of mitochondrial transcripts indicate cox24 null mutants to be blocked in processing of introns aI2 and aI3. The dependence of intron aI3 excision on Cox24p is also supported by the growth properties of the cox24 mutant harboring mitochondrial DNA with different intron compositions. The intermediate phenotype of the cox24 mutant in the background of intronless mitochondrial DNA, however, suggests that in addition to its role in splicing of the COX1 pre-mRNA, Cox24p still has another function. Based on the analysis of a cox14-cox24 double mutant, we propose that the other function of Cox24p is related to translation of the COX1 mRNA. © 2006 by The American Society for Biochemistry and Molecular Biology, Inc.

ZmPUMP encodes a fully functional monocot plant uncoupling mitochondrial protein whose affinity to fatty acid is increased with the introduction of a His pair at the second matrix loop

Uncoupling proteins (UCPs) are specialized mitochondrial transporter proteins that uncouple respiration from ATP synthesis. In this study, cDNA encoding maize uncoupling protein (ZmPUMP) was expressed in Escherichia coli and recombinant ZmPUMP reconstituted in liposomes. ZmPUMP activity was associated with a linoleic acid (LA)-mediated H+ efflux with Km of 56.36 ± 0.27 μM and Vmax of 66.9 μmol H+ min-1 (mg prot)-1. LA-mediated H+ fluxes were sensitive to ATP inhibition with Ki of 2.61 ± 0.36 mM (at pH 7.2), a value similar to those for dicot UCPs. ZmPUMP was also used to investigate the importance of a histidine pair present in the second matrix loop of mammalian UCP1 and absent in plant UCPs. ZmPUMP with introduced His pair (Lys155His and Ala157His) displayed a 1.55-fold increase in LA-affinity while its activity remained unchanged. Our data indicate conserved properties of plant UCPs and suggest an enhancing but not essential role of the histidine pair in proton transport mechanism. © 2006 Elsevier Inc. All rights reserved.