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.     

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.     

Directed Molecular Evolution

We propose the existence of a relationship of stereochemical complementarity between gene sequences that code for interacting components: nucleic acid-nucleic acid, protein-protein and protein-nucleic acid. Such a relationship would impose evolutionary constraints on the DNA sequences themselves, thus retaining these sequences and governing the direction of the evolutionary process. Therefore, we propose that prebiotic, template-directed autocatalytic synthesis of mutally cognate peptides and polynucleotides resulted in their amplification and evolutionary conservation in contemporary prokaryotic and eukaryotic organisms as a genetic regulatory apparatus. If this proposal is correct, then the relationships between the sequences in DNA coding for these interactions constitute a life code of which the genetic code is only one aspect of the many related interactions encoded in DNA.     

An evolutionary conserved pattern of 18S rRNA sequence complementarity to mRNA 5′ UTRs and its implications for eukaryotic gene translation regulation

There are several key mechanisms regulating eukaryotic gene expression at the level of protein synthesis. Interestingly, the least explored mechanisms of translational control are those that involve the translating ribosome per se, mediated for example via predicted interactions between the ribosomal RNAs (rRNAs) and mRNAs. Here, we took advantage of robustly growing large-scale data sets of mRNA sequences for numerous organisms, solved ribosomal structures and computational power to computationally explore the mRNA–rRNA complementarity that is statistically significant across the species. Our predictions reveal highly specific sequence complementarity of 18S rRNA sequences with mRNA 5′ untranslated regions (UTRs) forming a well-defined 3D pattern on the rRNA sequence of the 40S subunit. Broader evolutionary conservation of this pattern may imply that 5′ UTRs of eukaryotic mRNAs, which have already emerged from the mRNA-binding channel, may contact several complementary spots on 18S rRNA situated near the exit of the mRNA binding channel and on the middle-to-lower body of the solvent-exposed 40S ribosome including its left foot. We discuss physiological significance of this structurally conserved pattern and, in the context of previously published experimental results, propose that it modulates scanning of the 40S subunit through 5′ UTRs of mRNAs.     

Long noncoding RNA complementarity and target transcripts abundance

Eukaryotic mRNA metabolism regulates its stability, localization, and translation using complementarity with counter-part RNAs. To modulate their stability, small and long noncoding RNAs can establish complementarity with their target mRNAs. Although complementarity of small interfering RNAs and microRNAs with target mRNAs has been studied thoroughly, partial complementarity of long noncoding RNAs (lncRNAs) with their target mRNAs has not been investigated clearly. To address that research gap, our lab investigated whether the sequence complementarity of two lncRNAs, lincRNA-p21 and OIP5-AS1, influenced the quantity of target RNA expression. We predicted a positive correlation between lncRNA complementarity and target mRNA quantity. We confirmed this prediction using RNA affinity pull down, microarray, and RNA-sequencing analysis. In addition, we utilized the information from this analysis to compare the quantity of target mRNAs when two lncRNAs, lincRNA-p21 and OIP5-AS1, are depleted by siRNAs. We observed that human and mouse lincRNA-p21 regulated target mRNA abundance in complementarity-dependent and independent manners. In contrast, affinity pull down of OIP5-AS1 revealed that changes in OIP5-AS1 expression influenced the amount of some OIP5-AS1 target mRNAs and miRNAs, as we predicted from our sequence complementarity assay. Altogether, the current study demonstrates that partial complementarity of lncRNAs and mRNAs (even miRNAs) assist in determining target RNA expression and quantity.     

Characterization of MicA interactions suggests a potential novel means of gene regulation by small non-coding RNAs

MicA is a small non-coding RNA that regulates ompA mRNA translation in Escherichia coli. MicA has an inhibitory function, base pairing to the translation initiation region of target mRNAs through short sequences of complementarity, blocking their ribosome-binding sites. The MicA structure contains two stem loops, which impede its interaction with target mRNAs, and it is thought that the RNA chaperone protein Hfq, known to be involved in MicA regulation of ompA, may structurally remodel MicA to reveal the ompA-binding site for cognate pairing. To further characterize these interactions, we undertook biochemical and biophysical studies using native MicA and a ‘stabilized’ version, modified to mimic the conformational state of MicA where the ompA-binding site is exposed. Our data corroborate two proposed roles for Hfq: first, to bring both MicA and ompA into close proximity, and second, to restructure MicA to allow exposure of the ompA-binding site for pairing, thereby demonstrating the RNA chaperone function of Hfq. Additionally, at accumulated MicA levels, we identified a Mg²⁺-dependent self-association that occludes the ompA-recognition region. We discuss the potential contribution of an Mg²⁺-mediated conformational switch of MicA for the regulation of MicA function.     

