Discovery Studio 讲义中文

附录：准备一个PDB文件作为同源性

建模模板

同源建模的基本原理是，你映射的一个未知的蛋白质序列一种已知蛋白质的结构。因此，如果你没有已知的蛋白质，或模板，你将无法建立模型。模板的共同来源是蛋白质在结构生物信息学研究的实验室数据银行（目标）。该网站RCSB是HTTP：/ /www.rcsb。org /。

蛋白质数据库（PDB）可能是世界上领先的公共源三维生物分子数据（1）。截至七月2006，超过37000项可在PDB。每个月都有更多的人加入。X射线衍射仪和其它固态技术占大多数的结构。然而，超过5500的核磁共振结构还可用。这些沉积的结构包括蛋白质，肽，核酸，碳水化合物，这些分子的配合物。

在发现工作室环境中工作的这些分子是一个关键的过程给你的建模工作。这个练习如何准备一个PDB文件作为一个模板同源建模项目。你将在课程中学习如何：

?加载PDB文件直接从蛋白质数据银行，

?生产和检验蛋白质报告，

?清除晶体单元电池，

?分裂分子，

?删除不需要的组件，和保存已完成的文件。

1、开始发现工作室

我们必须有发现工作室开始行使。

推出发现工作室客户端，如果它没有运行。

如果发现工作室已经在运行，从窗口菜单，选择关闭所有命令

如果提示保存任何分子或数据，请选择“不”。

2、加载PDB文件

现在，我们将直接从目标通过Discovery Studio界面加载PDB文件。注意：文件|打开网址…命令只会获得文件通过网络连接一个数据库服务器。如果连接不可能返回错误。检查你的导师或系统管理员，如果一个连接是可从您的车间位置。

从文件下拉菜单，选择命令打开网址。在对话框中，网址应该是指目标网站。更换与1t64 URL的最后四个字。

点击打开按钮。

如果一个PDB不能连接，所需的文件中的数据文件是可用的目录

1t64_original.pdb。

在3D窗口中显示的结构是两人HDAC8蛋白的晶体结构在复杂的抑制剂，684个晶体的水，和几个离子。抑制剂是曲古抑菌素A，一个很好的研究组蛋白去乙酰化酶抑制剂（2）。

从“视图”菜单中选择“层次”命令。

注意层次视图中组件的分布。我们将重新安排组件为了方便工作的结构。

3、检查蛋白质报告

在深入了解蛋白质，我们应该花一点时间来看看有什么。一种快速的方法完成这项任务是使用蛋白质报告。

在蛋白质报告和公用事业工具面板，去蛋白质信息部分，然后点击蛋白质报告命令。

这将一个文本文档总结关键信息的HTML窗口打开1t64蛋白结构。正如我们有多个分子在视图中，该报告从一个分子列表。

滚动下来的报告，你会注意到，首先有2条链，一个和B。

此外，还有其他的构象和缺失的残基，这两个链。还有更多链中的问题的残留量比在乙链。进一步向下滚动，你会看到有几个配体和离子包括在结构和晶体的水。

根据数据的报告，我们可以得出结论，乙链是一个更好的模板比连锁。但我们要清楚我们几个非原子以及。

4、删除单元电池

晶体单元单元在分子力学计算中的意义不大。有时它只是更容易删除单元电池。

从结构的菜单，选择水晶细胞pullright菜单。选择删除单元格命令。

请注意，在层次视图中的更改为单元单元格信息被删除。

扩大对象1t64验证蛋白质和水分子仍然存在。

5、拆分分子系统的组成部分

这个文件是由多个部分组成的。我们可以通过层次结构查看组件其他练习窗口。然而，我们实际上可以分开的组件的PDB文件更快的分割命令。

访问工具浏览器。然后访问蛋白质报告和公用事业工具。

找分裂工具组。列出了三个选项，并将在一个瞬间描述。

单击所有命令。

注意在层次视图中的更改。PDB文件的四个组成部分已

隔离到单独的对象

四的成分是蛋白质（1t64_a和1t64_b），配体（1t64_nonprotein），和结晶水（1t64_water）。拆分命令的三个选项允许在如何处理对象的操作中有一定的灵

活性。

所有

打出的层次结构视图中所有链和非蛋白结构为独立的对象，并列出每个氨基酸序列为序列单独序列。

蛋白

打出所有的蛋白链结构为独立的对象层次结构中的每个视图和列表的氨基酸序列为序列单独序列。

