同义密码子使用模式对蛋白产物表达及构象形成的影响*

鉴于遗传密码子的简并性能够将基因遗传信息的容量提升,同义密码子使用偏嗜性得以在生物体的基因组中广泛存在。虽然同义密码子之间碱基的变化并不能导致氨基酸种类的改变,在研究mRNA半衰期、编码多肽翻译效率及肽链空间构象正确折叠的准确性和翻译等这一系列过程中发现,同义密码子使用的偏嗜性在某种程度上通过精微调控翻译机制体现其遗传学功能。同义密码子指导tRNA在翻译过程中识别核糖体的速率变化是由氨基酸的特定顺序决定,并且在新生多肽链合成时,蛋白质共翻译转运机制同时调节其空间构象的正确折叠从而保证蛋白的正常生物学功能。某些同义密码子使用偏嗜性与特定蛋白结构的形成具有显著相关性,密码子使用偏嗜性一旦改变将可能导致新生多肽空间构象出现错误折叠。结合近些年来国内外在此领域的研究成果,阐述同义密码子使用偏嗜性如何发挥精微调控翻译的生物学功能与作用。     

A Synonymous Single Nucleotide Polymorphism in ΔF508 CFTR Alters the Secondary Structure of the mRNA and the Expression of the Mutant Protein

Recent advances in our understanding of translational dynamics indicate that codon usage and mRNA secondary structure influence translation and protein folding. The most frequent cause of cystic fibrosis (CF) is the deletion of three nucleotides (CTT) from the cystic fibrosis transmembrane conductance regulator (CFTR) gene that includes the last cytosine (C) of isoleucine 507 (Ile507ATC) and the two thymidines (T) of phenylalanine 508 (Phe508TTT) codons. The consequences of the deletion are the loss of phenylalanine at the 508 position of the CFTR protein (ΔF508), a synonymous codon change for isoleucine 507 (Ile507ATT), and protein misfolding. Here we demonstrate that the ΔF508 mutation alters the secondary structure of the CFTR mRNA. Molecular modeling predicts and RNase assays support the presence of two enlarged single stranded loops in the ΔF508 CFTR mRNA in the vicinity of the mutation. The consequence of ΔF508 CFTR mRNA “misfolding” is decreased translational rate. A synonymous single nucleotide variant of the ΔF508 CFTR (Ile507ATC), that could exist naturally if Phe-508 was encoded by TTC, has wild type-like mRNA structure, and enhanced expression levels when compared with native ΔF508 CFTR. Because CFTR folding is predominantly cotranslational, changes in translational dynamics may promote ΔF508 CFTR misfolding. Therefore, we propose that mRNA “misfolding” contributes to ΔF508 CFTR protein misfolding and consequently to the severity of the human ΔF508 phenotype. Our studies suggest that in addition to modifier genes, SNPs may also contribute to the differences observed in the symptoms of various ΔF508 homozygous CF patients.     

Manipulating the genetic code for membrane protein production: What have we learnt so far?

With synthetic gene services, molecular cloning is as easy as ordering a pizza. However choosing the right RNA code for efficient protein production is less straightforward, more akin to deciding on the pizza toppings. The possibility to choose synonymous codons in the gene sequence has ignited a discussion that dates back 50years: Does synonymous codon use matter? Recent studies indicate that replacement of particular codons for synonymous codons can improve expression in homologous or heterologous hosts, however it is not always successful. Furthermore it is increasingly apparent that membrane protein biogenesis can be codon-sensitive. Single synonymous codon substitutions can influence mRNA stability, mRNA structure, translational initiation, translational elongation and even protein folding. Synonymous codon substitutions therefore need to be carefully evaluated when membrane proteins are engineered for higher production levels and further studies are needed to fully understand how to select the codons that are optimal for higher production. This article is part of a Special Issue entitled: Protein Folding in Membranes.     

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

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

Mutation-selection models of codon substitution and their use to estimate selective strengths on codon usage

