Dissecting the role of low-complexity regions in the evolution of vertebrate proteins

Low-complexity regions (LCRs) in proteins are tracts that are highly enriched in one or a few aminoacids. Given their high abundance, and their capacity to expand in relatively short periods of time through replication slippage, they can greatly contribute to increase protein sequence space and generate novel protein functions. However, little is known about the global impact of LCRs on protein evolution. We have traced back the evolutionary history of 2,802 LCRs from a large set of homologous protein families from H.sapiens, M.musculus, G.gallus, D.rerio and C.intestinalis. Transcriptional factors and other regulatory functions are overrepresented in proteins containing LCRs. We have found that the gain of novel LCRs is frequently associated with repeat expansion whereas the loss of LCRs is more often due to accumulation of amino acid substitutions as opposed to deletions. This dichotomy results in net protein sequence gain over time. We have detected a significant increase in the rate of accumulation of novel LCRs in the ancestral Amniota and mammalian branches, and a reduction in the chicken branch. Alanine and/or glycine-rich LCRs are overrepresented in recently emerged LCR sets from all branches, suggesting that their expansion is better tolerated than for other LCR types. LCRs enriched in positively charged amino acids show the contrary pattern, indicating an important effect of purifying selection in their maintenance. We have performed the first large-scale study on the evolutionary dynamics of LCRs in protein families. The study has shown that the composition of an LCR is an important determinant of its evolutionary pattern.

Comparative genomic analysis of the odorant-binding protein family in 12 Drosophila genomes: purifying selection and birth-and-death evolution

Background: Chemoreception is a widespread mechanism that is involved in critical biologic processes, including individual and social behavior. The insect peripheral olfactory system comprises three major multigene families: the olfactory receptor (Or), the gustatory receptor (Gr), and the odorant-binding protein (OBP) families. Members of the latter family establish the first contact with the odorants, and thus constitute the first step in the chemosensory transduction pathway.Results: Comparative analysis of the OBP family in 12 Drosophila genomes allowed the identification of 595 genes that encode putative functional and nonfunctional members in extant species, with 43 gene gains and 28 gene losses (15 deletions and 13 pseudogenization events). The evolution of this family shows tandem gene duplication events, progressive divergence in DNA and amino acid sequence, and prevalence of pseudogenization events in external branches of the phylogenetic tree. We observed that the OBP arrangement in clusters is maintained across the Drosophila species and that purifying selection governs the evolution of the family; nevertheless, OBP genes differ in their functional constraints levels. Finally, we detect that the OBP repertoire evolves more rapidly in the specialist lineages of the Drosophila melanogaster group (D. sechellia and D. erecta) than in their closest generalists.Conclusion: Overall, the evolution of the OBP multigene family is consistent with the birth-and-death model. We also found that members of this family exhibit different functional constraints, which is indicative of some functional divergence, and that they might be involved in some of the specialization processes that occurred through the diversification of the Drosophila genus.

Identification of potential small molecule binding pockets on Rho family GTPases

Rho GTPases are conformational switches that control a wide variety of signaling pathways critical for eukaryotic cell development and proliferation. They represent attractive targets for drug design as their aberrant function and deregulated activity is associated with many human diseases including cancer. Extensive high-resolution structures (.100) and recent mutagenesis studies have laid the foundation for the design of new structure-based chemotherapeutic strategies. Although the inhibition of Rho signaling with drug-like compounds is an active area of current research, very little attention has been devoted to directly inhibiting Rho by targeting potential allosteric non-nucleotide binding sites. By avoiding the nucleotide binding site, compounds may minimize the potential for undesirable off-target interactions with other ubiquitous GTP and ATP binding proteins. Here we describe the application of molecular dynamics simulations, principal component analysis, sequence conservation analysis, and ensemble small-molecule fragment mapping to provide an extensive mapping of potential small-molecule binding pockets on Rho family members. Characterized sites include novel pockets in the vicinity of the conformationaly responsive switch regions as well as distal sites that appear to be related to the conformations of the nucleotide binding region. Furthermore the use of accelerated molecular dynamics simulation, an advanced sampling method that extends the accessible time-scale of conventional simulations, is found to enhance the characterization of novel binding sites when conformational changes are important for the protein mechanism.

