Anatomy of noncovalent interactions between the nucleobases or ribose and π-containing amino acids in RNA–protein complexes

A set of >300 nonredundant high-resolution RNA–protein complexes were rigorously searched for π-contacts between an amino acid side chain (W, H, F, Y, R, E and D) and an RNA nucleobase (denoted π–π interaction) or ribose moiety (denoted sugar–π). The resulting dataset of >1500 RNA–protein π-contacts were visually inspected and classified based on the interaction type, and amino acids and RNA components involved. More than 80% of structures searched contained at least one RNA–protein π-interaction, with π–π contacts making up 59% of the identified interactions. RNA–protein π–π and sugar–π contacts exhibit a range in the RNA and protein components involved, relative monomer orientations and quantum mechanically predicted binding energies. Interestingly, π–π and sugar–π interactions occur more frequently with RNA (4.8 contacts/structure) than DNA (2.6). Moreover, the maximum stability is greater for RNA–protein contacts than DNA–protein interactions. In addition to highlighting distinct differences between RNA and DNA–protein binding, this work has generated the largest dataset of RNA–protein π-interactions to date, thereby underscoring that RNA–protein π-contacts are ubiquitous in nature, and key to the stability and function of RNA–protein complexes.     

From Nonspecific DNA–Protein Encounter Complexes to the Prediction of DNA–Protein Interactions

DNA–protein interactions are involved in many essential biological activities. Because there is no simple mapping code between DNA base pairs and protein amino acids, the prediction of DNA–protein interactions is a challenging problem. Here, we present a novel computational approach for predicting DNA-binding protein residues and DNA–protein interaction modes without knowing its specific DNA target sequence. Given the structure of a DNA-binding protein, the method first generates an ensemble of complex structures obtained by rigid-body docking with a nonspecific canonical B-DNA. Representative models are subsequently selected through clustering and ranking by their DNA–protein interfacial energy. Analysis of these encounter complex models suggests that the recognition sites for specific DNA binding are usually favorable interaction sites for the nonspecific DNA probe and that nonspecific DNA–protein interaction modes exhibit some similarity to specific DNA–protein binding modes. Although the method requires as input the knowledge that the protein binds DNA, in benchmark tests, it achieves better performance in identifying DNA-binding sites than three previously established methods, which are based on sophisticated machine-learning techniques. We further apply our method to protein structures predicted through modeling and demonstrate that our method performs satisfactorily on protein models whose root-mean-square Cα deviation from native is up to 5 Å from their native structures. This study provides valuable structural insights into how a specific DNA-binding protein interacts with a nonspecific DNA sequence. The similarity between the specific DNA–protein interaction mode and nonspecific interaction modes may reflect an important sampling step in search of its specific DNA targets by a DNA-binding protein.     

The yeast Pif1p DNA helicase preferentially unwinds RNA DNA substrates

Pif1p is the prototypical member of the PIF1 family of DNA helicases, a subfamily of SFI helicases conserved from yeast to humans. Baker's yeast Pif1p is involved in the maintenance of mitochondrial, ribosomal and telomeric DNA and may also have a general role in chromosomal replication by affecting Okazaki fragment maturation. Here we investigate the substrate preferences for Pif1p. The enzyme was preferentially active on RNA–DNA hybrids, as seen by faster unwinding rates on RNA–DNA hybrids compared to DNA–DNA hybrids. When using forked substrates, which have been shown previously to stimulate the enzyme, Pif1p demonstrated a preference for RNA–DNA hybrids. This preferential unwinding could not be correlated to preferential binding of Pif1p to the substrates that were the most readily unwound. Although the addition of the single-strand DNA-binding protein replication protein A (RPA) stimulated the helicase reaction on all substrates, it did not diminish the preference of Pif1p for RNA–DNA substrates. Thus, forked RNA–DNA substrates are the favored substrates for Pif1p in vitro. We discuss these findings in terms of the known biological roles of the enzyme.     

