仓库源文站点原文


layout: post title: Amber构建糖蛋白的方法 categories:


【按】以下翻译自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文件

  1. 删除PDB文件中所有的HETATM原子. 这些原子通常包括结合的配体, 非结晶的水分子和非配位的金属离子. 如果存在氢原子也要全部删除.
  2. 通常, 我们要检查蛋白质文件以确保其中没有重复原子. 这可以通过使用LEaP加载蛋白质后, 检查对应的警告来完成. 对我们要处理的这个特殊例子, 第119位残基(HIS)包含了重复的侧链原子. 删除其中任意一组重复原子即可.
  3. 检查PDB文件的开头部分看是否存在二硫键(SSBOND). 3RN3蛋白中有四个二硫键, 存在于下面几对半胱氨酸残基之间: 26-84, 40-95, 58-110和65-72. 将这8个半胱氨酸残基的名称由CYS改为CYX.
  4. 目前, 可以将糖分子与丝氨酸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节中的库文件列表.
  5. 准备一个包含蛋白质和糖分子的PDB文件, 其中糖分子相对蛋白质表面的取向适当. 进行这一步操作有几种方法, 包括:

13.6.3.2 示例: 将枝状糖连接到3RN3(N-连接的糖基化)

在这个例子中, 我们假定上面章节中已经生成的糖分子(branch.pdb)在蛋白质文件中相对于ASN 34的取向已经适当, 并且复合物已保存为新的PDB文件(例如, 3rn3_nlink.pdb). 最后一个氨基酸残基应该为VAL 124, 糖分子的名称应该为4YB 125, 4YB 126, VMB 127, OMA 128OMA 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>