Mechanism and thermodynamics of binding of the polypyrimidine tract binding protein to RNA

Reviewing PTBP1 Domain Modularity in the Pre-Genomic Era: A Foundation to Guide the Next Generation of Exploring PTBP1 Structure-Function Relationships

Polypyrimidine tract binding protein 1 (PTBP1) is one of the most well-described RNA binding proteins, known initially for its role as a splicing repressor before later studies revealed its numerous roles in RNA maturation, stability, and translation. While PTBP1's various biological roles have been well-described, it remains unclear how its four RNA recognition motif (RRM) domains coordinate these functions. The early PTBP1 literature saw extensive effort placed in detailing structures of each of PTBP1's RRMs, as well as their individual RNA sequence and structure preferences. However, limitations in high-throughput and high-resolution genomic approaches (i.e., next-generation sequencing had not yet been developed) precluded the functional translation of these findings into a mechanistic understanding of each RRM's contribution to overall PTBP1 function. With the emergence of new technologies, it is now feasible to begin elucidating the individual contributions of each RRM to PTBP1 biological functions. Here, we review all the known literature describing the apo and RNA bound structures of each of PTBP1's RRMs, as well as the emerging literature describing the dependence of specific RNA processing events on individual RRM domains. Our goal is to provide a framework of the structure-function context upon which to facilitate the interpretation of future studies interrogating the dynamics of PTBP1 function.

Allosteric mechanism for SL RNA recognition by polypyrimidine tract binding protein RRM1: An atomistic MD simulation and network-based study

Polypyrimidine tract-binding protein (PTB), an RNA-binding protein, is involved in the regulation of diverse processes in mRNA metabolism. However, the allosteric modulation of its binding with RNA remains unclear. We explore the dynamic characteristics of PTB RNA recognition motif 1 (RRM1) in its RNA-free and wild-type/mutant RNA-bound states to understand the issues using molecular dynamics (MD) simulation, perturbation response scanning (PRS) and protein structure network (PSN) models. It is found that RNA binding strengthens RRM1 stability, while L151G mutation in α3 helix far away from the interface makes the complex unstable. The latter is caused by long-distance dynamic couplings, which makes intermolecular electrostatic and entropy energies unfavorable. The weakened couplings between interface β sheets and C-terminal parts upon mutation reveal RNA recognition is co-regulated by these regions. Interestingly, PRS analysis reveals the allostery caused by the perturbation on α3 helix has already been pre-encoded in the equilibrium dynamics of the protein structure. PSN analysis shows the details of the allosteric signal transmission, revealing the necessity of strong couplings between α3 helix and interface for maintaining the high binding affinity. This study sheds light on the mechanisms of PTB allostery and RNA recognition and can provide important information for drug design.

Role of salt-bridging interactions in recognition of viral RNA by arginine-rich peptides

Interactions between RNA molecules and proteins are critical to many cellular processes and are implicated in various diseases. The RNA-peptide complexes are good model systems to probe the recognition mechanism of RNA by proteins. In this work, we report studies on the binding-unbinding process of a helical peptide from a viral RNA element using nonequilibrium molecular dynamics simulations. We explored the existence of various dissociation pathways with distinct free-energy profiles that reveal metastable states and distinct barriers to peptide dissociation. We also report the free-energy differences for each of the four pathways to be 96.47 ± 12.63, 96.1 ± 10.95, 91.83 ± 9.81, and 92 ± 11.32 kcal/mol. Based on the free-energy analysis, we further propose the preferred pathway and the mechanism of peptide dissociation. The preferred pathway is characterized by the formation of sequential hydrogen-bonding and salt-bridging interactions between several key arginine amino acids and the viral RNA nucleotides. Specifically, we identified one arginine amino acid (R8) of the peptide to play a significant role in the recognition mechanism of the peptide by the viral RNA molecule.

RNA Structural Dynamics As Captured by Molecular Simulations: A Comprehensive Overview