Dimeric assembly of human Suv3 helicase promotes its RNA unwinding function in mitochondrial RNA degradosome for RNA decay

Human Suv3 is a unique homodimeric helicase that constitutes the major component of the mitochondrial degradosome to work cooperatively with exoribonuclease PNPase for efficient RNA decay. However, the molecular mechanism of how Suv3 is assembled into a homodimer to unwind RNA remains elusive. Here, we show that dimeric Suv3 preferentially binds to and unwinds DNA–DNA, DNA–RNA, and RNA–RNA duplexes with a long 3′ overhang (≥10 nucleotides). The C‐terminal tail (CTT)‐truncated Suv3 (Suv3ΔC) becomes a monomeric protein that binds to and unwinds duplex substrates with ~six to sevenfold lower activities relative to dimeric Suv3. Only dimeric Suv3, but not monomeric Suv3ΔC, binds RNA independently of ATP or ADP, and is capable of interacting with PNPase, indicating that dimeric Suv3 assembly ensures its continuous association with RNA and PNPase during ATP hydrolysis cycles for efficient RNA degradation. We further determined the crystal structure of the apo‐form of Suv3ΔC, and SAXS structures of dimeric Suv3 and PNPase–Suv3 complex, showing that dimeric Suv3 caps on the top of PNPase via interactions with S1 domains, and forms a dumbbell‐shaped degradosome complex with PNPase. Overall, this study reveals that Suv3 is assembled into a dimeric helicase by its CTT for efficient and persistent RNA binding and unwinding to facilitate interactions with PNPase, promote RNA degradation, and maintain mitochondrial genome integrity and homeostasis.     

Putative intermediary stages for the molecular evolution from a ribozyme to a catalytic RNP

A hypothetical evolutionary pathway from a ribozyme to a catalytic RNA–protein complex (RNP) is proposed and examined. In this hypothesis for an early phase of molecular evolution, one RNA–RNA interaction in the starting ribozyme is replaced with an RNA–protein interaction via two intermediary stages. At each stage, the original RNA–RNA interaction and a newly introduced RNA–protein interaction are designed to coexist. The catalytic RNPs corresponding to the intermediary stages were constructed by employing the Tetrahymena ribozyme together with molecular modeling. Analyses of the RNPs indicate that the protein can fully replace the original role of the RNA–RNA interaction in the starting ribozyme and that the association of a protein with a ribozyme might be beneficial for improving the ribozymatic activity.     

Visualization and characterization of RNA–protein interactions in living cells

RNA–protein interactions are the structural and functional basis of significant numbers of RNA molecules. RNA–protein interaction assays though, still mainly depend on biochemical tests in vitro. Here, we establish a convenient and reliable RNA fluorescent three-hybrid (rF3H) method to detect/interrogate the interactions between RNAs and proteins in cells. A GFP tagged highly specific RNA trap is constructed to anchor the RNA of interest to an artificial or natural subcellular structure, and RNA–protein interactions can be detected and visualized by the enrichment of RNA binding proteins (RBPs) at these structures. Different RNA trapping systems are developed and detection of RNA–protein complexes at multiple subcellular structures are assayed. With this new toolset, interactions between proteins and mRNA or noncoding RNAs are characterized, including the interaction between a long noncoding RNA and an epigenetic modulator. Our approach provides a flexible and reliable method for the characterization of RNA–protein interactions in living cells.     

Insights into structure, dynamics and hydration of locked nucleic acid (LNA) strand-based duplexes from molecular dynamics simulations