非蛋白

打出配体，如为单独的对象的水域在层次视图其他非蛋白结构和列表的任何非蛋白多肽序列在序列视图分离的序列。

6、除去水

对于这次演习，我们将删除晶体的水分子。有几种方法水可以被移除。在这种情况下，当我们拆分系统和开发新的层次，我们将使用层次视图。

在层次视图中，选择进入1t64_water。输入应选择在两个层次视图和三维窗口。

点击键盘上的删除键。

从两个层次视图和三维窗口中的水分子被删除。

7、除去配位体和离子

无论是配位体，也没有离子所需的我们的同源性建模实验。所以我们可以删除它们。

在层次视图中，选择进入1t64_nonprotein。

点击键盘上的删除键。

从两个层次视图和三维窗口中除去的配位体和离子分子。

8、删除一个蛋白质单位

双蛋白链在本单位细胞中显示。然而，只有一个将需要作为模板。一个人的意志，将取决于个体链的质量。在这种情况下，我们在蛋白质报告中看到，有更多的残留物与问题（缺少的原子和替代的构象），因此，我们将保留只有乙链。

在层次视图中，选择进入1t64_a。

点击键盘上的删除键。

只有乙链仍然可见。

9、清洁蛋白质

发现工作室可以自动完成所需的任务，以适当的准备一个蛋白质能量计算。在进行清洁操作之前，我们可以指定操作通过偏好对话框进行。

从“编辑”菜单中选择“首选项”命令。

在“首选项”对话框中，展开“蛋白质公用事业”页。选择干净

的蛋白质页面。

在这一点上，我们有几个可供选择的选项。这些选项将共同解决可能在PDB文件存在的问题。例如，X射线晶体结构通常会没有氢原子，所以必须添加。此外，链端必须设置正确的化学物质和残留的残留的原子必须放置。描述的选项下面。

对的问题非标准名称

检查原子名称是否符

合标准如果需要的话，他

们的名字和改正。

障碍（保留一个设置）

检查无序原子和保留

第一套。

不完整的残基

增加了缺失的侧链原

子的氨基酸。这操作不会

填充缺少的循环区域。

所需的PH值允许的可离子化的氨基的质子化状态酸和末端被控制使用标准的pKa值。质子化状态之后的任何修改的氢调通过下列选项总站。

氢如果修改氢选检

查，氢会增加按需要。

所有的氢原子：添加所有氢原子。

只有那些极性氢：氢原子可能会参与氢键将加入。

无氢：没有添加氢原子。

总站

如果修改总站选项被选中，然后总站将被添加或删除，如。

固定连接和键订单

确保氨基酸有正确的键秩序。此选项将不会影响非标准氨基酸或配体。

现在我们将设置干净的蛋白质偏好。

打开喜好如图1所示（在下一页）。单击“确定”按钮偏好对话框并返回到三维窗口。

现在，我们实际上可以把蛋白质清洗干净。

在层次视图中，选择进入1t64_A

图1：清洁蛋白偏好板

访问工具浏览器。然后访问蛋白质报告和公用事业工具。

单击“干净”命令。

几分钟后，蛋白质结构，是目前在三维窗口的变化。最值得注意的是氢原子现在加入到结构中。其他操作也已执行也如完成一个缺失的残基和校正原子名。

10、拯救新的复杂

我们将保存这种蛋白质用于以后的使用。

确保三维窗口是活动的。

点击“适合屏幕”按钮。

从文件下拉菜单，选择“另存为…”命令。

对象名称输入名称1t64。确保类型的文件：参数设置为查看（*。MSV）。

单击“保存”按钮。

11、完成课程

关上发现工作室的所有窗口。

从“窗口”菜单中选择“关闭”命令。

或者，你可以简单地退出发现工作室。

如果要求保存任何数据，请选择No。

有了这个教训，我们采取了由各种分子实体PDB文件。我们有打破这些分子组成和重组他们的需要。然后我们保存结构。使用发现工作室的同源建模，我们需要尽我们的严格进行能量计算时。因此，它是没有必要的类型的分子，以只要我们不做一个能量计算。同源建模只需要一个有效的蛋白质结构为模板，这就是我们在这里所做的。

References:

1. Berman HM,Westbrook J, Feng Z, Gilliland G, Bhat TN, et al. 2000. The Protein Data Bank. Nucleic Acids Res 28: 235

