Bio. DBC.SASA模块
Bio. DBC实体的溶剂可及表面积的计算。
使用Shrake & Rupley算法开发的“滚球”算法,该算法使用球体(与溶剂分子半径相等)来探测分子表面。
- 参考:
Shrake,A; Rupley,JA.(1973)。J Mol Biol“蛋白质原子的环境和溶剂暴露。溶菌酶和胰岛素”。
- class Bio.PDB.SASA.ShrakeRupley(probe_radius=1.40, n_points=100, radii_dict=None)
基类:
object
使用Shrake-Rupley算法计算SASA。
- __init__(probe_radius=1.40, n_points=100, radii_dict=None)
初始化课程。
- 参数:
probe_radius (float) -- 探头半径A。默认值为1.40,大约是水分子的半径。
n_points (int) -- 每个原子表面的分辨率。默认为100。更多的点会导致更精确的测量,但会减慢计算速度。
radii_dict (dict) -- 用户提供的原子半径字典用于计算。值将替换/补充默认ATOMIC_RALDI字典中的值。
>>> sr = ShrakeRupley() >>> sr = ShrakeRupley(n_points=960) >>> sr = ShrakeRupley(radii_dict={"O": 3.1415})
- compute(entity, level='A')
计算实体的表面可达性表面积。
生成的原子表面可访问性值附加到每个实体(或原子)的.sasa属性,具体取决于级别。例如,如果level=“R”,则所有残基都将具有.sasa属性。Atoms将始终被分配具有其个人值的.sasa属性。
- 参数:
entity (Bio.PDB.Entity) -- 输入实体。
level -- 分配ASA值的级别,可以是“A”(Atom)、“R”(剩余)、“C”(链)、“M”(模型)或“S”(结构)之一。实体的ASA值是其子实体所有ASA值的总和。别名为“A”。
>>> from Bio.PDB import PDBParser >>> from Bio.PDB.SASA import ShrakeRupley >>> p = PDBParser(QUIET=1) >>> # This assumes you have a local copy of 1LCD.pdb in a directory called "PDB" >>> struct = p.get_structure("1LCD", "PDB/1LCD.pdb") >>> sr = ShrakeRupley() >>> sr.compute(struct, level="S") >>> print(round(struct.sasa, 2)) 7053.43 >>> print(round(struct[0]["A"][11]["OE1"].sasa, 2)) 9.64
- __firstlineno__ = 73
- __static_attributes__ = ('_sphere', 'n_points', 'probe_radius', 'radii_dict')