Locked nucleic acid (LNA) is a chemically modified nucleic acid with its sugar ring locked in an RNA-like (C3′-endo) conformation. LNAs show extraordinary thermal stabilities when hybridized with DNA, RNA or LNA itself. We performed molecular dynamics simulations on five isosequential duplexes (LNA–DNA, LNA–LNA, LNA–RNA, RNA–DNA and RNA–RNA) in order to characterize their structure, dynamics and hydration. Structurally, the LNA–DNA and LNA–RNA duplexes are found to be similar to regular RNA–DNA and RNA–RNA duplexes, whereas the LNA–LNA duplex is found to have its helix partly unwound and does not resemble RNA–RNA duplex in a number of properties. Duplexes with an LNA strand have on average longer interstrand phosphate distances compared to RNA–DNA and RNA–RNA duplexes. Furthermore, intrastrand phosphate distances in LNA strands are found to be shorter than in DNA and slightly shorter than in RNA. In case of induced sugar puckering, LNA is found to tune the sugar puckers in partner DNA strand toward C3′-endo conformations more efficiently than RNA. The LNA–LNA duplex has lesser backbone flexibility compared to the RNA–RNA duplex. Finally, LNA is less hydrated compared to DNA or RNA but is found to have a well-organized water structure.     

The role of disorder in RNA binding affinity and specificity

Most RNA-binding modules are small and bind few nucleotides. RNA-binding proteins typically attain the physiological specificity and affinity for their RNA targets by combining several RNA-binding modules. Here, we review how disordered linkers connecting RNA-binding modules govern the specificity and affinity of RNA–protein interactions by regulating the effective concentration of these modules and their relative orientation. RNA-binding proteins also often contain extended intrinsically disordered regions that mediate protein–protein and RNA–protein interactions with multiple partners. We discuss how these regions can connect proteins and RNA resulting in heterogeneous higher-order assemblies such as membrane-less compartments and amyloid-like structures that have the characteristics of multi-modular entities. The assembled state generates additional RNA-binding specificity and affinity properties that contribute to further the function of RNA-binding proteins within the cellular environment.     

Directed evolution of orthogonal RNA–RBP pairs through library-vs-library in vitro selection

RNA-binding proteins (RBPs) and their RNA ligands play many critical roles in gene regulation and RNA processing in cells. They are also useful for various applications in cell biology and synthetic biology. However, re-engineering novel and orthogonal RNA–RBP pairs from natural components remains challenging while such synthetic RNA–RBP pairs could significantly expand the RNA–RBP toolbox for various applications. Here, we report a novel library-vs-library in vitro selection strategy based on Phage Display coupled with Systematic Evolution of Ligands by EXponential enrichment (PD-SELEX). Starting with pools of 1.1 × 10¹² unique RNA sequences and 4.0 × 10⁸ unique phage-displayed L7Ae-scaffold (LS) proteins, we selected RNA–RBP complexes through a two-step affinity purification process. After six rounds of library-vs-library selection, the selected RNAs and LS proteins were analyzed by next-generation sequencing (NGS). Further deconvolution of the enriched RNA and LS protein sequences revealed two synthetic and orthogonal RNA–RBP pairs that exhibit picomolar affinity and >4000-fold selectivity.     

Selective Glutaraldehyde-Mediated Coupling of Proteins to the 3′-Adenine Terminus of Polymerase Chain Reaction Products

Attachment of proteins to the 3′ end of DNA increases stability of the DNA in serum and retards clearance of DNA by major organs, thereby enhancing in vivo half-life and therapeutic potential of DNA. Unfortunately, the length of DNA molecules that can be produced with 3 ′ modifications by solid-phase synthesis for protein attachment is limited to 45–60 nucleotides due to uncertainties about sequence fidelity for longer oligonucleotides. Here we describe selective covalent coupling of proteins or other molecules to the 3′-adenine overhang of unlabeled and fluorophore-labeled double-stranded polymerase chain reaction products putatively at the N⁶ position of adenine using 2.5% glutaraldehyde at pH 6.0 and 4°C for at least 16 h. Gel mobility shift analyses and fluorescence analyses of the shifted bands supported conjugate formation between double-stranded polymerase chain reaction products and β2-microglobulin. In addition, blunt-ended DNA ladder fragments treated with glutaraldehyde at 4°C showed no evidence of DNA–DNA or DNA–protein conjugate formation. With the present cold glutaraldehyde technique, longer DNA–3′-protein conjugates might be easily mass-produced. The protein portion of a DNA–3′-protein conjugate could possess functionality as well, such as receptor binding for cell entry, cytotoxicity, or opsonization.     

