The UlaG protein familly defines novel structural and functional motifs grafted on an ancient RNase fold

Background: Bacterial populations are highly successful at colonizing new habitats and adapting to changing environmental conditions, partly due to their capacity to evolve novel virulence and metabolic pathways in response to stress conditions and to shuffle them by horizontal gene transfer (HGT). A common theme in the evolution of new functions consists of gene duplication followed by functional divergence. UlaG, a unique manganese-dependent metallo-b-lactamase (MBL) enzyme involved in L-ascorbate metabolism by commensal and symbiotic enterobacteria, provides a model for the study of the emergence of new catalytic activities from the modification of an ancient fold. Furthermore, UlaG is the founding member of the so-called UlaG-like (UlaGL) protein family, a recently established and poorly characterized family comprising divalent (and perhaps trivalent)metal-binding MBLs that catalyze transformations on phosphorylated sugars and nucleotides. Results: Here we combined protein structure-guided and sequence-only molecular phylogenetic analyses to dissect the molecular evolution of UlaG and to study its phylogenomic distribution, its relatedness with present-day UlaGL protein sequences and functional conservation. Phylogenetic analyses indicate that UlaGL sequences are present in Bacteria and Archaea, with bona fide orthologs found mainly in mammalian and plant-associated Gramnegative and Gram-positive bacteria. The incongruence between the UlaGL tree and known species trees indicates exchange by HGT and suggests that the UlaGL-encoding genes provided a growth advantage under changing conditions. Our search for more distantly related protein sequences aided by structural homology has uncovered that UlaGL sequences have a common evolutionary origin with present-day RNA processing and metabolizing MBL enzymes widespread in Bacteria, Archaea, and Eukarya. This observation suggests an ancient origin for the UlaGL family within the broader trunk of the MBL superfamily by duplication, neofunctionalization and fixation. Conclusions: Our results suggest that the forerunner of UlaG was present as an RNA metabolizing enzyme in the last common ancestor, and that the modern descendants of that ancestral gene have a wide phylogenetic distribution and functional roles. We propose that the UlaGL family evolved new metabolic roles among bacterial and possibly archeal phyla in the setting of a close association with metazoans, such as in the mammalian gastrointestinal tract or in animal and plant pathogens, as well as in environmental settings. Accordingly, the major evolutionary forces shaping the UlaGL family include vertical inheritance and lineage-specific duplication and acquisition of novel metabolic functions, followed by HGT and numerous lineage-specific gene loss events.

Splitting statistical potentials into meaningful scoring functions: testing the prediction of near-native structures from decoy conformations

Background: Recent advances on high-throughput technologies have produced a vast amount of protein sequences, while the number of high-resolution structures has seen a limited increase. This has impelled the production of many strategies to built protein structures from its sequence, generating a considerable amount of alternative models. The selection of the closest model to the native conformation has thus become crucial for structure prediction. Several methods have been developed to score protein models by energies, knowledge-based potentials and combination of both.Results: Here, we present and demonstrate a theory to split the knowledge-based potentials in scoring terms biologically meaningful and to combine them in new scores to predict near-native structures. Our strategy allows circumventing the problem of defining the reference state. In this approach we give the proof for a simple and linear application that can be further improved by optimizing the combination of Zscores. Using the simplest composite score () we obtained predictions similar to state-of-the-art methods. Besides, our approach has the advantage of identifying the most relevant terms involved in the stability of the protein structure. Finally, we also use the composite Zscores to assess the conformation of models and to detect local errors.Conclusion: We have introduced a method to split knowledge-based potentials and to solve the problem of defining a reference state. The new scores have detected near-native structures as accurately as state-of-art methods and have been successful to identify wrongly modeled regions of many near-native conformations.

Divergent evolution and purifying selection of the flaA gene sequences in Aeromonas

