ECFPs_by_rdkit.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 58,
   "metadata": {},
   "outputs": [],
   "source": [
    "# 参考http://www.rdkit.org/docs/GettingStartedInPython.html\n",
    "from rdkit import Chem\n",
    "from rdkit import DataStructs\n",
    "from rdkit.Chem import AllChem\n",
    "from rdkit.Chem import Draw"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# 分子指纹生成"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "* 初始化"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 17,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAcIAAACWCAIAAADCEh9HAAAO40lEQVR4nO3da2wU5R4G8Ge2paVyL6QeQYMo9IYgHIO3D6wGFAwgB2ugJiqCEixqbYsg2iY9KjUBJC2IShCj+AECKglgRJAQmoZgiCAgEluUFG25VEoppSyl7L7nw3rW2t1ud+d9553Z9vmFNLM7szN/+Hce5rYzhhACRERklsvuAoiIYhtjlIhICmOUiEgKY5SISApjlIhICmOUiEgKY5SISApjlIhICmOUiEgKY5SISApjlIhICmOUiEgKY5SISApjlIhICmOUiEgKY5SISApjlIhICmOUiEgKY5SISApjlIhICmOUiEgKY5SISApjlIhICmOUiEgKY5SISApjlIhICmOUiEgKY5SISApjlIhICmOUiEgKY5SISApjlIhICmOUiEhKvLmPGcZfA0L8PRzQ9k0hzC2BrMUOEqlicmvUv2q1+9l2fWs7ihyIHSRSxeTWaISuXUNtrcnPulz1Pt8lc58dOnRofLy1f7Vugh0k6pwwq+1HA8P+gcDP8nIBmPzjdheY/kuVlpZWV1eb/qt1E+wgkRJS/+EHH1NrJykJd9xhcuYDBw66w9SHGxoa8vPzDxw4sHnzZpPL7jbYQSJ5hjB79Msw/j5wFhj2D7T9qV9tbW1aWlpzc/O+ffvcbrcNFcQIdpBICUsueAp58lebIUOGLF68GEBeXp7X67WtjljGDhJFzpKtUfx/b9Gu87wejyczM7O6unrdunXz5s2zpwjHYweJlDAZo51edQjYtkvot2XLllmzZqWkpFRWVvbv39+2OpyKHSRSxfx1o4HLDIXA8ePo1QuFhf+49tDeSw5nzpzpdrvr6uqWLl1qZx1OxQ4SqaLm2OjChWhuRmOjkpkpU1ZWFhcXt3r16srKSrtrcTp2kMg0BTG6bRt27cKAASgulp+ZSmPGjJk7d25ra+vChQvtrsXR2EEiGeZPMfldv45Ro1BVhfffx8svq6pKmbq6urS0tEuXLn3zzTePPfaY3eU4ETtIJEl2a7S0FFVVyMjA/PlK6lEsJSWlqKgIQEFBQWtrq93lOBE7SCRL5itQ586Jfv0EIL79VsU3qqzR0tKSlpYGoLS01O5aHIcdJJInFaNz5ghATJ+uqBbL7NixA0D//v3r6ursrsVZ2EEieeZj9PBh4XKJhARRWamwHqv4D6u9+OKLdhfiIOwgkRImTzEJAbcbFRVYvBjLlik8xmCVX375ZfTo0T6f79ChQ3fffbfd5diPHSRSxeQppk2bUFGBlBS8+abaeqySnp6+YMECr9ebn59vdy2OwA4SqWJma9TjQUYGTp/GJ59g7lwrqrJEQ0NDamrqhQsXvvzyy6ysLLvLsRM7SKSQma3RZctw+jTGjsVzz6kux0oDBgx4++23ARQUFFy9etXucuzEDhKpFO3B1D/+EL16/XVf9Jhz48YN/2G1pUuX2l2LbdhBIrWijtFZswQgnnrKimJ02Lt3L4Cbbrrp999/t7sWe7CDRGpFF6P79wvDEElJIqYfk/PEE08AePbZZ+0uxAbsIJFyUcSo1+udMuUyIIqLLStHi1OnTvXs2dMwjIqKCrtr0YodJLJCFKeYPv300927/zV9+s7Fi604SKvPsGHDCgoKhBB5eXk+n8/ucvRhB4ksEWHcXr58+ZZbbgGwadMmC1Ndl6ampsGDBwPYsGGD3bVowg4SWSTSGH3ttdcAPPjggz6fz9KCtNmwYQOAm2++ubGx0e5adGAHiSwSUYz++uuviYmJLpfr4MGDVhekjc/nu++++wAUFhbaXYvl2EEi60QUo1OnTgXwwgsvWF2NZj/88IPL5UpMTDx58qTdtViLHSSyTucxumfPHgB9+vQ5c+aMhoI0e+aZZwBkZWXZXYiF2EEiS3USo62trXfddReAFStW6ClIs7Nnz/bt2xfA7t277a7FEuwgkdU6idGysjIAd95557Vr1/QUpJ//+b0jR45sbW21uxb12EEiq4WL0fr6+oEDBwLYvn27toL0u3bt2vDhwwF89NFHdteiGDtIpEG4GM3JyQEwYcIEbdXY5YsvvgCQnJx84cIFu2tRiR0k0qDDGD1+/Hh8fHx8fPxPP/2ksyC7PPLIIwDy8vLsLkQZdpBIjw5j1P9L+eqrr2osxk5dL3TYQSI9Qt/9fuvWrVlZWcnJyVVVVf6Da91BTk7O2rVrJ06c+N1332letHHYCAyLf4u27/hfRosd1LbQt6YNbvuyeMcZbYsmpwhO1paWlhEjRgD48MMPtce6nerr65OTkwF8/fXXOpeLQwg5HPwyQuyg5g7+d+otIYepmwhxh6f33nvv5MmTmZmZ8+bN0xvpNktOTi4qKgKQl5fX0tKiZ6HGYcPc9mYY7KDODhK136kXQrjd7oqKij179kyYMMGusuzS2tqampp69uxZEyuhO9ddPrvcxELDxKiJkGUHTXfw6f88MvzG8Wg/NXF24Z4NJYGX/p36wJ5+8Y4z/mHu7Hdh8XYXoJsQ4sqVK3369LG7EDLJmR0MpORb0wYX7zgTCFP/S1tLI+sF7+e/++67ADIzM7vkd0I2bdqUkpKycePGkGNXrlwJYPjw4dq+8xP+6Ke5Y6PsoM4Oio6PjfI4aTcR4tjowoULR4wYceLEiXXr1umNdMt5PJ4lS5bU1dWFfDzvxYsXS0pKAKxatSoxMVF7dcqwg7HeQYoxIcN169at6IrfCSkuLgYwduxYr9cbPHb+/PkAJk6cqLkq5WfqBTuoF7dGu7nQ140CmDRp0u7du3Nzc1etWqUt0y1VU1OTnp7e3NxcXl4+fvz4dmOPHj16zz33GIZx5MiRkSNHaq4t+CrRtleSwtTVo+ygHiGvG+XFpN1LR/n6888/+78TcuzYMY2xbqHs7GwA2dnZIcc+9NBDAPLz8zVXZR12kEiPcPuMCxYsQFe5scX+/fsNw0hKSqoO9YD2LVu2AEhOTq6vr9dfm3XYQSINIrpR3rZt27QVZAWv1ztu3DgAxaEe0O7xeG6//XYAa9eu1V6atdhBIg06OYPhP6wW6zf9/fjjjwHceuutV65cCR77zjvvoOve9JcdJLJa5w8RGTVqFIDly5frKUi58M9nr6mp6d27N7ruIyjYQSKrdf1H2i1atAjAAw88EPL57E8//TSAJ598Un9h2rCDRJaK6LLEadOmAXj++eetrka58M9nP3DggGEY3eHxvOwgkXUiitHwv8pOFiY+fD7fvffeC6CoqEh/YZqxg0TWifRLMuF3rJwp/M7sZ599BmDIkCFNTU36a9OPHSSySKQxGjjM39EtIZwm8Hz2kKdWmpqaBg8eDODzzz/XX5st2EEii0Txle3169f36NFj6fTpornZuoJUCX+hzxtvvAHg/vvvj6FNM3nsIJEVoohRr9dbNWWKAESoS6AdJfxl57/99lvPnj0Nw/j+++/112YjdpDIClHeQGj/fmEYIilJhPpCnnOE/xLkjBkzAMyePVtvUc7ADhKpFv192LKzBSA6uD2EE4S/JcfevXsB9O7du7a2Vn9tjsAOEikVfYz+8Yfo1UsAorzcgnoUePTRRwHk5uYGj7px48bo0aMBlJSU6C/MKdhBIqVM3RW4uFgAYuxYEerWufby3654wIABIW9X/MEHHwAYNmyYx+PRX5uDsINE6piK0atXxdChAhDr16uuR0rg+exr1qwJHnvx4sVBgwYB+Oqrr/TX5izsIJE6Jp9RITZuFIBISRGXLimtR0r4R7nl5uYCePjhh/UX5kTsIJEiZmNUCDF+vADEokXqipFy7ty5vn37Ati1a1fw2BMnTvTo0SMuLu7o0aP6a3ModpBIBYkYPXxYuFwiIUFUVqqrx7zZs2cDmDFjRsixkydPBpCTk6O5KkdjB4lUkIhRIcTcuQIQjz+uqBjzDh065HK5EhISqqqqgsdu377df9bizz//1F+bo7GDRNLkYvT8edGvnwDEzp2K6jFp0qRJAF5//fXgUS0tLampqQDKysr0F+Z07CCRNLkYFUIsWyYAkZEhrl9XUY9J58+fz8/Pv3z5cvCoFStWAMjIyLhua4XOxQ4SyenwOfWRun4do0ahqgqrV+OVV6RmZYG6urrU1NTGxsadO3f6D65Re+wgkRyX7AwSErBiBQAUF+PCBfmC1CosLGxsbJw6dSrXwA6xg0RypLdG/SZPxq5deOklrFmjYG6K/Pjjj+PGjYuLizt27FhaWprd5TgbO0hklvTWqN/KlejVC/37q5mbInl5eV6vNzc3l2tg59hBIrMUbY0CuHTJUSvh5s2bs7OzU1JSqqqq+vXrZ3c5sYAdJDJF0dYo4Kg10OPxLFmyBEBJSQnXwEixg0SmqIvRAMOAYaifbTSWL19eXV09ZsyYOXPm2FtJTGIHiaKhbqf+r/kZ8M8wMKBdTU1Nenp6c3Pzvn373G63LTXEMHaQKErxKmfWdsUTAoaBgweRnW1uZt+MGfPKkSMmPtjQ0NDc3Dxz5kyugVFjB4mipzRGg3k8OHXK3Ee9t912yuxnS0tL/Y/rIVnsIFFnlO7Ut9sNNAx4PKitNTeziy5Xg89n7rNDhw6Nj7f4f4guiR0kip7FMWrTwTUyiR0kip4FZ+qJiLoTpTHqPynhxw2ZWMQOEkVP9QVPwF/rIdfA2MUOEkXDghglIupOeGyUiEgKY5SISApjlIhICmOUiEgKY5SISApjlIhICmOUiEgKY5SISApjlIhICmOUiEgKY5SISApjlIhICmOUiEgKY5SISApjlIhICmOUiEgKY5SISApjlIhICmOUiEhKvN0F/M0IPJMSEEIEXrZ9WlS7aQJvhpnG/7LdDIPnE7y4SJZFASH/uSKcjB2kmOaUGG37yx1YbfDPNSF4mnYTdDSftgOBn22nafdOyGmCl0Vttf1H6yiq2EHqkpy4Ux/JxkJgmjATd7SBY8WyKFrsIHUZTtkapa4heE+i7bakvbURWcSJW6NWaHc0jSwihGiXm/6da/kMZQfJsbpLjJJduBFKXZ4TYzSSjQ4T03R00kPJsiha7CB1GQ66+KPdBSshV6EwF7V09E67q2RMzyfkOxQQyRVIId9kBynWOShGiYhikRN36omIYghjlIhICmOUiEjK/wBw57lLwgWdsgAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<PIL.PngImagePlugin.PngImageFile image mode=RGB size=450x150 at 0x7AB9848>"
      ]
     },
     "execution_count": 17,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# 初始化rdkit mol对象m1和m2，并查看\n",
    "m1 = Chem.MolFromSmiles('ClC1=COCNC1')\n",
    "m2 = Chem.MolFromSmiles('BrC1=COCNC1')\n",
    "Draw.MolsToGridImage([m1, m2], subImgSize=(150, 150), legends=['ClC1CNCOC1', 'BrC1CNCOC1'])"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "* 以SparseBitVects方式生成ECFPs指纹，需要传入mol和radius"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 18,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "0.7"
      ]
     },
     "execution_count": 18,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "fp1 = AllChem.GetMorganFingerprint(m1, 2)\n",
    "fp2 = AllChem.GetMorganFingerprint(m2, 2)\n",
    "DataStructs.DiceSimilarity(fp1, fp2)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "* 以ExplicitBitVects方式生成ECFPs指纹，需要传入mol，radius和长度nBits"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 19,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "0.6842105263157895"
      ]
     },
     "execution_count": 19,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "fp1 = AllChem.GetMorganFingerprintAsBitVect(m1, 2, nBits=1024)\n",
    "fp2 = AllChem.GetMorganFingerprintAsBitVect(m2, 2, nBits=1024)\n",
    "DataStructs.DiceSimilarity(fp1, fp2)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "* 通过设置useFeatures=True生成FCFPs <br /> \n",
    "可以看到ffp1和ffp2完全一致，即Cl和Br在特征上被归为一类"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 20,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "1.0"
      ]
     },
     "execution_count": 20,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "ffp1 = AllChem.GetMorganFingerprint(m1, 2, useFeatures=True)\n",
    "ffp2 = AllChem.GetMorganFingerprint(m2, 2, useFeatures=True)\n",
    "DataStructs.DiceSimilarity(ffp1, ffp2)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 21,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "1.0"
      ]
     },
     "execution_count": 21,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "ffp1 = AllChem.GetMorganFingerprintAsBitVect(m1, 2, useFeatures=True)\n",
    "ffp2 = AllChem.GetMorganFingerprintAsBitVect(m2, 2, useFeatures=True)\n",
    "DataStructs.DiceSimilarity(ffp1, ffp2)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# 分子指纹提取"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "* 通过向参数bitInfo传入字典，可以获取所有非空指纹信息 <br /> \n",
    "字典info键表示位的索引，值为原子索引、半径构成的元组，出现多个元组时，表示子结构出现了多次<br />\n",
    "键为2968968094的值中含有多个元组，表示原子索引为4,6，半径为0的子结构相同，都记录在这一位上"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 24,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "{39328034: ((1, 1),),\n",
       " 211414882: ((5, 2),),\n",
       " 362715007: ((6, 1),),\n",
       " 397705891: ((4, 1),),\n",
       " 718785834: ((1, 2),),\n",
       " 1016841875: ((0, 0),),\n",
       " 1078999752: ((3, 1),),\n",
       " 1289643292: ((5, 1),),\n",
       " 2132511834: ((5, 0),),\n",
       " 2626911012: ((4, 2),),\n",
       " 2968968094: ((4, 0), (6, 0)),\n",
       " 3145053355: ((2, 1),),\n",
       " 3189457552: ((3, 0),),\n",
       " 3217380708: ((1, 0),),\n",
       " 3218693969: ((2, 0),),\n",
       " 3452535345: ((0, 1),),\n",
       " 3916672621: ((6, 2),),\n",
       " 4196496269: ((2, 2),),\n",
       " 4217154750: ((3, 2),)}"
      ]
     },
     "execution_count": 24,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "info = {}\n",
    "fp_explain = AllChem.GetMorganFingerprint(m1, 2, bitInfo=info)\n",
    "info"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## 方法一\n",
    "* 以提取出211414882的结构为例，首先提取m1中，半径为2，第5个原子的结构 <br />\n",
    "submol接收提取出的子结构 <br />\n",
    "amap用于接收原子索引的映射关系，键为原始分子中的原子索引，值为子结构中的原子索引<br />\n",
    "env是被提取出的键的索引"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 25,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "{1: 4, 3: 0, 4: 1, 5: 2, 6: 3}\n",
      "[4, 5, 3, 6]\n"
     ]
    }
   ],
   "source": [
    "amap = {}\n",
    "env = Chem.FindAtomEnvironmentOfRadiusN(m1, 2, 5)\n",
    "submol=Chem.PathToSubmol(m1, env, atomMap=amap)\n",
    "print(amap)\n",
    "print(list(env))"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "* 用SMILES表示该子结构"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 26,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "'CCNCO'"
      ]
     },
     "execution_count": 26,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "Chem.MolToSmiles(submol)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "* 以第5个原子为中心，用SMILES表示该子结构，可以看出是氮原子"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 28,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "'N(CO)CC'"
      ]
     },
     "execution_count": 28,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "Chem.MolToSmiles(submol, rootedAtAtom=amap[5], canonical=False)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## 方法二\n",
    "* 遍历env中的所有键，获取键连接的所有原子，保存在atoms中 <br />\n",
    "通过给定分子，及要提取的原子信息、键信息，获取子结构"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 29,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "'N(CC)CO'"
      ]
     },
     "execution_count": 29,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "atoms=set()\n",
    "for bidx in env: \n",
    "    atoms.add(m1.GetBondWithIdx(bidx).GetBeginAtomIdx())\n",
    "    atoms.add(m1.GetBondWithIdx(bidx).GetEndAtomIdx())\n",
    "Chem.MolFragmentToSmiles(m1, atomsToUse=list(atoms), bondsToUse=env, rootedAtAtom=5)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# 分子指纹可视化"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "* 绘制211414882代表的子结构"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 68,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAJYAAACWCAIAAACzY+a1AAAP9UlEQVR4nO2da2wTVxbHz9gev+LYTmKSYB5JKCGUR2EJAspDW9ryoUV9ZpetupB2VYXSlqoNXbXbdlGaQrVdCRqkqqBoVaksVFu6VLRF6gfEUrVQWCHCRrukEWHz4pU4fr899szc/TC7Xnc8Nk58x/aN5vcB4bkz9x77nztz7uOcoRBCoEAyqmIboJAvioTEo0hIPIqExKNISDyKhMSjSEg8ioTEo0hIPIqExKNISDyKhMSjSEg8ioTEo0hIPIqExKNISDyKhMSjSEg8ioTEo0hIPIqExKNISDyKhMSjSEg8ioTEo0hIPIqExKMpTDOxGB8Oc2o1ZTKpNRqqMI0KIIQYhuF5nqZpmqYL2XRhkFfCYJD717/CIyOxUIgTjlAU2Gz0ggWGpiajrFoihDwej9vtDoVCyegtmqYtFkt1dbXBYJCv6QJDyRSchhD84x+hy5dDPC9dv8mkvu8+q92ulaP1SCQyPDwci8UynVBdXT179myKKuj9QCZkkZDn4cwZ79BQxl9QQKWCn//c2tiIuUMEAoHBwUGe57OfZjKZGhsbVSrivQFZvsDf/x64o34AwPPw3Xe+sbE4xqZjsVgu+gFAKBQaGRnB2HSxwC/h+Hj8ypVwjicLKnIctjvByMhILvoJeL1er9eLq+ligV/Cy5dDkzo/EOCuXYtiaToQCITDuf71CNy+fRtL00UEp0f6P+egsrt77PnnZ4pKUw92d4+lFv3739GFC435G+DxeCZ7SSwWi0QiRiOG1osFzl4oOEaCPMl/U9VKLUrF4UhgcapCocndAPK5qnQo0NBexKefWv75T13qkd27IRcPf8mS31y5cipTaSKRyFT0xhtvLFmyxGazpRfF4zj9qcJTIAmFu2iy/4VCKp9PnXqCz5dTPbNmuaf29Nq9e/e6devef//9KVxb4uCXMP0pmM7Wrf5f/SqQ/KhSwZYtM9TqO3dDleoTns/o+/T390t2RLfbvX379tOnT1++fHnFihWiUtJn3fBLmOxqWbQsK/uJ32+z0XPn5jhRUpmlLBaL+aS6c3V1dWtra3d39759+44ePSoazhPty4Acg4qqKuk/akk3VWDePD2WpisqKjIVtba22u32gYGBr776KvU4TdMmkwlL68UCv4TLl5dN6nytVnX33Xj6QUVFhVYrPemq0+lefvllADh06FAwGEwer6mpIX2mFKeEwk8xf75hzhyd0OGef35m8j8gNZwAgHvvLdfp8JhBUVRdXV2m0k2bNq1YscLj8Xz88cfCEaPRWF1djaXpIiLLNDfD8CdPum/c4HfvnhEOq3bu9CxdykieuXRp2b33mvG27nA4bt68KVk0MDCwdetWlUp17NixxsbGpqamTL2WIGSZ5tbpVI88UnX6tDUcVt19NyOpn0pFrVljxq4fANTU1NTX10suQSxYsOCxxx5jWfbDDz+cHvqBfOuFfX2wfDkgBF1dPp3uJ8MAjYZqaND/7Gcmq1XGUWk8Hh8fH/d4PBzHpR6PxWKbN2/2+/3ffPPNQw89JJ8BhQPJw333IQDU3o4QQuEwNzoa6+8PX70auX2bSSR4mRpNh+f5UCjkcrmcTqfX62UYBiG0b98+AFi4cGE8Hi+YJfIhi4R//SsCQJWVyO2Wo/p8YRimqakJALq6uoptCwbwSxiNooYGBIAOHcJeNzZOnjwJAFardWJioti25At+CffsQQBo8WKUSGCvGyfCg3DHjh3FNiRfMEt48yYymRAAOnUKb8X46e/vp2larVb39vYW25a8wCzh1q0IAP3iF3hrlYtXXnkFADZu3FhsQ/ICp4QXLiCKQjodunYNY60y4vF4hBXE48ePF9uWqYNNQp5Hq1cjAPT227iqLAQHDx4EgLlz54bD4WLbMkWwSfjJJwgA1dYivx9XlYWAZdlly5YBwN69e4ttyxTBI2EwiOx2BID+/Gcs9RWUM2fOAIDRaLx+/XqxbZkKeCR8800EgNasQXzhJl5w8uSTTwJAa2trsQ2ZChgkHBoaevDB74xGdOFC/pUVh6GhIb1eT1HU2bNni23LpMEgYUtLCwDs2EGUG5PGW2+9BQDNzc0cxxXblsmRr4TCg8RkMt26dQuLQcUiGAza7XYAOHz4cLFtmRx5Sciy7D333AMA7733Hi6Disjhw4cBoKamxk+UV52XhMKgqqGhIRqN4jKoiPA8v3r1agB4m6ix7dQlTE5tfPHFFxgNKi4XLlygKEqn010jZYYpdwljsVggEPD7/ZFIhOd5NF0mGNPZtm0bALS0tCCEeJ6PRqN+v9/v90ejUb4kx0x32HiRSCQcDofH40ndJa1SqVwu1+bNm3me7+npEWY3pg3j4+NNTU2BQODIkSNLly5lWTZZpFarKyoqamtrdTpdlhoKTLbtTy6X68qVKw6HQ7TLnef5d955J5FIPP3004I7M52ora198cUXAaCjo4NhfrJxi+M4l8vV19c3Pj5eJOskyCjh7du3R0dHJSNmz549e/78+fLy8meeeWZwcDB7PyaO0dHRRx99dM6cOUNDQ19++WX6CQihW7dulU6Qt7SEbrd7bExi2y4AJBKJAwcOAEBbW5vVavX7/Tdu3JDPvgIzPj7ucrm0Wu3OnTsB4KOPPvL7/ZJnZvmJCoyEhCzLZlHls88+Gx0dra+v37Jli3DE6XSSHmUpwDBMMvLtgQceWL16dSAQSG79TmdsbEx0py0KEhI6nU7R3sskyd3su3bt0mj+vwu0pJ4NU8bhcKQ+FF577TW1Wn3s2LHBwUHJ8xFCBfjinZ2d2U+QkDBLxPrBgwdDodCGDRvWrl2bejwQCGRSnRQQQqLsF/PmzXv88cc5jtu/f3+mq7xeb9FdAbGEHMdlSpo0MDDw9ddfazSaV199VVSEEJpsqolSg2GY1PGDwEsvvWQ2my9evHju3DnJq7L8XHnS2dmZ2v9EH1MRb4nPEq7+ww8/IISeeOIJyeihLBcSgWTEvdlsbmtr279/f29vb/LZL8Ln80nGpeYDRVEdHR3wv7toZ2en8FGS4qRLKCl4nvd6vZkyeT377LOrVq3asGFDga3KHfHsDMdxvb29kqemRnald8TGxkazGX+YUgEQZg1pmh4dHRUVqVQqIYbU4/FkumEuWrRIjryKQv/r6OgQumDyY/qZ4l6oVqv1er2kuUJk14kTJ7q6uoShYRKKosrKJhfcWyKwLBuJRADAYrFoNBrR49BsNlMUxTBMJv2En0sOw0RqZbmRStw9KiszpiR44YUXysvLz507d/78+dTjZrNZrVZnuqqUCQQCAFBWVkbTtChUX6PRGAwGhJBwjiQVFRVFj/OWkHDGjBmZ9KisrHzuuecA4IMPPkj9g62trZXJPlmJxWIMw6hUKiFjgiju3mKxUBQVDoczeWoURRXgi2fpfwISEmo0mjlz5mS64KmnnqqrqxsZGfn888+FIzNmzCAxaQRCSMibUF5eLvgyOp1O2HsBAAaDQafTCeGJmWqYOXNmKSxZSLthVVVVM2dKJxjRaDTt7e0A0N3d7fF4LBZLFr1LmXA4zLKsRqNJzTtTW1trs9koihJcs2AwmCk1ZpafqMBkXKmw2+11dXWSrvb69evXrl0bDoePHDly1113Ff1hMAWS3Svdi66rq7Pb7Wq1OpFICJ6OCIqiZs2aVV9fXwA7c2GKS75ut/vhhx/mOO7ixYvNzc3y24kZn88XjUb1en16tiGO45xOJ0KIZVmPx1P6S765pktgGEaITKdpWtg1297efuDAgfXr13///fdkdcREIuFyuSiKstlsqZP1AqnqIoSELw4AWq1Wp9OV4DedesYLn8+3YMECp9N5/PhxYTcwKbjd7ng8bjKZysvLRUVJdbO45aXG1PPOWK3Wd999FwB27dol+cwoTaLRaDweTw4kRAgLvGVlZaToB3mmDmpra1u2bNn169e7urpwGSQryYGEMO0iKo1Go4lEQq1WkzVGyjd10Lfffnv//fcbjcb+/v65c+fiMksmgsFgKBSiaTo9PTBCSFjrtlqtZL1LJt8EXhs3bmxpaYlEIl1/+AMWg+SD4zhhUVNyOj4UCnEcR9M0WfoBlgRew8PDf9m+fdf58/q//Q3WrMFilhx4vd5YLGYwGKxWq6goOZCw2WzEJQnGkEavoaHhrVWr9JEI7NwJOb/mo8DE4/FYLEZRVLoXCgCBQAAhZDAYiNMPsGVCfPNNsNuhpweOHsVTIW6E1QaTyZTuamZXt/TBJKHJBMKz8PXXIfPSTLGIRCKCqym5qJlFXSLAl4902zZYtw4cDvjjH7HViQOe57MMJLKrSwRY85H29MCqVaDRQF8fzJ+Prdr8EF7kpNVqq6qqREU8zzudTp7nKyoqZFp8LwBYswI3N8Ovfw3xOPzudzirzYM7DiR4ntdqteTqB/izAo+PQ1MTBAJw6hRs2oSz5qkSj8eFGVHRcZZlnU4nAJA4kEgFd27u2lp4/XUAgPZ2SNtZWxS0Wq3khJngxRiNRqL1A1nSq//2tzB/PvT1wZ/+hL9yTCR3zRA6kEhFBgl1OhDebvX734Pbjb/+vElOdptMJuVdvhloaYFNm8Djgb17ASCRSPj9fpfLJbyiPPcXtcpEJBIRds2QO5BIRa6XHPz3LQcAIydOuH+6TUilUgnbF4riB/I8L+TjrqysLKn9E1NGrtsIt3Ch/5e/BJat3LNHVMTzvNvt/vHHHx0Oh0ytZyEYDCKE9Hr99NAPZJKQZdmrV6+OtLWxFov54kXL2bPp5yCEbt68menNPDIhbEojdzpUElkim4aHh6PRKJjNY21txoGByKJFmc50OBx6vV7y/axyEI1GAcBoNKbveiIX/M9Cj8czPDyc+/lqtXrJkiUF+02F3WkluBFtyuC/kUo+4ZpXrmxeuVLyfCGZC3YzRCRDZA0Gw3TSD7BLyDBM+m625pUrey5d6rl0KZOKoiB3hUmB+faVHnEv6Cf8X1Cx59Il0YBaiHPH3jkoitqzZw+kBAdlCbQkF8wS5hhxb7FYRLuMhBlnvNA0nXvEOrlgljDHnsTzvGiORqPRYO+FhK7CTxbMEub4Xk4hPWTyo0qlWr58uRxeRu4R6+SCWcKpzTqWlZXJ5CXmHrFOLpg9UpqmRevjqY5oqmuTSpbofoU7gn9oHwqFrl69KjooqCipn06nW7x48TQbqxUSWVYqbty4MTExkVPzFNXY2DidZiwLjyzT3LNnz06Pnk2Hoqi6ujpFvzyRbb0QYGxsbGxsLFP9Wq22vr5e0S9/ZJQQABiGmZiY8Pl8ySR1FEUZjcbKykqbzTYNNj2UAvJKmIRl2UQiQVGUVqtVlMNLgSRUkA+lQxCPIiHxKBISjyIh8SgSEo8iIfEoEhKPIiHxKBISjyIh8SgSEo8iIfEoEhKPIiHxKBISjyIh8SgSEo8iIfEoEhKPIiHxKBISjyIh8SgSEo8iIfEoEhKPIiHxKBISjyIh8fwHCdKJLmRUeUUAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<PIL.PngImagePlugin.PngImageFile image mode=RGB size=150x150 at 0x8C881C8>"
      ]
     },
     "execution_count": 68,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "Draw.DrawMorganBit(m1, 211414882, info)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "* 在ECFPs中的颜色的含义：<br />\n",
    "蓝色：表示中心原子<br />\n",
    "黄色：表示带有芳香性质的原子<br />\n",
    "灰色：表示在脂肪环中的原子"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 70,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/svg+xml": [
       "<svg baseProfile=\"full\" height=\"150px\" version=\"1.1\" viewBox=\"0 0 150 150\" width=\"150px\" xml:space=\"preserve\" xmlns=\"http://www.w3.org/2000/svg\" xmlns:rdkit=\"http://www.rdkit.org/xml\" xmlns:xlink=\"http://www.w3.org/1999/xlink\">\n",
       "<!-- END OF HEADER -->\n",
       "<rect height=\"150\" style=\"opacity:1.0;fill:#FFFFFF;stroke:none\" width=\"150\" x=\"0\" y=\"0\"> </rect>\n",
       "<ellipse cx=\"55.9881\" cy=\"92.3437\" rx=\"5.77835\" ry=\"5.77835\" style=\"fill:#E5E533;fill-rule:evenodd;stroke:#E5E533;stroke-width:2px;stroke-linecap:butt;stroke-linejoin:miter;stroke-opacity:1\"/>\n",
       "<ellipse cx=\"35.4002\" cy=\"46.7329\" rx=\"5.77835\" ry=\"5.77835\" style=\"fill:#E5E533;fill-rule:evenodd;stroke:#E5E533;stroke-width:2px;stroke-linecap:butt;stroke-linejoin:miter;stroke-opacity:1\"/>\n",
       "<ellipse cx=\"58.8609\" cy=\"63.5951\" rx=\"5.77835\" ry=\"5.77835\" style=\"fill:#E5E533;fill-rule:evenodd;stroke:#E5E533;stroke-width:2px;stroke-linecap:butt;stroke-linejoin:miter;stroke-opacity:1\"/>\n",
       "<ellipse cx=\"85.1943\" cy=\"51.7087\" rx=\"5.77835\" ry=\"5.77835\" style=\"fill:#9999E5;fill-rule:evenodd;stroke:#9999E5;stroke-width:2px;stroke-linecap:butt;stroke-linejoin:miter;stroke-opacity:1\"/>\n",
       "<ellipse cx=\"108.655\" cy=\"68.5709\" rx=\"5.77835\" ry=\"5.77835\" style=\"fill:#CCCCCC;fill-rule:evenodd;stroke:#CCCCCC;stroke-width:2px;stroke-linecap:butt;stroke-linejoin:miter;stroke-opacity:1\"/>\n",
       "<ellipse cx=\"137.403\" cy=\"71.4437\" rx=\"5.77835\" ry=\"5.77835\" style=\"fill:#CCCCCC;fill-rule:evenodd;stroke:#CCCCCC;stroke-width:2px;stroke-linecap:butt;stroke-linejoin:miter;stroke-opacity:1\"/>\n",
       "<ellipse cx=\"120.541\" cy=\"94.9043\" rx=\"5.77835\" ry=\"5.77835\" style=\"fill:#CCCCCC;fill-rule:evenodd;stroke:#CCCCCC;stroke-width:2px;stroke-linecap:butt;stroke-linejoin:miter;stroke-opacity:1\"/>\n",
       "<path class=\"bond-0\" d=\"M 54.7995,89.7104 L 30.7148,100.582\" style=\"fill:none;fill-rule:evenodd;stroke:#E5E5E5;stroke-width:2px;stroke-linecap:butt;stroke-linejoin:miter;stroke-opacity:1\"/>\n",
       "<path class=\"bond-0\" d=\"M 57.1767,94.977 L 33.092,105.848\" style=\"fill:none;fill-rule:evenodd;stroke:#E5E5E5;stroke-width:2px;stroke-linecap:butt;stroke-linejoin:miter;stroke-opacity:1;stroke-dasharray:6,6\"/>\n",
       "<path class=\"bond-7\" d=\"M 55.9881,92.3437 L 58.8609,63.5951\" style=\"fill:none;fill-rule:evenodd;stroke:#000000;stroke-width:2px;stroke-linecap:butt;stroke-linejoin:miter;stroke-opacity:1\"/>\n",
       "<path class=\"bond-7\" d=\"M 50.6693,87.4569 L 52.6802,67.3329\" style=\"fill:none;fill-rule:evenodd;stroke:#000000;stroke-width:2px;stroke-linecap:butt;stroke-linejoin:miter;stroke-opacity:1;stroke-dasharray:6,6\"/>\n",
       "<path class=\"bond-1\" d=\"M 12.5042,60.2377 L 36.5889,49.3663\" style=\"fill:none;fill-rule:evenodd;stroke:#E5E5E5;stroke-width:2px;stroke-linecap:butt;stroke-linejoin:miter;stroke-opacity:1\"/>\n",
       "<path class=\"bond-1\" d=\"M 10.1269,54.971 L 34.2116,44.0996\" style=\"fill:none;fill-rule:evenodd;stroke:#E5E5E5;stroke-width:2px;stroke-linecap:butt;stroke-linejoin:miter;stroke-opacity:1;stroke-dasharray:6,6\"/>\n",
       "<path class=\"bond-2\" d=\"M 35.4002,46.7329 L 58.8609,63.5951\" style=\"fill:none;fill-rule:evenodd;stroke:#000000;stroke-width:2px;stroke-linecap:butt;stroke-linejoin:miter;stroke-opacity:1\"/>\n",
       "<path class=\"bond-2\" d=\"M 35.5469,53.9544 L 51.9693,65.7579\" style=\"fill:none;fill-rule:evenodd;stroke:#000000;stroke-width:2px;stroke-linecap:butt;stroke-linejoin:miter;stroke-opacity:1;stroke-dasharray:6,6\"/>\n",
       "<path class=\"bond-3\" d=\"M 58.8609,63.5951 L 85.1943,51.7087\" style=\"fill:none;fill-rule:evenodd;stroke:#000000;stroke-width:2px;stroke-linecap:butt;stroke-linejoin:miter;stroke-opacity:1\"/>\n",
       "<path class=\"bond-4\" d=\"M 85.1943,51.7087 L 108.655,68.5709\" style=\"fill:none;fill-rule:evenodd;stroke:#000000;stroke-width:2px;stroke-linecap:butt;stroke-linejoin:miter;stroke-opacity:1\"/>\n",
       "<path class=\"bond-5\" d=\"M 108.655,68.5709 L 137.403,71.4437\" style=\"fill:none;fill-rule:evenodd;stroke:#000000;stroke-width:2px;stroke-linecap:butt;stroke-linejoin:miter;stroke-opacity:1\"/>\n",
       "<path class=\"bond-8\" d=\"M 108.655,68.5709 L 120.541,94.9043\" style=\"fill:none;fill-rule:evenodd;stroke:#000000;stroke-width:2px;stroke-linecap:butt;stroke-linejoin:miter;stroke-opacity:1\"/>\n",
       "<path class=\"bond-6\" d=\"M 137.403,71.4437 L 120.541,94.9043\" style=\"fill:none;fill-rule:evenodd;stroke:#E5E5E5;stroke-width:2px;stroke-linecap:butt;stroke-linejoin:miter;stroke-opacity:1\"/>\n",
       "<text dominant-baseline=\"central\" style=\"font-size:9px;font-style:normal;font-weight:normal;fill-opacity:1;stroke:none;font-family:sans-serif;fill:#7F7F7F\" text-anchor=\"end\" x=\"30.9403\" y=\"105.675\"><tspan>*</tspan></text>\n",
       "<text dominant-baseline=\"central\" style=\"font-size:9px;font-style:normal;font-weight:normal;fill-opacity:1;stroke:none;font-family:sans-serif;fill:#7F7F7F\" text-anchor=\"end\" x=\"10.3525\" y=\"60.0639\"><tspan>*</tspan></text>\n",
       "</svg>"
      ],
      "text/plain": [
       "<IPython.core.display.SVG object>"
      ]
     },
     "execution_count": 70,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "mol = Chem.MolFromSmiles('c1ccccc1CC1CC1')\n",
    "bi = {}\n",
    "fp = AllChem.GetMorganFingerprintAsBitVect(mol, radius=2, bitInfo=bi)\n",
    "Draw.DrawMorganBit(mol, 872, bi, useSVG=True)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "* 展示RDKit的分子指纹<br />\n",
    "黄色表示带有芳香性质的原子"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 73,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAJYAAACWCAIAAACzY+a1AAAMpklEQVR4nO2ce1BTVx7HfzcJIeGR8pRHS0ARK66iCFQ6VC0uPlrRaQe0jsrsDLrU7QOr7RTH6ZZZnVnZGay0Xacy26nFR3cH6UPsA98rKqPWigqVFUGMgWCIQEgg5Hnv/nFtlg0JULj3xtP5fSaTSbjn3POdfPjdV+4JxTAMICQj8nYAZKKgQuJBhcSDCokHFRIPKiQeVEg8qJB4UCHxoELiQYXEgwqJBxUSDyokHlRIPKiQeFAh8aBC4kGFxIMKiQcVEg8qJB5USDyokHhQIfGgQuJBhcSDCokHFRIPKiQeVEg8qJB4UCHxoELiQYXEgwqJBxUSDyokHsnYmtEOx0MAh0gUSlFSfhMhv5JRFFosdSZTlcVylWFMAAAglkp/J5cvl8tXoMvHBMrTTwfRtF6vL7ZYLrpdKpEog4L+6uMznc9syJhwr5Cme7u7/2i33xupJ+UXEvJ3qTSJr2jI2HB7OMPo9X8e2R8AMIypt/cdmtZzHwr5NbhRaDafs1gueepQXz/Q3GxmX9N0T3//P/iKhowNNwpNpiOeWhsMjk2b7i5d2vTpp13sBthkOsYwFv7yIaPiqpBhrFbrNU+tpVJqyZIgi4V+/311fn5rT4+dYUxW6w2eQyIj4arQ4XjAMDZPrWUy0c6dMfv2TVEoxMeP67OybtXVGR0ONc8hkZEYviH16M/JypXBZ87MSE8P0GptAwP0CMoRAXA9tReJwsbSLTpaeuTItNpaY2amQiwO5yEYMlZcq1AkekIiiXX5Y2Vlt07nWmpiMZWZqQCgfHxm8RgQGQ03R6Ry+fKhby9cMG7dei8rq+nsWcPwxlJpqlg8yfn26NGjer2e65DISLhR6O+/WiQKdb59+mnZ888/odPZ1q+/8+67KrOZHto9MPBPzjdXr15dvXr17Nmza2treYyM/D9uFFJUQFDQTuduMjzc5+DBqTt2xEgk1KFDD1944T9NTYPsooCADUMvsIWHh6elpd2/f3/RokXbtm2z2R6LwxyHw3HhwoWOjg6Hw+HtLPzAeGBw8Fxn53yNJsX5qKlJnDJFBgC+vqIdO2L6+vYyDO3Sy2q1bt++XSwWA0BGRoZKpfK0fgEwGo3l5eXTpk2jKGrKlCnp6emtra1ezMMTHhUyDGO3d/T2btdonnFavHs3ecOGR3u+l1566eHDh2471tXVTZ48OTg42FsKOzo6tm3bFhwczEaNjY0NCwsDgODg4CNHjnglEn+MpJDF4egdHDzV319hNH5uMh2z2e5XVVWFhIQAQERExA8//OC2V3d398WLF7lOOzo3b94sKCiQyWSsvIyMjMrKSrvd3tvb+8orr7B/zMvLMxqNwmfjidEVukWlUi1YsAAAKIoqLCy0WCzcxhoH58+fz87OpigKAEQiUXZ2dl1dnUubiooKf39/AJg8efLwpYQyToUMw9jt9pKSEh8fHwBITU1tbm4euT1N04cPH3Y4HOMe0S0Wi6WysjItLY2tsICAgIKCgtu3b3tqf+vWrTlz5gCARCIpLi7mPI/wjF8hy+XLl+Pj4wFALpeXlZWN0LK0tBQAMjMz1Wr1BAdl6evrKysri4mJYeVFREQUFxd3d3eP2tFsNhcVFYlEIgBYtGhRe3s7J3m8xUQVMgzT19e3fv169nPMycnp6elx26ympiYyMhIAQkNDv/rqq4mM2NbWVlRUFBQUxA6alJRUXl4+ODj4q1Zy4sSJqKgoAAgLCzt69OhE8ngXDhSyVFZWsp+pUqmsra1120ar1S5f/ujST0FBnk73aU/P2zpd3sOHG/T6v5hMx2l6lH3qtWvX8vLyJBKJ82ilurqapl3PbcbI0Dx5eXkDAwPjW4934UwhwzBtbW0ZGRkAIBaLi4qKrFbr8DY0Te/b97Fc7pOWFnD//tyh550aTYpW+8Lg4InhvRwOR3V1dVZWFvtxS6XSVatWXblyZeKZaZouKyvz9fUFgMTExPr6+omvU2C4VMgwjM1mKy4uZk/t582b19LS4tLA4ejW6dacPj3j8uVZLv6cD4PhY2d7s9lcUVExY8YMVp5CoSgsLORqb+qksbFx1qxZACCTycrKysZd1l6BY4UsZ8+efeqpp9hP/NChQ0OW2HS6P3gy98sjVaNJ6e//V1dXV0lJSXR0NCsvLi6upKREr9fzEZhhGJPJVFhYyI61ZMkSjUbD00Ccw4tChmH0ev2aNWtcTqX7+w84Jbk8Ghpmf//9dPZ1Xd3MjRuj/PzkbPe5c+dWVFTYbDaeog7l66+/Dg0NBYBJkyZ99913Aow4cfhSyOI8lY6Li7t4sfbBg8VuK6+jI2Xt2rAbN5K++ebpFSuCxWL29JzKzs4+efIkrwmH09nZuXTpUudVC7PZLHCAXwu/ChmGaWpqSk5OBgCJRPz221FqteshjEaTolanlJbGpqT4/3K0QuXmhl64kMl3Nk+wxzhSqRQAZs6c2dDQ4K0kY4F3hQzDWK3W4uJikYgCgIyMwGvXkob6O3w4ISbm0fSMsDDJO+9ENzbOZhfRdL8A8Tzx448/JiQkjOWqhXcRQiFLdfWWiAgfAFAoxPv2TXEqPH16BkVBXJzvjh0xra3JQ+3a7Z2CxXOLwWAoKChg/71efvnlsVz6ER6P02I4Z2Dgn21tf9u69d7Jk30AkJsbWlKi9PMTAcCVK/2pqQGiYV8/R0b+m6IChIk3AlVVVQUFBb29vTExMQcPHly4cKG3E/0fwk0R9fFJCA2VfP751J07Y3x9RVVV3S++2KTX2wHgmWfc+BOLox8HfwCQm5t7/fr1+fPnq9XqzMzMzZs3W61Wb4f6H8IplErniETBFAUbNkw6fjwxMVGekCAPCvI4wVEme16wbKOiVCrPnDlTUlIikUg++uij5557rqWlxduhHiHchhQA+vv3G4172ddmM22zMYGBYvexKGl4eJVYHC1YtjFy6dKldevW3b17NzAwsLS01Lmn9CKCzrX391/nnFUqk4k8+QOAwMDXHkN/AJCenl5fX7927Vqj0fjqq6/m5ubevHnTYvHmxCBBqxAAHA5dT89rdnsbAAXgfmh//7UKxVYhU42DAwcOvP766xRFGY1GAJDJZNHR0VFRUcOflUplYGAgf0mEVggADNNvMOwxmY4B0C6LRKIQhWKzy73Ijy1NTU3PPvtsX1+fRCKx2+0jtAwLC4uMjPTkmP2eZNx4QSGL3a42m49brTdoWg8gkkhifH3TZbLfU5TcK3nGx/bt23ft2rVq1ary8nKNRtPZ2Tn8Wa1Wj3xL7QQr2GsKfxtotdq4uDibzXb79m32BhRPzbRabXt7u9vnkQWHhoZGRkY++eSTkZGR7733HnvBaCiocKLk5+fv37+/sLDwww8/HN8aent7PVVwe3v70HPQhoaGmTNnunRHhROlsbExKSnJz89PpVKxX1Rxi1ar7erqUqvVXV1dOTk5brar3ru299th2bJlALBr1y6vjI5VyAGnTp1avHhxdHR0W1sb+xWVkODP6HFAVlZWcnKyRqP54osvhB8dFXLDli1bAGD37t3Cb9VwQ8oNNpstPj5erVbX1NSw920IBlYhN/j4+LzxxhsAsHv3boGHxirkDIPBoFQq+/r66uvr2Zk3woBVyBkKhSI/Px8A9uzZI+S4WIVcolKppk6dSlFUa2urc8oV32AVcklsbGxubq7NZtu7d69gg2IVcsxPP/2UmpqqUCjUarVCoRBgRKxCjklJSVm4cKHBYPjss8+EGRGrkHuOHTu2cuXK2NjYlpYW51RI/sAq5J7s7OzExESVSvXll18KMBwq5B6Kot566y0AYH9egPfhcEPKBxaLJS4u7sGDB+fOnWN/3YU/sAp5wdfXd9OmTSDI9TasQr7o6elRKpUmk+nnn39OTEzkbyCsQr4ICQnJy8tjGGbc99SMEaxCHrlz58706dOlUum9e/ciIiJ4GgWrkEcSEhJWrFhhNps/+eQT/kbBKuSX8+fPL1iwIDw8XKVSyeW83OWMVcgv8+fPT09P1+l0Bw4c4GkIVMg77G01H3zwAU27ziHhBFTIOzk5OfHx8c3Nzd9++y0f60eFvCMWi998803g7TQfD2eEYGBgIDY2tru7+9KlS/PmzeN25ViFQuDv779x40YAqKys5HzlWIUC0dnZef369WXLlrG/Hs4hqJB4cENKPKiQeFAh8aBC4kGFxIMKiQcVEg8qJB5USDyokHhQIfGgQuJBhcSDCokHFRIPKiQeVEg8qJB4UCHxoELiQYXEgwqJBxUSDyokHlRIPKiQeFAh8aBC4kGFxIMKiQcVEg8qJB5USDyokHhQIfGgQuJBhcSDCokHFRLPfwEbmdyLWCYmCAAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<PIL.PngImagePlugin.PngImageFile image mode=RGB size=150x150 at 0x89B5208>"
      ]
     },
     "execution_count": 73,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "rdkinfo = {}\n",
    "rdkfp = Chem.RDKFingerprint(mol, maxPath=5, bitInfo=rdkinfo)\n",
    "Draw.DrawRDKitBit(mol, list(rdkinfo.keys())[0], rdkinfo)"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.7.4"
  },
  "toc-autonumbering": true
 },
 "nbformat": 4,
 "nbformat_minor": 4
}