核酸结构分析及建设3DNA

摘要

追梦的软件包是一个流行和通用的生物信息学工具与功能分析，构建，可视化三维核酸结构。本文提出了详细的协议的一个子集3DNA提供新的和流行的特点，适用于个别结构和相关结构的合奏。

摘要

追梦的软件包是一个流行和通用的生物信息学工具与功能分析，构建，可视化三维核酸结构。本文提出了详细的协议的一个子集3DNA提供新的和流行的特点，适用于个别结构和相关结构的合奏。协议1列出下载并安装该软件所需要的指令的集合。其次，协议2，通过分析的核酸结构，包括碱基对的分配和确定刚体参数描述结构，协议3，由一个原子的重建的说明模型的结构从它的刚体参数。追梦，版本2.1，最新版本的结构歌舞团的分析和操作的新功能，如核磁共振（NMR）测量和分子动力学（MD推导出）模拟，这些功能在协议中的4和5。除了 ​​追梦的独立软件包，w3DNA位于http://w3dna.rutgers.edu ， web服务器，提供了一个友好的用户界面功能的软件选择。协议6演示了一个新的功能的网站建设模型长的DNA结合蛋白的分子修饰在用户指定的位置。

引言

了解三维结构的DNA，RNA和蛋白，药物和其他配体的配合物，用于解密其不同的生物功能是至关重要的，允许设计合理的治疗学。这种结构的探索需要三个单独的，但密切相关的部分组成：分析（提取型态的形状和相互作用），建模（评估能量和分子动力学），和可视化。结构分析和模型建设基本上同一枚硬币的两面，他们两个补充和可视化。

3DNA一套计算机程序，是一个日益流行的结构与功能的生物信息学工具包，分析，构建和可视化三维核酸结构。早期出版物概述软件^1，提供的食谱执行所选任务²的功能，介绍了基于Web的界面^软件的流行特点，结构特点数据库收集使用 3DNA ^4，5和示出的软件的实用程序在分析中的DNA和RNA结构^6，7。

本文目的是把的3DNA软件套件实验室科学家和其他人的利益和/或需要调查的DNA和RNA的空间组织与国家的最先进的计算工具。这里提出的协议包括一步一步的说明（一）下载并安装一个Mac OS X系统上的软件，（II-III）DNA碱基对组成的步骤，水平结构分析和修改（ IV-V），以分析和相关的DNA结构调整套，及（vi）与用户友好的w3DNA Web界面蛋白质修饰过的DNA链构建模型。软件具有分析能力，解决个别结构，利用X射线晶体学方法以及大型利用核磁共振（NMR）方法测定结构或计算机仿真技术所产生的合奏。

这里的结构检查包括（i）的高分辨率晶体结构的DNA势必到HBB蛋白从伯氏疏螺旋体 ^（蜱传播的细菌，导致莱姆病在人类^{9，10），（}二）两个大组顺序相关的DNA分子与分子模拟¹¹ - 4500快照d的（GGCAAAATTTTGCC）的_，2和_{d（CCGTTTTAAAACGG）2}在100皮秒为单位收集，在计算过程中，及（iii）O3 DNA操作员基于NMR的结构的一个小的合奏绑定的头饰在大肠杆菌 lac阻遏蛋白^12。下面的说明包括有关如何访问这些结构，以及如何使用3DNA（此文件的一个副本被发现相关文件的原子坐标信息3DNA在http://forum.x3dna.org/jove论坛）审查和修改这些结构。

研究方案

1。安装软件包

