Comparative mutational analysis of cis-acting RNA signals for translational frameshifting in HIV-1 and HTLV-2

Human immunodeficiency virus type 1 (HIV-1) and human T cell leukemia virus type II (HTLV-2) use a similar mechanism for –1 translational frameshifting to overcome the termination codon in viral RNA at the end of the gag gene. Previous studies have identified two important RNA signals for frameshifting, the slippery sequence and a downstream stem–loop structure. However, there have been somewhat conflicting reports concerning the individual contributions of these sequences. In this study we have performed a comprehensive mutational analysis of the cis-acting RNA sequences involved in HIV-1 gag–pol and HTLV-2 gag–pro frameshifting. Using an in vitro translation system we determined frameshifting efficiencies for shuffled HIV-1/HTLV-2 RNA elements in a background of HIV-1 or HTLV-2 sequences. We show that the ability of the slippery sequence and stem–loop to promote ribosomal frameshifting is influenced by the flanking upstream sequence and the nucleotides in the spacer element. A wide range of frameshift efficiency rates was observed for both viruses when shuffling single sequence elements. The results for HIV-1/HTLV-2 chimeric constructs represent strong evidence supporting the notion that the viral wild-type sequences are not designed for maximal frameshifting activity but are optimized to a level suited to efficient viral replication.

RegA proteins from phage T4 and RB69 have conserved helix–loop groove RNA binding motifs but different RNA binding specificities

The RegA proteins from the bacteriophage T4 and RB69 are translational repressors that control the expression of multiple phage mRNAs. RegA proteins from the two phages share 78% sequence identity; however, in vivo expression studies have suggested that the RB69 RegA protein binds target RNAs with a higher affinity than T4 RegA protein. To study the RNA binding properties of T4 and RB69 RegA proteins more directly, the binding sites of RB69 RegA protein on synthetic RNAs corresponding to the translation initiation region of two RB69 target genes were mapped by RNase protection assays. These assays revealed that RB69 RegA protein protects nucleotides –9 to –3 (relative to the start codon) on RB69 gene 44, which contains the sequence GAAAAUU. On RB69 gene 45, the protected site (nucleotides –8 to –3) contains a similar purine-rich sequence: GAAAUA. Interestingly, T4 RegA protein protected the same nucleotides on these RNAs. To examine the specificity of RNA binding, quantitative RNA gel shift assays were performed with synthetic RNAs corresponding to recognition elements (REs) in three T4 and three RB69 mRNAs. Comparative gel shift assays demonstrated that RB69 RegA protein has an ∼7-fold higher affinity for T4 gene 44 RE RNA than T4 RegA protein. RB69 RegA protein also binds RB69 gene 44 RE RNA with a 4-fold higher affinity than T4 RegA protein. On the other hand, T4 RegA exhibited a higher affinity than RB69 RegA protein for RB69 gene 45 RE RNA. With respect to their affinities for cognate RNAs, both RegA proteins exhibited the following hierarchy of affinities: gene 44 > gene 45 > regA. Interestingly, T4 RegA exhibited the highest affinity towards RB69 gene 45 RE RNA, whereas RB69 RegA protein had the highest affinity for T4 gene 44 RE RNA. The helix–loop groove RNA binding motif of T4 RegA protein is fully conserved in RB69 RegA protein. However, homology modeling of the structure of RB69 RegA protein reveals that the divergent residues are clustered in two areas of the surface, and that there are two large areas of high conservation near the helix–loop groove, which may also play a role in RNA binding.

Identification in mice of two isoforms of the cysteinyl leukotriene 1 receptor that result from alternative splicing

Two classes of human G protein-coupled receptors, cysteinyl leukotriene 1 (CysLT1) and CysLT2 receptors, recently have been characterized and cloned. Because the CysLT1 receptor blockers are effective in treating human bronchial asthma and the mouse is often used to model human diseases, we isolated the mouse CysLT1 receptor from a mouse lung cDNA library and found two isoforms. A short isoform cDNA containing two exons encodes a polypeptide of 339 aa with 87.3% amino acid identity to the human CysLT1 receptor. A long isoform has two additional exons and an in-frame upstream start codon resulting in a 13-aa extension at the N terminus. Northern blot analysis revealed that the mouse CysLT1 receptor mRNA is expressed in lung and skin; and reverse transcription–PCR showed wide expression of the long isoform with the strongest presence in lung and skin. The gene for the mouse CysLT1 receptor was mapped to band XD. Leukotriene (LT) D4 induced intracellular calcium mobilization in Chinese hamster ovary cells stably expressing either isoform of the mouse CysLT1 receptor cDNA. This agonist effect of LTD4 was fully inhibited by the CysLT1 receptor antagonist, MK-571. Microsomal membranes from each transformant showed a single class of binding sites for [3H]LTD4; and the binding was blocked by unlabeled LTs, with the rank order of affinities being LTD4 >> LTE4 = LTC4 >> LTB4. Thus, the dominant mouse isoform with the N-terminal amino acid extension encoded by an additional exon has the same ligand response profile as the spliced form and the human receptor.

