Nucleic Acids Research - Advance Access

The Ontology Lookup Service: more data and better tools for controlled vocabulary queries

Cote, R. G., Jones, P., Martens, L., Apweiler, R., Hermjakob, H. — 2008-05-08

The Ontology Lookup Service (OLS) (http://www.ebi.ac.uk/ols) provides interactive and programmatic interfaces to query, browse and navigate an ever increasing number of biomedical ontologies and controlled vocabularies. The volume of data available for querying has more than quadrupled since it went into production and OLS functionality has been integrated into several high-usage databases and data entry tools. Improvements have been made to both OLS query interfaces, based on user feedback and requirements, to improve usability and service interoperability and provide novel ways to perform queries.

PaLS: filtering common literature, biological terms and pathway information

Alibes, A., Canada, A., Diaz-Uriarte, R. — 2008-05-08

Many biological experiments and their subsequent analysis yield lists of genes or proteins that can potentially be important to the prognosis or diagnosis of certain diseases (e.g. cancer). Nowadays, information about the function of those genes or proteins may be already gathered in some databases, but it is essential to understand if some of the members of those lists have a function in common or if they belong to the same metabolic pathway. To help researchers filter those genes or proteins that have such information in common, we have developed PaLS (pathway and literature strainer, http://pals.bioinfo.cnio.es). PaLS takes a list or a set of lists of gene or protein identifiers and shows which ones share certain descriptors. Four publicly available databases have been used for this purpose: PubMed, which links genes with those articles that make reference to them; Gene Ontology, an annotated ontology of terms related to the cellular component, biological process or molecular function where those genes or proteins are involved; KEGG pathways and Reactome pathways. Those descriptors among these four sources of information that are shared by more members of the list (or lists) are highlighted by PaLS.

DNA binding of dinuclear iron(II) metallosupramolecular cylinders. DNA unwinding and sequence preference

Malina, J., Hannon, M. J., Brabec, V. — 2008-05-08

[Fe₂L₃]⁴⁺ (L = C₂₅H₂₀N₄) is a synthetic tetracationic supramolecular cylinder (with a triple helical architecture) that targets the major groove of DNA and can bind to DNA Y-shaped junctions. To explore the DNA-binding mode of [Fe₂L₃]⁴⁺, we examine herein the interactions of pure enantiomers of this cylinder with DNA by biochemical and molecular biology methods. The results have revealed that, in addition to the previously reported bending of DNA, the enantiomers extensively unwind DNA, with the M enantiomer being the more efficient at unwinding, and exhibit preferential binding to regular alternating purine–pyrimidine sequences, with the M enantiomer showing a greater preference. Also, interestingly, the DNA binding of bulky cylinders [Fe₂(L-CF₃)₃]⁴⁺ and [Fe₂(L-Ph)₃]⁴⁺ results in no DNA unwinding and also no sequence preference of their DNA binding was observed. The observation of sequence-preference in the binding of these supramolecular cylinders suggests that a concept based on the use of metallosupramolecular cylinders might result in molecular designs that recognize the genetic code in a sequence-dependent manner with a potential ability to affect the processing of the genetic code.

MultiBind and MAPPIS: webservers for multiple alignment of protein 3D-binding sites and their interactions

Shulman-Peleg, A., Shatsky, M., Nussinov, R., Wolfson, H. J. — 2008-05-08

Analysis of protein–ligand complexes and recognition of spatially conserved physico-chemical properties is important for the prediction of binding and function. Here, we present two webservers for multiple alignment and recognition of binding patterns shared by a set of protein structures. The first webserver, MultiBind (http://bioinfo3d.cs.tau.ac.il/MultiBind), performs multiple alignment of protein binding sites. It recognizes the common spatial chemical binding patterns even in the absence of similarity of the sequences or the folds of the compared proteins. The input to the MultiBind server is a set of protein-binding sites defined by interactions with small molecules. The output is a detailed list of the shared physico-chemical binding site properties. The second webserver, MAPPIS (http://bioinfo3d.cs.tau.ac.il/MAPPIS), aims to analyze protein–protein interactions. It performs multiple alignment of protein–protein interfaces (PPIs), which are regions of interaction between two protein molecules. MAPPIS recognizes the spatially conserved physico-chemical interactions, which often involve energetically important hot-spot residues that are crucial for protein–protein associations. The input to the MAPPIS server is a set of protein-protein complexes. The output is a detailed list of the shared interaction properties of the interfaces.

TRE-dependent transcription activation by JDP2-CHOP10 association

Weidenfeld-Baranboim, K., Bitton-Worms, K., Aronheim, A. — 2008-05-07

The c-Jun dimerization protein 2, JDP2, is a member of the activating protein 1 (AP-1) family of transcription factors. Overexpression of JDP2 has been shown to result in repression of AP-1-dependent transcription and inhibition of cellular transformation. Other studies suggested that JDP2 may function as an oncogene. Here we describe the identification of CHOP10, a member of the CCAAT enhancer binding proteins, as a protein associating with JDP2. In contrast to the inhibition of transcription by JDP2, JDP2–CHOP complex strongly enhances transcription from promoters containing TPA response elements (TRE), but not from those containing cyclic AMP response elements (CRE). The association between JDP2 and CHOP10 involves the leucine zipper motifs of both proteins, whereas, the basic domain of CHOP10 contributes to the association of the JDP2–CHOP10 complex with the DNA. DNA binding of JDP2–CHOP complex is observed both in vitro and in vivo. Finally, overexpression of JDP2 results in increased cell viability following ER stress and counteracts CHOP10 pro-apoptotic activity. JDP2 expression may determine the threshold for cell sensitivity to ER stress. This is the first report describing TRE-dependent activation of transcription by JDP2 and thus may provide an explanation for the as yet unexplored oncogenic properties of JDP2.

OREST: the online resource for EST analysis

Waegele, B., Schmidt, T., Mewes, H. W., Ruepp, A. — 2008-05-07

The generation of expressed sequence tag (EST) libraries offers an affordable approach to investigate organisms, if no genome sequence is available. OREST (http://mips.gsf.de/genre/proj/orest/index.html) is a server-based EST analysis pipeline, which allows the rapid analysis of large amounts of ESTs or cDNAs from mammalia and fungi. In order to assign the ESTs to genes or proteins OREST maps DNA sequences to reference datasets of gene products and in a second step to complete genome sequences. Mapping against genome sequences recovers additional 13% of EST data, which otherwise would escape further analysis. To enable functional analysis of the datasets, ESTs are functionally annotated using the hierarchical FunCat annotation scheme as well as GO annotation terms. OREST also allows to predict the association of gene products and diseases by Morbid Map (OMIM) classification. A statistical analysis of the results of the dataset is possible with the included PROMPT software, which provides information about enrichment and depletion of functional and disease annotation terms. OREST was successfully applied for the identification and functional characterization of more than 3000 EST sequences of the common marmoset monkey (Callithrix jacchus) as part of an international collaboration.

The Jpred 3 secondary structure prediction server

Cole, C., Barber, J. D., Barton, G. J. — 2008-05-07

Jpred (http://www.compbio.dundee.ac.uk/jpred) is a secondary structure prediction server powered by the Jnet algorithm. Jpred performs over 1000 predictions per week for users in more than 50 countries. The recently updated Jnet algorithm provides a three-state (-helix, β-strand and coil) prediction of secondary structure at an accuracy of 81.5%. Given either a single protein sequence or a multiple sequence alignment, Jpred derives alignment profiles from which predictions of secondary structure and solvent accessibility are made. The predictions are presented as coloured HTML, plain text, PostScript, PDF and via the Jalview alignment editor to allow flexibility in viewing and applying the data. The new Jpred 3 server includes significant usability improvements that include clearer feedback of the progress or failure of submitted requests. Functional improvements include batch submission of sequences, summary results via email and updates to the search databases. A new software pipeline will enable Jnet/Jpred to continue to be updated in sync with major updates to SCOP and UniProt and so ensures that Jpred 3 will maintain high-accuracy predictions.

SEQATOMS: a web tool for identifying missing regions in PDB in sequence context

Brandt, B. W., Heringa, J., Leunissen, J. A. M. — 2008-05-07

With over 46 000 proteins, the Protein Data Bank (PDB) is the most important database with structural information of biological macromolecules. PDB files contain sequence and coordinate information. Residues present in the sequence can be absent from the coordinate section, which means their position in space is unknown. Similarity searches are routinely carried out against sequences taken from PDB SEQRES. However, there no distinction is made between residues that have a known or unknown position in the 3D protein structure. We present a FASTA sequence database that is produced by combining the sequence and coordinate information. All residues absent from the PDB coordinate section are masked with lower-case letters, thereby providing a view of these residues in the context of the entire protein sequence, which facilitates inspecting ‘missing’ regions. We also provide a masked version of the CATH domain database. A user-friendly BLAST interface is available for similarity searching. In contrast to standard (stand-alone) BLAST output, which only contains upper-case letters, our output retains the lower-case letters of the masked regions. Thus, our server can be used to perform BLAST searching case-sensitively. Here, we have applied it to the study of missing regions in their sequence context. SEQATOMS is available at http://www.bioinformatics.nl/tools/seqatoms/.

Microarray retriever: a web-based tool for searching and large scale retrieval of public microarray data

2008-05-07

The major public microarray repositories Gene Expression Omnibus and ArrayExpress are growing rapidly. This enables meta-analysis studies, in which expression data from multiple individual studies are combined. To facilitate these types of studies, we developed Microarray Retriever for searching and retrieval of data from GEO and ArrayExpress. The tool allows access to the two repositories simultaneously, to search in the repositories using complex queries, to retrieve microarray data for published articles and to download data in one structured archive. The tool is available on the web at: http://www.lgtc.nl/MaRe/

Molecular basis for the lack of enantioselectivity of human 3-phosphoglycerate kinase

2008-05-07

Non-natural l-nucleoside analogues are increasingly used as therapeutic agents to treat cancer and viral infections. To be active, l-nucleosides need to be phosphorylated to their respective triphosphate metabolites. This stepwise phosphorylation relies on human enzymes capable of processing l-nucleoside enantiomers. We used crystallographic analysis to reveal the molecular basis for the low enantioselectivity and the broad specificity of human 3-phosphoglycerate kinase (hPGK), an enzyme responsible for the last step of phosphorylation of many nucleotide derivatives. Based on structures of hPGK in the absence of nucleotides, and bound to l and d forms of MgADP and MgCDP, we show that a non-specific hydrophobic clamp to the nucleotide base, as well as a water-filled cavity behind it, allows high flexibility in the interaction between PGK and the bases. This, combined with the dispensability of hydrogen bonds to the sugar moiety, and ionic interactions with the phosphate groups, results in the positioning of different nucleotides so to expose their diphosphate group in a position competent for catalysis. Since the third phosphorylation step is often rate limiting, our results are expected to alleviate in silico tailoring of l-type prodrugs to assure their efficient metabolic processing.

NetMHC-3.0: accurate web accessible predictions of human, mouse and monkey MHC class I affinities for peptides of length 8-11

Lundegaard, C., Lamberth, K., Harndahl, M., Buus, S., Lund, O., Nielsen, M. — 2008-05-07

NetMHC-3.0 is trained on a large number of quantitative peptide data using both affinity data from the Immune Epitope Database and Analysis Resource (IEDB) and elution data from SYFPEITHI. The method generates high-accuracy predictions of major histocompatibility complex (MHC): peptide binding. The predictions are based on artificial neural networks trained on data from 55 MHC alleles (43 Human and 12 non-human), and position-specific scoring matrices (PSSMs) for additional 67 HLA alleles. As only the MHC class I prediction server is available, predictions are possible for peptides of length 8–11 for all 122 alleles. artificial neural network predictions are given as actual IC₅₀ values whereas PSSM predictions are given as a log-odds likelihood scores. The output is optionally available as download for easy post-processing. The training method underlying the server is the best available, and has been used to predict possible MHC-binding peptides in a series of pathogen viral proteomes including SARS, Influenza and HIV, resulting in an average of 75–80% confirmed MHC binders. Here, the performance is further validated and benchmarked using a large set of newly published affinity data, non-redundant to the training set. The server is free of use and available at: http://www.cbs.dtu.dk/services/NetMHC.

FFPred: an integrated feature-based function prediction server for vertebrate proteomes

Lobley, A. E., Nugent, T., Orengo, C. A., Jones, D. T. — 2008-05-07

One of the challenges of the post-genomic era is to provide accurate function annotations for large volumes of data resulting from genome sequencing projects. Most function prediction servers utilize methods that transfer existing database annotations between orthologous sequences. In contrast, there are few methods that are independent of homology and can annotate distant and orphan protein sequences. The FFPred server adopts a machine-learning approach to perform function prediction in protein feature space using feature characteristics predicted from amino acid sequence. The features are scanned against a library of support vector machines representing over 300 Gene Ontology (GO) classes and probabilistic confidence scores returned for each annotation term. The GO term library has been modelled on human protein annotations; however, benchmark performance testing showed robust performance across higher eukaryotes. FFPred offers important advantages over traditional function prediction servers in its ability to annotate distant homologues and orphan protein sequences, and achieves greater coverage and classification accuracy than other feature-based prediction servers. A user may upload an amino acid and receive annotation predictions via email. Feature information is provided as easy to interpret graphics displayed on the sequence of interest, allowing for back-interpretation of the associations between features and function classes.

Cooperation between EZH2, NSPc1-mediated histone H2A ubiquitination and Dnmt1 in HOX gene silencing

Wu, X., Gong, Y., Yue, J., Qiang, B., Yuan, J., Peng, X. — 2008-05-06

An intricate interplay between DNA methylation and polycomb-mediated gene silencing has been highlighted recently. Here we provided evidence that Nervous System Polycomb 1 (NSPc1), a BMI1 homologous polycomb protein, plays important roles in promoting H2A ubiquitination and cooperates with DNA methylation in HOX gene silencing. We showed that NSPc1 stimulates H2A ubiquitination in vivo and in vitro through direct interaction with both RING2 and H2A. RT-PCR analysis revealed that loss of NSPc1, EZH2 or DNA methyltransferase 1 (Dnmt1), or inhibition of DNA methylation in HeLa cells de-represses the expression of HOXA7. Chromatin immunoprecipitation (ChIP) assays demonstrated that NSPc1, EZH2 and Dnmt1 bind to the promoter of HOXA7, which is frequently hypermethylated in tumors. Knockdown of NSPc1 results in significant reduction of H2A ubiquitination and DNA demethylation as well as Dnmt1 dissociation in the HOXA7 promoter. Meanwhile Dnmt1 deficiency affects NSPc1 recruitment and H2A ubiquitination, whereas on both cases EZH2-mediated H3K27 trimethylation remains unaffected. When EZH2 was depleted, however, NSPc1 and Dnmt1 enrichment was abolished concomitant with local reduction of H3K27 trimethylation, H2A ubiquitination and DNA methylation. Taken together, our findings indicated that NSPc1-mediated H2A ubiquitination and DNA methylation, both being directed by EZH2, are interdependent in long-term target gene silencing within cancer cells.

ProfCom: a web tool for profiling the complex functionality of gene groups identified from high-throughput data

Antonov, A. V., Schmidt, T., Wang, Y., Mewes, H. W. — 2008-05-06

ProfCom is a web-based tool for the functional interpretation of a gene list that was identified to be related by experiments. A trait which makes ProfCom a unique tool is an ability to profile enrichments of not only available Gene Ontology (GO) terms but also of ‘complex functions’. A ‘Complex function’ is constructed as Boolean combination of available GO terms. The complex functions inferred by ProfCom are more specific in comparison to single terms and describe more accurately the functional role of genes. ProfCom provides a user friendly dialog-driven web page submission available for several model organisms and supports most available gene identifiers. In addition, the web service interface allows the submission of any kind of annotation data. ProfCom is freely available at http://webclu.bio.wzw.tum.de/profcom/.

GraphWeb: mining heterogeneous biological networks for gene modules with functional significance

Reimand, J., Tooming, L., Peterson, H., Adler, P., Vilo, J. — 2008-05-06

Deciphering heterogeneous cellular networks with embedded modules is a great challenge of current systems biology. Experimental and computational studies construct complex networks of molecules that describe various aspects of the cell such as transcriptional regulation, protein interactions and metabolism. Groups of interacting genes and proteins reflect network modules that potentially share regulatory mechanisms and relate to common function. Here, we present GraphWeb, a public web server for biological network analysis and module discovery. GraphWeb provides methods to: (1) integrate heterogeneous and multispecies data for constructing directed and undirected, weighted and unweighted networks; (ii) discover network modules using a variety of algorithms and topological filters and (iii) interpret modules using functional knowledge of the Gene Ontology and pathways, as well as regulatory features such as binding motifs and microRNA targets. GraphWeb is designed to analyse individual or multiple merged networks, search for conserved features across multiple species, mine large biological networks for smaller modules, discover novel candidates and connections for known pathways and compare results of high-throughput datasets. The GraphWeb is available at http://biit.cs.ut.ee/graphweb/.

Selection of a novel class of RNA-RNA interaction motifs based on the ligase ribozyme with defined modular architecture

Ohuchi, S. P., Ikawa, Y., Nakamura, Y. — 2008-05-06

To develop molecular tools for the detection and control of RNA molecules whose functions rely on their 3D structures, we have devised a selection system to isolate novel RNA motifs that interact with a target RNA structure within a given structural context. In this system, a GAAA tetraloop and its specific receptor motif (11-ntR) from an artificial RNA ligase ribozyme with modular architecture (the DSL ribozyme) were replaced with a target structure and random sequence, respectively. Motifs recognizing the target structure can be identified by in vitro selection based on ribozyme activity. A model selection targeting GAAA-loop successfully identified motifs previously known as GAAA-loop receptors. In addition, a new selection targeting a C-loop motif also generated novel motifs that interact with this structure. Biochemical analysis of one of the C-loop receptor motifs revealed that it could also function as an independent structural unit.

RAPIDO: a web server for the alignment of protein structures in the presence of conformational changes

Mosca, R., Schneider, T. R. — 2008-05-06

Rapid alignment of proteins in terms of domains (RAPIDO) is a web server for the 3D alignment of crystal structures of different protein molecules in the presence of conformational change. The structural alignment algorithm identifies groups of equivalent atoms whose interatomic distances are constant (within a defined tolerance) in the two structures being compared and considers these groups of atoms as rigid bodies. In addition to the functionalities provided by existing tools, RAPIDO can identify structurally equivalent regions also when these consist of fragments that are distant in terms of sequence and separated by other movable domains. Furthermore, RAPIDO takes the variation in the reliability of atomic coordinates into account in the comparison of distances between equivalent atoms by employing weighting-functions based on the refined B-values. The regions identified as equivalent by RAPIDO furnish reliable sets of residues for the superposition of the two structures for subsequent detailed analysis. The RAPIDO server, with related documentation, is available at http://webapps.embl-hamburg.de/rapido.

Phylogeny.fr: robust phylogenetic analysis for the non-specialist

2008-05-06

Phylogenetic analyses are central to many research areas in biology and typically involve the identification of homologous sequences, their multiple alignment, the phylogenetic reconstruction and the graphical representation of the inferred tree. The Phylogeny.fr platform transparently chains programs to automatically perform these tasks. It is primarily designed for biologists with no experience in phylogeny, but can also meet the needs of specialists; the first ones will find up-to-date tools chained in a phylogeny pipeline to analyze their data in a simple and robust way, while the specialists will be able to easily build and run sophisticated analyses. Phylogeny.fr offers three main modes. The ‘One Click’ mode targets non-specialists and provides a ready-to-use pipeline chaining programs with recognized accuracy and speed: MUSCLE for multiple alignment, PhyML for tree building, and TreeDyn for tree rendering. All parameters are set up to suit most studies, and users only have to provide their input sequences to obtain a ready-to-print tree. The ‘Advanced’ mode uses the same pipeline but allows the parameters of each program to be customized by users. The ‘A la Carte’ mode offers more flexibility and sophistication, as users can build their own pipeline by selecting and setting up the required steps from a large choice of tools to suit their specific needs. Prior to phylogenetic analysis, users can also collect neighbors of a query sequence by running BLAST on general or specialized databases. A guide tree then helps to select neighbor sequences to be used as input for the phylogeny pipeline. Phylogeny.fr is available at: http://www.phylogeny.fr/

Crystal structure of trioxacarcin A covalently bound to DNA

Pfoh, R., Laatsch, H., Sheldrick, G. M. — 2008-05-03

We report a crystal structure that shows an antibiotic that extracts a nucleobase from a DNA molecule ‘caught in the act’ after forming a covalent bond but before departing with the base. The structure of trioxacarcin A covalently bound to double-stranded d(AACCGGTT) was determined to 1.78 Å resolution by MAD phasing employing brominated oligonucleotides. The DNA–drug complex has a unique structure that combines alkylation (at the N7 position of a guanine), intercalation (on the 3'-side of the alkylated guanine), and base flip-out. An antibiotic-induced flipping-out of a single, nonterminal nucleobase from a DNA duplex was observed for the first time in a crystal structure.

Structural and functional insights into human Tudor-SN, a key component linking RNA interference and editing

Li, C.-L., Yang, W.-Z., Chen, Y.-P., Yuan, H. S. — 2008-05-03

Human Tudor-SN is involved in the degradation of hyper-edited inosine-containing microRNA precursors, thus linking the pathways of RNA interference and editing. Tudor-SN contains four tandem repeats of staphylococcal nuclease-like domains (SN1–SN4) followed by a tudor and C-terminal SN domain (SN5). Here, we showed that Tudor-SN requires tandem repeats of SN domains for its RNA binding and cleavage activity. The crystal structure of a 64-kD truncated form of human Tudor-SN further shows that the four domains, SN3, SN4, tudor and SN5, assemble into a crescent-shaped structure. A concave basic surface formed jointly by SN3 and SN4 domains is likely involved in RNA binding, where citrate ions are bound at the putative RNase active sites. Additional modeling studies provide a structural basis for Tudor-SN's preference in cleaving RNA containing multiple I·U wobble-paired sequences. Collectively, these results suggest that tandem repeats of SN domains in Tudor-SN function as a clamp to capture RNA substrates.

ConTra: a promoter alignment analysis tool for identification of transcription factor binding sites across species

Hooghe, B., Hulpiau, P., van Roy, F., De Bleser, P. — 2008-05-03

Transcription factors (TFs) are key components in signaling pathways, and the presence of their binding sites in the promoter regions of DNA is essential for their regulation of the expression of the corresponding genes. Orthologous promoter sequences are commonly used to increase the specificity with which potentially functional transcription factor binding sites (TFBSs) are recognized and to detect possibly important similarities or differences between the different species. The ConTra (conserved TFBSs) web server provides the biologist at the bench with a user-friendly tool to interactively visualize TFBSs predicted using either TransFac (1) or JASPAR (2) position weight matrix libraries, on a promoter alignment of choice. The visualization can be preceded by a simple scoring analysis to explore which TFs are the most likely to bind to the promoter of interest. The ConTra web server is available at http://bioit.dmbr.ugent.be/ConTra/index.php.

Spatial effects on the speed and reliability of protein-DNA search

Wunderlich, Z., Mirny, L. A. — 2008-05-03

Strong experimental and theoretical evidence shows that transcription factors (TFs) and other specific DNA-binding proteins find their sites using a two-mode search: alternating between three-dimensional (3D) diffusion through the cell and one-dimensional (1D) sliding along the DNA. We show that, due to the 1D component of the search process, the search time of a TF can depend on the initial position of the TF. We formalize this effect by discriminating between two types of searches: global and local. Using analytical calculations and simulations, we estimate how close a TF and binding site need to be to make a local search likely. We then use our model to interpret the wide range of experimental measurements of this parameter. We also show that local and global searches differ significantly in average search time and the variability of search time. These results lead to a number of biological implications, including suggestions of how prokaryotes achieve rapid gene regulation and the relationship between the search mechanism and noise in gene expression. Lastly, we propose a number of experiments to verify the existence and quantify the extent of spatial effects on the TF search process in prokaryotes.

Purine twisted-intercalating nucleic acids: a new class of anti-gene molecules resistant to potassium-induced aggregation

2008-05-02

Sequence-specific targeting of genomic DNA by triplex-forming oligonucleotides (TFOs) is a promising strategy to modulate in vivo gene expression. Triplex formation involving G-rich oligonucleotides as third strand is, however, strongly inhibited by potassium-induced TFO self-association into G-quartet structures. We report here that G-rich TFOs with bulge insertions of (R)-1-O-[4-(1-pyrenylethynyl)-phenylmethyl] glycerol (called twisted intercalating nucleic acids, TINA) show a much lower tendency to aggregate in potassium than wild-type analogues do. We designed purine-motif TINA–TFOs for binding to a regulatory polypurine-polypyrimidine (pur/pyr) motif present in the promoter of the KRAS proto-oncogene. The binding of TINA–TFOs to the KRAS target has been analysed by electrophoresis mobility shift assays and DNase I footprinting experiments. We discovered that in the presence of potassium the wild-type TFOs did not bind to the KRAS target, differently from the TINA analogues, whose binding was observed up to 140 mM KCl. The designed TINA–TFOs were found to abrogate the formation of a DNA–protein complex at the pur/pyr site and to down-regulate the transcription of CAT driven by the murine KRAS promoter. Molecular modelling of the DNA/TINA–TFO triplexes are also reported. This study provides a new and promising approach to create TFOs to target in vivo the genome.

OGtree: a tool for creating genome trees of prokaryotes based on overlapping genes

Jiang, L.-W., Lin, K.-L., Lu, C. L. — 2008-05-02

OGtree is a web-based tool for constructing genome trees of prokaryotic species based on a measure of combining overlapping-gene content and overlapping-gene order in their whole genomes. The overlapping genes (OGs) are defined as adjacent genes whose coding sequences overlap partially or entirely. In fact, OGs are ubiquitous in microbial genomes and more conserved between species than non-OGs. Based on these properties, it has been suggested that OGs can serve as better phylogenetic characters than non-OGs for reconstructing the evolutionary relationships among microbial genomes. OGtree takes the accession numbers of prokaryotic genomes as its input. It then downloads their complete genomes from the National Centre for Biotechnology Information and identifies OGs in each genome and their orthologous OGs in other genomes. Next, OGtree computes an overlapping-gene distance between each pair of input genomes based on a combination of their OG content and orthologous OG order. Finally, it utilizes distance-based methods of building tree to reconstruct the genome trees of input prokaryotic genomes according to their pairwise OG distance. OGtree is available online at http://bioalgorithm.life.nctu.edu.tw/OGtree/.

Ribosomal RNAs are tolerant toward genetic insertions: evolutionary origin of the expansion segments

Yokoyama, T., Suzuki, T. — 2008-05-02

Ribosomal RNAs (rRNAs), assisted by ribosomal proteins, form the basic structure of the ribosome, and play critical roles in protein synthesis. Compared to prokaryotic ribosomes, eukaryotic ribosomes contain elongated rRNAs with several expansion segments and larger numbers of ribosomal proteins. To investigate architectural evolution and functional capability of rRNAs, we employed a Tn5 transposon system to develop a systematic genetic insertion of an RNA segment 31 nt in length into Escherichia coli rRNAs. From the plasmid library harboring a single rRNA operon containing random insertions, we isolated surviving clones bearing rRNAs with functional insertions that enabled rescue of the E. coli strain (7rrn) in which all chromosomal rRNA operons were depleted. We identified 51 sites with functional insertions, 16 sites in 16S rRNA and 35 sites in 23S rRNA, revealing the architecture of E. coli rRNAs to be substantially flexible. Most of the insertion sites show clear tendency to coincide with the regions of the expansion segments found in eukaryotic rRNAs, implying that eukaryotic rRNAs evolved from prokaryotic rRNAs suffering genetic insertions and selections.

Structural and evolutionary classification of Type II restriction enzymes based on theoretical and experimental analyses

Orlowski, J., Bujnicki, J. M. — 2008-05-02

For a very long time, Type II restriction enzymes (REases) have been a paradigm of ORFans: proteins with no detectable similarity to each other and to any other protein in the database, despite common cellular and biochemical function. Crystallographic analyses published until January 2008 provided high-resolution structures for only 28 of 1637 Type II REase sequences available in the Restriction Enzyme database (REBASE). Among these structures, all but two possess catalytic domains with the common PD-(D/E)XK nuclease fold. Two structures are unrelated to the others: R.BfiI exhibits the phospholipase D (PLD) fold, while R.PabI has a new fold termed ‘half-pipe’. Thus far, bioinformatic studies supported by site-directed mutagenesis have extended the number of tentatively assigned REase folds to five (now including also GIY-YIG and HNH folds identified earlier in homing endonucleases) and provided structural predictions for dozens of REase sequences without experimentally solved structures. Here, we present a comprehensive study of all Type II REase sequences available in REBASE together with their homologs detectable in the nonredundant and environmental samples databases at the NCBI. We present the summary and critical evaluation of structural assignments and predictions reported earlier, new classification of all REase sequences into families, domain architecture analysis and new predictions of three-dimensional folds. Among 289 experimentally characterized (not putative) Type II REases, whose apparently full-length sequences are available in REBASE, we assign 199 (69%) to contain the PD-(D/E)XK domain. The HNH domain is the second most common, with 24 (8%) members. When putative REases are taken into account, the fraction of PD-(D/E)XK and HNH folds changes to 48% and 30%, respectively. Fifty-six characterized (and 521 predicted) REases remain unassigned to any of the five REase folds identified so far, and may exhibit new architectures. These enzymes are proposed as the most interesting targets for structure determination by high-resolution experimental methods. Our analysis provides the first comprehensive map of sequence-structure relationships among Type II REases and will help to focus the efforts of structural and functional genomics of this large and biotechnologically important class of enzymes.

Microarray-based global mapping of integration sites for the retrotransposon, intracisternal A-particle, in the mouse genome

2008-05-01

Mammalian genomes contain numerous evolutionary harbored mobile elements, a part of which are still active and may cause genomic instability. Their movement and positional diversity occasionally result in phenotypic changes and variation by causing altered expression or disruption of neighboring host genes. Here, we describe a novel microarray-based method by which dispersed genomic locations of a type of retrotransposon in a mammalian genome can be identified. Using this method, we mapped the DNA elements for a mouse retrotransposon, intracisternal A-particle (IAP), within genomes of C3H/He and C57BL/6J inbred mouse strains; consequently we detected hundreds of probable IAP cDNA–integrated genomic regions, in which a considerable number of strain-specific putative insertions were included. In addition, by comparing genomic DNAs from radiation-induced myeloid leukemia cells and its reference normal tissue, we detected three genomic regions around which an IAP element was integrated. These results demonstrate the first successful genome-wide mapping of a retrotransposon type in a mammalian genome.

Enhanced identification and biological validation of differential gene expression via Illumina whole-genome expression arrays through the use of the model-based background correction methodology

Ding, L.-H., Xie, Y., Park, S., Xiao, G., Story, M. D. — 2008-05-01

Despite the tremendous growth of microarray usage in scientific studies, there is a lack of standards for background correction methodologies, especially in single-color microarray platforms. Traditional background subtraction methods often generate negative signals and thus cause large amounts of data loss. Hence, some researchers prefer to avoid background corrections, which typically result in the underestimation of differential expression. Here, by utilizing nonspecific negative control features integrated into Illumina whole genome expression arrays, we have developed a method of model-based background correction for BeadArrays (MBCB). We compared the MBCB with a method adapted from the Affymetrix robust multi-array analysis algorithm and with no background subtraction, using a mouse acute myeloid leukemia (AML) dataset. We demonstrated that differential expression ratios obtained by using the MBCB had the best correlation with quantitative RT–PCR. MBCB also achieved better sensitivity in detecting differentially expressed genes with biological significance. For example, we demonstrated that the differential regulation of Tnfr2, Ikk and NF-kappaB, the death receptor pathway, in the AML samples, could only be detected by using data after MBCB implementation. We conclude that MBCB is a robust background correction method that will lead to more precise determination of gene expression and better biological interpretation of Illumina BeadArray data.

Sequence-specific DNA cleavage mediated by bipyridine polyamide conjugates

2008-05-01

The design of molecules that damage a selected DNA sequence provides a formidable opportunity for basic and applied biology. For example, such molecules offer new prospects for controlled manipulation of the genome. The conjugation of DNA-code reading molecules such as polyamides to reagents that induce DNA damages provides an approach to reach this goal. In this work, we showed that a bipyridine conjugate of polyamides was able to induce sequence-specific DNA breaks in cells. We synthesized compounds based on two polyamide parts linked to bipyridine at different positions. Bipyridine conjugates of polyamides were found to have a high affinity for the DNA target and one of them produced a specific and high-yield cleavage in vitro and in cultured cells. The bipyridine conjugate studied here, also presents cell penetrating properties since it is active when directly added to cell culture medium. Harnessing DNA damaging molecules such as bipyridine to predetermined genomic sites, as achieved here, provides an attractive strategy for targeted genome modification and DNA repair studies.

Chemical mapping of cytosines enzymatically flipped out of the DNA helix

Daujotyte, D., Liutkeviciute, Z., Tamulaitis, G., Klimasauskas, S. — 2008-05-01

Haloacetaldehydes can be employed for probing unpaired DNA structures involving cytosine and adenine residues. Using an enzyme that was structurally proven to flip its target cytosine out of the DNA helix, the HhaI DNA methyltransferase (M.HhaI), we demonstrate the suitability of the chloroacetaldehyde modification for mapping extrahelical (flipped-out) cytosine bases in protein–DNA complexes. The generality of this method was verified with two other DNA cytosine-5 methyltransferases, M.AluI and M.SssI, as well as with two restriction endonucleases, R.Ecl18kI and R.PspGI, which represent a novel class of base-flipping enzymes. Our results thus offer a simple and convenient laboratory tool for detection and mapping of flipped-out cytosines in protein–DNA complexes.

Human and mouse introns are linked to the same processes and functions through each genome's most frequent non-conserved motifs

Tsirigos, A., Rigoutsos, I. — 2008-05-01

We identified the most frequent, variable-length DNA sequence motifs in the human and mouse genomes and sub-selected those with multiple recurrences in the intergenic and intronic regions and at least one additional exonic instance in the corresponding genome. We discovered that these motifs have virtually no overlap with intronic sequences that are conserved between human and mouse, and thus are genome-specific. Moreover, we found that these motifs span a substantial fraction of previously uncharacterized human and mouse intronic space. Surprisingly, we found that these genome-specific motifs are over-represented in the introns of genes belonging to the same biological processes and molecular functions in both the human and mouse genomes even though the underlying sequences are not conserved between the two genomes. In fact, the processes and functions that are linked to these genome-specific sequence-motifs are distinct from the processes and functions which are associated with intronic regions that are conserved between human and mouse. The findings show that intronic regions from different genomes are linked to the same processes and functions in the absence of underlying sequence conservation. We highlight the ramifications of this observation with a concrete example that involves the microsatellite instability gene MLH1.

MassNet: a functional annotation service for protein mass spectrometry data

2008-04-29

Although mass spectrometry has been frequently used to identify proteins, there are no web servers that provide comprehensive functional annotation of those identified proteins. It is necessary to provide such web service due to a rapid increase in the data. We, therefore, introduce MassNet, which provides (i) physico-chemical analysis information, (ii) KEGG pathway assignment (iii) Gene Ontology mapping and (iv) protein–protein interaction (PPI) prediction for the data from MASCOT, Prospector and Profound. MassNet provides the prediction information for PPIs using both 3D structural interaction and experimental interaction deposited in PSIMAP, BIND, DIP, HPRD, IntAct, MINT, CYGD and BioGrid. The web service is freely available at http://massnet.kr or http://sequenceome.kobic.re.kr/MassNet/.

Designating eukaryotic orthology via processed transcription units

Ho, M.-R., Jang, W.-J., Chen, C.-h., Ch'ang, L.-Y., Lin, W.-c. — 2008-04-29

Orthology is a widely used concept in comparative and evolutionary genomics. In addition to prokaryotic orthology, delineating eukaryotic orthology has provided insight into the evolution of higher organisms. Indeed, many eukaryotic ortholog databases have been established for this purpose. However, unlike prokaryotes, alternative splicing (AS) has hampered eukaryotic orthology assignments. Therefore, existing databases likely contain ambiguous eukaryotic ortholog relationships and possibly misclassify alternatively spliced protein isoforms as in-paralogs, which are duplicated genes that arise following speciation. Here, we propose a new approach for designating eukaryotic orthology using processed transcription units, and we present an orthology database prototype using the human and mouse genomes. Currently existing programs cover less than 69% of the human reference sequences when assigning human/mouse orthologs. In contrast, our method encompasses up to 80% of the human reference sequences. Moreover, the ortholog database presented herein is more than 92% consistent with the existing databases. In addition to managing AS, this approach is capable of identifying orthologs of embedded genes and fusion genes using syntenic evidence. In summary, this new approach is sensitive, specific and can generate a more comprehensive and accurate compilation of eukaryotic orthologs.

MolAxis: a server for identification of channels in macromolecules

Yaffe, E., Fishelovitch, D., Wolfson, H. J., Halperin, D., Nussinov, R. — 2008-04-29

MolAxis is a freely available, easy-to-use web server for identification of channels that connect buried cavities to the outside of macromolecules and for transmembrane (TM) channels in proteins. Biological channels are essential for physiological processes such as electrolyte and metabolite transport across membranes and enzyme catalysis, and can play a role in substrate specificity. Motivated by the importance of channel identification in macromolecules, we developed the MolAxis server. MolAxis implements state-of-the-art, accurate computational-geometry techniques that reduce the dimensions of the channel finding problem, rendering the algorithm extremely efficient. Given a protein or nucleic acid structure in the PDB format, the server outputs all possible channels that connect buried cavities to the outside of the protein or points to the main channel in TM proteins. For each channel, the gating residues and the narrowest radius termed ‘bottleneck’ are also given along with a full list of the lining residues and the channel surface in a 3D graphical representation. The users can manipulate advanced parameters and direct the channel search according to their needs. MolAxis is available as a web server or as a stand-alone program at http://bioinfo3d.cs.tau.ac.il/MolAxis.

NOBAI: a web server for character coding of geometrical and statistical features in RNA structure

Knudsen, V., Caetano-Anolles, G. — 2008-04-29

The Numeration of Objects in Biology: Alignment Inferences (NOBAI) web server provides a web interface to the applications in the NOBAI software package. This software codes topological and thermodynamic information related to the secondary structure of RNA molecules as multi-state phylogenetic characters, builds character matrices directly in NEXUS format and provides sequence randomization options. The web server is an effective tool that facilitates the search for evolutionary history embedded in the structure of functional RNA molecules. The NOBAI web server is accessible at ‘http://www.manet.uiuc.edu/nobai/nobai.php’. This web site is free and open to all users and there is no login requirement.

Snap-to-it probes: chelate-constrained nucleobase oligomers with enhanced binding specificity

Morgan, J. R., Lyon, R. P., Maeda, D. Y., Zebala, J. A. — 2008-04-29

We describe snap-to-it probes, a novel probe technology to enhance the hybridization specificity of natural and unnatural nucleic acid oligomers using a simple and readily introduced structural motif. Snap-to-it probes were prepared from peptide nucleic acid (PNA) oligomers by modifying each terminus with a coordinating ligand. The two coordinating ligands constrain the probe into a macrocyclic configuration through formation of an intramolecular chelate with a divalent transition metal ion. On hybridization with a DNA target, the intramolecular chelate in the snap-to-it probe dissociates, resulting in the probe ‘snapping-to’ and binding the target nucleic acid. Thermal transition analysis of snap-to-it probes with complementary and single-mismatch DNA targets revealed that the transition between free and target-bound probe conformations was a reversible equilibrium, and the intramolecular chelate provided a thermodynamic barrier to target binding that resulted in a significant increase in mismatch discrimination. A 4–6°C increase in specificity (T_m) was observed from snap-to-it probes bearing either terminal iminodiacetic acid ligands coordinated with Ni²⁺, or terminal dihistidine and nitrilotriacetic acid ligands coordinated with Cu²⁺. The difference in specificity of the PNA oligomer relative to DNA was more than doubled in snap-to-it probes. Snap-to-it probes labeled with a fluorophore-quencher pair exhibited target-dependent fluorescence enhancement upon binding with target DNA.

Enzymatic synthesis of structure-free DNA with pseudo-complementary properties

2008-04-29

Long single-stranded DNAs and RNAs possess considerable secondary structure under conditions that support stable hybrid formation with oligonucleotides. Consequently, different oligomeric probes can hybridize to the same target with efficiencies that vary by several orders of magnitude. The ability to enzymatically generate structure-free single-stranded copies of any nucleic acid without impairing Watson–Crick base pairing to short probes would eliminate this problem and significantly improve the performance of many oligonucleotide-based applications. Synthetic nucleic acids that exhibit these properties are defined as pseudo-complementary. Previously, we described a pseudo-complementary A-T couple consisting of 2-aminoadenine (nA) and 2-thiothymine (sT) bases. The nA-sT couple is a mismatch even though nA-T and A-sT are stable base pairs. Here we show that 7-alkyl-7-deazaguanine and N⁴-alkylcytosine (where alkyl = methyl or ethyl) can be used in conjunction with nA and sT to render DNA largely structure-free and pseudo-complementary. The deoxynucleoside triphosphates (dNTPs) of these bases are incorporated into DNA by selected mesophilic and thermophilic DNA polymerases and the resulting primer extension products hybridize with good specificity and stability to oligonucleotide probes composed of the standard bases. Further optimization and characterization of the synthesis and properties of pseudo-complementary DNA should lead to an ideal target for use with oligonucleotide probes that are <25 nt in length.

Base-pair neutral homozygotes can be discriminated by calibrated high-resolution melting of small amplicons

2008-04-29

Genotyping by high-resolution melting analysis of small amplicons is homogeneous and simple. However, this approach can be limited by physical and chemical components of the system that contribute to intersample melting variation. It is challenging for this method to distinguish homozygous G::C from C::G or A::T from T::A base-pair neutral variants, which comprise ~16% of all human single nucleotide polymorphisms (SNPs). We used internal oligonucleotide calibrators and custom analysis software to improve small amplicon (42–86 bp) genotyping on the LightScanner®. Three G/C (PAH c.1155C>G, CHK2 c.1-3850G>C and candidate gene BX647987 c.261+22,290C>G) and three T/A (CPS1 c.3405-29A>T, OTC c.299-8T>A and MSH2 c.1511-9A>T) human single nucleotide variants were analyzed. Calibration improved homozygote genotyping accuracy from 91.7 to 99.7% across 1105 amplicons from 141 samples for five of the six targets. The average T_m standard deviations of these targets decreased from 0.067°C before calibration to 0.022°C after calibration. We were unable to generate a small amplicon that could discriminate the BX647987 c.261+22,290C>G (rs1869458) SNP, despite reducing standard deviations from 0.086°C to 0.032°C. Two of the sites contained symmetric nearest neighbors adjacent to the SNPs. Unexpectedly, we were able to distinguish these homozygotes by T_m even though current nearest neighbor models predict that the two homozygous alleles would be identical.

Insights into anti-termination regulation of the hut operon in Bacillus subtilis: importance of the dual RNA-binding surfaces of HutP

2008-04-29

The anti-termination protein, HutP, regulates the gene expression of the hut (histidine utilization) operon of Bacillus subtilis, by destabilizing the hut terminator RNA located upstream of the coding region encoding l-histidine degradation enzymes. On the basis of biochemical, in vivo and X-ray structural analyses, we now report that HutP uses its dual RNA-binding surfaces to access two XAG-rich regions (sites I and II) within the terminator RNA to mediate the destabilization process. In this process, HutP initiates destabilization at the 5'-end of its mRNA by binding to the first XAG-rich region (site I) and then accesses the second XAG-rich region (site II), located downstream of the stable G-C-rich segment of the terminator stem. By this action, HutP appears to disrupt the G-C-rich terminator stem, and thus prevents premature termination of transcription in the RNA segment preceding the regions encoding for the histidine degradation enzymes.

High-throughput functional annotation and data mining with the Blast2GO suite

2008-04-29

Functional genomics technologies have been widely adopted in the biological research of both model and non-model species. An efficient functional annotation of DNA or protein sequences is a major requirement for the successful application of these approaches as functional information on gene products is often the key to the interpretation of experimental results. Therefore, there is an increasing need for bioinformatics resources which are able to cope with large amount of sequence data, produce valuable annotation results and are easily accessible to laboratories where functional genomics projects are being undertaken. We present the Blast2GO suite as an integrated and biologist-oriented solution for the high-throughput and automatic functional annotation of DNA or protein sequences based on the Gene Ontology vocabulary. The most outstanding Blast2GO features are: (i) the combination of various annotation strategies and tools controlling type and intensity of annotation, (ii) the numerous graphical features such as the interactive GO-graph visualization for gene-set function profiling or descriptive charts, (iii) the general sequence management features and (iv) high-throughput capabilities. We used the Blast2GO framework to carry out a detailed analysis of annotation behaviour through homology transfer and its impact in functional genomics research. Our aim is to offer biologists useful information to take into account when addressing the task of functionally characterizing their sequence data.

A new pheromone trail-based genetic algorithm for comparative genome assembly

Zhao, F., Zhao, F., Li, T., Bryant, D. A. — 2008-04-29

Gap closing is considered one of the most challenging and time-consuming tasks in bacterial genome sequencing projects, especially with the emergence of new sequencing technologies, such as pyrosequencing, which may result in large amounts of data without the benefit of large insert libraries for contig scaffolding. We propose a novel algorithm to align contigs with more than one reference genome at a time. This approach can successfully overcome the limitations of low degrees of conserved gene order for the reference and target genomes. A pheromone trail-based genetic algorithm (PGA) was used to search globally for the optimal placement for each contig. Extensive testing on simulated and real data sets shows that PGA significantly outperforms previous methods, especially when assembling genomes that are only moderately related. An extended version of PGA can predict additional candidate connections for each contig and can thus increase the likelihood of identifying the correct arrangement of each contig. The software and test data sets can be accessed at http://sourceforge.net/projects/pga4genomics/.

Crystal structure of the 25 kDa subunit of human cleavage factor Im

Coseno, M., Martin, G., Berger, C., Gilmartin, G., Keller, W., Doublie, S. — 2008-04-29

Cleavage factor I_m is an essential component of the pre-messenger RNA 3'-end processing machinery in higher eukaryotes, participating in both the polyadenylation and cleavage steps. Cleavage factor I_m is an oligomer composed of a small 25 kDa subunit (CF I_m25) and a variable larger subunit of either 59, 68 or 72 kDa. The small subunit also interacts with RNA, poly(A) polymerase, and the nuclear poly(A)-binding protein. These protein–protein interactions are thought to be facilitated by the Nudix domain of CF I_m25, a hydrolase motif with a characteristic /β/ fold and a conserved catalytic sequence or Nudix box. We present here the crystal structures of human CF I_m25 in its free and diadenosine tetraphosphate (Ap₄A) bound forms at 1.85 and 1.80 Å, respectively. CF I_m25 crystallizes as a dimer and presents the classical Nudix fold. Results from crystallographic and biochemical experiments suggest that CF I_m25 makes use of its Nudix fold to bind but not hydrolyze ATP and Ap₄A. The complex and apo protein structures provide insight into the active oligomeric state of CF I_m and suggest a possible role of nucleotide binding in either the polyadenylation and/or cleavage steps of pre-messenger RNA 3'-end processing.

Identification of transcriptional regulatory cascades in retinoic acid-induced growth arrest of HepG2 cells

Nakanishi, M., Tomaru, Y., Miura, H., Hayashizaki, Y., Suzuki, M. — 2008-04-29

All-trans retinoic acid (ATRA) is a potent inducer of cell differentiation and growth arrest. Here, we investigated ATRA-induced regulatory cascades associated with growth arrest of the human hepatoma cell line HepG2. ATRA induced >2-fold changes in the expression of 402 genes including 55 linked to cell-cycle regulation, cell growth or apoptosis during 48 h treatment. Computational search predicted that 250 transcriptional regulatory factors (TRFs) could recognize the proximal upstream regions of any of the 55 genes. Expression of 61 TRF genes was significantly changed during ATRA incubation, providing many potential regulatory edges. We focused on six TRFs that could regulate many of the 55 genes and found a total of 160 potential edges in which the expression of each of the genes was changed later than the expression change of the corresponding regulator. RNAi knockdown of the selected TRFs caused perturbation of the respective potential targets. The genes showed an opposite regulation pattern by ATRA and specific siRNA treatments were selected as strong candidates for direct TRF targets. Finally, 36 transcriptional regulatory edges were validated by chromatin immunoprecipitation. These analyses enabled us to depict a part of the transcriptional regulatory cascades closely linked to ATRA-induced cell growth arrest.

CRISPRcompar: a website to compare clustered regularly interspaced short palindromic repeats

Grissa, I., Vergnaud, G., Pourcel, C. — 2008-04-28

Clustered regularly interspaced short palindromic repeat (CRISPR) elements are a particular family of tandem repeats present in prokaryotic genomes, in almost all archaea and in about half of bacteria, and which participate in a mechanism of acquired resistance against phages. They consist in a succession of direct repeats (DR) of 24–47 bp separated by similar sized unique sequences (spacers). In the large majority of cases, the direct repeats are highly conserved, while the number and nature of the spacers are often quite diverse, even among strains of a same species. Furthermore, the acquisition of new units (DR + spacer) was shown to happen almost exclusively on one side of the locus. Therefore, the CRISPR presents an interesting genetic marker for comparative and evolutionary analysis of closely related bacterial strains. CRISPRcompar is a web service created to assist biologists in the CRISPR typing process. Two tools facilitates the in silico investigation: CRISPRcomparison and CRISPRtionary. This website is freely accessible at http://crispr.u-psud.fr/CRISPRcompar/.

GLUE-IT and PEDEL-AA: new programmes for analyzing protein diversity in randomized libraries

Firth, A. E., Patrick, W. M. — 2008-04-28

There are many methods for introducing random mutations into nucleic acid sequences. Previously, we described a suite of programmes for estimating the completeness and diversity of randomized DNA libraries generated by a number of these protocols. Our programmes suggested some empirical guidelines for library design; however, no information was provided regarding library diversity at the protein (rather than DNA) level. We have now updated our web server, enabling analysis of translated libraries constructed by site-saturation mutagenesis and error-prone PCR (epPCR). We introduce GLUE-Including Translation (GLUE-IT), which finds the expected amino acid completeness of libraries in which up to six codons have been independently varied (according to any user-specified randomization scheme). We provide two tools for assisting with experimental design: CodonCalculator, for assessing amino acids corresponding to given randomized codons; and AA-Calculator, for finding degenerate codons that encode user-specified sets of amino acids. We also present PEDEL-AA, which calculates amino acid statistics for libraries generated by epPCR. Input includes the parent sequence, overall mutation rate, library size, indel rates and a nucleotide mutation matrix. Output includes amino acid completeness and diversity statistics, and the number and length distribution of sequences truncated by premature termination codons. The web interfaces are available at http://guinevere.otago.ac.nz/stats.html.

HOMCOS: a server to predict interacting protein pairs and interacting sites by homology modeling of complex structures

Fukuhara, N., Kawabata, T. — 2008-04-28

As protein–protein interactions are crucial in most biological processes, it is valuable to understand how and where protein pairs interact. We developed a web server HOMCOS (Homology Modeling of Complex Structure, http://biunit.naist.jp/homcos) to predict interacting protein pairs and interacting sites by homology modeling of complex structures. Our server is capable of three services. The first is modeling heterodimers from two query amino acid sequences posted by users. The server performs BLAST searches to identify homologous templates in the latest representative dataset of heterodimer structures generated from the PQS database. Structure validity is evaluated by the combination of sequence similarity and knowledge-based contact potential energy as previously described. The server generates a sequence-replaced model PDB file and a MODELLER script to build full atomic models of complex structures. The second service is modeling homodimers from one query sequence. The third service is identification of potentially interacting proteins for one query sequence. The server searches the dataset of heterodimer structures for a homologous template, outputs the candidate interacting sequences in the Uniprot database homologous for the interacting partner template proteins. These features are useful for wide range of researchers to predict putative interaction sites and interacting proteins.

The RosettaDock server for local protein-protein docking

Lyskov, S., Gray, J. J. — 2008-04-28

The RosettaDock server (http://rosettadock.graylab.jhu.edu) identifies low-energy conformations of a protein–protein interaction near a given starting configuration by optimizing rigid-body orientation and side-chain conformations. The server requires two protein structures as inputs and a starting location for the search. RosettaDock generates 1000 independent structures, and the server returns pictures, coordinate files and detailed scoring information for the 10 top-scoring models. A plot of the total energy of each of the 1000 models created shows the presence or absence of an energetic binding funnel. RosettaDock has been validated on the docking benchmark set and through the Critical Assessment of PRedicted Interactions blind prediction challenge.

CoPub: a literature-based keyword enrichment tool for microarray data analysis

2008-04-28

Medline is a rich information source, from which links between genes and keywords describing biological processes, pathways, drugs, pathologies and diseases can be extracted. We developed a publicly available tool called CoPub that uses the information in the Medline database for the biological interpretation of microarray data. CoPub allows batch input of multiple human, mouse or rat genes and produces lists of keywords from several biomedical thesauri that are significantly correlated with the set of input genes. These lists link to Medline abstracts in which the co-occurring input genes and correlated keywords are highlighted. Furthermore, CoPub can graphically visualize differentially expressed genes and over-represented keywords in a network, providing detailed insight in the relationships between genes and keywords, and revealing the most influential genes as highly connected hubs. CoPub is freely accessible at http://services.nbic.nl/cgi-bin/copub/CoPub.pl.

MassTRIX: mass translator into pathways

Suhre, K., Schmitt-Kopplin, P. — 2008-04-28

Recent technical advances in mass spectrometry (MS) have brought the field of metabolomics to a point where large numbers of metabolites from numerous prokaryotic and eukaryotic organisms can now be easily and precisely detected. The challenge today lies in the correct annotation of these metabolites on the basis of their accurate measured masses. Assignment of bulk chemical formula is generally possible, but without consideration of the biological and genomic context, concrete metabolite annotations remain difficult and uncertain. MassTRIX responds to this challenge by providing a hypothesis-driven approach to high precision MS data annotation. It presents the identified chemical compounds in their genomic context as differentially colored objects on KEGG pathway maps. Information on gene transcription or differences in the gene complement (e.g. samples from different bacterial strains) can be easily added. The user can thus interpret the metabolic state of the organism in the context of its potential and, in the case of submitted transcriptomics data, real enzymatic capacities. The MassTRIX web server is freely accessible at http://masstrix.org

The HIV-1 transcriptional activator Tat has potent nucleic acid chaperoning activities in vitro

2008-04-28

The human immunodeficiency virus type 1 (HIV-1) is a primate lentivirus that causes the acquired immunodeficiency syndrome (AIDS). In addition to the virion structural proteins and enzyme precursors, that are Gag, Env and Pol, HIV-1 encodes several regulatory proteins, notably a small nuclear transcriptional activator named Tat. The Tat protein is absolutely required for virus replication since it controls proviral DNA transcription to generate the full-length viral mRNA. Tat can also regulate mRNA capping and splicing and was recently found to interfere with the cellular mi- and siRNA machinery. Because of its extensive interplay with nucleic acids, and its basic and disordered nature we speculated that Tat had nucleic acid-chaperoning properties. This prompted us to examine in vitro the nucleic acid-chaperoning activities of Tat and Tat peptides made by chemical synthesis. Here we report that Tat has potent nucleic acid-chaperoning activities according to the standard DNA annealing, DNA and RNA strand exchange, RNA ribozyme cleavage and trans-splicing assays. The active Tat(44–61) peptide identified here corresponds to the smallest known sequence with DNA/RNA chaperoning properties.

PVS: a web server for protein sequence variability analysis tuned to facilitate conserved epitope discovery

Garcia-Boronat, M., Diez-Rivero, C. M., Reinherz, E. L., Reche, P. A. — 2008-04-27

We have developed PVS (Protein Variability Server), a web-based tool that uses several variability metrics to compute the absolute site variability in multiple protein-sequence alignments (MSAs). The variability is then assigned to a user-selected reference sequence consisting of either the first sequence in the alignment or a consensus sequence. Subsequently, PVS performs tasks that are relevant for structure-function studies, such as plotting and visualizing the variability in a relevant 3D-structure. Neatly, PVS also implements some other tasks that are thought to facilitate the design of epitope discovery-driven vaccines against pathogens where sequence variability largely contributes to immune evasion. Thus, PVS can return the conserved fragments in the MSA—as defined by a user-provided variability threshold—and locate them in a relevant 3D-structure. Furthermore, PVS can return a variability-masked sequence, which can be directly submitted to the RANKPEP server for the prediction of conserved T-cell epitopes. PVS is freely available at: http://imed.med.ucm.es/PVS/.

A role for Caf1 in mRNA deadenylation and decay in trypanosomes and human cells

Schwede, A., Ellis, L., Luther, J., Carrington, M., Stoecklin, G., Clayton, C. — 2008-04-27

The eukaryotic Ccr4/Caf1/Not complex is involved in deadenylation of mRNAs. The Caf1 and Ccr4 subunits both potentially have deadenylating enzyme activity. We investigate here the roles of Ccr4 and Caf1 in deadenylation in two organisms that separated early in eukaryotic evolution: humans and trypanosomes. In Trypanosoma brucei, we found a complex containing CAF1, NOT1, NOT2 and NOT5, DHH1 and a possible homologue of Caf130; no homologue of Ccr4 was found. Trypanosome CAF1 has deadenylation activity, and is essential for cell survival. Depletion of trypanosome CAF1 delayed deadenylation and degradation of constitutively expressed mRNAs. Human cells have two isozymes of Caf1: simultaneous depletion of both inhibited degradation of an unstable reporter mRNA. In both species, depletion of Caf1 homologues inhibited deadenylation of bulk RNA and resulted in an increase in average poly(A) tail length.

Rapid unwinding of triplet repeat hairpins by Srs2 helicase of Saccharomyces cerevisiae

Dhar, A., Lahue, R. S. — 2008-04-25

Expansions of trinucleotide repeats cause at least 15 heritable human diseases. Single-stranded triplet repeat DNA in vitro forms stable hairpins in a sequence-dependent manner that correlates with expansion risk in vivo. Hairpins are therefore considered likely intermediates during the expansion process. Unwinding of a hairpin by a DNA helicase would help protect against expansions. Yeast Srs2, but not the RecQ homolog Sgs1, blocks expansions in vivo in a manner largely dependent on its helicase function. The current study tested the idea that Srs2 would be faster at unwinding DNA substrates with an extrahelical triplet repeat hairpin embedded in a duplex context. These substrates should mimic the relevant intermediate structure thought to occur in vivo. Srs2 was faster than Sgs1 at unwinding several substrates containing triplet repeat hairpins or another structured loop. In contrast, control substrates with an unstructured loop or a Watson–Crick duplex were unwound equally well by both enzymes. Results with a fluorescently labeled, three-way junction showed that Srs2 unwinding proceeds unabated through extrahelical triplet repeats. In summary, Srs2 maintains its facile unwinding of triplet repeat hairpins embedded within duplex DNA, supporting the genetic evidence that Srs2 is a key helicase in Saccharomyces cerevisiae for preventing expansions.

Software.ncrna.org: web servers for analyses of RNA sequences

2008-04-25

We present web servers for analysis of non-coding RNA sequences on the basis of their secondary structures. Software tools for structural multiple sequence alignments, structural pairwise sequence alignments and structural motif findings are available from the integrated web server and the individual stand-alone web servers. The servers are located at http://software.ncrna.org, along with the information for the evaluation and downloading. This website is freely available to all users and there is no login requirement.

OligoHeatMap (OHM): an online tool to estimate and display hybridizations of oligonucleotides onto DNA sequences

Croce, O., Chevenet, F., Christen, R. — 2008-04-25

The efficiency of molecular methods involving DNA/DNA hybridizations depends on the accurate prediction of the melting temperature (T_m) of the duplex. Many softwares are available for T_m calculations, but difficulties arise when one wishes to check if a given oligomer (PCR primer or probe) hybridizes well or not on more than a single sequence. Moreover, the presence of mismatches within the duplex is not sufficient to estimate specificity as it does not always significantly decrease the T_m. OHM (OligoHeatMap) is an online tool able to provide estimates of T_m for a set of oligomers and a set of aligned sequences, not only as text files of complete results but also in a graphical way: T_m values are translated into colors and displayed as a heat map image, either stand alone or to be used by softwares such as TreeDyn to be included in a phylogenetic tree. OHM is freely available at http://bioinfo.unice.fr/ohm/, with links to the full source code and online help.

ProdoNet: identification and visualization of prokaryotic gene regulatory and metabolic networks

2008-04-25

ProdoNet is a web-based application for the mapping of prokaryotic genes and the corresponding proteins to common gene regulatory and metabolic networks. For a given list of genes, the system detects shared operons, identifies co-expressed genes and deduces joint regulators. In addition, the contribution to shared metabolic pathways becomes visible on KEGG maps. Furthermore, the co-occurrence of genes of interest in gene expression profiles can be added to the visualization of the global network. In this way, ProdoNet provides the basis for functional genomics approaches and for the interpretation of transcriptomics and proteomics data. As an example, we present an investigation of an experimental membrane subproteome analysis of Pseudomonas aeruginosa with ProdoNet. The ProdoNet dataset on transcriptional regulation is based on the PRODORIC Prokaryotic Database of Gene Regulation and the Virtual Footprint tool. ProdoNet is accessible at http://www.prodonet.tu-bs.de.

GenoWatch: a disease gene mining browser for association study

2008-04-25

A human gene association study often involves several genomic markers such as single nucleotide polymorphisms (SNPs) or short tandem repeat polymorphisms, and many statistically significant markers may be identified during the study. GenoWatch can efficiently extract up-to-date information about multiple markers and their associated genes in batch mode from many relevant biological databases in real-time. The comprehensive gene information retrieved includes gene ontology, function, pathway, disease, related articles in PubMed and so on. Subsequent SNP functional impact analysis and primer design of a target gene for re-sequencing can also be done in a few clicks. The presentation of results has been carefully designed to be as intuitive as possible to all users.

The GenoWatch is available at the website http://genepipe.ngc.sinica.edu.tw/genowatch