最会设计的科研狗

文章目录一、绘制png：MolDraw2DCairo二、绘制svg：MolDraw2DSVG三、绘制png设置(svg设置类似)一、绘制png：MolDraw2DCairo先导入所要用到的库和MolDraw2DCairo模块。from rdkit import Chemfrom rdkit.Chem.Draw import rdMolDraw2Dfrom IPython.display import Image首先，创建一个Cairo drawer：Chem.Draw.rdMolDraw2D

文章目录一、分子指纹计算二、相似性搜索三、自定义搜索函数一、分子指纹计算本文介绍在windows环境下，使用rdkit函数在postgresql数据库中进行相似性搜索。环境搭建、数据表准备不再赘述，可以参考这篇文章。在上述工作基础上，继续进行指纹计算、建立索引。操作之前先看看在postgresql中支持的指纹函数：layered_fp(mol)：另一种rdkit原创指纹，官方文档的解释是它一种子结构指纹，与rdkit拓扑子图的生成步骤一致，但根据子图生成指纹向量的过程有所不同。在子结构指纹类别中，l

文章目录一、环境搭建二、数据表准备三、结构搜索与查询1.smiles子结构搜索2.smarts子结构搜索3.立体信息的结构搜索4.带取代基的结构搜索一、环境搭建本文以windows为例进行操作和演示。假设都已经安装好了postgresql，下面依次进行用户切换、启动服务端、创建数据库、导入数据、连接数据库、加载插件、函数测试。在windows上安装带rdkit插件的postgresql以及基本操作可以参考这篇文章。在cmd中切换到postgres用户runas /user:postgres cm

文章目录一、环境配置1.windows下的安装与初始化2.用户设置与服务启动二、操作使用1.cmd操作2.python操作一、环境配置PostgreSQL是一个开源、可扩展的关系型数据库，rdkit官网文档里也是以postgresql作为案例进行演示。本文介绍带rdkit插件的postgresql在windows系统上的环境配置，以及在python环境下的基本操作。1.windows下的安装与初始化直接在cmd下使用conda安装C:\Users\Administrator> conda

文章目录一、初级篇1.氢原子显示与隐藏2.芳香键与kekule式转换二、高级篇1.Atom和Bond对象的编辑功能2.RWMol类的编辑功能一、初级篇1.氢原子显示与隐藏正常情况下，分子在rdkit中存储时，氢以隐式氢的形式存储，即不会在图片中显示出来。当需要加入氢原子时，例如要生成和优化立体结构，可以通过函数加上氢原子。加氢：AddHs()去氢：RemoveHs()>>> from rdkit import Chem>>> m = Chem.MolFr

文章目录分子基础操作与药效团查找1.原子操作2.键操作3.环操作4.手动实现氧族药效团查找分子基础操作与药效团查找1.原子操作在rdkit中，分子中的每一个原子都是对象，可以通过原子对象的属性和函数来获取各种信息。对原子进行遍历：m.GetAtoms()获取原子索引：GetIdx()获取原子序号：GetAtomicNum()获取原子符号：GetSymbol()获取原子连接数（受H是否隐藏影响）：GetDegree()获取原子总连接数（与H是否隐藏无关）：GetTotalDegree()

文章目录Recap和BRICS对分子拆解与合成1.Recap拆解2.BRICS拆解3.BRICS合成Recap和BRICS对分子拆解与合成1.Recap拆解另一个与化学反应相关的功能是Recap，Recap可以模仿实验室中的正向合成过程来进行逆向操作，对分子进行一系列的转换与分解，最终得到一组合理的分子片段。Recap可以对拆解过程进行追踪，并形成类似树的数据结构。原始分子记为根节点（root），被拆解的分子记为父节点（parent），拆解得到的分子记为子节点(children)，某节点下的所有分支节

文章目录在化学反应中保护原子1.rxn文件创建反应2.保护目标原子在化学反应中保护原子本文是化学反应的进阶操作，关于使用rdkit进行化学反应的操作可以参考这篇文章。有时在使用rxn文件时，很难准确表达或全面考虑到每个细节，导致不需要的原子发生反应，并产生副产物。而rdkit也提供了保护原子的操作，避免不需要的原子参与到反应中。还是以形成酰胺键为例。1.rxn文件创建反应通过rxn文件创建一个反应模式：ReactionFromRxnFile()检查一下该反应：ReactionToImage()