Recognition and discrimination of target mRNAs by Sib RNAs,a cis-encoded sRNA family

Five Sib antitoxin RNAs, members of a family of cis-encoded small regulatory RNAs (sRNAs) in Escherichia coli, repress their target mRNAs, which encode Ibs toxins. This target repression occurs only between cognate sRNA–mRNA pairs with an exception of ibsA. We performed co-transformation assays to assess the ability of SibC derivatives to repress ibsC expression, thereby revealing the regions of SibC that are essential for ibsC mRNA recognition. SibC has two target recognition domains, TRD1 and TRD2, which function independently. The target site for TRD1 is located within the ORF of ibsC, whereas the target site for TRD2 is located in the translational initiation region. The TRD1 sequence is sufficient to repress ibsC expression. In contrast, TRD2 requires a specific structure in addition to the recognition sequence. An in vitro structural probing analysis showed that the initial interactions at these two recognition sites allowed base-pairing to progress into the flanking sequences. Displacement of the TRD1 and TRD2 domains of SibC by the corresponding domains of SibD changed the target specificity of SibC from ibsC to ibsD, suggesting that these two elements modulate the cognate target recognition of each Sib RNA by discriminating among non-cognate ibs mRNAs.     

Protein folding in membranes: insights from neutron diffraction studies of a membrane beta-sheet oligomer

Studies of the assembly of the hexapeptide Acetyl-Trp-Leu₅ (AcWL₅) into β-sheets in membranes have provided insights into membrane protein folding. Yet, the exact structure of the oligomer in the lipid bilayer is unknown. Here we use neutron diffraction to study the disposition of the peptides in bilayers. We find that pairs of adjacent deuterium-labeled leucines have no well-defined peak or dip in the transmembrane distribution profiles, indicative of heterogeneity in the depth of membrane insertion. At the same time, the monomeric homolog AcWL₄ exhibits a homogeneous, well-defined, interfacial location in neutron diffraction experiments. Thus, although the bilayer location of monomeric AcWL₄ is determined by hydrophobicity matching or complementarity within the bilayer, the AcWL₅ molecules in the oligomer are positioned at different depths within the bilayer because they assemble into a staggered transmembrane β-sheet. The AcWL₅ assembly is dominated by protein-protein interactions rather than hydrophobic complementarity. These results have implications for the structure and folding of proteins in their native membrane environment and highlight the importance of the interplay between hydrophobic complementarity and protein-protein interactions in determining the structure of membrane proteins.     

Concentration and Localization of Coexpressed ELAV/Hu Proteins Control Specificity of mRNA Processing

Neuronally coexpressed ELAV/Hu proteins comprise a family of highly related RNA binding proteins which bind to very similar cognate sequences. How this redundancy is linked to in vivo function and how gene-specific regulation is achieved have not been clear. Analysis of mutants in Drosophila ELAV/Hu family proteins ELAV, FNE, and RBP9 and of genetic interactions among them indicates that they have mostly independent roles in neuronal development and function but have converging roles in the regulation of synaptic plasticity. Conversely, ELAV, FNE, RBP9, and human HuR bind ELAV target RNA in vitro with similar affinities. Likewise, all can regulate alternative splicing of ELAV target genes in nonneuronal wing disc cells and substitute for ELAV in eye development upon artificially increased expression; they can also substantially restore ELAV''s biological functions when expressed under the control of the elav gene. Furthermore, ELAV-related Sex-lethal can regulate ELAV targets, and ELAV/Hu proteins can interfere with sexual differentiation. An ancient relationship to Sex-lethal is revealed by gonadal expression of RBP9, providing a maternal fail-safe for dosage compensation. Our results indicate that highly related ELAV/Hu RNA binding proteins select targets for mRNA processing through alteration of their expression levels and subcellular localization but only minimally by altered RNA binding specificity.     

Identification of a novel cognate cytosolic Hsp70 gene (MnHsc70-2) from oriental river prawn Macrobrachium nipponense and comparison of its expressions with the first cognate Hsc70 (MnHsc70-1) under different stresses

