Deciphering the rules by which dynamics of mRNA secondary structure affect translation efficiency in Saccharomyces cerevisiae

Messenger RNA (mRNA) secondary structure decreases the elongation rate, as ribosomes must unwind every structure they encounter during translation. Therefore, the strength of mRNA secondary structure is assumed to be reduced in highly translated mRNAs. However, previous studies in vitro reported a positive correlation between mRNA folding strength and protein abundance. The counterintuitive finding suggests that mRNA secondary structure affects translation efficiency in an undetermined manner. Here, we analyzed the folding behavior of mRNA during translation and its effect on translation efficiency. We simulated translation process based on a novel computational model, taking into account the interactions among ribosomes, codon usage and mRNA secondary structures. We showed that mRNA secondary structure shortens ribosomal distance through the dynamics of folding strength. Notably, when adjacent ribosomes are close, mRNA secondary structures between them disappear, and codon usage determines the elongation rate. More importantly, our results showed that the combined effect of mRNA secondary structure and codon usage in highly translated mRNAs causes a short ribosomal distance in structural regions, which in turn eliminates the structures during translation, leading to a high elongation rate. Together, these findings reveal how the dynamics of mRNA secondary structure coupling with codon usage affect translation efficiency.     

The relationship among gene expression, folding free energy and codon usage bias in Escherichia coli

Taking advantage of microarray data in Escherichia coli genome, the relationship among mRNA expression levels, folding free energy and codon usage bias are investigated. Our results indicate that mRNA expression is correlated to the stability of mRNA secondary structure and the codon usage bias. The decrease of the stability of mRNA structure contributes to the increase of mRNA expression. There is a negative correlation between codon adaptation index (CAI) and mRNA expression in genes with less stable structure. The relationship between the stability of mRNA structure and mRNA half-life indicates the stability of mRNA structure is different from mRNA half-life.     

Determining the relationship of gene expression and global mRNA stability in Drosophila melanogaster and Escherichia coli using linear models

There are several sequence-dependent factors regulating gene expression. Some of them have been extensively studied, among the most prominent are GC content and codon usage bias. Other factors hypothesized to have an impact on gene expression are gene length and the thermodynamic stability of mRNA secondary structure. In this work, we analyzed two different microarray datasets of Drosophila melanogaster gene expression and one dataset of Escherichia coli. To investigate the relationship between gene expression, codon usage bias and GC content of first, second and third codon position, gene length and mRNA stability we employed a multiple regression analysis using a comprehensive linear model. It is shown that codon usage bias and GC content of the first, second and third codon position show a significant influence on gene expression, whereas no significant effect of mRNA secondary structure stability is observed.     

mRNA的序列、结构以及翻译速率与蛋白质结构的关系

mRNA所包含的核苷酸序列通过三联体密码子决定了蛋白质的氨基酸序列。但是, 由于对氨基酸同义密码使用频率上的差异, 密码子与反密码子相互作用效率上的不同, 以及密码子上下文关系和mRNA 不同区域二级结构上的差异, 造成了核糖体对mRNA 不同区域翻译速度上的差异, 加之共翻译折叠的作用, 使得mRNA 的序列和结构影响着蛋白质空间结构的形成。     

Translation in Bacillus subtilis: roles and trends of initiation and termination, insights from a genome analysis.

We analysed the Bacillus subtilis protein coding sequences termini, and compared it to other genomes. The analysis focused on signals, com-positional biases of nucleotides, oligonucleotides, codons and amino acids and mRNA secondary structure. AUG is the preferred start codon in all genomes, independent of their G+C content, and seems to induce less stable mRNA structures. However, it is not conserved between homologous genes neither is it preferred in highly expressed genes. In B.subtilis the ribosome binding site is very strong. We found that downstream boxes do not seem to exist either in Escherichia coli or in B.subtilis. UAA stop codon usage is correlated with the G+C content and is strongly selected in highly expressed genes. We found less stable mRNA structures at both termini, which we related to mRNA-ribosome and mRNA-release-factor interactions. This pattern seems to impose a peculiar A-rich nucleotide and codon usage bias in these regions. Finally the analysis of all proteins from B.subtilis revealed a similar amino acid bias near both termini of proteins consisting of over-representation of hydrophilic residues. This bias near the stop codon is partially release-factor specific.     

大肠杆菌mRNA编码区长度、形成二级结构倾向与密码子偏好性的关系

从GenBank获得大肠杆菌K-12MG1655株的全基因组序列,计算了与基因密码子偏好性相关的多个参数(Nc、CAI、GC、GC3s),对其mRNA编码区长度、形成二级结构倾向与密码子偏好性之间的关系进行了统计学分析,发现虽然翻译效率(包括翻译速度和翻译精度)是制约大肠杆菌高表达基因的密码子偏好性的主要因素,同时,mRNA编码区长度及其形成二级结构的倾向也是形成这种偏好性的不可忽略的原因,而且对偏好性有一定程度的削弱。另外对mRNA编码区形成二级结构倾向的生物学意义进行了讨论分析。     

