Structural mimicry in the phage phi21 N peptide-boxB RNA complex

RPflex: A Coarse-Grained Network Model for RNA Pocket Flexibility Study

RNA regulates various biological processes, such as gene regulation, RNA splicing, and intracellular signal transduction. RNA’s conformational dynamics play crucial roles in performing its diverse functions. Thus, it is essential to explore the flexibility characteristics of RNA, especially pocket flexibility. Here, we propose a computational approach, RPflex, to analyze pocket flexibility using the coarse-grained network model. We first clustered 3154 pockets into 297 groups by similarity calculation based on the coarse-grained lattice model. Then, we introduced the flexibility score to quantify the flexibility by global pocket features. The results show strong correlations between the flexibility scores and root-mean-square fluctuation (RMSF) values, with Pearson correlation coefficients of 0.60, 0.76, and 0.53 in Testing Sets I–III. Considering both flexibility score and network calculations, the Pearson correlation coefficient was increased to 0.71 in flexible pockets on Testing Set IV. The network calculations reveal that the long-range interaction changes contributed most to flexibility. In addition, the hydrogen bonds in the base–base interactions greatly stabilize the RNA structure, while backbone interactions determine RNA folding. The computational analysis of pocket flexibility could facilitate RNA engineering for biological or medical applications.

Early antitermination in the atypical coliphage mEp021 mediated by the Gp17 protein

The coliphage mEp021 belongs to a phage group with a unique immunity repressor, and its life cycle requires the host factor Nus. mEp021 has been classified as non-lambdoid based on its specific characteristics. The mEp021 genome carries a gene encoding an N_λ-like antiterminator protein, termed Gp17, and three nut sites (nutL, nutR1, and nutR2). Analysis of plasmid constructs containing these nut sites, a transcription terminator, and a GFP reporter gene showed high levels of fluorescence when Gp17 was expressed, but not in its absence. Like lambdoid N proteins, Gp17 has an arginine-rich motif (ARM), and mutations in its arginine codons inhibit its function. In infection assays using the mutant phage mEp021ΔGp17::Kan (where gp17 has been deleted), gene transcripts located downstream of transcription terminators were obtained only when Gp17 was expressed. In contrast to phage lambda, mEp021 virus particle production was partially restored (>1/3 relative to wild type) when nus mutants (nusA1, nusB5, nusC60, and nusE71) were infected with mEp021 and Gp17 was overexpressed. Our results suggest that RNA polymerase reads through the third nut site (nutR2), which is more than 7.9 kbp downstream of nutR1.

The glucocorticoid receptor DNA-binding domain recognizes RNA hairpin structures with high affinity

The glucocorticoid receptor (GR) binds the noncoding RNA Gas5 via its DNA-binding domain (DBD) with functional implications in pro-apoptosis signaling. Here, we report a comprehensive in vitro binding study where we have determined that GR-DBD is a robust structure-specific RNA-binding domain. GR-DBD binds to a diverse range of RNA hairpin motifs, both synthetic and biologically derived, with apparent mid-nanomolar affinity while discriminating against uniform dsRNA. As opposed to dimeric recognition of dsDNA, GR-DBD binds to RNA as a monomer and confers high affinity primarily through electrostatic contacts. GR-DBD adopts a discrete RNA-bound state, as assessed by NMR, distinct from both free and DNA-bound. NMR and alanine mutagenesis suggest a heightened involvement of the C-terminal α-helix of the GR-DBD in RNA-binding. RNA competes for binding with dsDNA and occurs in a similar affinity range as dimer binding to the canonical DNA element. Given the prevalence of RNA hairpins within the transcriptome, our findings strongly suggest that many RNAs have potential to impact GR biology.

An intricate balance of hydrogen bonding, ion atmosphere and dynamics facilitates a seamless uracil to cytosine substitution in the U-turn of the neomycin-sensing riboswitch

