Sounds of silence: synonymous nucleotides as a key to biological regulation and complexity

Messenger RNA is a key component of an intricate regulatory network of its own. It accommodates numerous nucleotide signals that overlap protein coding sequences and are responsible for multiple levels of regulation and generation of biological complexity. A wealth of structural and regulatory information, which mRNA carries in addition to the encoded amino acid sequence, raises the question of how these signals and overlapping codes are delineated along non-synonymous and synonymous positions in protein coding regions, especially in eukaryotes. Silent or synonymous codon positions, which do not determine amino acid sequences of the encoded proteins, define mRNA secondary structure and stability and affect the rate of translation, folding and post-translational modifications of nascent polypeptides. The RNA level selection is acting on synonymous sites in both prokaryotes and eukaryotes and is more common than previously thought. Selection pressure on the coding gene regions follows three-nucleotide periodic pattern of nucleotide base-pairing in mRNA, which is imposed by the genetic code. Synonymous positions of the coding regions have a higher level of hybridization potential relative to non-synonymous positions, and are multifunctional in their regulatory and structural roles. Recent experimental evidence and analysis of mRNA structure and interspecies conservation suggest that there is an evolutionary tradeoff between selective pressure acting at the RNA and protein levels. Here we provide a comprehensive overview of the studies that define the role of silent positions in regulating RNA structure and processing that exert downstream effects on proteins and their functions.     

Nucleotide-dependent displacement and dynamics of the α-1 helix in kinesin revealed by site-directed spin labeling EPR

In kinesin X-ray crystal structures, the N-terminal region of the α-1 helix is adjacent to the adenine ring of the bound nucleotide, while the C-terminal region of the helix is near the neck-linker (NL). Here, we monitor the displacement of the α-1 helix within a kinesin monomer bound to microtubules (MTs) in the presence or absence of nucleotides using site-directed spin labeling EPR. Kinesin was doubly spin-labeled at the α-1 and α-2 helices, and the resulting EPR spectrum showed dipolar broadening. The inter-helix distance distribution showed that 20% of the spins have a peak characteristic of 1.4–1.7 nm separation, which is similar to what is predicted from the X-ray crystal structure, albeit 80% were beyond the sensitivity limit (>2.5 nm) of the method. Upon MT binding, the fraction of kinesin exhibiting an inter-helix distance of 1.4–1.7 nm in the presence of AMPPNP (a non-hydrolysable ATP analog) and ADP was 20% and 25%, respectively. In the absence of nucleotide, this fraction increased to 40–50%. These nucleotide-induced changes in the fraction of kinesin undergoing displacement of the α-1 helix were found to be related to the fraction in which the NL undocked from the motor core. It is therefore suggested that a shift in the α-1 helix conformational equilibrium occurs upon nucleotide binding and release, and this shift controls NL docking onto the motor core.     

Insights into Eukaryotic Primer Synthesis from Structures of the p48 Subunit of Human DNA Primase

DNA replication in all organisms requires polymerases to synthesize copies of the genome. DNA polymerases are unable to function on a bare template and require a primer. Primases are crucial RNA polymerases that perform the initial de novo synthesis, generating the first 8–10 nucleotides of the primer. Although structures of archaeal and bacterial primases have provided insights into general priming mechanisms, these proteins are not well conserved with heterodimeric (p48/p58) primases in eukaryotes. Here, we present X-ray crystal structures of the catalytic engine of a eukaryotic primase, which is contained in the p48 subunit. The structures of p48 reveal that eukaryotic primases maintain the conserved catalytic prim fold domain, but with a unique subdomain not found in the archaeal and bacterial primases. Calorimetry experiments reveal that Mn^2 + but not Mg^2 + significantly enhances the binding of nucleotide to primase, which correlates with higher catalytic efficiency in vitro. The structure of p48 with bound UTP and Mn^2 + provides insights into the mechanism of nucleotide synthesis by primase. Substitution of conserved residues involved in either metal or nucleotide binding alter nucleotide binding affinities, and yeast strains containing the corresponding Pri1p substitutions are not viable. Our results reveal that two residues (S160 and H166) in direct contact with the nucleotide were previously unrecognized as critical to the human primase active site. Comparing p48 structures to those of similar polymerases in different states of action suggests changes that would be required to attain a catalytically competent conformation capable of initiating dinucleotide synthesis.     