Differential Selective Constraints Shaping Codon Usage Pattern of Housekeeping and Tissue-specific Homologous Genes of Rice and Arabidopsis

Intra-genomic variation between housekeeping and tissue-specific genes has always been a study of interest in higher eukaryotes. To-date, however, no such investigation has been done in plants. Availability of whole genome expression data for both rice and Arabidopsis has made it possible to examine the evolutionary forces in shaping codon usage pattern in both housekeeping and tissue-specific genes in plants. In the present work, we have taken 4065 rice–Arabidopsis homologous gene pairs to study evolutionary forces responsible for codon usage divergence between housekeeping and tissue-specific genes. In both rice and Arabidopsis, it is mutational bias that regulates error minimization in highly expressed genes of both housekeeping and tissue-specific genes. Our results show that, in comparison to tissue-specific genes, housekeeping genes are under strong selective constraint in plants. However, in tissue-specific genes, lowly expressed genes are under stronger selective constraint compared with highly expressed genes. We demonstrated that constraint acting on mRNA secondary structure is responsible for modulating codon usage variations in rice tissue-specific genes. Thus, different evolutionary forces must underline the evolution of synonymous codon usage of highly expressed genes of housekeeping and tissue-specific genes in rice and Arabidopsis.Key words: error minimization, housekeeping, mRNA folding energy, synonymous rates, tissue specific, tRNA copy number     

Reinvestigating the codon and amino acid usage of S. cerevisiae genome: a new insight from protein secondary structure analysis

Biased usage of synonymous codons has been elucidated under the perspective of cellular tRNA abundance for quite a long time now. Taking advantage of publicly available gene expression data for Saccharomyces cerevisiae, a systematic analysis of the codon and amino acid usages in two different coding regions corresponding to the regular (helix and strand) as well as the irregular (coil) protein secondary structures, have been performed. Our analyses suggest that apart from tRNA abundance, mRNA folding stability is another major evolutionary force in shaping the codon and amino acid usage differences between the highly and lowly expressed genes in S. cerevisiae genome and surprisingly it depends on the coding regions corresponding to the secondary structures of the encoded proteins. This is obviously a new paradigm in understanding the codon usage in S. cerevisiae. Differential amino acid usage between highly and lowly expressed genes in the regions coding for the irregular protein secondary structure in S. cerevisiae is expounded by the stability of the mRNA folded structure. Irrespective of the protein secondary structural type, the highly expressed genes always tend to encode cheaper amino acids in order to reduce the overall biosynthetic cost of production of the corresponding protein. This study supports the hypothesis that the tRNA abundance is a consequence of and not a reason for the biased usage of amino acid between highly and lowly expressed genes.     

同义密码子的反常蛋白质二级结构偏好性

统计分析了 119种人蛋白质和 92种大肠杆菌蛋白质的mRNA序列和蛋白质二级结构的关系 .从二肽频数出发 ,研究了同义密码子使用对蛋白质二级结构的影响 ,证明其影响在 10 %到 2 0 %的量级 .对于人和大肠杆菌 ,在 90 %置信水平上 ,4 0 0对二肽中分别有 79对和 6 0对 ,在 95 %置信水平上 ,分别有 4 5对和 36对二肽的相应密码子二联体具有不同于氨基酸的反常二级结构偏好性 ,并且这种反常不能归因于随机涨落     

枯草芽孢杆菌同义密码子使用偏性对蛋白质折叠速率的影响

目前,有关同义密码子使用偏性对蛋白质折叠的影响研究中,样本蛋白均来源于不同的物种。考虑到同义密码子使用偏性的物种差异性,选取枯草杆菌的核蛋白为研究对象。首先,将每条核蛋白按二级结构截取为α螺旋片段、β折叠片段和无规卷曲（α-β混合）片段,并计算其蛋白质折叠速率。然后,整理每个片段相应的核酸序列信息,计算其同义密码子使用度。在此基础上,分析枯草芽孢杆菌核蛋白的同义密码子使用偏性与蛋白质折叠速率的相关性。发现对于不同二级结构的肽链片段,都有部分密码子的使用偏性与其对应的肽链折叠速率显著相关。进一步分析发现,与肽链片段折叠速率显著相关的密码子绝大部分为枯草杆菌全序列或核蛋白序列的每一组同义密码子中使用度最高的密码子。结果表明,在蛋白质的折叠过程中,枯草芽孢杆菌的同义密码子使用偏性起着重要作用。     

Translational features of human alpha 2b interferon production in Escherichia coli

The yield of human alpha 2b interferon in Escherichia coli was optimized by replacement of low-usage arginine codons located in the mRNA 5' end. The differences observed among the various gene variants suggest that codon usage, Shine-Dalgarno-like sequences, and mRNA secondary structure contribute to the performance of E. coli translation machinery.     

