Last active
October 25, 2020 14:46
-
-
Save adelenelai/59a8794e1f030941c19bcb50aa8adf3f to your computer and use it in GitHub Desktop.
20201024_inchikey_cis_trans
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
{ | |
"cells": [ | |
{ | |
"cell_type": "markdown", | |
"metadata": {}, | |
"source": [ | |
"# Issue: Nitrogen sp2 isomers get the same InChI Key\n", | |
"\n", | |
"https://sourceforge.net/p/rdkit/mailman/message/37135337/" | |
] | |
}, | |
{ | |
"cell_type": "code", | |
"execution_count": 1, | |
"metadata": {}, | |
"outputs": [], | |
"source": [ | |
"import rdkit.Chem as Chem\n", | |
"import rdkit.Chem.Draw as Draw" | |
] | |
}, | |
{ | |
"cell_type": "code", | |
"execution_count": 2, | |
"metadata": {}, | |
"outputs": [], | |
"source": [ | |
"m1_cis = Chem.MolFromSmiles(\"C/N=C(/NC#N)NCCSCc1nc[nH]c1C\")\n", | |
"m1_trans = Chem.MolFromSmiles(\"C/N=C(\\\\NC#N)NCCSCc1nc[nH]c1C\")" | |
] | |
}, | |
{ | |
"cell_type": "code", | |
"execution_count": 3, | |
"metadata": {}, | |
"outputs": [ | |
{ | |
"data": { | |
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAlgAAADICAIAAAC7/QjhAAAn6UlEQVR4nO3daVwUV7oG8LfZRZRFFFFxRcXdAFGJu6JJhFFvIppxj0s7WcTEMelohmBcZsCog2bmZojXOLhExWRikBgzRESN4tKIigqCimxugCBrQzf93g/VtmhQWaq7wXr+v3yAprrOiynqOefUqWoZMxMAAIBUmZm6AAAAAFNCEAIAgKQhCAEAQNIQhAAAIGkIQgAAkDQEIQAASBqCEAAAJA1BCAAAkoYgBAAASUMQAgCApCEIAQBA0hCEAAAgaQhCAACQNAQhAABIGoIQAAAkDUEIAACShiAEAABJQxACAICkIQgBAEDSEIQAACBpCEIAAJA0BCEAAEgaghAAACQNQQgAAJKGIAQAAElDEAIAgKQhCAEAQNIQhAAAIGkIQgAAkDQEIQAASBqCEAAAJA1BCAAAkoYgBAAASUMQAgCApCEIAQBA0hCEAAAgaQhCAACQNAQhAABIGoIQAAAkDUEIAACShiAEAABJQxACAICkIQgBAEDSEIQAACBpCEIAAJA0BCEAAEgaghAAACQNQQgAAJKGIAQAAElDEAIAgKQhCAEAQNIQhAAAIGkIQgAAkDQEIQAASBqCEAAAJA1BCAAAkoYgBAAASUMQAgCApCEIAQBA0hCEAAAgaXUOQpWKSksNUQkAAIAJPCsIS0qovJyISKOhkhLdi+vWUc+etH07MRu+OgAAAAN7VhCuW0d+fkREKSkUHExExEyHD1NODs2ZQ2PH0qVLRqkRAADAYJ4zNWphQdu3P/pWJqO4OIqIoDZt6MgReuklWrSI8vMNWyIAAIDhPCcIV6ygLVt0URcfT8ePk0xGs2dTSgotXkxE9PXX1Ls37d59W6vVGr5aAAAAkT0nCG1saO1aUihIq6VFi2jECPrDHygzkxwdafNmSkqi8eOprIw//NBv0KBB8fHxxikaXjw3bpBabeoiAECSnr9qdMQI6tePmCkggGxtKTqa+vSh0FCqrCQPD/rlF9qxI9nSMjchIWHYsGELFiy4d++eEeqGpuv3i7CqqugPf6Du3bEICwBMgesiO5tnzWIiJmJ3d46O1r1eWloaHBxsY2NDRA4ODiEhIRUVFXXaM0hHUBCPHs3MnJTES5cyM6ens4eH7rgaNIiPHDFpfQAgMbW6j1AupyVL6MEDat+etm+nX3+l3r3p2jXy96dJk+j6dba1tV25cmVSUpKfn19hYeEnn3wyaNCgiooKA4c4NFVPLMLq3JkuX6aICGrXjs6codGjadw4On/eZOUBgKQ8PwgzMmjbNtq8mTw8dDNXY8fShQsUFkb29hQVRTNmBC9ZsqS4uNjd3T06OjomJqZXr15Dhw61trY2wi8ATVH1RVh799K5c2RmRrNnU1oarVlDLVvSr7+SlxctW/ZNVlaWqYsFgBddbYaNiYk8fLhu5srbm+Pjda/fucMKxU2ZTEZEHTt2jIyMFF5XqVRFRUUGGsNCUxcUxPHxfPQoDx7M77zD9vYsk3FAAKel6TbIz2eFgn187hKRlZWVXC6/d++eSUsGgBdZraZGBw6ko0dp505q146USho6lBYtorw8tYsLhYR0Onv2rI+PT2Zm5tSpU0ePHp2UlGRtbd2iRQvDBrgk5efnb9u2zdRViENYhGVmRnI5WVvTvn3Upw998AHl5ZGTE4WE0K5dZdOnT1er1V9//XXPnj2/+OILlUpl6qpfTFqt9sKFC6auAsB06hSbJSUcHMzW1tyunaZdu+5hYWFqtZqZtVptREREmzZtiMjCwkIul+fl5RkmuaVLo9EMGTKEiJYvX67Vak1djpiyslguZ3NzJmI7O1YoWD+hcOnSJX9/f+FY7dChQ3h4uEajMWmxL6CPPvrIwsLiq6++MnUhAKZRtyAUXL3KCxfqljoMGDDg2LFjwuv3799///33LSwsiMjFxWX37t1VVVWiVit1O3futLS0JKI5c+YIXZAmbcECVii4sFD37cWLPGGCbga+XTvevj1b/zvGxMR4enoSkYeHxwvwizc2y5YtE/6cX7w+FkBt1CcIBVFRUV27dhX+fvz9/TMyMoTXk5OTx48f37dvX5lM5uXldfLkSZFKBWbm6OhoW1tbIpo0aVJZWZmpy6m/69fZzIyJuE0b/vJLrqzUvX7yJA8bxlZW7OY2smfPnpGRkcKpuaqqaufOnTExMaYs+sX1f//3f0If68033ywtLTV1OQBGVf8gZOby8vJVq1YJ52U7O7vNmzcLr2u12n379nXo0IGIzMzM5s+ff/fuXTGqBWbm06dPt2rViohGjRr14MEDU5dTf6dO8ciRuiFgp04cHs7CDIJWyz/+mN69e3ehm/XKK6/89ttvpi72xRcTE+Pg4EBEgwcPvn37tqnLATCeBgWhIDs7e9asWUT0xRdfVH8dd9kbzuXLl4V+hpeXV1PvZOzfz7166eLw5Zf52LFs4fXKysrw8PC2bdsKcejr63vhwgXTlvrCS01NFfof7du3T0xMNHU5AEYiQhAK4uLiasy5tLQ0/WKHnj17Hjp0SKwWJS49Pb1Hjx5E1LVr12vXrpm6nAapquLISO7UiV1cNLa2zr6+vvqzcHFxcXBwsJ2dnbAOa+HChdnZ2SYt9gWXl5c3fPhwYY4nWv/sKIAXmmhB+GwHDhxwd3cnok6d3P/4R3VmpnGafcHduXPnpZdeIiJXV9cXYLRUWspffnnc3t5emFGfM2eO/sJzbm6uQqEQHtHg7e2tUCgKCgpMWuyLTKVSzZgxg4jMzc2//PJLsXZbXMzCRW21mouLuaSEhZ6zWs246xhMy0hByMwqlWrt2rXDhh0k4ubNee1aVqmM1vgLq7i4eNy4cUTk6Ogo1oU0056w8vPzFQqFMKP+xN30qampixcvFh7g4Ozs/Pe//x2T7Qai1WqDg4OFf+rAwEBRbll54hmzq1bxL78wMyckcGBgw3cPUH/GC0JBTg7PmsUyGRNxt2788Fk0UH8qlWrKlClEZGtre/DgwYbvsDGcsG7cuDF9+nQzMzPh9sHqgXfmzJnRo0cLk+1ubm64s9Bw9uzZI/RIXnvttYYsy6qq4vPnedkyHjeOIyIeHVc//cSVlXzqFIIQTKxWT5YRUbt2tH07HTlC/fvT9es0dSqNG0fJyUau4oVibW29Z8+ehQsXlpWVTZw4sSGPnsnIoJgYot89FFujIbXaqJ8X2KVLl127dl28eNHf33/evHlWVlb6H7388suxsbExMTEDBw7MyspatGjRgAED9u3bZ7ziJGPatGmHDx9u3br1oUOHhg8fnpmZWfv3ajSUkECbNtHUqdS6NQ0cSCkpjz1jlojWraPZsykoyCDFA9SBqRJYreZNm9jBgYnY2ppXrcL1gwYR5rKISCaTbdiwoZbvysnhqCgODmZ/f3ZxYSK2sODly/nIER42jOPidD33kSP5rbd43DjT9NyfNuDTaDTbtm1zc3MTjuQFC2LOnTNyaZJw7do1Dw8PImrXrp1SqXzGlqWlpXFxxStX8pgxbGurWwks/NelC/v5PXrGLKZGoVExWRAK8vM5MJDNzdnPz/TTcS+ATZs2CdOJCoWixg1u3bp14MCBv/1ts78/t2372KmKiJ2d+bXX+IMPmtIJq6ysLDQ0dNAgX5lMKzy8u4kvoW2M8vPzR40aRUTNmzffv39/9R8VFxfHxMQEBwf7+vpaW1uPGPEP/eHUtSvPmsXh4XzjBvPDh60z84IF/OGHPHt2oz6uQFJk3Ag+Efz8efr+ezp9mmbOJE9P2raNHBzIy4vGjaNz5+jbb2nTJlOX2HTs2LFj/vz5arX6nXfe+cc//lFUVHTp0qWEh65cuUJEZmZmtraFJSUtWrakfv3Iy0v3X+/eJJPRZ5/RhAk0ZAgtXEjm5mRvT2PH0vjxdO4cRUQ00v8XBQUUGkqbN1N5OVla0ttv06pV5OJi6rJeIJWVlXK5PCIiwtzcfPny5QMGDDh+/PjRo0eTkpK0Wq2wjbm5+ZQpy1xdQ0aMoGHDqHXrmnfFTKNG0fHjdPgwPbzaC2BKjSIIieizz2jMGAoKojVrKCqKHBzo8GFydaX8fOrVq5GefButqKiot956q7y83NHRsaCgoPqP7O3tPT09vb29fXwU/fu36tbtWfvZv59mzaL/+Z/Hrhc2ZtnZtHo1bd1KVVVkZ0fvvUeffkr4HBSxMPNf//rXoKAge3v7wsJC4UULC4sBAwYMHTp02LBhvr6+jo6OtdnV2rX0l79Qjx508SLhc0vB5BpREE6YQJWV9PHHNHQoOTjQ4MGNfRTSmB0+fPiNN95QqVRWVlYDBgzweqhXr17C3GltpKdT795UUUGnT9PLLxu0XjFdvEiffEI//0xE1K4d7dtHL71EzZqRRkMqFclkZGlJVlak0VB5OWKyzqZNmxYZGdm1a9e5c+eOHDly0KBBwsrSOqmspIEDKTmZ1q6lFSsMUSZAHViYuoDHCJ9RBw3XqVOnoqKijh07pqen1z75ntClCy1ZQqGhtGQJnThBMpm4NRpK//508CCdPEkKBaWn04ED9Je/UGwspaToZt2FPtbFi+hj1Yfw0LvAwMAlS5bUeydWVvSvf9GoUbRmDb31Fj18ej+AaRj79omnEa7oODnRuXO0YQMFBdH48UREnp44VdWHsNK9S5cu9U5BwaefUtu2FB9P338vUmXG8sordOwYnThBlpamvxvkRZKXl0dEzs7ODdzPiBE0YwaVl9O774pRFkADNJYgJCK1mgoKqKTE1HW8ELKysohIuLWguLg4MTExX3/3Vl20aEErVxIRffwxVVSIWaERyGTUqRMR4fY1MVUPwm+//Xb+/PmxsbH129X69eToSL/8Qv/5j5gVAtRVIwpCYelZwwYwoCMEYceOHYno9OnTnp6eU6dOrd+uFiygfv0oPZ02bxazQmOysaG1a0mh0H27YgXt3k0hISatqcmqHoTHjh375ptvUlNT67crFxf661+JiJYsoeJirXg1AtRNI4qdqioiInNzU9fxQqg+IhSmSfU3nteVuTmFhRERrV5Nd++KVaBRlZSQnR0uRIlDCELhEzEbPk0ql9Orr1Z07vyvVasUz98awDAaURBiRCii6uFXfXRYP2PG0IQJ1K1bQVjYFrEqNJpVqyg5mby8qFUrXH4WgTDHLoRf9VCsHzMzWrfu6qlTi8PCws6fPy9SjQB104hiByNCEVUPv+qjw3rbsCH70qU269e/e/nyZVEqNCbhTrUmd42zESovLy8tLbWxsRE+IbJ6KNZb//793333XY1Gs2jRIv29+QDG1IiCECNCEYk4NSrw8OiwaNEijUbTkEXzpoIgFMsTc6FirSBds2ZN+/btz5w588033zRwVwD10IhiByNCsRQWFhYVFdnZ2Tk4OJBII0IiWrVqlZOT0+HDhw8dOtTwIo0JQSiW6snHzMKIsCFTo4IWLVps2LCBiBQKRW5uboPLBKibRhSEGBGK5YmLgtnZ2SRGEDo5Oa1YsYKIli5dqtFoGrg3YxKCsLLS1HU0fdWD8MGDB2q12t7evvqHZNXbtGnTJkyYcP/+fYUCq2bA2BpR7Gi18TKZmZ3dZFMX0uRVD8L79++XlJQ4ODi0bNmy4XtevHhx9+7dk5OTt2xpSqtmMCIUS/UgFGVeVKPR6CcYNm3aZGNj8+9///vIkSMNrhSgDhpREFZVVTFzRUV97vuG6qpfFBTlAiE9XBZhZWUVEhJCRMHBwQ8ePGhoocaCIBSL6EG4ePHi119/fd26dUTk7u6+cOFCMzOzlJQUMYoFqK3GFYREZI6LhA1W/aJgw++dIKKvvvqqR48e8fHxRPTGG28MHz48Nzd3/fr1YhRrDAhCsVS/KNjwINy4ceO//vUvGxub4cOHCzs/ePCgcB4AMKZGFITCyukGPhsT6PHwa/iIcO/eve+//35BQYHwABFmdnBwsLOzmzBhgkj1GlyzZqpevcpatWoyQ9hG6/cjwnqvlDl48ODHH38sk8m2bt3q4+OjVqsDAgKuX7/u6ek5e/ZsEWsGeK5GlDoYEYrl93fT1zsIDx8+PGfOHK1WGxoaOmfOHCJatmzZgQMHZDKZpaWleCUblqVlanJy8/T0EaYupMmrHoQNuYkwMTFx2rRpVVVVa9asmT59OjPPnz//yJEj7dq1+/HHH5s3by5u2QDP1og+hgkjQrGIdTf9mTNnJk+eXFFRoVAoPvroIyJas2bNxo0brays9u3b5+3tLWrVBmRtbU1EFZgbbbDqo8B63ztx6xa9807LNm06vPnmYGEd8urVq3fs2GFnZ/fTTz916NBB7KoBnqMRBSFGhKJg5pycHJlM1r59e2rANcIrV2jlynOlpaXz5s3729/+RkTh4eFBQUFmZmY7d+589dVXRa/ccIT1/QjChqs+CiwuLqa6jwhLSsjfnxITu/n5nfj6azsi2rdv38qVK83MzHbt2jVw4EDxiwZ4Lm40oqOjicjPz8/UhYhPq1WpVKnGaevWrVtE1KZNG+HbqKio0NDQu3fv1mknN29yhw5MxO+996tarWbmH374wdzcXCaThYeHi1+0geXk5BCRq6urqQtp8oThWmZmpvCtSqWqqKio/durqnjSJCbibt343j1m5lOnUoQPuN+0aZMhCgaojcYShOXl5cIVcjc3t5s3bxq6uZKSksOHDx89erSsrMzQbZWWJqSmvpaZuTg19VVDt8XMp0+fJiIvL6967yE3lz08mIhHjuTycmbm2NhYYXZx7dq1ohVqRMKEnpOTk6kLEZ8x+1jMLIRWvf9qPviAidjJia9eZWZOT+c2bbQ+Pqvff3+xmFUC1JHpg1Cr1e7du7dz585E1KxZMyKytbX9/PPPy4VzsAGai4yM7NSpk42Nja2tbYcOHSIiIgzRkN7du5vu3Fmn1aqLin7VaO4btK2Kioq5c+cSUceOHS9evFiPPRQVsZcXE3G/fnz/PjPz2bP88sunnJxavffeeyKXa0ilpaX6wYowide8eXPTliQ6Y/axKisrN2zY4ODgYGlpGRAQcOLEibru4dQpJmJraz56lJn5wQPu25eJ+LXXWK0Wv2CA2jNxEJ49e1a4hYiIevXqtWPHjlmzZslkMiIyRETFxcV5enoKzfXt27dXr17C1+PGjbty5Yq4bTHzgwc/Z2V9pNEUZWV9mJzsc/361MrK29nZH2s0haK3xczfffddt27diEh4xCgRDR06NCoqSqvV1n4n8+YxEbu78507zMwpKdy6NRPx+++n12k/phUVFdW5c+c1a9YI31ZWVhKRpaWlcVrXarXC0yEMzWh9rKNHj/bu3Vs4qIQ/TyIaPHjwv//97zp1WLdt4x07mJkrK9nXl4m4d28uNMhfA0AdmCwIs7Oz5XK5sEbU2dk5LCxMo9EIPzpz5sygQYOEP7YxY8ZcunSp4c2lpaUFBAQI+2zfvn14eHhVVVVVVVVERETr1q2Fs2RgYGChSH+UZWWXUlNfVSpJqaSsrD+r1feYOStraWrqa0olXbjgkpf3DbNo58qEhISRI0cKv13Pnj23bt2qUCjs7e2FVwYOHBgREaGuXa/79m2eMIGvX2dmzs7mzp2ZiP38uLJSrGIN6+bNm5MmTdL3A/ThLRxpRsino0eP9u/fXy6X+/r6njt3zkCtaLXq3NytanW+vo+l1aqzsj4qLj4mbkP5+flyuVwIP3d390OHDqWlpSkUCv0aGXt7e7lcfvny5Trt9k9/YiJu25YzMsStF6A+TBCExcW8bt1OYRbUxsZm+fLlRUVFT2yj0Wj++c9/Ojk5EZGVldXy5Z//bpPaKigoXLp0qbBu0M7ObvXq1U9c4bh//35gYKCwWrVt27ZCRtazMWa1Oi8zMzAhwUKppMREh9u3Q4qLj6WkDMvMfD8lZVhx8bHU1PFCQF654lVScrLeDQlycnLkcrlQfKtWrcLCwvSBl5+fv2bNGhcXF+GE1blz502bNpWUlNZyz3l53Ls3E/GQIVxS0sAyn6pSvIBVq9VhYWEtWrQQzs76rpXwevPmzVu2bLlgwYLbt2+L1eIT7ty5M3PmTCEzbG1ticjMzGzWrFmiX/MuLU1ITn5ZqaRr1ybr+1i3b68Rjqtr1yaWl9ctlmqk1WojIiKEwLOxsQkODq4++CsrK9u6dav+FhozM7MFC7b++CM/7M0+S2goE3GzZhwf3/AyAURg1CCsquKICHZ15a5dky0tLf39/W/cuCH86ObNmzNmzEhLS6u+vT6iRo067OrKERFcp8k5tZrDw9nD427LlvbCKenWrVtP2zgxMVE/Sevt7R1f979Rrbbi7t2wxER7pZISEixu3pSr1fq1mlq1Ole/ZWFh1MWLnZRKUipl6emzKivv1LUtZi4tLQ0JCRHO+88YzlZUVERERAjzWjY2Dl27agIDOTv7uTvnoUOZiPv25fz8elRXK9evX3d1dQ0JCWn49eBjx3jUqMzmze2I6K233tL/j/7tt9/69eunP1kTUYsWLT7//PPi4uIGl/+IkBnCHXXNmjULDg6+deuWQqEQlpZYWVnJ5fK6LtytkUZzPyNjkVJpplRSUlKXO3c26vtYKtW127dDhMNPqTS7fj2goiK93g1duHDhlVdeEf7dRo8enZyc/LQtlUqlXC63tbX18sol4nbtWKHgrKyn7vn779nMjGUy3r273tUBiMx4QRgbywMHMpFukHH6tK6bXFRUtGLFCmGAOG3atN+/UalMGD1aI7xx2DCu5WzTgQO6pY9EvGDBnsTExNq8KyoqSrjlTgjOWp68tFrt1as/JiV1EbrkaWl+5eVPPXEIqqpKsrM/SUiwFgaOd+5sruXUJVdb7yOcp/z9/a9du/a85qr2798/d2608A9iY8MLF3JKSs0bV1by668zEbu58cN18gaxcuVK4Vfo0qXLrl276jcQv3eP585lmYyJ+M03t/33v/8VXhd6UUL4devW7eeff05OTtZPjzs7O4eEhKhUqob/FomJiUOGDBF26+fnp+/bMXNmZqZ+vG5nZxccHNyQVcr370eeP99GqaSEBMvMzMCqKiHLH+tjqdV52dkK4bg6d65ZdrZCoymoUyulpaXBwcHCDIqrq2str9Pn5+dv2KB1d9f9xVlZ8R//yMdqmqZdupSJODS0TkUBGJYxgjA1lQMCdH8hHTo8GtgJl+jatm0rXIEPCAiofhKpTqvlyEh2c2MiNjPjWbN0NyHV6MoVnjBB11z37hwZWbdqS0pKgoODhbsFHB0dq0821kipVI4YMaJjR4fERKdLl3oVFkbXvi2VKi0tzV+ppI0bh3p4ePzyyy/PfUt8fLyPj49w2vXy8joqrMCrNaWSZ81ic3Pdv6S/P8fEPLnN4cNsbs4uLvz4+NwgYmJi9PdQ9+nT58CBA7V/r1bLERG6tTxWVqxQ6G72EMZn+ku/CoWi+ogzPj5+xAjds9Z69OgRGRlZ70VADx7wJ5+sFnLOzc3tP//5T42bXb58+Ynr07Xv9AhUqtTUVF+hm5WSMqKs7DlXzVWqtOvXpyqVMqWSzp9vFRf3dS0jv3pHUC6XP3jwoE51VlXxoUM8aZLuABPWHv/v/3J+PgsdALWai4u5Foc5gFGJH4TFxY8O+qwsVijY2pqJuHlzDg5mfYc4Job/+EfdZ9q98sorp06deu6eS0o4OFi3N0dHDgt78oJEXh4HBur+CB0dOSSE693jT01N1T9UesCAAcdq6txmZmbqrwm5uLicPLlTq63FFZLfKSw8MGzYYKGtadOmZT1lXikjI0O/pFa/3qcezTFzaiovWsQ2NrqzlY8Px8Y+dqrav7+2I++Gq6qqioyM7NKli/Av4Ovrm5CQUJs3rl6tq//VVx9ldlKSSj/F7evre1W4Ye1xWq32+++/79mzp7DZjBk36tidYGaOimI3N/bxibSwsAgMDPz9de4n/Prrr15eXkKLHh4ekbXroJWWli5fvvzMGc+Ha6y2M9c2tktLz169OiYuzsnevqWbm1t4eLjm6Vfwrl+/rj/gPT09z5w5U8tWapSTwyEhup5rhw786ac8ejQzc1ISL13akB0DGIT4QRgU9OigFwaCwhhOv0bh6lXd61ZWKh+fsXv27KlTlzw5mceN050BP/nk0ek7J4cdHZmILS05MFCcK1tRUVFdu3bVz0BmPFziJlyfs7OzE64ANXy5aWVlpX6hh62t7RMLEwoKChQKhTBItbW1VSgUolziunePg4PZ2ZmJeN48E5+qysrKQkNDHR0dicjRscvbb6vT05/zlrw87teP9VN3paUcHMxWVjxs2Ly2bdtGREQ8+7hSq9VfffXVqFErhGPJ359ruTw5NfXRETh0KCclPWdeWk+Y03Z3dxeOqCFDhtTYwdI7cOCAcH+tn1/PjIx36neDxLlz/9VfJR04cOChQ4ee2EA49oTnXDs4OFRfv91AFRW8Zw/v2MFBQTxuHEdEIAihkTJIEFY/6D//nC9c0P0oP//JAWK9F0lERfGAAfzhh4+dvt99l319a3s6q6WKioqwsDAh82xtbT/77LNdu3bpH93p7+9/XbjVQAw5OTn6MV+3bt0iIyPVanV4eHibNm30s8cZYq83LynhvXsby6kqPz//z3/+89ixZ4Q7rz/4gHNzn7W9Pul++IE7dmQiNjfnZcsKat8vKS3lkBC2t9f12AIC+CnT88zMlZUcEqIbSQtzEvUYk1dWVoaHh+tX8/r6+iYlJT2xjXAk6Cck6rF0qzphzK3v0o0dO1apVAo/iouL098gGBAQIMqKnt8LCuIjR3jYMI6LQxBCY2SQIKzxoFeruUsX3alq0SIW5S/uidO34Z5PkZGRMWXKFKpm8ODBJ0829OaHGsXGxvbt21doRf9o/9GjRxvupjRuZKeqrCyWy3VT3HZ2rFDwM+Yds7N5yhTd+Oyll7gWU+w1yMt71EWzsuLAQC4oeHKb2Fjd8iuZ7DlXqWujpKREv+jXwsJCLpcLK12r3wTSvHnzkJAQscZn5eXlX3zxhXBLkkwmmzx58uTJk4Wjq1evXnFxcaK0UqOgII6P56NHefBg0x9dAL9nkCB82kG/fj2PHcvnz4vZljFP37/++uu8efMCAwN37dpl0MesCANBZ2fniRMndu/evZbXkxqiEZ6qLl9+tMaqXTsOD3+yo1NZyWFhbGfHRGxvX8M147q6do2nTdMtQHVy4uPHH028f/edrpK+ffn48Qa1Ut2tW7fkcrmFhYUQe2+//XafPn2EcJo6dWpOTo5oLT1UVFQUHBwsLNK2srISbvYQZfXsMwhHFzMvWNBYji6A6gwVhFzTQS/6Yz0a4elbRIWFhcXFxXV6un+9NdpTVWwsv/yyLoT69+f79x+F0/z5utenT2cR75I/e5bHjuUuXXjFikcT7x9+yEOHcnAwG+L/RkpKSkBAgEwmE56NJ9zsIX4z1WRmZm7evHnv3r3pz70SCyABMmYmsW3ZQjExtHAhjRsn+r4f89lnNGECDRlCCxdSy5a0YYNhmwOTYKbvvqNPP6Vx46hVK/rtN4qNpUuX6O9/p7Nn6YsvyBCfjXjvHv3jH3TqFM2cSZ6etG0brV9PD5+yaRAnT54sLCzMysqaM2eOcCc+ABiHQT6Y98wZ2rePxo83xL4fo1CQpSUR0VdfUXm5wZsDk5DJKCCAJk+m8nJav54sLGj7dvL0JAcHunjRUI22aUNEtGIFBQXRmjW6MgxK/yQXADAyM0PsVKslIjIzyL4fs3EjxcUREV28SH/5i8GbAxOytKSWLYmIVqygLVsoP98YjdrY0Nq1pFAYoy0AMBWDhFVVFZFRgpCINBpSq0mtNkZb0BgYOZxGjKCHt+EBwIvJIFOjwojQ3NwQ+37SunW0Ywfl59PDzxaEF59xwgkT7wASYcAgNM6IcMUKGj+ezp2jiAhjNAemtWqV7ostWwze1saNNHgwjR9PFy9SRARt2mTwFgHAJAw4NWqcESGA4WDiHUAKmvaIMChI94WnJ3l6Grw5kBpMvANIAUaEAE+1YgXt3k0hIaauAwAMySAjwtatv/L2vm1j8wbRQEPsHwAAQCwGCcJbtw4plVFVVd4IQmi6MPEOIBEGuqFeS0Rmxlk2CgAA0AAGukZYRUTmuEgIAACNHkaEAAAgaQhCAACQNEyNAgCApGFECAAAkoYRIQAASJpBPqFeqVQWFhZ6e3s7ODiIvnMAAAARGSQIAQAAmgpcxgMAAElDEAIAgKQhCAEAQNIQhAAAIGkIQgAAkDQEIQAASBqCEAAAJA1BCAAAkoYgBAAASUMQAgCApCEIAQBA0hCEAAAgaQhCAACQNAQhAABIGoIQAAAkDUEIAACShiAEAABJQxACAICkIQgBAEDSEIQAACBpCEIAAJA0BCEAAEgaghAAACQNQQgAAJKGIAQAAElDEAIAgKQhCAEAQNIQhAAAIGkIQgAAkDQEIQAASBqCEAAAJA1BCAAAkoYgBAAASUMQAgCApCEIAQBA0hCEAAAgaQhCAACQNAQhAABIGoIQAAAkDUEIAACShiAEAABJQxACAICkIQgBAEDSEIQAACBpCMLa+umnn9q2bWvqKqApUalU//znP0tKSkxdCAA8i4yZTV1D05CXl3f58uWRI0eauhBoMnJzc9u0aZORkdGxY0dT1wIAT4UR4VPt3LnT3d3d1dV14cKF5eXlqampoaGhRLRy5coePXq4ubl9+umn6EbAM0yaNImIXnvttZ9//tnf33/VqlW9e/cmor179/r4+Li6uk6ZMiUjI4OIfvvtt0mTJi1fvtzNza1///7Hjh0jHGkARsNQk3PnzslkspCQkP3793ft2nXLli3R0dGtW7c+duyYs7NzXFzcN998Y2dnd/z4cVNXCo1XVFQUEX377bcHDx4kInd3982bN6tUqtatW69bt+7EiRPDhw+Xy+XMHB0dTUSvv/56dHS0t7f3kCFDcKQBGI2FiXO4sYqIiBg0aJBCoSAiCwuLmzdvCq+np6eXlZWVlZXNnDmzf//+Li4upqwSGrchQ4YQ0dChQ1NTU4nowIEDHh4eubm5YWFh06dPv3PnTo8ePc6fPy9sbGlpuWfPnpYtW2o0mnnz5uFIAzAaTI3W7ObNm3369BG+9vPze++994SvZ8yY4e/vP3HixNatW4eHhzs6OpquRmhKHB0dPTw8iKhFixYJCQmOjo6dOnX6+eef9Rs4OTm1bNmSiOzt7VUqFY40AKNBENbMw8MjJydH+PrEiRO7d+8Wvr558+bGjRvv3Lnz5ZdfRkVFbd261XQ1QlMik8mEL3744Yfw8PCDBw+WlZUtXbpUv4G5uXn17XGkARgNgrBmEydOPHHiRGxs7L179/70pz8VFhYKr+/bt2/mzJlENGPGjM6dO5eXl5uySmjchGy7f/9+9Rdv3LjRoUMHb2/ve/fuRUREVFVV1fheHGkARoMgrNmQIUP8/PzGjh3bsWNHV1fXuXPnEpFMJps9e3ZaWlr79u2dnZ2Zed68eaauFBovJyenQYMGvfTSS7m5ufoR4ezZsysqKlq1atWzZ88xY8ZcuXJl06ZNT7wRRxqAMeE+wme5d++ehYWFk5NT9RdVKlVqamqrVq3at29vqsKgCcnPz2/VqlX1V5j56tWr3bp1s7S0LCgosLa2trW1/f0bcaQBGAeCEAAAJA1TowAAIGkIQgAAkDQEIQAASBqCEAAAJA1BCAAAkvb/QK83n1NJPB4AAAAASUVORK5CYII=\n", | |
"text/plain": [ | |
"<PIL.PngImagePlugin.PngImageFile image mode=RGB size=600x200 at 0x12023FA90>" | |
] | |
}, | |
"execution_count": 3, | |
"metadata": {}, | |
"output_type": "execute_result" | |
} | |
], | |
"source": [ | |
"Draw.MolsToGridImage([m1_cis,m1_trans],legends=[\"cis\",\"trans\"])" | |
] | |
}, | |
{ | |
"cell_type": "markdown", | |
"metadata": {}, | |
"source": [ | |
"## Compare InChIKeys" | |
] | |
}, | |
{ | |
"cell_type": "code", | |
"execution_count": 4, | |
"metadata": {}, | |
"outputs": [ | |
{ | |
"data": { | |
"text/plain": [ | |
"'AQIXAKUUQRKLND-UHFFFAOYSA-N'" | |
] | |
}, | |
"execution_count": 4, | |
"metadata": {}, | |
"output_type": "execute_result" | |
} | |
], | |
"source": [ | |
"inchikey_cis = Chem.inchi.MolToInchiKey(m1_cis)\n", | |
"inchikey_cis" | |
] | |
}, | |
{ | |
"cell_type": "code", | |
"execution_count": 5, | |
"metadata": {}, | |
"outputs": [ | |
{ | |
"data": { | |
"text/plain": [ | |
"'AQIXAKUUQRKLND-UHFFFAOYSA-N'" | |
] | |
}, | |
"execution_count": 5, | |
"metadata": {}, | |
"output_type": "execute_result" | |
} | |
], | |
"source": [ | |
"inchikey_trans = Chem.inchi.MolToInchiKey(m1_trans)\n", | |
"inchikey_trans" | |
] | |
}, | |
{ | |
"cell_type": "code", | |
"execution_count": 6, | |
"metadata": {}, | |
"outputs": [ | |
{ | |
"data": { | |
"text/plain": [ | |
"True" | |
] | |
}, | |
"execution_count": 6, | |
"metadata": {}, | |
"output_type": "execute_result" | |
} | |
], | |
"source": [ | |
"inchikey_cis == inchikey_trans\n", | |
"#True" | |
] | |
}, | |
{ | |
"cell_type": "markdown", | |
"metadata": {}, | |
"source": [ | |
"## Compare InChIs" | |
] | |
}, | |
{ | |
"cell_type": "code", | |
"execution_count": 7, | |
"metadata": {}, | |
"outputs": [ | |
{ | |
"data": { | |
"text/plain": [ | |
"'InChI=1S/C10H16N6S/c1-8-9(16-7-15-8)5-17-4-3-13-10(12-2)14-6-11/h7H,3-5H2,1-2H3,(H,15,16)(H2,12,13,14)'" | |
] | |
}, | |
"execution_count": 7, | |
"metadata": {}, | |
"output_type": "execute_result" | |
} | |
], | |
"source": [ | |
"inchi_cis = Chem.inchi.MolToInchi(m1_cis)\n", | |
"inchi_cis" | |
] | |
}, | |
{ | |
"cell_type": "code", | |
"execution_count": 8, | |
"metadata": {}, | |
"outputs": [ | |
{ | |
"data": { | |
"text/plain": [ | |
"'InChI=1S/C10H16N6S/c1-8-9(16-7-15-8)5-17-4-3-13-10(12-2)14-6-11/h7H,3-5H2,1-2H3,(H,15,16)(H2,12,13,14)'" | |
] | |
}, | |
"execution_count": 8, | |
"metadata": {}, | |
"output_type": "execute_result" | |
} | |
], | |
"source": [ | |
"inchi_trans = Chem.inchi.MolToInchi(m1_trans)\n", | |
"inchi_trans" | |
] | |
}, | |
{ | |
"cell_type": "code", | |
"execution_count": 9, | |
"metadata": {}, | |
"outputs": [ | |
{ | |
"data": { | |
"text/plain": [ | |
"True" | |
] | |
}, | |
"execution_count": 9, | |
"metadata": {}, | |
"output_type": "execute_result" | |
} | |
], | |
"source": [ | |
"inchi_cis == inchi_trans\n", | |
"#True" | |
] | |
}, | |
{ | |
"cell_type": "markdown", | |
"metadata": {}, | |
"source": [ | |
"## Compare Auxilliary Info" | |
] | |
}, | |
{ | |
"cell_type": "code", | |
"execution_count": 10, | |
"metadata": {}, | |
"outputs": [ | |
{ | |
"data": { | |
"text/plain": [ | |
"('InChI=1S/C10H16N6S/c1-8-9(16-7-15-8)5-17-4-3-13-10(12-2)14-6-11/h7H,3-5H2,1-2H3,(H,15,16)(H2,12,13,14)',\n", | |
" 'AuxInfo=1/1/N:17,1,8,9,11,5,14,16,12,3,6,2,7,4,15,13,10/rA:17CNCNCNNCCSCCNCNCC/rB:s1;d+2;s3;s4;t5;s3;s7;s8;s9;s10;s11;s12;d13;s14;d12s15;s16;/rC:;;;;;;;;;;;;;;;;;')" | |
] | |
}, | |
"execution_count": 10, | |
"metadata": {}, | |
"output_type": "execute_result" | |
} | |
], | |
"source": [ | |
"inchi_aux_cis = Chem.inchi.MolToInchiAndAuxInfo(m1_cis)\n", | |
"inchi_aux_cis" | |
] | |
}, | |
{ | |
"cell_type": "code", | |
"execution_count": 11, | |
"metadata": {}, | |
"outputs": [ | |
{ | |
"data": { | |
"text/plain": [ | |
"('InChI=1S/C10H16N6S/c1-8-9(16-7-15-8)5-17-4-3-13-10(12-2)14-6-11/h7H,3-5H2,1-2H3,(H,15,16)(H2,12,13,14)',\n", | |
" 'AuxInfo=1/1/N:17,1,8,9,11,5,14,16,12,3,6,2,7,4,15,13,10/rA:17CNCNCNNCCSCCNCNCC/rB:s1;d-2;s3;s4;t5;s3;s7;s8;s9;s10;s11;s12;d13;s14;d12s15;s16;/rC:;;;;;;;;;;;;;;;;;')" | |
] | |
}, | |
"execution_count": 11, | |
"metadata": {}, | |
"output_type": "execute_result" | |
} | |
], | |
"source": [ | |
"inchi_aux_trans = Chem.inchi.MolToInchiAndAuxInfo(m1_trans)\n", | |
"inchi_aux_trans" | |
] | |
}, | |
{ | |
"cell_type": "code", | |
"execution_count": 12, | |
"metadata": {}, | |
"outputs": [ | |
{ | |
"data": { | |
"text/plain": [ | |
"False" | |
] | |
}, | |
"execution_count": 12, | |
"metadata": {}, | |
"output_type": "execute_result" | |
} | |
], | |
"source": [ | |
"inchi_aux_cis == inchi_aux_trans" | |
] | |
}, | |
{ | |
"cell_type": "markdown", | |
"metadata": {}, | |
"source": [ | |
"It seems the isomers differ in their aux info : d+2 vs. d-2. \n", | |
"\n", | |
"Their InChIs and InChIKeys as generated by RDKit are otherwise identical. Not sure why this discrepancy arises." | |
] | |
}, | |
{ | |
"cell_type": "markdown", | |
"metadata": {}, | |
"source": [ | |
"### Practical Application\n", | |
"\n", | |
"Originally, Gustavo wanted to generate all stereoisomers, then filter out duplicates on InChIKey.\n", | |
"\n", | |
"Why not filter out duplicates using canonical SMILES? " | |
] | |
}, | |
{ | |
"cell_type": "code", | |
"execution_count": 13, | |
"metadata": {}, | |
"outputs": [ | |
{ | |
"data": { | |
"text/plain": [ | |
"'C/N=C(/NC#N)NCCSCc1nc[nH]c1C'" | |
] | |
}, | |
"execution_count": 13, | |
"metadata": {}, | |
"output_type": "execute_result" | |
} | |
], | |
"source": [ | |
"m1_cis_can = Chem.CanonSmiles('C/N=C(/NC#N)NCCSCc1nc[nH]c1C')\n", | |
"m1_cis_can" | |
] | |
}, | |
{ | |
"cell_type": "code", | |
"execution_count": 14, | |
"metadata": {}, | |
"outputs": [ | |
{ | |
"data": { | |
"text/plain": [ | |
"'C/N=C(\\\\NC#N)NCCSCc1nc[nH]c1C'" | |
] | |
}, | |
"execution_count": 14, | |
"metadata": {}, | |
"output_type": "execute_result" | |
} | |
], | |
"source": [ | |
"m2_trans_can = Chem.CanonSmiles('C/N=C(\\\\NC#N)NCCSCc1nc[nH]c1C')\n", | |
"m2_trans_can" | |
] | |
}, | |
{ | |
"cell_type": "code", | |
"execution_count": 15, | |
"metadata": {}, | |
"outputs": [ | |
{ | |
"data": { | |
"text/plain": [ | |
"False" | |
] | |
}, | |
"execution_count": 15, | |
"metadata": {}, | |
"output_type": "execute_result" | |
} | |
], | |
"source": [ | |
"m1_cis_can == m2_trans_can" | |
] | |
}, | |
{ | |
"cell_type": "markdown", | |
"metadata": {}, | |
"source": [ | |
"Above, canonical SMILES (within RDKit) are able to distinguish the stereoisomers. \n", | |
"\n", | |
"However, note that different toolkits have different forms of canonical SMILES, so you'd need to stay within RDKit for your deduplication task!\n", | |
"\n", | |
"See Greg's post on canonical SMILES:\n", | |
"https://github.com/rdkit/rdkit/issues/2747\n", | |
"\n", | |
"See Paolo's gist on comparing molecule identity:\n", | |
"https://gist.github.com/ptosco/1e1c23ad24c90444993fa1db21ccb48b" | |
] | |
}, | |
{ | |
"cell_type": "code", | |
"execution_count": null, | |
"metadata": {}, | |
"outputs": [], | |
"source": [] | |
}, | |
{ | |
"cell_type": "markdown", | |
"metadata": {}, | |
"source": [ | |
"## Possible related issue\n", | |
"\n", | |
"Example from \n", | |
"\n", | |
"https://sourceforge.net/p/inchi/mailman/inchi-discuss/thread/CA%2BZ3Zfc0nLLRJPYwQ2EJL-ePPH8RgzedbFxkRMKi6ifcsp4FHw%40mail.gmail.com/#msg31493868" | |
] | |
}, | |
{ | |
"cell_type": "code", | |
"execution_count": 16, | |
"metadata": {}, | |
"outputs": [], | |
"source": [ | |
"m2_cis = Chem.MolFromSmiles(\"C1CCCCCC(=O)OCC/C=C\\C1\") \n", | |
"m2_trans = Chem.MolFromSmiles(\"C1CCCCCC(=O)OCC/C=C/C1\") " | |
] | |
}, | |
{ | |
"cell_type": "code", | |
"execution_count": 17, | |
"metadata": {}, | |
"outputs": [ | |
{ | |
"data": { | |
"text/plain": [ | |
"'FZKPUQQWULXMCD-ALCCZGGFSA-N'" | |
] | |
}, | |
"execution_count": 17, | |
"metadata": {}, | |
"output_type": "execute_result" | |
} | |
], | |
"source": [ | |
"inchikey_cis = Chem.inchi.MolToInchiKey(m2_cis)\n", | |
"inchikey_cis" | |
] | |
}, | |
{ | |
"cell_type": "code", | |
"execution_count": 18, | |
"metadata": {}, | |
"outputs": [ | |
{ | |
"data": { | |
"text/plain": [ | |
"'FZKPUQQWULXMCD-FNORWQNLSA-N'" | |
] | |
}, | |
"execution_count": 18, | |
"metadata": {}, | |
"output_type": "execute_result" | |
} | |
], | |
"source": [ | |
"inchikey_trans = Chem.inchi.MolToInchiKey(m2_trans)\n", | |
"inchikey_trans" | |
] | |
}, | |
{ | |
"cell_type": "code", | |
"execution_count": 19, | |
"metadata": {}, | |
"outputs": [ | |
{ | |
"data": { | |
"text/plain": [ | |
"False" | |
] | |
}, | |
"execution_count": 19, | |
"metadata": {}, | |
"output_type": "execute_result" | |
} | |
], | |
"source": [ | |
"inchikey_cis == inchikey_trans" | |
] | |
}, | |
{ | |
"cell_type": "markdown", | |
"metadata": {}, | |
"source": [ | |
"InChiKey function of PubChem seems to be working OK for cis/trans isomers in this case.\n", | |
"\n", | |
"Check their aux info:" | |
] | |
}, | |
{ | |
"cell_type": "code", | |
"execution_count": 20, | |
"metadata": {}, | |
"outputs": [ | |
{ | |
"data": { | |
"text/plain": [ | |
"('InChI=1S/C12H20O2/c13-12-10-8-6-4-2-1-3-5-7-9-11-14-12/h5,7H,1-4,6,8-11H2/b7-5-',\n", | |
" 'AuxInfo=1/0/N:1,2,14,3,13,4,12,5,11,6,10,7,8,9/rA:14CCCCCCCOOCCCCC/rB:s1;s2;s3;s4;s5;s6;d7;s7;s9;s10;s11;d-12;s1s13;/rC:;;;;;;;;;;;;;;')" | |
] | |
}, | |
"execution_count": 20, | |
"metadata": {}, | |
"output_type": "execute_result" | |
} | |
], | |
"source": [ | |
"inchi_aux_cis = Chem.inchi.MolToInchiAndAuxInfo(m2_cis)\n", | |
"inchi_aux_cis" | |
] | |
}, | |
{ | |
"cell_type": "code", | |
"execution_count": 21, | |
"metadata": {}, | |
"outputs": [ | |
{ | |
"data": { | |
"text/plain": [ | |
"('InChI=1S/C12H20O2/c13-12-10-8-6-4-2-1-3-5-7-9-11-14-12/h5,7H,1-4,6,8-11H2/b7-5+',\n", | |
" 'AuxInfo=1/0/N:1,2,14,3,13,4,12,5,11,6,10,7,8,9/rA:14CCCCCCCOOCCCCC/rB:s1;s2;s3;s4;s5;s6;d7;s7;s9;s10;s11;d+12;s1s13;/rC:;;;;;;;;;;;;;;')" | |
] | |
}, | |
"execution_count": 21, | |
"metadata": {}, | |
"output_type": "execute_result" | |
} | |
], | |
"source": [ | |
"inchi_aux_cis = Chem.inchi.MolToInchiAndAuxInfo(m2_trans)\n", | |
"inchi_aux_cis" | |
] | |
}, | |
{ | |
"cell_type": "code", | |
"execution_count": 22, | |
"metadata": {}, | |
"outputs": [ | |
{ | |
"data": { | |
"text/plain": [ | |
"False" | |
] | |
}, | |
"execution_count": 22, | |
"metadata": {}, | |
"output_type": "execute_result" | |
} | |
], | |
"source": [ | |
"inchi_aux_cis == inchi_aux_trans" | |
] | |
}, | |
{ | |
"cell_type": "markdown", | |
"metadata": {}, | |
"source": [ | |
"As before, the aux info differs: d-12 vs. d+12." | |
] | |
}, | |
{ | |
"cell_type": "code", | |
"execution_count": 23, | |
"metadata": {}, | |
"outputs": [ | |
{ | |
"data": { | |
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAlgAAADICAIAAAC7/QjhAAAegklEQVR4nO3de1hUdf4H8M/MwHARGUBEBw0QxAQ1r6viPcW8S1lupuB98XF/q+XzbI3tr8TULbZsZc1dFx9dU9Mt7Fk3HzOLVBATpcwLSqEyMiOIoDhg6HCZmfP74/wiUy5zPd85c96vv3hmZ+a87Xnv+Zwz5ybjOI4AAACkSs46AAAAAEsYhAAAIGkYhAAAIGkYhAAAIGkYhAAAIGkYhAAAIGkYhAAAIGkYhAAAIGkYhAAAIGkYhAAAIGkYhAAAIGkYhAAAIGkYhAAAIGkYhAAAIGkYhAAAIGkYhAAAIGkYhAAAIGkYhAAAIGkYhAAAIGkYhAAAIGkYhAAAIGkYhAAAIGkYhAAAIGkYhAAAIGkYhAAAIGkYhAAAIGlerAMwYDAYdD8rLS3V6/U6nW7KlCndunWbM2eOSqViHRBECb0CECkZx3GsM7iKwWDQarVarfbmzZsVFRX83yUlJTU1NY+/2dfXt76+3sfHZ+bMmSkpKVOnTlUoFIJHBhFArwA8jOgHoclkKisre2QzXKfT6fX6hoaGFj8SGBgYERERFRUV+TO1Wl1SUrJ3795jx45ZLBYiioyMXLBgwYIFC6Kjo4X9B4FbQK8ApEM0g7ChoaG8vPyRzXCtVqvX600mU4sfCQ4OVqvV4eHh0Q/hX2ltKeXl5R999NH27duvXbvGvzJ48ODU1NS5c+cGBAS45B8GTKFXAODWg7C0tDQnJ2fr1q06na6ysrLF9ygUCrVa3bwZHhERwf8RFRXl5+dn96LPnj27bdu2ffv21dXVEZGfn9/06dNTU1MnTJggk8ns/lpwB+gVADzMfQfh/v37ly1bFhwcrNVqiUipVHbv3v2RLXG1Wt2jRw9/f38XZTAajYcOHdq2bdvRo0f5/1C9evV66aWXFi5cGBUV5aKFgkuhVwDwKM79VFZWJiUl8fGSkpKys7PLy8stFgvDSMXFxWlpac0rKblcnpiYuGvXrvv37zNMBTZBrwCgRW43CA8fPqxWq4lIpVJlZmayjvMrZrM5Ozs7JSWleV8hKCgoNTU1Ly+PdTRoB3oFAK1xo0FYW1ubmprKrwgmTJig1+tZJ2qVwWDIzMwcOXJk8451fHx8enp6ZWUl62jwKPQKANrmLoPwm2++6dmzJxH5+vqmp6ebzWbWiaxy+fJljUYTFhbGr7YUCkViYmJWVlZjYyPraMBx6BUAWIH9IKyvr9doNPxVxk899dSFCxdYJ7KZyWTKzs6ePXu2t7c3v+bq2rXrypUrz58/zzqadKFXAGAlxoOwsLBwwIABROTl5aXRaBoaGtjmcVB1dXVmZib/L+INHjw4IyPjzp07rKNJC3oFANZjNgjNZnNGRoaPjw8RRUdHnzx5klUSVzh9+vSyZcuCgoL41Zafn9/rr79eU1PDOpfnQ68AwFZsBqFWqx0zZgwRyWSy1NTUurq6x99z9+7dpUuX3rhxQ/h4zlJfX5+VlTV9+nQvL6+EhIS4uLjq6mrWoTwZegUAdmAwCHft2sXfVqpr166HDh1q8T1fffVV9+7diejZZ58VOJ4rFBUVxcfHE9Ho0aONRiPrOJ4JvWIdB0CsBB2Et27dmjFjBv+rzuzZs1vcjH3w4IFGo5HL5UQ0fPjw4uJiIRO6Tnl5eUREBBHNnDnTZDKxjuNR0Cv0CsARwg3C/fv3h4aG8hcL79mzp8X3FBQUPPnkk0Tk7e2dlpbmYf/Hvnz5cnBwMBEtX76cdRbPgV6hVwAOEmIQ1tTUNF/R/Mwzz5SVlT3+nqampvT0dP408T59+pw9e1aAYMLLzc3lz+PYuHEj6yyih141Q68AHOHyQZidnf3EE0/wZ7hlZGS0eGvHoqKiIUOG8Oc4rFy5sr6+3tWpGPrkk0/kcrlMJtu9ezfrLCKGXj0CvQKwmwsHodFobD4qM2zYsBaPylgslszMzA4dOhBRZGTk8ePHXZfHfbz33ntEpFQqs7OzWWcRH/SqNegVgH1cNQgLCgp69+7NX9Hc2lEZnU43fvz45nMcDAaDi8K4oVdeeYWIAgMDcZcQm6BXbUOvAOzg/EHIH5VRKpX8LYNbOyqTlZUVEhJCRGFhYQcOHHB6DDdnNptfeOEFIurWrZtOp2MdRwTQK2ugVwB2cPIg/OGHH5qPyqSmprb4WDWDwTBv3jx+g33KlCk3b950bgaxMBqNo0aN4k/iuHv3Lus4bg29sh56BWArpw1CK4/KHDlyJDw8nP/1xt0eCye86urquLi40JCQ88nJnMjvh+ki6JUd0CsAmzhtEO7evZs/gXvx4sW1tbWPv+H+/fsrV66UyWRENGLEiGvXrjlr0aJ2/fr14qFDOSJu7lyO6dPS3RN6ZR/0CsB6zhmEBQUFMpksMjKytaMy+fn5sbGxJLbHwgnk4kUuKIgj4l57jXUU94JeOQS9ArCOcwbhxIkTiWj16tWP/0+NjY3Nj4UbNGjQpUuXnLJET3PsGKdUckRcRgbrKG4EvXIUegVgBScMwtzcXCJSqVQt3uPRbDaPGzdOoVB4wGPhXGvvXk4m4+Ry7tNPWUdxC+iVc6BXAO1xwiBMSEggog0bNrT2huvXr+fn5zu+IM+3YQNHxPn6cnl5rKOwh145DXoF0CYZx3HkgIMHDyYlJXXu3LmkpKRjx46OfBUQEa1YQVu2UKdOdPIk9e7NOg0z6JWToVcArZM78mGLxZKWlkZEb7zxBtZWzpGRQc89R9XVNGUK3brFOg0b6JXzoVcArXNoj3Dv3r3JyckRERFXrlzhz3EHJzAaKTGRTp2ifv0oL49UKtaBhIZeuYTkewXQGvv3CE0m0/r164koLS0Naytn8vOjzz6jJ5+kwkKaNYsaG1kHEhR65SrS7hVAG+wfhDt27CguLu7Vq9f8+fOdGAiIiEJD6fPPKSyMjh2j995jnUZQ6JULSbhXAG2w86fR+vr62NjYsrKyrKys2bNnOz0WEBEVFND779O//kUdOrCOIhD0SgjS6xVA2+wchBs3bnz11Vefeuqpc+fO8U+GA3AcegUAwrNnENbV1cXExFRVVX3xxReTJ092RSx4VGEhHTlCV69Sz56UmEiDBrEO5HzoFQMS6BVAu+zZ6N64cWNVVdWoUaOwthLIvn00ZAgdOUKBgXTsGA0bRh98wDqT86FXQpNGrwDaZ+sV+Ldv3w4MDCSinJwcJ1/cDy26eZPz8eH+/OdfXlm3jlMouCtX2GVyPvRKaNLoFYA1bN4jfOedd+7duzdt2rSxY8c6fyzD4z75hIKD6bXXfnll6VLiODp4kF0m50OvhCaNXgFYw7ZBePPmza1bt8pkMv5KLxBCSQn17k1eXr+8olbTwIF04wa7TE6GXjEggV4BWMm2Qbh27Vqj0fjb3/524MCBLgoEjwoIoOvXf/WKxULXr1NoKKNAzodeMSCBXgFYyYZBePXq1Q8//FChUKxdu9ZleeAxffuSTkeFhb+88v33dPcuDRnCLpMzoVdseHqvAKxnwyBMS0trampatGhRb9y9XkizZ1OPHvTyy6TXExHp9bRiBQ0dSp5yaiV6xYan9wrAetZeR1hYWDhgwAClUllcXBwREeHqWPArJSWUnEynT5OfH9XXU1wcffMNBQWxjuUE6BVLntsrAJtYu0f4pz/9yWKxLF++HGsrBmJiKD+fqqroyy9JLqerV8nbm3Um50CvWPLcXgHYxKo9woKCguHDh3fo0OHatWtdunQRIBa0asgQOnuWvv6aJkxgHcVR6JUb8aBeAdjKqj3C1atXcxy3atUqrK3Y4y+zy81lncMJ0Cs34kG9ArBV+4Pwyy+/PH78eHBw8KpVqwQIBO3gV1g5OYxjOAy9ci+e0isAO7QzCDmOS0tLI6LVq1cHBwcLEgnaNHYsKRRUUEBGI+so9kOv3I5H9ArAPu0Mwv/85z9nzpxRq9V/+MMfhAkE7VCpqF8/amigM2dYR7EfeuV2PKJXAPZpaxCazeY1a9YQ0Ztvvunv7y9UJGiPyA/noFduSuS9ArBbW4Nwz549RUVFUVFRS5YsESwQtE/kKyz0yk2JvFcAdmv18onGxsa4uDitVrtnz57k5GSBY0Fb7t6lzp1JqSSDgXx9WaexDXrlvsTcKwBHtLpHmJmZqdVq+/TpM3fuXCEDQftCQig+nurr6bvvWEexGXrlvsTcKwBHtDoIb926FRYWlpaWJpfb8xR7cC3RnuyOXrk10fYKwBGtroy2b99eVVU1evRoIdOAtUR7OAe9cmui7RWAI9rZKpfJZMLkANuMG0cyGZ06RY2NrKPYA71yUyLvFYB98POUOHXubOnfv6ZXr7Jz51hHAQ+CXoEkYRCK1YqEhODz5z86fpx1EPAo6BVIEAahWI0eM4aIcnE4B5wKvQIJwiAUq7FjxxLRyZMnTSYT6yzgOdArkCAMQrFSq9WxsbF1dXXncDgHnAe9AgnCIBQxfuMdv2KBc6FXIDUYhCKGFRa4AnoFUoNBKGLjxo0jory8PLPZzDoLeA70CqQGg1DEunfvHh0dXVtbe+HCBdZZwHOgVyA1GITihl+xwBXQK5AUL9YBwCFjx47duXPnl19+2Wj1PbHCfX1T6uutXYCfH61caWc4EC30CiSl1ecRdunSpaqq6tatW126dBE4E1ivtLS0R48egYGB9+7ds/Ijk6KijpSWWruAkBCqriYiKiykI0fo6lXq2ZMSE2nQIHviolciIbpeATgCe4TiFhUVFRkZqdPpFi1a1LlzZ2s+EuHvTw8eWLuADh2IiPbto0WLaNQoGjiQjh2j//1f+utfacUKe1ODu0OvQFq4VoSFhRHRrVu3WnsDuImUlBQi2rx5s6sWcPMm5+PD/fnPv7yybh2nUHBXrtjxZeiVWIirVwCOwMkyopeQkCCXy7VarasW8MknFBxMr732yytLlxLH0cGDrloiuAH0CqQDgxDaU1JCvXuT10O/oqvVNHAg3bjBLhOIH3oFbgODUPTy8/MtFkt0dLSrFhAQQNev/+oVi4WuX6fQUFctEdwAegXSgZNlRO/EiRNEdO7cOY1GY837I/z9/8emkxr69iWdjgoLqV+//3/x++/p7l0aMsSeuCAS6BVIBy6fEDchTnOvqKDevSkqij78kCIiSK+nF18ki4XOnLEjMHolCqLrFYAjsEcobvy9PxISEp5++mkrPxLu60s2XfisVFJ2NiUnU2Qk+flRQwNNnUp79tgXGEQBvQJJwSAUN36FNWnSpFWrVrlwMTExlJ9Pt2/TjRsUG0sdO7pwWeAG0CuQFAxCceNXWPydIV2uc2ey7tpqEDv0CiQFZ42KWFlZmVarValU/fv3Z50FPAd6BVKDQShiOTk5RDRq1CiFQsE6C3gO9AqkBoNQxAT9/QokA70CqcEgFDF+hcU/TxzAWdArkBoMQrGqqKi4evVqQEDAwIEDWWcBz4FegQRhEIoVv9k+atQoLy+c+gtOg16BBGEQihUO5IAroFcgQRiEYoUVFrgCegUShEEoSrdv3/7xxx/9/f0HDx7MOgt4DvQKpKmdQdjaLbmBrdzcXI7jRowYoVQqWWexB3rlnsTeKwD7tDoIly5dGhYWlpeXJ2QasJJ4f79Cr9yZeHsF4IhWB2HXrl2rqqreeusti8UiZCCwhnhXWOiVOxNvrwAc0erzCBsbG+Pi4rRa7Z49e5KTkwWOBW24e/du586dlUqlwWDw9fVlHcc26JXbEnWvABzR6h6hUql88803iejNN99sbGwUMBK0Izc312KxJCQkiHFthV65LVH3CsARbZ0sk5KSEh8fX1paumPHDsECQbvE/vsVeuWexN4rALu1NQgVCsW6deuIaP369Q8ePBAqErRD7Css9Mo9ib1XAHZr5/KJWbNmDRs2rKKiYsuWLcIEgrbV1tYWFhb6+PgMGzaMdRb7oVfuxjN6BWCfdgahTCZ76623iCg9Pd1gMAgSCdpy4sQJs9k8dOhQPz8/1lnsh165G8/oFYB92r+zzKRJk55++mmDwbBp0yYBAkHbPOb3K/TKrXhMrwDsYNUt1tLT02Uy2aZNmyorK10dCNrmSSss9Mp9eFKvAGxl1SAcOnTotGnT6urq/vKXv7g6ELThp59+On/+vLe3d0JCAussToBeuQkP6xWAray96fbbb78tl8u3bt2q1+tdGgjakJeXZzKZfvOb33To0IF1FudAr9yB5/UKwCbWDsJ+/fq9+OKL9fX169evd2kgaIPn/X6FXrkDz+sVgE1seAzTW2+95e3tvXPnzh9//NF1gaANHrnCQq+Y88heAVjPhkEYGxu7cOFCs9m8du1al+WBVtXV1X3//fdeXl4jRoxgncWZ0Cu2PLVXANaz7cG8a9eu9fPzy8rKOnfunIsCQWtOnTrV1NQ0aNCgjh07ss7iZOgVQx7cKwAr2TYIw8PDly9fznEcf99kEExTU9OGDRtGjhw5a9Ys1lmcD71ixbN7BWClVh/D1Jo7d+7ExMTcu3cvJycHBxWEwXHckiVLdu7c2aVLl6tXr3rkljt6JTwp9ArAGrbtERJRaGjoqlWriOiNN95wQR5owZo1a3bu3Onn53fgwAFPXVuhV8KTQq8ArGHzHiER1dXVxcTEVFVVffHFF5MnT3ZFLGi2ffv23/3udwqF4tNPP3322WdZx3Eh9EpI0ukVQLts3iMkooCAgFdffZWINBqNxWJxdiT4xeHDh5cvX05Ef/vb3zx+bYVeCUZSvQJoH2cXo9HYvXt3IsrKyrLvG6Bd3377bUBAABGtWbOGdRaBoFcCkGCvANpm5yDkOO6f//wnEfXq1aupqcmJgYBXUlLSpUsXIpo3b57FYmEdRzjolUtJtlcAbbDnGCHPZDL17du3uLh4x44dixcvds7+KRAR0Z07d0aOHHnlypXx48d/8cUXSqWSdSLhoFeuI+VeAbTB/kFIRHv37k1OTo6IiLhy5YqPj48TY0mZ0WicMGFCfn7+oEGDcnNz+V+xJAW9cgX0CqA19pws0+yll14aMGCAXq/PzMx0ViCJM5vN8+bNy8/P79Gjx+effy7NtRV65XToFUAbHNojJKKDBw8mJSV17ty5pKQElyI5bsWKFVu2bOnUqdPJkyd79+7NOg4z6JVzoVcAbXBoj5CIZs6cmZCQcPv27c2bNzslkJRt2LBhy5Ytvr6+n332mcTXVuiVE6FXAG1zdI+QiE6cODF27FiVSqXVakNCQpwSS4L27duXnJwsk8mysrKef/551nHYQ6+cAr0CaJeje4RENGbMmIkTJ9bW1r733nuOf5s0HT9+fNGiRRzHbdq0CWsrHnrlOPQKwCpOuQijoKBAJpNFRkYeOHDAKV8oKRcvXlSpVESk0WhYZ3Ev6JUj0CsAKzlnEHIct3v3bv5M98WLF9fW1jrraz3ejRs3+HupzJkzx2w2s47jdtAr+6BXANZz2iC0WCyZmZkdOnQgosjIyOPHjzvrmz1YdXU1f/LCuHHj6uvrWcdxR+iVHdArAJs4bRDyfvjhhyFDhhCRTCZLTU29f/++c7/fkxiNxlGjRhFRnz59DAYD6zhuDb2yHnoFYCsnD0KO45qamtLT0/m7N8XHx589e9bpi/AAZrP5hRdeIKJu3brp9XrWcUQAvbIGegVgB+cPQl5BQQH/44yXl1daWprJZHLRgkTq5ZdfJqLAwMALFy6wziIm6FXb0CsAO7hqEHIcZzQaNRqNXC4nomHDhhUXF7tuWeLy7rvvEpFSqfz6669ZZxEf9Ko16BWAfVw4CHnZ2dlPPPEEEfn5+WVkZODJLx9//LFcLpfJZLt372adRcTQq0egVwB2c/kg5DiupqYmNTWVv2zxmWeeKSsrE2Ch7iknJ4e/GOD9999nnUX00Ktm6BWAI4QYhLz9+/eHhoYSUVBQ0J49ewRbrvu4dOlScHAwES1fvpx1Fs+BXqFXAA4SbhByHHfr1q0ZM2bwm/CzZ8+urq4WculslZeXR0REENHMmTNxiodzoVfoFYAjBB2EvF27dvGPQ+vateuhQ4eEDyC8oqKi+Ph4Iho9erTRaGQdxzOhV6zjAIiVE266bav58+dfvHhxzJgx/Ib8smXL7t+/L3wMATQ0NOzfv3/GjBn9+vVTqVRxcXH//e9/fX19WefyTOgV61wAosVqApvN5oyMDP4If3R09MmTJ1klcYXTp08vW7YsKCiI/4/s5+f3+uuv19TUsM7l+dArALAVs0HIKywsHDBgABF5eXlpNJqGhga2eRxUUVGRkZHRv3//5u2MwYMHZ2Rk3Llzh3U0aUGvAMB6jAchx3H19fUajUahUBDRU089JcY7YphMpuzs7NmzZ3t7e/Prqa5du65cufL8+fOso0kXegUAVmI/CHnffPNNz549icjX1zc9PV0sD465fPmyRqMJCwvj11MKhSIxMTErK6uxsZF1NOA49AoArOAug5DjuNra2ubroydMmODOtww2GAyZmZkjR45s/qkqPj4+PT29srKSdTR4FHoFAG1zo0HIO3z4sFqtJiKVSpWZmck6zq+Yzebs7OyUlBR/f39+PRUUFJSampqXl8c6GrQDvQKA1rjdIOQ4rrKyMikpiV8jJCUlZWdnl5eXs72ZZHFxcVpaWlRUFJ9KLpcnJibu2rULD8YTEfQKAFok4ziu7esrWNm/f/+yZctCQkJKSkqISKlUdu/eXa1Wh4eHR/9MrVb36NGjeTva6YxG46FDh7Zt23b06FH+P1SvXr1eeumlhQsXNq+8QFzQKwB4hPsOQiIqLS3NycnZunWrTqerrKxs8T0KhUKtVkdFRUVGRkZGRkZERPB/REVF+fn52b3os2fPbtu2bd++fXV1dUQUGBiYlJQ0f/78CRMmyGQyu78W3AF6BQAPc+tB+LCGhoby8nKtVnvz5s2Kigrtz/R6vclkavEjwcHBj2zp8xv74eHhrS2lvLz8o48+2r59+7Vr1/hXBg8enJqaOnfuXP72XeBh0CsAEM0gbI3JZCorK9PpdDqdrrS0VK/X83/r9fqGhoYWPxIYGBgREdG8sR8ZGalWq0tKSvbu3Xvs2DGLxUJEkZGRCxYsWLBgQXR0tLD/IHAL6BWAdIh+ELbBYDDwW/cPb+yXlJTU1NQ8/mZfX9/6+npfX98ZM2akpKRMnTqVvxYb4BHoFYCH8eRB2BqDwaD7WfPG/pQpU7p16zZnzhyVSsU6IIgSegUgUlIchAAAAM0YPIYJAADAfWAQAgCApGEQAgCApGEQAgCApGEQAgCApGEQAgCApGEQAgCApGEQAgCApGEQAgCApGEQAgCApGEQWuv06dOzZs1inQLExGw2FxUVNTU1sQ4CAG3BILRWQEDAk08+yToFiMndu3f79OlTUVHBOggAtMWLdQD3VVlZuXPnzgcPHkyePHnEiBEqlapfv35EdO3atQMHDhiNxueee45/BaBFe/fuJaLdu3e/8MIL33333eDBgw8ePKjRaCoqKg4fPqzX64cPH56YmOjt7a3T6U6dOtWnT5/PPvtMpVItXLgwMDAQTQMQBvYIW1ZdXT1w4MDDhw+XlpYmJibm5ORcvHjxlVde+fHHH/v166fX67/99tvhw4eXlpayTgru6+rVq0T0ww8/FBcXL1u2bOjQobm5uffu3Rs6dOjOnTvv3bu3ZMmSv//970R06dKlJUuWTJ8+XavVvv7660uXLkXTAASDPcKWbdmyxd/f/+jRo97e3n369Dl8+PDYsWOJ6MiRI0FBQWvXru3UqdPmzZtbfBYrAG/t2rX/+Mc/3nnnnStXrjx48ODQoUPTpk0rLi6eNGnStm3b5HJ5QEDAxx9//MorrxBRQ0PD0aNHY2Njx48fv3r1ajQNQDDYI2zZpUuXEhISvL29iUij0bz77rv864sXL+7YsWN4ePi4ceNMJlPfvn2ZxgTR6Nix49SpU4koNja2f//+EydOjI6Ofvvtt5sfCBoSEhIbG0tETzzxRE1NDZoGIBgMwpap1eqffvqJ/7uoqOjo0aP832azuaCg4MSJE+PHj1+3bh1/EAigXd7e3jKZjIg+//zzP/7xjwsXLvzqq682btzY/AalUvnw+9E0AMFgELbs+eefz8nJOXPmTE1Nzfz58wsLC/nXN2/ePGfOnLi4uDVr1vTt21er1bLNCe5MLpcT0b179x5+8dy5czExMfPmzQsJCfn3v/9tsVha/CyaBiAYDMKWjR49etKkScOHD+/WrVtISEhqaioRyWSypUuX6nS6rl27dunSpamp6fe//z3rpOC+OnXq1L9//wEDBty+fZvfHSSiBQsWGI3G0NDQ6Ojo4cOHFxUVffDBB0TU/AaZTIamAQhJ1nyIAh5XVVXl5eUVEhLy8Iscx50/fz4sLKxbt26sgoFYmEymmpqa0NDQh1/kOK64uDgmJsbb29tgMPj4+Pj7+z/+WTQNQBgYhAAAIGn4aRQAACQNgxAAACQNgxAAACQNgxAAACQNgxAAACTt/wBNNniCLw9bAwAAAABJRU5ErkJggg==\n", | |
"text/plain": [ | |
"<PIL.PngImagePlugin.PngImageFile image mode=RGB size=600x200 at 0x12084DCD0>" | |
] | |
}, | |
"execution_count": 23, | |
"metadata": {}, | |
"output_type": "execute_result" | |
} | |
], | |
"source": [ | |
"Draw.MolsToGridImage([m2_cis,m2_trans],legends=[\"cis\",\"trans\"])" | |
] | |
}, | |
{ | |
"cell_type": "markdown", | |
"metadata": {}, | |
"source": [ | |
"Not sure why they are depicted identically, they shouldn't be. Issues with depiction..." | |
] | |
}, | |
{ | |
"cell_type": "code", | |
"execution_count": null, | |
"metadata": {}, | |
"outputs": [], | |
"source": [] | |
} | |
], | |
"metadata": { | |
"kernelspec": { | |
"display_name": "rdkit-kernel", | |
"language": "python", | |
"name": "my-rdkit-env" | |
}, | |
"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.7" | |
} | |
}, | |
"nbformat": 4, | |
"nbformat_minor": 4 | |
} |
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment