Computational analysis of amino acids and their sidechain analogs in crowded solutions of RNA nucleobases with implications for the mRNA–protein complementarity hypothesis

Many critical processes in the cell involve direct binding between RNAs and proteins, making it imperative to fully understand the physicochemical principles behind such interactions at the atomistic level. Here, we use molecular dynamics simulations and 15 μs of sampling to study the behavior of amino acids and amino acid sidechain analogs in high-concentration aqueous solutions of standard RNA nucleobases. Structural and energetic analysis of simulated systems allows us to derive interaction propensity scales for different amino acid/nucleobase combinations. The derived scales closely match and greatly extend the available experimental data, providing a comprehensive foundation for studying RNA–protein interactions in different contexts. By using these scales, we demonstrate a statistically significant connection between nucleobase composition of human mRNA coding sequences and nucleobase interaction propensities of their cognate protein sequences. For example, pyrimidine density profiles of mRNAs match uracil-propensity profiles of their cognate proteins with a median Pearson correlation coefficient of R = −0.70. Our results provide support for the recently proposed hypotheses that mRNAs and their cognate proteins may be physicochemically complementary to each other and bind, especially if unstructured, with the complementarity level being negatively influenced by mRNA adenine content. Finally, we utilize the derived scales to refine the complementarity hypothesis and closely examine its physicochemical underpinnings.     

Sequence signatures of direct complementarity between mRNAs and cognate proteins on multiple levels

A potential connection between physico-chemical properties of mRNAs and cognate proteins, with implications concerning both the origin of the genetic code and mRNA–protein interactions, is unexplored. We compare pyrimidine content of naturally occurring mRNA coding sequences with the propensity of cognate protein sequences to interact with pyrimidines. The latter is captured by polar requirement, a measure of solubility of amino acids in aqueous solutions of pyridines, heterocycles closely related to pyrimidines. We find that the higher the pyrimidine content of an mRNA, the stronger the average propensity of its cognate protein’s amino acids to interact with pyridines. Moreover, window-averaged pyrimidine profiles of individual mRNAs strongly mirror polar-requirement profiles of cognate protein sequences. For example, 4953 human proteins exhibit a correlation between the two with |R| > 0.8. In other words, pyrimidine-rich mRNA regions quantitatively correspond to regions in cognate proteins containing residues soluble in pyrimidine mimetics and vice versa. Finally, by studying randomized genetic code variants we show that the universal genetic code is highly optimized to preserve these correlations. Overall, our findings redefine the stereo-chemical hypothesis concerning code’s origin and provide evidence of direct complementary interactions between mRNAs and cognate proteins before development of ribosomal decoding, but also presently, especially if both are unstructured.     

Absolute binding-free energies between standard RNA/DNA nucleobases and amino-acid sidechain analogs in different environments

Despite the great importance of nucleic acid–protein interactions in the cell, our understanding of their physico-chemical basis remains incomplete. In order to address this challenge, we have for the first time determined potentials of mean force and the associated absolute binding free energies between all standard RNA/DNA nucleobases and amino-acid sidechain analogs in high- and low-dielectric environments using molecular dynamics simulations and umbrella sampling. A comparison against a limited set of available experimental values for analogous systems attests to the quality of the computational approach and the force field used. Overall, our analysis provides a microscopic picture behind nucleobase/sidechain interaction preferences and creates a unified framework for understanding and sculpting nucleic acid–protein interactions in different contexts. Here, we use this framework to demonstrate a strong relationship between nucleobase density profiles of mRNAs and nucleobase affinity profiles of their cognate proteins and critically analyze a recent hypothesis that the two may be capable of direct, complementary interactions.     

Genome-wide analysis of mRNAs targeted to yeast mitochondria

It is agreed that nuclear-encoded mitochondrial proteins are post-translationally targeted to mitochondria, even if, in some cases, a co-translational phase can assist the import of precursor proteins. We used yeast DNA microarrays to analyse the mRNA populations associated with free and mitochondrion-bound polysomes. As expected, many mRNAs, known to encode mitochondrial proteins, are localized to free cytoplasmic polysomes, but many are localized to mitochondrion-bound polysomes. Furthermore, the 3′-UTR of six randomly chosen mitochondrion-bound mRNAs contains sufficient information to target, in vivo, non-translatable RNA to the vicinity of mitochondria. Interestingly, genes producing mRNAs that are targeted to mitochondria are mainly of ancient bacterial origin, whereas those producing mRNAs that are translated in the cytoplasm are mainly of eukaryotic origin. These observations, which support the recent hypotheses concerning the dual origin of the mitochondrial proteome, provide new insights into the biogenesis of mitochondria.     

The 3' untranslated region of human vimentin mRNA interacts with protein complexes containing eEF-1gamma and HAX-1

Previously, we have shown that the vimentin 3′ untranslated region (3′UTR) contains a highly conserved region, which is sufficient for the perinuclear localization of a reporter mRNA. This region was shown to specifically bind protein(s) by band shift analyses. UV-cross-linking studies suggest these proteins are 46- and 35-kDa in mass. Here, we have used this sequence as ‘bait’ to isolate RNA binding proteins using the yeast three-hybrid method. This technique relies on a functional assay detecting bona fide RNA–protein interaction in vivo. Three cDNA isolates, HAX-1, eEF-1γ and hRIP, code for proteins of a size consistent with in vitro cross- linking studies. In all cases, recombinant proteins were capable of binding RNA in vitro. Although hRIP is thought to be a general mRNA binding protein, this represents an unreported activity for eEF-1γ and HAX-1. Moreover, HAX-1 binding appears to be specific to vimentin’s 3′UTR. Both in vivo synthesized eEF-1γ and HAX-1 proteins were ‘pulled out’ of HeLa whole cell extracts by binding to a RNA affinity column comprised of vimentin’s 3′UTR. Moreover, size-fractionation of extracts results in the separation of large complexes containing either eEF-1γ or HAX-1. Thus, in addition to their known functions, both eEF-1γ and HAX-1 are RNA binding proteins, which suggests new roles in mRNA translation and/or perinuclear localization.     

Mutational and Structural Analysis of the Tandem Zinc Finger Domain of Tristetraprolin

Tristetraprolin (TTP), the best known member of a class of tandem (R/K)YKTELCX₈CX₅CX₃H zinc finger proteins, can destabilize target mRNAs by first binding to AU-rich elements (AREs) in their 3′-untranslated regions (UTRs) and subsequently promoting deadenylation and ultimate destruction of those mRNAs. This study sought to determine the roles of selected amino acids in the RNA binding domain, known as the tandem zinc finger (TZF) domain, in the ability of the full-length protein to bind to AREs within the tumor necrosis factor α (TNF) mRNA 3′-UTR. Within the CX₈C region of the TZF domain, mutation of some of the residues specific to TTP, not found in other members of the TTP protein family, resulted in decreased binding to RNA as well as inhibited mRNA deadenylation and decay. Evaluation of simulation solution models revealed a distinct structure in the second zinc finger of TTP that was induced by the presence of these TTP-specific residues. In addition, mutations within the lead-in sequences preceding the first C of highly conserved residues within the CX₅C or CX₃H regions or within the linker region between the two fingers also perturbed both RNA binding and the simulation model of the TZF domain in complex with RNA. We conclude that, although the majority of conserved residues within the TZF domain of TTP are required for productive binding, not all residues at sequence-equivalent positions in the two zinc fingers of the TZF domain of TTP are functionally equivalent.     

In vitro selection of an RNA sequence that interacts with high affinity with thymidylate synthase

Previous studies have shown that the repressive effect of thymidylate synthase (TS) mRNA translation is mediated by direct binding of TS itself to two cis-acting elements on its cognate mRNA. To identify the optimal RNA nucleotides that interact with TS, we in vitro synthesized a completely degenerate, linear RNA pool of 25 nt and employed in vitro selection to isolate high affinity RNA ligands that bind human TS protein. After 10 rounds of selection and amplification, a single RNA molecule was selected that bound TS protein with nearly 20-fold greater affinity than native, wild-type TS RNA sequences. Secondary structure analysis of this RNA sequence predicted it to possess a stem–loop structure. Deletion and/or modification of the UGU loop element within the RNA sequence decreased binding to TS by up to 1000-fold. In vivo transfection experiments revealed that the presence of the selected RNA sequence resulted in a significant increase in the expression of a heterologous luciferase reporter construct in human colon cancer H630 and TS-overexpressing HCT-C:His-TS+ cells, but not in HCT-C18 cells expressing a functionally inactive TS. In addition, the presence of this element in H630 cells leads to induced expression of TS protein. An immunoprecipitation method using RT–PCR confirmed a direct interaction between human TS protein and the selected RNA sequence in transfected human cancer H630 cells. This study identified a novel RNA sequence from a degenerate RNA library that specifically interacts with TS.