2. Yoshida M, Kijima M, Akita M, Beppu T. 1990. Potent and specific inhibition of mammalian histone deacetylase both in vivo and in vitro by trichostatin A. J Biol Chem 265: 17174 265：17174

手在练习：序列搜索

课程1：爆炸搜索

在ε-组蛋白赖氨酸的氨基乙酰化具体发挥关键作转录和基因调控（11）。乙酰基的去除是通过，组蛋白去乙酰化酶（HDAC），deactylate 核小体组蛋白。最近，抑制剂组蛋白去乙酰化酶已确定诱导生长抑制、分化和凋亡癌细胞死亡（9），因此具有一定的药物潜力。然而，有HDACs的几类；和，在每个基因的表达有显著作用（7）。开发类HDAC 抑制剂是一种策略被用于抗癌的发展药物。

为了设计这些特异性抑制剂，对不同的目标酶的结构将所需。然而，而不是等待晶体结构，工作已经开始创造HDACs同源模型（15）。在王和同事的研究中，同源模型建立了四个I类组蛋白去乙酰化酶抑制剂的结合模式探讨。我们在这个工作坊中使用这项工作作为练习的基础。

在开始任何同源建模实验之前，你必须要有一个同源蛋白建立您的模型。发现工作室提供访问爆炸和PSI-BLAST搜索蛋白质序列数据库。爆炸搜索被用来识别目标序列的近同源序列（1）。重复BLAST

搜索用于识别更遥远的同系物（2）。

在本教程中，您将运行一个序列的BLAST搜索人类组蛋白去乙酰化酶（HDAC1）在PDB数据库中找到同源序列。我们会用点击来帮助确定一个合适的模板建立HDAC1同源模型。

完成这门课之后，你会很熟悉：

?从外部来源阅读到发现工作室的序列，

?使用爆炸协议进行搜索，并

?加载和分析所获得的点击量。

1、获得目标序列

注意：如果你的互联网连接和浏览器可在您的指导，请检查

训练场地。如果没有，一个准备好的文件可供行使。

我们将开始获得我们的目标或未知结构的序列。

在你的电脑上启动您的互联网浏览器。

输入网址https://www.360docs.net/doc/353619579.html,/database/去NCBI数据库。

在搜索框中，改变Entrez蛋白质。

然后输入GenBank登录号q13547。点击去进行搜索。

你会看到人类组蛋白去乙酰化酶HDAC1的入口（14）。该项的格式为无法装入发现工作室。我们必须改变格式。

更改显示设置FASTA。

HDAC的序列是现在显示在FASTA格式或皮尔森（12），可阅读发现工作室。FASTA格式是一种常用的简单的ASCII文本文件转移序列。但是，我们需要把这个序列插入一个文件中。

接下来的显示设置下拉菜单标记发送到。更改发送到参数文件。

下一步是依赖于特定的浏览器和设置。因此，指令是相当普遍。

一个对话框应该问你想做的文件。选择保存。如果提示一个文件名，输入sequences.fasta。如果提示一个位置，选择桌面。

保存文件。

该文件将保存在您的桌面上。如果您有麻烦定位或下载该文件，不关心。序列在车间的数据文件是可用的，通常在数据文件桌面上的文件夹。

此外，还有无数的地方，我们可以得到的序列。例如，我们可以使用Accelrys GCG（俗称威斯康星州包）和搜索数据库。或者我们可以去其他的数据库，如互联网上可用SWISS-PROT（HTTP：/

/www.expasy。org /运动/）。无论我们得到的序列，它必须在被发现工作室阅读的文件格式。合适的格式的列表是讲座介绍。

2、开始发现工作室

我们必须有发现工作室开始行使。

推出发现工作室客户端，如果它没有运行。

如果发现工作室已经在运行，从窗口菜单中，选择关闭所有命令。如果提示保存任何分子或数据，请选择“不”。

3、负荷的HDAC1序列

我们将负载的人类蛋白deacteylase序列文件。

从文件下拉菜单，选择命令打开。

在类型的文件：下拉，选择序列格式过滤掉非顺序文件。

单击桌面图标以跳转到桌面。选择刚才下载的序列文件sequences.fasta或任何你已经命名的文件。如果您无法找到该文件，你可以浏览到桌面/文件目录和使用文件hdac1.fasta。

点击打开按钮。

序列窗口将出现我们未知的氨基酸序列。

Discovery Studio给出序列来自NCBI指数名称。我们将它重命名为更容易理解的东西。

