| 查看: 566 | 回复: 1 | |||
| 当前主题已经存档。 | |||
cad_0木虫 (著名写手)
玛雅城的小石头
|
[交流]
rasmol使用方法
|
||
|
rasmol使用方法 作者:swallowx 整理时间:2005-08-31 目前在结构生物学领域有许多图形显示的程序，每个都有自己不同的特点。可能很多人知道rasmol，除了图形界面中的一些功能外，该程序的命令行方式有着很强大的功能。下面将介绍一些常用的使用方法。 Rasmol（http://www.openrasmol.org/）程序有多种版本，有unix, windows, Mac等。另外还有一个windows和linux版本http://www.geneinfinity.org/rastop/，该版本将一些原来rasmol菜单下没有的命令加入菜单，比原来的版本方便。 在unix系统下，读入结构文件可以直接用命令的方式，如 rasmol 1crn.pdb. 在windows下，可以先打开raswin, 然后在File的菜单下读入结构文件。Rasmol 所识别的文件有下面几种： pdb： Brookhaven Protein Databank，来源于 www.rcsb.org nmrpdb NMR multi-pdb file format mopac mopac file format; either cartesian or z-matrix format。 mdl Molecular Design Limited's MOL file format mol2 Tripos' Sybyl Mol2 file format。 xyz MSC's XMol XYZ file format。 alchemy Tripos' Alchemy file format。 charmm CHARMm file format。 如果想读入Charmm软件包的结构文件，用命令行的方式是：rasmol –charmm 1crn.crd. 在打开rasmol后，会出现两个窗口，一个是图形窗口，另外一个是命令行窗口。可以在图形窗口中进行一系列的结构操作，但是有些的操作还需要命令行来补充。下面就将常用的命令进行一下总结。 restrict protein： 在图形窗口中去除所有的非蛋白质原子。 restrict lys: 在图形窗口中去除所有的非lys残基。 select all: 选择所有的原子。 select protein： 选择蛋白质原子。 select hetero： 选择非蛋白质，非DNA原子。 一、常用的选择命令： 1. 结构文件中链的选择： 每条链都有一个字母或数字表示。选择一条链时，必须用:或*说明字母代表链。比如： select :d 选择d链的所有原子 select *d 选择d链的所有原子 select :d,:e 选择d或e链的所有原子 select glu:2 选择2号链的所有glu。 2. 通过残基名称选择 PDB文件中每个残基都有1-3个字符串的名称。所有的氨基酸用3字符表示，DNA，RNA用单字符表示。水分子用HOH表示。其他的配体的名称可以用文本编辑器打开PDB文件搜索，配体原子对应的坐标在文本中由HETATM开头，而蛋白质或DNA的原子是以ATOM开头，见下面的例子： ATOM 1902 O GLY R 62 -32.180 -32.765 46.907 1.00 38.84 HETATM 1955 O HOH 1 -26.069 -22.429 17.059 1.00 53.88 有些基团的名称含有数字，如SO3，PO4。在选择这些残基时，残基需要加上方括号，如select [SO3]. 钙原子或其他金属原子一般用2字符表示，如CA, MN,MG, ZN. 如果需要选择钙原子，可以用”select ca”. 钙原子ca的表示与蛋白质的a碳原子表示冲突，因此如果PDB结构中有钙原子，只想选择蛋白质的a碳原子，可以用select protein and *.ca. 另外一些例子有： select lys:a 选择A链的所有lys. Select (lys,arg) and :b 选择B链的lys或arg。 3. 通过残基的数字选择 每个残基都有一个数字相对应，下面是一些选择的例子： Select 32 选择每条链的32号残基以及32号杂原子 Select 19-32 选择每条链的19-32号残基以及19-32号杂原子 select 19-32 and not hetero 选择每条链的19-32号残基 select 19-32 and hetero 选择19-32 号杂原子 select 19-32:y 选择Y链的19-32号残基 select asp47 选择所有链的47号位置的asp 4. 原子的选择 PDB文件中每个原子都有一个序列号对应，可以在图形界面上点击，查看序列号。如果想选择原子可以用： Select atomno = 131 选择第131个原子 Select atomno = 217, atomno = 1426 选择第217和1426号原子 select atomno >= 195 and atomno <= 277 选择195-277号原子 PDB原子名称，PDB文件中所有的原子的命名采用标准命名，CA表示a－碳原子CB表示b－碳原子，以此类推。CG, CD, CE, CZ, CH (gamma, delta, epsilon, zeta, eta)。N7 (7th nitrogen in a residue), O2P (second oxygen on a phosphorus), OE2 (second oxygen on an epsilon carbon), HD1 (1st hydrogen on a delta carbon). 可以在图形界面中点击原子，然后在命令行窗口上读取原子信息。其他实例： select *.cg 选择所有g位的碳原子。 select lys.cg 选择所有赖氨酸的g位的碳原子 select :a.cg 选择A链所有g位的碳原子 select lys:a.cg 选择A链所有赖氨酸的g位的碳原子 select 27-42:a.cg 选择A链从27到42位氨基酸的g位的碳原子 select *.h? 选择所有2字符的氢原子 select *.h??? 选择所有的氢原子。 元素名称，键入全名，如magnesium, iron, sulfur等。 另外一些选择实例： 用户可以自定义一些区域进行操作，比如： define activesite（15,67,109） select activesite color green 另外可以通过下面的命令选择某个残基周围的其他原子： select within（4.5, ser72) 选择ser72周围4.5?内的原子。 如果不想选择该范围内的某个残基，可以用: select within（4.5, ser72) and not lys80 二、进一步的操作： 在选择了一些基团 后，就可以对它们进行进一步的操作，比如修改残基的表示方法，可以变成球棍模型，空间堆积模型等等。也可以对它们进行不同的着色。rasmol提供了多种图形显示方法，对原子的显示有wireframe，spacefill，sticks，ball and stick，对于二级结构的显示有ribbons，strands和cartoons。另外还有backbone的显示方法。 在选择了一个残基后，可以有下列操作： select 172:A color green 将该残基着绿色 wireframe 0.5 数字表示一个相对值，此时可以看见该残基变粗了。 此时可以在图形界面中的display菜单中选择sticks ，spacefill或ball and stick的方式。然后在命令行窗口可以进行参数设定，比如，在将残基变成ball and stick方式后，在键入 spacefill 0.3， 可以发现原子的表示比原来要小，如果在键入wireframe 0.1，可以发现，化学键的表示变细了。通过这种方法，可以很容易区分所感兴趣的残基。 |
回复此楼» 猜你喜欢
08工学 309分求调剂
已经有7人回复
-
各位老师好，求调剂，本科211，一志愿天津大学生物与医药学硕，差两名录取。
已经有9人回复
-
考研调剂
已经有12人回复
-
材料相关专业344求调剂双非工科学校或课题组
已经有18人回复
-
297求调剂
已经有24人回复
-
085600材料与化工329分求调剂
已经有19人回复
-
310求调剂
已经有17人回复
-
求调剂，985材料与化工348分
已经有14人回复
-
材料工程085601，270求调剂
已经有42人回复
-
279求调剂
已经有11人回复
» 本主题相关商家推荐: (我也要在这里推广)
1
 |
2楼2006-10-11 16:52:14