BACKGROUND: The bacterial flagellum is the most important organelle of motility in bacteria and plays a key role in many bacterial lifestyles, including virulence. The flagellum also provides a paradigm of how hierarchical gene regulation, intricate protein-protein interactions and controlled protein secretion can result in the assembly of a complex multi-protein structure tightly orchestrated in time and space. As if to stress its importance, plants and animals produce receptors specifically dedicated to the recognition of flagella. Aside from motility, the flagellum also moonlights as an adhesion and has been adapted by humans as a tool for peptide display. Flagellar sequence variation constitutes a marker with widespread potential uses for studies of population genetics and phylogeny of bacterial species. RESULTS: We sequenced the complete flagellin gene (flaA) in 18 different species and subspecies of Aeromonas. Sequences ranged in size from 870 (A. allosaccharophila) to 921 nucleotides (A. popoffii). The multiple alignment displayed 924 sites, 66 of which presented alignment gaps. The phylogenetic tree revealed the existence of two groups of species exhibiting different FlaA flagellins (FlaA1 and FlaA2). Maximum likelihood models of codon substitution were used to analyze flaA sequences. Likelihood ratio tests suggested a low variation in selective pressure among lineages, with an omega ratio of less than 1 indicating the presence of purifying selection in almost all cases. Only one site under potential diversifying selection was identified (isoleucine in position 179). However, 17 amino acid positions were inferred as sites that are likely to be under positive selection using the branch-site model. Ancestral reconstruction revealed that these 17 amino acids were among the amino acid changes detected in the ancestral sequence. CONCLUSION: The models applied to our set of sequences allowed us to determine the possible evolutionary pathway followed by the flaA gene in Aeromonas, suggesting that this gene have probably been evolving independently in the two groups of Aeromonas species since the divergence of a distant common ancestor after one or several episodes of positive selection. REVIEWERS: This article was reviewed by Alexey Kondrashov, John Logsdon and Olivier Tenaillon (nominated by Laurence D Hurst).

Meta-structure correlation in protein space unveils different selection rules for folded and intrinsically disordered proteins

The number of existing protein sequences spans a very small fraction of sequence space. Natural proteins have overcome a strong negative selective pressure to avoid the formation of insoluble aggregates. Stably folded globular proteins and intrinsically disordered proteins (IDP) use alternative solutions to the aggregation problem. While in globular proteins folding minimizes the access to aggregation prone regions IDPs on average display large exposed contact areas. Here, we introduce the concept of average meta-structure correlation map to analyze sequence space. Using this novel conceptual view we show that representative ensembles of folded and ID proteins show distinct characteristics and responds differently to sequence randomization. By studying the way evolutionary constraints act on IDPs to disable a negative function (aggregation) we might gain insight into the mechanisms by which function - enabling information is encoded in IDPs.

Vacuole membrane protein 1 is an endoplasmic reticulum protein required for organelle biogenesis, protein secretion, and development

Vacuole membrane protein 1 (Vmp1) is membrane protein of unknown molecular function that has been associated with pancreatitis and cancer. The social amoeba Dictyostelium discoideum has a vmp1-related gene that we identified previously in a functional genomic study. Loss-of-function of this gene leads to a severe phenotype that compromises Dictyostelium growth and development. The expression of mammalian Vmp1 in a vmp1 Dictyostelium mutant complemented the phenotype, suggesting a functional conservation of the protein among evolutionarily distant species and highlights Dictyostelium as a valid experimental system to address the function of this gene. Dictyostelium Vmp1 is an endoplasmic reticulum protein necessary for the integrity of this organelle. Cells deficient in Vmp1 display pleiotropic defects in the secretory pathway and organelle biogenesis. The contractile vacuole, which is necessary to survive under hypoosmotic conditions, is not functional in the mutant. The structure of the Golgi apparatus, the function of the endocytic pathway and conventional protein secretion are also affected in these cells. Transmission electron microscopy of vmp1 cells showed the accumulation of autophagic features that suggests a role of Vmp1 in macroautophagy. In addition to these defects observed at the vegetative stage, the onset of multicellular development and early developmental gene expression are also compromised.

Clathrin Assembly Lymphoid Myeloid Leukemia (CALM) Protein: Localization in Endocytic-coated Pits, Interactions with Clathrin, and the Impact of Overexpression on Clathrin-mediated Traffic

The clathrin assembly lymphoid myeloid leukemia (CALM) gene encodes a putative homologue of the clathrin assembly synaptic protein AP180. Hence the biochemical properties, the subcellular localization, and the role in endocytosis of a CALM protein were studied. In vitro binding and coimmunoprecipitation demonstrated that the clathrin heavy chain is the major binding partner of CALM. The bulk of cellular CALM was associated with the membrane fractions of the cell and localized to clathrin-coated areas of the plasma membrane. In the membrane fraction, CALM was present at near stoichiometric amounts relative to clathrin. To perform structure-function analysis of CALM, we engineered chimeric fusion proteins of CALM and its fragments with the green fluorescent protein (GFP). GFP-CALM was targeted to the plasma membrane-coated pits and also found colocalized with clathrin in the Golgi area. High levels of expression of GFP-CALM or its fragments with clathrin-binding activity inhibited the endocytosis of transferrin and epidermal growth factor receptors and altered the steady-state distribution of the mannose-6-phosphate receptor in the cell. In addition, GFP-CALM overexpression caused the loss of clathrin accumulation in the trans-Golgi network area, whereas the localization of the clathrin adaptor protein complex 1 in the trans-Golgi network remained unaffected. The ability of the GFP-tagged fragments of CALM to affect clathrin-mediated processes correlated with the targeting of the fragments to clathrin-coated areas and their clathrin-binding capacities. Clathrin-CALM interaction seems to be regulated by multiple contact interfaces. The C-terminal part of CALM binds clathrin heavy chain, although the full-length protein exhibited maximal ability for interaction. Altogether, the data suggest that CALM is an important component of coated pit internalization machinery, possibly involved in the regulation of clathrin recruitment to the membrane and/or the formation of the coated pit.