The Crithidia fasciculata RNH1 gene encodes both nuclear and mitochondrial isoforms of RNase H

The Crithidia fasciculata RNH1 gene encodes an RNase H, an enzyme that specifically degrades the RNA strand of RNA–DNA hybrids. The RNH1 gene is contained within an open reading frame (ORF) predicted to encode a protein of 53.7 kDa. Previous work has shown that RNH1 expresses two proteins: a 38 kDa protein and a 45 kDa protein which is enriched in kinetoplast extracts. Epitope tagging of the C-terminus of the RNH1 gene results in localization of the protein to both the kinetoplast and the nucleus. Translation of the ORF beginning at the second in-frame methionine codon predicts a protein of 38 kDa. Insertion of two tandem stop codons between the first ATG codon and the second in-frame ATG codon of the ORF results in expression of only the 38 kDa protein and the protein localizes specifically to the nucleus. Mutation of the second methionine codon to a valine codon prevents expression of the 38 kDa protein and results in exclusive production of the 45 kDa protein and localization of the protein only in the kinetoplast. These results suggest that the kinetoplast enzyme results from processing of the full-length 53.7 kDa protein. The nuclear enzyme appears to result from translation initiation at the second in-frame ATG codon. This is the first example in trypanosomatids of the production of nuclear and mitochondrial isoforms of a protein from a single gene and is the only eukaryotic gene in the RNase HI gene family shown to encode a mitochondrial RNase H.

Pseudoknots in prion protein mRNAs confirmed by comparative sequence analysis and pattern searching

The human prion gene contains five copies of a 24 nt repeat that is highly conserved among species. An analysis of folding free energies of the human prion mRNA, in particular in the repeat region, suggested biased codon selection and the presence of RNA patterns. In particular, pseudoknots, similar to the one predicted by Wills in the human prion mRNA, were identified in the repeat region of all available prion mRNAs available in GenBank, but not those of birds and the red slider turtle. An alignment of these mRNAs, which share low sequence homology, shows several co-variations that maintain the pseudoknot pattern. The presence of pseudoknots in yeast Sup35p and Rnq1 suggests acquisition in the prokaryotic era. Computer generated three-dimensional structures of the human prion pseudoknot highlight protein and RNA interaction domains, which suggest a possible effect in prion protein translation. The role of pseudoknots in prion diseases is discussed as individuals with extra copies of the 24 nt repeat develop the familial form of Creutzfeldt–Jakob disease.

Mapping the early steps in the pH-induced conformational conversion of the prion protein

Under certain conditions, the prion protein (PrP) undergoes a conformational change from the normal cellular isoform, PrPC, to PrPSc, an infectious isoform capable of causing neurodegenerative diseases in many mammals. Conversion can be triggered by low pH, and in vivo this appears to take place in an endocytic pathway and/or caveolae-like domains. It has thus far been impossible to characterize the conformational change at high resolution by experimental methods. Therefore, to investigate the effect of acidic pH on PrP conformation, we have performed 10-ns molecular dynamics simulations of PrPC in water at neutral and low pH. The core of the protein is well maintained at neutral pH. At low pH, however, the protein is more dynamic, and the sheet-like structure increases both by lengthening of the native β-sheet and by addition of a portion of the N terminus to widen the sheet by another two strands. The side chain of Met-129, a polymorphic codon in humans associated with variant Creutzfeldt–Jakob disease, pulls the N terminus into the sheet. Neutralization of Asp-178 at low pH removes interactions that inhibit conversion, which is consistent with the Asp-178–Asn mutation causing human prion diseases.

Selenoprotein oxidoreductase with specificity for thioredoxin and glutathione systems

