SelenoDB 1.0 : a database of selenoprotein genes, proteins and SECIS elements

Selenoproteins are a diverse group of proteinsusually misidentified and misannotated in sequencedatabases. The presence of an in-frame UGA (stop)codon in the coding sequence of selenoproteingenes precludes their identification and correctannotation. The in-frame UGA codons are recodedto cotranslationally incorporate selenocysteine,a rare selenium-containing amino acid. The developmentof ad hoc experimental and, more recently,computational approaches have allowed the efficientidentification and characterization of theselenoproteomes of a growing number of species.Today, dozens of selenoprotein families have beendescribed and more are being discovered in recentlysequenced species, but the correct genomic annotationis not available for the majority of thesegenes. SelenoDB is a long-term project that aims toprovide, through the collaborative effort of experimentaland computational researchers, automaticand manually curated annotations of selenoproteingenes, proteins and SECIS elements. Version 1.0 ofthe database includes an initial set of eukaryoticgenomic annotations, with special emphasis on thehuman selenoproteome, for immediate inspectionby selenium researchers or incorporation into moregeneral databases. SelenoDB is freely available athttp://www.selenodb.org.

Characterization of alternatively spliced products and tissue-specific isoforms of USP28 and USP25

Background: The ubiquitin-dependent protein degradation pathway is essential for the proteolysis of intracellular proteins and peptides. Deubiquitinating enzymes constitute a complex protein family involved in a multitude of cellular processes. The ubiquitin-specific proteases (UBP) are a group of enzymes whose predicted function is to reverse the ubiquitinating reaction by removing ubiquitin from a large variety of substrates. We have lately reported the characterization of human USP25, a specific-ubiquitin protease gene at 21q11.2, with a specific pattern of expression in murine fetal brains and adult testis. Results: Database homology searches at the DNA and protein levels and cDNA library screenings led to the identification of a new UBP member in the human genome, named USP28, at 11q23. This novel gene showed preferential expression in heart and muscle. Moreover, cDNA, expressed sequence tag and RT-PCR analyses provided evidence for alternatively spliced products and tissue-specific isoforms. Concerning function, USP25 overexpression in Down syndrome fetal brains was shown by real-time PCR. Conclusions: On the basis of the genomic and protein sequence as well as the functional data, USP28 and USP25 establish a new subfamily of deubiquitinating enzymes. Both genes have alternatively spliced exons that could generate protein isoforms with distinct tissue-specific activity. The overexpression of USP25 in Down syndrome fetal brains supports the gene-dosage effects suggested for other UBP members related to aneuploidy syndromes.

Insertion of a malE B-Galactosidase fusion protein into the envelope of Escherichia coli disrupts biogenesis of outer membrane proteins and processing of inner membrane proteins

The synthesis of a membrane-bound MalE ,B-galactosidase hybrid protein, when induced by growth of Escherichia coli on maltose, leads to inhibition of cell division and eventually a reduced rate of mass increase. In addition, the relative rate of synthesis of outer membrane proteins, but not that of inner membrane proteins, was reduced by about 50%o. Kinetic experiments demonstrated that this reduction coincided with the period of maximum synthesis of the hybrid protein (and another maltose-inducible protein, LamB). The accumulation of this abnormal protein in the envelope therefore appeared specifically to inhibit the synthesis, the assembly of outer membrane proteins, or both, indicating that the hybrid protein blocks some export site or causes the sequestration of some limiting factor(s) involved in the export process. Since the MalE protein is normally located in the periplasm, the results also suggest that the synthesis of periplasmic and outer membrane proteins may involve some steps in common. The reduced rate of synthesis of outer membrane proteins was also accompanied by the accumulation in the envelope of at least one outer membrane protein and at least two inner membrane proteins as higher-molecular-weight forms, indicating that processing (removal of the N-terminal signal sequence) was also disrupted by the presence of the hybrid protein. These results may indicate that the assembly of these membrane proteins is blocked at a relatively late step rather than at the level of primary recognition of some site by the signal sequence. In addition, the results suggest that some step common to the biogenesis of quite different kinds of envelope protein is blocked by the presence of the hybrid protein.

Bcl-2 family proteins and cytoskeleton changes involved in DM-1 cytotoxic effect on melanoma cells