在“序列”窗口中，双击“序列”或“序列”的名称。确保整个序列被高亮显示。

右键单击并选择“重命名”顺序。在打开输入的对话框中名称HDAC1。点击确定。取消一切。

3、运行一次爆炸搜索

你可以通过协议访问爆炸程序。

从协议资源管理器中，打开蛋白质建模协议组。

选择爆炸协议。

参数浏览器将显示为一个爆炸搜索可用的设置。

序列的参数，输入序列窗口名称序列：HDAC1。

以类似的方式指定结构，特定的序列被附加到该名称的序列窗口。您只能选择一个序列，在搜索中使用。

确保数据库参数设置为pdb_nr95。

默认的数据库是一个高质量的、非冗余设置来自PDB（3）。其他支持数据库的nrdB NCBI、EBI、PDB nrdb90，和Swiss-Prot中描述演讲。如果你将一个新的格式化数据库添加到数据库目录中，它也是上市。

确保矩阵参数设置为BLOSUM62。

你必须使用BLOSUM选项（6）或（4）PAM矩阵。该BLOSUM62矩阵一个合理的开始在大多数搜索。

确保间隙的调整设置为true。

这个选项让你选择是否使用跳空爆炸或进行计算搜索可以运行。空位BLAST是PSI-BLAST一迭代的等效（2）。

我们将把剩下的选项设置为默认值如下：

缺口成本存在： 11延伸：1

期望值 10

滤波器低复杂度 Ture

输出序列数 250

这四个参数允许您指定搜索将继续进行。它超越了本次研讨会的范围，以解释所有的选项可能。通过爆炸的讨论搜索在Accelrys GCG

车间是可用的。爆炸选项也彻底在最近的一本书（8）中讨论。

单击“运行”按钮。

运行将开始和完成在几秒钟。

4、查看地图视图中的输出

当运行完成时，任务完成对话框出现。点击“确定”按钮。

当运行完成时，出现了“爆炸”窗口（图1）。

使爆炸窗口活动。

图1：用地图视图窗口

地图视图给出了一个图形表示的点击量，显示为一行每序列。酒吧的颜色根据点击的位。着色的一个传说是显示在地图视图的顶部。分数越高，分数越高。

将查询序列显示在视图的顶部，作为一个标尺，用残基索引。视图可以可以放大和缩小使用视图|变换命令，鼠标上的滚动按钮，或控制键和+或键盘上的键。当视图滚动到它的最大尺寸，的标尺被改变来显示单个残基。

你可以将鼠标悬停在地图上的任何一个点击查看和点击的信息

显示。

把鼠标放在第一次成功，等待弹出出现。

弹出的文本框中给出了描述序列，序列号，该开始和结束位置的查询序列和开始和结束的命中序列，和命中得分。

注意，第一次打击是组蛋白去乙酰化酶8的PDB文件1t64代表（13）。

这是让我们想建立一个组蛋白去乙酰化酶模型。在原文件中Wang et al.，结构1t64不可用。

把鼠标移到二次。注意，这个打击是组蛋白去乙酰化酶样蛋白（HDLP）通过PDB文件1c3p代表（5）。

HDLP的发现是用原来的建模研究的蛋白质几乎相同王和同事，

1c3r（5）。所以，我们可以说，我们已经找到了相同的模板在他们的学习中。我们没有得到同样的打击1c3r由于使用的非冗余数据库。把鼠标移到第三个命中。请注意，这个打也是一种组蛋白去乙酰化酶

样蛋白但细菌的同源组蛋白去乙酰化酶水解酶。它是由PDB文件表示（10）1zz1。

这个打击，是1t64，不可在原王等人的时间。研究。

点击鼠标左键可以选择点击。给了一个具体的打击当前选择，按住鼠标移动或控制键。

图2：在表视图中显示的爆炸搜索结果

5、检查表视图中的点击量

数据的另一个视图可以作为表视图（图2）。表格视图显示每行一行每行附加信息。

确保在地图视图中没有选择。

在地图视图的底部是标签。单击“表格”选项卡。现在的表格视图是显示。

数据报告基本上是相同的，在地图视图中发现，除了表格形式。

在表中的点击率可以通过左键点击列标题来排序。什么时候该表被排序，根据顺序在地图视图中的点击顺序改变表格视图。

在表视图中，可以选择左键单击。

这些点击率也在地图视图中选择。反向也是正确的；在地图视图中选择一个命中也选择在表视图中的命中。