1. 连接到的追梦3DNA在http://x3dna.org网站和单击链接到3DNA论坛。在论坛内，选择“注册”链接，并按照说明进行操作，以创建一个新的帐户。
2. 以下说明详细安装一个基于OS X的Macintosh电脑，一个默认的bash'壳上的软件。 Linux或Windows（Cygwin的的MinGW / MSYS）系统常用shell（包括'tcsh的'）的程序内发现3DNA论坛。
3. 一旦你已经建立了用户名和密码，登录进入论坛。欢迎部分，在新页面上选择“下载”链接选择'3 DNA的下载“。这会带给你追梦的下载页面，可以下载各种版本的软件。
4. 选择Mac OS X的英特尔链接下载一个压缩的文件（tar.gz的）包含软件。
5. 双集团公司k对（tar.gz）文件来创建一个文件夹名为-V2.1 x3dna的。这个文件夹，其中包含的3DNA软件套件，可以移动到任何位置。对于这个食谱的目的，我们将它拖放到“应用程序”目录。在这个阶段，你已完成的tar.gz文件，可以删除它。
6. 打开终端应用程序，通常发现在“实用工具”下，更改到V2.1 x3dna的第5步中创建的目录，使用以下命令（在这里结束，下面的说明一个回车）：
   CD / Applications/x3dna-v2.1的
7. 3DNA正常打字需要运行安装脚本：
   ./bin/x3dna_setup
8. 将步骤7中印两线出口设置。如果用户已经在“应用程序”目录中放置软件，两行应显示如下：
   出口X3DNA = / Applications/x3dna-v2.1号
   出口PATH = / Applications/x3dna-v2.1/bin：$ PATH
   一个不同的'应用'字，除了将APPEAR在出口命令，如果用户选择安装的软件在不同的目录。此备用目录名称会取代在上述两个命令的'应用'。请注意，如果默认的shell是tcsh的'，而不是'庆典'，两行的内容将是：的setenv X3DNA / Applications/x3dna-v2.1和SETENV路径/ Applications/x3dna-v2.1/bin：$ PATH中。
9. 请在“搜索”看，如果一个文件名为'。bashrc'中已经存在你的电脑上。
10. 打开TextEdit应用程序在“应用程序”目录下的“格式”菜单中选择“纯文本”。
11. 如果'。bashrc'中存在文件从第9步，在打开文件TextEdit和导出设置粘贴步骤8至文件末尾。如果没有'。bashrc中的文件退出，粘贴导出设置到一个空白文件，并从第8步保存的文件名'。bashrc中'在你的home目录。 （主目录是由“房子”在Macintosh上的图标表示。）
12. 为了让新的设置到bê，你必须运行新保存'。bashrc'中文件内的终端程序，使用下面的命令：
    bashrc中源〜/。
13. 追梦的套房现在已经安装在你的系统内的终端程序，您可以执行协议2-5。协议6，一个网络接口使用w3DNA追梦软件选择的特征。

2。晶体结构分析

1. 我们说明使用莱姆病螺旋体 ⁸ HBB蛋白结合的DNA的核酸结构的分析。与这个结构的原子坐标文件，可以发现蛋白质数据银行^13（PDB http://www.rcsb.org ），它是在结构标识符2np2的分配。随附的补充资料包括一份这个文件，也可以找到http://forum 3DNA论坛。x3dna.org/jove。
2. 追梦的结构分析的第一步是创建一个文件，该文件列出了所有的配对基地。这是通过“find_pair程序通过键入以下内容：
   find_pair 2np2.pdb 2np2.bps
   这里2np2.pdb文件的原子坐标上述2np2.bps的是输入输出的文本文件，其中包含的碱基配对信息。后者碱基对文件由用户如果有必要，可以检查和编辑。
3. 在分析中的下一个步骤是确定所特有的结构的几何参数。这些包括标准的化学扭转角^{14，pseudorotation}角度描述的褶皱的糖环^15，刚体参数指定基地和连续的碱基对¹⁶等措施，立体化学结构的安排。命令'分析'作为输入碱基对步骤2中生成的文件，并创建几个OU的TPUT文件，出现在当前工作目录。计算这些值是通过输入以下命令：
   分析2np2.bps，
   用在这里的输出文件包括的2np2.out，其中包含计算的参数的概述，bp_step.par，其中包含刚体上述参数的列表。后者的文件可以方便地阅读到电子表格程序，作进一步的分析和策划。代表结果（ 图1）的一个例子。

3。从硬体参数构建的DNA结构