The neomycin sensing riboswitch is the smallest biologically functional RNA riboswitch, forming a hairpin capped with a U-turn loop-a well-known RNA motif containing a conserved uracil. It was shown previously that a U→C substitution of the eponymous conserved uracil does not alter the riboswitch structure due to C protonation at N3. Furthermore, cytosine is evolutionary permitted to replace uracil in other U-turns. Here, we use molecular dynamics simulations to study the molecular basis of this substitution in the neomycin sensing riboswitch and show that a structure-stabilizing monovalent cation-binding site in the wild-type RNA is the main reason for its negligible structural effect. We then use NMR spectroscopy to confirm the existence of this cation-binding site and to demonstrate its effects on RNA stability. Lastly, using quantum chemical calculations, we show that the cation-binding site is altering the electronic environment of the wild-type U-turn so that it is more similar to the cytosine mutant. The study reveals an amazingly complex and delicate interplay between various energy contributions shaping up the 3D structure and evolution of nucleic acids.

LoaP is a broadly conserved antiterminator protein that regulates antibiotic gene clusters in Bacillus amyloliquefaciens

A valuable resource available in the search for new natural products is the diverse microbial life that spans the planet. A large subset of these microorganisms synthesize complex specialized metabolites exhibiting biomedically important activities. A limiting step to the characterization of these compounds is an elucidation of the genetic regulatory mechanisms that oversee their production. Although proteins that control transcription initiation of specialized metabolite gene clusters have been identified, those affecting transcription elongation have not been broadly investigated. In this study, we analysed the phylogenetic distribution of the large, widespread NusG family of transcription elongation proteins and found that it includes a cohesive outgroup of paralogues (herein coined LoaP), which are often positioned adjacent or within gene clusters for specialized metabolites. We established Bacillus amyloliquefaciens LoaP as a paradigm for this protein subgroup and showed that it regulated the transcriptional readthrough of termination sites located within two different antibiotic biosynthesis operons. Both of these antibiotics have been implicated in plant-protective activities, demonstrating that LoaP controls an important regulon of specialized metabolite genes for this microorganism. These data therefore reveal transcription elongation as a point of regulatory control for specialized metabolite pathways and introduce a subgroup of NusG proteins for this purpose.

HK022 Nun Requires Arginine-Rich Motif Residues Distinct from λ N

Bacteriophage λ N protein binds boxB RNA hairpins in the nut (N utilization) sites of immediate early λ transcripts and interacts with host factors to suppress transcriptional termination at downstream terminators. In opposition to λ N, the Nun protein of HK022 binds the boxBs of coinfecting λ transcripts, interacts with a similar or identical set of host factors, and terminates transcription to suppress λ replication. Comparison of N-boxB and Nun-boxB nuclear magnetic resonance (NMR) structural models suggests similar interactions, though limited mutagenesis of Nun is available. Here, libraries of Nun's arginine-rich motif (ARM) were screened for the ability to exclude λ coinfection, and mutants were assayed for Nun termination with a boxB plasmid reporter system. Several Nun ARM residues appear to be immutable: Asp26, Arg28, Arg29, Arg32, Trp33, and Arg36. Asp26 and Trp33 appear to be unable to contact boxB and are not found at equivalent positions in λ N ARM. To understand if the requirement of Asp26, Trp33, and Arg36 indicated differences between HK022 Nun termination and λ N antitermination complexes, the same Nun libraries were fused to the activation domain of λ N and screened for clones able to complement N-deficient λ. Mutants were assayed for N antitermination. Surprisingly, Asp26 and Trp33 were still essential when Nun ARM was fused to N. Docking suggests that Nun ARM contacts a hydrophobic surface of the NusG carboxy-terminal domain containing residues necessary for Nun function. These findings indicate that Nun ARM relies on distinct contacts in its ternary complex and illustrate how protein-RNA recognition can evolve new regulatory functions.

IMPORTANCE

λ N protein interacts with host factors to allow λ nut-containing transcripts to elongate past termination signals. A competing bacteriophage, HK022, expresses Nun protein, which causes termination of λ nut transcripts. λ N and HK022 Nun use similar arginine-rich motifs (ARMs) to bind the same boxB RNAs in nut transcripts. Screening libraries of Nun ARM mutants, both in HK022 Nun and in a λ N fusion, revealed amino acids essential to Nun that could bind one or more host factors. Docking suggests that NusG, which is present in both Nun termination and N antitermination, is a plausible partner. These findings could help understand how transcription elongation is regulated and illustrate how subtle differences allow ARMs to evolve new regulatory functions.

