如何更改 pdb 文件的链名称?

2023-12-08

我想将 PDB 文件的链重命名为“6gch”-https://www.rcsb.org/struct/6GCH.

我检查了Biopython手册,似乎找不到任何东西。 任何输入都会有很大帮助!


你确实可以改变id链元素的属性。之后您可以使用PDBIO保存修改后的结构。

但请注意,此过程会相当多地修改 PDB。 PDBIO 不存储诸如REMARK、SHEET 和SSBOND 之类的条目。如果您知道自己需要这些,则必须小心。 此外,此过程将 HETATM 移动到相应链的末尾,而原始 PDB 将它们定位在文件的末尾。

由于 6GCH 有 3 个链,所以我使用字典renames配置旧链名称到新链名称的映射。如果该字典中不包含链名称,则不会进行重命名。

from Bio.PDB import PDBList, PDBIO, PDBParser

pdbl = PDBList()

io = PDBIO()
parser = PDBParser()
pdbl.retrieve_pdb_file('6gch', pdir='.', file_format="pdb")

# pdb6gch.ent is the filename when retrieved by PDBList
structure = parser.get_structure('6gch', 'pdb6gch.ent')

renames = {
    "E": "A",
    "F": "B",
    "G": "C"
}

for model in structure:
    for chain in model:
        old_name = chain.get_id()
        new_name = renames.get(old_name)
        if new_name:
            print(f"renaming chain {old_name} to {new_name}")
            chain.id = new_name
        else:
            print(f"keeping chain name {old_name}")

io.set_structure(structure)
io.save('6gch_renamed.pdb')
本文内容由网友自发贡献,版权归原作者所有,本站不承担相应法律责任。如您发现有涉嫌抄袭侵权的内容,请联系:hwhale#tublm.com(使用前将#替换为@)

如何更改 pdb 文件的链名称? 的相关文章

随机推荐

  • Eclipse 如何将 .java 文件作为小程序运行?

    我一直在尝试运行一个从命令行创建的简单小程序 我尝试这样做 C java Applet 显然这不起作用 但是 我注意到 如果我选择该类并选择作为 java applet 运行 则 Eclipse 允许我运行该 applet Eclipse
  • 如何将文档内容动态加载到iframe中?

    我想做这样的东西 A B B页面在A页面里面 我想把所有的js css都插入到A页面中 我想把B页面分开 我设计它使用一个iFrame 动态获
  • 检查单词是否存在两次? (正则表达式)[关闭]

    很难说出这里问的是什么 这个问题模棱两可 含糊不清 不完整 过于宽泛或言辞激烈 无法以目前的形式合理回答 如需帮助澄清此问题以便重新打开 访问帮助中心 我对正则表达式很糟糕 我想要的是检查一个字符串是否有两次 http 一词 例如 http
  • modelica:计算随时间变化的连续变量的最小值/最大值

    如上所述 我希望计算连续变量随时间变化的最小值 和 或最大值 这是一个演示的最小示例 model MinMaxTest Real u Real u min start 10 Real u max start 10 equation u si
  • 用JS发送文件到PHP服务器

    我可以使用拖放将文件从浏览器发送到服务器 并且此代码工作正常 var temp new FormData temp append file content e originalEvent dataTransfer files 0 Ajax
  • Kotlin 术语“可变”是错误的吗?

    我的理解是 变量 一词指的是重新分配引用的能力 常量 意味着不能重新分配引用 Java中final和not的本质区别 var something new obj gt reference can be re assigned val som
  • 从嵌套 json 文件中删除 python dict 项

    我有一个 JSON 文件 当我尝试删除 python 字典中的项目时 我从返回 KeyError 0 的 API 中获取该文件 我认为这是我缺乏技能和 json 格式的结合 我的目标是删除所有实例192 168 1 1 from ip ad
  • Hbase 错误“错误:KeeperErrorCode = NoNode for /hbase/master”

    在 hbase shell 中执行任何命令时 我在 hbase shell 中收到以下错误 ERROR KeeperErrorCode NoNode for hbase master 启动HBASE HOSTCHND hbase 2 0 0
  • 实现可以将提示视为实际语句吗?

    在C中 register存储限定符是hint对于实现来说 应该尽可能快地访问此类标识符 例如 存储在 CPU 寄存器中 6 7 1具有存储类说明符寄存器的对象的标识符声明表明对该对象的访问应尽可能快 这些建议的有效程度取决于实施情况 and
  • FTP从服务器直接下载文件到客户端

    我尝试将文件从 FTP 服务器下载到客户端 如果我使用ftp get 文件被下载到PHP服务器 PHP服务器可以将输出写入浏览器 所以下载过程是 FTP server gt PHP server gt client 这会使流量加倍 这对于下
  • 如何在服务方法中完全运行后台进程?

    您好 任何人都可以告诉我如何在没有活动的情况下在服务中运行此代码 我已经在活动中完成了此代码 但我不希望它成为一个应用程序 我需要它在服务中只是为了在服务上显示它 谢谢我已尝试过 但我的活动每 30 分钟显示一次 这是我的代码 public
  • SwiftUI:使用切换删除单元格时索引超出范围

    我在删除包含切换的单元格时遇到问题 我的模型看起来像这样 class Model ObservableObject Published var items Item init items Item self items items stru
  • relativelayout正在全屏显示wrap_content

    Why does FOOBARZ get layed out all the way at the bottom when no elements are layout height fill parent in other words a
  • javascript 代码在 JSF xhtml 页面中不起作用

    这是运行代码fiddle您将看到它在这里工作得很好 但是当我在 xhtml 页面中使用 glassfish server 3 2 1 在 eclipse 中运行此代码时 它会给出此错误 javax servlet ServletExcept
  • 在 Web 视图中重复使用 MSAL 的访问令牌

    在我的 iOS 应用程序中 我目前正在通过MSAL iOS SDK并且运行良好 在应用程序的某个位置 我还需要显示一个网站 通过 UIWebView 打开 该网站也需要相同的 Azure AD 身份验证 由于我已经通过 SDK 进行了身份验
  • 将 OpenMP 与 Windows SDK 结合使用

    我知道 VC2010 Express Edition 不包含 OpenMP 支持 因此会报告 omp h 文件丢失 因此 我在Windows中安装了Windows SDK v7 1 64位版本 然而 即使我跑了 set DISTUTIL U
  • 检查字符串中的特殊字符(/*-+_@&$#%)?

    如何检查字符串以确保它仅包含数字 字母或空格 在 C 中这很简单 private bool HasSpecialChars string yourString return yourString Any ch gt char IsLette
  • Prolog:为什么我的谓词返回 false?

    所以我写了一个谓词来计算一个元素在列表列表中出现的次数 count 0 base case count Elem Rest OtherLists Elem Count Elem is the head of sublist count Re
  • grep for String 并在相应行打开

    我在通过以下方式获得的几行中的多个文件中多次出现特定字符串grep grep rn include cpp mystring lib mlib actionbuttonrule cpp 300 mystring Foobar lib mli
  • 如何更改 pdb 文件的链名称?

    我想将 PDB 文件的链重命名为 6gch https www rcsb org struct 6GCH 我检查了Biopython手册 似乎找不到任何东西 任何输入都会有很大帮助 你确实可以改变id链元素的属性 之后您可以使用PDBIO保