Influenza A virus protein PB1-F2 from different strains shows distinct structural signatures

The proapoptotic influenza A virus PB1-F2 protein contributes to viral pathogenicity and is present in most human and avian influenza isolates. The structures of full-length PB1-F2 of the influenza strains Pandemic flu 2009 H1N1, 1918 Spanish flu H1N1, Bird flu H5N1 and H1N1 PR8, have been characterized by NMR and CD spectroscopy. The study was conducted using chemically synthesized full-length PB1-F2 protein and fragments thereof. The amino acid residues 30–70 of PR8 PB1-F2 were found to be responsible for amyloid formation of the protein, which could be assigned to formation of β-sheet structures, although α-helices were the only structural features detected under conditions that mimic a membranous environment. At membranous conditions, in which the proteins are found in their most structured state, significant differences become apparent between the PB1-F2 variants investigated. In contrast to Pandemic flu 2009 H1N1 and PR8 PB1-F2, which exhibit a continuous extensive C-terminal α-helix, both Spanish flu H1N1 and Bird flu H5N1 PB1-F2 contain a loop region with residues 66–71 that divides the C-terminus into two shorter helices. The observed structural differences are located to the C-terminal ends of the proteins to which most of the known functions of these proteins have been assigned. A C-terminal helix–loop–helix motif might be a structural signature for PB1-F2 of the highly pathogenic influenza viruses as observed for 1918 Spanish flu H1N1 and Bird flu H5N1 PB1-F2. This signature could indicate the pathological nature of viruses emerging in the future and thus aid in the recognition of these viruses.     

Analysis of the local dynamics of human insulin and a rapid-acting insulin analog by hydrogen/deuterium exchange mass spectrometry

Human insulin and insulin lispro (lispro), a rapid-acting insulin analog, have identical primary structures, except for the transposition of a pair of amino acids. This mutation results in alterations in their higher order structures, with lispro dissociating more easily than human insulin. In our previous study performed using hydrogen/deuterium exchange mass spectrometry (HDX/MS), differences were observed in the rates and levels of deuteration among insulin analog products, which were found to be related to their self-association stability. In this study, we carried out peptide mapping of deuterated human insulin and lispro to determine the regions responsible for these deuteration differences and to elucidate the type of structural changes that affect their HDX reactivity. We identified A3–6 and B22–24 as the 2 regions that showed distinct differences in the number of deuterium atoms incorporated between human insulin and lispro. These regions contain residues that are thought to participate in hexamerization and dimerization, respectively. We also determined that over time, the differences in deuteration levels decreased in A3–6, whereas they increased in B22–24, suggesting a difference in the dynamics between these 2 regions. This article is part of a Special Issue entitled: Mass spectrometry in structural biology.     

Unique arrangement of coding sequences for 5 S, 5.8 S, 18 S and 25 S ribosomal RNA in Saccharomyces cerevisiae as determined by R-loop and hybridization analysis.