The structure of the SOLE element of oskar mRNA

mRNA localization by active transport is a regulated process that requires association of mRNPs with protein motors for transport along either the microtubule or the actin cytoskeleton. oskar mRNA localization at the posterior pole of the Drosophila oocyte requires a specific mRNA sequence, termed the SOLE, which comprises nucleotides of both exon 1 and exon 2 and is assembled upon splicing. The SOLE folds into a stem-loop structure. Both SOLE RNA and the exon junction complex (EJC) are required for oskar mRNA transport along the microtubules by kinesin. The SOLE RNA likely constitutes a recognition element for a yet unknown protein, which either belongs to the EJC or functions as a bridge between the EJC and the mRNA. Here, we determine the solution structure of the SOLE RNA by Nuclear Magnetic Resonance spectroscopy. We show that the SOLE forms a continuous helical structure, including a few noncanonical base pairs, capped by a pentanucleotide loop. The helix displays a widened major groove, which could accommodate a protein partner. In addition, the apical helical segment undergoes complex dynamics, with potential functional significance.

The structure of the SOLE element of oskar mRNA

mRNA localization by active transport is a regulated process that requires association of mRNPs with protein motors for transport along either the microtubule or the actin cytoskeleton. oskar mRNA localization at the posterior pole of the Drosophila oocyte requires a specific mRNA sequence, termed the SOLE, which comprises nucleotides of both exon 1 and exon 2 and is assembled upon splicing. The SOLE folds into a stem-loop structure. Both SOLE RNA and the exon junction complex (EJC) are required for oskar mRNA transport along the microtubules by kinesin. The SOLE RNA likely constitutes a recognition element for a yet unknown protein, which either belongs to the EJC or functions as a bridge between the EJC and the mRNA. Here, we determine the solution structure of the SOLE RNA by Nuclear Magnetic Resonance spectroscopy. We show that the SOLE forms a continuous helical structure, including a few noncanonical base pairs, capped by a pentanucleotide loop. The helix displays a widened major groove, which could accommodate a protein partner. In addition, the apical helical segment undergoes complex dynamics, with potential functional significance.

Randomized codon mutagenesis reveals that the HIV Rev arginine-rich motif is robust to substitutions and that double substitution of two critical residues alters specificity

The binding of the arginine-rich motif (ARM) of HIV Rev protein to its high-affinity site in stem IIB in the Rev response element (RRE) initiates assembly of a ribonucleoprotein complex that mediates the export of essential, incompletely spliced viral transcripts. Many biochemical, genetic, and structural studies of Rev-RRE IIB have been published, yet the roles of many peptide residues in Rev ARM are unconfirmed by mutagenesis. Rev aptamer I (RAI) is an optimized RRE IIB that binds Rev with higher affinity and for which mutational data are sparse. Randomized-codon libraries of Rev ARM were assayed for their ability to bind RRE IIB and RAI using a bacterial reporter system based on bacteriophage λ N-nut antitermination. Most Rev ARM residues tolerated substitutions without strong loss of binding to RRE IIB, and all except arginine 39 tolerated substitution without strong loss of binding to RAI. The pattern of critical Rev residues is not the same for RRE IIB and RAI, suggesting important differences between the interactions. The results support and aid the interpretation of existing structural models. Observed clinical variation is consistent with additional constraints on Rev mutation. By chance, we found double mutants of two highly critical residues, arginine 35 (to glycine) and asparagine 40 (to valine or lysine), that bind RRE IIB well, but not RAI. That an apparently distinct binding mode occurs with only two mutations highlights the ability of ARMs to evolve new recognition strategies and supports the application of neutral theories of evolution to protein-RNA recognition.

Understanding the physical basis of the salt dependence of the electrostatic binding free energy of mutated charged ligand-nucleic acid complexes