Melanoma is one of the most aggressive types of skin cancer and its incidence rate is still increasing. All existing treatments are minimally effective. Consequently, new therapeutic agents for melanoma treatment should be developed. The DM-1 compound is a curcumin analog that possesses several curcumin characteristics, such as antiproliferative, antitumor, and anti-metastatic properties. The aim of this study was to evaluate the different signaling pathways involved in the cytotoxic effect of DM-1 on melanoma cells. The apoptotic process and cytoskeletal changes were evaluated by immunoblotting and immunofluorescence, respectively, in melanoma cells. After DM-1 treatment, SK-MEL-5 melanoma cells showed actin filament disorganization with spicule formation throughout the cytoskeleton and significant reduction of focal adhesion as well as they were present only at cell extremities, conferring a poor connection between the cell and the substrate. Besides this, there was significant filopodium retraction and loss of typical cytoskeleton scaffold. These modifications contributed to cell detachment followed by cell death. Furthermore, DM-1-induced apoptosis was triggered by multiple Bcl-2 proteins involved in both the extrinsic and the intrinsic apoptotic pathways. SK-MEL-5 cells showed a death mechanism mainly by Bcl-2/Bax ratio decrease, whereas A375 cells presented apoptosis induction by Mcl-1 and Bcl-xL downregulation. In SK-MEL-5 and A375 melanoma cells, there was a significant increase in the active form of caspase 9, and the inactive form of the effector caspase 3 was decreased in both cell lines. Expression of cleaved poly ADP ribose polymerase was increased after DM-1 treatment in these melanoma cell lines, demonstrating that the apoptotic process occurred. Altogether, these data elucidate the cellular and molecular mechanisms involved in the cytotoxicity induced by the antitumor agent DM-1 in melanoma cells.

rDock: A Fast, Versatile and Open Source Program for Docking Ligands to Proteins and Nucleic Acids

Identification of chemical compounds with specific biological activities is an important step in both chemical biology and drug discovery. When the structure of the intended target is available, one approach is to use molecular docking programs to assess the chemical complementarity of small molecules with the target; such calculations provide a qualitative measure of affinity that can be used in virtual screening (VS) to rank order a list of compounds according to their potential to be active. rDock is a molecular docking program developed at Vernalis for high-throughput VS (HTVS) applications. Evolved from RiboDock, the program can be used against proteins and nucleic acids, is designed to be computationally very efficient and allows the user to incorporate additional constraints and information as a bias to guide docking. This article provides an overview of the program structure and features and compares rDock to two reference programs, AutoDock Vina (open source) and Schrodinger's Glide (commercial). In terms of computational speed for VS, rDock is faster than Vina and comparable to Glide. For binding mode prediction, rDock and Vina are superior to Glide. The VS performance of rDock is significantly better than Vina, but inferior to Glide for most systems unless pharmacophore constraints are used; in that case rDock and Glide are of equal performance. The program is released under the Lesser General Public License and is freely available for download, together with the manuals, example files and the complete test sets, at http://rdock.sourceforge.net/

Combinatorial Synthesis and Biological Evaluation of Inhibitors of D-Ala-D-Ala-Ligase

Report for the scientific sojourn carried out at the Max Planck Institut of Molecular Phisiology, Germany, from 2006 to 2008.The work carried out during this postdoctoral stage was focused on two different projects. Firstly, identification of D-Ala D-Ala Inhibitors and the development of new synthethic approaches to obtain lipidated peptides and proteins and the use of these lipidated proteins in biological and biophysical studies. In the first project, new D-Ala D-Ala inhibitors were identified by using structural alignments of the ATP binding sites of the bacterial ligase DDl and protein and lipid kinases in complex with ATP analogs. We tested a series of commercially available kinase inhibitors and found LFM-A13 and Tyrphostine derivatives to inhibit DDl enzyme activity. Based on the initial screening results we synthesized a series of malononitrilamide and salicylamide derivatives and were able to confirm the validity of these scaffolds as inhibitors of DDl. From this investigation we gained a better understanding of the structural requirements and limitations necessary for the preparation of ATP competitive DDl inhibitors. The compounds in this study may serve as starting points for the development of bi-substrate inhibitors that incorporate both, an ATP competitive and a substrate competitive moiety. Bisubstrate inhibitors that block the ATP and D-Ala binding sites should exhibit enhanced selectivity and potency profiles by preferentially inhibiting DDl over kinases. In the second project, an optimized synthesis for tha alkylation of cysteins using the thiol ene reaction was establisehd. This new protocol allowed us to obtain large amounts of hexadecylated cysteine that was required for the synthesis of differently lipidated peptides. Afterwards the synthesis of various N-ras peptides bearing different lipid anchors was performed and the peptides were ligated to a truncated N-ras protein. The influence of this differently lipidated N-ras proteins on the partioning and association of N-Ras in model membrane subdomains was studied using Atomic Force Microscopy.

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.

