python中利用biopython和dssp提取pdb文件中的二级结构

最新推荐文章于 2023-12-09 18:32:35 发布
最新推荐文章于 2023-12-09 18:32:35 发布
阅读量4k 收藏 8
点赞数 2
分类专栏： python 生物信息
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/yangbobor/article/details/108102335
版权
以CASP11数据集为例，选取domain中的pdb数据作为转换：在python中，安装biopython包和dssp包pip3 install biopythonconda install -c salilab dssppython脚本如下，可以提取所有sequence和ss的信息到同一个文档：from Bio.PDB import DSSP, PDBParserimport osdir='/***/casp11.domains/'list=os.listdir('/***/casp
摘要由CSDN通过智能技术生成

以CASP11数据集为例，选取domain中的pdb数据作为转换：
CASP11数据集
在python中，安装biopython包和dssp包

pip3 install biopython
conda install -c salilab dssp

python脚本如下，可以提取所有sequence和ss的信息到同一个文档：

from Bio.PDB import DSSP, PDBParser
import os

dir='/***/casp11.domains/'
list=os.listdir('/***/casp11.domains')
最低0.47元/天 解锁文章
yangbobor
关注
  • 2
    点赞
  • 踩
  • 8
    收藏
    觉得还不错? 一键收藏
  • 1
    评论
专栏目录
使用DSSP计算二级结构以及计算a-helix个数
weixin_39417324的博客
12-03 1395
安装DSSP使用anaconda： 解析氨基酸层面的dssp信息： 官网说明：https://osgeo.cn/biopython/Bio.PDB.DSSP.html https://www.docin.com/p-1973044677.html（这是资料来源，写的很清楚，很不错了）计算a-helix的个数：
Nothing
Mystra
03-31 1056
Pthon批量处理将pdb文件生成dssp文件
09-21
主要介绍了Pthon批量处理将pdb文件生成dssp文件,通过本例主要学习遍历目录下文件的方法,需要的朋友可以参考下
Python生物信息学②从PDB文件中提取蛋白序列
DrugAI
03-24 8467
环境 OS version : Win10 x64 python_version : Python 3.6.5 实例代码 aa_codes = { 'ALA':'A', 'CYS':'C', 'ASP':'D', 'GLU':'E', 'PHE':'F', 'GLY':'G', 'HIS':'H', 'LYS':'K', 'ILE':'I', '...
溶剂可及性（三）dssp和biopython
only my railgun
07-24 2347
溶剂可及性，DSSP， biopython
如何统计PDB文件二级结构的多少及位置并和作图
qq_23341021的博客
04-19 1056
很多时候我们想要展示蛋白质二级结构的多少和位置，如下简单介绍了两种可以实现此要求的方法。
SecStAnT:二级结构分析工具，用于数据选择和统计-开源
04-29
SecStAnT是一种工具，用于自动从蛋白质数据库（PDB）创建具有用户定义的结构组成的结构数据集，并计算其内部变量分布。 SecStAnT能够1.根据用户指定的二级结构（基于内部PDB分类或DSSP定义）和/或序列动机从PDB数据集中选择结构。 2.建立具有不同分辨率级别的数据集（所有原子，仅骨架，仅Cα，...）3.评估内部变量的统计分布： 单一变量分布（包括与原子表示形式最相关的分布，例如PHI和PSI，以及一些用于基于Cα的表示形式的分布）b。 两个变量的相关性（包括PHI-PSI Ramachandran图及其在基于Cα的表示中的等价关系）c。 三个变量的相关性
毕业设计 蛋白质二级结构预测 pyhton
最新发布
05-12
蛋白质二级结构预测是一个涉及生物信息学、机器学习和计算生物学的研究领域。以下是一个可能的毕业设计项目框架： ### 1. 需求分析 - **功能需求**：确定系统需要实现的基本功能，如序列输入、预测结果输出、预测...
基于机器学习的二级结构空间特征回归分析.pdf
09-24
传统的二级结构归属算法，如DSSP和STRIDE，依赖于氢键模式，但本文提出了一种新的方法，利用蛋白质相邻肽平面中心间距离这一几何特征。 在实验中，作者采用了SVM-RFE（支持向量机-递归特征消除）和McTwo（最大信息...
biopandas:在熊猫DataFrames中处理分子结构
05-01
在熊猫DataFrames中处理分子结构 链接 文档： : 源代码存储库： : PyPI： ：//pypi.python.org/pypi/biopandas 如何贡献： : //rasbt.github.io/biopandas/contributing/ 变更日志： ./ docs / sources / CHANGELOG.md 如果您是计算生物学家，则很有可能对蛋白质结构文件进行了太多次诅咒。 是的，我说的是Goode Olde蛋白质数据库格式，又名“ PDB文件”。 毫无疑问，PDB是一种结构整齐的格式（如果正确部署的话）。 但是，使用“现代”编程语言处理PDB文件有点麻烦-我很确定我们都同意这一点。 作为机器学习和“数据科学”的人，我爱上了熊猫DataFrames来处理几乎所有可以加载到内存中的内容。 那么，为什么我们不把熊猫带到结构生物学世界呢？ 在熊猫
python毕业设计 蛋白质二级结构预测+源代码+文档说明+论文
abc6838的博客
12-09 587
文件说明train.npy为训练集，test.npy为测试集main.py为主程序，运行它，就能得到我预设的效果，修改它就可以得到更多的功能mytools.py 为我写的工具箱，方便程序编写saved_model.h5 保存训练后的模型，方便多次训练a.yaml 保存了训练中 误差率 loss和准确率q3 随训练代数epoch的变化1.mp3 为训练完成的通知铃声。
利用python从蛋白文件中提取特定类型的原子
recher_He1107的博客
11-16 559
从蛋白文件中提取某种原子并写出
Pymol & BioPython | PDB文件中氨基酸序列的提取
热门推荐
DrugAI
11-01 1万+
1. Pymol 当前目录下有一个PDB文件，利用pymol的命令模式： pymol receptor.pdb -c -d "save receptor.fasta" 2.BioPython from Bio import PDB parser = PDB.PDBParser() structure = parser.get_structure('2FH7', '2FH7.p...
从蛋白质结构文件pdb中提取序列和骨架原子坐标
yang605439282的博客
01-19 5638
pdb文件是蛋白质结构存储的标准文件，用txt打开其主要内容如下所示 上图部分是model1下的蛋白原子分布，序号1所在列是原子序号，序号2所在列是氨基酸种类，序号3所在列是氨基酸序号，序号4所在三列是原子的三维坐标，序号5所在列是原子种类 用python实现提取序列： from Bio import PDB '''本文件用于提取单个pdb文件的氨基酸''' parser = PDB.PDBParser() io = PDB.PDBIO() fileName_save = 'HHH_b1_0
PDB文件：每个开发人员都必须知道的
weixin_34355559的博客
08-15 2564
PDB Files: What Every Developer Must Knowhttp://www.wintellect.com/CS/blogs/jrobbins/archive/2009/05/11/pdb-files-what-every-developer-must-know.aspx PDB文件：每个开发人员都必须知道的 一 什么是PDB文件大部分的开发人员应该都知道PDB文件是用来...
python读、写、拷贝等基本文件操作
liyali2020的博客
09-29 671
本篇简单介绍文件的读、写、拷问文件、修改名称、删除等基本操作。 python操作文件时，可直接open(“文件路径”,mode) open之后，需得close()。因为文件对象会占用操作系统的资源，并且操作系统同一时间能打开的文件数量也是有限的 python引入了with语句自动调用close()方法，这样可不用考虑底层的东西，只关注操作即可。 1 .文件权限 文件无外乎读r(read)和写w(write)，r和w文件指针默认都在从文件开头，故r+会在文件开头追加内容，w,w+会删除原内容。所以还有一个追加
Linux下安装TMalign、IUpred2a和DSSP
续起大胡子
01-09 1316
Linux中安装TMalign、IUpred2a和DSSP
蛋白质结构预测与分析
qq_40459859的博客
06-05 1万+
蛋白质的结构 一级结构 氨基酸序列 二级结构 周期性的结构构象 三级结构 整条多肽链的三维空间结构 四级结构 多个亚基形成的复合体结构 蛋白质的二级结构 蛋白质经过折叠后会形成规则的片段，这些规则的片段构成了蛋白质的二级结构单元。 图 1蛋白质的二级结构 DSSP指认 DSSP（Definition of Secondary Structure of Proteins），即蛋白质二级结构定义词典。DSSP并不预测二级结构，而是根据二级结构的定义对已经测定三级结构的蛋白质的各.
用python写一段代码，用dssp方法从pdb文件中计算蛋白质球度
03-24
以下是用Python编写的代码，使用DSSP方法从PDB文件中计算蛋白质球度： ```python import Bio.PDB import os # 设置PDB文件路径 pdb_file = "protein.pdb" # 创建PDB解析器 parser = Bio.PDB.PDBParser() # 读取PDB文件 structure = parser.get_structure("protein", pdb_file) # 创建DSSP计算器 model = structure[0] dssp = Bio.PDB.DSSP(model, pdb_file) # 计算蛋白质球度 q3 = Bio.PDB.Polypeptide.three_to_one(dssp[0][1]) for i in range(1, len(dssp)): aa = Bio.PDB.Polypeptide.three_to_one(dssp[i][1]) if aa == "X": continue q3 += aa q3 = q3.replace("G", "") q3 = q3.replace("P", "") q3 = q3.replace("A", "") q3 = q3.replace("V", "") q3 = q3.replace("L", "") q3 = q3.replace("I", "") q3 = q3.replace("M", "") q3 = q3.replace("F", "") q3 = q3.replace("Y", "") q3 = q3.replace("W", "") q3 = q3.replace("S", "C") q3 = q3.replace("T", "C") q3 = q3.replace("N", "C") q3 = q3.replace("Q", "C") q3 = q3.replace("H", "C") q3 = q3.replace("K", "C") q3 = q3.replace("R", "C") q3 = q3.replace("C", "H") q3 = q3.replace("D", "E") q3 = q3.replace("E", "E") q3 = q3.replace("B", "X") q3 = q3.replace("Z", "X") q3 = q3.replace("J", "X") q3 = q3.replace("O", "X") q3 = q3.replace("U", "X") q3 = q3.replace("X", "-") print(q3) ``` 请注意，此代码需要安装Biopython库才能运行。

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 1
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值