The predictions of the derivative of the electrostatic binding free energy of a biomolecular complex, ΔG(el), with respect to the logarithm of the 1:1 salt concentration, d(ΔG(el))/d(ln[NaCl]), SK, by the Poisson-Boltzmann equation, PBE, are very similar to those of the simpler Debye-Hückel equation, DHE, because the terms in the PBE's predictions of SK that depend on the details of the dielectric interface are small compared to the contributions from long-range electrostatic interactions. These facts allow one to obtain predictions of SK using a simplified charge model along with the DHE that are highly correlated with both the PBE and experimental binding data. The DHE-based model developed here, which was derived from the generalized Born model, explains the lack of correlation between SK and ΔG(el) in the presence of a dielectric discontinuity, which conflicts with the popular use of this supposed correlation to parse experimental binding free energies into electrostatic and nonelectrostatic components. Moreover, the DHE model also provides a clear justification for the correlations between SK and various empirical quantities, like the number of ion pairs, the ligand charge on the interface, the Coulomb binding free energy, and the product of the charges on the complex's components, but these correlations are weak, questioning their usefulness.

Conformational distributions at the N-peptide/boxB RNA interface studied using site-directed spin labeling

In bacteriophage λ, interactions between a 22-amino acid peptide (called the N-peptide) and a stem-loop RNA element (called boxB) play a critical role in transcription anti-termination. The N-peptide/boxB complex has been extensively studied, and serves as a paradigm for understanding mechanisms of protein/RNA recognition. Particularly, ultrafast spectroscopy techniques have been applied to monitor picosecond fluorescence decay behaviors of 2-aminopurines embedded at various positions of the boxB RNA. The studies have led to a model in which the bound N-peptide exists in dynamic equilibrium between two states, with peptide C-terminal fragment either stacking on (i.e., the stacked state) or peeling away from (i.e., the unstacked state) the RNA loop. The function of the N-peptide/boxB complex seems to correlate with the fraction of the stacked state. Here, the N-peptide/boxB system is studied using the site-directed spin labeling technique, in which X-band electron paramagnetic resonance spectroscopy is applied to monitor nanosecond rotational behaviors of stable nitroxide radicals covalently attached to different positions of the N-peptide. The data reveal that in the nanosecond regime the C-terminal fragment of bound N-peptide adopts multiple discrete conformations within the complex. The characteristics of these conformations are consistent with the proposed stacked and unstacked states, and their distributions vary upon mutations within the N-peptide. These results suggest that the dynamic two-state model remains valid in the nanosecond regime, and represents a unique mode of function in the N-peptide/boxB RNA complex. It also demonstrates a connection between picosecond and nanosecond dynamics in a biological complex.

Binding of the bacteriophage P22 N-peptide to the boxB RNA motif studied by molecular dynamics simulations

Protein-RNA interactions are important for many cellular processes. The Nut-utilization site (N)-protein of bacteriophages contains an N-terminal arginine-rich motif that undergoes a folding transition upon binding to the boxB RNA hairpin loop target structure. Molecular dynamics simulations were used to investigate the dynamics of the P22 N-peptide-boxB complex and to elucidate the energetic contributions to binding. In addition, the free-energy changes of RNA and peptide conformational adaptation to the bound forms, as well as the role of strongly bound water molecules at the peptide-RNA interface, were studied. The influence of peptide amino acid substitutions and the salt dependence of interaction were investigated and showed good agreement with available experimental results. Several tightly bound water molecules were found at the RNA-binding interface in both the presence and absence of N-peptide. Explicit consideration of the waters resulted in shifts of calculated contributions during the energetic analysis, but overall similar binding energy contributions were found. Of interest, it was found that the electrostatic field of the RNA has a favorable influence on the coil-to-alpha-helix transition of the N-peptide already outside of the peptide-binding site. This result may have important implications for understanding peptide-RNA complex formation, which often involves coupled folding and association processes. It indicates that electrostatic interactions near RNA molecules can lead to a shift in the equilibrium toward the bound form of an interacting partner before it enters the binding pocket.

The RNA-binding domain of bacteriophage P22 N protein is highly mutable, and a single mutation relaxes specificity toward lambda