1. 除了一个核酸结构的分析，追梦的软件提供的能力，建立一个结构模型刚体使用'重建'命令的参数。在这个协议中，我们将展示如何构建两个DNA螺旋模型，首先用刚体参数计算协议2 HBB-DNA复合物，然后使用修改后的参数设置在w一族选定的滚动角被改变。侧倾角测量连续的碱基对的DNA进入槽¹⁶的弯曲程度。轧辊的正的值相关联的与缩小的大沟和负值的小沟变窄。的是标准基准帧¹⁷所通 ​​过3DNA和其他流行的核酸分析软件， 例如的曲线+ ^18，可确保数值的一致性Watson-Crick碱基对中的计算参数。
2. 3DNA“重建”功能默认情况下，将创建只包含含氮碱基中的原子坐标的精确的原子模型。为了创造一个近似的原子模型，其中包括骨干和基础原子的坐标，输入以下命令：
   x3dna_utils cp_std BDNA
   此命令指示3DNA刚体参数引入一个标准的B型DNA骨干构建设模式。
3. 裴LD的双螺旋模型，刚体HBB-DNA结构中找到的参数，键入以下内容：
   重建原子bp_step.par 2np2_rbld_org.pdb
   在这里，'原子'指定的全原子模型的建设与所有重型（非氢）原子的原子坐标，'bp_step.par'议定书2第3步中生成的文件，其中包含了输入的刚体参数，'2 np2_rbld_org.pdb'是输出文件创建的结构的原子坐标。后者按照同规格的PDB文件格式，从而终止。PDB'。
4. 可编辑'bp_step.par的文件，生成一个修改过的结构。在Macintosh上的变化可以用TextEdit应用程序中所描述的协议1步骤10。
5. 在这里，我们修改两个站点的最大弯曲极端辊形成的夹角。值，以度表示，从64.95的值改为0，第17行和600.93第26行。我们用一个新名称保存修改后的文件，'bp_step_roll0.par'，关闭的TextEdit程序。
6. 我们建立的结构，在步骤3中，通过打字输入文件与修改后的一套步骤参数：
   重建原子bp_step_roll0.par 2np2_rbld_roll0.pdb
7. 可以被视为在一个标准的分子观者如PyMOL的（ www.pymol.org ），这两个模型，使用生成的坐标文件（的。pdb）作为输入。请参阅图像重建结构的代表性的结果（ 图2）。

4。多模型结构文件分析

1. 协议2，在此我们分析一个单一的含核酸结构，在这里我们将展示如何分析结构大合奏。我们研究了随着时间的推移沿着一系列的14个碱基对的DNA链得到的分子动力学（MD）模拟¹¹的刚体参数的变化。分析认为，4500结构，间距为100皮秒的增量在MD计算，4,500单独的文件保存在。随附的补充数据包括这些文件的副本，以及反向的DNA序列（见碱基的确切顺序的介绍性文字）模拟了一套4500的其他文件。两套文件和一个较小的数据子集，也可以发现3DNA在http://forum.x3dna.org/jove论坛上。
2. 在分析中的第一个步骤是通过运行的'find_pair程序进行协议2，识别配对的碱基相同。由于整个组的结构股相同的碱基配对方案一需要运行'find_pair'，只有一次有代表性的文件。在这里，我们选择的第1001个结构在DNA链含中央步骤，的文件名md_AT_set1.pdb.1001的模拟。下面的指令属韭碱基配对信息，并将其存储在文件中md_AT_set1.bps：
   find_pair md_AT_set1.pdb.1001 md_AT_set1.bps
   用户可以检查所确定的碱基配对和知识结构的基础上在文件中进行手动更改。
3. 下一个步骤是确定整组使用“x3dna_ensemble分析”命令的结构的几何参数。该程序使用作为输入碱基对生成的文件在步骤2中（md_AT_set1.bps）以及一个文件包含要分析的结构的列表。后面的文件，这里称为md_AT_set1.list，包含4500个文件进行分析，在个别行上输入的文件名中的顺序生成的分子动力学模拟。这个文件的副本包含在补充材料，也提供3DNA论坛。用户应注意，由于大量的文件，计算的时间会很长（约20分钟AMD羿龙II X4 965处理器，每1000个文件）。钍e命令运行通过键入以下内容：
   x3dna_ensemble分析 - B md_AT_set1.bps的-L md_AT_set1.list邻md_AT_set1.out
   在上面的例子'-B'涉及到碱基对在第2步中创建的文件，“-L”指定的文件列表，'邻'允许用户指定一个名称到输出文件中。也有很多选择'x3dna_ensemble分析“命令，可以配合使用，如具体型号的选择。用户可以键入'x3dna_ensemble分析-H“来获取有关这些选项的更多信息。