Current models of codon substitution are formulated at the levels of nucleotide substitution and do not explicitly consider the separate effects of mutation and selection. They are thus incapable of inferring whether mutation or selection is responsible for evolution at silent sites. Here we implement a few population genetics models of codon substitution that explicitly consider mutation bias and natural selection at the DNA level. Selection on codon usage is modeled by introducing codon-fitness parameters, which together with mutation-bias parameters, predict optimal codon frequencies for the gene. The selective pressure may be for translational efficiency and accuracy or for fine-tuning translational kinetics to produce correct protein folding. We apply the models to compare mitochondrial and nuclear genes from several mammalian species. Model assumptions concerning codon usage are found to affect the estimation of sequence distances (such as the synonymous rate d(S), the nonsynonymous rate d(N), and the rate at the 4-fold degenerate sites d(4)), as found in previous studies, but the new models produced very similar estimates to some old ones. We also develop a likelihood ratio test to examine the null hypothesis that codon usage is due to mutation bias alone, not influenced by natural selection. Application of the test to the mammalian data led to rejection of the null hypothesis in most genes, suggesting that natural selection may be a driving force in the evolution of synonymous codon usage in mammals. Estimates of selection coefficients nevertheless suggest that selection on codon usage is weak and most mutations are nearly neutral. The sensitivity of the analysis on the assumed mutation model is discussed.     

Codon Usage in Signal Sequences Affects Protein Expression and Secretion Using Baculovirus/Insect Cell Expression System

By introducing synonymous mutations into the coding sequences of GP64sp and FibHsp signal peptides, the influences of mRNA secondary structure and codon usage of signal sequences on protein expression and secretion were investigated using baculovirus/insect cell expression system. The results showed that mRNA structural stability of the signal sequences was not correlated with the protein production and secretion levels, and FibHsp was more tolerable to codon changes than GP64sp. Codon bias analyses revealed that codons for GP64sp were well de-optimized and contained more non-optimal codons than FibHsp. Synonymous mutations in GP64sp sufficiently increased its average codon usage frequency and resulted in dramatic reduction of the activity and secretion of luciferase. Protein degradation inhibition assay with MG-132 showed that higher codon usage frequency in the signal sequence increased the production as well as the degradation of luciferase protein, indicating that the synonymous codon substitutions in the signal sequence caused misfolding of luciferase instead of slowing down the protein production. Meanwhile, we found that introduction of more non-optimal codons into FibHsp could increase the production and secretion levels of luciferase, which suggested a new strategy to improve the production of secretory proteins in insect cells.     

Codon-by-Codon Modulation of Translational Speed and Accuracy Via mRNA Folding

Rapid cell growth demands fast protein translational elongation to alleviate ribosome shortage. However, speedy elongation undermines translational accuracy because of a mechanistic tradeoff. Here we provide genomic evidence in budding yeast and mouse embryonic stem cells that the efficiency–accuracy conflict is alleviated by slowing down the elongation at structurally or functionally important residues to ensure their translational accuracies while sacrificing the accuracy for speed at other residues. Our computational analysis in yeast with codon resolution suggests that mRNA secondary structures serve as elongation brakes to control the speed and hence the fidelity of protein translation. The position-specific effect of mRNA folding on translational accuracy is further demonstrated experimentally by swapping synonymous codons in a yeast transgene. Our findings explain why highly expressed genes tend to have strong mRNA folding, slow translational elongation, and conserved protein sequences. The exquisite codon-by-codon translational modulation uncovered here is a testament to the power of natural selection in mitigating efficiency–accuracy conflicts, which are prevalent in biology.     

Genome-wide selection on codon usage at the population level in the fungal model organism Neurospora crassa

