A thermodynamic analysis of the sequence-specific binding of RNA by bacteriophage MS2 coat protein

A reaction-diffusion model to study RNA motion by quantitative fluorescence recovery after photobleaching

Fluorescence recovery after photobleaching (FRAP) is a powerful technique to study molecular dynamics inside living cells. During the past years, several laboratories have used FRAP to image the motion of RNA-protein and other macromolecular complexes in the nucleus and cytoplasm. In the case of mRNAs, there is growing evidence indicating that these molecules assemble into large ribonucleoprotein complexes that diffuse throughout the nucleus by Brownian motion. However, estimates of the corresponding diffusion rate yielded values that differ by up to one order of magnitude. In vivo labeling of RNA relies on indirect tagging with a fluorescent probe, and here we show how the binding affinity of the probe to the target RNA influences the effective diffusion estimates of the resulting complex. We extend current reaction-diffusion models for FRAP by allowing for diffusion of the bound complex. This more general model can be used to fit any fluorescence recovery curve involving two interacting mobile species in the cell (a fluorescent probe and its target substrate). The results show that interpreting FRAP data in light of the new model reconciles the discrepant mRNA diffusion-rate values previously reported.

Molecular design and delivery of siRNA

Double stranded short interfering RNAs (siRNAs) mediate gene silencing in a sequence specific manner. By virtue of their specific gene silencing activity and owing to the recent discoveries on their plasmid and virus driven expression, siRNAs are being widely adopted in research and therapeutics. Efforts were made to optimize the siRNA expression system for the application in therapy. One major obstacle in developing RNA interference (RNAi) therapy is the delivery of siRNAs to the target cells. Combination of novel molecular targeting technologies, such as recombinant protein technology and ribosome display technology, will enable to deliver gene silencing agents to target cells specifically and efficiently.

Eukaryotic Transcription: What Does It Mean for a Gene to Be ‘on’?

Until recently, transcription could only be observed by measuring mRNA production of cell populations, thus obscuring the kinetics at the level of individual transcription events. A new study now shows that eukaryotic transcription, visualised in individual living cells, occurs in bursts -- much as it does in prokaryotes.

Real-time kinetics of gene activity in individual bacteria

Protein levels have been shown to vary substantially between individual cells in clonal populations. In prokaryotes, the contribution to such fluctuations from the inherent randomness of gene expression has largely been attributed to having just a few transcripts of the corresponding mRNAs. By contrast, eukaryotic studies tend to emphasize chromatin remodeling and burst-like transcription. Here, we study single-cell transcription in Escherichia coli by measuring mRNA levels in individual living cells. The results directly demonstrate transcriptional bursting, similar to that indirectly inferred for eukaryotes. We also measure mRNA partitioning at cell division and correlate mRNA and protein levels in single cells. Partitioning is approximately binomial, and mRNA-protein correlations are weaker earlier in the cell cycle, where cell division has recently randomized the relative concentrations. Our methods further extend protein-based approaches by counting the integer-valued number of transcript with single-molecule resolution. This greatly facilitates kinetic interpretations in terms of the integer-valued random processes that produce the fluctuations.

Alanine scanning of MS2 coat protein reveals protein-phosphate contacts involved in thermodynamic hot spots

The co-crystal structure of the MS2 coat protein dimer with its RNA operator reveals eight amino acid side-chains contacting seven of the RNA phosphates. These eight amino acids and five nearby control positions were individually changed to an alanine residue and the binding affinities of the mutant proteins to the RNA were determined. In general, the data agreed well with the crystal structure and previous RNA modification data. Interestingly, amino acid residues that are energetically most important for complex formation cluster in the middle of the RNA binding interface, forming thermodynamic hot spots, and are surrounded by energetically less relevant amino acids. In order to evaluate whether or not a given alanine mutation causes a global change in the RNA-protein interface, the affinities of the mutant proteins to RNAs containing one of 14 backbone modifications spanning the entire interface were determined. In three of six protein mutations tested, thermodynamic coupling between the site of the mutation and RNA groups that can be even more than 16 A away was detected. This suggests that, in some cases, the mutation may subtly alter the entire protein-RNA interface.