Podocalyxin is a novel polysialylated neural adhesion protein with multiple roles in neural development and synapse formation

Neural development and plasticity are regulated by neural adhesion proteins, including the polysialylated form of NCAM (PSA-NCAM). Podocalyxin (PC) is a renal PSA-containing protein that has been reported to function as an anti-adhesin in kidney podocytes. Here we show that PC is widely expressed in neurons during neural development. Neural PC interacts with the ERM protein family, and with NHERF1/2 and RhoA/G. Experiments in vitro and phenotypic analyses of podxl-deficient mice indicate that PC is involved in neurite growth, branching and axonal fasciculation, and that PC loss-of-function reduces the number of synapses in the CNS and in the neuromuscular system. We also show that whereas some of the brain PC functions require PSA, others depend on PC per se. Our results show that PC, the second highly sialylated neural adhesion protein, plays multiple roles in neural development.

RPS4Y gene family evolution in primates

Background: The RPS4 gene codifies for ribosomal protein S4, a very well-conserved protein present in all kingdoms. In primates, RPS4 is codified by two functional genes located on both sex chromosomes: the RPS4X and RPS4Y genes. In humans, RPS4Y is duplicated and the Y chromosome therefore carries a third functional paralog: RPS4Y2, which presents a testis-specific expression pattern. Results: DNA sequence analysis of the intronic and cDNA regions of RPS4Y genes from species covering the entire primate phylogeny showed that the duplication event leading to the second Y-linked copy occurred after the divergence of New World monkeys, about 35 million years ago. Maximum likelihood analyses of the synonymous and non-synonymous substitutions revealed that positive selection was acting on RPS4Y2 gene in the human lineage, which represents the first evidence of positive selection on a ribosomal protein gene. Putative positive amino acid replacements affected the three domains of the protein: one of these changes is located in the KOW protein domain and affects the unique invariable position of this motif, and might thus have a dramatic effect on the protein function.Conclusion: Here, we shed new light on the evolutionary history of RPS4Y gene family, especially on that of RPS4Y2. The results point that the RPS4Y1 gene might be maintained to compensate gene dosage between sexes, while RPS4Y2 might have acquired a new function, at least in the lineage leading to humans.

ORMDL proteins are a conserved new family of endoplasmic reticulum membrane proteins

Background: Annotations of completely sequenced genomes reveal that nearly half of the genes identified are of unknown function, and that some belong to uncharacterized gene families. To help resolve such issues, information can be obtained from the comparative analysis of homologous genes in model organisms. Results: While characterizing genes from the retinitis pigmentosa locus RP26 at 2q31-q33, we have identified a new gene, ORMDL1, that belongs to a novel gene family comprising three genes in humans (ORMDL1, ORMDL2 and ORMDL3), and homologs in yeast, microsporidia, plants, Drosophila, urochordates and vertebrates. The human genes are expressed ubiquitously in adult and fetal tissues. The Drosophila ORMDL homolog is also expressed throughout embryonic and larval stages, particularly in ectodermally derived tissues. The ORMDL genes encode transmembrane proteins anchored in the endoplasmic reticulum (ER). Double knockout of the two Saccharomyces cerevisiae homologs leads to decreased growth rate and greater sensitivity to tunicamycin and dithiothreitol. Yeast mutants can be rescued by human ORMDL homologs. Conclusions: From protein sequence comparisons we have defined a novel gene family, not previously recognized because of the absence of a characterized functional signature. The sequence conservation of this family from yeast to vertebrates, the maintenance of duplicate copies in different lineages, the ubiquitous pattern of expression in human and Drosophila, the partial functional redundancy of the yeast homologs and phenotypic rescue by the human homologs, strongly support functional conservation. Subcellular localization and the response of yeast mutants to specific agents point to the involvement of ORMDL in protein folding in the ER.

Podocalyxin is a novel polysialylated neural adhesion protein with multiple roles in neural development and synapse formation