The arrangement of the coding sequences for the 5 S, 5.8 S, 18 S and 25 S ribosomal RNA from Saccharomyces cerevisiae was analyzed in λ-yeast hybrids containing repeating units of the ribosomal DNA. After mapping of restriction sites, the positions of the coding sequences were determined by hybridization of purified rRNAs to restriction fragments, by R-loop analysis in the electron microscope, and by electrophoresis of S₁ nuclease-treated rRNA/rDNA hybrids in alkaline agarose gels. The R-loop method was improved with respect to the length calibration of RNA/DNA duplexes and to the spreading conditions resulting in fully extended 18 S and 25 S rRNA R-loops. The qualitative results are: (1) the 5 S rRNA genes, unlike those in higher eukaryotes, alternate with the genes of the precursor for the 5.8 S, 18 S and 25 S rRNA; (2) the coding sequence for 5.8 S rRNA maps, as in higher eukaryotes, between the 18 S and 25 S rRNA coding sequences. The quantitative results are: (1) the tandemly repeating rDNA units have a constant length of 9060 ± 100 nucleotide pairs with one SstI, two HindIII and, dependent on the strain, six or seven EcoRI sites; (2) the 18 S and 25 S rRNA coding regions consist of 1710 ± 80 and 3360 ± 80 nucleotide pairs, respectively; (3) an 18 S rRNA coding region is separated by a 780 ± 70 nucleotide pairs transcribed spacer from a 25 S rRNA coding region. This is then followed by a 3210 ± 100 nucleotide pairs mainly non-transcribed spacer which contains a 5 S rRNA gene.     

Hfq assists small RNAs in binding to the coding sequence of ompD mRNA and in rearranging its structure

The bacterial protein Hfq participates in the regulation of translation by small noncoding RNAs (sRNAs). Several mechanisms have been proposed to explain the role of Hfq in the regulation by sRNAs binding to the 5′-untranslated mRNA regions. However, it remains unknown how Hfq affects those sRNAs that target the coding sequence. Here, the contribution of Hfq to the annealing of three sRNAs, RybB, SdsR, and MicC, to the coding sequence of Salmonella ompD mRNA was investigated. Hfq bound to ompD mRNA with tight, subnanomolar affinity. Moreover, Hfq strongly accelerated the rates of annealing of RybB and MicC sRNAs to this mRNA, and it also had a small effect on the annealing of SdsR. The experiments using truncated RNAs revealed that the contributions of Hfq to the annealing of each sRNA were individually adjusted depending on the structures of interacting RNAs. In agreement with that, the mRNA structure probing revealed different structural contexts of each sRNA binding site. Additionally, the annealing of RybB and MicC sRNAs induced specific conformational changes in ompD mRNA consistent with local unfolding of mRNA secondary structure. Finally, the mutation analysis showed that the long AU-rich sequence in the 5′-untranslated mRNA region served as an Hfq binding site essential for the annealing of sRNAs to the coding sequence. Overall, the data showed that the functional specificity of Hfq in the annealing of each sRNA to the ompD mRNA coding sequence was determined by the sequence and structure of the interacting RNAs.     

An intron in a ribosomal protein gene from Tetrahymena

We have cloned and sequenced a single copy gene encoding a ribosomal protein from the ciliate Tetrahymena thermophila. The gene product was identified as ribosomal protein S25 by comparison of the migration in two-dimensional polyacrylamide gels of the protein synthesized by translation in vitro of hybrid-selected mRNA and authentic ribosomal proteins. The proteins show strong homology to ribosomal protein S12 from Escherichia coli. The coding region of the gene is interrupted by a 979-bp intron 68 bp downstream of the translation start. This is the first intron in a protein encoding gene of a ciliate to be described at the nucleotide sequence level. The intron obeys the GT/AG rule for splice junctions of nuclear mRNA introns from higher eukaryotes but lacks the pyrimidine stretch usually found in the immediate vicinity of the 3' splice junction. The structure of the intron and the fact that it is found together with the well described self-splicing rRNA intron is discussed in relation to the evolution of RNA splicing.     

Secondary Structures for 5′ Regions of R2 Retrotransposon RNAs Reveal a Novel Conserved Pseudoknot and Regions that Evolve under Different Constraints

