Java教程

化学分子溶解度预测模型(Rdkit构建)

本文主要是介绍化学分子溶解度预测模型(Rdkit构建),对大家解决编程问题具有一定的参考价值,需要的程序猿们随着小编来一起学习吧!

各位朋友好,今天我讲述如何用Rdkit构建化学分子的溶解度预测模型。首先我们要了解一个非常重要的包Rdkit。

v2-7b72a3113620fc081da6b223feeaa13d_720w.jpg

 

Rdkit介绍

•Rdkit是开源的化学信息工具包

•采用了进行封装,提供Python2/Python3的接口

•通过SWIG技术提供Java 和C# 接口

•提供了大量对化学分子2D/3D的计算操作

•生成用于机器学习的分子描述商业友好的BSD协议

•核心数据结构和算法用C++实现

•通过Boost.Python技术对RDKit符

•基于PostgreSQL搭建分子数据库

•KNIME中的化学信息计算支持(https://www.knime.com/rdkit)

RDKIT可以处理化学分子内容如下图

v2-7076c80e668b0e38d0971a2b53310e60_720w.jpg

watermark,size_14,text_QDUxQ1RP5Y2a5a6i,color_FFFFFF,t_100,g_se,x_10,y_10,shadow_20,type_ZmFuZ3poZW5naGVpdGk=

接下来,我们要探索如何用Rdkit构建化学分子的溶解度预测模型。

我们来看看术语‘溶解度solubility’是什么意思。

溶解度,符号S,在一定温度下,某固态物质在100g溶剂中达到饱和状态时所溶解的溶质的质量,叫做这种物质在这种溶剂中的溶解度。物质的溶解度属于物理性质。

溶解度的单位是克(或者是克/100克溶剂),而不是没有单位.

例如,在20°C的时候,100克水里溶解0.165克氢氧化钙,溶液就饱和了,氢氧化钙在20°C的溶解度就是0.165克,也可以写成0.165克/100克水。

v2-825ae88334e96139c6721e0b55b7f4e7_720w.jpg

每一种化学分子都有相应溶解度,那么我们可否根据已知化学分子溶解度来预测其它化学分子溶解度呢?

v2-a5a33d8d55cc2607cd70290fb3795ee6_720w.jpgv2-70b53534fde640d0c6133746f06b4ad2_720w.jpg

答案是可以的。我们可以用python的rdkit,sklearn包搭建分子溶解度预测模型。

v2-3d0981eef18c609426510f5ff1d38288_720w.jpg

建模过程中,我们需要解决一些难题。

难点1:化学分子如何转换为数字?用于机器学习计算。

v2-072ffdbe51bafc9ad3eda805d02a911d_720w.jpg

难点2—如何构建机器学习模型?

v2-b56b59c479def8d2b47f7c9c4284e2bc_720w.jpg

欢迎给位朋友关注<python机器学习-乳腺癌细胞挖掘>教程:https://edu.51cto.com/sd/7036f

里面有详细讲述如何用Rdkit构建化学分子的溶解度预测模型,并解决上述难题。

教程计算逻辑如下图

v2-8b9c65d0d2130058c13eeb3a6741548a_720w.jpg

教程有完整数据和调试好的python脚本

v2-69328fa0584b2d91adb5fdf9f046c29a_720w.jpg

这是化学分子式和溶解度数据

v2-2e8aaa71e4e782516ecbdd2750a316f2_720w.jpg

这是通过python脚本,把dat数据保存到Excel文件。

v2-0e460702af350e182cb7175ea7347248_720w.jpg

部分python脚本把smiles字符串形式转换为分子式,并可视化展示截图

v2-a92155c834e2c6b95fa2252babf59e7e_720w.jpg

模型校验可视化

v2-e419c2461024879c9967b95e4af1d0b1_720w.jpg

 

术语补充知识

SMILES(Simplified molecular input line entry specification),简化分子线性输入规范,是一种用ASCII字符串明确描述分子结构的规范。SMILES由Arthur Weininger和David Weininger于20世纪80年代晚期开发,并由其他人,尤其是日光化学信息系统有限公司(Daylight Chemical Information Systems Inc.),修改和扩展。

由于SMILES用一串字符来描述一个三维化学结构,它必然要将化学结构转化成一个生成树,此系统采用纵向优先遍历树算法。转化时,先要去掉氢,还要把环打开。表示时,被拆掉的键端的原子要用数字标记,支链写在小括号里。

SMILES字符串可以被大多数分子编辑软件导入并转换成二维图形或分子的三维模型。转换成二维图形可以使用Helson的“结构图生成算法”(Structure Diagram Generation algorithms)。

v2-15a8a91edadd1eba8e6d33e6fab28492_720w.jpg

 

这篇关于化学分子溶解度预测模型(Rdkit构建)的文章就介绍到这儿,希望我们推荐的文章对大家有所帮助,也希望大家多多支持为之网!