Distribution of events of positive selection and population differentiation in a metabolic pathway: the case of asparagine N-glycosylation

Background: Asparagine N-Glycosylation is one of the most important forms of protein post-translational modification in eukaryotes. This metabolic pathway can be subdivided into two parts: an upstream sub-pathway required for achieving proper folding for most of the proteins synthesized in the secretory pathway, and a downstream sub-pathway required to give variability to trans-membrane proteins, and involved in adaptation to the environment andinnate immunity. Here we analyze the nucleotide variability of the genes of this pathway in human populations, identifying which genes show greater population differentiation and which genes show signatures of recent positive selection. We also compare how these signals are distributed between the upstream and the downstream parts of the pathway, with the aim of exploring how forces of population differentiation and positive selection vary among genes involved in the same metabolic pathway but subject to different functional constraints. Results:Our results show that genes in the downstream part of the pathway are more likely to show a signature of population differentiation, while events of positive selection are equally distributed among the two parts of the pathway. Moreover, events of positive selection arefrequent on genes that are known to be at bifurcation points, and that are identified as beingin key position by a network-level analysis such as MGAT3 and GCS1.Conclusions: These findings indicate that the upstream part of the Asparagine N-Glycosylation pathway has lower diversity among populations, while the downstream part is freer to tolerate diversity among populations. Moreover, the distribution of signatures of population differentiation and positive selection can change between parts of a pathway, especially between parts that are exposed to different functional constraints. Our results support the hypothesis that genes involved in constitutive processes can be expected to show lower population differentiation, while genes involved in traits related to the environment should show higher variability. Taken together, this work broadens our knowledge on how events of population differentiation and of positive selection are distributed among different parts of a metabolic pathway.

K(+) channel expression distinguishes subpopulations of parvalbumin- and somatostatin-containing neocortical interneurons

Kv3.1 and Kv3.2 K+ channel proteins form similar voltage-gated K+ channels with unusual properties, including fast activation at voltages positive to −10 mV and very fast deactivation rates. These properties are thought to facilitate sustained high-frequency firing. Kv3.1 subunits are specifically found in fast-spiking, parvalbumin (PV)-containing cortical interneurons, and recent studies have provided support for a crucial role in the generation of the fast-spiking phenotype. Kv3.2 mRNAs are also found in a small subset of neocortical neurons, although the distribution of these neurons is different. We raised antibodies directed against Kv3.2 proteins and used dual-labeling methods to identify the neocortical neurons expressing Kv3.2 proteins and to determine their subcellular localization. Kv3.2 proteins are prominently expressed in patches in somatic and proximal dendritic membrane as well as in axons and presynaptic terminals of GABAergic interneurons. Kv3.2 subunits are found in all PV-containing neurons in deep cortical layers where they probably form heteromultimeric channels with Kv3.1 subunits. In contrast, in superficial layer PV-positive neurons Kv3.2 immunoreactivity is low, but Kv3.1 is still prominently expressed. Because Kv3.1 and Kv3.2 channels are differentially modulated by protein kinases, these results raise the possibility that the fast-spiking properties of superficial- and deep-layer PV neurons are differentially regulated by neuromodulators. Interestingly, Kv3.2 but not Kv3.1 proteins are also prominent in a subset of seemingly non-fast-spiking, somatostatin- and calbindin-containing interneurons, suggesting that the Kv3.1–Kv3.2 current type can have functions other than facilitating high-frequency firing.

Biana: a software framework for compiling biological interactions and analyzing networks