Sequences from the 5′ region of R2 retrotransposons of four species of silk moth are reported. In Bombyx mori, this region of the R2 messenger RNA contains a binding site for R2 protein and mediates interactions critical to R2 element insertion into the host genome. A model of secondary structure for a segment of this RNA is proposed on the basis of binding to oligonucleotide microarrays, chemical mapping, and comparative sequence analysis. Five conserved secondary structures are identified, including a novel pseudoknot. There is an apparent transition from an entirely RNA structure coding function in most of the 5′ segment to a protein coding function near the 3′ end. This suggests that local regions evolved under separate functional constraints (structural, coding, or both).     

RNA structure-based ribosome recruitment: lessons from the Dicistroviridae intergenic region IRESes

In eukaryotes, the canonical process of initiating protein synthesis on an mRNA depends on many large protein factors and the modified nucleotide cap on the 5' end of the mRNA. However, certain RNA sequences can bypass the need for these proteins and cap, using an RNA structure-based mechanism called internal initiation of translation. These RNAs are called internal ribosome entry sites (IRESes), and the cap-independent initiation pathway they support is critical for successful infection by many viruses of medical and economic importance. In this review, we briefly describe and compare mechanistic and structural groups of viral IRES RNAs, focusing on those IRESes that are capable of direct ribosome recruitment using specific RNA structures. We then discuss in greater detail some recent advances in our understanding of the intergenic region IRESes of the Dicistroviridae, which use the most streamlined ribosome-recruitment mechanism yet discovered. By combining these findings with knowledge of canonical translation and the behavior of other IRESes, mechanistic models of this RNA structure-based process are emerging.     

Dissection of Conformational Conversion Events during Prion Amyloid Fibril Formation Using Hydrogen Exchange and Mass Spectrometry

A molecular understanding of prion diseases requires an understanding of the mechanism of amyloid fibril formation by the prion protein. In particular, it is necessary to define the sequence of the structural events describing the conformational conversion of monomeric PrP to aggregated PrP. In this study, the sequence of the structural events in the case of amyloid fibril formation by recombinant mouse prion protein at pH 7 has been characterized by hydrogen–deuterium exchange and mass spectrometry. The observation that fibrils are substantially more stable to hydrogen–deuterium exchange than is native monomer allows both forms to be quantified during the course of the aggregation reaction. Under the aggregation conditions utilized, native monomeric protein and amyloid fibrils are the only forms of the protein detectable during the course of the fibril formation reaction, suggesting that monomer directly adds on to the fibril template. Conformational conversion is shown to occur in two steps after the binding of monomer to fibril, with helix 1 unfolding only after helices 2 and 3 transform into β-sheet. Local stability in the β-sheet core region (residues ~ 159–225) of the fibrils is shown to be sequence dependent in that it varies along the length of the core, and local stability in protein molecules that are ordered in the structurally heterogeneous sequence segment 109–132 is shown to be similar to that in the core. This new understanding of the structural events during prion protein aggregation has important bearing on our comprehension of the molecular basis of prion pathogenesis.     

Influence of the three nucleotides upstream of the initiation codon on expression of the Escherichia coli lacZ gene in Saccharomyces cerevisiae.

By introducing synthetic oligonucleotides into a lacZ-yeast expression vector a set of 47 plasmids (out of 64 possible) was generated, differing only in the three bases immediately upstream of the AUG initiation codon of the Escherichia coli lacZ gene. Expression of the beta-galactosidase fusion protein encoded by the different plasmids was determined in Saccharomyces cerevisiae by immunogel electrophoresis. Among the clones tested we found a factor 3 difference in expression. A slight nucleotide preference was found in positions -3(A > G > C = U) and -2 (G > C = U > A). The choice of the nucleotide at position -1 immediately 5' of the AUG did not effect translation efficiency. Increasing homology to the yeast consensus sequence (AAAAAAAUGUCU) was not concomitant with an increased translation efficiency. Our results indicate that the choice of nucleotides immediately preceding the initiation codon in yeast does not dramatically influence translation efficiency, as in prokaryotes or higher eukaryotes.     