A T42A Ran Mutation: Differential Interactions with Effectors and Regulators, and Defect in Nuclear Protein Import

Ran, the small, predominantly nuclear GTPase, has been implicated in the regulation of a variety of cellular processes including cell cycle progression, nuclear-cytoplasmic trafficking of RNA and protein, nuclear structure, and DNA synthesis. It is not known whether Ran functions directly in each process or whether many of its roles may be secondary to a direct role in only one, for example, nuclear protein import. To identify biochemical links between Ran and its functional target(s), we have generated and examined the properties of a putative Ran effector mutation, T42A-Ran. T42A-Ran binds guanine nucleotides as well as wild-type Ran and responds as well as wild-type Ran to GTP or GDP exchange stimulated by the Ran-specific guanine nucleotide exchange factor, RCC1. T42A-Ran·GDP also retains the ability to bind p10/NTF2, a component of the nuclear import pathway. In contrast to wild-type Ran, T42A-Ran·GTP binds very weakly or not detectably to three proposed Ran effectors, Ran-binding protein 1 (RanBP1), Ran-binding protein 2 (RanBP2, a nucleoporin), and karyopherin ß (a component of the nuclear protein import pathway), and is not stimulated to hydrolyze bound GTP by Ran GTPase-activating protein, RanGAP1. Also in contrast to wild-type Ran, T42A-Ran does not stimulate nuclear protein import in a digitonin permeabilized cell assay and also inhibits wild-type Ran function in this system. However, the T42A mutation does not block the docking of karyophilic substrates at the nuclear pore. These properties of T42A-Ran are consistent with its classification as an effector mutant and define the exposed region of Ran containing the mutation as a probable effector loop.

T-Coffee: a web server for the multiple sequence alignment of protein and RNA sequences using structural information and homology extension

This article introduces a new interface for T-Coffee, a consistency-based multiple sequence alignment program. This interface provides an easy and intuitive access to the most popular functionality of the package. These include the default T-Coffee mode for protein and nucleic acid sequences, the M-Coffee mode that allows combining the output of any other aligners, and template-based modes of T-Coffee that deliver high accuracy alignments while using structural or homology derived templates. These three available template modes are Expresso for the alignment of protein with a known 3D-Structure, R-Coffee to align RNA sequences with conserved secondary structures and PSI-Coffee to accurately align distantly related sequences using homology extension. The new server benefits from recent improvements of the T-Coffee algorithm and can align up to 150 sequences as long as 10 000 residues and is available from both http://www.tcoffee.org and its main mirror http://tcoffee.crg.cat.

Protein oligomers studied by solid-state NMR the case of the full-length nucleoid-associated protein histone-like nucleoid structuring protein

Members of the histone-like nucleoid structuring protein (H-NS) family play roles both as architectural proteins and as modulators of gene expression in Gram-negative bacteria. The H-NS protein participates in modulatory processes that respond to environmental changes in osmolarity, pH, or temperature. H-NS oligomerization is essential for its activity. Structural models of different truncated forms are available. However, high-resolution structural details of full-length H-NS and its DNA-bound state have largely remained elusive. We report on progress in characterizing the biologically active H-NS oligomers with solid-state NMR. We compared uniformly ((13)C,(15)N)-labeled ssNMR preparations of the isolated N-terminal region (H-NS 1-47) and full-length H-NS (H-NS 1-137). In both cases, we obtained ssNMR spectra of good quality and characteristic of well-folded proteins. Analysis of the results of 2D and 3D (13)C-(13)C and (15)N-(13)C correlation experiments conducted at high magnetic field led to assignments of residues located in different topological regions of the free full-length H-NS. These findings confirm that the structure of the N-terminal dimerization domain is conserved in the oligomeric full-length protein. Small changes in the dimerization interface suggested by localized chemical shift variations between solution and solid-state spectra may be relevant for DNA recoginition.