Background: The analysis and usage of biological data is hindered by the spread of information across multiple repositories and the difficulties posed by different nomenclature systems and storage formats. In particular, there is an important need for data unification in the study and use of protein-protein interactions. Without good integration strategies, it is difficult to analyze the whole set of available data and its properties.Results: We introduce BIANA (Biologic Interactions and Network Analysis), a tool for biological information integration and network management. BIANA is a Python framework designed to achieve two major goals: i) the integration of multiple sources of biological information, including biological entities and their relationships, and ii) the management of biological information as a network where entities are nodes and relationships are edges. Moreover, BIANA uses properties of proteins and genes to infer latent biomolecular relationships by transferring edges to entities sharing similar properties. BIANA is also provided as a plugin for Cytoscape, which allows users to visualize and interactively manage the data. A web interface to BIANA providing basic functionalities is also available. The software can be downloaded under GNU GPL license from http://sbi.imim.es/web/BIANA.php.Conclusions: BIANA's approach to data unification solves many of the nomenclature issues common to systems dealing with biological data. BIANA can easily be extended to handle new specific data repositories and new specific data types. The unification protocol allows BIANA to be a flexible tool suitable for different user requirements: non-expert users can use a suggested unification protocol while expert users can define their own specific unification rules.

Assembly and disassembly of the Golgi Complex: two prodesses arranged in a cis to trans direction

We have studied the disassembly and assembly of two morphologically and functionally distinct parts of the Golgi complex, the cis/middle and trans cisterna/trans network compartments. For this purpose we have followed the redistribution of three cis/middle- (GMPc-1, GMPc-2, MG 160) and two trans- (GMPt-1 and GMPt-2) Golgi membrane proteins during and after treatment of normal rat kidney (NRK) cells with brefeldin A (BFA). BFA induced complete disassembly of the cis/middle- and trans-Golgi complex and translocation of GMPc and GMPt to the ER. Cells treated for short times (3 min) with BFA showed extensive disorganization of both cis/middle- and trans-Golgi complexes. However, complete disorganization of the trans part required much longer incubations with the drug. Upon removal of BFA the Golgi complex was reassembled by a process consisting of three steps: (a) exist of cis/middle proteins from the ER and their accumulation into vesicular structures scattered throughout the cytoplasm; (b) gradual relocation and accumulation of the trans proteins in the vesicles containing the cis/middle proteins; and (c) assembly of the cisternae, and reconstruction of the Golgi complex within an area located in the vicinity of the centrosome from which the ER was excluded. Reconstruction of the cis/middle-Golgi complex occurred under temperature conditions inhibitory of the reorganization of the trans-Golgi complex, and was dependent on microtubules. Reconstruction of the trans-Golgi complex, disrupted with nocodazole after selective fusion of the cis/middle-Golgi complex with the ER, occurred after the release of cis/middle-Golgi proteins from the ER and the assembly of the cis/middle cisternae.

PGC-1α induces mitochondrial and myokine transcriptional programs and lipid droplet and glycogen accumulation in cultured human skeletal muscle cells

The transcriptional coactivator peroxisome proliferator-activated receptor-gamma coactivator 1 alpha (PGC-1α) is a chief activator of mitochondrial and metabolic programs and protects against atrophy in skeletal muscle (skm). Here we tested whether PGC-1α overexpression could restructure the transcriptome and metabolism of primary cultured human skm cells, which display a phenotype that resembles the atrophic phenotype. An oligonucleotide microarray analysis was used to reveal the effects of PGC-1α on the whole transcriptome. Fifty-three different genes showed altered expression in response to PGC-1α: 42 upregulated and 11 downregulated. The main gene ontologies (GO) associated with the upregulated genes were mitochondrial components and processes and this was linked with an increase in COX activity, an indicator of mitochondrial content. Furthermore, PGC-1α enhanced mitochondrial oxidation of palmitate and lactate to CO2, but not glucose oxidation. The other most significantly associated GOs for the upregulated genes were chemotaxis and cytokine activity, and several cytokines, including IL-8/CXCL8, CXCL6, CCL5 and CCL8, were within the most highly induced genes. Indeed, PGC-1α highly increased IL-8 cell protein content. The most upregulated gene was PVALB, which is related to calcium signaling. Potential metabolic regulators of fatty acid and glucose storage were among mainly regulated genes. The mRNA and protein level of FITM1/FIT1, which enhances the formation of lipid droplets, was raised by PGC-1α, while in oleate-incubated cells PGC-1α increased the number of smaller lipid droplets and modestly triglyceride levels, compared to controls. CALM1, the calcium-modulated δ subunit of phosphorylase kinase, was downregulated by PGC-1α, while glycogen phosphorylase was inactivated and glycogen storage was increased by PGC-1α. In conclusion, of the metabolic transcriptome deficiencies of cultured skm cells, PGC-1α rescued the expression of genes encoding mitochondrial proteins and FITM1. Several myokine genes, including IL-8 and CCL5, which are known to be constitutively expressed in human skm cells, were induced by PGC-1α.