Incorporation of thio-pseudoisocytosine into triplex-forming peptide nucleic acids for enhanced recognition of RNA duplexes

Peptide nucleic acids (PNAs) have been developed for applications in biotechnology and therapeutics. There is great potential in the development of chemically modified PNAs or other triplex-forming ligands that selectively bind to RNA duplexes, but not single-stranded regions, at near-physiological conditions. Here, we report on a convenient synthesis route to a modified PNA monomer, thio-pseudoisocytosine (L), and binding studies of PNAs incorporating the monomer L. Thermal melting and gel electrophoresis studies reveal that L-incorporated 8-mer PNAs have superior affinity and specificity in recognizing the duplex region of a model RNA hairpin to form a pyrimidine motif major-groove RNA₂–PNA triplex, without appreciable binding to single-stranded regions to form an RNA–PNA duplex or, via strand invasion, forming an RNA–PNA₂ triplex at near-physiological buffer condition. In addition, an L-incorporated 8-mer PNA shows essentially no binding to single-stranded or double-stranded DNA. Furthermore, an L-modified 6-mer PNA, but not pseudoisocytosine (J) modified or unmodified PNA, binds to the HIV-1 programmed −1 ribosomal frameshift stimulatory RNA hairpin at near-physiological buffer conditions. The stabilization of an RNA₂–PNA triplex by L modification is facilitated by enhanced van der Waals contacts, base stacking, hydrogen bonding and reduced dehydration energy. The destabilization of RNA–PNA and DNA–PNA duplexes by L modification is due to the steric clash and loss of two hydrogen bonds in a Watson–Crick-like G–L pair. An RNA₂–PNA triplex is significantly more stable than a DNA₂–PNA triplex, probably because the RNA duplex major groove provides geometry compatibility and favorable backbone–backbone interactions with PNA. Thus, L-modified triplex-forming PNAs may be utilized for sequence-specifically targeting duplex regions in RNAs for biological and therapeutic applications.     

Inferring primase-DNA specific recognition using a data driven approach

DNA–protein interactions play essential roles in all living cells. Understanding of how features embedded in the DNA sequence affect specific interactions with proteins is both challenging and important, since it may contribute to finding the means to regulate metabolic pathways involving DNA–protein interactions. Using a massive experimental benchmark dataset of binding scores for DNA sequences and a machine learning workflow, we describe the binding to DNA of T7 primase, as a model system for specific DNA–protein interactions. Effective binding of T7 primase to its specific DNA recognition sequences triggers the formation of RNA primers that serve as Okazaki fragment start sites during DNA replication.     

Controlling RNA self-assembly to form filaments

Fundamental control over supra-molecular self-assembly for organization of matter on the nano-scale is a major objective of nanoscience and nanotechnology. ‘RNA tectonics’ is the design of modular RNA units, called tectoRNAs, that can be programmed to self-assemble into novel nano- and meso-scopic architectures of desired size and shape. We report the three-dimensional design of tectoRNAs incorporating modular 4-way junction (4WJ) motifs, hairpin loops and their cognate loop–receptors to create extended, programmable interaction interfaces. Specific and directional RNA–RNA interactions at these interfaces enable conformational, topological and orientational control of tectoRNA self-assembly. The interacting motifs are precisely positioned within the helical arms of the 4WJ to program assembly from only one helical stacking conformation of the 4WJ. TectoRNAs programmed to assemble with orientational compensation produce micrometer-scale RNA filaments through supra-molecular equilibrium polymerization. As visualized by transmission electron microscopy, these RNA filaments resemble actin filaments from the protein world. This work emphasizes the potential of RNA as a scaffold for designing and engineering new controllable biomaterials mimicking modern cytoskeletal proteins.     

Statistical prediction of single-stranded regions in RNA secondary structure and application to predicting effective antisense target sites and beyond