Antitermination in bacteriophage P22, a lambdoid phage, uses the arginine-rich domain of the N protein to recognize boxB RNAs in the nut site of two regulated transcripts. Using an antitermination reporter system, we screened libraries in which each nonconserved residue in the RNA-binding domain of P22 N was randomized. Mutants were assayed for the ability to complement N-deficient virus and for antitermination with P22 boxB(left) and boxB(right) reporters. Single amino acid substitutions complementing P22 N(-) virus were found at 12 of the 13 positions examined. We found evidence for defined structural roles for seven nonconserved residues, which was generally compatible with the nuclear magnetic resonance model. Interestingly, a histidine can be replaced by any other aromatic residue, although no planar partner is obvious. Few single substitutions showed bias between boxB(left) and boxB(right), suggesting that the two RNAs impose similar constraints on genetic drift. A separate library comprising only hybrids of the RNA-binding domains of P22, lambda, and phi21 N proteins produced mutants that displayed bias. P22 N(-) plaque size plotted against boxB(left) and boxB(right) reporter activities suggests that lytic viral fitness depends on balanced antitermination. A few N proteins were able to complement both lambda N- and P22 N-deficient viruses, but no proteins were found to complement both P22 N- and phi21 N-deficient viruses. A single tryptophan substitution allowed P22 N to complement both P22 and lambda N(-). The existence of relaxed-specificity mutants suggests that conformational plasticity provides evolutionary transitions between distinct modes of RNA-protein recognition.

Recognition of a flipped base in a hairpinloop DNA by a small peptide

Two tiny hairpin DNAs, CORE (dAGGCTTCGGCCT) and AP2 (dAGGCTXCGGCCT; X: abasic nucleotide), fold into almost the same tetraloop hairpin structure with one exception, that is, the sixth thymine (T6) of CORE is exposed to the solvent water (Kawakami, J. et al., Chem. Lett. 2001, 258-259). In the present study, we selected small peptides that bind to CORE or AP2 from a combinatorial pentapeptide library with 2.5 x 10(6) variants. On the basis of the structural information, the selected peptide sequences should indicate the essential qualifications for recognition of the hairpin loop DNA with and without a flipped base. In the selected DNA binding peptides, aromatic amino acids such as histidine for CORE and glutamine/aspartic acid for AP2 were found to be abundant amino acids. This amino acid preference suggests that CORE-binding peptides use pi-pi stacking to recognize the target while hydrogen bonding is dominant for AP2-binding peptides. To investigate the binding properties of the selected peptide to the target, surface plasmon resonance was used. The binding constant of the interaction between CORE and a CORE-binding peptide (HWHHE) was about 1.1 x 10(6) M(-1) at 25 degrees C and the resulting binding free energy change at 25 degrees C (DeltaG degrees (25)) was -8.2 kcal mol(-1). The binding of the peptide to AP2 was also analyzed and the resulting binding constant and DeltaG degrees (25) were about 4.2 x 10(4) M(-1) and -6.3 kcal mol(-1), respectively. The difference in the binding free energy changes (DeltaDeltaG degrees (25)) of 1.9 kcal mol(-1) was comparable to the values reported in other systems and was considered a consequence of the loss of pi-pi stacking. Moreover, the stabilization effect by stacking affected the dissociation step as well as the association step. Our results suggest that the existence of an aromatic ring (T6 base) produces new dominant interactions between peptides and nucleic acids, although hydrogen bonding is the preferable mode of interaction in the absence of the flipping base. These findings regarding CORE and AP2 recognition are expected to give useful information in the design of novel artificial DNA binding peptides.

Bacteriophage P22 antitermination boxB sequence requirements are complex and overlap with those of lambda

Transcription antitermination in phages lambda and P22 uses N proteins that bind to similar boxB RNA hairpins in regulated transcripts. In contrast to the lambda N-boxB interaction, the P22 N-boxB interaction has not been extensively studied. A nuclear magnetic resonance structure of the P22 N peptide boxB(left) complex and limited mutagenesis have been reported but do not reveal a consensus sequence for boxB. We have used a plasmid-based antitermination system to screen boxBs with random loops and to test boxB mutants. We find that P22 N requires boxB to have a GNRA-like loop with no simple requirements on the remaining sequences in the loop or stem. U:A or A:U base pairs are strongly preferred adjacent to the loop and appear to modulate N binding in cooperation with the loop and distal stem. A few GNRA-like hexaloops have moderate activity. Some boxB mutants bind P22 and lambda N, indicating that the requirements imposed on boxB by P22 N overlap those imposed by lambda N. Point mutations can dramatically alter boxB specificity between P22 and lambda N. A boxB specific for P22 N can be mutated to lambda N specificity by a series of single mutations via a bifunctional intermediate, as predicted by neutral theories of evolution.