4. 下一步骤是提取所需的刚体参数从概述的输出文件，md_AT_set1.out，在步骤3中产生。这是使用“x3dna_ensemble提取物”命令执行。一上市可用的几何参数，可以发现输入'x3dna_ensemble提取L'。在这里，我们提取的滚动角，并创建一个以逗号分隔的文本文件（CSV）的侧倾角度用以下命令：
   x3dna_ensemble提取物 - 对滚-S，-F md_AT_set1.out邻md_AT_set1_roll.csv
   '-P卷在此命令中指定要提取的参数，“-S”表示一个逗号分隔的文件，'-F'表示合奏分析命令的输出文件，和'邻'允许用户在所有分析的结构的DNA命名为CSV输出文件中的侧滚角沿。这个文件可以读取到其它程序中作进一步的分析和策划。请参阅代表结果（ 图3），在上述两个MD生成的数据集和用于生成这些图像在MATLAB脚本3DNA论坛辊的变化进行分析。

5。多模型结构的叠加到一个共同的参考框架

1. 追梦（2.1版）的最新版本有一个新功能，它允许用户在多个结构从一个普通的角度看。 'x3dna_ensemble重新定位“命令相关结构的叠加集合上课在一个共同的碱基对或碱基对的步骤。以下协议适用于这种能力O3 DNA运营商绑定在Lac阻遏蛋白¹²的头饰核磁共振派生结构。下的标识符2kek的多模式结构的文件，该文件包含了所有的结构模型在一个文件中，而不是单独的文件中，例如用于每个协议4的4500型号，这些结构被存储在PDB中。随附的补充资料包括一份这个文件，也可以找到3DNA在http://forum.x3dna.org/jove论坛。
2. 正如在先前的协议中，第一步是计算的DNA的碱基配对可与'find_pair“程序获得。在这种情况下，我们使用输入文件kek.pdb '2'，并键入以下内容：
   find_pair 2kek.pdb 2kek.bps
   碱基配对产生PDB文件使用的第一款车型。手动更正的输出bASE对的文件，'2 kek.bps“，可以进行知识结构的基础上。
3. 该方案的x3dna_ensemble重新定向'对齐在一个多模型结构的一个共同的参考帧的单独的模型。此程序需要pdb文件和碱基配对的文件作为输入，并从以前的步骤，运行以下命令：
   重新调整-B 2kek.bps x3dna_ensemble F 1-E 2kek.pdb邻2kek_fr1.pdb
   在这里'-B'选项指定的碱基配对文件，参数“-f 1'表示的碱基对上所有的结构是一致的，在这种情况下，碱基对的1序列息链的5'-末端'-E'指定的输入集合文件，“-o”参数允许用户提供一个输出文件名，在这种情况下，'kek_fr1.pdb'。有关命令选项的更多信息，可以通过键入'x3dna_ensemble重新调整H'。请参阅代表结果（ 图4）对齐的结构的讨论。

6。工程业CTION的蛋白质修饰过的DNA分子