查看表视图中的前两个点击。在值列，期望值击中可能被发现。为1t64和1c3p，很低的期望值，表明这些是好的点击。打1zz1还具有低的期望值却不低，作为第一个。通知还宽的差异之间的能量值为第一三安打，其余的桌子。

值是由同一位得分机会得分的期望。低值表示一个真正的打击。然而，单独的值可能是不够的。你也应该考虑对准长度。

在表视图中，检查列的序列长度和对齐长度。注意前两打（1t64和1c3p）是在他们的一个重要部分HDAC1对齐长度。第三打1zz1对齐在更短的长度。最后，在剩下的条目表盖HDAC1更小的部分。

最后一个标签是文本视图。这是原爆输出为ASCII文本文件。6、载荷对准

点击几个方式可以加载到发现工作室。

返回地图视图。

同时按住移位键，左键点击三次点击。

在地图视图区域中，单击右键单击右键。

选择命令加载序列和对齐方式。还注意到五个加载命令可用。

几分钟后，一个新的序列窗口打开的三个点击和原始基因序列（与原来的名字）。

从地图视图或表视图中，五个加载选项可供选择序列，序列，和结构的序列窗口和三维窗口。这些简要说明选项。

负载选择的序列加载的序列的选定的点击进入一个序列窗口

负载选定的结构加载的蛋白质结构的选定的点击进入一个三

维窗口。序列被提取到一个序列中窗口，如果序列

窗口优先选项自动创建序列窗口，它是关闭的默认

。

负载序列和调整加载序列和序列比对的序列，选择的点击率，

并将它们对齐根据爆破或PSI-BLAST查询序列校准

。序列区域以外的命中对齐通过爆炸或报道是对齐

的方法序列中的空白。

负载结构和调整负载的蛋白质结构的选定的点击，提取他们的

序列为一序列窗口，并将其对齐序列的查询序列，

根据爆炸或PSI-BLAST对齐。

负载对齐的区域加载选定的命中率从爆炸或PSI-BLAST输出。

加载顺序和对齐方式和负载结构和对齐命令加载查询序列在其全部进入序列窗口。对齐查询序列中的空白当点击加载时选择点击。

任何命令，负荷结构只有在数据库中搜索是pdb_nr95或PDB，我们这里。对于其他的数据库，执行命令将带来一个错误消息指出，PDB 标识码丢失。的PDB文件的位置指定编辑|偏好…|文件资源管理器| PDB的位置。

7、保存序列比对

我们将为未来的使用保存这一校准。

从“文件”菜单中选择“保存为。

设置类型的文件：参数序列文件。

设置文件名为blast_search。

单击“保存”按钮。

在正常的工作，你可能想保存三PDB文件尽可能好我们希望使用的模板。

8、检查序列比对

花一刻时间检查序列比对。

激活序列窗口序列结构对齐-序列窗口。

如果你向下滚动序列窗口，你会发现大量的残留物不匹配。在你看到任何一个之前，你需要在序列标尺上滚动1500实质对齐。通过寻找背景颜色的对齐方式将被发现序列。

序列窗口显示未知的命中的对齐方式。着色是基于命中与目标序列相似性的研究。

9、为下一个练习做准备

在进行下一个练习之前，我们将关闭所有的窗口。

从“窗口”菜单中选择“关闭”命令。

如果要求保存任何数据，请选择NO。

爆炸算法是一种有用的算法来寻找未知的模板。这是一个敏感从序列数据库中查找近同源序列的方法。如果序列数据库是从结构数据库PDB等，我们可以使用模板结构建立同源建模。识别更远程的同系物，PSI-BLAST可能在下一个练习中所描述的。

参考文献：

1. Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ. 1990.

Basic local alignment

search tool. J Mol Biol 215: 403

2. Altschul SF, Madden TL, Schaffer AA, Zhang J, Zhang Z, et al. 1997. Gapped BLAST

and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids

Res 25: 3389

3. Berman HM, Westbrook J, Feng Z, Gilliland G, Bhat TN, et al. 2000. The Protein Data

Bank. Nucleic Acids Res 28: 235

4. Dayhoff MO, Schwartz RM, Orcutt BC. 1978. A model of evolutionary change in

proteins. In Atlas of Protein Sequence and Structure, ed. MO Dayhoff, pp. 345.

Washington, DC: National Biomedical Research Foundation

5. Finnin MS, Donigian JR, Cohen A, Richon VM, Rifkind RA, et al. 1999. Structures of a