With both catalytic and genetic functions, ribonucleic acid (RNA) is perhaps the most pluripotent chemical species in molecular biology, and its functions are intimately linked to its structure and dynamics. Computer simulations, and in particular atomistic molecular dynamics (MD), allow structural dynamics of biomolecular systems to be investigated with unprecedented temporal and spatial resolution. We here provide a comprehensive overview of the fast-developing field of MD simulations of RNA molecules. We begin with an in-depth, evaluatory coverage of the most fundamental methodological challenges that set the basis for the future development of the field, in particular, the current developments and inherent physical limitations of the atomistic force fields and the recent advances in a broad spectrum of enhanced sampling methods. We also survey the closely related field of coarse-grained modeling of RNA systems. After dealing with the methodological aspects, we provide an exhaustive overview of the available RNA simulation literature, ranging from studies of the smallest RNA oligonucleotides to investigations of the entire ribosome. Our review encompasses tetranucleotides, tetraloops, a number of small RNA motifs, A-helix RNA, kissing-loop complexes, the TAR RNA element, the decoding center and other important regions of the ribosome, as well as assorted others systems. Extended sections are devoted to RNA-ion interactions, ribozymes, riboswitches, and protein/RNA complexes. Our overview is written for as broad of an audience as possible, aiming to provide a much-needed interdisciplinary bridge between computation and experiment, together with a perspective on the future of the field.

Synergy between NMR measurements and MD simulations of protein/RNA complexes: application to the RRMs, the most common RNA recognition motifs

RNA recognition motif (RRM) proteins represent an abundant class of proteins playing key roles in RNA biology. We present a joint atomistic molecular dynamics (MD) and experimental study of two RRM-containing proteins bound with their single-stranded target RNAs, namely the Fox-1 and SRSF1 complexes. The simulations are used in conjunction with NMR spectroscopy to interpret and expand the available structural data. We accumulate more than 50 μs of simulations and show that the MD method is robust enough to reliably describe the structural dynamics of the RRM-RNA complexes. The simulations predict unanticipated specific participation of Arg142 at the protein-RNA interface of the SRFS1 complex, which is subsequently confirmed by NMR and ITC measurements. Several segments of the protein-RNA interface may involve competition between dynamical local substates rather than firmly formed interactions, which is indirectly consistent with the primary NMR data. We demonstrate that the simulations can be used to interpret the NMR atomistic models and can provide qualified predictions. Finally, we propose a protocol for 'MD-adapted structure ensemble' as a way to integrate the simulation predictions and expand upon the deposited NMR structures. Unbiased μs-scale atomistic MD could become a technique routinely complementing the NMR measurements of protein-RNA complexes.

PbWoxT1 mRNA from pear (Pyrus betulaefolia) undergoes long-distance transport assisted by a polypyrimidine tract binding protein

Little is known about the mechanisms by which mRNAs are transported over long distances in the phloem between the rootstock and the scion in grafted woody plants. We identified an mRNA in the pear variety 'Du Li' (Pyrus betulaefolia) that was shown to be transportable in the phloem. It contains a WUSCHEL-RELATED HOMEOBOX (WOX) domain and was therefore named Wox Transport 1 (PbWoxT1). A 548-bp fragment of PbWoxT1 is critical in long-distance transport. PbWoxT1 is rich in CUCU polypyrimidine domains and its mRNAs interact with a polypyrimidine tract binding protein, PbPTB3. Furthermore, the expression of PbWoxT1 significantly increased in the stems of wild-type (WT) tobacco grafted onto the rootstocks of PbWoxT1 or PbPTB3 co-overexpressing lines, but this was not the case in WT plants grafted onto PbWoxT1 overexpressing rootstocks, suggesting that PbPTB3 mediates PbWoxT1 mRNA long-distance transport. We provide novel information that adds a new mechanism with which to explain the noncell-autonomous manner of WOX gene function, which enriches our understanding of how WOX genes work in fruit trees and other species.

