Proteome-wide analysis reveals clues of complementary interactions between mRNAs and their cognate proteins as the physicochemical foundation of the genetic code.

Despite more than 50 years of effort, the origin of the genetic code remains enigmatic. Among different theories, the stereochemical hypothesis suggests that the code evolved as a consequence of direct interactions between amino acids and appropriate bases. If indeed true, such physicochemical foundation of the mRNA/protein relationship could also potentially lead to novel principles of protein-mRNA interactions in general. Inspired by this promise, we have recently explored the connection between the physicochemical properties of mRNAs and their cognate proteins at the proteome level. Using experimentally and computationally derived measures of solubility of amino acids in aqueous solutions of pyrimidine analogs together with knowledge-based interaction preferences of amino acids for different nucleobases, we have revealed a statistically significant matching between the composition of mRNA coding sequences and the base-binding preferences of their cognate protein sequences. Our findings provide strong support for the stereochemical hypothesis of genetic code's origin and suggest the possibility of direct complementary interactions between mRNAs and cognate proteins even in present-day cells.

Evidence of direct complementary interactions between messenger RNAs and their cognate proteins.

Recently, the ability to interact with messenger RNA (mRNA) has been reported for a number of known RNA-binding proteins, but surprisingly also for different proteins without recognizable RNA binding domains including several transcription factors and metabolic enzymes. Moreover, direct binding to cognate mRNAs has been detected for multiple proteins, thus creating a strong impetus to search for functional significance and basic physico-chemical principles behind such interactions. Here, we derive interaction preferences between amino acids and RNA bases by analyzing binding interfaces in the known 3D structures of protein-RNA complexes. By applying this tool to human proteome, we reveal statistically significant matching between the composition of mRNA sequences and base-binding preferences of protein sequences they code for. For example, purine density profiles of mRNA sequences mirror guanine affinity profiles of cognate protein sequences with quantitative accuracy (median Pearson correlation coefficient R = -0.80 across the entire human proteome). Notably, statistically significant anti-matching is seen only in the case of adenine. Our results provide strong evidence for the stereo-chemical foundation of the genetic code and suggest that mRNAs and cognate proteins may in general be directly complementary to each other and associate, especially if unstructured.

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.

Linker and/or transmembrane regions of influenza A/Group-1, A/Group-2, and type B virus hemagglutinins are packed differently within trimers.

Influenza virus hemagglutinin is a homotrimeric spike glycoprotein crucial for virions' attachment, membrane fusion, and assembly reactions. X-ray crystallography data are available for hemagglutinin ectodomains of various types/subtypes but not for anchoring segments. To get structural information for the linker and transmembrane regions of hemagglutinin, influenza A (H1-H16 subtypes except H8 and H15) and B viruses were digested with bromelain or subtilisin Carlsberg, either within virions or in non-ionic detergent micelles. Proteolytical fragments were analyzed by sodium dodecyl sulfate-polyacrylamide gel electrophoresis and matrix-assisted laser desorption/ionization time-of-flight mass spectrometry. Within virions, hemagglutinins of most influenza A/Group-1 and type B virus strains were more susceptible to digestion with bromelain and/or subtilisin compared to A/Group-2 hemagglutinins. The cleavage sites were always located in the hemagglutinin linker sequence. In detergent, 1) bromelain cleaved hemagglutinin of every influenza A subtype in the linker region; 2) subtilisin cleaved Group-2 hemagglutinins in the linker region; 3) subtilisin cleaved Group-1 hemagglutinins in the transmembrane region; 4) both enzymes cleaved influenza B virus hemagglutinin in the transmembrane region. We propose that the A/Group-2 hemagglutinin linker and/or transmembrane regions are more tightly associated within trimers than type A/Group-1 and particularly type B ones. This hypothesis is underpinned by spatial trimeric structure modeling performed for transmembrane regions of both Group-1 and Group-2 hemagglutinin representatives. Differential S-acylation of the hemagglutinin C-terminal anchoring segment with palmitate/stearate residues possibly contributes to fine tuning of transmembrane trimer packing and stabilization since decreased stearate amount correlated with deeper digestion of influenza B and some A/Group-1 hemagglutinins.