Thioredoxin (Trx) and glutathione (GSH) systems are considered to be two major redox systems in animal cells. They are reduced by NADPH via Trx reductase (TR) or oxidized GSH (GSSG) reductase and further supply electrons for deoxyribonucleotide synthesis, antioxidant defense, and redox regulation of signal transduction, transcription, cell growth, and apoptosis. We cloned and characterized a pyridine nucleotide disulfide oxidoreductase, Trx and GSSG reductase (TGR), that exhibits specificity for both redox systems. This enzyme contains a selenocysteine residue encoded by the TGA codon. TGR can reduce Trx, GSSG, and a GSH-linked disulfide in in vitro assays. This unusual substrate specificity is achieved by an evolutionary conserved fusion of the TR and glutaredoxin domains. These observations, together with the biochemical probing and molecular modeling of the TGR structure, suggest a mechanism whereby the C-terminal selenotetrapeptide serves a role of a protein-linked GSSG and shuttles electrons from the disulfide center within the TR domain to either the glutaredoxin domain or Trx.

Highly expressed and alien genes of the Synechocystis genome

Comparisons of codon frequencies of genes to several gene classes are used to characterize highly expressed and alien genes on the Synechocystis PCC6803 genome. The primary gene classes include the ensemble of all genes (average gene), ribosomal protein (RP) genes, translation processing factors (TF) and genes encoding chaperone/degradation proteins (CH). A gene is predicted highly expressed (PHX) if its codon usage is close to that of the RP/TF/CH standards but strongly deviant from the average gene. Putative alien (PA) genes are those for which codon usage is significantly different from all four classes of gene standards. In Synechocystis, 380 genes were identified as PHX. The genes with the highest predicted expression levels include many that encode proteins vital for photosynthesis. Nearly all of the genes of the RP/TF/CH gene classes are PHX. The principal glycolysis enzymes, which may also function in CO2 fixation, are PHX, while none of the genes encoding TCA cycle enzymes are PHX. The PA genes are mostly of unknown function or encode transposases. Several PA genes encode polypeptides that function in lipopolysaccharide biosynthesis. Both PHX and PA genes often form significant clusters (operons). The proteins encoded by PHX and PA genes are described with respect to functional classifications, their organization in the genome and their stoichiometry in multi-subunit complexes.

Comparison of intron-containing and intron-lacking human genes elucidates putative exonic splicing enhancers

Of the rules used by the splicing machinery to precisely determine intron–exon boundaries only a fraction is known. Recent evidence suggests that specific short sequences within exons help in defining these boundaries. Such sequences are known as exonic splicing enhancers (ESE). A possible bioinformatical approach to studying ESE sequences is to compare genes that harbor introns with genes that do not. For this purpose two non-redundant samples of 719 intron-containing and 63 intron-lacking human genes were created. We performed a statistical analysis on these datasets of intron-containing and intron-lacking human coding sequences and found a statistically significant difference (P = 0.01) between these samples in terms of 5–6mer oligonucleotide distributions. The difference is not created by a few strong signals present in the majority of exons, but rather by the accumulation of multiple weak signals through small variations in codon frequencies, codon biases and context-dependent codon biases between the samples. A list of putative novel human splicing regulation sequences has been elucidated by our analysis.

Nonsense-mediated Decay of mRNA for the Selenoprotein Phospholipid Hydroperoxide Glutathione Peroxidase Is Detectable in Cultured Cells but Masked or Inhibited in Rat Tissues

Previous studies of mRNA for classical glutathione peroxidase 1 (GPx1) demonstrated that hepatocytes of rats fed a selenium-deficient diet have less cytoplasmic GPx1 mRNA than hepatocytes of rats fed a selenium-adequate diet. This is because GPx1 mRNA is degraded by the surveillance pathway called nonsense-mediated mRNA decay (NMD) when the selenocysteine codon is recognized as nonsense. Here, we examine the mechanism by which the abundance of phospholipid hydroperoxide glutathione peroxidase (PHGPx) mRNA, another selenocysteine-encoding mRNA, fails to decrease in the hepatocytes and testicular cells of rats fed a selenium-deficient diet. We demonstrate with cultured NIH3T3 fibroblasts or H35 hepatocytes transiently transfected with PHGPx gene variants under selenium-supplemented or selenium-deficient conditions that PHGPx mRNA is, in fact, a substrate for NMD when the selenocysteine codon is recognized as nonsense. We also demonstrate that the endogenous PHGPx mRNA of untransfected H35 cells is subject to NMD. The failure of previous reports to detect the NMD of PHGPx mRNA in cultured cells is likely attributable to the expression of PHGPx cDNA rather than the PHGPx gene. We conclude that 1) the sequence of the PHGPx gene is adequate to support the NMD of product mRNA, and 2) there is a mechanism in liver and testis but not cultured fibroblasts and hepatocytes that precludes or masks the NMD of PHGPx mRNA.