Functional Relevance of Amino Acid Residues Involved in Interactions with Ordered Nucleic Acid in a Spherical Virus

In the spherical virion of the parvovirus minute virus of mice, several amino acid side chains of the capsid were previously found to be involved in interactions with the viral single-stranded DNA molecule. We have individually truncated by mutation to alanine many (ten) of these side chains and analyzed the effects on capsid assembly, stability and conformation, viral DNA encapsidation, and virion infectivity. Mutation of residues Tyr-270, Asp-273, or Asp-474 led to a drastic reduction in infectivity. Mutant Y270A was defective in capsid assembly; mutant D273A formed stable capsids, but it was essentially unable to encapsidate the viral DNA or to externalize the N terminus of the capsid protein VP2, a connected conformational event. Mutation of residues Asp-58, Trp-60, Asn-183, Thr-267, or Lys-471 led to a moderate reduction in infectivity. None of these mutations had an effect on capsid assembly or stability, or on the DNA encapsidation process. However, those five mutant virions were substantially less stable than the parental virion in thermal inactivation assays. The results with this model spherical virus indicate that several capsid residues that are found to be involved in polar interactions or multiple hydrophobic contacts with the viral DNA molecule contribute to preserving the active conformation of the infectious viral particle. Their effect appears to be mediated by the non-covalent interactions they establish with the viral DNA. In addition, at least one acidic residue at each DNA-binding region is needed for DNA packaging.

The crystal structure of a high affinity RNA stem-loop complexed with the bacteriophage MS2 capsid: further challenges in the modeling of ligand-RNA interactions

We have determined the structure to 2.8 A of an RNA aptamer (F5), containing 2'-deoxy-2-aminopurine (2AP) at the -10 position, complexed with MS2 coat protein by soaking the RNA into precrystallised MS2 capsids. The -10 position of the RNA is an important determinant of binding affinity for coat protein. Adenine at this position in other RNA stem-loops makes three hydrogen bonds to protein functional groups. Substituting 2AP for the -10 adenine in the F5 aptamer yields an RNA with the highest yet reported affinity for coat protein. The refined X-ray structure shows that the 2AP base makes an additional hydrogen bond to the protein compared to adenine that is presumably the principal origin of the increased affinity. There are also slight changes in phosphate backbone positions compared to unmodified F5 that probably also contribute to affinity. Such phosphate movements are common in structures of RNAs bound to the MS2 T = 3 protein shell and highlight problems for de novo design of RNA binding ligands.

Rotavirus virus-like particles as surrogates in environmental persistence and inactivation studies

Virus-like particles (VLPs) with the full-length VP2 and VP6 rotavirus capsid proteins, produced in the baculovirus expression system, have been evaluated as surrogates of human rotavirus in different environmental scenarios. Green fluorescent protein-labeled VLPs (GFP-VLPs) and particles enclosing a heterologous RNA (pseudoviruses), whose stability may be monitored by flow cytometry and antigen capture reverse transcription-PCR, respectively, were used. After 1 month in seawater at 20 degrees C, no significant differences were observed between the behaviors of GFP-VLPs and of infectious rotavirus, whereas pseudovirus particles showed a higher decay rate. In the presence of 1 mg of free chlorine (FC)/liter both tracers persisted longer in freshwater at 20 degrees C than infectious viruses, whereas in the presence of 0.2 mg of FC/liter no differences were observed between tracers and infectious rotavirus at short contact times. However, from 30 min of contact with FC onward, the decay of infectious rotavirus was higher than that of recombinant particles. The predicted Ct value for a 90% reduction of GFP-VLPs or pseudoviruses induces a 99.99% inactivation of infectious rotavirus. Both tracers were more resistant to UV light irradiation than infectious rotavirus in fresh and marine water. The effect of UV exposure was more pronounced on pseudovirus than in GFP-VLPs. In all types of water, the UV dose to induce a 90% reduction of pseudovirus ensures a 99.99% inactivation of infectious rotavirus. Recombinant virus surrogates open new possibilities for the systematic validation of virus removal practices in actual field situations where pathogenic agents cannot be introduced.