An allosteric path to transcription termination

Transcription termination signals in bacteria occur in RNA as a strong hairpin followed by a stretch of U residues at the 3' terminus. To release the transcript, RNA polymerase (RNAP) is thought to translocate forward without RNA synthesis. Here we provide genetic and biochemical evidence supporting an alternative model in which extensive conformational changes across the enzyme lead to termination without forward translocation. In this model, flexible parts of the RNA exit channel (zipper, flap, and zinc finger) assist the initial step of hairpin folding (nucleation). The hairpin then invades the RNAP main channel, causing RNA:DNA hybrid melting, structural changes of the catalytic site, and DNA-clamp opening induced by interaction with the G(trigger)-loop. Our results envision the elongation complex as a flexible structure, not a rigid body, and establish basic principles of the termination pathway that are likely to be universal in prokaryotic and eukaryotic systems.

A stem-loop structure in the wingless transcript defines a consensus motif for apical RNA transport

Although the subcellular localization of mRNA transcripts is a well-established mechanism for controlling protein localization, the basis for the recognition of mRNA localization elements is only now emerging. For example, although localization elements have been defined for many mRNAs that localize to apical cytoplasm in Drosophila embryos, no unifying properties have been identified within these elements. In this study, we identify and characterize an apical localization element in the 3'UTR of the Drosophila wingless mRNA. We show that this element, referred to as WLE3, is both necessary and sufficient for apical RNA transport. Full, unrestricted activity, however, requires the presence of one of several downstream potentiating elements. Comparison of WLE3 sequences within the Drosophila genus, and their predicted secondary structures, defines a highly conserved stem-loop structure. Despite these high levels of sequence and predicted structure conservation, however, mutagenesis shows significant leeway for both sequence and structure variation in the predicted stem-loop. Importantly, the features that emerge as crucial include an accessible distal helix sequence motif, which is also found in the predicted structures of other apical localization elements.

A novel stable RNA pentaloop that interacts specifically with a motif peptide of lambda-N protein

To achieve a novel specific peptide-nucleic acid binding model, we designed an in vitro selection procedure to decrease the energetic contribution of the electrostatic interaction in the total binding energy and to increase the contribution of hydrogen bonding and pi-pi stacking. After the selection of hairpin-loop RNAs that specifically bound to a model peptide of lambda N protein (N peptide), a new thermostable pentaloop RNA motif (N binding thermostable RNA hairpin: NTS RNA) was revealed. The obtained NTS RNA was able to bind to the N peptide with superior specificity to the boxB RNA, which is the naturally occurring partner of the lambda N protein.

Morphing molecular specificities between Arm-peptide and NUT-RNA in the antitermination complexes of bacteriophages lambda and P22

Bacteriophage lambda's N-protein includes a 17-amino-acid segment, Arm, rich in arginine and having specific affinity for a 15-nucleotide RNA stem-loop called BOX-B. Parallel but different Arm/BOX-B sequences in lambda's cousin, phage P22, account for some of the type specificity that distinguishes lambda from P22: the N of each works only with its cognate BOX-B in vivo. We find that the specificity of N(lambda) can be shifted gradually to that of N(22) by substituting sets of particular amino acids from Arm(22) into Arm of N(lambda). The determinative amino acids are generally those shown by nuclear magnetic resonance to contact BOX-B RNA; gain or loss of these contact amino acids is reasonably expected to contribute to the affinity of each amino acid sequence. Intermediate sequences may show no function with either BOX-B, or weak function with both BOX-B(lambda) and BOX-B(22), the latter suggesting possible evolutionary paths for specificity shifts.

Sequence and structure space of RNA-binding peptides

Studies of RNA-binding peptides, and recent combinatorial library experiments in particular, have demonstrated that diverse peptide sequences and structures can be used to recognize specific RNA sites. The identification of large numbers of sequences capable of binding to a particular site has provided extensive phylogenetic information used to deduce basic principles of recognition. The high frequency at which RNA-binding peptides are found in large sequence libraries suggests plausible routes to evolve sequence-specific binders, facilitating the design of new binding molecules and perhaps reflecting characteristics of natural evolution.