New bioinformatic tools for analysis of nucleotide modifications in eukaryotic rRNA

This report presents a valuable new bioinformatics package for research on rRNA nucleotide modifications in the ribosome, especially those created by small nucleolar RNA:protein complexes (snoRNPs). The interactive service, which is not available elsewhere, enables a user to visualize the positions of pseudouridines, 2'-O-methylations, and base methylations in three-dimensional space in the ribosome and also in linear and secondary structure formats of ribosomal RNA. Our tools provide additional perspective on where the modifications occur relative to functional regions within the rRNA and relative to other nearby modifications. This package of new tools is presented as a major enhancement of an existing but significantly upgraded yeast snoRNA database available publicly at http://people.biochem.umass.edu/sfournier/fournierlab/snornadb/. The other key features of the enhanced database include details of the base pairing of snoRNAs with target RNAs, genomic organization of the yeast snoRNA genes, and information on corresponding snoRNAs and modifications in other model organisms.     

Data mining of imperfect double-stranded RNA in 3' untranslated regions of eukaryotic mRNAs

Recent developments in the study of RNA silencing indicate that double-stranded RNA (dsRNA) can be used in eukaryotes to block expression of a corresponding cellular gene. There is also a large class of small non-coding RNAs having potential to form a distinct, stable stem-loop in numbers of eukaryotic genomes. We had reported that a large imperfect dsRNA structure with hundreds of base-pairs (bp) in the 3' untranslated region (3' UTR) of cytotoxic ribonuclease was correlated with the translation suppression. In this study, we search for such dsRNAs in a 3' UTR database. The occurrence rate of large dsRNA in 3' UTRs ranges from 0.01% in plant to 0.30% in vertebrate mRNAs. However, small imperfect dsRNAs of ~ 30 bp are much more prevalent than large ones. The small dsRNAs are statistically very significant and uniquely well-ordered. Most of them have the conserved structural features of pre-miRNAs. Our data mining of the dsRNAs in the 3' UTR database can be used to explore RNA-based regulation of gene expression.     

Sequence analysis and compositional properties of untranslated regions of human mRNAs

A detailed computer analysis of the untranslated regions, 5′-UTR and 3′- UTR, of human mRNA sequences is reported. The compositional properties of these regions, compared with those of the corresponding coding regions, indicate that 5′-UTR and 3′-UTR are less affected by the isochore compartmentalization than the corresponding third codon positions of mRNAs. The presence of higher functional constraints in 5′-UTR is also reported. Dinucleotide analysis shows a depletion of CpG and TpA in both sequences. A search for significant sequence motifs using the WORDUP algorithm reveals the patterns already known to have a functional role in the mRNA UTR, and several other motifs whose functional roles remain to be demonstrated. This type of analysis may be particularly useful for guiding site-directed mutagenesis experiments. In addition, it can be used for assessing the nature of anonymous sequences now produced in large amounts in megabase sequencing projects.     

UTRdb: a specialized database of 5' and 3' untranslated regions of eukaryotic mRNAs.