histone deacetylase homologue bound to the TSA and SAHA inhibitors. Nature 401: 188

6. Henikoff S, Henikoff JG. 1992. Amino acid substitution matrices from protein blocks.

Proc Natl Acad Sci USA 89: 10915

7. Khochbin S, Verdel A, Lemercier C, Seigneurin-Berny D. 2001. Functional significance

of histone deacetylase diversity. Curr Opin Genet Dev 11: 162

8. Korf I, Yandell M, Bedell J. 2003. BLAST. Sebastopol: O'Reilly & Associates, Inc. 339

pp.

9. Marks PA, Richon VM, Rifkind RA. 2000. Histone deacetylase inhibitors: inducers of

differentiation or apoptosis of transformed cells. J Natl Cancer Inst. 92: 1210

10. Nielsen TK, Hildmann C, Dickmanns A, Schwienhorst A, Ficner R. 2005. Crystal

structure of a bacterial class 2 histone deacetylase homologue. J Mol Biol 354: 107

11. Pazin MJ, Kadonaga JT. 1997. What's up and down with histone deacetylation and

transcription? Cell 89: 325

12. Pearson WR, Lipman DJ. 1988. Improved tools for biological sequence comparison.

Proc Natl Acad Sci USA 85: 2444

13. Somoza JR, Skene RJ, Katz BA, Mol C, Ho JD, et al. 2004. Structural snapshots of

human HDAC8 provide insights into the class I histone deacetylases. Structure 12: 1325

14. Taunton J, Hassig CA, Schreiber SL. 1996. A mammalian histone deacetylase related to

the yeast transcriptional regulator Rpd3p. Science 272: 408

15. Wang DF, Helquist P, Wiech NL, Wiest O. 2005. Toward selective histone deacetylase

inhibitor design: homology modeling, docking studies, and molecular dynamics

simulations of human class I histone deacetylases. J Med

Chem 48: 6936

手在练习：序列结构比对

教训1：多序列比对

当几个序列是可用的，可以执行一个多序列比对对齐他们。由此产生的校准可以用来确定哪些领域的参考蛋白是保守序列。因此，对齐表明，参考蛋白可能也被结构上保守。序列比对也被用来建立一个一对一参考蛋白的氨基酸与未知的对应关系结构上保守区域的蛋白质。这种通信是传输的基础从参考模型蛋白的坐标。并与同源建模未知和参考蛋白之间的适当序列比对是一个关键步骤（5，7，13）。

有2个多序列比对方法可在发现工作室。一种方法是align2d设有自动同源建模程序建模（9）。该方法使用一个全局动态规划算法的可变间隙罚函数这取决于插入或删除的结构上下文。align2d访问

通过序列结构比对和构建模型协议。

第二种方法，和一个用于这项运动，是align123。基于Clustal （12），一个二级结构的匹配项被添加到原来的Clustal W多序列比对得分。我们将通过多align123序列比对。

在这个练习中，我们从一组不同的序列从PSI-BLAST结果以前的练习。我们将执行序列的比对和识别的关键残基。

在这次演习结束时，你会熟悉：

?加载多个序列文件，

?进行二级结构预测，

?运行多序列比对协议，以及

?分析结果。

1开始发现工作室

我们必须有发现工作室开始行使。

推出发现工作室客户端，如果它没有运行。

如果发现工作室已经在运行，从窗口菜单，选择关闭所有命令。

如果提示保存任何分子或数据，请选择“不”。

一般情况下，我们没有看到在序列窗口中显示的任何二次结构信息。我们必须要求它显示，因为我们将在这项练习中使用该信息。

从“编辑”菜单中选择“首选项”命令。

在“首选项”对话框中，选择“序列窗口”页。

在面板的右下方显示二级结构的切换。点击确定。

2。负载未对齐的序列

根据以前的练习的结果，我们将加载一组选定的序列

与我们未知的关系最为密切。从PSI-BLAST搜索结果结合

从PDB文件1c3p序列的HDLP蛋白，和1zz1，细菌

水解酶。一个多序列比对运行（非常类似，我们将在这一运动）。的序列进行聚类分析，采用进化踪迹序列聚类分析工具。一个节点被选中HDAC1序列。冗余和非人类

序列被删除。添加远程同源1c3p和1zz1是。和

由此产生的序列集被保存。序列的完整列表是：

名称描述Swiss-Prot

加入