RNA dynamics in live Escherichia coli cells

We describe a method for tracking RNA molecules in Escherichia coli that is sensitive to single copies of mRNA, and, using the method, we find that individual molecules can be followed for many hours in living cells. We observe distinct characteristic dynamics of RNA molecules, all consistent with the known life history of RNA in prokaryotes: localized motion consistent with the Brownian motion of an RNA polymer tethered to its template DNA, free diffusion, and a few examples of polymer chain dynamics that appear to be a combination of chain fluctuation and chain elongation attributable to RNA transcription. We also quantify some of the dynamics, such as width of the displacement distribution, diffusion coefficient, chain elongation rate, and distribution of molecule numbers, and compare them with known biophysical parameters of the E. coli system.

RNA-binding protein-mediated translational repression of transgene expression in plants

We have demonstrated that RNA-binding proteins from coliphages and yeast can function as translational repressors in plants. RNA sequences called translational operators were inserted at a cap-proximal position in the 5'-UTR of mRNAs of two reporter genes, gus or aroA:CP4. Translation of the reporter mRNAs was efficiently repressed when the RNA binding protein that specifically binds to its cognate operator was co-expressed. The efficiency of translational repression by RNA-binding protein positively correlated with the amount of binding protein in transformed plant cells. Detailed studies on coliphage MS2 coat protein-mediated translational repression also suggested that the efficiency of translational repression was position-dependent. A translational operator situated at the cap-proximal position was more efficient in conferring repression than one that was placed cap-distal. Translational repression can be an efficient means for regulation of transgene expression, thereby broadening current approaches for transgene regulation in plants.

Interaction of C5 protein with RNA aptamers selected by SELEX

RNA aptamers binding to C5 protein, the protein component of Escherichia coli RNase P, were selected and characterized as an initial step in elucidating the mechanism of action of C5 protein as an RNA-binding protein. Sequence analyses of the RNA aptamers suggest that C5 protein binds various RNA molecules with dissociation constants comparable to that of M1 RNA, the RNA component of RNase P. The dominant sequence, W2, was chosen for further study. Interactions between W2 and C5 protein were independent of Mg2+, in contrast to the Mg2+ dependency of M1 RNA-C5 protein interactions. The affinity of W2 for C5 protein increased with increasing concentration of monovalent NH4+, suggesting interactions via hydrophobic attraction. W2 forms a fairly stable complex with C5 protein, although the stability of this complex is lower than that of the complex of M1 RNA with C5 protein. The core RNA motif essential for interaction with C5 protein was identified as a stem-loop structure, comprising a 5 bp stem and a 20 nt loop. Our results strongly imply that C5 protein is an interacting partner protein of some cellular RNA species apart from M1 RNA.

Substitution of an essential adenine in the U1A-RNA complex with a non-polar isostere

The RNA recognition motif (RRM) binds to single-stranded RNA target sites of diverse sequences and structures. A conserved mode of base recognition by the RRM involves the simultaneous formation of a network of hydrogen bonds with the base functional groups and a stacking interaction between the base and a highly conserved aromatic amino acid. We have investigated the energetic contribution of the functional groups involved in the recognition of an essential adenine, A6, in stem-loop 2 of U1 snRNA by the N-terminal RRM of the U1A protein. Previously, we found that elimination of individual hydrogen bond donors and acceptors on A6 destabilized the complex by 0.8-1.9 kcal/mol, while mutation of the aromatic amino acid (Phe56) that stacks with A6 to Ala destabilized the complex by 5.5 kcal/mol. Here we continue to probe the contribution of A6 to complex stability through mutation of both the RNA and protein. We have removed two hydrogen-bonding functional groups by introducing a U1A mutation, Ser91Ala, and replacing A6 with tubercidin, purine, or 1-deazaadenine. We find that the complex is destabilized an additional 1.2-2.6 kcal/mol by the elimination of the second hydrogen bond donor or acceptor. Surprisingly, deletion of all of the functional groups involved in hydrogen bonds with the U1A protein by substituting adenine with 4-methylindole reduced the binding free energy by only 2.0 kcal/mol. Experiments with U1A proteins containing mutations of Phe56 suggested that improved stacking interactions due to the greater hydrophobicity of 4-methylindole than adenine may be partly responsible for the small destabilization of the complex upon substitution of 4-methylindole for A6. The data imply that hydrophobic interactions can compensate energetically for the disruption of the complex hydrogen-bonding network between nucleotide and protein.