On homology searches by protein Blast and the characterization of the age of genes

Background: It has been shown in a variety of organisms, including mammals, that genes that appeared recently in evolution, for example orphan genes, evolve faster than older genes. Low functional constraints at the time of origin of novel genes may explain these results. However, this observation has been recently attributed to an artifact caused by the inability of Blast to detect the fastest genes in different eukaryotic genomes. Distinguishing between these two possible explanations would be of great importance for any studies dealing with the taxon distribution of proteins and the origin of novel genes. Results: Here we used simulations of protein sequences to examine the capacity of Blast to detect proteins of diverse evolutionary rates in the different species of an eukaryotic phylogenetic tree that included metazoans, fungi and plants. We simulated the evolution of protein genes with the same evolutionary rates than those observed in functional mammalian genes and with among-site rate heterogeneity. Under these conditions, we found that only a very small percentage of simulated ancestral eukaryotic proteins was affected by the Blast artifact. We show that the good detectability of Blast is due to the heterogeneity of protein evolutionary rates at different sites, since only a small conserved motif in a sequence suffices to detect its homologues. Our results indicate that Blast, at least when applied within eukaryotes, only misses homologues of extremely fast-evolving sequences, which are rare in the mammalian genome, as well as sequences evolving homogeneously or pseudogenes.Conclusion: Although great care should be exercised in the recognition of remote homologues, most functional mammalian genes can be detected in eukaryotic genomes by Blast. That is, the majority of functional mammalian genes are not as fast as for not being detected in other metazoans, fungi or plants, if they had been present in these organisms. Thus, the correlation previously found between age and rate seems not to be due to a pure Blast artifact, at least for mammals. This may have important implications to understand the mechanisms by which novel genes originate.

Comparative analysis of amino acid repeats in rodents and humans

Amino acid tandem repeats, also called homopolymeric tracts, are extremely abundant in eukaryotic proteins. To gain insight into the genome-wide evolution of these regions in mammals, we analyzed the repeat content in a large data set of rat-mouse-human orthologs. Our results show that human proteins contain more amino acid repeats than rodent proteins and that trinucleotide repeats are also more abundant in human coding sequences. Using the human species as an outgroup, we were able to address differences in repeat loss and repeat gain in the rat and mouse lineages. In this data set, mouse proteins contain substantially more repeats than rat proteins, which can be at least partly attributed to a higher repeat loss in the rat lineage. The data are consistent with a role for trinucleotide slippage in the generation of novel amino acid repeats. We confirm the previously observed functional bias of proteins with repeats, with overrepresentation of transcription factors and DNA-binding proteins. We show that genes encoding amino acid repeats tend to have an unusually high GC content, and that differences in coding GC content among orthologs are directly related to the presence/absence of repeats. We propose that the different GC content isochore structure in rodents and humans may result in an increased amino acid repeat prevalence in the human lineage.

High pressure induced changes on sarcoplasmic protein fraction and quality indicators

The combined effect of pressure and mild temperature treatments on bovine sarcoplasmic proteins and quality parameters was assessed. M. longissimus dorsi samples were pressurised in a range of 200–600 MPa and 10–30 °C. High Pressure Processing (HPP) induced a reduction of protein solubility (p < 0.001) compared to non-treated controls (NT), more pronounced above 200 MPa. HPP at pressures higher than 200 MPa induced a strong modification (p < 0.001) of meat colour and a reduction of water holding capacity (WHC). SDS–PAGE analysis demonstrated that HPP significantly modified the composition of the sarcoplasmic protein fraction. The pressurisation temperature mainly affected protein solubility and colour; a smaller effect was observed on protein profiles. Significant correlations (p < 0.001) between sarcoplasmic protein solubility and both expressible moisture (r = −0.78) and colour parameters (r = −0.81 to −0.91) suggest that pressure induced denaturation of sarcoplasmic proteins could influence to some extent WHC and colour modifications of beef. Changes in protein band intensities were also significantly correlated with protein solubility, meat lightness and expressible moisture. These results describe the changes induced by HPP on sarcoplasmic proteins and confirm a relationship between modification of the sarcoplasmic protein fraction and alteration of meat quality characteristics