代码

PDB ID参考文献

HDAC1人类组蛋白去乙酰化酶1 q13547无（11）

1c3p组蛋白去乙酰化酶

蛋白（HDLP）

o67135 1c3p（2）

1t64人类组蛋白去乙酰化酶8（10）q9by41 1t64

1zz1细菌水解酶q70i53 1zz1（6）

hdac2_human人类组蛋白去乙酰化酶2 q92769无（14）hdac3_human人类组蛋白去乙酰化酶3 o15379无（1）

这些序列被编译成与空白删除单BSML文件。最后的文件是你将用什么来开始这个练习。

从文件下拉菜单，选择命令打开。

在类型的文件：下拉，选择序列格式过滤掉非顺序文件。从桌面到数据文件目录。选择文件align1.bsml。

点击打开按钮。

序列窗口将出现六个氨基酸序列，包括我们的

未知的HDAC1。序列对齐时，很少有证据表明

任何同源性。

3。二级结构预测

的一个优点是align123能力包括二级结构的援助

校准。甚至包括二级结构的预测，已被证明增加

序列比对灵敏度（3，8）。发现工作室有可以使二

结构的预测，从而可以利用align123。

确保在序列窗口中没有选择任何选择。

从“序列”菜单中选择“二级结构预测”命令。

序列窗口视图的变化，包括二级结构预测。一个文本

窗口显示从计算的原始输出。它超出了这个范围

本次研讨会在这个窗口中解释所有的条目。

关闭dscprediction文本窗口只显示序列窗口。

表明二次预测是基于从国王和斯腾伯格的DSC算法（4）

预测每个序列的二级结构。如果一个结构可用于任何

序列，我们可以使用的实际结构，以确定二级结构元素。很显然，使用真实的结构比预测更准确。

序列窗口现在显示一个二级结构的卡通。红栏显示

预测α-螺旋。蓝色箭头表示预测的链。目前，序列

都没有对齐。但我们可以作出调整，包括二级结构

信息。

4。运行一个align123对齐

你访问align123程序通过协议。

从协议资源管理器中，打开蛋白质建模协议组。

选择多序列比对协议。

参数的浏览器会出现可align123搜索设置。

在这个练习中，所有的默认值将被用于对齐，并显式地表示下面的清晰度。

对于参数名称模式，选择参数值对齐序列。

参数名称序列窗口1、选择序列窗口的名称，align1。

序列窗口1中的参数名称序列将显示所有的序列

将用于校准。

序列窗口2和序列中的序列窗口2。

成对的对齐速度应设置为慢。

多序列比对的打分矩阵可以设置BLOSUM表明BLOSUM

评分矩阵将使用。

多个对准间隙开罚应设置为10。

设置为0.05的多个对齐间隙扩展罚

我们仍然需要指定的二级结构信息是要包括在内。

单击“显示全部”按钮。

几乎在参数管理器列表的底部，选择参数使用

二级结构。设置值为真。

单击“绿色运行”按钮启动作业。

运行将完成在不到一分钟。

对于重要的align123简要解释提供以下参数。

模式指定一个简单的多序列比对，

序列的对齐方式，或对齐

将进行配置文件。一个配置文件只是一个或多个

先前已对齐的序列。

序列窗口1设置包含序列的序列窗口的名称。当模式对齐序列，序列窗口1

将所有的序列保存为对齐。

当模式是对齐序列的配置文件，

序列窗口1保存预先对齐的序列（配置文件）。当模式对齐2配置文件，序列窗口1

持有第一个配置文件。

序列序列

窗口1

允许序列的选择在序列中对齐

窗口1。

序列窗口2此参数不在设置模式设置为对齐时使用序列。

当模式是对齐序列的配置文件，

序列窗口2保存序列与配置文件对齐

在序列窗口中指定的1。

当模式对齐2配置文件，序列窗口2

持有二次配置文件。

序列序列

窗口2

当使用配置文件，允许选择的序列是

顺序窗口2。

配对比较

速度

设置此参数以快速生成初始对准

FASTA算法。

慢速设置使用一个全动态规划方法

产生一个更好的整体调整。

多序列比对

评分矩阵

指定要使用的评分函数的家庭。在

迭代法的123，评分矩阵将调整

根据序列的序列特征

考虑。

多对准间隙

开罚

在对齐的最后阶段设置间隙的开罚

过程。它将不同的对齐迭代。差距更大

开罚会产生较少的差距。