Polypyrimidine tract-binding proteins of potato mediate tuberization through an interaction with StBEL5 RNA

Polypyrimidine tract-binding (PTB) proteins are a family of RNA-binding proteins that function in a wide range of RNA metabolic processes by binding to motifs rich in uracils and cytosines. A PTB protein of pumpkin was identified as the core protein of an RNA-protein complex that trafficks RNA. The biological function of the PTB-RNA complex, however, has not been demonstrated. In potato, six PTB proteins have been identified, and two, designated StPTB1 and StPTB6, are similar to the phloem-mobile pumpkin type. RNA binding assays confirmed the interaction of StPTB1 and StPTB6 with discrete pyrimidine-rich sequences of the 3'-untranslated regions of the phloem-mobile mRNA, StBEL5. The promoter of StPTB1 was active in companion cells of phloem in both stem and petioles. Expression of both types was evident in phloem cells of roots and in stolons during tuber formation. RNA accumulation of both PTB proteins was induced by short days in leaves in correlation with enhanced accumulation of StBEL5 RNA. StPTB suppression lines exhibited reduced tuber yields and decreased StBEL5 RNA accumulation, whereas StPTB overexpression lines displayed an increase in tuber production correlated with the enhanced production in stolons of steady-state levels of StBEL5 transcripts and RNA of key tuber identity genes. In StPTB overexpression lines, both the stability and long-distance transport of StBEL5 transcripts were enhanced, whereas in suppression lines stability and transport decreased. Using a transgenic approach, it is shown that the StPTB family of RNA-binding proteins regulate specific stages of development through an interaction with phloem-mobile transcripts of StBEL5.

Transcriptional analysis of phloem-associated cells of potato

Background

Numerous signal molecules, including proteins and mRNAs, are transported through the architecture of plants via the vascular system. As the connection between leaves and other organs, the petiole and stem are especially important in their transport function, which is carried out by the phloem and xylem, especially by the sieve elements in the phloem system. The phloem is an important conduit for transporting photosynthate and signal molecules like metabolites, proteins, small RNAs, and full-length mRNAs. Phloem sap has been used as an unadulterated source to profile phloem proteins and RNAs, but unfortunately, pure phloem sap cannot be obtained in most plant species.

Results

Here we make use of laser capture microdissection (LCM) and RNA-seq for an in-depth transcriptional profile of phloem-associated cells of both petioles and stems of potato. To expedite our analysis, we have taken advantage of the potato genome that has recently been fully sequenced and annotated. Out of the 27 k transcripts assembled that we identified, approximately 15 k were present in phloem-associated cells of petiole and stem with greater than ten reads. Among these genes, roughly 10 k are affected by photoperiod. Several RNAs from this day length-regulated group are also abundant in phloem cells of petioles and encode for proteins involved in signaling or transcriptional control. Approximately 22 % of the transcripts in phloem cells contained at least one binding motif for Pumilio, Nova, or polypyrimidine tract-binding proteins in their downstream sequences. Highlighting the predominance of binding processes identified in the gene ontology analysis of active genes from phloem cells, 78 % of the 464 RNA-binding proteins present in the potato genome were detected in our phloem transcriptome.

Conclusions

As a reasonable alternative when phloem sap collection is not possible, LCM can be used to isolate RNA from specific cell types, and along with RNA-seq, provides practical access to expression profiles of phloem tissue. The combination of these techniques provides a useful approach to the study of phloem and a comprehensive picture of the mechanisms associated with long-distance signaling. The data presented here provide valuable insights into potentially novel phloem-mobile mRNAs and phloem-associated RNA-binding proteins.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1844-2) contains supplementary material, which is available to authorized users.

Patterns of missplicing due to somatic U2AF1 mutations in myeloid neoplasms