Predicted highly expressed and putative alien genes of Deinococcus radiodurans and implications for resistance to ionizing radiation damage

Predicted highly expressed (PHX) and putative alien genes determined by codon usages are characterized in the genome of Deinococcus radiodurans (strain R1). Deinococcus radiodurans (DEIRA) can survive very high doses of ionizing radiation that are lethal to virtually all other organisms. It has been argued that DEIRA is endowed with enhanced repair systems that provide protection and stability. However, predicted expression levels of DNA repair proteins with the exception of RecA tend to be low and do not distinguish DEIRA from other prokaryotes. In this paper, the capability of DEIRA to resist extreme doses of ionizing and UV radiation is attributed to an unusually high number of PHX chaperone/degradation, protease, and detoxification genes. Explicitly, compared with all current complete prokaryotic genomes, DEIRA contains the greatest number of PHX detoxification and protease proteins. Other sources of environmental protection against severe conditions of UV radiation, desiccation, and thermal effects for DEIRA are the several S-layer (surface structure) PHX proteins. The top PHX gene of DEIRA is the multifunctional tricarboxylic acid (TCA) gene aconitase, which, apart from its role in respiration, also alerts the cell to oxidative damage.

Long term trends in the evolution of H(3) HA1 human influenza type A

We have studied the HA1 domain of 254 human influenza A(H3N2) virus genes for clues that might help identify characteristics of hemagglutinins (HAs) of circulating strains that are predictive of that strain’s epidemic potential. Our preliminary findings include the following. (i) The most parsimonious tree found requires 1,260 substitutions of which 712 are silent and 548 are replacement substitutions. (ii) The HA1 portion of the HA gene is evolving at a rate of 5.7 nucleotide substitutions/year or 5.7 × 10−3 substitutions/site per year. (iii) The replacement substitutions are distributed randomly across the three positions of the codon when allowance is made for the number of ways each codon can change the encoded amino acid. (iv) The replacement substitutions are not distributed randomly over the branches of the tree, there being 2.2 times more changes per tip branch than for non-tip branches. This result is independent of how the virus was amplified (egg grown or kidney cell grown) prior to sequencing or if sequencing was carried out directly on the original clinical specimen by PCR. (v) These excess changes on the tip branches are probably the result of a bias in the choice of strains to sequence and the detection of deleterious mutations that had not yet been removed by negative selection. (vi) There are six hypervariable codons accumulating replacement substitutions at an average rate that is 7.2 times that of the other varied codons. (vii) The number of variable codons in the trunk branches (the winners of the competitive race against the immune system) is 47 ± 5, significantly fewer than in the twigs (90 ± 7), which in turn is significantly fewer variable codons than in tip branches (175 ± 8). (viii) A minimum of one of every 12 branches has nodes at opposite ends representing viruses that reside on different continents. This is, however, no more than would be expected if one were to randomly reassign the continent of origin of the isolates. (ix) Of 99 codons with at least four mutations, 31 have ratios of non-silent to silent changes with probabilities less than 0.05 of occurring by chance, and 14 of those have probabilities <0.005. These observations strongly support positive Darwinian selection. We suggest that the small number of variable positions along the successful trunk lineage, together with knowledge of the codons that have shown positive selection, may provide clues that permit an improved prediction of which strains will cause epidemics and therefore should be used for vaccine production.

Cellular mechanisms of β-amyloid production and secretion