Investigation of a conserved stacking interaction in target site recognition by the U1A protein

Three highly conserved aromatic residues in RNA recognition motifs (RRM) participate in stacking interactions with RNA bases upon binding RNA. We have investigated the contribution of one of these aromatic residues, Phe56, to the complex formed between the N-terminal RRM of the spliceosomal protein U1A and stem-loop 2 of U1 snRNA. Previous work showed that the aromatic group is important for high affinity binding. Here we probe how mutation of Phe56 affects the kinetics of complex dissociation, the strength of the hydrogen bonds formed between U1A and the base that stacks with Phe56 (A6) and specific target site recognition. Substitution of Phe56 with Trp or Tyr increased the rate of dissociation of the complex, consistent with previously reported results. However, substitution of Phe56 with His decreased the rate of complex association, implying a change in the initial formation of the complex. Simultaneous modification of residue 56 and A6 revealed energetic coupling between the aromatic group and the functional groups of A6 that hydrogen bond to U1A. Finally, mutation of Phe56 to Leu reduced the ability of U1A to recognize stem-loop 2 correctly. Taken together, these experiments suggest that Phe56 contributes to binding affinity by stacking with A6 and participating in networks of energetically coupled interactions that enable this conserved aromatic amino acid to play a complex role in target site recognition.

Modifying the specificity of an RNA backbone contact

The interaction between the MS2 bacteriophage coat protein homodimer and its cognate RNA hairpin is facilitated by 21 different RNA-protein contacts. In one of these contacts, the 2'-hydroxyl group at ribose -5 of the RNA acts as a hydrogen bond donor to Glu63 in one subunit of the protein. Previous experiments showed that substitution of ribose -5 with deoxyribose resulted in a 24-fold decrease in binding affinity between RNA and protein. Using a protein where the two MS2 monomers were fused to increase stability, the contribution of this contact to the overall binding affinity was investigated by site-directed mutagenesis. When Glu63 was substituted with glutamine, aspartate, or alanine, the binding affinity of the hairpin for the protein was weakened by 12 to 100-fold, similar to that observed with deoxyribose at position -5. However, the specificity of the three mutant proteins for RNAs with various modifications at the 2'-position of ribose -5 differed dramatically. While the Glu63Asp protein resembled the wild-type protein in preferring the 2'-hydroxyl group over a proton or a bulky 2'-substituent, both the Glu63Ala and Glu63Gln proteins preferred bulky 2'-substituents over the 2'-hydroxyl group by more than 100-fold. These experiments emphasize the ease with which the specificity of a protein-nucleic acid interaction can be changed at thermodynamically important sites.

Structural basis of pyrimidine specificity in the MS2 RNA hairpin-coat-protein complex

We have determined the X-ray structures of six MS2 RNA hairpin-coat-protein complexes having five different substitutions at the hairpin loop base -5. This is a uracil in the wild-type hairpin and contacts the coat protein both by stacking on to a tyrosine side chain and by hydrogen bonding to an asparagine side chain. The RNA consensus sequence derived from coat protein binding studies with natural sequence variants suggested that the -5 base needs to be a pyrimidine for strong binding. The five -5 substituents used in this study were 5-bromouracil, pyrimidin-2-one, 2-thiouracil, adenine, and guanine. The structure of the 5-bromouracil complex was determined to 2.2 A resolution, which is the highest to date for any MS2 RNA-protein complex. All the complexes presented here show very similar conformations, despite variation in affinity in solution. The results suggest that the stacking of the -5 base on to the tyrosine side chain is the most important driving force for complex formation. A number of hydrogen bonds that are present in the wild-type complex are not crucial for binding, as they are missing in one or more of the complexes. The results also reveal the flexibility of this RNA-protein interface, with respect to functional group variation, and may be generally applicable to other RNA-protein complexes.