Single-stranded regions in RNA secondary structure are important for RNA–RNA and RNA–protein interactions. We present a probability profile approach for the prediction of these regions based on a statistical algorithm for sampling RNA secondary structures. For the prediction of phylogenetically-determined single-stranded regions in secondary structures of representative RNA sequences, the probability profile offers substantial improvement over the minimum free energy structure. In designing antisense oligonucleotides, a practical problem is how to select a secondary structure for the target mRNA from the optimal structure(s) and many suboptimal structures with similar free energies. By summarizing the information from a statistical sample of probable secondary structures in a single plot, the probability profile not only presents a solution to this dilemma, but also reveals ‘well-determined’ single-stranded regions through the assignment of probabilities as measures of confidence in predictions. In antisense application to the rabbit β-globin mRNA, a significant correlation between hybridization potential predicted by the probability profile and the degree of inhibition of in vitro translation suggests that the probability profile approach is valuable for the identification of effective antisense target sites. Coupling computational design with DNA–RNA array technique provides a rational, efficient framework for antisense oligonucleotide screening. This framework has the potential for high-throughput applications to functional genomics and drug target validation.     

Identification of antisense RNA stem–loops that inhibit RNA–protein interactions using a bacterial reporter system

Many well-characterized examples of antisense RNAs from prokaryotic systems involve hybridization of the looped regions of stem–loop RNAs, presumably due to the high thermodynamic stability of the resulting loop–loop and loop–linear interactions. In this study, the identification of RNA stem–loops that inhibit U1A protein binding to the hpII RNA through RNA–RNA interactions was attempted using a bacterial reporter system based on phage λ N-mediated antitermination. As a result, loop sequences possessing 7–8 base complementarity to the 5′ region of the boxA element important for functional antitermination complex formation, but not the U1 hpII loop, were identified. In vitro and in vivo mutational analysis strongly suggested that the selected loop sequences were binding to the boxA region, and that the structure of the antisense stem–loop was important for optimal inhibitory activity. Next, in an attempt to demonstrate the ability to inhibit the interaction between the U1A protein and the hpII RNA, the rational design of an RNA stem–loop that inhibits U1A-binding to a modified hpII was carried out. Moderate inhibitory activity was observed, showing that it is possible to design and select antisense RNA stem–loops that disrupt various types of RNA–protein interactions.     

RecJ-like protein from Pyrococcus furiosus has 3′–5′ exonuclease activity on RNA: implications for proofreading of 3′-mismatched RNA primers in DNA replication

Replicative DNA polymerases require an RNA primer for leading and lagging strand DNA synthesis, and primase is responsible for the de novo synthesis of this RNA primer. However, the archaeal primase from Pyrococcus furiosus (Pfu) frequently incorporates mismatched nucleoside monophosphate, which stops RNA synthesis. Pfu DNA polymerase (PolB) cannot elongate the resulting 3′-mismatched RNA primer because it cannot remove the 3′-mismatched ribonucleotide. This study demonstrates the potential role of a RecJ-like protein from P. furiosus (PfRecJ) in proofreading 3′-mismatched ribonucleotides. PfRecJ hydrolyzes single-stranded RNA and the RNA strand of RNA/DNA hybrids in the 3′–5′ direction, and the kinetic parameters (K_m and K_cat) of PfRecJ during RNA strand digestion are consistent with a role in proofreading 3′-mismatched RNA primers. Replication protein A, the single-stranded DNA–binding protein, stimulates the removal of 3′-mismatched ribonucleotides of the RNA strand in RNA/DNA hybrids, and Pfu DNA polymerase can extend the 3′-mismatched RNA primer after the 3′-mismatched ribonucleotide is removed by PfRecJ. Finally, we reconstituted the primer-proofreading reaction of a 3′-mismatched ribonucleotide RNA/DNA hybrid using PfRecJ, replication protein A, Proliferating cell nuclear antigen (PCNA) and PolB. Given that PfRecJ is associated with the GINS complex, a central nexus in archaeal DNA replication fork, we speculate that PfRecJ proofreads the RNA primer in vivo.