Bio.PDB.SASA模块¶
Bio.PDB实体溶剂可及表面积的计算
使用Shrake&Rupley算法开发的“滚球”算法,该算法使用一个(与溶剂分子半径相等的)球体来探测分子的表面。
- 参考资料:
题名/责任者:A.(1973)。J·摩尔生物:“环境与蛋白质原子、溶菌酶和胰岛素的溶剂暴露”。
- 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) -- 默认情况下,探针的半径是1.40,大致相当于水分子的半径。
n_points (int) -- 每个原子表面的分辨率。默认值为100。点数越多,测量越精确,但会减慢计算速度。
radii_dict (dict) -- 用户提供的用于计算的原子半径字典。值将替换/补充默认ATOM_RADII字典中的值。
示例:>sr=shrakerupley()>sr=shrakerupley(n_point=960)>sr=shrakerupley(Radii_dict={“O”:3.1415})
- compute(entity, level='A')¶
计算实体的表面辅助功能表面积。
生成的原子表面可访问性值附加到每个实体(或原子)的.sasa属性,具体取决于级别。例如,如果LEVEL=“R”,则所有残留物都将具有.sasa属性。原子将始终被分配一个.sasa属性及其各自的值。
- 参数:
entity (Bio.PDB.Entity) -- 输入实体。
level -- 指定ASA值的级别,可以是“A”(原子)、“R”(残渣)、“C”(链)、“M”(模型)或“S”(结构)之一。实体的ASA值是其子项的所有ASA值之和。默认为“A”。
示例:>从Bio.PDB导入PDBParser>p=PDBParser(Quiet=1)>struct=p.get_Structure(“1LCD”,“pdb/1LCD.pdb”)>sr=ShrakeRupley()>sr.computer(struct,level=“S”)>Print(round(struct.sasa,2))7053.43>打印(round(struct [0] ["A"] [11] [“OE1”] .sasa,2))9.64