全 文 ：·综述与专论· 2015, 31(12):50-55
生物技术通报
BIOTECHNOLOGY BULLETIN
膜蛋白在生物体的许多生命活动中起着非常重
要的作用，如细胞的增殖和分化、能量转换、信号
转导及物质运输等。目前有大约 60% 的药物作用
靶点是膜蛋白［1］。关于膜蛋白结构的研究对探究疾
病的发病机理和有效治疗有着重要意义，也是当前
的研究热点之一。研究膜蛋白结构的技术包括 X 射
线衍射、核磁共振波谱、电子显微镜、原子力显微
镜、红外光谱和圆二色谱等。其中 X 射线衍射和
核磁共振波谱技术是对膜蛋白三维结构进行研究的
主要方法。尤其利用固体核磁共振技术（Solid state
NMR）可在接近膜蛋白的天然环境的磷脂双分子层
中研究膜蛋白的三维结构信息和动力学特征，在膜
蛋白结构研究领域具有独特的优势和广阔的发展前
景［2，3］。采用固体核磁共振研究膜蛋白的结构与功
收稿日期 ：2015-03-03
作者简介 ：王世山，男，硕士研究生，研究方向 ：膜蛋白结构与功能 ；E-mail ：seebinghe@foxmail.com
通讯作者 ：杨俊，E-mail ：yangjun@wipm.ac.cn
膜蛋白原核表达的基因序列优化
王世山  陈艳可  杨俊
（中国科学院武汉物理与数学研究所 武汉磁共振中心 波谱与原子分子物理国家重点实验室，武汉 430071）
摘 要 ： 膜蛋白的结构与功能研究是当前的热点之一。目前获取膜蛋白的最主要途径是通过原核系统进行表达。但是膜蛋
白在原核系统中表达水平相对较低，通过对基因序列进行优化促使膜蛋白正确高效表达是一项非常重要的技术手段。归纳了近些
年来膜蛋白原核表达技术在基因序列优化研究上的最新进展，并从稀有密码子的优化、mRNA 的稳定性与翻译起始、mRNA 与核
糖体行为、翻译的效率与膜蛋白的折叠 4 个方面对基因序列上影响膜蛋白表达的因素进行了总结，旨在为膜蛋白的原核表达提供
一定的思路和参考。
关键词 ： 膜蛋白 ；超量表达 ；基因序列优化 ；密码子优化 ；mRNA 稳定性 ；翻译速率 ；膜蛋白折叠
DOI ：10.13560/j.cnki.biotech.bull.1985.2015.12.007
The Gene Sequence Optimization of Membrane Protein in Prokaryotic
Expression System
Wang Shishan Chen Yanke Yang Jun
（State Key Laboratory of Magnetic Resonance and Atomic and Molecular Physics，Wuhan Centre for Magnetic Resonance，Wuhan Institute of
Physics and Mathematics，Chinese Academy of Sciences，Wuhan 430071）
Abstract: Studies of the structure and function of membrane proteins are hotspots of life science in recent years. The most important
approach to acquire membrane protein is recombinant expression in prokaryotic system. However, the yield of membrane protein in prokaryotic
system is relatively low. Gene sequence optimization is a effective technique to enhance expression of membrane protein in prokaryotic system.
In this article, we summarize the latest progress on gene sequence optimization for the prokaryotic expression of membrane proteins, and discuss
the factors which affect the expression of membrane protein, including the optimization of rare codons, mRNA stability and translation starting,
mRNA and ribosome behaviour, translation rate and membrane protein folding. Our paper provides new possibilities for the over-expression of
membrane protein in prokaryotic system. .
Key words: membrane protein ；over-expression ；gene sequence optimization ；rare codons optimization ；mRNA stability ；translation
rate ；membrane protein folding
2015,31(12) 51王世山等：膜蛋白原核表达的基因序列优化
能，需要制备大量的同位素富集的膜蛋白样品，重
组表达就成为了获取膜蛋白样品的重要手段。常用
于重组膜蛋白的表达系统有真核表达系统、原核表
达系统和近些年来发展的无细胞表达系统。其中以
大肠杆菌（E.coli）为代表的原核表达系统因为操作
简单、成本相对低廉、遗传背景清楚、方便同位素
标记，以及有大量可利用的表达载体和宿主菌株等
原因，是目前获取重组膜蛋白的最主要途径。对于
一些膜蛋白而言，采用增加蛋白可溶性或者促使蛋
白形成包涵体的标签进行融合表达，是很好的增加
蛋白产量的办法，但是目前还没有普遍有效的融合
标签可用于所有膜蛋白的超量表达［4-6］。因此，如
何提高蛋白表达量一直是膜蛋白三维结构研究过程
中首先需要解决的问题。本文归纳了近些年来基因
序列优化在膜蛋白原核表达技术研究上的最新进展，
并从稀有密码子的优化、mRNA 的稳定性与翻译起
始、mRNA 与核糖体行为、翻译的效率与膜蛋白的
折叠 4 个方面对基因序列影响膜蛋白表达的优化因
素进行了总结，旨在对原核系统中膜蛋白的超量表
达提供一些思路和参考。
1 膜蛋白基因序列特点
基因的碱基序列上携带的遗传信息，不仅决定
了所合成肽链的氨基酸组成和顺序，还会影响基因
的转录和翻译过程，这些遗传信息是决定蛋白质丰
度的主要原因［7］。膜蛋白通常含有一个或几个疏水
性很强的跨膜片段，其组成以疏水性氨基酸残基为
主［8］。编码疏水性氨基酸的密码子中，其尿嘧啶含
量要明显高于亲水性氨基酸的密码子［9］，与可溶蛋
白相比膜蛋白在基因组转录中处于明显劣势［10，11］，
这种碱基特异性是可能原因之一。此外，在膜蛋白
基因序列上还存在影响蛋白折叠和插入的信号，它
们对膜蛋白的结构和功能的保持至关重要。对多药
耐药基因的研究发现，同义密码子的替换会影响其
产物 P-糖蛋白折叠和插入膜的时机，从而使 P-糖蛋
白底物特异性发生改变［12］。在啤酒酵母膜蛋白的跨
膜区后第 45 个密码子及第 75 个密码子左右经常会
出现一些翻译暂停信号，由于翻译过程中核糖体通
道正好可以容纳 30-72 个氨基酸，这些信号可能与
新合成的跨膜区的折叠有关［13，14］。
因此，对于外源膜蛋白在原核表达系统中出现
的表达水平低、可溶性差、易错误折叠和聚集的现
象［4，15］，在具体实践中结合表达宿主对目的基因进
行碱基序列优化是促使膜蛋白高效，且正确表达的
一项十分常用的解决手段［4，16，17］。
2 膜蛋白基因序列的优化分析
蛋白的表达是遗传信息从基因的 DNA 序列经过
转录、翻译和肽链折叠而形成蛋白质等一系列的过
程。因此膜蛋白基因序列的优化会涉及到 mRNA 的
转录水平、翻译的效率以及翻译过程中肽链的折叠
等等多个方面。下面将从稀有密码子优化、mRNA
的稳定性与翻译起始、mRNA 与核糖体行为、翻译
的效率与膜蛋白的折叠四个方面进行介绍。
2.1 稀有密码子的优化
除甲硫氨酸和色氨酸之外，其他的氨基酸至
少对应两个密码子。不同物种之间对密码子的偏好
有明显差异，使用频率低的密码子被称为稀有密码
子［18，19］。大量研究表明，稀有密码子的存在会影响
膜蛋白的最终表达水平，因此异源表达的膜蛋白通
常要在基因序列上对密码子进行优化［9，17，20-22］。许
多膜蛋白，如秀丽隐杆线虫氯离子通道蛋白（Caen-
orhabditis elegans GluCl ion channel）、人类细胞色素
P4504F11（Human cytochrome P450 4F11）、G 蛋 白
偶联受体（G-protein-coupled receptor）等，经过密码
子优化后在大肠杆菌中得到了良好的表达［22-25］。除
此之外，稀有密码子的出现在有些情况下还会造成
点突变甚至移码突变，从而表达出序列错误的多肽。
如有些情况下稀有密码子 AGA 编码的精氨酸（Arg）
会被错误翻译成赖氨酸（Lys）［26］，串联的稀有密码
子 如 AGG-AGG、AGA-AGA 及 3 个 连 续 的 CGG 可
能会造成阅读框的改变［27，28］，终止密码子 UGA 被
翻译成色氨酸（Trp）造成通读而产生后续的冗余序
列［29］。因此，应该避免在膜蛋白的基因序列中含有
大量稀有密码子，对此最直接的办法是对基因序列
中的密码子进行优化，目前有 Codon Usage Analyzer、
Graphical Codon Usage Analyzer 等 许 多 在 线 程 序 或
网站可以利用。另一种解决方式是采用带有稀有密
码子额外拷贝的表达菌株，如已经商业化的 BL21
Codon Plus 或 Rosetta 菌株，这在一定程度上也可以
生物技术通报 Biotechnology Bulletin 2015,Vol.31,No.1252
避免或减轻稀有密码子造成的不良影响［20，30］。
值得一提的是，大肠杆菌在不同的培养基中
对密码子的偏好性可能会发生变化，比如 2012 年，
Weissman 课题组发现丝氨酸密码子在 LB 培养基和
含葡萄糖的富营养培养基中表现有很大差异［31］。但
到目前为止，这种关于培养基成分对大肠杆菌密码
子偏好性影响的研究并不多见。
2.2 mRNA的稳定性与翻译起始
mRNA 的降解是影响膜蛋白基因高效表达的重
要因素。mRNA 的降解与其自身的碱基序列有很大
关联，一个碱基的突变就可能会使 mRNA 稳定性出
现巨大变化，从而导致膜蛋白的最终表达水平大受
影响。以大肠杆菌 OmpA 蛋白和人体多巴胺 D2 型
受体（DRD2）为例 ：当 OmpA 基因 5 端的密码子
被一些不常用的同义密码子取代后，mRNA 的半衰
期缩短了 4 倍，mRNA 水平和最终的蛋白产量都降
低了 10 倍［32］；至于 DRD2 基因，当第 957 位的胞
嘧 啶（C） 被 胸 腺 嘧 啶（T） 替 换 后，mRNA 的 半
衰期从 8 h 降至 4 h，蛋白产量也下降到了原来的
50%，但是同时将第 1101 碱基突变为腺嘌呤（A）
后，mRNA 稳定性和蛋白表达水平又恢复正常［33］。
mRNA 的降解还可能与其第 20-45 个密码子的裸露
有密切关联，在这一段序列前面或后面加一些低翻
译速率的密码子分别可以缩短或延长 mRNA 的半衰
期，原因可能是核糖体在这段序列上更长时间的停
留和包裹可以保护 mRNA 免于被降解［34］。
在 mRNA 起始密码子附近的 -4-+37 位碱基处
避免形成稳定的二级结构对于蛋白的表达非常重要，
这正是核糖体在翻译起始的部位［35］。mRNA 5 端
临近翻译起始部位疏松的二级结构有利于核糖体亚
基的快速识别和重组，提高翻译的起始效率［7］。高
效的翻译起始和延伸，一方面加快了肽链的合成速
率 ；另一方面，由于核糖体迅速的结合在一定程度
上抑制了 mRNA 自身的降解，使转录出来的 mRNA
保持一个比较高的水平［7］。在有些条件下，mRNA
5 端二级结构的优化甚至比全局的密码子优化对膜
蛋白基因的表达影响更为明显［17，20，35］。如对人体
神经细胞突触回蛋白和大肠杆菌 FtsH 蛋白的研究发
现，膜蛋白的表达水平与其 N 端氨基酸密码子选择
不同而导致的 mRNA 5 端的二级结构变化有着很大
关系［20，36］。事实上，很多高表达基因为了使 mRNA
5 端保持有利于翻译起始的二级结构而不得不在密
码子偏好性上作出妥协［37-39］。因此，运用聚合酶链
式反应（PCR）等技术对基因序列的 5 端进行改造
或者添加一些碱基降低翻译起始部位 mRNA 二级结
构的稳定性，是使膜蛋白获得高效表达的一个有效
手段［17，40］。
2.3 mRNA与核糖体行为
膜蛋白基因序列的优化还需要考虑到核糖体诱
捕和类似核糖体结合序列（SD 序列）等影响核糖体
行为的因素。在蛋白质的合成过程中，一条 mRNA
通常会同时结合多个核糖体，每个核糖体合成一条
肽链，它们的相对独立并同时进行使得翻译过程十
分高效。但是，mRNA 上密码子的使用不仅影响到
核糖体在某个位置的移动速率，还会关系到相邻核
糖体之间的距离，在整体翻译速率非常快的情况下，
后面有一个翻译速率慢的密码子有可能引起核糖体
的碰撞和形成队列，这通常会使这些核糖体被诱捕
而长时间不能移动，从而影响了翻译的效率，因此
在序列优化时需要防止翻译过程中核糖体的碰撞和
队列来节省翻译的时间［41］。另外，膜蛋白基因中类
似核糖体结合序列（SD 序列）的结构也是一个不能
忽略的因素。一方面，类 SD 序列的存在可能会使核
糖体在该位置上被长时间扣押而导致翻译停滞，极
大地降低肽链延伸速率甚至造成蛋白合成提前终止；
另一方面，原核表达是转录和翻译的耦合，终止密
码子附近的类 SD 序列导致核糖体在翻译即将结束
时移动减速，使其不会影响到终止子茎环结构的折
叠从而确保转录的正常终止［31］。
2.4 翻译的效率与膜蛋白的折叠
膜蛋白基因携带的遗传信息还会调节翻译的肽
链的折叠、复合物的形成及膜蛋白的插入等［12，42］。
适时的翻译速率的减缓对膜蛋白在大肠杆菌中正确
表达是非常重要的，慢的翻译速率可以让超量表达
的膜蛋白更好地折叠、定位并保真［43］。如对蛋白家
族数据库（Pfam）的统计显示，许多与膜相关的蛋
白质都有很长的稀有密码子簇，这些稀有密码子簇
大都在前 130 个氨基酸残基内，并且在靠近 N 端的
2015,31(12) 53王世山等：膜蛋白原核表达的基因序列优化
区域稀有密码子的数量会大幅攀升。在基因表达的
过程中，稀有密码子簇引起的翻译的减速或暂停，
有利于初始合成的肽链在核糖体通道里形成所需的
α 螺旋结构。离开核糖体后，α 螺旋上的疏水残基与
水环境相互作用，快速地重新排列形成更为复杂的
结构，所以存在于膜蛋白结构域之间的稀有密码子
簇的一个作用可能就是通过形成多个翻译暂停，对
蛋白结构进行分步包装，使蛋白正确地折叠和定位。
因此在基因编码区设置适当的翻译暂停位点对于膜
蛋白结构的正确形成，如跨膜螺旋的折叠，有很好
的促进作用［42］。
3 总结
我们从稀有密码子的优化、mRNA 的稳定性与
翻译起始、mRNA 与核糖体行为、翻译的效率与膜
蛋白的折叠 4 个方面对基因序列上影响膜蛋白表达
的优化因素进行了总结，并列举了一些优化方案和
成功实例。基因的序列优化工作是膜蛋白原核表达
的重要开端，在实际操作过程中多个因素之间经常
会相互影响，必要的时候可能需要作出一定程度的
妥协才能得到最优的序列，一段合适的基因序列可
以让后面的表达纯化工作达到事半功倍的效果。目
前对于 mRNA 降解、翻译暂停和膜蛋白折叠等过程
的认识还不够充分，所以膜蛋白基因序列的优化也
还需要未来更多工作的支持。膜蛋白的表达本身是
多个过程的综合考量，除了本文所讲的基因序列优
化外，后续的表达载体、融合标签以及宿主菌的选择、
培养条件的优化等也是膜蛋白高效表达时应该认真
考虑和设计的工作。
参 考 文 献
［1］Yildirim MA, Goh KI, Cusick ME, et al. Drug-target network［J］.
Nature Biotechnology, 2007, 25（10）：1119-1126.
［2］Cross TA, Ekanayake V, Paulino J, et al. Solid state NMR ：The
essential technology for helical membrane protein structural
characterization［J］. Journal of Magnetic Resonance, 2014, 239 ：
100-109.
［3］Chen YK, Zhang ZF, Tang XQ, et al. Conformation and topology of
diacylglycerol kinase in E. coli membranes revealed by solid-state
NMR spectroscopy［J］. Angewandte Chemie-International Edition,
2014, 53（22）：5624-5628.
［4］Su PC, Si W, Baker DL, et al. High-yield membrane protein
expression from E. coli using an engineered outer membrane protein
F fusion［J］. Protein Science, 2013, 22（4）：434-443.
［5］Dong GF, Wang CZ, Wu YH, et al. Tat peptide-mediated soluble
expression of the membrane protein LSECtin-CRD in Escherichia
coli［J］. PLoS One, 2013, 8（12）：e83579.
［6］Zuo X, Lie S, Hall J, et al. Enhanced expression and purification
of membrane proteins by SUMO fusion in Escherichia coli［J］.
Journal of Structural and Functional Genomics, 2005, 6（2-3）：
103-111.
［7］Guimaraes JC, Rocha M, Arkin AP. Transcript level and sequence
determinants of protein abundance and noise in Escherichia
coli［J］. Nucleic Acids Research, 2014, 42（8）：4791-4799.
［8］Prilusky J, Bibi E. Studying membrane proteins through the eyes
of the genetic code revealed a strong uracil bias in their coding
mRNAs［J］. Proceedings of the National Academy of Sciences of
the United States of America, 2009, 106（16）：6662-6666.
［9］Norholm MHH, Light S, Virkki MTI, et al. Manipulating the genetic
code for membrane protein production ：what have we learnt so
far?［J］. Biochimica Et Biophysica Acta-Biomembranes, 2012,
1818（4）：1091-1096.
［10］Greenbaum D, Jansen R, Gerstein M. Analysis of mRNA expression
and protein abundance data ：an approach for the comparison of
the enrichment of features in the cellular population of proteins and
transcripts［J］. Bioinformatics, 2002, 18（4）：585-596.
［11］Jansen R, Gerstein M. Analysis of the yeast transcriptome with
structural and functional categories ：characterizing highly
expressed proteins［J］. Nucleic Acids Research, 2000, 28（6）：
1481-1488.
［12］Kimchi-Sarfaty C, Oh JM, Kim IW, et al. A “silent” polymorphism
in the MDR1 gene changes substrate specificity［J］. Science,
2007, 315（5811）：525-528.
［13］Kepes F. The“+70 pause”：Hypothesis of a translational control
of membrane protein assembly［J］. Journal of Molecular Biology,
1996, 262（2）：77-86.
［14］Dessen P, Kepes F. The PAUSE software for analysis of
translational control over protein targeting：application to E-nidulans
membrane proteins［J］. Gene, 2000, 244（1-2）：89-96.
［15］Katzen F, Peterson TC, Kudlicki W. Membrane protein expression ：
生物技术通报 Biotechnology Bulletin 2015,Vol.31,No.1254
no cells required［J］. Trends in Biotechnology, 2009, 27（8）：
455-460.
［16］Wang Q, Mei C, Zhen HH, et al. Codon preference optimization
increases prokaryotic cystatin C expression［J］. Journal of Biom-
edicine and Biotechnology, 2012, doi ：org/10. 1155/2012/732017.
［17］Norholm MHH, Toddo S, Virkki MTI, et al. Improved production
of membrane proteins in Escherichia coli by selective codon
substitutions［J］. Febs Letters, 2013, 587（15）：2352-2358.
［18］Grosjean H, Fiers W. Preferential codon usage in prokaryotic genes-
the optimal codon anticodon interaction energy and the selective co-
don usage in efficiently expressed genes［J］. Gene, 1982, 18（3）：
199-209.
［19］Gouy M, Gautier C. Codon usage in bacteria-correlation with gene
expressivity［J］. Nucleic Acids Research, 1982, 10（22）：
7055-7074.
［20］Löw C, Jegerschöld C, Kovermann M, et al. Optimisation of over-
expression in E. coli and biophysical characterisation of human
membrane protein synaptogyrin 1［J］. PLoS One, 2012, 7（6）：
e38244.
［21］Hardt B, Volker C, Mundt S, et al. Human endo-alpha 1,
2-mannosidase is a golgi-resident type II membrane protein［J］.
Biochimie, 2005, 87（2）：169-179.
［22］Hassan KA, Xu ZQ, Watkins RE, et al. Optimized production and
analysis of the staphylococcal multidrug efflux protein QacA［J］.
Protein Expression and Purification, 2009, 64（2）：118-124.
［23］Slimko EM, Lester HA. Codon optimization of Caenorhabditis
elegans GluCl ion channel genes for mammalian cells dramatically
improves expression levels［J］. Journal of Neuroscience Methods,
2003, 124（1）：75-81.
［24］Sohl CD, Guengerich FP. Kinetic Analysis of the three-step steroid
aromatase reaction of human cytochrome P450 19A1［J］. Journal
of Biological Chemistry, 2010, 285（23）：17734-17743.
［25］Baneres JL, Martin A, Hullot P, et al. Structure-based analysis of
GPCR function ：conformational adaptation of both agonist and
receptor upon leukotriene B-4 binding to recombinant BLT1［J］.
Journal of Molecular Biology, 2003, 329（4）：801-814.
［26］Calderone TL, Stevens RD, Oas TG. High-level misincorporation of
lysine for arginine at AGA codons in a fusion protein expressed in
Escherichia coli［J］. Journal of Molecular Biology, 1996, 262（4）：
407-412.
［27］Gurvich OL, Baranov PV, Gesteland RF, et al. Expression levels
influence ribosomal frameshifting at the tandem rare arginine
codons AGG_AGG and AGA_AGA in Escherichia coli［J］.
Journal of Bacteriology, 2005, 187（12）：4023-4032.
［28］McNulty DE, Claffee BA, Huddleston MJ, et al. Mistranslational
errors associated with the rare arginine codon CGG in Escherichia
coli［J］. Protein Expression and Purification, 2003, 27（2）：
365-374.
［29］Vyas VV, Esposito D, Sumpter TL, et al. Clinical manufacturing
of recombinant human interleukin 15. I. Production cell line
development and protein expression in E. coli with stop codon
optimization［J］. Biotechnology Progress, 2012, 28（2）：497-
507.
［30］Whittaker MM, Whittaker JW. Expression and purification of
recombinant Saccharomyces cerevisiae mitochondrial carrier protein
YGR257Cp（Mtm1p）［J］. Protein Expression and Purification,
2014, 93 ：77-86.
［31］Li GW, Oh E, Weissman JS. The anti-Shine-Dalgarno sequence
drives translational pausing and codon choice in bacteria［J］.
Nature, 2012, 484（7395）：538-541.
［32］Deana A, Ehrlich R, Reiss C. Silent mutations in the Escherichia
coli ompA leader peptide region strongly affect transcription and
translation in vivo［J］. Nucleic Acids Research, 1998, 26（20）：
4778-4782.
［33］Duan JB, Wainwright MS, Comeron JM, et al. Synonymous
mutations in the human dopamine receptor D2（DRD2）affect
mRNA stability and synthesis of the receptor［J］. Human
Molecular Genetics, 2003, 12（3）：205-216.
［34］Pedersen M, Nissen S, Mitarai N, et al. The functional half-life
of an mRNA depends on the ribosome spacing in an early coding
region［J］. Journal of Molecular Biology, 2011, 407（1）：35-
44.
［35］Kudla G, Murray AW, Tollervey D, et al. Coding-Sequence
determinants of gene expression in Escherichia coli［J］. Science,
2009, 324（5924）：255-258.
［36］Makino S, Qu JN, Uemori K, et al. A silent mutation in the ftsH
gene of Escherichia coli that affects FtsH protein production and
colicin tolerance［J］. Molecular & General Genetics, 1997, 254
（5）：578-583.
［37］Hockenberry AJ, Sirer MI, Amaral LAN, et al. Quantifying
2015,31(12) 55王世山等：膜蛋白原核表达的基因序列优化
position-dependent codon usage bias［J］. Molecular Biology and
Evolution, 2014, 31（7）：1880-1893.
［38］Bentele K, Saffert P, Rauscher R, et al. Efficient translation
initiation dictates codon usage at gene start［J］. Molecular
Systems Biology, 2013, 9 ：675.
［39］ Goodman DB, Church GM, Kosuri S. Causes and effects of
N-terminal codon bias in bacterial genes［J］. Science, 2013, 342
（6157）：475-479.
［40］ Tang ZM, Salamanca-Pinzon SG, Wu ZL, et al. Human cytochrome
P450 4F11 ：Heterologous expression in bacteria, purification,
and characterization of catalytic function［J］. Archives of
Biochemistry and Biophysics, 2010, 494（1）：86-93.
［41］ Mitarai N, Pedersen S. Control of ribosome traffic by position-
dependent choice of synonymous codons［J］. Physical Biology,
2013, 10（5）：056011.
［42］Chartier M, Gaudreault F, Najmanovich R. Large-scale analysis
of conserved rare codon clusters suggests an involvement in co-
translational molecular recognition events［J］. Bioinformatics,
2012, 28（11）：1438-1445.
［43］Fluman N, Navon S, Bibi E, et al. mRNA-programmed translation
pauses in the targeting of E. coli membrane proteins. ［J］. Elife,
2014, 3. doi ：10. 7554/eLife. 03440.
（责任编辑 狄艳红）