layout: post
title: Amber构建糖蛋白的方法
categories:
- 2017年04月06日 18:29:59 翻译: 房筱; 校对: 李继存
【按】以下翻译自Amber 16手册.
13.6.3 使用LEaP构建糖蛋白
本节给出的LEaP命令假定你已经有一个包含糖分子和蛋白质的PDB文件, 并且二者之间的相对构象适当. 即便是为了将最简单的糖与最简单的蛋白质相连接, 也需要你对LEaP命令有透彻的理解, 但这超出了本节的目的. 下面会给出与生成相应PDB文件有关的几个选项(条目5a-5c).
本例中使用的蛋白质是牛核糖核酸酶A(PDBID: 3RN3). 我们使用上面章节的第二个示例中创建的枝状寡糖, 并将其(N-连接)连接到ASN 34以产生核糖核酸酶B.
13.6.3.1 在LEaP中设置糖基化蛋白质的PDB文件
- 删除PDB文件中所有的
HETATM
原子. 这些原子通常包括结合的配体, 非结晶的水分子和非配位的金属离子. 如果存在氢原子也要全部删除.
- 通常, 我们要检查蛋白质文件以确保其中没有重复原子. 这可以通过使用LEaP加载蛋白质后, 检查对应的警告来完成. 对我们要处理的这个特殊例子, 第119位残基(HIS)包含了重复的侧链原子. 删除其中任意一组重复原子即可.
- 检查PDB文件的开头部分看是否存在二硫键(
SSBOND
). 3RN3蛋白中有四个二硫键, 存在于下面几对半胱氨酸残基之间: 26-84, 40-95, 58-110和65-72. 将这8个半胱氨酸残基的名称由CYS
改为CYX
.
- 目前, 可以将糖分子与丝氨酸SER, 苏氨酸THR, 羟脯氨酸HYP和天冬酰胺ASN连接起来. 必须在蛋白质的PDB文件中手动重命名氨基酸后, 才能将其加载到LEaP中. 修饰的残基名称可取为
OLS
(O连接的SER), OLT
(O连接的THR), OLP
(O连接的羟脯氨酸HYP)和NLN
(N连接的ASN). 当source leaprc.GLYCAM_06j-1
时, 修改后的包含氨基酸残基的库文件会自动加载. 更多信息, 请参阅3.3节中的库文件列表.
- 准备一个包含蛋白质和糖分子的PDB文件, 其中糖分子相对蛋白质表面的取向适当. 进行这一步操作有几种方法, 包括:
- 实验用的PDB文件中通常存在一个或多个糖残基. 在这种情况下, 合理的方法是将GLYCAM生成的待连接糖分子, 与实验用PDB文件中的糖分子进行叠合, 然后保存更改后的坐标. 如果你使用这种方法, 记得从PDB文件中删除实验中的糖分子! 还必须保证PDB文件中每个碳水化合物残基以
TER
与其他残基分开. 还要记得删除糖分子末端的OH或OMe残基. 另外, 实验中的糖分子可以保留在PDB文件中, 只要它们的名称遵循GLYCAM 3字符编码, 并且其原子在PDB文件中的名称和顺序与GLYCAM的标准符合. 这些操作很麻烦, 但确实能成功达到目的. 再次, 请务必确保在蛋白质和碳水化合物, 以及碳水化合物残基和碳水化合物残基之间插入TER
.
- 使用分子建模程序将GLYCAM生成的糖分子与蛋白质对齐, 并将坐标保存在单个文件中. 记住删除糖分子末端的OH或OMe.
- 使用在线的糖蛋白构建工具http://www.glycam.org. 该工具允许用户上传蛋白质的坐标, 构建糖分子(或从库中选择), 并将其连接到蛋白质. 然后可以下载所需的AMBER文件. 该网站也能用于处理仅含蛋白质的文件, 只要上传到糖蛋白构建工具就可以了, 很方便.
13.6.3.2 示例: 将枝状糖连接到3RN3(N-连接的糖基化)
在这个例子中, 我们假定上面章节中已经生成的糖分子(branch.pdb
)在蛋白质文件中相对于ASN 34的取向已经适当, 并且复合物已保存为新的PDB文件(例如, 3rn3_nlink.pdb
). 最后一个氨基酸残基应该为VAL 124
, 糖分子的名称应该为4YB 125
, 4YB 126
, VMB 127
, OMA 128
和OMA 129
.
记得将残基ASN 34
的名称从ASN
更改为NLN
. 对于糖分子的结构, 请确保PDB文件中的每个糖残基都以TER
与其他残基分开. 不要在这里使用seq
命令, 并且所有的连接信息(糖分子之间, 以及糖与蛋白质之间的)都将单独指定.
在xleap
中输入以下命令(如果不需要图形界面, 则可以使用tleap
). 或者将命令保存到文件中然后source
.
<div class="highlight" style="background: #f8f8f8"><pre style="line-height: 125%"><span></span><span style="color: #AA22FF">source</span><span style="color: #bbbbbb"> </span>leaprc.GLYCAM_06j-1<span style="color: #bbbbbb"> </span><span style="color: #008800; font-style: italic">#</span><span style="color: #bbbbbb"> </span><span style="color: #008800; font-style: italic">加载</span><span style="color: #bbbbbb"> </span><span style="color: #008800; font-style: italic">ff14SB</span><span style="color: #bbbbbb"> </span><span style="color: #008800; font-style: italic">的</span><span style="color: #bbbbbb"> </span><span style="color: #008800; font-style: italic">GLYCAM-06</span><span style="color: #bbbbbb"> </span><span style="color: #008800; font-style: italic">leaprc文件</span><span style="color: #bbbbbb"> </span>
<span style="color: #AA22FF">source</span><span style="color: #bbbbbb"> </span>leaprc.protein.ff14SB<span style="color: #bbbbbb"> </span><span style="color: #008800; font-style: italic">#</span><span style="color: #bbbbbb"> </span><span style="color: #008800; font-style: italic">加载蛋白力场</span><span style="color: #bbbbbb"></span>
<span style="color: #B8860B">glyprot</span><span style="color: #bbbbbb"> </span><span style="color: #666666">=</span><span style="color: #bbbbbb"> </span>loadpdb<span style="color: #bbbbbb"> </span>3rn3_nlink.pdb<span style="color: #bbbbbb"> </span><span style="color: #008800; font-style: italic">#</span><span style="color: #bbbbbb"> </span><span style="color: #008800; font-style: italic">加载蛋白和糖的pdb文件</span><span style="color: #bbbbbb"></span>
bond<span style="color: #bbbbbb"> </span>glyprot.125.O4<span style="color: #bbbbbb"> </span>glyprot.126.C1<span style="color: #bbbbbb"> </span><span style="color: #008800; font-style: italic">#</span><span style="color: #bbbbbb"> </span><span style="color: #008800; font-style: italic">设定糖之间的键</span><span style="color: #bbbbbb"></span>
bond<span style="color: #bbbbbb"> </span>glyprot.126.O4<span style="color: #bbbbbb"> </span>glyprot.127.C1<span style="color: #bbbbbb"></span>
bond<span style="color: #bbbbbb"> </span>glyprot.127.O6<span style="color: #bbbbbb"> </span>glyprot.128.C1<span style="color: #bbbbbb"></span>
bond<span style="color: #bbbbbb"> </span>glyprot.127.O3<span style="color: #bbbbbb"> </span>glyprot.129.C1<span style="color: #bbbbbb"></span>
bond<span style="color: #bbbbbb"> </span>glyprot.34.SG<span style="color: #bbbbbb"> </span>glyprot.125.C1<span style="color: #bbbbbb"> </span><span style="color: #008800; font-style: italic">#</span><span style="color: #bbbbbb"> </span><span style="color: #008800; font-style: italic">设定糖和蛋白之间的键</span><span style="color: #bbbbbb"></span>
bond<span style="color: #bbbbbb"> </span>glyprot.26.SG<span style="color: #bbbbbb"> </span>glyprot.84.SG<span style="color: #bbbbbb"> </span><span style="color: #008800; font-style: italic">#</span><span style="color: #bbbbbb"> </span><span style="color: #008800; font-style: italic">设定二硫键</span><span style="color: #bbbbbb"></span>
bond<span style="color: #bbbbbb"> </span>glyprot.40.SG<span style="color: #bbbbbb"> </span>glyprot.95.SG<span style="color: #bbbbbb"></span>
bond<span style="color: #bbbbbb"> </span>glyprot.58.SG<span style="color: #bbbbbb"> </span>glyprot.110.SG<span style="color: #bbbbbb"></span>
bond<span style="color: #bbbbbb"> </span>glyprot.65.SG<span style="color: #bbbbbb"> </span>glyprot.72.SG<span style="color: #bbbbbb"></span>
addions<span style="color: #bbbbbb"> </span>glyprot<span style="color: #bbbbbb"> </span>CL<span style="color: #bbbbbb"> </span><span style="color: #666666">0</span><span style="color: #bbbbbb"> </span><span style="color: #008800; font-style: italic">#</span><span style="color: #bbbbbb"> </span><span style="color: #008800; font-style: italic">添加电荷抗衡离子,</span><span style="color: #bbbbbb"> </span><span style="color: #008800; font-style: italic">使体系电中性</span><span style="color: #bbbbbb"></span>
solvateBox<span style="color: #bbbbbb"> </span>glyprot<span style="color: #bbbbbb"> </span>TIP3P<span style="color: #bbbbbb"> </span>BOX<span style="color: #bbbbbb"> </span><span style="color: #666666">8</span><span style="color: #bbbbbb"> </span><span style="color: #008800; font-style: italic">#</span><span style="color: #bbbbbb"> </span><span style="color: #008800; font-style: italic">添加溶剂</span><span style="color: #bbbbbb"></span>
savepdb<span style="color: #bbbbbb"> </span>glyprot<span style="color: #bbbbbb"> </span>3nr3_glycan.pdb<span style="color: #bbbbbb"> </span><span style="color: #008800; font-style: italic">#</span><span style="color: #bbbbbb"> </span><span style="color: #008800; font-style: italic">保存pdb文件,</span><span style="color: #bbbbbb"> </span><span style="color: #008800; font-style: italic">top,</span><span style="color: #bbbbbb"> </span><span style="color: #008800; font-style: italic">crd文件</span><span style="color: #bbbbbb"></span>
saveamberparm<span style="color: #bbbbbb"> </span>glyprot<span style="color: #bbbbbb"> </span>3nr3_glycan.top<span style="color: #bbbbbb"> </span>3nr3_glycan.crd<span style="color: #bbbbbb"></span>
quit<span style="color: #bbbbbb"> </span><span style="color: #008800; font-style: italic">#</span><span style="color: #bbbbbb"> </span><span style="color: #008800; font-style: italic">退出</span><span style="color: #bbbbbb"></span>
</pre></div>