多对准间隙

延期罚款

在对齐的最后阶段设置间隙扩展罚

过程。它将不同的对齐迭代。差距更大

延期罚款会产生较短的差距。

在线帮助信息的剩余隐藏参数，

现在让我们看看结果。

5。检查结果

当“作业完成”对话框出现时，请单击“确定”按钮。

一个新的序列窗口将出现与对齐的序列。双功能应

立即明显。首先，二级结构元件在所有的

序列。二、指示序列相似性的残留显着改善。

新的显示可以作为建筑HDAC1的同源模型校准。

一些编辑的手可能是必要的，可以进行。例如，一个可以使

一个参数，大的C-末端间隙从400到450可以关闭。目前，我们

将离开的对齐方式为。

6。结束练习

在完成练习前，重置偏好。

从“编辑”菜单中选择“首选项”命令。

在“首选项”对话框中，选择“序列窗口”页。

点击重置按钮。点击确定。

现在，我们将关闭所有的窗口。

从“窗口”菜单中选择“关闭”命令。

如果要求保存任何数据，请选择号码。

在今天，我们可以使用这个信息的大量的序列信息

帮助了解蛋白质之间的关系。代表关系的一种方式是

多序列比对。一组序列的仔细校准可以高亮显示

保守的区域和可能的重要残基。结构比较

对齐，相似性和差异可以变得更加明显。

References:

1. Dangond F, Hafler DA, Tong JK, Randall J, Kojima R, et al. 1998. Differential display

cloning of a novel human histone deacetylase (HDAC3) cDNA from PHA-activated

immune cells. Biochem Biophys Res Commun 242: 648

2. Finnin MS, Donigian JR, Cohen A, Richon VM, Rifkind RA, et al. 1999. Structures of a

histone deacetylase homologue bound to the TSA and SAHA inhibitors. Nature 401: 188

3. Fischer D, Eisenberg D. 1996. Protein fold recognition using sequence-derived

predictions. Protein Science 5: 947

4. King RD, Sternberg MJ. 1996. Identification and application of the concepts important

for accurate and reliable protein secondary structure prediction. Prot Sci 5: 2298

5. Marti-Renom MA, Stuart AC, Fiser A, Sanchez R, Melo F, Sali A. 2000. Comparative

protein structure modeling of genes and genomes. Annu Rev Biophys Biomol Struct 29:

291

6. Nielsen TK, Hildmann C, Dickmanns A, Schwienhorst A, Ficner R. 2005. Crystal

structure of a bacterial class 2 histone deacetylase homologue. J Mol Biol 354: 107

7. Olszewski KA, Yan L, Edwards D, Yeh T. 2000. From fold recognition to homology

modeling: an analysis of protein modeling challenges at different levels of prediction

complexity. Comput Chem 24: 499

8. Russell RB, Copley RR, Barton GJ. 1996. Protein fold recognition by mapping predicted

secondary structures. J Mol Biol 259: 349

9. Sanchez R, Sali A. 1997. Evaluation of comparative protein structure modeling by

MODELLER-3. Proteins Suppl 1: 50

10. Somoza JR, Skene RJ, Katz BA, Mol C, Ho JD, et al. 2004. Structural snapshots of

human HDAC8 provide insights into the class I histone deacetylases. Structure 12: 1325

11. Taunton J, Hassig CA, Schreiber SL. 1996. A mammalian histone deacetylase related to

the yeast transcriptional regulator Rpd3p. Science 272:

408

12. Thompson JD, Higgins DG, Gibson TJ. 1994. CLUSTAL W: improving the sensitivity of

progressive multiple sequence alignment through sequence weighting, position-specific

gap penalties and weight matrix choice. Nucleic Acids Res 22: 4673

13. Venclovas C, Zemla A, Fidelis K, Moult J. 2003. Assessment of progress over the CASP

experiments. Proteins 53 Suppl 6: 585

14. Yang WM, Inouye C, Zeng Y, Bearss D, Seto E. 1996. Transcriptional repression by

YY1 is mediated by interaction with a mammalian homolog of the yeast global regulator

RPD3. Proc Natl Acad Sci U S A 93: 12845

手在练习：序列结构比对

Figure 1: 1FDT (3) with a

bound estradiol and NADP

2课：用3dma多重结构比对

在这个课程中，你将学习如何调整一组结构

利用结构比对算法3DMA。这个

NAD（P）结合Rossmann折叠超家族是一个多样化的