Neural development and plasticity are regulated by neural adhesion proteins, including the polysialylated form of NCAM (PSA-NCAM). Podocalyxin (PC) is a renal PSA-containing protein that has been reported to function as an anti-adhesin in kidney podocytes. Here we show that PC is widely expressed in neurons during neural development. Neural PC interacts with the ERM protein family, and with NHERF1/2 and RhoA/G. Experiments in vitro and phenotypic analyses of podxl-deficient mice indicate that PC is involved in neurite growth, branching and axonal fasciculation, and that PC loss-of-function reduces the number of synapses in the CNS and in the neuromuscular system. We also show that whereas some of the brain PC functions require PSA, others depend on PC per se. Our results show that PC, the second highly sialylated neural adhesion protein, plays multiple roles in neural development.

Synthesis and multi-target biological profiling of a novel family of rhein derivatives as disease-modifying anti-Alzheimer agents

We have synthesized a family of rheinhuprine hybrids to hit several key targets for Alzheimer"s disease. Biological screening performed in vitro and in Escherichia coli cells has shown that these hybrids exhibit potent inhibitory activities against human acetylcholinesterase butyrylcholinesterase, and BACE-1, dual Aβ42 and tau anti-aggregating activity, and brain permeability. Ex vivo studies with the leads (+)- and ()-7e in brain slices of C57bl6 mice have revealed that they efficiently protect against the Aβ-induced synaptic dysfunction , preventing the loss of synaptic proteins and/or have a positive effect on the induction of long term potentiation. In vivo studies in APP-PS1 transgenic mice treated i.p. for 4 weeks with (+)- and ()-7e have shown a central soluble Aβ lowering effect, accompanied by an increase in the levels of mature amyloid precursor protein (APP). Thus, (+)- and ()-7e emerge as very promising disease-modifying anti-Alzheimer drug candidates.

Synthesis and multi-target biological profiling of a novel family of rhein derivatives as disease-modifying anti-Alzheimer agents

We have synthesized a family of rheinhuprine hybrids to hit several key targets for Alzheimer"s disease. Biological screening performed in vitro and in Escherichia coli cells has shown that these hybrids exhibit potent inhibitory activities against human acetylcholinesterase butyrylcholinesterase, and BACE-1, dual Aβ42 and tau anti-aggregating activity, and brain permeability. Ex vivo studies with the leads (+)- and ()-7e in brain slices of C57bl6 mice have revealed that they efficiently protect against the Aβ-induced synaptic dysfunction , preventing the loss of synaptic proteins and/or have a positive effect on the induction of long term potentiation. In vivo studies in APP-PS1 transgenic mice treated i.p. for 4 weeks with (+)- and ()-7e have shown a central soluble Aβ lowering effect, accompanied by an increase in the levels of mature amyloid precursor protein (APP). Thus, (+)- and ()-7e emerge as very promising disease-modifying anti-Alzheimer drug candidates.

ORMDL proteins are a conserved new family of endoplasmic reticulum membrane proteins

Background: Annotations of completely sequenced genomes reveal that nearly half of the genes identified are of unknown function, and that some belong to uncharacterized gene families. To help resolve such issues, information can be obtained from the comparative analysis of homologous genes in model organisms. Results: While characterizing genes from the retinitis pigmentosa locus RP26 at 2q31-q33, we have identified a new gene, ORMDL1, that belongs to a novel gene family comprising three genes in humans (ORMDL1, ORMDL2 and ORMDL3), and homologs in yeast, microsporidia, plants, Drosophila, urochordates and vertebrates. The human genes are expressed ubiquitously in adult and fetal tissues. The Drosophila ORMDL homolog is also expressed throughout embryonic and larval stages, particularly in ectodermally derived tissues. The ORMDL genes encode transmembrane proteins anchored in the endoplasmic reticulum (ER). Double knockout of the two Saccharomyces cerevisiae homologs leads to decreased growth rate and greater sensitivity to tunicamycin and dithiothreitol. Yeast mutants can be rescued by human ORMDL homologs. Conclusions: From protein sequence comparisons we have defined a novel gene family, not previously recognized because of the absence of a characterized functional signature. The sequence conservation of this family from yeast to vertebrates, the maintenance of duplicate copies in different lineages, the ubiquitous pattern of expression in human and Drosophila, the partial functional redundancy of the yeast homologs and phenotypic rescue by the human homologs, strongly support functional conservation. Subcellular localization and the response of yeast mutants to specific agents point to the involvement of ORMDL in protein folding in the ER.

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.

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.

