PredictProtein--an open resource for online prediction of protein structural and functional features

in Nucleic Acids Research (2014), 42(8),

PredictProtein is a meta-service for sequence analysis that has been predicting structural and functional features of proteins since 1992. Queried with a protein sequence it returns: multiple sequence ... [more ▼]

PredictProtein is a meta-service for sequence analysis that has been predicting structural and functional features of proteins since 1992. Queried with a protein sequence it returns: multiple sequence alignments, predicted aspects of structure (secondary structure, solvent accessibility, transmembrane helices (TMSEG) and strands, coiled-coil regions, disulfide bonds and disordered regions) and function. The service incorporates analysis methods for the identification of functional regions (ConSurf), homology-based inference of Gene Ontology terms (metastudent), comprehensive subcellular localization prediction (LocTree3), protein–protein binding sites (ISIS2), protein–polynucleotide binding sites (SomeNA) and predictions of the effect of point mutations (non-synonymous SNPs) on protein function (SNAP2). Our goal has always been to develop a system optimized to meet the demands of experimentalists not highly experienced in bioinformatics. To this end, the PredictProtein results are presented as both text and a series of intuitive, interactive and visually appealing figures. The web server and sources are available at http://ppopen.rostlab.org. [less ▲]

The FurA regulon in Anabaena sp. PCC 7120: in silico prediction and experimental validation of novel target genes.

in Nucleic acids research (2014), 42(8), 4833-46

In the filamentous cyanobacterium Anabaena sp. PCC 7120, the ferric uptake regulator FurA functions as a global transcriptional regulator. Despite several analyses have focused on elucidating the FurA ... [more ▼]

In the filamentous cyanobacterium Anabaena sp. PCC 7120, the ferric uptake regulator FurA functions as a global transcriptional regulator. Despite several analyses have focused on elucidating the FurA-regulatory network, the number of target genes described for this essential transcription factor is limited to a handful of examples. In this article, we combine an in silico genome-wide predictive approach with experimental determinations to better define the FurA regulon. Predicted FurA-binding sites were identified upstream of 215 genes belonging to diverse functional categories including iron homeostasis, photosynthesis and respiration, heterocyst differentiation, oxidative stress defence and light-dependent signal transduction mechanisms, among others. The probabilistic model proved to be effective at discerning FurA boxes from non-cognate sequences, while subsequent electrophoretic mobility shift assay experiments confirmed the in vitro specific binding of FurA to at least 20 selected predicted targets. Gene-expression analyses further supported the dual role of FurA as transcriptional modulator that can act both as repressor and as activator. In either role, the in vitro affinity of the protein to its target sequences is strongly dependent on metal co-regulator and reducing conditions, suggesting that FurA couples in vivo iron homeostasis and the response to oxidative stress to major physiological processes in cyanobacteria. [less ▲]

uORFdb--a comprehensive literature database on eukaryotic uORF biology.

in Nucleic acids research (2014), 42(1), 60-7