Recently, recurrent mutations of spliceosomal genes were frequently identified in myeloid malignancies, as well as other types of cancers. One of these spliceosomal genes, U2AF1, was affected by canonical somatic mutations in aggressive type of myeloid malignancies. We hypothesized that U2AF1 mutations causes defects of splicing (missplicing) in specific genes and that such misspliced genes might be important in leukemogenesis. We analyzed RNA deep sequencing to compare splicing patterns of 201 837 exons between the cases with U2AF1 mutations (n = 6) and wild type (n = 14). We identified different alternative splicing patterns in 35 genes comparing cells with mutant and wild-type U2AF1. U2AF1 mutations are associated with abnormal splicing of genes involved in functionally important pathways, such as cell cycle progression and RNA processing. In addition, many of these genes are somatically mutated or deleted in various cancers. Of note is that the alternative splicing patterns associated with U2AF1 mutations were associated with specific sequence signals at the affected splice sites. These novel observations support the hypothesis that U2AF1 mutations play a significant role in myeloid leukemogenesis due to selective missplicing of tumor-associated genes.

Testing of the GROMOS Force-Field Parameter Set 54A8: Structural Properties of Electrolyte Solutions, Lipid Bilayers, and Proteins

The GROMOS 54A8 force field [Reif et al. J. Chem. Theory Comput.2012, 8, 3705–3723] is the first of its kind to contain nonbonded parameters for charged amino acid side chains that are derived in a rigorously thermodynamic fashion, namely a calibration against single-ion hydration free energies. Considering charged moieties in solution, the most decisive signature of the GROMOS 54A8 force field in comparison to its predecessor 54A7 can probably be found in the thermodynamic equilibrium between salt-bridged ion pair formation and hydration. Possible shifts in this equilibrium might crucially affect the properties of electrolyte solutions or/and the stability of (bio)molecules. It is therefore important to investigate the consequences of the altered description of charged oligoatomic species in the GROMOS 54A8 force field. The present study focuses on examining the ability of the GROMOS 54A8 force field to accurately model the structural properties of electrolyte solutions, lipid bilayers, and proteins. It is found that (i) aqueous electrolytes involving oligoatomic species (sodium acetate, methylammonium chloride, guanidinium chloride) reproduce experimental salt activity derivatives for concentrations up to 1.0 m (1.0-molal) very well, and good agreement between simulated and experimental data is also reached for sodium acetate and methylammonium chloride at 2.0 m concentration, while not even qualitative agreement is found for sodium chloride throughout the whole range of examined concentrations, indicating a failure of the GROMOS 54A7 and 54A8 force-field parameter sets to correctly account for the balance between ion–ion and ion–water binding propensities of sodium and chloride ions; (ii) the GROMOS 54A8 force field reproduces the liquid crystalline-like phase of a hydrated DPPC bilayer at a pressure of 1 bar and a temperature of 323 K, the area per lipid being in agreement with experimental data, whereas other structural properties (volume per lipid, bilayer thickness) appear underestimated; (iii) the secondary structure of a range of different proteins simulated with the GROMOS 54A8 force field at pH 7 is maintained and compatible with experimental NMR data, while, as also observed for the GROMOS 54A7 force field, α-helices are slightly overstabilized with respect to 3₁₀-helices; (iv) with the GROMOS 54A8 force field, the side chains of arginine, lysine, aspartate, and glutamate residues appear slightly more hydrated and present a slight excess of oppositely-charged solution components in their vicinity, whereas salt-bridge formation properties between charged residues at the protein surface, as assessed by probability distributions of interionic distances, are largely equivalent in the GROMOS 54A7 and 54A8 force-field parameter sets.

Thermodynamic basis of selectivity in guide-target-mismatched RNA interference