Uncovering the functional constraints underlying the genomic organisation of the Odorant-Binding Protein genes

Animal olfactory systems have a critical role for the survival and reproduction of individuals. In insects, the odorant-binding proteins (OBPs) are encoded by a moderately sized gene family, and mediate the first steps of the olfactory processing. Most OBPs are organized in clusters of a few paralogs, which are conserved over time. Currently, the biological mechanism explaining the close physical proximity among OBPs is not yet established. Here, we conducted a comprehensive study aiming to gain insights into the mechanisms underlying the OBP genomic organization. We found that the OBP clusters are embedded within large conserved arrangements. These organizations also include other non-OBP genes, which often encode proteins integral to plasma membrane. Moreover, the conservation degree of such large clusters is related to the following: 1) the promoter architecture of the confined genes, 2) a characteristic transcriptional environment, and 3) the chromatin conformation of the chromosomal region. Our results suggest that chromatin domains may restrict the location of OBP genes to regions having the appropriate transcriptional environment, leading to the OBP cluster structure. However, the appropriate transcriptional environment for OBP and the other neighbor genes is not dominated by reduced levels of expression noise. Indeed, the stochastic fluctuations in the OBP transcript abundance may have a critical role in the combinatorial nature of the olfactory coding process.

Genetic variation in the European hake, Merluccius merluccius. Description of protein loci and genetic divergence between Atlantic and Mediterranean populations = Anàlisi de la variabilitat genètica en el lluç europeu, Merluccius merluccius. Descripció dels loci electroforètics i divergència genètica entre les poblacions atlàntica i mediterrània

We examined the genetic population structure of the european hake (Merluccius merluccius) using electrophoretically detectable population markers in 35 protein loci. Samples were collected from 7 locations in the Atlantic Ocean and Mediterranean Sea. Six loci were polymorphic using the 0.05 criterion of polymorphism. Sample heterozigosities ranged from 0.052 to 0.072 and averaged 0.0625. In this study, significant allele frequency differences were detected between Atlantic and Mediterranean populations in three polymorphic loci: GAPDH-1*, GPI-2* and SOD-1*. Two major genetic groups were considered: a North-Atlantic stock and the Mediterranean stock. The Nei genetic distance, D, (based on 33 loci) between samples from these two groups ranged from 0.002 to 0.006. Genetic differenciation between these areas appears to reflect the barrier effect of Strait of Gibraltar. On average over loci, 96.92 % of the total gene diversity was contained within samples, 0.23 % expressed differences among locations within areas, and 2.64 % differences between regions. A review of morphological variation together with the genetic data presented here suggest that the populations of hake from these areas are subdivided into two different stocks: the North-Atlantic stock and the Mediterranean stock. The most conservative approach to the management of these stocks is to consider the Atlantic and Mediterranean stocks independently from oneanother

Structure and dynamics of the membrane attaching nitric oxide transporter nitrophorin 7

Nitrophorins represent a unique class of heme proteins that are able to perform the delicate transportation and release of the free-radical gaseous messenger nitric oxide (NO) in a pH-triggered manner. Besides its ability to bind to phospholipid membranes, the N-terminus contains an additional Leu-Pro-Gly stretch, which is a unique sequence trait, and the heme cavity is significantly altered with respect to other nitrophorins. These distinctive features encouraged us to solve the X-ray crystallographic structures of NP7 at low and high pH and bound with different heme ligands (nitric oxide, histamine, imidazole). The overall fold of the lipocalin motif is well preserved in the different X-ray structures and resembles the fold of other nitrophorins. However, a chain-like arrangement in the crystal lattice due to a number of head-to-tail electrostatic stabilizing interactions is found in NP7. Furthermore, the X-ray structures also reveal ligand-dependent changes in the orientation of the heme, as well as in specific interactions between the A-B and G-H loops, which are considered to be relevant for the biological function of nitrophorins. Fast and ultrafast laser triggered ligand rebinding experiments demonstrate the pH-dependent ligand migration within the cavities and the exit route. Finally, the topological distribution of pockets located around the heme as well as from inner cavities present at the rear of the protein provides a distinctive feature in NP7, so that while a loop gated exit mechanism to the solvent has been proposed for most nitrophorins, a more complex mechanism that involves several interconnected gas hosting cavities is proposed for NP7.