The 5' and 3' untranslated regions of eukaryotic mRNAs may play a crucial role in the regulation of gene expression controlling mRNA localization, stability and translational efficiency. For this reason we developed UTRdb (http://bigarea.area.ba.cnr.it:8000/BioWWW/#U TRdb), a specialized database of 5' and 3' untranslated sequences of eukaryotic mRNAs cleaned from redundancy. UTRdb entries are enriched with specialized information not present in the primary databases including the presence of nucleotide sequence patterns already demonstrated by experimental analysis to have some functional role. All these patterns have been collected in the UTRsite database so that it is possible to search any input sequence for the presence of annotated functional motifs. Furthermore, UTRdb entries have been annotated for the presence of repetitive elements.     

Statistical evidence for conserved, local secondary structure in the coding regions of eukaryotic mRNAs and pre-mRNAs

Owing to the degeneracy of the genetic code, protein-coding regions of mRNA sequences can harbour more than only amino acid information. We search the mRNA sequences of 11 human protein-coding genes for evolutionarily conserved secondary structure elements using RNA-Decoder, a comparative secondary structure prediction program that is capable of explicitly taking the known protein-coding context of the mRNA sequences into account. We detect well-defined, conserved RNA secondary structure elements in the coding regions of the mRNA sequences and show that base-paired codons strongly correlate with sparse codons. We also investigate the role of repetitive elements in the formation of secondary structure and explain the use of alternate start codons in the caveolin-1 gene by a conserved secondary structure element overlapping the nominal start codon. We discuss the functional roles of our novel findings in regulating the gene expression on mRNA level. We also investigate the role of secondary structure on the correct splicing of the human CFTR gene. We study the wild-type version of the pre-mRNA as well as 29 variants with synonymous mutations in exon 12. By comparing our predicted secondary structures to the experimentally determined splicing efficiencies, we find with weak statistical significance that pre-mRNAs with high-splicing efficiencies have different predicted secondary structures than pre-mRNAs with low-splicing efficiencies.     

UTRdb and UTRsite: specialized databases of sequences and functional elements of 5' and 3' untranslated regions of eukaryotic mRNAs

The 5' and 3' untranslated regions of eukaryotic mRNAs may play a crucial role in the regulation of gene expression controlling mRNA localization, stability and translational efficiency. For this reason we developed UTRdb, a specialized database of 5' and 3' untranslated sequences of eukaryotic mRNAs cleaned from redundancy. UTRdb entries are enriched with specialized information not present in the primary databases including the presence of nucleotide sequence patterns already demonstrated by experimental analysis to have some functional role. All these patterns have been collected in the UTRsite database so that it is possible to search any input sequence for the presence of annotated functional motifs. Furthermore, UTRdb entries have been annotated for the presence of repetitive elements. All internet resources implemented for retrieval and functional analysis of 5' and 3' untranslated regions of eukaryotic mRNAs are accessible at http://bigarea.area.ba.cnr.it:8000/EmbIT/UTRH ome/     

The proteome: structure, function and evolution

This paper reports two studies to model the inter-relationships between protein sequence, structure and function. First, an automated pipeline to provide a structural annotation of proteomes in the major genomes is described. The results are stored in a database at Imperial College, London (3D-GENOMICS) that can be accessed at www.sbg.bio.ic.ac.uk. Analysis of the assignments to structural superfamilies provides evolutionary insights. 3D-GENOMICS is being integrated with related proteome annotation data at University College London and the European Bioinformatics Institute in a project known as e-protein (http://www.e-protein.org/). The second topic is motivated by the developments in structural genomics projects in which the structure of a protein is determined prior to knowledge of its function. We have developed a new approach PHUNCTIONER that uses the gene ontology (GO) classification to supervise the extraction of the sequence signal responsible for protein function from a structure-based sequence alignment. Using GO we can obtain profiles for a range of specificities described in the ontology. In the region of low sequence similarity (around 15%), our method is more accurate than assignment from the closest structural homologue. The method is also able to identify the specific residues associated with the function of the protein family.     

Annealing function of GroEL: structural and bioinformatic analysis

The Escherichia coli chaperonin system, GroEL–GroES, facilitates folding of substrate proteins (SPs) that are otherwise destined to aggregate. The iterative annealing mechanism suggests that the allostery-driven GroEL transitions leading to changes in the microenvironment of the SP constitutes the annealing action of chaperonins. To describe the molecular basis for the changes in the nature of SP–GroEL interactions we use the crystal structures of GroEL (T state), GroEL–ATP (R state) and the GroEL–GroES–(ADP)₇ (R″ state) complex to determine the residue-specific changes in the accessible surface area and the number of tertiary contacts as a result of the T→R→R″ transitions. We find large changes in the accessible area in many residues in the apical domain, but relatively smaller changes are associated with residues in the equatorial domain. In the course of the T→R transition the microenvironment of the SP changes which suggests that GroEL is an annealing machine even without GroES. This is reflected in the exposure of Glu386 which loses six contacts in the T→R transition. We also evaluate the conservation of residues that participate in the various chaperonin functions. Multiple sequence alignments and chemical sequence entropy calculations reveal that, to a large extent, only the chemical identities and not the residues themselves important for the nominal functions (peptide binding, nucleotide binding, GroES and substrate protein release) are strongly conserved. Using chemical sequence entropy, which is computed by classifying aminoacids into four types (hydrophobic, polar, positively charged and negatively charged) we make several new predictions that are relevant for peptide binding and annealing function of GroEL. We identify a number of conserved peptide binding sites in the apical domain which coincide with those found in the 1.7 Å crystal structure of ‘mini-chaperone’ complexed with the N-terminal tag. Correlated mutations in the HSP60 family, that might control allostery in GroEL, are also strongly conserved. Most importantly, we find that charged solvent-exposed residues in the T state (Lys 226, Glu 252 and Asp 253) are strongly conserved. This leads to the prediction that mutating these residues, that control the annealing function of the SP, can decrease the efficacy of the chaperonin function.     

The structure and function of eukaryotic photosystem I

Eukaryotic photosystem I consists of two functional moieties: the photosystem I core, harboring the components for the light-driven charge separation and the subsequent electron transfer, and the peripheral light-harvesting complex (LHCI). While the photosystem I-core remained highly conserved throughout the evolution, with the exception of the oxidizing side of photosystem I, the LHCI complex shows a high degree of variability in size, subunits composition and bound pigments, which is due to the large variety of different habitats photosynthetic organisms dwell in. Besides summarizing the most current knowledge on the photosystem I-core structure, we will discuss the composition and structure of the LHCI complex from different eukaryotic organisms, both from the red and the green clade. Furthermore, mechanistic insights into electron transfer between the donor and acceptor side of photosystem I and its soluble electron transfer carrier proteins will be given. This article is part of a Special Issue entitled: Regulation of Electron Transport in Chloroplasts.     

Characterization of the 5' and 3' untranslated regions in murine mdm2 mRNAs

The murine double minute 2 (mdm2) gene is essential for embryogenesis in mice that express the p53 tumor suppressor protein. Mdm2 levels must be regulated tightly because overexpression of mdm2 contributes to tumorigenesis. We investigated whether the 5' and 3' untranslated regions (UTRs) of murine mdm2 affect the expression of MDM2 proteins. Induction of mdm2 expression by p53 results in synthesis of an mdm2 mRNA with a short 5' UTR. The long 5' UTR increases internal initiation of translation of a minor MDM2 protein, p76(MDM2), without affecting the efficiency of translation of the full-length p90(MDM2). We discovered two alternative 3' untranslated regions in murine mdm2 mRNA expressed in the testis. The longer 3' UTR contains a consensus instability element, but mdm2 mRNAs containing the long and short 3' UTRs have comparable half-lives. The 3' UTRs do not affect either initiation codon use or translation efficiency. Thus, the murine 5' UTR, but not the 3'UTR, influences the ratio of the two MDM2 proteins but neither UTR affects MDM2 abundance significantly.     

Serpins: structure,function and molecular evolution

The superfamily of serine proteinase inhibitors (serpins) are involved in a number of fundamental biological processes such as blood coagulation, complement activation, fibrinolysis, angiogenesis, inflammation and tumor suppression and are expressed in a cell-specific manner. The average protein size of a serpin family member is 350-400 amino acids, but gene structure varies in terms of number and size of exons and introns. Previous studies of all known serpins identified 16 clades and 10 orphan sequences. Vertebrate serpins can be conveniently classified into six sub-groups. We provide additional data that updates the phylogenetic analysis in the context of structural and functional properties of the proteins. From these, we can conclude that the functional classification of serpins relies on their protein structure and not on sequence similarity.     

Fe-only hydrogenases: structure,function and evolution

Hydrogenases are enzymes capable of catalyzing the oxidation of molecular hydrogen or its production from protons and electrons according to the reversible reaction: H(2)<==>2H(+)+2e(-). Most of these enzymes fall into to major classes: NiFe and Fe-only hydrogenases. Extensive spectroscopic, electrochemical and structural studies have shed appreciable light on the catalytic mechanism of hydrogenases. Although evolutionarily unrelated, NiFe and Fe-hydrogenases share a common, unusual feature: an active site low-spin Fe center with CO and CN coordination. We have recently focused our attention on Fe-hydrogenases because from structural studies by us and others, it appears to be a simpler system than the NiFe counterpart. Thus the primary hydrogen binding site has been identified and plausible, electron, proton and hydrogen pathways from and to the buried active site may be proposed from the structural data. The extensive genome sequencing effort currently under way has shown that eukaryotic organisms contain putatively gene coding sequences that display significant homology to Fe-hydrogenases. Here, we summarize the available evidence concerning the mechanism of these enzymes and carry out a structural comparison between Fe-hydrogenases and related proteins of unknown metal content from yeast, plant, worm, insect and mammals.     

Crocodile transthyretin: structure,function, and evolution

Structure and function were studied for Crocodylus porosus transthyretin (crocTTR), an important intermediate in TTR evolution. The cDNA for crocTTR mRNA was cloned and sequenced and the amino acid sequence of crocTTR was deduced. In contrast to mammalian TTRs, but similar to avian and lizard TTRs, the subunit of crocTTR had a long and hydrophobic NH(2)-terminal region. Different from the situation in mammals, triiodothyronine (T(3)) was bound by crocTTR with higher affinity than thyroxine (T(4)). Recombinant crocTTR and a chimeric construct, with the NH(2)-terminal region of crocTTR being replaced by that of Xenopus laevis TTR, were synthesized in the yeast Pichia pastoris. Analysis of the affinity of the chimeric TTRs showed that the NH(2)-terminal region modulates T(4) and T(3) binding characteristics of TTR. The structural differences of the NH(2)-terminal regions of reptilian and amphibian TTRs were caused by a shift in splice sites at the 5' end of exon 2. The comparison of crocodile and other vertebrate TTRs shows that TTR evolution is an example for positive Darwinian evolution and identifies its molecular mechanism.     

Protein kinase structure and function analysis with chemical tools

Protein kinases are the largest enzyme superfamily involved in cell signal transduction and represent therapeutic targets for a range of diseases. There have been intensive efforts from many labs to understand their catalytic mechanisms, discover inhibitors and discern their cellular functions. In this review, we will describe two approaches developed to analyze protein kinases: bisubstrate analog inhibition and phosphonate analog utilization. Both of these methods have been used in combination with the protein semisynthesis method expressed protein ligation to advance our understanding of kinase-substrate interactions and functional elucidation of phosphorylation. Previous work on the nature of the protein kinase mechanism suggests it follows a dissociative transition state. A bisubstrate analog was designed against the insulin receptor kinase to mimic the geometry of a dissociative transition state reaction coordinate distance. This bisubstrate compound proved to be a potent inhibitor against the insulin receptor kinase and occupied both peptide and nucleotide binding sites. Bisubstrate compounds with altered hydrogen bonding potential as well as varying spacers between the adenine and the peptide demonstrate the importance of the original design features. We have also shown that related bisubstrate analogs can be used to potently block serine/threonine kinases including protein kinase A. Since many protein kinases recognize folded protein substrates for efficient phosphorylation, it was advantageous to incorporate the peptide-ATP conjugates into protein structures. Using expressed protein ligation, a Src-ATP conjugate was produced and shown to be a high affinity ligand for the Csk tyrosine kinase. Nonhydrolyzable mimics of phosphoSer/phosphoTyr can be useful in examining the functionality of phosphorylation events. Using expressed protein ligation, we have employed phosphonomethylene phenylalanine and phosphonomethylene alanine to probe the phosphorylation of Tyr and Ser, respectively. These tools have permitted an analysis of the SH2-phosphatases (SHP1 and SHP2), revealing a novel intramolecular stimulation of catalytic activity mediated by the corresponding phosphorylation events. They have also been used to characterize the cellular regulation of the melatonin rhythm enzyme by phosphorylation.     

PFIT and PFRIT: bioinformatic algorithms for detecting glycosidase function from structure and sequence

The identification of the enzymes involved in the metabolism of simple and complex carbohydrates presents one bioinformatic challenge in the post-genomic era. Here, we present the PFIT and PFRIT algorithms for identifying those proteins adopting the alpha/beta barrel fold that function as glycosidases. These algorithms are based on the observation that proteins adopting the alpha/beta barrel fold share positions in their tertiary structures having equivalent sets of atomic interactions. These are conserved tertiary interaction positions, which have been implicated in both structure and function. Glycosidases adopting the alpha/beta barrel fold share more conserved tertiary interactions than alpha/beta barrel proteins having other functions. The enrichment pattern of conserved tertiary interactions in the glycosidases is the information that PFIT and PFRIT use to predict whether any given alpha/beta barrel will function as a glycosidase or not. Using as a test set a database of 19 glycosidase and 45 nonglycosidase alpha/beta barrel proteins with low sequence similarity, PFIT and PFRIT can correctly predict glycosidase function for 84% of the proteins known to function as glycosidases. PFIT and PFRIT incorrectly predict glycosidase function for 25% of the nonglycosidases. The program PSI-BLAST can also correctly identify 84% of the 19 glycosidases, however, it incorrectly predicts glycosidase function for 50% of the nonglycosidases (twofold greater than PFIT and PFRIT). Overall, we demonstrate that the structure-based PFIT and PFRIT algorithms are both more selective and sensitive for predicting glycosidase function than the sequence-based PSI-BLAST algorithm.