Silencing in RNAi is strongly affected by guide-strand/target-mRNA mismatches. Target nucleation is thought to occur at positions 2-8 of the guide ("seed region"); successful hybridization in this region is the primary determinant of target-binding affinity and hence target cleavage. To define a molecular basis for the target sequence selectivity in RNAi, we studied all possible distinct single mismatches in seven positions of the seed region-a total of 21 substitutions. We report results from soft-core thermodynamic integration simulations to determine changes in targeting binding-free energies to Argonaute due to single mismatches in the guide strand, which arise during binding of an imperfectly matched target mRNA. In agreement with experiment, most mismatches impair target binding, consistent with a prominent role for binding affinity changes in RNAi sequence selectivity. Individual Argonaute residues located near the mismatched base pair are found to contribute significantly to binding affinity changes. We also use this methodology to analyze the mismatch-dependent free energy changes for dissociation of a DNA•RNA hybrid from Argonaute, as a model for the escape of miRNAs from the silencing pathway. Several mismatched sequences of the miRNA have increased affinity to Argonaute, implying that some mismatches may reduce the probability for escape. Furthermore, calculations of base-substitution-dependent free energy changes for binding ssDNA reveal mild sequence sensitivity as expected for guide strand binding to Argonaute. Our findings give a thermodynamic basis for RNAi target sequence selectivity and suggest that miRNA mismatches may increase silencing effectiveness and thus could be evolutionarily advantageous.

Interactions between PTB RRMs induce slow motions and increase RNA binding affinity

Polypyrimidine tract binding protein (PTB) participates in a variety of functions in eukaryotic cells, including alternative splicing, mRNA stabilization, and internal ribosomal entry site-mediated translation initiation. Its mechanism of RNA recognition is determined in part by the novel geometry of its two C-terminal RNA recognition motifs (RRM3 and RRM4), which interact with each other to form a stable complex (PTB1:34). This complex itself is unusual among RRMs, suggesting that it performs a specific function for the protein. In order to understand the advantage it provides to PTB, the fundamental properties of PTB1:34 are examined here as a comparative study of the complex and its two constituent RRMs. Both RRM3 and RRM4 adopt folded structures that NMR data show to be similar to their structure in PRB1:34. The RNA binding properties of the domains differ dramatically. The affinity of each separate RRM for polypyrimidine tracts is far weaker than that of PTB1:34, and simply mixing the two RRMs does not create an equivalent binding platform. (15)N NMR relaxation experiments show that PTB1:34 has slow, microsecond motions throughout both RRMs including the interdomain linker. This is in contrast to the individual domains, RRM3 and RRM4, where only a few backbone amides are flexible on this time scale. The slow backbone dynamics of PTB1:34, induced by packing of RRM3 and RRM4, could be essential for high-affinity binding to a flexible polypyrimidine tract RNA and also provide entropic compensation for its own formation.

A polypyrimidine tract binding protein, pumpkin RBP50, forms the basis of a phloem-mobile ribonucleoprotein complex

RNA binding proteins (RBPs) are integral components of ribonucleoprotein (RNP) complexes and play a central role in RNA processing. In plants, some RBPs function in a non-cell-autonomous manner. The angiosperm phloem translocation stream contains a unique population of RBPs, but little is known regarding the nature of the proteins and mRNA species that constitute phloem-mobile RNP complexes. Here, we identified and characterized a 50-kD pumpkin (Cucurbita maxima cv Big Max) phloem RNA binding protein (RBP50) that is evolutionarily related to animal polypyrimidine tract binding proteins. In situ hybridization studies indicated a high level of RBP50 transcripts in companion cells, while immunolocalization experiments detected RBP50 in both companion cells and sieve elements. A comparison of the levels of RBP50 present in vascular bundles and phloem sap indicated that this protein is highly enriched in the phloem sap. Heterografting experiments confirmed that RBP50 is translocated from source to sink tissues. Collectively, these findings established that RBP50 functions as a non-cell-autonomous RBP. Protein overlay, coimmunoprecipitation, and cross-linking experiments identified the phloem proteins and mRNA species that constitute RBP50-based RNP complexes. Gel mobility-shift assays demonstrated that specificity, with respect to the bound mRNA, is established by the polypyrimidine tract binding motifs within such transcripts. We present a model for RBP50-based RNP complexes within the pumpkin phloem translocation stream.