Hepatocyte nuclear factor-4 alpha regulates the human apolipoprotein AV gene: Identification of a novel response element and involvement in the control by peroxisome proliferator-activated receptor-gamma coactivator-1 alpha, AMP-activated protein kinase, and mitogen-activated protein kinase pathway

The recently discovered apolipoprotein AV (apoAV) gene has been reported to be a key player in modulating plasma triglyceride levels. Here we identify the hepatocyte nuclear factor-4 (HNF-4 ) as a novel regulator of human apoAV gene. Inhibition of HNF-4 expression by small interfering RNA resulted in down-regulation of apoAV. Deletion, mutagenesis, and binding assays revealed that HNF-4 directly regulates human apoAV promoter through DR1 [a direct repeat separated by one nucleotide (nt)], and via a novel element for HNF-4 consisting of an inverted repeat separated by 8 nt (IR8). In addition, we show that the coactivator peroxisome proliferator-activated receptor- coactivator-1 was capable of stimulating the HNF-4 -dependent transactivation of apoAV promoter. Furthermore, analyses in human hepatic cells demonstrated that AMP-activated protein kinase (AMPK) and the MAPK signaling pathway regulate human apoAV expression and suggested that this regulation may be mediated, at least in part, by changes in HNF-4 . Intriguingly, EMSAs and mice with a liver-specific disruption of the HNF-4 gene revealed a species-distinct regulation of apoAV by HNF-4 , which resembles that of a subset of HNF-4 target genes. Taken together, our data provide new insights into the binding properties and the modulation of HNF-4 and underscore the role of HNF-4 in regulating triglyceride metabolism.

Lipid phosphate phosphatase 3 participates in transport carrier formation and protein trafficking in the early secretory pathway

The inhibition of phosphatidic acid phosphatase (PAP) activity by propanolol indicates that diacylglycerol (DAG) is required for the formation of transport carriers at the Golgi and for retrograde trafficking to the ER. Here we report that the PAP2 family member lipid phosphate phosphatase 3 (LPP3, also known as PAP2b) localizes in compartments of the secretory pathway from ER export sites to the Golgi complex. The depletion of human LPP3: (i) reduces the number of tubules generated from the ER-Golgi intermediate compartment and the Golgi, with those formed from the Golgi being longer in LPP3-silenced cells than in control cells; (ii) impairs the Rab6-dependent retrograde transport of Shiga toxin subunit B from the Golgi to the ER, but not the anterograde transport of VSV-G or ssDsRed; and (iii) induces a high accumulation of Golgi-associated membrane buds. LPP3 depletion also reduces levels of de novo synthesized DAG and the Golgi-associated DAG contents. Remarkably, overexpression of a catalytically inactive form of LPP3 mimics the effects of LPP3 knockdown on Rab6-dependent retrograde transport. We conclude that LPP3 participates in the formation of retrograde transport carriers at the ER-Golgi interface, where it transitorily cycles, and during its route to the plasma membrane.

Family Size Evolution in Drosophila Chemosensory Gene Families: A Comparative Analysis with a Critical Appraisal of Methods

Gene turnover rates and the evolution of gene family sizes are important aspects of genome evolution. Here, we use curated sequence data of the major chemosensory gene families from Drosophila-the gustatory receptor, odorant receptor, ionotropic receptor, and odorant-binding protein families-to conduct a comparative analysis among families, exploring different methods to estimate gene birth and death rates, including an ad hoc simulation study. Remarkably, we found that the state-of-the-art methods may produce very different rate estimates, which may lead to disparate conclusions regarding the evolution of chemosensory gene family sizes in Drosophila. Among biological factors, we found that a peculiarity of D. sechellia's gene turnover rates was a major source of bias in global estimates, whereas gene conversion had negligible effects for the families analyzed herein. Turnover rates vary considerably among families, subfamilies, and ortholog groups although all analyzed families were quite dynamic in terms of gene turnover. Computer simulations showed that the methods that use ortholog group information appear to be the most accurate for the Drosophila chemosensory families. Most importantly, these results reveal the potential of rate heterogeneity among lineages to severely bias some turnover rate estimation methods and the need of further evaluating the performance of these methods in a more diverse sampling of gene families and phylogenetic contexts. Using branch-specific codon substitution models, we find further evidence of positive selection in recently duplicated genes, which attests to a nonneutral aspect of the gene birth-and-death process.