文章目录化学反应1.SMARTS创建反应2.rxn文件创建反应3.产物后处理化学反应Rdkit中提供了基于SMARTS的化学反应操作，可以通过SMARTS或rxn反应文件构建反应模式，再对指定的反应物进行匹配，将匹配上的结构按反应模式进行合成。1.SMARTS创建反应从SMARTS创建一个反应模式：Chem.ReactionFromSmarts()下面这个反应表示一个羧基，与至少带有一个氢的氮原子反应，形成酰胺定义反应物：reactants进行一个反应：RunReactants(reacta

文章目录一、骨架分解1.Murcko Scaffold2.Generic Framework二、侧链分离1.rdRGroupDecomposition2.ReplaceCore一、骨架分解1.Murcko ScaffoldMurcko骨架由Murcko等人设计并用药物的形状、结构分析。他们将药物分子拆解成四种单元：环系结构（ring system）、接头（linker）、骨架（scaffold）、侧链（side chain），其中scaffold又由ring system和linker组成；scaff

文章目录一、分子片段生成二、片段指纹生成三、指纹重要性分析一、分子片段生成分子片段（Molecular Fragments）是一组相连的原子，并可能包含有相关官能团。在rdkit中提供了一系列用于分析、操作分子片段的工具。说起来比较抽象，操作起来也比较抽象。获取官能团库：RDConfig.RDDataDir目录下的’FunctionalGroups.txt’根据官能团库实例化一个参数器：FragmentCatalog.FragCatParams()>>> import os

文章目录一、化学特征和药效团提取二、化学特征文件介绍1.化学特征（chemical features）2.FDef文件语法三、2D药效团指纹1.编码原理2.参数设置3.生成2D药效团指纹4.修改FDef设置5.Gobbi 2D药效团指纹一、化学特征和药效团提取Rdkit中有一系列基于SMARTS定义的化学特征。要对分子进行化学特征分析和计算时，需要先导入一个特征库，创建一个特征工厂，并通过特征工厂计算化学特征。先来操作一下，原理及介绍在第二部分。获取特征库：RDConfig.RDDataDir目录下

文章目录一、描述符计算模块1.rdkit.Chem.Lipinski模块2.rdkit.Chem.Descriptors模块3.rdkit.ML.Descriptors.MoleculeDescriptors模块二、原子描述符可视化1.原子partial charge可视化2.原子logP可视化一、描述符计算模块1.rdkit.Chem.Lipinski模块rdkit中提供了许多描述符的计算方法，可用于分子筛选、成药性评估等。以lipinski类药的相关规则为例，可以通过rdkit.Chem.Lipi

文章目录一、分子指纹提取1.Topological Fingerprints2.MACCS3.Atom Pairs and Topological Torsions4.Morgan Fingerprints (Circular Fingerprints)二、相似性比较一、分子指纹提取1.Topological Fingerprints也叫RDKit topological fingerprint。该算法受Daylight fingerprint启发而产生,它会计算所有介于minPath和maxPath

文章目录最大公共子结构1.搜索方法2.bondCompare参数3.atomCompare参数4.maximizeBonds参数5.matchValences参数6.ringMatchesRingOnly参数7.completeRingsOnly参数8.timeout参数最大公共子结构rdkit中提供了在多个分子间查找最大公共子结构（Maximum Common Substructure, MCS）的功能。该功能可以应用在相同骨架搜索、分子库分析、寻找药效团等场景中。1.搜索方法FindMCS(m

文章目录一、摩根分子指纹计算1.简介2.SparseIntVects3.ExplicitBitVects4.FCFPs5.更多泛化二、摩根分子指纹提取1.提取方法一2.提取方法二三、指纹可视化一、摩根分子指纹计算1.简介摩根分子指纹（Morgan Fingerprints），是一种圆形指纹，也属于拓扑型指纹，是通过对标准的摩根算法进行改造后得到。可以大致等同于扩展连通性指纹（Extended...

文章目录分子子结构操作1.删除子结构2.替换子结构3.切掉侧链4.切掉母核分子子结构操作RDKit包含了一些修改分子的函数，这些函数可以方便地对分子进行子结构删除/替换等操作。更复杂的操作可以看Chemical Reactions中相关的功能。1.删除子结构先初始化一下定义一个苯丙氨酸分子，要把苯甲基删掉，变成甘氨酸>>> from rdkit import Ch...

文章目录一、SMARTS简介二、原子属性三、键属性四、逻辑操作符五、递归SMARTS六、组合匹配七、反应SMARTS一、SMARTS简介是什么SMARTS（SMiles ARbitrary Target Specification）是一种用于描述分子模式和属性的语言。SMILES所有的符号和属性在SMARTS中同样适用，因此它也是SMILES的延伸。此外，SMARTS还包括了逻辑操作符和...

文章目录一、简单子结构搜索1. 初始化2. 获取所有匹配的结构3. 查看匹配的结构二、带条件的子结构搜索1. 定义条件2. 对侧链设置条件3. 使用条件进行筛选4. 其他条件举例一、简单子结构搜索1. 初始化m表示待搜索的原始分子p表示带侧链的四元环结构>>> from rdkit import Chem>>> from rdkit.Chem i...

文章目录一、单个子结构搜索与展示1.判断是否包含子结构2.搜索并返回子结构3.获取匹配的原子和索引4.绘制子结构二、多个子结构搜索与展示1.搜索并返回所有子结构2.获取匹配的原子和索引3.设置原子和键的颜色4.绘制多个子结构一、单个子结构搜索与展示首先初始化一个分子，并设置一个要搜索的SMARTS/SMILES>>> from rdkit import Chem>&...

文章目录一、简介二、构象生成算法三、四、一、简介让计算机识别分子结构是计算化学码农的必备技能，也是对分子进行后续操作的基础。本文整理和总结了rdkit进行读取、输出和可视化的一些方法，包含对SMILES、SDF、MOL、MOL2、CSV等文件的处理，以及分子的结构展示。二、构象生成算法先来了解两种rdkit生成3D构象的方法：距离几何法距离几何法步骤：1.通过分子的连接表信息和一套规...

读取分子：读SMILES、读DataFrame、读取sdf、mol、mol2...输出分子：输出SMILES、输出sdf、mol、pdb...可视化：单独可视化、批量可视化、3D可视化...

文章目录一、简介二、应用三、特性四、与路径指纹的对比五、表示方法1.整数标识列（list of integer identifiers）2.定长比特串（fixed-length bit string）六、生成流程1.初始化原子标识符2.标识符的迭代更新3.标识符去重七、参数设置1.主要参数2.原子属性八、功能基指纹九、特征可视化十、rdkit计算ECFPs一、简介本文参考自ChemAxon的介...