博客
关于我
pysmiles:一个用于读写SMILES表达式的python库
阅读量:529 次
发布时间:2019-03-08

本文共 870 字,大约阅读时间需要 2 分钟。

SMILES 表达式与 pysmiles 的应用

技术背景

SMILES(SMILES Notation)是化学领域中用于描述分子结构的简洁文本表达方式。通过特定的符号规则,SMILES能够准确表示分子中原子的连接关系及相关属性。例如:

  • 氢原子通常省略,只显示原子符号(如 C 表示碳)。
  • 双键和三键分别用 =# 表示。
  • 分支关系用圆括号表示,主链不加标记。
  • 不成键的分子用 . 分隔,环状结构用数字标记开头和结尾。

这种表达方式简洁高效,广泛应用于化学信息处理中。

pysmiles 的安装

pysmiles 是一个基于 Python 的 SMILES 解析与生成工具,结合 networkx 框架存储分子结构信息。安装过程如下:

pip install pysmiles

安装完成后,可以通过 Python 环境调用 pysmiles API。

pysmiles 的使用

以下是一个简单的 SMILES 解析示例:

from pysmiles import read_smilesimport networkx as nxsmiles = 'N#CC#N'mol = read_smiles(smiles)print(mol.nodes)  # 打印节点信息print(mol.edges)  # 打印边信息print(nx.to_numpy_matrix(mol))  # 打印邻接矩阵elements = nx.get_node_attributes(mol, name="element")nx.draw(mol, with_labels=True, labels=elements)plt.savefig('pysmiles.png')

运行该脚本会生成一个分子结构图,图中标注原子元素。

总结

pysmiles 提供了便捷的 SMILES 解析与生成功能,结合 networkx 实现分子结构可视化。它适合用于化学信息处理中的结构分析与可视化。

如需进一步探索 SMILES 的应用,可参考相关工具的官方文档或示例代码。

转载地址:http://yjgnz.baihongyu.com/

你可能感兴趣的文章
ORCHARD 是什么?
查看>>
Struts2中使用Session的两种方法
查看>>
Stream API:filter、map和flatMap 的用法
查看>>
STM32工作笔记0032---编写跑马灯实验---寄存器版本
查看>>
order by rand()
查看>>
SSM(Spring+SpringMvc+Mybatis)整合开发笔记
查看>>
Orderer节点启动报错解决方案:Not bootstrapping because of 3 existing channels
查看>>
org.apache.axis2.AxisFault: org.apache.axis2.databinding.ADBException: Unexpected subelement profile
查看>>
sql查询中 查询字段数据类型 int 与 String 出现问题
查看>>
org.apache.commons.beanutils.BasicDynaBean cannot be cast to ...
查看>>
org.apache.dubbo.common.serialize.SerializationException: com.alibaba.fastjson2.JSONException: not s
查看>>
sqlserver学习笔记(三)—— 为数据库添加新的用户
查看>>
org.apache.http.conn.HttpHostConnectException: Connection to refused
查看>>
org.apache.ibatis.binding.BindingException: Invalid bound statement错误一例
查看>>
org.apache.ibatis.exceptions.PersistenceException:
查看>>
org.apache.ibatis.exceptions.TooManyResultsException: Expected one result (or null) to be returned
查看>>
org.apache.ibatis.type.TypeException: Could not resolve type alias 'xxxx'异常
查看>>
org.apache.poi.hssf.util.Region
查看>>
org.apache.xmlbeans.XmlOptions.setEntityExpansionLimit(I)Lorg/apache/xmlbeans/XmlOptions;
查看>>
org.apache.zookeeper.KeeperException$ConnectionLossException: KeeperErrorCode = ConnectionLoss for /
查看>>