Many organisms exhibit biased codon usage in their genome, including the fungal model organism Neurospora crassa. The preferential use of subset of synonymous codons (optimal codons) at the macroevolutionary level is believed to result from a history of selection to promote translational efficiency. At present, few data are available about selection on optimal codons at the microevolutionary scale, that is, at the population level. Herein, we conducted a large-scale assessment of codon mutations at biallelic sites, spanning more than 5,100 genes, in 2 distinct populations of N. crassa: the Caribbean and Louisiana populations. Based on analysis of the frequency spectra of synonymous codon mutations at biallelic sites, we found that derived (nonancestral) optimal codon mutations segregate at a higher frequency than derived nonoptimal codon mutations in each population; this is consistent with natural selection favoring optimal codons. We also report that optimal codon variants were less frequent in longer genes and that the fixation of optimal codons was reduced in rapidly evolving long genes/proteins, trends suggestive of genetic hitchhiking (Hill-Robertson) altering codon usage variation. Notably, nonsynonymous codon mutations segregated at a lower frequency than synonymous nonoptimal codon mutations (which impair translational efficiency) in each N. crassa population, suggesting that changes in protein composition are more detrimental to fitness than mutations altering translation. Overall, the present data demonstrate that selection, and partly genetic interference, shapes codon variation across the genome in N. crassa populations.     

Isolation of Ribosome Bound Nascent Polypeptides in vitro to Identify Translational Pause Sites Along mRNA

The rate of translational elongation is non-uniform. mRNA secondary structure, codon usage and mRNA associated proteins may alter ribosome movement on the message^{for review see 1}. However, it''s now widely accepted that synonymous codon usage is the primary cause of non-uniform translational elongation rates¹. Synonymous codons are not used with identical frequency. A bias exists in the use of synonymous codons with some codons used more frequently than others². Codon bias is organism as well as tissue specific^2,3. Moreover, frequency of codon usage is directly proportional to the concentrations of cognate tRNAs⁴. Thus, a frequently used codon will have higher multitude of corresponding tRNAs, which further implies that a frequent codon will be translated faster than an infrequent one. Thus, regions on mRNA enriched in rare codons (potential pause sites) will as a rule slow down ribosome movement on the message and cause accumulation of nascent peptides of the respective sizes^5-8. These pause sites can have functional impact on the protein expression, mRNA stability and protein folding^{for review see 9}. Indeed, it was shown that alleviation of such pause sites can alter ribosome movement on mRNA and subsequently may affect the efficiency of co-translational (in vivo) protein folding^1,7,10,11. To understand the process of protein folding in vivo, in the cell, that is ultimately coupled to the process of protein synthesis it is essential to gain comprehensive insights into the impact of codon usage/tRNA content on the movement of ribosomes along mRNA during translational elongation.Here we describe a simple technique that can be used to locate major translation pause sites for a given mRNA translated in various cell-free systems^6-8. This procedure is based on isolation of nascent polypeptides accumulating on ribosomes during in vitro translation of a target mRNA. The rationale is that at low-frequency codons, the increase in the residence time of the ribosomes results in increased amounts of nascent peptides of the corresponding sizes. In vitro transcribed mRNA is used for in vitro translational reactions in the presence of radioactively labeled amino acids to allow the detection of the nascent chains. In order to isolate ribosome bound nascent polypeptide complexes the translation reaction is layered on top of 30% glycerol solution followed by centrifugation. Nascent polypeptides in polysomal pellet are further treated with ribonuclease A and resolved by SDS PAGE. This technique can be potentially used for any protein and allows analysis of ribosome movement along mRNA and the detection of the major pause sites. Additionally, this protocol can be adapted to study factors and conditions that can alter ribosome movement and thus potentially can also alter the function/conformation of the protein.     

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

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

New insights into the interplay between codon bias determinants in plants

Codon bias is the non-random use of synonymous codons, a phenomenon that has been observed in species as diverse as bacteria, plants and mammals. The preferential use of particular synonymous codons may reflect neutral mechanisms (e.g. mutational bias, G|C-biased gene conversion, genetic drift) and/or selection for mRNA stability, translational efficiency and accuracy. The extent to which these different factors influence codon usage is unknown, so we dissected the contribution of mutational bias and selection towards codon bias in genes from 15 eudicots, 4 monocots and 2 mosses. We analysed the frequency of mononucleotides, dinucleotides and trinucleotides and investigated whether the compositional genomic background could account for the observed codon usage profiles. Neutral forces such as mutational pressure and G|C-biased gene conversion appeared to underlie most of the observed codon bias, although there was also evidence for the selection of optimal translational efficiency and mRNA folding. Our data confirmed the compositional differences between monocots and dicots, with the former featuring in general a lower background compositional bias but a higher overall codon bias.     