1. 的w3DNA Web界面包括一些流行的功能的3DNA软件包。在这里，我们提请注意的能力，构建三维模型与Web服务器的蛋白质修饰过的DNA。结合的DNA片段采用刚体参数所选择的蛋白质-DNA复合物，这里的值HBB 2协议分析的蛋白质的DNA结合。未绑定的DNA区域，可以采取三个用户选定的螺旋形式（A，B或C）之一。 A，B和C是指有三个典型的DNA模型，从早期的光纤衍射实验，分别为11，10，每回合的双螺旋^19-21和9个碱基对。
2. 构建蛋白质修饰过的DNA模型的第一步是位于http://w3dna.rutgers.edu，参观w3DNA网站。从菜单中选择“重建”的顶部左侧的页面激活vates在线模式建设能力。
3. 下一个步骤是选择在浅色阴影的框中，然后在新的页面的中间发现的'束缚的蛋白质与DNA模板“链接。此选择将激活一个下拉菜单，该菜单允许用户指定结合的蛋白质的数量。在这里，我们选择了两个结合蛋白，并单击“继续”按钮。请注意，模型的大小受限制，以1,000个碱基对，或2000个核苷酸，和第三结合蛋白。
4. 结合的蛋白质的数量的选择导致在图5中所示的类似的规范页面。在该页面允许用户指定的DNA序列，螺旋形式的未结合的DNA结合蛋白的身份和位置的。
5. 用户类型或粘贴在中间的规格页中的文本框（ 图5，标签1）按照期望的顺序。注意唯一标准的残留物 - A，T，C，G，U - 被接受。在这里，我们进入了81个碱基对的坎opolymer，完全的A·T对组成序列息链中的A。
6. 有一个下拉菜单（ 图5，标签2），下面的文本框中，选择未绑定的螺旋构象的DNA区域。在这个例子中，我们选择的B型构象。
7. 底部进入的规格页的装订位置（s）和文件名 ​​（S）结合蛋白（S）（ 图5，标签3）附近有一个文本框。的结合位置指定的蛋白质结合的片段的DNA序列上的中心的位置。如果结合的DNA包含奇数个碱基对，如六溴代二苯DNA结构中的结合位置表示的位置的中间的碱基对的片段。 HBB-DNA结合的情况下，这中间的碱基对是T·T mispair。在这个例子中，在这里我们将两个拷贝的贫血症状的蛋白质的DNA模板，放置在位置20和63的81个碱基对的DNA序列中的T·T对序列。请注意，绑定的片段的序列不具有配合在步骤5输入。如果蛋白质结合的区域包含偶数数目的碱基对中的结合位置表示的位置上的DNA碱基对的中间步骤。也就是说，中央结合的片段的碱基对的步骤，n和n +1，其中n是沿段的指定数目的碱基对之间。另外请注意，用户可以与任何蛋白质修饰的DNA，只要其与DNA的复合物的结构是已知的并存储在标准PDB格式^13。
8. 在规范的页的底部是一个盒子，它允许用户生成的蛋白结合的DNA（ 图5，标签4）的预览。在这里，我们请选中该复选框并单击“继续”按钮。这个动作会产生一个查看页面，列出所选的参数，以及可能出现的任何错误，从选定的结合位点的重叠。用户可以LECT“返回”按钮，如有需要作出改变或'建立'按钮继续。
9. 翻页显示的静态图像的DNA-蛋白质复合物在其大多数扩展的安排，并允许上线交互式可视化通过Jmol的Webmol的。用户还可以下载一个坐标（PDB）文件，按照同规格的PDB格式。这是通过选择“重建pdb文件生成的图像下面的链接下载。到一个分子的进一步可视化的图形程序，后者可以方便地阅读文件。代表结果的的贫血症状装饰DNA的例子（ 图6）如上所述，一个稍微长链（86个碱基对的）在位置20和67的中心的HBB。

结果

常规使用的的3DNA软件工具来分析核酸结构。例如，碱基对的身份和特征在双螺旋片段的碱基的DNA和RNA结构的安排刚体参数自动计算，为每一个新的条目存储在核酸数据库^22，世界范围内的存储库核酸的结构信息。刚体参数确定协议2的值轻易表露在三维结构的扭曲，如极端的DNA大沟，大正滚动角（64.95°，60.93°）弯曲成两个网站，发现在AT·AT与伯氏疏螺旋体 HBB蛋白^8（ ?...

讨论

在这篇文章中提出的协议集，只触及3DNA程序套件的能力。的工具，可以适用核糖核酸结构识别非经典的碱基对，以确定的二级结构的上下文中，会发生这样的配对，量化螺旋片段的空间配置，测量沿链主链的碱基重叠，等等。 rebuild命令允许用户构建简单和翔实的基地和基地对像中的插图所示， 图1块表示。建设工具还包括在一个给定的结构化模板功能，以“线”不同的序列，产生无...