The 70-kDa family of heat-shock proteins (Hsp70) plays an important role in the host immunity, which is widely expressed in eukaryotic cells as a major chaperone protein. In the present study, the full-length complementary DNA (cDNA) of a second cognate cytosolic Hsp70 family member (MnHsc70-2) was cloned and characterized from Macrobrachium nipponense, which is an economically and nutritionally important crustacean. The cDNA was 2,717 bp, containing an open reading frame (ORF) of 1,950 bp, which encodes a protein of 649 amino acids with a theoretical molecular weight of 71.1 kDa and an isoelectric point of 5.27. Sequence alignment showed that the MnHsc70-2 shared 75–97 % identity with other heat-shock proteins. Compared to the previously identified cognate Hsp70 (MnHsc70-1) in M. nipponense, MnHsc70-2 showed quite different expression profiles under unstressed conditions in all tested tissues, including the hemocytes, heart, hepatopancreas, gill, intestine, nerve, and muscle. The phylogenetic analysis demonstrated that MnHsc70-2 showed the closest relationship with MnHsc70-1. Heat-inducibility assays showed that two isolated messenger RNAs (mRNAs) displayed different expression profiles in both the hepatopancreas and gill tissues. MnHsc70-1 mRNA expression level decreased at first and then increased to the normal level, whereas MnHsc70-2 mRNA level increased at first and then decreased. The expressions of two MnHsc70s showed substantial obvious heat-inducible regulation in both the hepatopancreas and gill. Under bacterial challenge by Aeromonas hydrophila, both MnHsc70-1 and MnHsc70-2 mRNA level was up-regulated moderately. The results suggested that two cognate Hsc70s may play essential functions in mediating responses to heat-shock and bacterial challenge.     

ARC-1, a sequence element complementary to an internal 18S rRNA segment, enhances translation efficiency in plants when present in the leader or intercistronic region of mRNAs

The sequences of different plant viral leaders with known translation enhancer ability show partial complementarity to the central region of 18S rRNA. Such complementarity might serve as a means to attract 40S ribosomal subunits and explain in part the translation-enhancing property of these sequences. To verify this notion, we designed β-glucuronidase (GUS) mRNAs differing only in the nature of 10 nt inserts in the center of their 41 base leaders. These were complementary to consecutive domains of plant 18S rRNA. Sucrose gradient analysis revealed that leaders with inserts complementary to regions 1105–1114 and 1115–1124 (‘ARC-1’) of plant 18S rRNA bound most efficiently to the 40S ribosomal subunit after dissociation from 80S ribosomes under conditions of high ionic strength, a treatment known to remove translation initiation factors. Using wheat germ cell-free extracts, we could demonstrate that mRNAs with these leaders were translated more than three times more efficiently than a control lacking such a complementarity. Three linked copies of the insert enhanced translation of reporter mRNA to levels comparable with those directed by the natural translation enhancing leaders of tobacco mosaic virus and potato virus Y RNAs. Moreover, inserting the same leaders as intercistronic sequences in dicistronic mRNAs substantially increased translation of the second cistron, thereby revealing internal ribosome entry site activity. Thus, for plant systems, the complementary interaction between mRNA leader and the central region of 18S rRNA allows cap-independent binding of mRNA to the 43S pre-initiation complex without assistance of translation initiation factors.     

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.     

Changes in gene expression during myogenic differentiation: II. Identification of the proteins encoded by myotube-specific complementary DNA sequences

Changes in the pattern of protein synthesis during terminal differentiation of a mouse, myogenic cell line have been examined by two-dimensional gel analysis. In addition the the increase in messenger RNAs coding for the major contractile proteins, several other new proteins are expressed after cell fusion, together with the diminution or loss of proteins expressed in mononucleate cells. In the accompanying paper (Affara et al., 1980), analysis of rnRNA changes using fractionated complementary DNA probes showed that a new group of mRNAs enters the polysomes after cell fusion. To identify which proteins are coded by mRNAs represented in this myotube-specific complementary DNA, we have vised a combination of in vitro translation and sulphydryl chromatography. The results indicate a correspondence between many of the new proteins appearing after cell fusion and the mRNAs encoded by myotube-specific complementary DNA sequences. Included amongst these proteins are some of the contractile polypeptides.     

Intrinsic disorder in protein sense‐antisense recognition

In addition to the well‐established sense‐antisense complementarity abundantly present in the nucleic acid world and serving as a basic principle of the specific double‐helical structure of DNA, production of mRNA, and genetic code‐based biosynthesis of proteins, sense‐antisense complementarity is also present in proteins, where sense and antisense peptides were shown to interact with each other with increased probability. In nucleic acids, sense‐antisense complementarity is achieved via the Watson‐Crick complementarity of the base pairs or nucleotide pairing. In proteins, the complementarity between sense and antisense peptides depends on a specific hydropathic pattern, where codons for hydrophilic and hydrophobic amino acids in a sense peptide are complemented by the codons for hydrophobic and hydrophilic amino acids in its antisense counterpart. We are showing here that in addition to this pattern of the complementary hydrophobicity, sense and antisense peptides are characterized by the complementary order‐disorder patterns and show complementarity in sequence distribution of their disorder‐based interaction sites. We also discuss how this order‐disorder complementarity can be related to protein evolution.