同义密码子使用模式对多肽链共翻译折叠的影响研究进展

同义密码子使用模式作为核苷酸与氨基酸的纽带,其多样性介导了核糖体扫描速率,同时扩充了基因的遗传信息存储量。随着新型技术的应用,发现特异性密码子和密码子结合力可调节核糖体扫描速率并影响蛋白质构象。同义密码子使用模式通过多种方式在不同环节影响着核糖体扫描速率,同时还影响着自身mRNA的稳定性。本文简述了密码子使用模式如何在核糖体扫描翻译mRNA的过程中实现对多肽链翻译延伸的调控,为今后生物工程学领域如何优化蛋白高效表达提供可参考的思路与理念。     

The relationship between synonymous codon usage and protein structure in Escherichia coli and Homo sapiens

The role of silent position in the codon on the protein structure is an interesting and yet unclear problem. In this paper, 563 Homo sapiens genes and 417 Escherichia coli genes coding for proteins with four different folding types have been analyzed using variance analysis, a multivariate analysis method newly used in codon usage analysis, to find the correlation between amino acid composition, synonymous codon, and protein structure in different organisms. It has been found that in E. coli, both amino acid compositions in differently folded proteins and synonymous codon usage in different gene classes coding for differently folded proteins are significantly different. It was also found that only amino acid composition is different in different protein classes in H. sapiens. There is no universal correlation between synonymous codon usage and protein structure in these two different organisms. Further analysis has shown that GC content on the second codon position can distinguish coding genes for different folded proteins in both organisms.     

In vivo introduction of unpreferred synonymous codons into the Drosophila Adh gene results in reduced levels of ADH protein

The evolution of codon bias, the unequal usage of synonymous codons, is thought to be due to natural selection for the use of preferred codons that match the most abundant species of isoaccepting tRNA, resulting in increased translational efficiency and accuracy. We examined this hypothesis by introducing 1, 6, and 10 unpreferred codons into the Drosophila alcohol dehydrogenase gene (Adh). We observed a significant decrease in ADH protein production with number of unpreferred codons, confirming the importance of natural selection as a mechanism leading to codon bias. We then used this empirical relationship to estimate the selection coefficient (s) against unpreferred synonymous mutations and found the value (s >or= 10(-5)) to be approximately one order of magnitude greater than previous estimates from population genetics theory. The observed differences in protein production appear to be too large to be consistent with current estimates of the strength of selection on synonymous sites in D. melanogaster.     

Human Coding Synonymous Single Nucleotide Polymorphisms at Ramp Regions of mRNA Translation

According to the ramp model of mRNA translation, the first 50 codons favor rare codons and have slower speed of translation. This study aims to detect translational selection on coding synonymous single nucleotide polymorphisms (sSNP) to support the ramp theory. We investigated fourfold degenerate site (FFDS) sSNPs with A↔G or C↔T substitutions in human genome for distribution bias of synonymous codons (SC), grouped by CpG or non-CpG sites. Distribution bias of sSNPs between the 3^rd ∼50^th codons and the 51^st ∼ remainder codons at non-CpG sites were observed. In the 3^rd ∼50^th codons, G→A sSNPs at non-CpG sites are favored than A→G sSNPs [P = 2.89×10⁻³], and C→T at non-CpG sites are favored than T→C sSNPs [P = 8.50×10⁻³]. The favored direction of SC usage change is from more frequent SCs to less frequent SCs. The distribution bias is more obvious in synonymous substitutions CG(G→A), AC(C→T), and CT(C→T). The distribution bias of sSNPs in human genome, i.e. frequent SCs to less frequent SCs is favored in the 3^rd ∼50^th codons, indicates translational selection on sSNPs in the ramp regions of mRNA templates.