Synonymous codon usage in different protein secondary structural classes of human genes: implication for increased non-randomness of GC3 rich genes towards protein stability

The relationship between the synonymous codon usage and different protein secondary structural classes were investigated using 401 Homo sapiens proteins extracted from Protein Data Bank (PDB). A simple Chi-square test was used to assess the significance of deviation of the observed and expected frequencies of 59 codons at the level of individual synonymous families in the four different protein secondary structural classes. It was observed that synonymous codon families show non-randomness in codon usage in four different secondary structural classes. However,when the genes were classified according to their GC3 levels there was an increase in non-randomness in high GC3 group of genes. The non-randomness in codon usage was further tested among the same protein secondary structures belonging to four different protein folding classes of high GC3 group of genes. The results show that in each of the protein secondary structural unit there exist some synonymous family that shows class specific codon-usage pattern. Moreover, there is an increased non-random behaviour of synonymous codons in sheet structure of all secondary structural classes in high GC3 group of genes. Biological implications of these results have been discussed.     

人β神经生长因子基因稀有密码子及mRNA二级结构对其在大肠杆菌中表达的影响

目的:研究人β神经生长因子(β-NGF)基因中稀有密码子及其mRNA二级结构对其在大肠杆菌中表达量的影响。方法:根据对人β-ngf中稀有密码子及其mRNA二级结构的研究,同义突变人β-ngf基因,通过PCR得到人β-ngf的5’端同义突变基因rh-β-ngfp32和全同义突变基因rh-β-ngfmu,将这2个序列克隆入载体pET3a中,得到重组质粒pET3a-ngfp32和pET3a-ngfmu,分别转化大肠杆菌BL21(DE3)感受态,IPTG诱导表达,收集菌体,SDS-PAGE检测其表达量的改变。结果:构建的pET3a-ngfp32和pET3a-ngfmu表达载体酶切和测序结果正确,SDS-PAGE结果显示,与在重组菌pET3a-NGF总蛋白中的表达量相比,目的蛋白rh-β-NGF在重组菌pET3a-NGFP32和pET3a-NGFmu中的表达量均明显增高,并且在重组菌pET3a-NGFmu中的表达量高于重组菌pET3a-NGFP32。结论:目的蛋白rh-β-NGF在重组菌pET3a-NGFP32和pET3a-NGFmu中表达量的增高,说明人β-ngf基因中稀有密码子和mRNA的二级结构对其在大肠杆菌中的表达有较为明显的影响,结果为构建rh-β-NGF的大肠杆菌工程菌株奠定了基础。     

On the origin of synonymous codon usage divergence between thermophilic and mesophilic prokaryotes

Synonymous codon usage analysis between thermophilic and mesophilic prokaryotes has gained wide attention in recent years. Although it is known that thermophilic and mesophilic prokaryotes use different subset of synonymous codons, no reason for this difference is known so far. In the present communication, by analyzing a large number of thermophilic and mesophilic prokaryotes, we provide evidence that bias in the selection of synonymous codons between thermophilic and mesophilic prokaryotes is related to differential folding pattern of mRNA secondary structures. Moreover, we observe that error-minimizing property has significant influence in differentiating the synonymous codon usage between thermophilic and mesophilic prokaryotes. Biological implications of these results are discussed.     

AT2-AT3-profiling: a new look at synonymous codon usage

The teleology of synonymous codon usage (SCU) still awaits a unifying concept. Here the 2nd codon letter of human mRNA-codons was graphically, aided by a computer program, put in relation to the 3rd codon letter, the carrier of SCU: AT2, the density of A+T in 2nd codon position, behaves to AT3, the analogous density of the 3rd codon position, mostly in an inverse fashion that can be expressed as typical figures: mRNAs with an overall AT-density below 50% have a tendency to produce bulky figures called "red dragons" (when redness is attributed to graph-areas, where AT3< AT2), while mRNAs with an AT-density above 50% produce a pattern called "harlequin" consisting of alternating red and blue (blueness, in analogy, when AT3>AT2) diamonds. With more diversion of AT3 from AT2, the harlequin patterns can assume the pattern of a "blue dragon". By analysing the mRNA of known proteins, these patterns can be correlated with certain functional regions: proteins with multiple transmembrane passages show bulky "red dragons", structural proteins with a high glycine- and proline content such as collagen result in "blue dragons". Non-coding mRNAs tend to show a balance between AT2 and AT3 and hence "harlequin patterns". Signal peptides usually code red due to a low AT3 with an AT2-density at the expectance level. With this technique DNA-sequences of as yet unknown functional meaning were scanned. When stretches of harlequin patterns appear interrupted by red or blue dragons, closer scrutiny of these stretches can reveal ORFs which deserve to be looked at more closely for their protein-informational content. At least in humans, SCU appears to follow protein-dependent AT2-density in a reciprocal fashion and does not seem to serve the purpose of influencing mRNA secondary structure which is discussed in depth.