The major constituent of senile plaques in Alzheimer’s disease is a 42-aa peptide, referred to as β-amyloid (Aβ). Aβ is generated from a family of differentially spliced, type-1 transmembrane domain (TM)-containing proteins, called APP, by endoproteolytic processing. The major, relatively ubiquitous pathway of APP metabolism in cell culture involves cleavage by α-secretase, which cleaves within the Aβ sequence, thus precluding Aβ formation and deposition. An alternate secretory pathway, enriched in neurons and brain, leads to cleavage of APP at the N terminus of the Aβ peptide by β-secretase, thus generating a cell-associated β-C-terminal fragment (β-CTF). A pathogenic mutation at codons 670/671 in APP (APP “Swedish”) leads to enhanced cleavage at the β-secretase scissile bond and increased Aβ formation. An inhibitor of vacuolar ATPases, bafilomycin, selectively inhibits the action of β-secretase in cell culture, suggesting a requirement for an acidic intracellular compartment for effective β-secretase cleavage of APP. β-CTF is cleaved in the TM domain by γ-secretase(s), generating both Aβ 1–40 (90%) and Aβ 1–42 (10%). Pathogenic mutations in APP at codon 717 (APP “London”) lead to an increased proportion of Aβ 1–42 being produced and secreted. Missense mutations in PS-1, localized to chromosome 14, are pathogenic in the majority of familial Alzheimer’s pedigrees. These mutations also lead to increased production of Aβ 1–42 over Aβ 1–40. Knockout of PS-1 in transgenic animals leads to significant inhibition of production of both Aβ 1–40 and Aβ 1–42 in primary cultures, indicating that PS-1 expression is important for γ-secretase cleavages. Peptide aldehyde inhibitors that block Aβ production by inhibiting γ-secretase cleavage of β-CTF have been discovered.

Evolution of RNA editing in trypanosome mitochondria

Two different RNA editing systems have been described in the kinetoplast-mitochondrion of trypanosomatid protists. The first involves the precise insertion and deletion of U residues mostly within the coding regions of maxicircle-encoded mRNAs to produce open reading frames. This editing is mediated by short overlapping complementary guide RNAs encoded in both the maxicircle and the minicircle molecules and involves a series of enzymatic cleavage-ligation steps. The second editing system is a C34 to U34 modification in the anticodon of the imported tRNATrp, thereby permitting the decoding of the UGA stop codon as tryptophan. U-insertion editing probably originated in an ancestor of the kinetoplastid lineage and appears to have evolved in some cases by the replacement of the original pan-edited cryptogene with a partially edited cDNA. The driving force for the evolutionary fixation of these retroposition events was postulated to be the stochastic loss of entire minicircle sequence classes and their encoded guide RNAs upon segregation of the single kinetoplast DNA network into daughter cells at cell division. A large plasticity in the relative abundance of minicircle sequence classes has been observed during cell culture in the laboratory. Computer simulations provide theoretical evidence for this plasticity if a random distribution and segregation model of minicircles is assumed. The possible evolutionary relationship of the C to U and U-insertion editing systems is discussed.

Molecular mechanisms of translation initiation in eukaryotes

Translation initiation is a complex process in which initiator tRNA, 40S, and 60S ribosomal subunits are assembled by eukaryotic initiation factors (eIFs) into an 80S ribosome at the initiation codon of mRNA. The cap-binding complex eIF4F and the factors eIF4A and eIF4B are required for binding of 43S complexes (comprising a 40S subunit, eIF2/GTP/Met-tRNAi and eIF3) to the 5′ end of capped mRNA but are not sufficient to promote ribosomal scanning to the initiation codon. eIF1A enhances the ability of eIF1 to dissociate aberrantly assembled complexes from mRNA, and these factors synergistically mediate 48S complex assembly at the initiation codon. Joining of 48S complexes to 60S subunits to form 80S ribosomes requires eIF5B, which has an essential ribosome-dependent GTPase activity and hydrolysis of eIF2-bound GTP induced by eIF5. Initiation on a few mRNAs is cap-independent and occurs instead by internal ribosomal entry. Encephalomyocarditis virus (EMCV) and hepatitis C virus epitomize distinct mechanisms of internal ribosomal entry site (IRES)-mediated initiation. The eIF4A and eIF4G subunits of eIF4F bind immediately upstream of the EMCV initiation codon and promote binding of 43S complexes. EMCV initiation does not involve scanning and does not require eIF1, eIF1A, and the eIF4E subunit of eIF4F. Initiation on some EMCV-like IRESs requires additional noncanonical initiation factors, which alter IRES conformation and promote binding of eIF4A/4G. Initiation on the hepatitis C virus IRES is even simpler: 43S complexes containing only eIF2 and eIF3 bind directly to the initiation codon as a result of specific interaction of the IRES and the 40S subunit.