Approximately half of all human transcripts contain at least one upstream translational initiation site that precedes the main coding sequence (CDS) and gives rise to an upstream open reading frame (uORF ... [more ▼]

Approximately half of all human transcripts contain at least one upstream translational initiation site that precedes the main coding sequence (CDS) and gives rise to an upstream open reading frame (uORF). We generated uORFdb, publicly available at http://cbdm.mdc-berlin.de/tools/uorfdb, to serve as a comprehensive literature database on eukaryotic uORF biology. Upstream ORFs affect downstream translation by interfering with the unrestrained progression of ribosomes across the transcript leader sequence. Although the first uORF-related translational activity was observed >30 years ago, and an increasing number of studies link defective uORF-mediated translational control to the development of human diseases, the features that determine uORF-mediated regulation of downstream translation are not well understood. The uORFdb was manually curated from all uORF-related literature listed at the PubMed database. It categorizes individual publications by a variety of denominators including taxon, gene and type of study. Furthermore, the database can be filtered for multiple structural and functional uORF-related properties to allow convenient and targeted access to the complex field of eukaryotic uORF biology. [less ▲]

Predicting missing expression values in gene regulatory networks using a discrete logic modeling optimization guided by network stable states

in Nucleic Acids Research (2013), 41(1), 8

The development of new high-throughput technologies enables us to measure genome-wide transcription levels, protein abundance, metabolite concentration, etc. Nevertheless, these experimental data are ... [more ▼]

The development of new high-throughput technologies enables us to measure genome-wide transcription levels, protein abundance, metabolite concentration, etc. Nevertheless, these experimental data are often noisy and incomplete, which hinders data analysis, modeling and prediction. Here, we propose a method to predict expression values of genes involved in stable cellular phenotypes from the expression values of the remaining genes in a literature-based gene regulatory network. The consistency between predicted and known stable states from experimental data is used to guide an iterative network pruning that contextualizes the network to the biological conditions under which the expression data were obtained. Using the contextualized network and the property of network stability we predict gene expression values missing from experimental data. The prediction method assumes a Boolean model to compute steady states of networks and an evolutionary algorithm to iteratively prune the networks. The evolutionary algorithm samples the probability distribution of positive feedback loops or positive circuits and individual interactions within the subpopulation of the best-pruned networks at each iteration. The resulting expression inference is based not only on previous knowledge about local connectivity but also on a global network property (stability), providing robustness in the predictions. [less ▲]

Interplay of microRNAs, transcription factors and target genes: Linking dynamic expression changes to function

in Nucleic Acids Research (2013), 41(5), 2817-2831

MicroRNAs (miRNAs) are ubiquitously expressed small non-coding RNAs that, in most cases, negatively regulate gene expression at the posttranscriptional level. miRNAs are involved in fine-tuning ... [more ▼]

MicroRNAs (miRNAs) are ubiquitously expressed small non-coding RNAs that, in most cases, negatively regulate gene expression at the posttranscriptional level. miRNAs are involved in fine-tuning fundamental cellular processes such as proliferation, cell death and cell cycle control and are believed to confer robustness to biological responses. Here, we investigated simultaneously the transcriptional changes of miRNA and mRNA expression levels over time after activation of the Janus kinase/Signal transducer and activator of transcription (Jak/STAT) pathway by interferon-c stimulation of melanoma cells. To examine global miRNA and mRNA expression patterns, time-series microarray data were analysed. We observed delayed responses of miRNAs (after 24-48 h) with respect to mRNAs (12-24 h) and identified biological functions involved at each step of the cellular response. Inference of the upstream regulators allowed for identification of transcriptional regulators involved in cellular reactions to interferon-c stimulation. Linking expression profiles of transcriptional regulators and miRNAs with their annotated functions, we demonstrate the dynamic interplay of miRNAs and upstream regulators with biological functions. Finally, our data revealed network motifs in the form of feed-forward loops involving transcriptional regulators, mRNAs and miRNAs. Additional information obtained from integrating time-series mRNA and miRNA data may represent an important step towards understanding the regulatory principles of gene expression. © The Author(s) 2013. [less ▲]

The non-autonomous retrotransposon SVA is trans-mobilized by the human LINE-1 protein machinery

in Nucleic Acids Research (2012), 40(4), 1666-1683

SINE-VNTR-Alu (SVA) elements are non-autonomous, hominid-specific non-LTR retrotransposons and distinguished by their organization as composite mobile elements. They represent the evolutionarily youngest ... [more ▼]

SINE-VNTR-Alu (SVA) elements are non-autonomous, hominid-specific non-LTR retrotransposons and distinguished by their organization as composite mobile elements. They represent the evolutionarily youngest, currently active family of human non-LTR retrotransposons, and sporadically generate disease-causing insertions. Since preexisting, genomic SVA sequences are characterized by structural hallmarks of Long Interspersed Elements 1 (LINE-1, L1)-mediated retrotransposition, it has been hypothesized for several years that SVA elements are mobilized by the L1 protein machinery in trans. To test this hypothesis, we developed an SVA retrotransposition reporter assay in cell culture using three different human-specific SVA reporter elements. We demonstrate that SVA elements are mobilized in HeLa cells only in the presence of both L1-encoded proteins, ORF1p and ORF2p. SVA trans-mobilization rates exceeded pseudogene formation frequencies by 12-to 300-fold in HeLa-HA cells, indicating that SVA elements represent a preferred substrate for L1 proteins. Acquisition of an AluSp element increased the trans-mobilization frequency of the SVA reporter element by ∼25-fold. Deletion of (CCCTCT)n repeats and Alu-like region of a canonical SVA reporter element caused significant attenuation of the SVA trans-mobilization rate. SVA de novo insertions were predominantly full-length, occurred preferentially in G+C-rich regions, and displayed all features of L1-mediated retrotransposition which are also observed in preexisting genomic SVA insertions. © 2012 The Author(s). [less ▲]

Hierarchical regulation of the NikR-mediated nickel response in Helicobacter pylori.

in Nucleic acids research (2011), 39(17), 7564-75

Nickel is an essential metal for Helicobacter pylori, as it is the co-factor of two enzymes crucial for colonization, urease and hydrogenase. Nickel is taken up by specific transporters and its ... [more ▼]

Nickel is an essential metal for Helicobacter pylori, as it is the co-factor of two enzymes crucial for colonization, urease and hydrogenase. Nickel is taken up by specific transporters and its intracellular homeostasis depends on nickel-binding proteins to avoid toxicity. Nickel trafficking is controlled by the Ni(II)-dependent transcriptional regulator NikR. In contrast to other NikR proteins, NikR from H. pylori is a pleiotropic regulator that depending on the target gene acts as an activator or a repressor. We systematically quantified the in vivo Ni(2+)-NikR response of 11 direct NikR targets that encode functions related to nickel metabolism, four activated and seven repressed genes. Among these, four targets were characterized for the first time (hpn, hpn-like, hydA and hspA) and NikR binding to their promoter regions was demonstrated by electrophoretic mobility shift assays. We found that NikR-dependent repression was generally set up at higher nickel concentrations than activation. Kinetics of the regulation revealed a gradual and temporal NikR-mediated response to nickel where activation of nickel-protection mechanisms takes place before repression of nickel uptake. Our in vivo study demonstrates, for the first time, a chronological hierarchy in the NikR-dependent transcriptional response to nickel that is coherent with the control of nickel homeostasis in H. pylori. [less ▲]

A simple and efficient method to search for selected primary transcripts: non-coding and antisense RNAs in the human pathogen Enterococcus faecalis.

in Nucleic acids research (2011), 39(7), 46

Enterococcus faecalis is a commensal bacterium and a major opportunistic human pathogen. In this study, we combined in silico predictions with a novel 5'RACE-derivative method coined '5'tagRACE', to ... [more ▼]

Enterococcus faecalis is a commensal bacterium and a major opportunistic human pathogen. In this study, we combined in silico predictions with a novel 5'RACE-derivative method coined '5'tagRACE', to perform the first search for non-coding RNAs (ncRNAs) encoded on the E. faecalis chromosome. We used the 5'tagRACE to simultaneously probe and characterize primary transcripts, and demonstrate here the simplicity, the reliability and the sensitivity of the method. The 5'tagRACE is complementary to tiling arrays or RNA-sequencing methods, and is also directly applicable to deep RNA sequencing and should significantly improve functional studies of bacterial RNA landscapes. From 45 selected loci of the E. faecalis chromosome, we discovered and mapped 29 novel ncRNAs, 10 putative novel mRNAs and 16 antisense transcriptional organizations. We describe in more detail the oxygen-dependent expression of one ncRNA located in an E. faecalis pathogenicity island, the existence of an ncRNA that is antisense to the ncRNA modulator of the RNA polymerase, SsrS and provide evidences for the functional interplay between two distinct toxin-antitoxin modules. [less ▲]

Structural and mechanistic insights into Helicobacter pylori NikR activation.

in Nucleic acids research (2010), 38(9), 3106-18

NikR is a transcriptional metalloregulator central in the mandatory response to acidity of Helicobacter pylori that controls the expression of numerous genes by binding to specific promoter regions. NikR ... [more ▼]

NikR is a transcriptional metalloregulator central in the mandatory response to acidity of Helicobacter pylori that controls the expression of numerous genes by binding to specific promoter regions. NikR/DNA interactions were proposed to rely on protein activation by Ni(II) binding to high-affinity (HA) and possibly secondary external (X) sites. We describe a biochemical characterization of HpNikR mutants that shows that the HA sites are essential but not sufficient for DNA binding, while the secondary external (X) sites and residues from the HpNikR dimer-dimer interface are important for DNA binding. We show that a second metal is necessary for HpNikR/DNA binding, but only to some promoters. Small-angle X-ray scattering shows that HpNikR adopts a defined conformation in solution, resembling the cis-conformation and suggests that nickel does not trigger large conformational changes in HpNikR. The crystal structures of selected mutants identify the effects of each mutation on HpNikR structure. This study unravels key structural features from which we derive a model for HpNikR activation where: (i) HA sites and an hydrogen bond network are required for DNA binding and (ii) metallation of a unique secondary external site (X) modulates HpNikR DNA binding to low-affinity promoters by disruption of a salt bridge. [less ▲]

Sequence-structure relationships in RNA loops: establishing the basis for loop homology modeling.

in Nucleic Acids Research (2010), 38(3), 970-80

The specific function of RNA molecules frequently resides in their seemingly unstructured loop regions. We performed a systematic analysis of RNA loops extracted from experimentally determined three ... [more ▼]

The specific function of RNA molecules frequently resides in their seemingly unstructured loop regions. We performed a systematic analysis of RNA loops extracted from experimentally determined three-dimensional structures of RNA molecules. A comprehensive loop-structure data set was created and organized into distinct clusters based on structural and sequence similarity. We detected clear evidence of the hallmark of homology present in the sequence-structure relationships in loops. Loops differing by <25% in sequence identity fold into very similar structures. Thus, our results support the application of homology modeling for RNA loop model building. We established a threshold that may guide the sequence divergence-based selection of template structures for RNA loop homology modeling. Of all possible sequences that are, under the assumption of isosteric relationships, theoretically compatible with actual sequences observed in RNA structures, only a small fraction is contained in the Rfam database of RNA sequences and classes implying that the actual RNA loop space may consist of a limited number of unique loop structures and conserved sequences. The loop-structure data sets are made available via an online database, RLooM. RLooM also offers functionalities for the modeling of RNA loop structures in support of RNA engineering and design efforts. [less ▲]