Mechanistic Role of Structurally Dynamic Regions in Dicistroviridae IGR IRESs

Dicistroviridae intergenic region (IGR) internal ribosome entry site(s) (IRES) RNAs drive a cap-independent pathway of translation initiation, recruiting both small and large ribosomal subunits to viral RNA without the use of any canonical translation initiation factors. This ability is conferred by the folded three-dimensional structure of the IRES RNA, which has been solved by X-ray crystallography. Here, we report the chemical probing of Plautia stali intestine virus IGR IRES in the unbound form, in the 40S-subunit-bound form, and in the 80S-ribosome-bound form. The results, when combined with an analysis of crystal structures, suggest that parts of the IRES RNA change structure as the preinitiation complex forms. Using mutagenesis coupled with native gel electrophoresis, preinitiation complex assembly assays, and translation initiation assays, we show that these potentially structurally dynamic elements of the IRES are involved in different steps in the pathway of ribosome recruitment and translation initiation. Like tRNAs, it appears that the IGR IRES undergoes local structural changes that are coordinated with structural changes in the ribosome, and these are critical for the IRES mechanism of action.     

Conserved Element of the Dicistrovirus IGR IRES that Mimics an E-site tRNA/Ribosome Interaction Mediates Multiple Functions

The internal ribosome entry site within the intergenic region (IGR IRES) of the Dicistroviridae family mimics a tRNA to directly assemble 80 S ribosomes and initiate translation at a non-AUG codon from the ribosomal A-site. A comparison of IGR IRESs within this viral family reveals structural similarity but little sequence similarity. However, a few specific conserved elements exist, which likely have important roles in IRES function. In this study, we have generated a battery of mutations to characterize the role of a conserved loop (L1.1) region of the IGR IRES. Mutating specific nucleotides within the L1.1 region inhibited IGR IRES-mediated translation in rabbit reticulocyte lysates. By assaying different steps in IRES function, we found that the mutant L1.1 IRESs had reduced affinity for 80 S ribosomes but not 40 S subunits, indicating that the L1.1 region mediated either binding to preformed 80 S or 60 S joining. Furthermore, mutations in L1.1 altered the position of the ribosome on the mutant IRES, indicating that the tRNA-like anticodon/codon mimic within the ribosomal P-site is disrupted. Structural studies have revealed that the L1.1 region interacts with the L1 stalk of the 60 S subunit, which is similar to the interactions between the T-loop of the E-site tRNA and ribosomal protein rpL1. Our results demonstrate that the conserved L1.1 region directs multiple steps in IGR IRES-mediated translation including ribosome binding and positioning, which are functions that the E-site tRNA may normally mediate during translation.     

The translation of recombinant proteins in E. coli can be improved by in silico generating and screening random libraries of a −70/+96 mRNA region with respect to the translation initiation codon

Recombinant protein translation in Escherichia coli may be limited by stable (i.e. low free energy) secondary structures in the mRNA translation initiation region. To circumvent this issue, we have set-up a computer tool called ‘ExEnSo’ (Expression Enhancer Software) that generates a random library of 8192 sequences, calculates the free energy of secondary structures of each sequence in the −70/+96 region (base 1 is the translation initiation codon), and then selects the sequence having the highest free energy. The software uses this ‘optimized’ sequence to create a 5′ primer that can be used in PCR experiments to amplify the coding sequence of interest prior to sub-cloning into a prokaryotic expression vector. In this article, we report how ExEnSo was set-up and the results obtained with nine coding sequences with low expression levels in E. coli. The free energy of the −70/+96 region of all these coding sequences was increased compared to the non-optimized sequences. Moreover, the protein expression of eight out of nine of these coding sequences was increased in E. coli, indicating a good correlation between in silico and in vivo results. ExEnSo is available as a free online tool.