Evaluation of methylphosphonates as analogs for detecting phosphate contacts in RNA-protein complexes

The well-studied interaction between the MS2 coat protein and its cognate hairpin was used to test the utility of the methylphosphonate linkage as a phosphate analog. A nitrocellulose filter binding assay was used to measure the change in binding affinity upon introduction of a single methylphosphonate stereoisomer at 13 different positions in the RNA hairpin. Comparing these data to the available crystal structure of the complex shows that all phosphates that are in proximity to the protein show a weaker binding affinity when substituted with a phosphorothioate and control positions show no change. However, in two cases, a methylphosphonate isomer either increased or decreased the binding affinity where no interaction can be detected in the crystal structure. It is possible that methylphosphonate substitutions at these positions affect the structure or flexibility of the hairpin. The utility of the methylphosphonate substitution is compared to phosphate ethylation and phosphorothioate substitution experiments previously performed on the same system.

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 approximately 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.

Deletion of a single hydrogen bonding atom from the MS2 RNA operator leads to dramatic rearrangements at the RNA-coat protein interface

The MS2 coat protein binds specifically to an RNA hairpin formed within the viral genome. By soaking different RNA fragments into crystals of MS2 coat protein capsids it is possible to determine the X-ray structure of the RNA-protein complexes formed. Here we present the structure to 2.85 A resolution of a complex between a chemically modified RNA hairpin variant and the MS2 coat protein. This RNA variant has a substitution at the -5 base position, which has been shown previously to be pyrimidine-specific and is a uracil in the wild-type RNA. The modified RNA hairpin contains a pyridin-4-one base (4one) at this position that lacks the exocyclic 2-oxygen eliminating the possibility of forming a hydrogen bond to asparagine A87 in the protein. The 4one complex structure shows an unprecedented major conformational change in the loop region of the RNA, whereas there is almost no change in the conformation of the protein.

Interactions of Escherichia coli RNA with bacteriophage MS2 coat protein: genomic SELEX

Genomic SELEX is a method for studying the network of nucleic acid-protein interactions within any organism. Here we report the discovery of several interesting and potentially biologically important interactions using genomic SELEX. We have found that bacteriophage MS2 coat protein binds several Escherichia coli mRNA fragments more tightly than it binds the natural, well-studied, phage mRNA site. MS2 coat protein binds mRNA fragments from rffG (involved in formation of lipopolysaccharide in the bacterial outer membrane), ebgR (lactose utilization repressor), as well as from several other genes. Genomic SELEX may yield experimentally induced artifacts, such as molecules in which the fixed sequences participate in binding. We describe several methods (annealing of oligonucleotides complementary to fixed sequences or switching fixed sequences) to eliminate some, or almost all, of these artifacts. Such methods may be useful tools for both randomized sequence SELEX and genomic SELEX.

Themes in RNA-protein recognition

Atomic resolution structures are now available for more than 20 complexes of proteins with specific RNAs. This review examines two main themes that appear in this set of structures. A "groove binder" class of proteins places a protein structure (alpha-helix, 310-helix, beta-ribbon, or irregular loop) in the groove of an RNA helix, recognizing both the specific sequence of bases and the shape or dimensions of the groove, which are sometimes distorted from the normal A-form. A second class of proteins uses beta-sheet surfaces to create pockets that examine single-stranded RNA bases. Some of these proteins recognize completely unstructured RNA, and in others RNA secondary structure indirectly promotes binding by constraining bases in an appropriate orientation. Thermodynamic studies have shown that binding specificity is generally a function of several factors, including base-specific hydrogen bonds, non-polar contacts, and mutual accommodation of the protein and RNA-binding surfaces. The recognition strategies and structural frameworks used by RNA binding proteins are not exotically different from those employed by DNA-binding proteins, suggesting that the two kinds of nucleic acid-binding proteins have not evolved independently.