Skip to content

Instantly share code, notes, and snippets.

@iwatobipen
Created September 24, 2021 02:44
Show Gist options
  • Save iwatobipen/303f196615e9082863f87c3209d4d290 to your computer and use it in GitHub Desktop.
Save iwatobipen/303f196615e9082863f87c3209d4d290 to your computer and use it in GitHub Desktop.
Display the source blob
Display the rendered blob
Raw
{
"cells": [
{
"cell_type": "code",
"execution_count": 1,
"id": "comparative-cream",
"metadata": {},
"outputs": [],
"source": [
"from rdkit import Chem\n",
"from rdkit.Chem import AllChem\n",
"from rdkit.Chem.Draw import IPythonConsole\n",
"from rdkit.Chem import Draw\n",
"IPythonConsole.drawOptions.addAtomIndices = True\n",
"IPythonConsole.drawOptions.minFontSize = 16"
]
},
{
"cell_type": "code",
"execution_count": 2,
"id": "massive-columbia",
"metadata": {},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAcIAAACWCAIAAADCEh9HAAAABmJLR0QA/wD/AP+gvaeTAAAgAElEQVR4nO3deVxUVf8H8M8wgCCoIIi4C2oiKkjkFi48iQo4CKi4ULhkmVk/LAuw9BfZ83sS0gy3x6Q0KUUUNRgQREBEVEhBEcHcBWIXBJWdGc7vj4sjwrAIM/cyw3m//EPuvTPzceHLvfec+z08QggoiqKojlLhOgBFUZRio2WUoiiqU1S5DkBRyiMzM/Pq1avM73v37j137lxu81DsoGejFCV7mZmZfn5+XKegWMKjQ0wUJXPr16+3s7OztbXlOgjFBlpGKUrGSktL58yZc+XKFR6Px3UWig30op6iZGz//v0ff/wxraHdBz0bpShZqqurmzp16sWLFzU0NLjOQrGEno1SlCwdPXrU0dGR1tBuhZ6NUhRFdQo9G6UoiuoUWkYpiqI6hZZRipKB7Ozs3bt3u7i4nD59uqCggOs4FKvovVGK6riHDx+GhYUFBwdfvnyZ+VZSU1MbOnTo2bNnjY2NuU5HsYSWUYp6PfX19YmJiSEhISEhIffv32c2amtr29nZ2djYHDhw4MqVK4aGhpGRkRMmTOA2KsUOWkYpql1EIlFSUlJwcPCJEyfy8vKYjXp6evb29g4ODvb29lpaWgAqKioWLlwYFRWlo6MjFAqnT5/OaWq2PXnypLa21tDQkOsgrKJllKJaU1lZGRsbGxwcLBQKnz59ymwcPnz4/PnzHRwcrK2tVVWbtkmrra11c3M7fvx4jx49AgMDFyxYwHpqbnh4eFy/fl1DQ2Ps2LG+vr5cx2EPLaMUJUVxcXFERERwcHB0dHRNTQ2z0dTU1MHBQSAQWFlZtf6sp1gs/uSTT/bv38/n8/fv37969WpWUnOpsrJywoQJt2/fVlFRsbS0PH36dPc5J6X9RilFlZeX98MPP2RnZ7u6ui5atEgm7/no0SOhUBgcHJyYmFhfXw+Az+dbWVk5ODgsWLBg1KhR7XwfPp//888/GxkZbdy48cMPP3zy5ImHh4dMEnZZqqqq9fX1FRUVampqPXv2vHfvHi2jFNWl1dXV2dnZ7d6928zMrLq6upPvlpGRERwcHB4enpKSwmzR0NCwsbFxcHCYP39+h8uBl5eXlpbW+vXrPT09S0pKtm7dqsT9StTV1bdv375kyZLevXtXVlZqampynYg99KKeUkhnz54NCgr66quvxGKxiYlJB95BLBYnJiYGBwefOnUqJyeH2airq2tjYyMQCJydnXv16iWTqIGBgStXrqyrq1uxYsWvv/7a/F6qkhGJRObm5klJSbL6C+z6lPxflFJW2dnZiYmJ/v7+2dnZRkZGPj4+7XxhVVVVTExMeHh4SEhIUVERs3HIkCF2dnYCgcDW1lZNTU22UV1dXQ0MDJydnQMCAsrKyoKCgpS1ccmuXbtqamrOnTv38ccfd58aCno2SimooKCg69ev+/r61tfXGxkZZWVltX58SUnJ6dOnw8PDIyMjy8vLmY3GxsYCgcDFxaXNIaPOu3Llir29fUlJibW1dWhoaO/eveX6cZzIycm5deuWiYnJ0KFDuc7CKlpGKYWUm5trb29//vz5lJSUzZs3JyUlST0sKysrJCQkPDw8Pj6+rq4OgIqKioWFhUAgWLp0acfuBrTi2rVr5ubmfD5f6t6MjAxbW9ucnBxLS8uIiAgDAwPZfrrc1dXh6lVkZYHHw8iRsLBA4z/prVsoLsa0aVB59RHznBw8fAgzM+josJyXPYSiFJNQKHRwcFi1alVOTk6TXenp6T4+Po3PMVVVVa2srPz8/HJzc+WU5/z58xoaGo6OjlVVVS0d8+jRozfeeAOAsbHx/fv35ZRELvbtI/r6BHj5a9AgcvToywOcnAhAKiqavtDXlwAkJobNsCyjZZRSKh4eHsOHD5ecJfTp02fZsmXHjh179uyZvD86KSlJT08PgLW19dOnT1s6rKCgwMLCAsCAAQNu3Lgh71Sy8f33BCATJ5KICFJQQPLzSXAwGT2a8Hjkl18ajqFllKKUwLZt28aMGQNAX1/fzc1NKBRWV1ezGSA9PX3QoEEALC0tCwsLWzrs+fPnNjY2AHR1dS9dusRmwo64d4+oqhJz86Yl8vFjMngw0dIijx8TQssoRSk+sVjcp08fAGFhYWKxmKsYjx49Ymbpt37ZXl1dzTwyoKWlFRERwWbC17ZpEwHIqVNSdu3ZQwCyaxch3bqM0n6jlJLIyMh4+vTp8OHDBQKBigpn/7GHDx+ekJBgYWHx8OHD6dOnp6WlST2sR48eQUFBH3zwQUVFhZOTU1BQEMs5X0NiIgBMnixll43NywO6MTpvlFISiYmJAKZOncp1EPTv3//ChQvOzs4xMTHW1tbh4eFvv/1288P4fL6/v7+enp6vr6+rq2t+fv7nn3/Oftq25eVBXR1SH+UaMqThAAlLSzSZOlZSIs9wXQI9G6WURNcpowC0tbXDw8MXLlxYWlo6Z86cM2fOSD2Mx+P5+Pj4+fkB2LBhw8aNG9mN2T6EgM9vWhwZzIQnsfjllrffxvTpr/wyMmIpJ4e4vqtAUbIxevRoAFevXuU6yEsikeiDDz4AoK6ufrTx3KBmAgICmIdE161bx+GN3ZdKS8mRI+TrrwkhxMqKAKSoSMphmZkEIAsWEELvjVKUgispKbl7966mpqaZmRnXWV5iLts9PT1ra2tdXV1/+umnlo5cvnz5qVOnNDU1//vf/y5atEjSmq8lRE5PzTx+jN9/h4MD+vfHu+/Cxwd5eXjzTQBIT5dy/OXLAPDWW3IJo0C4ruMUJQPh4eEAZsyYwXUQ6fz8/JgHAby8vFo57Pz588xkg3feeaelia51dXXLli1zcnKaOXNmeHi4bPLdvk22biWTJxMer2FqvaoqeecdsmsXKS0lly8TgDg6Nn2VWEwmTyZ8Pnn4kJBufTZKyyilDDZt2tRSkaqoqKivr2c/UhPtvGy/efPmwIEDAUycOPExMx/zVampqbNnzyaEpKWl2dradipTejrx9iaWli8fTNLUJAIB2b+fFBS8cuSSJQQgn3xCJMX98WOydCkByIYNDVtoGe0mWnm2hFJo77zzDoCQkJDmuzw8PHR1dQMCAthP1URoaCjThdPZ2bmV5wIePnw4cuRIAGPGjMnOzm6yt7Kycvz48T4+Pp999tmRI0deO4RIRBISiLs7GTz4ZfXs25e4uJCAANLSs16VlQ1FU0ODmJuT8eOJmhpRUSGffELq6hqOoWW0O0hMTOzbt++JEye4DkLJmEgkYhom5efnN99rZWUFIDIykv1gzUku22fNmtXK86n5+fnm5uYAhg0bdvv27ca7Hj9+PG/evB07dkydOnXnzp3t/Ny6589JcDB5912io/Oyeg4bRtzdyblzRCRq17tcvUq8vcmKFWTVKvL99yQ9/ZW9J0+SrVtfVlWJpCSydSvJzGxnVEXUjcqop6cnAFVV1YMHD3KdhZKl1NRUAMbGxs131dbWampq8ni8kpIS9oNJ1eZlO6O0tHTatGkA+vbtm5iYKNnu7++/ZcsWQkhtba2BgUHrn1VcXBwQEODi4jLa0JCoqzdUT1NT4uVFEhJIF7jXoRyUtoyKxeJ79+41uYpnmvsyk/W4CkbJ3L59+wC8++67zXcxDfRMTU3ZT9WKBw8etHLZLlFRUWFvbw9AW1s7KiqK2fjXX39NnDgxLi7uv//97/z581t6/x9//HH69OmSln0qKiq5rq5k+3aiWG2lFIRyltGqqqqZM2euWbNm4sSJQqGw8a49e/YwTwq6u7t3hZEHqvOWL18OYM+ePc13MXOMVq9ezX6q1rVy2d5YXV3dqlWrAKirqx87dozZmJaW5ufnd/To0SYd+Zq3B+zRo4eNjY1c2wNSRFnL6OHDhz/77DNCSG5urrm5eZO9gYGBzEIRy5cvr2t+K4dSNEwrkJSUlOa7Fi9eDODXX39lP1WbSktLmfu2TS7bm6ivr//iiy8AMGs1N9krEokSEhK8vLyYNqYMLS0tgUAQEBBAx1TZoZxl1MvLSzKIaWRk1LxWxsTEMGvFODg4VFZWsh6Qkpni4mIej6elpSX1J+KQIUMAZGRksB+sPaRetkslWWyKmdRVVVUlFArXrFnTeNXSfv36Me0Ba2pq2PoTUIQoaxn19PQMDAxkfj9ixIja2trmx1y5ckVfXx/AzJkzy8rK2A0oRyKRqLy8nOsU7BEKhQCsra2b78rNzQWgo6PTJR6vbEHjy/bjx4+3cuS+ffuYe51jxozp2bOnpHqamJhs3LgxKSmJ3qTiinI+DDpmzBimQVlJSYmmpqbUtR4nTpx44cKFIUOGxMfHT5s2La9xlxqFFRwcPGXKlCVLljg4OIhEIq7jsKGVjiSXL18GMHnyZA775rVJVVX1wIEDX3zxRW1t7bJly/z9/Vs6cu3atSdOnNDU1CwsLKysrDQ1NfX29k5OTv7777+3bt06efJkea/KR7WI6zouF+Xl5ZMmTdq8ebO1tbXktFSqzMxMpqWFkZHRvXv3WEsoJ6dPn37+/DkhZNGiRXFxcVzHYYO1tTWAJgOJjA0bNgD49ttv2U/VAd999x0AHo/XymzQnJwcANra2v/88w+b2ajWKWcZJYRUVVVdunQpKyurzSNLSkqmTJkCwNDQMDU1lYVsLBAIBNevX+c6hdyJRCJtbW0ABU0eXiSEEMKcorZ+z7FL2bdvn46OjtSxMsaxY8cA2NnZsZmKalPXvdjpJA0Njbfffrs962X37ds3Ojp6zpw5BQUF1tbWFy9eZCGeXIWGhhJCJkyYwHUQuUtLSysvLx85cmT//v2b7Kqpqbl+/TqPx5s4cSIn2Tpg7dq19+/ff5PpqCRNl+qpSkkobRl9Ldra2mFhYS4uLmVlZXPnzo2IiOA6UcedOnVq9+7dR44c4ToIG1opK9euXauurh47dqyuri7ruTqOWVu0JbSMdk20jDZgGut+8MEHlZWV336boKBV6Jtvvvn0008dHR2PHTuWkZHBdRy5a6WsKF/FqampSU1NVVFRUaDz626CrsX0EtNkd9iwt319Vy5fjpISuLt35H2uXbv2448/VldXnzx5UtYZ2zB16tTx48czv+/RowfLn86+blVGU1JSampqxo8fzzQ3oboOWkZfwePxNm9e1bs3Pv8c69cjLw8vZj2/hpKSkg0bNqxfv14OAdtgZ2fH/odypaio6MGDB1paWuPGjWu+l3maXpnKqPL9YFAa9KJeCnd3BARATQ2+vli1Cq87/3L27NnDhw+XSzKqEaasTJ48mWmH3Ng///yTk5Ojo6PT+BFJRUfLaJdFy6h0772HU6egqYlDh+DigupqrgO1RSQSnTt37sqVK/n5+VxnYUl7rui78sT716V859dKQ3n+k8mcQIBz56Cnh5AQ2Nvj2TOuA0lTVVUVFhb20UcfDR48eNasWWvXrn3jjTeioqK4zsWGbnVjNCsrKzc3V09PT5nOr5UGLaOtmTIF8fEYNAhxcZg1C48ft+tVmZmZ8fHxZWVlMTExRD4rOJaUlBw6dMjJyUlPT2/+/Pn+/v6FhYWmpqYAysvLHR0dT5w4IY/P7TpEIlFKSgqPx5s8eXLzvcpXRpk/0ZQpU+gTn10QLaNtGDsWCQkYNQrJyZgxA9nZbb+ksLDw3r17bm5uzOMoMgyTnZ3t7+/v4OAwYMCAVatWhYaGVlVVMc9WZ2RkZGRkpKSkfPnllzU1NUuXLt2/f78MP7qruXHjRkVFxahRo5j+Mo0p5cQg5fvBoFQ4fopKQeTnkwkTCEAGDCBpaWx/evN2vHw+38rKys/PLycnp/nxPj4+7VnOV6Ht3r0bwIoVK5rvYp5DMzMzYz2UHDE/EmJjY7kOQklBy2h7lZaS6dMJQHR1yeXLcv84sVicnJzs7e3NdE5h9OzZk2nH22Znv0OHDjHj159++mlXbhPXYa6urgB+/vnn5ru2bdsG4KOPPmI/lZxUVlaqq6vz+XzahrlromX0NVRXkwULCEC0tIicFpqsqiJhYeTDD0VDhw6TVE9DQ8M1a9ZERES0sipvcyEhIRoaGgBcXV2ldlxVaEZGRgBu3LjRfNeCBQsAHDp0iP1UcpKQkACg+ToOVBdBy+jrEYnI6tUEIOrqJChIZm9bXk6EQuLmRnr3bli9cfJkRyMjI3d39+jo6A6vdBIXF8esPGxjY8M00FMOhYWFAHr16iWStjLwoEGDANy5c4f9YHLyww8/AFi7di3XQSjp6FNMr4fPxy+/oE8f7NgBV1doaeH2bYjFWLkSjXsMicXYvh0mJnB0bO3dcnMhFOLPP3H+POrqAIDHg6UlnJ3h7HzY1FS7k2mtra3PnTtnZ2cXExMza9asiIiI1jtfKAqmH/OkSZMkK19KSCYGMQs0KQc6vtTVcV3HFZWfH5kyhZSXN6z+7er6yt7qagKQJUukv/bBA+LnR6ysiIpKw7knn0+srIiPD5FH5+j79++PGDECgKmpqXK0+/X09ATwv//7v813HT16FMC8efPYTyU/zIJLStBWXFnRs9EOWr8e69aBWZ1EXR2BgVi5ErNnt/aSjAwEByM4GLduNWzR1MSsWXBwgKMjmjXMlJkRI0YkJCTMnTv35s2b06dPj4qKUpQp3Fu3br148aKWltaECRO+/vpryXbmbLSbTLx/9OhRQUGBvr4+87OQ6oJoGe04yQpPc+fi1i2sW4ebN6Gh8coxYjESExEcjJMnkZvbsLFvX8ybBwcH2NlBu7MX7u0yYMCA+Ph4gUBw+fLlGTNmREZGWlhYsPHBnZOenr59+/YxY8Y03lhXV3ft2jUejzdp0qTmL1G+Mir5E9GJ910WnX4vA+rq+P573L8vpR3UkSOYPh27diE3F8OHY/16xMWhqAi//w4XF5ZqKENXV/fs2bO2traFhYUzZsyIiYlh77M7qqCgIDo6eu/evUVFRZKNqamplZWVo0ePbn6ft6qq6saNG3w+n068p9hEy6hsLF6MmTPh44M7d17ZLhDAwgLffINr1/DoEfz8YG2NZuMiLNHS0goNDV26dGl5eblAIOj6D4xu2rTJ3Ny8f//+c+fOJS+eB2ulrJSVlc2bN2/WrFm9evViNag80TKqALi+Oavw1NXJwoWEEJKeTtTUyMyZpL6+jSEmbtXX13/++ecA+Hz+L7/8wnWcdpk0aZJkcGzp0qUA/P39uY3EjsrKSjU1NT6fr0zz1ZQPPRuVmbFj4e6O+HgEBXEdpVU8Hm/Hjh0+Pj5isXjNmjW+vr5cJ5Kurq5u9erVUVFRv//+u1gsHjBgALO9W52dXb16ta6uzszMTJvNG0DUa6JlVJa8vTFoEDw9UVHBdZS2eHl57d27l8fjbdy4cf369UQ+nag6Q01NzdPTk7kTGh0dzUwRzc/Pz8rK6t27d5NxJ2XFzEl4++23uQ5CtYaO1MtSr1746ScsXowdO7iO0g7r1q3T1dVdsWLFrl27SktLDx482LyNPLdGjx7t5eUl+bK4uPj7778HYGZm1nzivVLqVqfeiouejcqYiwvmzMFPP3Gdo32WLVsWGRnZq1evP/74Y8GCBVVVVVwnkkLSHnDgwIF79uwZPHhwamrqpUuXuM7Fhr/++gu0jHZ5tIzK3t69qK/nOkS7zZo1KzY2Vl9fPywszM7O7unTp1wnapCWlvbdd9+9+eabw4YN++ijj8LDw3k83pw5cwYOHFheXj5nzpzIyEiuM8rXgwcPCgsLDQwMjI2Nuc5CtYrrMS6FJxmpb2zz5q47Ui/VrVu3hgwZAmDcuHG5ublcxWhPe0CRSPThhx8CUFVVPXjwIFdRWfDHH38AcHR05DoI1QZaRjsrJYU8eNB0Y1UVSU6Wsr0ry8zMZIqXkZERy49vV1VVRUdHu7u7S4bjAejr67u5uQmFwubtAevr6729vQHweLzt27ezGZVN69atA+Dj48N1EKoNtIx2Smkp2byZCIVSdvn4kL17WQ/UOSUlJVOmTAFgaGiYmpoq748rLS09fvy4m5sb082P0f72gDt37mQW/lTWJv/MA7vx8fFcB6HaQMtopzx8SAAitc+6gQFRxDa7z58/nzNnDgAdHZ2EhAR5fERRUVFAQIBAIFBXV5dUT2ZFqeTk5Nd6qz/++ENNTQ3Axx9/rGRN/svLy1VVVdXU1CoqKrjOQrWBltFOUb4ySgipqalxcXFhbkqePn1aVm/74MEDPz+/5itK+fj4dOYeQlhYmKamJgAnJ6eqqipZpeXcuXPnALz11ltcB6HaRstopyhlGSWEiESiNWvWyGQYJz093dvbm1n8maGpqSkQCPbv319YWCiTtPHx8X369AHwr3/969mzZzJ5T8795z//AfA///M/XAeh2kbLaKcoaxklnRvGEYlECQkJ7u7uzHoejL59+7q5uR0/flwej4ffvHlz4MCBzOlbUVGRzN+ffQKBAEBgYCDXQai20TLaKUwZfe898s8/TX/p6yt2GWXs2rWr/cM4FRUVQqHQzc1NR0dHUj2HDRu2Zs0aoVAo72X1Hj58OHLkSAAmJiZZWVly/azOa/1MvL6+vl+/fgAePXrEViKq42gZ7RSmjLb0SwnKKGk0jLNy5Uqpo+ePHz9mhox69OjReMjIy8srISGhvr6etaj5+fkTJkwAMHDgwLS0NNY+t/0ePXrk5+dnY2Ojqqp6/fr1lg67e/cugP79+7OZjeqwrvUMtYKaMwdr1jTduHo1F1Hk4L333tPR0Vm8ePGhQ4fKysqOHj3KrNucmZkZGhoaHh5+/vx5kUgEQEVFxdLSUiAQLFu2rPH8edYYGhrGxcU5ODhcvHjR2to6PDy8izxGmZKSEhISEhISkp6ezmzR0NBIT09nin5zzKP0tCOJwuC6jis2Jb43WlVV9ffff5eXlzNfJiYmMt3mJ0+e7OHhMW7cOMl/IQ0NDQcHhwMHDnSRm5LV1dXMUvVaWlpnzpzhKgZzg9jLy6vxGqU6OjouLi4BAQFPnz5t5bVr164F8MMPP7CWluoMWkY7RVnLaFpa2rRp0z799FMLC4vExERmY2pqqr6+vqSxvK6ubnsqAidEItH7778PQF1d/dixY2x+dGVlpVAoXLNmTf9GixT269ePeSKrpqamPW9ibm4OQE7zdimZo2W0U5S1jObl5eXl5RFCDhw48PXXX0u2M/NJR48eHR0dLe8ho06qr6//8ssvAfD5/H379sn740pKSpgnshqvX2JsbOzu7v66N4ifPXvG5/PV1NQqKyvlF5iSIXpvlJJiwIABpaWlu3btioqK8vPzYzbW1NRER0cDcHZ2trGx4TRg23g83rZt2/T19b/66quPP/44MzPTp/mKg52WnZ195syZsLCwqKiouro6ZqOpqamLi8vixYsbz5ZtvytXrojFYktLS+axAqrro2WUkk5dXX3kyJEXLlxITk5m7u6dPn26rKwMQOMGIl2cl5eXoaHhBx984OvrW1FRIXkMv5MyMjLCw8PDwsIuX75MCMGLJ7JcXFwWLVrUeLbs60pNTWUm3nesBFPc4Pp0WLHV1JCUFJKdLWVXWhr5+2/WA8lIbW0t84h6fHy8q6srs3HhwoXM/5njx49zmu61hYSEMLMLXF1dO3wvQtLEz8TERPLtI2niV1pa2uF4YrGYGYx64403mLc1MDDQ09O7du1ah9+TYhMto5QUUVFRAoFg586dM2fOFAqFhJCnT59KrjEvXrzIdcDXFhcXx/SRsrGxea3HqKqrq5s38dPT02OeyJLMZOgAZjDq/fff19fXl7zzgAEDVq9ezfTZ0tbWjo6O7vD7U6yhZZSS7p9//omMjMzMzGS+PHjwoORb/YFiNVJ9ITk5mXk0aNKkScXFxe18leQMEcCIESO+/PLLixcvdqab1JMnT1pqD5iQkMC8c01NDbOOdI8ePU6cONHhz6LYQcso1S6Nx5QUdwT5/v37I0aMAGBqaipZ+L51a9eu7VgTvyZetz1gfX39Z599BoDP5//yyy+d+WhK3mgZpdqWl5cnWYlTR0eH6zidkpeXZ2ZmBmD48OF37txp83iRSNSZj2upPaCfnx9Tx2NiYrZs2cIcnJOT4+Xl5eHhIbkIYGYX8Hg82gO/K6NllGrbjkYLRpuamnIdp7OePHnCPGfZv39/OQ3jMO0Bx4wZI/l7k9oe8LffftuyZcuIESMIIfX19RYWFmfOnImNjTUzM5OU7z179jCzC9zd3dlsUEC1Hy2jVNveeustSTmwsbHhOo4MlJeX29raynYYp8PtAU1MTAgh169ft7e3Z7Y4OzsnJSVJDggMDGS6wyxfvrzNtVUo9tF5o1Qb7t+/n5ycLPlSgSaNtkJLSys0NHTFihVBQUECgeDw4cOLFi3q2FtVVlbGxsYGBwcLhULJ8tTDhg2bO3euQCCwtbVlKmCbCgoKJH+3gwYNys/Pl+xatmyZgYGBs7Pz77//XlpaeuzYMTozv0uhZZRqw+HDhxt/yXRHVgLq6uqBgYEDBgz46aefli5dum/fPmbd5nYqLi6OiIgIDg6Ojo6uqalhNpqamjo4OAgEgsY3Q9tJV1eXeboBQGlpqa6ubuO9s2bNio2Ntbe3DwsLs7OzCw0NZRr+U10C16fDVFfXeMYPAD8/P64TydhrDeM8fPhQ0jCU+QtRUVFhVpRqz4CVVMxFfXV19ejRo/Py8oqKikxMTKSuZHfr1q0hQ4YAGDduXG5ubsc+jpI5Wkap1ly5cqXJz12WGyaxY+/eva0P4zBDRpaWlpK/Bw0NDRsbGz8/v/z8/A5/7vnz521sbCZNmmRjY5OTk5OQkGBnZ2draxsXF9fSSzIzM5lerkZGRp1ZCpCSIVpGqdYwUxcbU9bubZJhHDc3N2YYR9IwlFmbhNEV2gMWFhYyBX2qiYn45k2uYlASPEKIbO4OUEqnvr5+6NChubm5jTdKZrArn9jYWGdn5+fPn0+dOtXY2DgyMvLJkyfMrsGDBzs6Ojo5Oc2cObOdQ0Zy9ezZs5VLlhzKz++dlYWwMEybxnWibmPlJ5QAAAgVSURBVI2WUapFMTEx/rNn9wSSgDsvNpaXl2tpaXEZS56uXr1qb2/fq1evR48eATA2NhYIBC4uLh0YMpK72lq89x6Cg9GzJ4KDYW/PdaDui5ZRqkVPpkzp+9dfzO9LgETg3716/fXsGbep5O3OnTvl5eWxsbFOTk5Nhte6HLEY69bB3x+qqvD3x6pVXAfqpmgZpVpQXY2+fVFV1XjbzX79xhcVcZWIkoIQbNmCLVvA42HbNnzxBdeBuiMZtLCllFN4eJMaCmDc2LGcZKFaxOPh22+xcyd4PHz5JTZu5DpQd0TLKNWCwMDm23iDB7MfRN5KSkok8+dzc3Pv3r2reJdo7u4ICICaGnx9sWoVRCKuA3UvtIxS0pSVISJCynaleBJUori4eOXKlWZmZhEREQA8PDw++ugjHx8fgUCgeJX0vfdw6hQ0NXHoEFxcUF3NdaBuhJZRSpoTJ/DiBO0VylVGe/fu/eOPP65YsYL5cuHCheHh4QcPHiwtLWVG6hWMQIBz56Cnh5AQ2NtD2QcDuw5aRilpjh6Vvl1ZHqhnqKur6+npSb5klu4ghJSWljJ98hXPlCmIj8egQYiLw6xZePyY60DdAi2jVDN5eYiPl75Luc5Gpfr6668XL17ceLl5BTN2LBISMGoUkpMxYways7kOpPxoGaWaOXoUYrH0Xcp1NtqESCTy9PSsqKj49ttvuc7SOUZGuHABEybg9m1MnYqbN7kOpORoGaWakTZG38DQkMUcbNi9e3diYuKff/4ZGxs7d+7cs2fPDho06Icffnjw4AHX0TrH0BBxcZg2DXl5cHJCXR3XgZQZnX5Pver2bTRa+uIVffrgRUNMpXHjxg2RSATAwMDg8ePHkm+HkSNHKkNDz6oqLF+Odevwr381bLl9GwEBSE3F06fQ08PUqXj//Zc/He/dg7c3Fi3CggVN38rVFebm8PJiL7zioGWUetU33+Df/5a+a/Ro3L7NbhpKprZvx1dfQVUVU6dCXx9ZWUhORs+eOHwYjo4AcPkyrKzw7bfw9n7lhYRARQVz5+LMGU6Cd3G0+z31qqCgFnc1WmKIUjwnT8LDA1On4sSJl/e4r1/H/PlYsgRXr2L8eE7zKTB6b5Rq5K+/cO9ei3u7wTC90iIEmzZBWxunTr0yTmhhgSNHUFODLVu4C6fwaBmlGmllcAlKPkyv5P7+G3fuYN48KYOEM2Zg7FhERKC2lotkyoBe1FMviMU4fry1A+jZqOJKTweAUaOk733zTWRk4P79hi+zslqcOExJQ8so9UJsLAoKWjuAllHFxUyxaOl6om9fACgtBdOa+rff8NtvbCVTBrSMUi/8+WcbB9CLesXFrGvf0nw1po9Jz54NrRHXrsWaNU2PefNN+aVTdLSMUi84OeHnn1s7gJ6NKq5hwwC0+GDo/fvg8TBsWMOENkNDWFi8cgCdFtkqWkapF+bOxd27SEpCYiIuX0Z6etNHQmkZVVxvvQVtbcTFgRA0WVTqyRNcugQzs4ZLe+r10ZF66oXqasTE4ORJJCZCTQ1OTti0Cd99B4EA/ftj9278+CMcHFBf3/SFQUFwcMDff3MRmmqfnj3x/vu4cwd+fq9sJwQbNqC6GuvXc5RMGdCzUQoA8PAhbG1x7x7MzGBhAZEIly7h5ElYWUEobDhPsbdHZKSU67u7dxEeTpev6Or+7/9w/jy++AI3bmDxYvTvj4cP8euvOHsWCxZg5Uqu8ykwWkYpoLYWzs7IzERgIJYta9goFuM//4G3N1auhFDIaT5KFnr1Qnw8Nm9GQAACAho2Ghjg++/h4dH0Sp96HbSMUsDx40hLw8aNL2soAD4f33yDq1cRFoarVzFxInf5KBnR0cGePdixA3fuoKwM/frhjTeg0nBnr7i4+PSdO25Pnqj07FlXVxcWFmZpaTmMGZvi8VBWBjU1LsN3YfTeKAWEhQHA8uVSdjETX5gDKOWgro7x4zF9OkxMJDX09OnTn376qefGjSItraKnT5cuXbpt27akpKSXr+rTBz17chO4y6NnoxSQkQEeD8bGUnYxE1+YZ2AYp09Lvvca3L0rz3AUG+bNmzdv3jxzc3MABgYGJ0+eVPje1SyiZZRCQ+vJHj2k7NLRAV6dts10VKMo6gVaRimgZ08UFkIsBp/fdFdlJQBoab3ccudO07PRPXuwc6ecI1JU10XLKAUYGeHuXeTnY/DgpruYC/bG1/sjRjSttnTaNtW90SEmCg0rTERFSdkVEgIA77zDah6KdXfu3Jk9e3ZxcbG9vX1KSsqyZcuCg4O3bdu2a9curqMpALqICAUUFWHUKPTujaSkV1rcX78OKysMGYKMDKiqNky/F4mano1+9x28vXHxIqysWA5OUV0BvainAAMDHDiAZctgaYnPPoOlJUQiXLiAPXugro7Dh6FK/59QVIvotwcFAFi0CIaG2LQJmzc3dCTp0QO2tvD1xejRXIejqC6NXtRTryorQ24uVFQwdOgrA/QACGlYIbK5+nrp2ymqG6BllKIoqlPoGQRFUVSn0DJKURTVKbSMUhRFdQoto1SD0NDQPn36HDlyhOsgFKVgaBmlAODZs2e+vr5r167lOghFKR5aRikA8PT09PLy0tbW5joIRSkeWkYpXLhwITc315F2wKOoDqHzRiksX768oKBAR0fn1q1b2tra4eHh+vr6XIeiKIVByyj10r///W9jY+N3332X6yAUpUhoGaVeYs5Ghw4dynUQilIktIxSFEV1Ch1ioiiK6hRaRimKojrl/wG9kEftmNMigwAAARt6VFh0cmRraXRQS0wgcmRraXQgMjAyMS4wMy41AAB4nHu/b+09BiDgAWJGBgjgheIGRjaGBJA4MzuYZoLymZjYFEzAfBYOiAJGNgcNIM3MwuaQAaKZGZEYEBkOBrAAIzYBqFJuoAsYmRiYmIGyLIzMrAysbBlMbOwJ7BwZTBycChxcCVzcClw8GUw8rAkijGysPFwc7GziMFcz8GrNXXZg1q9qexBnqUvzAcPKX3tA7FsftA78871mB2IvqD+xP3W5xH4QO9/v036VO8wHQOxvpVz7ijeJ7QOxaxhc7Dz6VoDN+bjhvP0a7kdgdlOnnsMXDwcw+5dAq8O9hTPBbOlFkg4WDWZgM1WDExz81O+C2SvVKuwb9k4Bs8UAg/5BwiTV7i8AAAEaelRYdE1PTCByZGtpdCAyMDIxLjAzLjUAAHicnVPBTgUhDLzvV/QHJJ0WCpxdT8Zn4sF/8O7/x0J9ZA+avEAaMrPdTqYFDhrr43z9+qa15DwOgo4g/jN67/QpzHz4708lCVsdlZx6yxrIs0zP9J/ENaZKTmDBQF6rbBeV2+MqkqyCozaj2Z4XpGKWw0vpttkRUuvK4YqNryrvj6twghb5RbnKlhevQMfoiJAU1vZUkEzEAlUz3lORVKUiamsB9lT8vpSm4cWv4t503UtmLjHdCul7XjQVLogz997anhd/PdY1VGCyNV03IXMfxJEu4igv4qhMgiC2Mo7q/DSIo7aIo75Ip/HKsAjuGUdYDgbRldGpPskb0cvtPH4AMPOxZ3e0dgEAAAEMelRYdFNNSUxFUyByZGtpdCAyMDIxLjAzLjUAAHicTY6xTkMxDEV/hbFIqeVrx06ctyB1YSrsiAE9MaIi1LEfjwNDOsXn6jo+p/Pp7fT0/H54edyx7/t8/+aH2+FoJOwehSl6ldCyHSuBBVqOTKxszWqGQt44LHsV3XxGIHPXmD0LV/7PeihLyT47KzJigprAZw+1wa1sOQUq9wJSeM9NkItYZNDcopdNqElDzYPNgNTZ0su6SlaEWWxWKrPJ/LhBeuQ3SsbGraRIrjefp8JD60zgwloey8f18vX6c/keTHM8X66fxAMLMGSBDF2goy6owxbY8AU+2oI2+oI+YkEM3CkgHe4lMHCnAbn9AiePdCMSK+CYAAAAAElFTkSuQmCC\n",
"text/plain": [
"<rdkit.Chem.rdchem.Mol at 0x7f664716f7b0>"
]
},
"execution_count": 2,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"epinephrine = Chem.MolFromSmiles('CNC[C@H](O)c1ccc(O)c(O)c1')\n",
"epinephrine"
]
},
{
"cell_type": "code",
"execution_count": 3,
"id": "official-jumping",
"metadata": {},
"outputs": [],
"source": [
"bi = {}\n",
"ecfp6 = AllChem.GetMorganFingerprintAsBitVect(epinephrine, radius=3, nBits=1024, bitInfo=bi)"
]
},
{
"cell_type": "code",
"execution_count": 4,
"id": "wooden-romania",
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"{1: ((3, 0),), 18: ((7, 3),), 33: ((0, 0),), 80: ((2, 0),), 128: ((1, 0),), 202: ((3, 3),), 227: ((4, 1),), 244: ((5, 3),), 266: ((8, 3),), 315: ((8, 2),), 338: ((3, 1),), 341: ((5, 1),), 343: ((1, 2),), 356: ((5, 0), (8, 0), (10, 0)), 406: ((0, 1),), 410: ((5, 2),), 426: ((6, 3),), 451: ((7, 2),), 578: ((9, 1), (11, 1)), 583: ((8, 1), (10, 1)), 589: ((10, 2),), 606: ((2, 1),), 632: ((3, 2),), 726: ((6, 1), (7, 1)), 760: ((12, 2),), 807: ((4, 0), (9, 0), (11, 0)), 813: ((12, 3),), 829: ((6, 2),), 844: ((2, 2),), 849: ((6, 0), (7, 0), (12, 0)), 875: ((12, 1),), 899: ((2, 3),), 971: ((1, 1),)}\n"
]
}
],
"source": [
"#bit info has index of bit as key and (atom index, radius) as value\n",
"print(bi)"
]
},
{
"cell_type": "code",
"execution_count": 5,
"id": "limited-roads",
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"0 1 ((3, 0),)\n",
"1 18 ((7, 3),)\n",
"2 33 ((0, 0),)\n",
"3 80 ((2, 0),)\n",
"4 128 ((1, 0),)\n",
"5 202 ((3, 3),)\n",
"6 227 ((4, 1),)\n",
"7 244 ((5, 3),)\n",
"8 266 ((8, 3),)\n",
"9 315 ((8, 2),)\n",
"10 338 ((3, 1),)\n",
"11 341 ((5, 1),)\n",
"12 343 ((1, 2),)\n",
"13 356 ((5, 0), (8, 0), (10, 0))\n",
"14 406 ((0, 1),)\n",
"15 410 ((5, 2),)\n",
"16 426 ((6, 3),)\n"
]
}
],
"source": [
"for idx, onbits in enumerate(ecfp6.GetOnBits()):\n",
" print(idx, onbits, bi[onbits])\n",
" if idx > 15:\n",
" break"
]
},
{
"cell_type": "code",
"execution_count": 6,
"id": "working-accreditation",
"metadata": {},
"outputs": [],
"source": [
"def getSubmolRadN(mol, radius, rootedAtAtom):\n",
" submols=[]\n",
" env=Chem.FindAtomEnvironmentOfRadiusN(mol, radius, rootedAtAtom)\n",
" amap={}\n",
" submol=Chem.PathToSubmol(mol, env, atomMap=amap)\n",
" subsmi=Chem.MolToSmiles(submol, rootedAtAtom=amap[rootedAtAtom], canonical=False)\n",
" submols.append(Chem.MolFromSmiles(subsmi, sanitize=False))\n",
" return submols"
]
},
{
"cell_type": "code",
"execution_count": 7,
"id": "limiting-highway",
"metadata": {},
"outputs": [],
"source": [
"# 266 ((8, 3),)\n",
"res = getSubmolRadN(epinephrine, 3, 8)"
]
},
{
"cell_type": "code",
"execution_count": 8,
"id": "young-mixture",
"metadata": {},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAcIAAACWCAIAAADCEh9HAAAABmJLR0QA/wD/AP+gvaeTAAAYPklEQVR4nO3de1BTVx4H8F/CQ14C8tL6QCsIKKICCq6KAtLiI9jHlnF2urHtrNLHzIKP7mKn2rhj14Euaqx1W1brNtp1KtXZWRQfBRRfVXlYEXmqqCiKvBXBhpic/ePSiBKR5CbcAN/PHx1Izj33lw7z9Z5zzz0RMcYIAAAMJRa6AACAvs1S6AIABqibN2+q1Wrtr0OHDrW3txewHjCYCIN6AEFER0c3NDRwPxcXF+/fv3/BggXClgSGQYwCCOzGjRszZsyorKy0sbERuhYwBOZGAQS2efPmDz/8EBnad+FqFEBITU1Nvr6+JSUlbm5uQtcCBsLVKICQUlNTf//73yND+zRcjQIIRqVSeXt7HzlyZPz48ULXAobDgifgJTc39+bNm9zPwcHBY8eOFbaevmXPnj1TpkxBhvZ1iFHgZf369fb29mPGjCEiHx8focvpYwoLCxMTE4WuAvjCoB54CQoK+vbbbwMDA4UuBEAwuMUEvNTU1Ny5c+f7778vKSkRuhYAYSBGgZfo6Ojs7Ozy8vK5c+fu2rVL6HIABIBBPRjHvn37ZDJZcXGx0IWYH6WSWltp8GCyshK6FDAJXI0CL48fP+Z+sLW11Wg0whZjXq5dow8+oFGjyMaGXF3J2prGjaNVq+ju3aeaff01ubrSkiW6O8nJIVdXmjatF+oFg+FOPRiutLRUIpGsXLnS0dFx/fr1K1asELois7F7Ny1dSu3tZGNDYWE0ZAjV1lJ+Pm3aRNu3048/UnR0R8tHj6ixkR480N1Pezs1NtKQIb1WOBgAV6NguPHjx//nP/+pr68vKCj45z//+eGHHwpdkXnIzqZ336X2dlq+nO7do5Mn6X//o7Nn6fZtio2llhZ67TUqKhK6SjAaXI0CL9OnT58+fbrQVZgTjYY++IA0GoqPp82bn3pr6FD64Qdqa6OMDPrznyknR5gKwdhwNQpgVJmZdPUq2dnR+vU63hWLaeNGEonoxAkqLe314sAkEKMARsVdY0ZEkKOj7ga+vjRx4pOW0PdhUA9gVGVlRNQRlM8TEEBFRU9djTY00KlTOlpevmzU4sAkEKPAS1ZW1po1a4KCgh4+fJiSkuLh4SF0RUJrbiYicnburo2LCxFRU9OTV06fptmzTVkWmBBiFAxUVla2cuXKw4cPE9GFCxdUKtXBgwfXrVv30UcfWVoO4L8rCwsiou6X0HLfZGdt/eQVb2964w0dLW/coB9/NGJ1YAoD+M8dDNXc3JyUlCSXy5VKpbOz8+rVqxcsWLB69epDhw4lJCRs27Zt06ZNCxcuFLpMgXBrPOvru2tTV0f02zUpx9+fvvhCR8uffkKMmj/cYgI9aDSaXbt2+fr6Jicnq1QqqVRaVlaWmJgYEBCQkZGRmZk5YcKEiooKiUTyyiuvDNDNSiZNIiIqLOyuzS+/PGkJfR9iFHoqJycnKCjonXfeqa2tnTNnzoULF3bt2jV06FBtg6ioqIsXL8rlcicnp6ysrClTpiQkJNy/f1/AmnsVtz1FVBQR0alTdO+e7mb5+XTtGllYUERE79UGpoQYhRe7ffv2kiVLIiMjCwsLR44cqVAojh8/Pnny5K4traysEhISrl27Fh8fr9FovvzySy8vry1btqi52cD+qrmZVq+m994jIvrd7ygoiFQqWrWKuu77o1LRypVERG+8QSNH9nadYCIM4PlaW1tlMhn33b92dnYymaytra2HxxYUFISFhXF/ZoGBgSdPnjRpqcJ4/Jh9/TVzc2NEzNKSVVYyxlh+PrO2ZkTs7bdZdfWTxleusFdfZURsyBB261bHixs3MiL22mu6+z96lBExLy8TfwzgBTEKumk0mrS0NE9PTyISiUSxsbE3btwwoJ/09HTuK0aISCKRXL9+3diVCuf4cTZ5MiNiRCw8nF28+OStw4eZoyMjYmIxmziRRUQwP7+OlsOHs4KCJy0Ro30fYhR0yMvLmzlzJpd9wcHBp0+f5tNbW1tbUlKSg4MDEdna2iYmJra0tBirVGFUVTGptCMWR41iCoWONvfusU8+YQEBzMaGETE7OxYSwjZsYM989u+/Z4GBbOVK3Sc6e5YFBrJFi4z/EcB4EKPwlDt37sTFxYnFYiJ66aWXUlNT1Wq1UXq+ffu2VCoViURENGLECIVCodFojNJzr2ptZTLZk2SUydijRy8+Sqk0fWUgGMQodFAqlXK53NHRkYisra3j4+Pv379v9LOcP39euyPUtGnTfv75Z6OfwlQ0GpaWxjw9GRETiVhsLLt5U+iawCwgRoExxtLT07VfMS+RSK5du2a6c2k0GoVCMWzYMG7WVSqV3r1713SnM468PDZjRscofupUxm+WA/oZxOhAV1paOm/ePC5A/fz8Dh8+3DvnffjwoUwmGzRoEBHZ29vLZLJHPRkd977qahYXx8TijrtDqanMSLMc0G8gRgeuhoaG+Ph47vl3FxcXuVyuUqn4dNjc3KzvIVeuXImNjeVC3NvbOy0tjU8BRqZUMrmcDR7MiJi1NYuPZyaY5YB+ADE6EKlUqtTUVDc3NyKytLSMi4urq6vj2WdDQ4Obm5tUKq2pqdH32Ozs7ICAAC5MIyMjL126xLMYI0hPZ2PHdoziJRJmylkO6OsQowNOVlbWxN92w5w7d25RUZFRut2/f7+VlRUROTk5paSkKPW8N80lu7u7uzbZa2trDavkyJEjCQkJK1asOHXqlCHHl5Sw6OiOAB0/nh05YlgZMHAgRgeQiooK7Qh63LhxRh9B8++/sbFRO88wZMgQA+YZ1q9f7+vr++9//zstLe3o0aN6HVtXV/fjp58yCwtGxNzc2Ndfs8eP9eoBBibE6IDQ0tLyzP2cX3/91UTnyszM5Hm1W1paOn/+fO1dr0OHDvXwwIaGBltb24qKCn3P2N7erp3lqAsJYXFxjPcsBwwciNF+Tq1WKxQKbh8msVhs2Nylvowy95qenu7l5aVdg3X16tUXHpKTk+Pn53fs2LE1a9Zs3bq1h7f+Dx8+7Ofnx51o3rx55aWl+pYKAxxitD87d+5caGgoFxAhISFnz57tzbNzKwEsLCwMXgnQ3t6ufSLAysrqhU8E7Nmzx9bW9v3339+1a9eCBQtef/317vs39SwHDBCIUbMWExPj7u6em5ur74G3bt0ykycvS0pKoqOj+axL5Z5P5eKYez718XOmLDMyMqZOncr9XF1dLRKJnrcfVedZDgcHB5POckC/hxg1X9999110dHRAQIBeT0y2trZq9wGxs7Mzk31A+D8llZ+fP2vWLK6HoKAgnXfhq6qqBg8e3NTUxBg7c+bM4MGDuwauILMc0L8hRs1UXV2dp6fn1atX9YrR9PT00aNHm+eudPyf2ef27uM+4PP27ktISAgMDFy1atXIkSNTUlKeebfzLEdoaOi5c+d4fSQAxhhi1GwtXrx4w4YNjLEexmhBQUHnizWz3SO5urqa5w5S3E7Stra2z7vcPnXq1I4dOwo67+n59CwHt4F/n9xfCswSYtQcZWRk+Pv7t7e3sx7EaF1dnfZOjpubm1wuf97UofnovJ/p1KlTDdjPtOexaJ6zHNCfIEbNkb+//4gRI4KDg4ODg21tbf38/Cq5b6d4mr43ss1K1931b+q/71z3g/TOkwDcLIdhG/gDdA8xao5u3rx57Te+vr779u3r+mxl52WVUVFRJSUlgpTKU9fvetJ3n6fn3TLqyS0pAKNAjJq7mTNnPrPgqbq6Oor7Fl8if3//n376SajajKWqqkoqlXKfyLCJy6amphUrVmgf6p8xYwY3/Tps2LCdO3caawN/AJ0Qo33PDz/84Obm5uzsnJSUpO8OIOas85c2h4eHX+z8DXE9o11OP3z4cEtLy741ywF9F2K07/H19SWifrlYhxuhe3h4aEfo9+7d07eTiIgIItq8ebMpKgToSkzQNzk5OQldgvGJxeIlS5aUl5cnJiZaWlru3r3b19c3OTlZqVT2vBNunpT7L0AvQIyC2eHmK4qKiiQSSXNz8+rVqydNmnTw4EGh6wLQDTEKZsrHx+fAgQOZmZn+/v4VFRUxMTGvvPJKcXGx0HUBPAsxCmYtKirqwoULKSkpTk5OWVlZb7zxhkajEboogKcgRsHcWVtbr1q1qqKiYtmyZRs2bOBWMgGYD0uhCwDoEQ8Pj3/9619CVwGgA/5h74cePXrU0tKiVquFLgRgQECM9kMSicTR0TEnJ0foQgAGBMQoAAAviFEAAF4QowAAvCBGAQB4QYwCAPCCGAUA4AUxCgDAC2IUAIAXxCgAAC+IUQAAXhCjAAC8IEYBAHhBjAIA8IIYBQDgBTEKAMALYhQAgBfEKAAAL4hRAABeEKMAALwgRgEAeEGMAgDwghgFAOAFMQoAwAtiFACAF8QoAAAviFEAAF4QowAAvCBGAQB4QYwCAPCCGAUA4AUxCgDAC2IUAIAXxCgAAC+IUQAAXhCjAAC8IEYBAHhBjAIA8GIpdAFgfGPHjq2trXVwcBC6EIABATHaD23fvl3oEgAGEAzqoc9Qq9UqlUroKgCehRiFviEnJyc4OHjTpk1CFwLwLMQomLvr16+/9dZbERERhYWFe/bs0Wg0QlcE8BTEKJivtra2devWTZgwYf/+/XZ2djKZ7Pz582Ix/mjBvOAWE5gjxti+ffs+/vjjqqoqkUgUGxv7j3/8Y/To0ULXBaADYhTMTn5+/vLly8+cOUNEwcHBW7ZsmTlzptBFATwXxkd91S+//CJ0CcZ39+7d999/PzQ09MyZMy+99FJqampubq5eGdra2lpSUkJE58+fN1mZAE9jYK5u37595cqVx48fP/N6bm4ut7ReIpFcv35diNKMT6lUyuVyR0dHIrK2to6Pj79//75ePWg0mrS0NG7gP2TIECIKDQ09d+6ciQoG0EKMmqPW1tZFixZNnDhxxowZPj4+lZWVz7yblJTEJamdnV1iYmJLS4tQpRpFenr62LFjuX/XJRLJtWvX9O0hPz9/1qxZXA9BQUGffvrp0KFDiUgsFkul0pqaGlOUDcBBjJqj4uLir776ivt56dKlH3/8cdc2t27dkkqlIpGIiEaOHKlQKDQaTe+WaQQlJSXz5s3j4s/Pz+/w4cP69nDnzp24uDgLCwsi4uYBuOv3lpYWmUw2aNAgInJwcJDJZL/++qsJPgEAYtS8qdXqBQsWbNy48XkNzp07FxoaysVQaGjo2bNne7M8PhoaGuLj47n4c3FxkcvlKpVKrx7a29u18wBWVlY65wEqKipiY2O5/z/jxo1LS0sz3icA6IAYNV8rVqx4+eWX586d232+qNVqhULRh8awKpUqNTXVzc2NiCwtLePi4urq6vTtJD093cvLSzsPcPXq1W4aZ2ZmTpw4kWs8d+7coqIiHuUDPAsxar4ePHhQVFQ0e/bslJSUFzbuK2PYrKwsnolWWlo6f/58rgdfX99Dhw51frelpWXHjh2fffZZRkZG59eNkt0AOiFGzd13330XERHRw8bmPIblX1tjY2N8fLylpSV3L17nPMCsWbPefPPN7du3jx8/Pjk5+Zl3+c8kAHSFGDVHJ0+e/OMf/1hZWVlTUzN//vzly5frdTj/Kz7j4n+lzF1Luru7a68la2truzY7derUiBEjuFtMubm5Q4cO1ZmSJSUl0dHRfO5rGRFyvB9AjJojlUr1+eefBwYG+vv7r1y5srW11YAezGEMa5R52+zs7ICAAC71IiMjL1269LyWX3755eLFi7W/2tradrOulv8qK55qamo+++wzd3f3NWvW9PKpwbgQo/2ZsGPYzqsIQkJCDFhFcOXKFe08gLe39wvnAWQy2bJly7S/Dhs2LD8/v5v2/Nf887F3794dO3Z88MEHq1ev7rWTgikgRvs//msz9cV/TevDhw+18wD29vYymezRo0cvPGrbtm1vvfUW97NarR40aNCtW7deeNTt27e11Q4fPvz03r2sF1fgrlmzBjHa1yFGB4reGcPyf8JKo9EoFIphw4YRkUgkkkqld+/e7eGxBQUF7u7u3MRrTk7O6NGj1Wp1D4/Ny8ubOXOmg7W1cswYFhzMTp/Wq2yDIUb7AcToAGLqMWx6erp2LzvDnvc/f/789OnTuR6mTZv2888/69tDTExMeHj4559/PmrUqB07duh1rEajKU5LY8OHMyImEjGplFVX61uAvhCj/QBidMCprq6Oi4vjNj/mnp7s+SVbNxQKBRd/U6dOPa3/pVznkfWIESMMfrZVpVKlpaWlpKQY/kBXayuTyZiNDSNidnZMJmM9mE8wGGK0H0CMDlDcGJYLvuDgYAOC7xlKpTIsLGznzp36hnJbW5t2HsDW1jYxMfHBgwc8izGCqiomlTIiRsRGjWIKhdHPoFKpGhsbV61atXz58sbGxr64JQJwEKMDF7eznKenJzcLGRsbe+PGjV6uIT09fcyYMXzmAUzr+HE2eXJHmIaHs4sXjdh3Xl7e2E6ampqM2Dn0JsToQNfa2iqTyWxsbLibQj28J87fhQsXwsLCuAANDAw8ceJEL5zUEGo1UyiYhwcjYmIxk0rZvXtC1wTmBTEKjDFWVVUllUq5UDP1tnv19fXa1ayurq5yubzr1tRmp6mJJSYya2tGxJydWVISUyqFrgnMBWIUnjh+/PjkyZO5MA0PD79o1DEs+21rOycnJ+3Wds3NzcY9hWmVl7OFCzvG+D4+7OBBoQsCs4AYhadwj296eHhoH9+8Z6QxbGZm5oQJE7iMjoqKKi4uNkq3AsjMZP7+HWEaFcUuXxa6IBAYYhR0aGpqSkxM5B4icnZ2TkpK4rPtXllZ2cKFC7kA9fHxOdgPLuLa25lczpycGBGzsmLx8azrZXVNDfvkEzZ5MrO3Z0Rs8GAWGsr+/vdnW+7cycaNYx99pPtEp0+zceNYVJRJPgUYCWIUnqu8vFwikWjj78CBA/r20NjYmJiYaG1trY1jZX+aUqyvZ/HxzMKCETFXVyaXM+0k74EDzMGh467UhAlsxgw2blzHBayHB+v8WMHGjYyIvfaa7lMcPcqImJeXyT8L8IAYhRfIzMz09/fXDsYv92wMa7rJAbNTUMDCwhgRs7Fh3Iqxs2eZlRUjevY5qPJy9uqrjIg5OLArVzpeRIz2fYhReLGut4a6X+R47NixSZMmcckbERFRWFjYa6UKQ6NhP/zANm/u+HniREbE3n1XR8v2djZnDiNi8+d3vIIY7fsQo9BTnRcqcdvudV2o1Hnh1KhRoxQmePjH3B0/zojYoEGsvl53g8LCjmf2ue+PQoz2fWIC6BlXV9ctW7bk5eXNnj27sbFx+fLlU6dOPXHiBPdua2vrunXrfHx8du/ezW1tV1FRsWTJEmFrFkB2NhFReDi5uupuMGkS+fkRY3TsWG/WBaZjKXQB0MdwTxwdOHAgISHh4sWL4eHhEolk3rx5X3zxRVVVFfdQaUpKCveM6UBUXExE9Nuchm5TplBZWUdLTksLlZToaFlVZdTiwCQQo2CImJiYqKiolJSU5OTkgwcPHjt2rK2tbfr06Vu2bAkJCRG6OkE1NRERubh014a7UG1sfPLKsWP023086HMQo2AgW1vbtWvXvvfee8nJyYGBgdbW1m+//Ta3092AJhYTETHWXRvuXQuLJ694elJUlI6W1dV09KgRqwNTQIwCLyNHjty6davQVZgTZ2eip680u2poIHr6ijUwkL79VkfLn35CjJo/3GICMCruq60vXequTWHhk5bQ9yFGAYwqIoKI6MSJ516QXr5MZWUkFlN4eC+WBSaEGAUwqjlzaMIEUirp0091vKvR0F/+QkQ0fz69/HIvlwYmghgFXrKzs998883w8PC//e1vSqVS6HLMgEhE33xDYjF98w3Fx1NLy5O3amvpD3+gI0fI3p7kcuFKBCPDLSYwXHl5+bJly7Zt2zZmzJilS5eKxeK1a9cKXZQZCAujPXvonXdo61basYNCQsjdne7cobw8UqloyBD673/J21voKsFoEKNgOF9f38uXL9vZ2RHRn/70p7S0NKErMhuLF9O0aZScTAcP0m8PepGXF73+Ov31r+Th8aSloyN5ej71Sme2tuTpSSNGmLxg4EHEul/gBtAza9euraur++abb4QuxPy0tVFzM7m4kI2N0KWASSBGwQiKi4sjIyNPnjzp6+srdC0AvQ23mICvysrKmJiYr776ChkKAxNiFHg5dOhQZGTk5s2bY2Njha4FQBgY1IPhMjIyFi1a5O3tze3n5OLisnfvXqGLAuhtiFEw3IMHD+rr67W/WlhYjB49WsB6AASBGAUA4AVzowAAvPwfxeM5Fmxgn+oAAACaelRYdHJka2l0UEtMIHJka2l0IDIwMjEuMDMuNQAAeJx7v2/tPQYg4AFiRgYI4IDiBkY2Bw2QAAsHQwaIZgQKgBnM+BlomrhBJjN5MDExOzCzeDCxsDqwsnkwsbEzsDE4iDCyMbCxsjAzicNsZ+B46KZmz8DgsB/ECcr6CWQfALOB4vthbKD4fiQ1ILY9khoY2x4mDjUHzBYDAPAHIqMgZqE5AAAAuXpUWHRNT0wgcmRraXQgMjAyMS4wMy41AAB4nKVSSwpCMQzc9xRzgVfSZz+va+tKVHDhHdx7f0w/pMpTKDSEdqYpw4REIcc9nZ8vSKxJKWArST8zxojHSkQq/yftbXQMFtLBEdU3rhKO+CfxmUXF6Gg3n1WMdvSlchtWWbqXncq4l6V7meioq2BGRTrCREd9RjNe+oz2XoZnZHjBytnIAbCVVNQIX04IIy+EURASypYaIaZVLsDpmtQb3w1xfgaidKgAAACPelRYdFNNSUxFUyByZGtpdCAyMDIxLjAzLjUAAHicXcxBDgIhDAXQq7jUBBo6UBjqHQavYIhLgzEsObzDKNOMXfX1/zRlzH0wndoZIbqZvNIIpK4GvIuEUWkDYbWWw5brvf7Lh//rx2/fdHR7dlH3Wp63d3mxgb4upT4AGQWGJ8HEVmDZCRyTgNgLPAdBaB+RJkHTCCo2swAAAABJRU5ErkJggg==\n",
"text/plain": [
"<rdkit.Chem.rdchem.Mol at 0x7f664718e8f0>"
]
},
"execution_count": 8,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"res[0]"
]
},
{
"cell_type": "code",
"execution_count": 9,
"id": "japanese-provision",
"metadata": {},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAJYAAACWCAIAAACzY+a1AAAZdUlEQVR4nO2deVxTV9rHn+RmIUDYBSMggguvgqCgr4KKOEopkqjUvaC1dbQure1HO2MXbWundqRaax3rSBdbRt9apYIRFQHBBQRl+RQRXBBBVMKekJCEbPee949gDAF3QnKd+/3wh3nOuec811/OvoSGEAIKMkO3tAMULwolIemhJCQ9lISkh5KQ9FASkh5KQtJDSUh6KAlJDyUh6aEkJD2UhKSHkpD0UBKSHkpC0kNJSHooCUkPJSHpoSQkPZSEpIeSkPRQEpIeSkLSQ0lIeigJSQ8lIemhJCQ9lISkh5KQ9FASkh5KQtLz3yWhWCzOzc1tamqytCN9yX+JhEijKb91a098/NTy8gNz5kTX1NRY2qU+g2FpB8yOWl0kk+3Q6Wrs7eGHH1gMxtU//xTdvr3ay2s7ixVkae/6gJe8FCqVaWLxOzpdV5lTKIjVq2tv31b5+7e2ta3s7MyyrHt9Au0lPmuv0Vxpa1sBQJjY9+5t9PBgzZ3rQqMxXV2TmcwRFnGvr3iZS6FMtstYv/v3Nbt3N1ZVqa5e7fT0ZAEAQtqOjt0W86+PeGlLoU53t6XlNRNjfn5HWZkiJMQuPJz7wEbz8Mig09362b0+5KUthVrtNeOPYrHu999bBwxgzJnjYqQfACCt9kY/+9a3vLQ9UoSUAHD/vubcOVl2dvu5czKtFo0cyWlt1R08OGz0aFtDTIJQWM7NPuDllLCsrOzo0fS0tOuVlUq9hcWiRUY6tLToWlq08+ZV/fLLUENZxDBXy3naC6mpqZmZmQCwbNmysLCwJ8a3uraQINo6OzO12usE0Y5hA5jMMRzODBrN9ikeJAoKCk6cOHHs2LGbN2/qjRwOffJkLp/vHB3t5OCAabVo3bo7QqGYxaLt2ePL57vQaGwPj6ynSb/fWL58+YYNG3g8np2dHYvFemJ8KyqF7723rrY2Xyqt2bLF09+f88B8vKNjt4PD3zic6F6fUqlU+fn56enpKSkpDQ0NeqObm1tMTMzMmeqwsGoWi2aIzGTS9u719fBg/vBD0+rVtVIpvmLFKqvSDwAaGhouXrxob2+/YMGCp4lvLaWQIPDc3JUBAVfy8jrOnpV++qmXSQQHh/V2dq8bPkokkjNnzqSnpwuFQplMpjf6+voKBIL58+eHh4fT6XQcb2ltjScIcc/svv++cevWehoNNm/+cMuWf5rvvZ6DkpISNpudm5vb0tLy5ZdfPjG+tZRClSo1IOAKAFRXdw4ezO4RTpPJdjGZgVIpLyMjIyUlJSsrS6PR6MNGjRo1f/58gUAQGhpq/AyGDXBx2SORvIfjLSbJrV07kMt1+fjjyi++2Nbervz222/pdEt0zjs64No1UKvBywv8/PS2MWPGMHQ6V5lsxT/+YRq/thZsbcHDo5sRWQEEoWlsjBKJQnNyRkVEONTWjhWJQo3/CgsDv/jCe8KEgTRaV62IYdikSZN27dp19+7dxyeO4+1S6TeNjdNEotD6+lCRKLSxcbpMtpsgOlJTU9lsNgAsWbJEq9X2z8t2IZGgpUsRg4EAuv4CAlBODkEQ8fHxa197LQbg3Natpk8NHIiWLDGxWYWEKlWhSBT6yy9D+Xzn69fHGJQ7e3bUhg28ESNsDF84DseGz+cnJSU1NTU9Yya6K1dOjxnjX1mZgxBusObk5HC5XAAQCARKpbJv3+uRqNVo/HjE5aKkJNTUhDo70cWLKCwMMZno7FmEUNupUzoAJBSaPmi1Esrl/1dZGRwR4RAb6ywQOG/a5FlfHxoR4WBQzsWFsWiR26+/Dm1vL3juXObOnQsA7u7upaWlxvaioiI3NzcAiIiIaG9vf+G3eQr27EEA6Pjxbka5HPn5ocBAhBDKz0fkkrCj4xeTmvPVV53c3Znu7syEBLfk5GF1dSF6u0p16blzUSgUMTExAGBvb5+VlWUcdO3aNW9vbwAIDAysr69/4Rd6EhERaMiQXuxffYUAUGXlM0loFd0ZDOMZf0QIiorkYrGuqGi0lxfrMTGfCVtb2+PHj69YseLXX38VCAQHDhyYP3++PmjkyJF5eXnR0dEVFRWTJ0/OysoaNmzYc2f0ZCorYdy4XuwjRwIAVFSApycAQEEBPOiydaFS9XzIKiRks8cDMAB0+o+1tSqxWOfhweyhnyeD4f0iGTEYjP3797u4uOzcuXPx4sXt7e0rVqzQB/n4+BQUFMTGxl66dGnKlCmnT58ODg5+kbweR0cHDBjQi13f1ZTJuiT87jvAsG4RlMqeD1nFNDed7mJrO9PwsaREAQChoXYm0eztEwBo8GLQaLRvvvlm27ZtOI6//fbbn3/+uSHIxcUlOzv7lVdeaWxsjIyMzM/P7zWFkpKSRwU9gYoK2LoV5HLgcqG5uZcI+k09Dg86AYcPg1ze7c9kOAEAVlIKAYDLfU+tLsXxegAoLdVLaP8gkAaA2OwJtrZz+yq7jRs3Ojg4vPPOO1u2bJFIJIZxob29fXp6ekJCQkpKSnR09JEjR2JjYwFAocDr6tQSiba5uXHHjo/s7dmZmSeeKieCgD//hPR0OHIErl8HAPD3h4AAuHWrl8j6ecHAQJBInuFlzNNePw86nailZbFIFDpyJAcAhEJ/Q+/mP/+ZVV5eguP4k1N5FlJTU21sbAAgISFBo9EYeaJbuXIlAGAY9uOPP128KP3xR1FSkigpSTRu3KyPP84ICoo6ebJNJtM9MmmVCp06hVauRB4eD0d+7u7or39FpaVo924EgE6fNn1kxAgUGIgIgnw9UiO0TU0HMYzGZNJqasY2NExua3tHLM5ksVh0Ol0qlfZ5foZxIZ/PNx4XEgTx2WefAQCNRps7d/MPP4iSkkQrVyYJBB/s3n0rKCgqKUmUnNwoFnebEJBLpejQIbRgAeJyHyrn54fWr0cXLiDdA8k7O9Ho0cjVFf3xB9J/dW7eRNHRCMNQdjZCJBxUGJOdnQ0AEyZMIIgOvSUvLw8AgoODzZRjUVHRgAEDAGDixImtra3GQevWJdJodACIjl6blCRyduaFhS0ID1/o4uL5978Lk5JEv//erNUSLS0tycnJfD6fzWZrvb27lBs1Cm3ciPLyEEH0kmtTE5o1CwEgJhPZ2yMA5O2N0tO7Qkk3qDCmsLAQAMLCwmg0exOLmXIcP378+fPno6OjL126NHXq1MzMTE9PTwCQSHQBAUvefNM2OXl9Zub3HR2tX311mU5nqNWKn35aM3To+KammszM0199lVVRUUIQBABgGHY1OnpsYCDMmQM+Po/L1d0dhEK4fx/KykCthsGDISTkYf9z7FgoKwNfX9OncnKAyzU19tFXuc/Qj74PHz5ssMTFxQFAcnKyWfOtr68PDAwEAF9f36qqKoRQcbFM3/6tXZvMZNoAwJgxr+7ZU7t585nw8EU+Pg/3oNrYdE37NTQ0mNXJXrGWxSY9CCE3NzexWHz37l39dAkAeHp6ikSiqqqq4cOHmzV3sVjM5/MLCwsHDhyYkZHR1ORTV9c1lL516/LevW8olTIbGzuVqmujhp2d0+jRUSEhMdu3L7S1tdiio3VJeOPGjZEjRw4aNKi+vl5vuXPnjq+vr6ura0tLi2GZwnx0dHTExcXl5OQ4OTlt2nSYyx1tCLp3r3L79tkYxmKxOKNHTw8KigoIiMQwJgCsWMEzv2uPxLraQn2zFx4ebmIJCwvrB/0AgMvlnjp1asmSJRcvXvT0HPhgLRkAAMMYarXSyckxMbHU+BFbW7oF9QPrlNC452LuvkxPWCzWb7/9JhKJOjvdzp1rN9hrakoBYNiw8SbxBw3quUDdr1jFBJuBnoIVFBRA93LZD2AY5u3t7ednY2f3cIqypqYEAPz8TKenAwNNJwL7GSuSUCaTXb9+nc1mjx07Vm/p7OwsLy/HMGxcr/P6ZobBoE2Z4mioJPWlcOjQbns7goLs3N2Z/e+bMVYk4aVLl3AcDwkJ0U96AUBxcbFWqw0KCrK3t3/8s2Zi8GD2tGlODAZNqZQ1Nt5mMtleXoGG0IAAuwkTHB7zeP9gRW2hNTSEPRk2jOPuzkxKuogQMXhwEIPBpNNh0CB2cLCdp6eFW0E9lIRPxsGBIZeXAwCfH7F4sbutLR3DLNoH7Y61SIgQKioqAoCJEycajJcvXwYrkBAefJkiIsK5XOyJkfub/p8Q6pXKykoA8Pb2Nliqq6sBwN3d3YJe6cFx3NHREQDu379vaV96wVq6M1ZbiwLAtWvXpFLpkCFD9NPf1obFKlKCIIqKcgICqlSq8zh+Lze3BgBCQgiCENPpLmBNEurHptbgSa9YphTW1tYuXRozb55ALv9Fp6tBSFtSIgOAwMCK5ubXVKpzYE0SWo8nvWIZCXm8u9u3S3x82ABRAOtkMp/qahWLRQsMtEVILpH8XSzOvnr1KpPJtMig3gQrl9ACFSlByKTSbwDmAfwBsA4ASksZBJEaHMx9cJCMqKr6fNSo/7GxsbXgIo4eiURSVVXF4XCCgqz0kpr+lhDH8fb26wSxC4ADINQbS0tLAWDcuIfLgT4+6oKCj9jspzpgZ1YKCwsRQuPGjXua05oWof8qUq1WK5FImpubNZrBCNmUl5fL5fLy8nKtVltSUgIAoaGRxvHV6osMhuWHrVZei0L/SKhSqdra2lpbW1UP9pMjhMrLyxMSEsrLy3EcxzAMwzAXl1CAh1PGOC7qB9+eiPVLaMZVe4RQZ2enQqHQ6XSPj7l///5NmzY5OzsfODA1JOSe3shgDBkw4A8z+faUEATh7Owsk8kaGhoGDhxoWWcehVlKIUEQcrm8ublZKpU+UT8ASEhIiI2NlUgkCxZknD3btVLOYPiZw7dnoqKiQiaT+fr6Wq1+8CgJcRzX77t91uRwHJfJZM3NzR0dHfp9eU8Di8Xat2/f4sWLlUr1smXVQqEYAGxs/vKsufc5FllwflZ6l7C9vT0kJMQgYVJS0oQJE7RarXGc6urqCRMmZGV1XSZo6K0oFIrn0B7DsB07PlqzxkurRWvW1P78M87hRD1rIn2O9TeE0FNClUr1/vvvSySSoKCg8vLyjz/+GACkUum9e/dMYmo0mnv37ikUip69lWdHC3CLRvt+0yaPL77wBoDNm8s+/PCT503thcDxho6Ova2ty5qbY/PzUwEgNJQFYEX7/Eww7bVrNJrAwMDt27dXVlYePHgwODhYY3JKsTtKpVLyTMdwAACATqczGDocv4DjxQB3AOoBcH3Q6tWRXl5jVq7cmJiYKJfLd+/e3Q93UWg0mra2Nh6Pp1Ac6uj4F0IaAJBIdLW1cg6H7um5r62tyMnpKwzr7VCgpTGV0MHBISEhIScn5+23375z5058fPzj/wefcujGeACTydSfcQEAgIVqtbdKdUGnGwBAYzAG29hMZbPD3nyT7uo6bNGiRd9//71YLE5OTmYyzbI/RalU5uTkpKSkpKenT5ky5bff5nR07DWElpQoEIIxY+wYDFp7e2lOzuwBAz6ZOjXGMtebPJpeBDh+/DiPx1u/fv2cOXPKy8vHjBmjt8vlcmPBlEolADCZTBrNdGRCp9P1ajGZTL1yj9gFSmezJ7PZk3sGzJo1KyMjY/bs2YcOHWpubk5LS+P2PEvwvDQ3NwuFQqFQmJOTY6j8m5rqpNJ9xuqUlsoBIDTUTqkk3nijet481/r6PVOmRJNAwlmzZo0YMUKhUBw9etRYs4CAgJ6RaTQai8XCcdxYMwzrg6XtqVOn5ubmxsTE5OTkTJ8+/dSpU/p7KZ6bO3fuCIXCEydOnDt3Tj/UodPpoaGhfD5/8eLF7u4/qVRnjePrT6qOG2d3/Lh45kznyEgHN7cWABHA4Bdxo88xlVAul8fGxsbExGRkZBw8eNBwsAEA0tLSjBWtra1dt24dALi4uJjJuZCQkMLCwldeeaW4uDgiIiIrK8vLy/RirydSWVmZkpJy4sQJ/UwsANjY2ERGRvL5/AULFvB4+vsXdI2NBcZP4TgqK1MAwNixdj//3FxR0dnYqCkpUfz++6nBg1e9+Kv1IaYSCoXCWbNmbdiwYcSIEcnJyZs2bTIEhYSEGLdJHA4HzI+fn19eXl5MTMyVK1cmT56cmZnp7+//xKdwHC8sLDxx4kRqauqtB0einZycoqKi+Hz+nDlzHBwcusdvRkit/7dUil+4IDt2TKzTgZMT5uTEsLPD4uPdYmKc/vWvxvz8i6+/bt0S1tXVjRo1CgD8/PwyMjIs4ZIpPB7v7NmzAoHg4sWL4eHhJ0+eNN4iZYxGozl58qRQKExPTxeLu27P8/T0nD179pw5cyIjIx/TLRKJNJmZ0tOn2wsLO3Q6BAAcDr29HV+6tHrjRs+vv663taXn5koTEp65GjA3phI6OjrqBwnt7e3Ozs6WcKkXnJ2ds7Ky5s2bl5GRERUVlZqaGhVlPPDHCUJBo7F1OjwhIUHfz/Lz8+Pz+fPnz580adJjjtTU1NSkp6enpBwpKLiq75NhGG38eHuBwNnf32bdujvnz8ukUvyzz7xqa1XffjvEx2f0o5KyGCbboWpqaqZNm3bz5s0FCxbk5+frjYmJiTwez/hCAYRQZWUlj8dLTU018wath2i12jfffBMA2Gz2kSNHEEKdnWfb2lY3NISJRKEi0bjm5vnr18f+85//uHHjxmPSwXG8pKTks88+G6m/qgcAAGxsGFFRjl9/7VNeHmy4puHSpcAhQ9gAMHy4TUnJaJFovE53r79e92npZaXi8uXLGRkZkyZNMnzTi4uLi4uLV61aZdyfbmtrO3z4cFRUlLlPbhqDEPrggw927tyJYdiOHdMXLmzrGQfDPF1cdjIYQ03s+gYyJSXl6NGjhvOLrq6uM2fOFAgEM2YMU6lWARAmEzHNzdr4+OrKSqWHBzMt7a2wsH1merXnxrqOiD4liYnbPvzwIxoN1q8ftGFDL1d60ekubm4HMMwDjMbvx48fl0ql+gg+Pj6zZ88WCARTp041NJAKxUGZbFfP1GQyfOnS6qIiuYuL84kTJ61typSUEiqVx3fvfveTT+4SBCxf7r5li3fP0bZCMb6gICwlJSU7O1ut7uptjho1SiAQ8Pn8RzWQSmWqTPYtQp09QkLXrKk/duyknZ3d0aNHo6N7v2XaIpBSwpaWhTrd7YyM9jVratVqYu5cl507hzCZNAC4e1edlSXNzpYaOpZ0Oj0sLEwgEMTFxY0Y8eTf9iGINqUyTa0uwfFGOp3DYPhzONFsdphWq12+fPmBAwdcudzaAwe4s2eb/T2fEks2xM8FjrcYuhtHjoywt8cAIDycu3btQKNL2cHGhh4bO27//v0tLS19lTVBEOvff/9eWBjCMLRvX18l+4KQrxRqtZWtrW8YPv75pyIhoRrDaK2tWgBwdMQiIhxmzHCMiXEaOHA5l/tO33uQmAgffQQIwcaNsG1b36f/jFh+i9iz0214Pnas3bFj/lIpnpYmjolxmjjRnsHoauRoNPNsG9y4ETw8YMUKSEwEhQK++w4sOvFNvlKIkLqpaZp+SQ8ApFK8okI5bJiNh4fpzIuz89dm3L0hFMKiRaBSweuvw6+/gnmWw54G61o3eRpoNDabPVX/b5kMT0i4VVenfvfd2hs3unUj6XQum23O3v/s2ZCRAQ4O8NtvMHMmyOVmzOuxkE9CAOByV+orSYUCZzBo8+e7BgfbtbYatsrRAMDefjmNZuaJ+MhIyM2FAQPgzBmYPh3aepln6AdIKSGD4evouAmAPnAgi8vFFi26delSx//+r+FKBcThvGL80zJmJDQU8vJg8GAoKoIZM6D7DrH+gZQSAgCHM9PF5bvz57GAANujR0fExbns2dMIADQam8t928npy/57NX9/uHQJgoLgrbe6WsTaWvjgA5gyBcaOhRkzYNs2eDArBGo1zJwJhw6ZJvLpp7B27fPlT8YeaRdsdpin57a6us8lkskSyVlfX3dHx9dtbGboT5j2KzweFBUBmw0AkJkJr70Grq6wcCG4ucHt27B1K/z733DmDAwfDjodZGRAz8Wy4mJ4MG37rJBYQgCIiJim06HU1IIpUzbOtux0iV6/1lZYvBiCgiA7Gwx35fztbzBpEsTHw+XLZsna0nMLLxfffIMAUEmJqf277xAAystDcjkCQFu2mEZ49VU0evTz5UnWttBKyc8HW1vo/gNuAACRkV2hZoDcFanV0dgI3r39GIr+kucHv5AJO3ZAUlK3CG1t8BRT8L1CSdin0OnQ60kuvdGw/2/6dBAIukXYteu586Qk7FO8vKCkBBACk8XIurquUD1jx8Jbb3WLkJLy3D1Sqi3sUyIjQa2GCxdM7ZmZAAB/McuELSVhnxIfD4MGwbvvQmvrQ2NZGSQmQnR0M4+XkpaWApBSWXnt2rW+ypOqSPsULhdSUkAggJEjIS4O3N3h9m1ITYXhw2H/fkOsk1VVyuJi/X7dF4d8i00koKUFkpKgoACkUvDwgOhoWLasa+yv0cCqVdFlZcKCAsPNuQAAu3aBRAJbtjxHbpSE/U1WVtbly5c3b97cVwlSEvY3cXFxP/744wue0jKG6s70KxUVFTwerw/1A6oU9jM3btxwdHR8cCKub6AkJD1URUp6KAlJDyUh6aEkJD2UhKSHkpD0UBKSHkpC0kNJSHooCUkPJSHpoSQkPZSEpIeSkPRQEpIeSkLSQ0lIeigJSQ8lIemhJCQ9lISkh5KQ9FASkh5KQtJDSUh6KAlJDyUh6aEkJD2UhKSHkpD0/D82/63X90kjDQAAAABJRU5ErkJggg==\n",
"text/plain": [
"<PIL.PngImagePlugin.PngImageFile image mode=RGB size=150x150 at 0x7F664719CE90>"
]
},
"execution_count": 9,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"Draw.DrawMorganBit(epinephrine, 266, bi)"
]
},
{
"cell_type": "code",
"execution_count": null,
"id": "dress-juvenile",
"metadata": {},
"outputs": [],
"source": []
}
],
"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.9"
}
},
"nbformat": 4,
"nbformat_minor": 5
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment