Skip to content

Instantly share code, notes, and snippets.

@greglandrum
Created January 5, 2021 06:32
Show Gist options
  • Save greglandrum/93a209c0155d82f0fde07d14b4a14fbe to your computer and use it in GitHub Desktop.
Save greglandrum/93a209c0155d82f0fde07d14b4a14fbe to your computer and use it in GitHub Desktop.
github3669.ipynb
Display the source blob
Display the rendered blob
Raw
{
"cells": [
{
"metadata": {
"trusted": true
},
"cell_type": "code",
"source": "from rdkit import Chem\nfrom rdkit.Chem.Draw import IPythonConsole\nfrom rdkit.Chem import Draw",
"execution_count": 14,
"outputs": []
},
{
"metadata": {
"trusted": true
},
"cell_type": "code",
"source": "m = Chem.MolFromSmiles('CCCCn1cccc1')\nIPythonConsole.drawOptions.addBondIndices = True\nm",
"execution_count": 11,
"outputs": [
{
"output_type": "execute_result",
"execution_count": 11,
"data": {
"text/plain": "<rdkit.Chem.rdchem.Mol at 0x7fb9de0981c0>",
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAcIAAACWCAIAAADCEh9HAAAABmJLR0QA/wD/AP+gvaeTAAAgAElEQVR4nO3deVxTV9oH8CcJS9gXcUFWBTd0BLF1Q3FDeVWqVg1qFbVjS61tsbWvxs+nWtqp00987cxQp7UTtxEVsSBF0bpRVEBc6m5RtAgoIhQbWcKSQELO+8ehkSJLIDe5N/B8/+hHmOTcBzv8enLPuc/hEUIAIYRQZ5mxXQBCiKN27NiRmJjY7Jve3t47d+5kpR7OwhhFCLXM3d195MiRTb+ze/dujUbDVj2cxcMP9QghXWRnZw8fPvzw4cPz5s1juxZuwRhFCOnkr3/96/nz53NzcwUCAdu1cAuf7QIQQiagtLQ0Pj5+7dq1mKEvwxhFCLVv27ZtVlZWK1asYLsQLsIYRQi1o7a2ViqVrlq1ytbWlu1auAhjFCHUjj179sjl8tWrV7NdCEfhEhNCqC0ajWbQoEFBQUF79+5luxaOwn2jCKG2HDly5OHDhy/vw0daOBtFCLUlKCjIxsbmzJkzbBfCXTgbRQi1SqFQLF++PCgoiO1COA1nowghpBdcqUcIIb1gjCKEkF4wRhFCSC8YowghpBeMUYQQ0gvGKEII6QX3jSKEDKKqqkqtVldXV6tUqpqamvr6+tra2rq6OoVCoVQq6+rqamtr6+vra2pqVCpVdXW1Wq2uqqpqaGiQy+WEkIqKCk9PT6lUambG9ZjCfaMIdWs01PT/Z7OhKisr9T9uxM7OzsrKKjc3197enomf1VC4HvMIoWbu379fVFRE04rGVmszO41GU1lZSWd2AFBeXg4AFRUVhBBGYq5ttra25ubmNjY2FhYW1tbWlpaWVlZWQqHQ0tLS2trawsLCxsbG3Nzc1tbWzMzMzs5OIBDY29vz+XwHBwcAcHJyWrdu3a1bt3744QeO9znF2ShCpuSnn35q8cDOThMKhTTdOv3PFgdxcHDg8/Vdetm1a9fbb789derUn376iYmf1VAwRhEyJQEBAbdv3w4MDOzbt6+VlVXbMzsej+fo6AgATk5OAODo6Mjj8WjA2dvbc/84ELlc3qdPn7q6ukePHnl4eLBdTqswRhEyGcXFxe7u7jY2NjKZzNLSku1yjCE8PDwxMXHLli3r169nu5ZW4YYnhEzG8ePHCSHTpk3rJhkKAEuXLgWAffv2sV1IWzBGETIZJ0+eBIAZM2awXYjxzJgxw8XF5e7du3fu3GG7llZhjCJkGurr69PS0gAgNDSU7VqMx9zcPDw8HAAOHDjAdi2twhhFyDRkZGRUVVX5+/t7enqyXYtR0c/1cXFxDQ0NbNfSMoxRhEzDiRMnAGDmzJlsF2JsY8eOHThwYHFx8blz59iupWUYowiZhm4bowCwePFiAIiLi2O7kJbhhieETEBBQUH//v2dnJyePXvG/WfMGZeXlzdgwAAbG5vffvvNxsaG7XKaw9koQibg+PHjABAaGtoNMxQAfHx8Ro8eXV1dnZKSwnYtLcAYRcgEdMOtTs3QhSZurtfjh3qEuE6hULi4uCiVyuLi4t69e7NdDjvKyspcXV01Gk1RURHX/hK64wcEhEzL2bNna2trR40axbX4aLF3lFwub2hoaLvZKG1MFRoaunLlSh2v5ezsHBoaeuzYse+//z4qKsqAP1XHYYwixHX6r9Ez1VS0xTajnebo6Kh7jALA0qVLjx07duDAAa7FKH6oR4jrfHx88vPzf/7551dffVXHt5w4cWLJkiV0qmhQTXtH0a5RdnZ2ZmZmrTUbpW30aGMqb29v3X8iAFAqla6urhUVFTk5OYMHDzbYz9RhOBtFiNPu3buXn5/fs2fPkSNH6v4uPp/fNEP1byr68j/pHwzwE7dKKBTOmzdvz549cXFxX3zxhTEv3TacjSLEaV999dW6deuWL1++d+9e3d9Fu9/TqWJXcu7cuSlTpnh5eRUUFPB4PLbLaYQbnhDitM5tdTIzM+t6GQoAkyZN8vb2fvz4cVZWFtu1vIAxihB3yeXyCxcuCASCadOmsV0LJ/B4vEWLFgHHNpBijCLEXampqfX19ePGjXN2dma7Fq6IiIgAgISEhLq6OrZraYQxihB30U/03bMdSWv8/PxGjBhRXl5O94FxAcYoQhxFCDl16hRgjL5kyZIlwKXP9bhSjxBH3bx5MzAw0N3dvbCwkDur0lxQUlLi4eEhEAiKi4t79OjBdjk4G0WIq7QPL2GGNuPq6jp16tT6+vrDhw936I1ZWVlvvvlmUFBQSEjIpk2bnj9/zkg9GKMIcRR2dWqD9mQR3d/yr3/9a8KECbm5ucHBwWPHjr169aq5uTkjxeCHeoS4qKysrHfv3nw+XyaT2dnZsV0O59TU1PTp06empubhw4f9+/dv9/XZ2dkjRoxYsWLFzp07GS+G+dloenp6VVUV48Mi1K2cOnVKrVYHBwdjhrbIxsZmzpw5hJCDBw/q8vo9e/bw+XwDPULKcIzeuHFj5syZ48aNe/ToEbMjI2SKCgsLv/vuu88++2znzp0ymUz3N+JWp3bRz/X79u3T5cWXL18eOHCgtbV1YmKiRCJJSkqqr69nrBTCqLy8PD8/PwBwcXHJyMhgdnCETEtiYqKlpeWQIUPmz5/fr18/Ozs7HX8pGhoaevXqBQD37983dJGmS61W9+nTBwB+/vnndl/s5eXl7e3dt29fDw+PgIAAgUAQGBhYWVnJSCUMxyghRC6Xz5o1CwAsLS3379/P+PgImYq+fftOmTKloaGBEKJUKv38/F599VVd3njp0iUA6Nevn4ELNHkffvghAKxZs6bdV3p5ebm7u2dmZtIvU1NT+Xz+xo0bGSmD+XujdnZ2R48e/eCDD+rq6pYtW7ZhwwaNRsP4VRDiOHrmR0BAAJ/PBwBLS8uAgAAdb3ZdvPjE2bkHnY6gNtDP9fHx8SqVqu1X9u3bt3///uPHj6dfhoSEDBo06MKFC4yUYZANTwKBYNu2bVKpVCAQbNmyZdGiRQqFwhAXQoizhEKhv79/bGxsamoqAJSVlZ09e3b27Nm6vDc+XlRZWfraa5sNXKPJGzly5NChQ589e0b/ktvQr1+/nJycpmmrUCisra2ZqYOROW1rTp8+7eDgAABjxowpKSkx6LUQ4prs7GwPDw8A8Pf39/PzCw8PVygU7b6rtJTw+cTKitTUGKFGk/f3v/8dABYvXtz2y5KTkwFgy5Yt9Mtdu3YBwHfffcdIDYaNUULIL7/84u3tDQBubm43btww9OUQ4o5//vOfdnZ2X3zxxeLFi62srBwcHA4dOtTuu/bsIQAkLMwIBXYFjx8/5vP5VlZW9EC91mg0muXLlwPAoEGDAgICAGD+/PlqtZqRGgweo4QQmUwWHBwMALa2tkePHjXCFRFi3dWrV3k83vbt2+mXv//+e0hIiIWFRUFBQdtvFIkIAPn2W4NX2GXQeNm7d2+7r0xLS/vss8+io6PPnDmj0WiYKsAYMUoIUSqVtEugQCCQSCTGuShCLPr6668BIC8vT/ud9PR0AEhMTGzjXSoVcXQkAKTJ+1A7duzYAQAhISFsFWCkZ+otLS1jY2NpgG7YsOHtt99ud2UNIZNG72Vdu3ZN+51bt24BgLu7exvvysqCigrw8wMdnm9EjcLDw62srM6ePVtUVMROBUaO7cTERLo6FhISUl5ebuSrI2Q0tbW1f/nLX+zt7Tdu3BgXF7dhwwZra+tZs2a1/S6xmACQ//1f49TYdSxYsAAAtm7dysrVWWhNcuXKlblz5/7222++vr7Hjx8fNGiQkQtAyDiqq6u//fbbtLS04uLinj17zp07d/Xq1W13FRo+HH75BdLSYMoUo5XZFRw5cuT111/39/enU34jY6fD09OnT2fPnn3jxg1nZ+ekpKRJkyYZvwaEuObJE/DyAltbkMnAwoLtakxKfX29m5ubTCa7ffv28OHDjXx1dvqNurm5paenz5kzp6ysLDQ0tEMHcCPUVf34IxAC06ZhhnaYhYWFSCSCDnYgZQprbZttbW2TkpLEYnF9ff2bb765Zs0afGYUdXMnTwIAYFOnzqEPhh44cKChocHIl2a/bfOuXbtWr16tUqlef/31/fv329jYsFsPQqyoqwMXF6ipgSdPwM2N7WpM08CBA3Nzc9PS0qYY99Yy+4eIvPXWW2fPnnVxcUlOTg4KCiosLGS7IoRYkJ4O1dUQEIAZ2nmLFy8GNk4MZT9GAWD8+PGXLl0aNGjQ7du3x4wZ03SrHULdBH6i119ERASPx0tKSqqtrTXmdTkRowDg6+t78eLFyZMnl5SUTJgwIT4+nu2KEDKqEycAAPD8On34+vqOGjVKLpenpKQY87pciVEAcHZ2PnPmzLvvvqtUKpcsWfLZZ5+xXRFCRpKfD7/+Ck5OMHo026WYuE6cGKo/DsUoAJiZmW3fvj0mJobH433++eeLFy9WKpVsF4WQwR0/DgDwP/8DZmZsl2LiFi5caG5ufurUqdLSUqNdlFsxSq1Zs+b48eP29vaHDh2aOnXqs2fP2K4IIcOin+jxxqj+evbsGRoaqlarExISjHZRLsYoAMyYMSMzM9PT0/PixYtjx47NyclhuyKEDEWhgIwM4PNh+nS2S+kStBtIjXZF9veNtqGkpGTOnDlXr161s7OLj4/Ho2k6pLy8PCEhQaFQuLm5zZs3TyAQsF0RatmxYzB7NowZA5cusV1Kl6BUKvv06VNZWXn//n3jtOzg6GyUcnV1zcjIeOONN6qqqubMmfPNN9+wXZEJqKurO3bs2LJly9zd3VetWvXJJ5+Eh4e/9tprcrmc7dJQy+hWJ1yjZ4pQKJw3bx4AHDx40EiXZKWvVIdoNJro6GgejwcAkZGRKpWK7Yq4SKPRZGZmvvPOO05OTvTfrEAgmD59+qeffkpPPB82bFi7fdcRK7y9CQC5epXtOrqQtLQ0APD29mawxX0bTCBGqfj4eKFQCAChoaFtH7rS3eTk5ERHR/v6+mr/0+jn5yeRSIqLi+kL8vPz/fz8AKBHjx7p6ensVouayc4mAKRXL9LQwHYpXUhDQwM9TPDChQtGuJzJxCghJCsrCydWWs+fP5dKpUFBQXSeDgBubm5RUVEtnhsol8vDwsIAwNLSct++fcavFrXm//6PAJAVK9iuo8tZv349AKxatcoI1zKlGCWE5OXl0YmVi4tLRkYG2+WwQKFQpKSkiEQibfdfBweHiIiIlJSUto85VKvVUVFRAMDj8cRicQNOfrhh0iQCQBIS2K6jy7l79y4AODk5KZVKQ1/LxGKUECKXy+mSvaWl5f79+9kux0gaGhoyMzMjIyPt7Oy0tz5DQkJiY2Orq6t1H0cqldL8XbBgQQ2eg862ykpibk7MzAiep2MI/v7+AJCcnGzoC5lejBJC1Gr1Bx980E0mVtnZ2dHR0fR8NGrkyJExMTGlpaWdG/D06dOOjo4AMHr06JKSEmarRR2SmEgAyMSJbNfRRW3duhUA5s+fb+gLcShGKyoqqqqqdH+9VCo1MzMDAJFIVFtba7jCWFFUVBQTExMYGKhNT09PT7FY/ODBA/0Hf/DgwYABA+jt1OvXr+s/IOqcv/6VABA8cdxAnj59KhAILCwsnj9/btALcSJGf/zxR/pbDQCvvPLKVZ23fpw+fdrBwaErTaxqa2sTEhLCwsLM/ni42snJKSIiIjU1ldmtGzKZbOLEiQBga2t75MgRBkdGOtJoSN++BIDcucN2KV3XtGnTAGDHjh0GvQr7MXrjxg0LCwvacjQtLW3QoEHOzs66Z+Ivv/xCP/Ca9MRKrVanpqZGRERom/9bWlqGhYUlJCTU1dUZ6KJKpXLZsmX03kh0dLSBroJac/06ASDu7sQoWxu7KXrOW3BwsEGvwn6MhoeH29raas+sv3fvHgBs3LhR9xFkMllwcDCdWB09etQwZRpKdna2WCzu3bs3TU8+nx8UFBQTEyOTyYxTQExMDJ/PB4CVK1fW19cb56KIECKVEh6PREayXUeXVl1dbWtry+Px8vPzDXcV9mO0Z8+eCxcubPqdESNGBAUFdWgQpVIZERFB168lpnCrqbCwUCKRDBw4UHvrc8iQIdHR0Xl5ecYvJjEx0draGgDGjx//+++/G7+Abqu4mDx6xHYRXd0bb7wBAJs3bzbcJViOUYVCwePxPvnkk6bfFIlEXl5eHR1Ko9FIJBI6sXrrrbe4ObEqLy+PjY0NCQnR7pl3dXWNiorKzMxkt7CbN2+6u7sDgK+vb05ODrvFIMSgEydOAMDgwYMNdwmWY/TJkycAsHXr1qbfjIyMtLKy6tyA2olVSEhIWVkZEzUyQKlU0j3zFn8cQG5lZSUSiVJSUrjTIqCoqGjkyJEA4OzsfPbsWbbLMXn+/gSAAJAePUjbU/xevRpf2eX2m3CCSqXq06cPAFy7ds1Al2A5RunJU5s2bWr6zYULF3p4eHR6zMuXL9O/NV9f3/v37+tdo16uXbsWFRXl4uLS9NanVCrt0NYuo6murp4zZw78cQwB2+WYNm2MApC33mrrlRijulOpVB999NHx48c79K41a9YAwIcffmigqti/N+rk5LR8+fKm3xk3btyYMWP0GbOoqIjuuGRrYkXbhfj4+DRrF8L9XVm0nxatOSoqqms/2mBQTWOUxyPnzrX6SoxRHRUXF0+YMAEAevfu3aFn8K5evQoAvXr1MtCHP/ZjdNasWS4uLtptPXTH7Nq1a/UclrYoBQALC4v//ve/+lapG227EG16uru7R0VF3bx50zgFMGX37t30/sOMGTMqKyvZLsckaWOUxyMAZNgw0trteoxRXWRkZLi6utKtjRcvXuzo24cOHQoAJ06cMERt7MdoWloaj8dbtGjR8+fPi4qKgoODraysGNmdYLSJlUKhoHvmm7ULYXzPvDFlZmb27NkTAIYPH/748WO2yzE92hgViRr/sGVLy6/EGG2XthfExIkTO/eRbvPmzQCwZMkSxmsjXIhRQsh//vMfW1tbGkC9e/f+8ccfGRx8165d9F/AzJkz5XI5gyO30S6ka3T9ePjw4eDBg+l2gp9//pntckyMNkZzcoidHQEg1takxekBxmgbqqqqwsPD6UMiUVFRnf5U/ujRIx6PZ21tzWwIUJyIUUJITU3N5cuXr1+/boiNSpmZmXSRh6mJFd0zTz9iMNIuhLPKysqmTJkCAEKhMC4uju1yTIk2Rp8/J5991vjnsLAWXokx2pr79+/TD+P29vaHDx/WZyiNRuPj42NhYZGWlsZUeVpciVFDy83N1X9iRduFjBgxQpueXl5eYrH4119/ZbZaTlGpVKtXr9Y+M2q6tymMrGmM1tY2nhQCQH74ofkrMUZblJycTDtmDBky5N69e/oMJZfLFyxYAADDhg0zxCpTd4lRQsjz588nT55MJ1YHDx7U/Y0VFRWxsbFhYWHawzWdnJwiIyMzMzO7T6ZonxldtGhR1+unZQhNY5QQkpDQ+KWrK2l2CA7GaDMqlUosFtPftSVLlnSoo+7LGJzStqYbxSghRKVSvfvuuzo243i5XYhQKKTtQhi881BTU5Oammq0J+j1ceLECXt7ewAYO3Zs17t9wbhmMUoImTy58TsfffSnV2KMNlVaWkqnO2ZmZvo/2M3glLYN3StGqaYTK4VC8fIL6J55eu5T0z3zzG79KS0tjY6Opjdtv/76awZHNpw7d+54eXkBQP/+/e/evct2OZz2cow+eEAsLQkAEQhI0xOzMEa10tPT6ZKDu7t7J3Y1NcXslLZt3TFGSZOJ1bhx47QTq8ePHzdrF+Ln5xcdHW2I3jBqtdrNzW3KlClbtmwxoRglhBQXF48aNQoA7OzsOvowSbfycowSQsTixm+OGvXiKFCMUUKIRqOJiYmhm2omTZr022+/6TPa06dP6fZtS0vLmJgYpopsTTeNUULI7du3PT09AcDb23vjxo1jx47Vpmffvn0//vjjW7duGbQAetJWdXW1acUoIUShUNCuOQKBwLQqN6YWY7Sm5sVa065djd/EGJXL5SKRCP44FqjtwxnbxeCUVkfdN0YJIU+fPn3llVf4fD5dO2KlXYgpxij549EG2qcqMjKSOw1WWKdddGwxRgkhSUmN3+/Vq3GtqZvHaE5ODj3u197ePikpSZ+hmJ3S6q5bxyghpKam5ujRo5s2bYqLi2Nlz7yJxih16NAhKysrAJg+fXpFs+Xn7ic7m4jFxN2d0KdsWotRQsisWY3/07p1hHTvGD148CBdwvX393/48KE+QzE7pe2Q7h6jrDPpGCWEXLx4kbbuHzp0qEEbjHPWo0fk738nQ4a8aERCD/5pI0YfPiRCIQEgFhYkL6+bxmjTJaClS5fqOYlhcErbCRijLDP1GCWE5Ofn0315PXr0SE9PZ7scI6moILGxJCSksfMIAHF2JpGRJDOz8XN9GzFKyIvnmhYt6o4xWlRUNG7cOKaWgBic0nYOxijLukCMEkLkcnlYWBj9rYiNjWW7HANSq0lqKomIINbWjdknFJKwMJKQ0LyBU9sxqlAQX9/G/k/m5n+K0X37iLHu6bEjPT2ddgR2d3e/dOmSPkMxO6XtNIxRlnE2Ruvr6wsLC1vcV9sitVpN/w9N70x1vUal166RqKgXM0c+nwQFEamUtNbpou0YJYScPv3iPoA2RrOyiEBAHB3J/v2G+1FY03QJaPLkyXouATE7pdUHxihrbty4IZVKt23bBgALFy6USqXnz59nuyhCCKmtrX3//ffpWSx8Pn/69OkFBQU6vlfb0GzBggVdo83Vo0dEIiEDBrwIOz8/Eh1N2v0raTdGCSHz5zeP0eJiEhb2osNeVzpgUPtgO1O7mpia0uoPY5Q127dv7/9nYrGY7aIIISQ8PLxHjx4HDhzIzc09cuSIi4vL6NGjde8ecObMGUdHRwAICAh48uSJQUs1nLIyIpWSoKAXtz7d3EhUFLl+XdcRdInRJ0+Ire2fYpQQotEQqbSxt17v3iQ5mYEfh3VNl4B+eLk7S0cwO6VlBMYoau7IkSMZGRnaL7/88ksA6NAq/K+//kofBnNzczPcOWKGoFCQlBQiEr24X2llRUQikpJCOro1VpcYJYRs2dI8RqmCghfP4ItEhDPHM3ZGXFwcXQIKCAjQf1cTg1NapmCMonbExcUBQFZWVofeJZPJJk6cCAC2trZHjhwxUG1MaWggmZkkMpLY2zcml0BAQkJIbCzp9KPYly6R1FSSmtrq2SFUfX3jy1JTSbP7yXRaamPT2BfKFJ+8NdyuJj2ntMzCGEXtWLdunUAgeN72nKolSqVy2bJlOvbTYsvduyQ6mvTr96dbnxIJh9bKHz4kEyY0rulHRhJOHinbMmaXgBic0jIOYxS1paSkxNnZOTIystMjaPtprVy50hBHG3TO06ckJoYEBb1ITw8PIhYTtg/kbplKRSSSxu5Q/fq1dcgod5w/f54uAXl4eFy+fFmfoZpOaSMiIji4dIkxilpVU1MTHBzs6en57NkzfcY5fPgwXfcfP37876yuPcvl8tjY2KVLZQLBiz3zq1aRCxcI9xtw//ILCQx8MS3lXpg0arYEpGdrWu7samoDxihqmUwmmzRpkqurKyNHpNy8edPDwwMAfH19c3Jy9B+wQ7QduOnJicHBiZaWjXvm/zjY2zTQaSld/hoyhHDwmEG5XD5//nymloAYnNIaFMYoasG1a9e8vb2HDh2am5vL1JhPnz4dOXIkADg7O589e5apYdt29erVNWvW0Kf+6e92cHDwnj2J5eXGub5B3LrVuA3AzIyIxUSpZLugP9y7d2/IkCGG2NXE8dMWMEZRc7t37xYKhUOHDj1//vy1PzCyO6+6unru3Ln0fIhvv/1W/wFbU1hYGBMT4+/vr+0hO3jw4OjoaK4tTXSaQkHEYkJvTQwb1oHdrIZz4MAB7RJQXl6ePkMxO6U1AoxR9Cc1NTXQkn/84x+MjE8bldIxo6KimH1mtLy8PDY2NiQkhDZCpTNfevggg1fhjosXycCBBICYmxOxuJ2dVYbD7BIQg1Nao8EYRX+i0WjKWqL7w/W62L17t4WFBQDMmDFD/xOulEplSkpKREQEXcgCAKFQSDtwc2dvgIHU1hKxmPD5jQeTGP22M3ny5Ak9OYKRJSAGp7TGhDGK2HHhwoWePXsCwPDhwx8/fty5Qejhg3QcaHL4oLy1fiFdVGoq8fRsbDclkRCjfQg+d+4cve/cHXY1tQFjFLHm4cOHgwcPBgBXV9efO7Lq/OjRI4lE4uvrq73nQA8f1L2FStdTWUkiIxs7AIwbR5jYXtEWjUYjkUjo6TszZszoxNMZTTE7pTU+jFHEprKysilTptCP4XFxcW2/+Pnz51KpNCgoSHvr083NLSoq6jpzKyx5eXnx8fGHDh0y/q4sRpw8SdzcGlsBSCTNny5lSmVlZdMlID1vcDM4pWULxihimUqleu+997TPjL7cSkqhUKSkpIhEIrr9hS4+REREpKSkMLiGq1aro6Ki+Hy+q6uri4sLj8dbuXKl7n2tuKO8nERGNj5cMG0aKSxkePy8vDz6OcDZ2fnEiRP6DMXslJZFGKOIE2JiYuiv08KFC2trawkhDQ0NmZmZkZGR9vb2ND0FAkFISEhsbGx1p/uFtG7v3r1OTk60A4tGo9m6dSsAGP9UH6YkJBAXFwJAHBzIvn1MPmOgUCgCAwNHjBih5xIQs1NadmGMIq44efIkTcwRI0a89957bm5u2lufY8aM+fe//23QB0kbGhqaPmtQWloKAH/7298Md0VDKy0lr79OAMiECStnzJhRVFTE1MhPnz7Vc+fGrVu3fHx86JT25MmTTBXGFoxRxCF37tzx8vLSTj89PT3FYvGDBw+MX0lsbCwAnD592viXZtahQyVOTk4A4OLi8v3337NdDiGE7N+/n25N031Kq1Aomu3AU3W0/6shYYwibikuLr5y5cratWuzsrKMf2uyoKBg3bp1c+fOdXBw+Oqrr4x8dQMpKSl57bXX6H+ZRCIRi91hlEplVFSUdldTrc5Hoa5du7bZ8yCMtGyjDQ4AAASqSURBVHpgilmLj6wgxBZXV1dXV9dRo0axcnW1Wi2Tyaqrq1Uq1YMHD+rq6iwtLVmphEF9+vRJSUnZt2/f+++/n5iYmJ6eLpVK6VO5xlRUVCQSiS5fviwUCrdt2/b222/r/t6SkpKRI0dKpVLtd/r162eAGjuL7RxHiIsuX74sEAg2b97MdiGtKiwslEqlUqlU9zaGBQUFkydPpr/4IpHImCvj2l1Nnp6eV65c6ejbg4ODFy1aZIjCGIExilDLBgwYMH36dLaraJlGowkJCbGzswOADj25oNFopFIpfeDS1dX12LFjhitSe0XtrqaZM2d2LrsHDBjw0UcfMV4bU/jszYMR4hCZTJacnKz9srq6+tmzZ3RxhoP27t17/vz5TZs2dfSNPB4vMjLyzp07EyZMKCkpmT179jvvvFNdXW2IIgGAnkC3YcMGjUYjFouPHTvm7OzciXFKSkqysrLGjh3r6+s7c+bMc+fOMV6qXtjOcYQ4Yfv27Twe75133jly5MjBgwfHjRtnbm7e9IRU7vj9999dXFzWrVt35swZ6OBsVEutVkskEnrn19vb+5wBTia5efMm3dXUo0cPfXY1aTSa9evXf/rppz/88MPevXvHjBnD5/NTU1MZLFVPGKMINUpISAgKCrKzs3NwcJgyZYohkoURIpHI09OzqqpKnxilsrOzaS9tOktl8LkG7a6mwMDADp3O3a6amho3N7epU6cyOKaeMEYRMiXHjx8HgJSUFEKI/jFKCFGpVBKJhPYt9PHx0b83a6d3Nelu5syZ3t7ejA/baRijCJmMyspKd3f3BQsW0C8ZiVHq9u3b9LAAMzMzsVis7OzJJIWFhWPGjAEAoVC4c+dO/QujZDKZ9s8qlcrHx2fSpElMDa4/jFGETMaqVassLS0zMjLy8vLy8vLoo1bJycmMbF1SKBRisZguqQ8bNqwTfbPOnj3bq1cvAPD09GQk3KmkpCRHR8etW7fevHnz9u3bb7zxBo/Ho/NxjsAYRcg0qNVqPr/lrTXr169n6ioXL14cOHAgAJibm4vFYh2PD2i6q2nWrFllZWVM1UMIqa2t/fjjjx0cHOgP6+XlFR8fz+D4+uMRQgyzBQAhxLD8/PymX164cGH58uXJycnBwcGd20jUIoVC8fnnn2/dulWj0YwaNWrv3r30cKTWyOXyFStWJCcn83i89evXf/nll63FvT7UanVBQYGFhYWnp6e24SxXsJ3jCKFOYvDe6MsyMzPpdiWhUCiRSFpr7Xrz5s3+/fsDQI8ePU6dOmWISrgPt98jZKp69+4tEokYnIc2NX78+Bs3bkRGRtbV1W3YsGHChAm5ubnNXrN///6goKD8/PzAwMCrV6+GhoYaohITwHaOI4Q47eTJk7T3q5WVlUQiof2VjbCryYTgvVGEUDsqKirEYvGOHTsAICQkZPPmzWvWrLly5YpQKPzmm29WrlzJdoEswxhFCOnk8OHD7777rkwms7CwqK+v9/HxSUpKortNuzmMUYSQrkpLS1etWvXKK6/cuHFj9+7djo6ObFfECRijCCGkF1ypRwghvWCMIoSQXjBGEUJILxijCCGkF4xRhBDSy/8Dk7wFm+4swgQAAADPelRYdHJka2l0UEtMIHJka2l0IDIwMjAuMDkuMQAAeJx7v2/tPQYg4AFiRgYI4ITiBkY2hgSQODOEZmJCp9kdNIA0MzObQwaYZsTD4Aaaz8jEwMTMwMzCwMKawMqWwcTGnsDOkcHEwZIgwsjKwsrGziEOcwQD56NZXQd8d4nagzh9V70P3DpyeB+IfVzo8f6mb4x2IPbWJQr7peTD94PY37s97CXqVoLVnJjC6OCtOQ0sLp/s63AttxEsLnNax+GBzRSwmUfXz7Q35ukDs8UAdC8uguTUHX0AAADgelRYdE1PTCByZGtpdCAyMDIwLjA5LjEAAHicnZIxDsIwDEX3nsIXwHJsJ05myoTowMAd2Lm/cJIqdAAJxbKq/+X65UftArXu6/X5glG8LgtAaU1fu5QCDyYifw9OiiJB6yZhzFq68inBGX4hjt0ogqHvnggls85RApop9d3AaTILYWKLXWUlmaL4hmXR/UYS6EDZ/qcwEqv1uwWTMJdFkKlwz8KR8xyFMdVP7+VZYp7MEjwL2U6hYjMUP5rbsxpXMowrHcZVHMZVGsaVtc1qXOUxye3H3yflA7gBXLZ1eQMr2n4m6AeEfAAAAJZ6VFh0U01JTEVTIHJka2l0IDIwMjAuMDkuMQAAeJwdjTEOAzEIBL+SMpE4xAI2IJd+yvV5wT0+OHSMZnf33vuLuw+v5305m8FBwiM9w2ldxkgv0CVsqTbRDBzhcjToLEQj4akxYMdLF0PSEo40y/nPGqTrlEU9ug2M6Claxiqnoh0dGr3U0qwqamUk+kcrMvMAqRj0eX7gICWRFBH5kwAAAABJRU5ErkJggg==\n"
},
"metadata": {}
}
]
},
{
"metadata": {
"trusted": true
},
"cell_type": "code",
"source": "fm = Chem.FragmentOnBonds(m,[3],addDummies=False)\nIPythonConsole.drawOptions.addBondIndices = False\nfm",
"execution_count": 12,
"outputs": [
{
"output_type": "execute_result",
"execution_count": 12,
"data": {
"text/plain": "<rdkit.Chem.rdchem.Mol at 0x7fb9de47a670>",
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAcIAAACWCAIAAADCEh9HAAAABmJLR0QA/wD/AP+gvaeTAAAd4klEQVR4nO3deVQUV74H8B+NQoPI4gJEMxgXIhCjNEsAgRyNYEC6I4gFggEzOqBxAbPiyyQSjZngOXlPcXsBdTLgBpSoYTUiGmVziNC4ACoaXIKigIAgDShd748yNT6jgN1N3+rm9zn5w3OnqfqeOZ6vfYtb9+owDAMIIYQUJSAdACGENBvWKEIIKQVrFCGElII1ihBCSsEafQmdnZ2kIyCEeAdrtG8NDQ3JyckSicTW1nbt2rV3794lnQghxCM6uODpRa5fv37o0KH09PQzZ87I5XIA0NXV7enp8fHxycnJ0dHRIR0QIcQLWKPPqq2tzcjIoGm6uLiY/T9HKBR6eXlRFOXm5ubp6Xn37t3vv//+k08+IZ0UIcQLWKNPVFZWZmVlZWZmFhUVsSOGhobvvPMORVEBAQHDhw9nB3Nzc/38/IYMGVJQUODi4kIuL0KILwZ7jVZWVtI0nZaWVl1dzY6YmZmJxWKJRDJnzpxhw4b9+UdWr14dHx8/ceJEqVTK1StCaNAapDXKtueBAweuXLnCjowcOXLOnDkURb377rt6enq9/GxXV5ebm5tUKv3ggw9+/PFHteRFCPHXIKpRuVxeXFxM03R6enpdXR07OHr0aB8fH4qifHx8hg4d2s9LXb161cHBoa2tbe/evQsXLhywyAghDaD9NdrT01NSUkLTNE3Td+7cYQetrKz8/f0pipo+fbpAoMiqr127dkVERBgZGZWVlb3++usqjYwQ0iRaW6NdXV0FBQWZmZkpKSn37t1jB8ePHy+RSCiKcnd3V37FUmho6IEDB5ycnIqKinp/DoAQ0mLaVqMymez48eM0TWdkZLS2trKDdnZ2EolELBZ7eHio8F6tra0ikai2tnbNmjXfffedCq+MENIgWlKjHR0d+fn5NE0fOXKkra2NHbSzs6MoKjg42NbWdoDuW1pa6uHh0dPTc/ToUW9v7wG6C0KIzzS7RpubmzMzM7OysnJych4+fMgOsu25cOFCa2trNWTYsGHDV199ZWFhUVFRYWlpqYY7IoR4RSNrtKmpKTs7m6bpY8eOdXd3A4BAIHBzc6Moav78+WPHjlVnGLlcPnv27Pz8fHxJFKHBSZNqtKGhITc3l6bpo0ePPn78GAB0dXVdXV0pigoKCnrllVdIBaurq7O3t29sbNy0adPq1atJxUAIEaEBNXrz5s3Dhw/TNF1SUsJuEaKvr+/p6SkWi0NCQszNzUkHBADIyckRi8V6enolJSUikYh0HISQ+vC3Rv+8RYiBgcGsWbMoivL39zc2NiYd8FmrVq3atm3bpEmTysvL8SVRhAYP3tVoL1uEzJs3z8jIiGy8XnR1dbm6ulZUVCxevHj37t2k4yCE1IQvNdrLFiF+fn6GhoZk4/VTVVWVs7NzR0fHvn37QkNDScdBCKkD4Rpl23P//v01NTXsSP+3COGnxMTEpUuXmpiYSKXS8ePHk46DEBpwBGq09y1CfH19hwwZouZIf9bc3GxmZqbYz4aEhKSkpDg7OxcVFfV/uxOEkIZSX41yW4SkpaXV19ezg+PGjZs7d64yW4SoFrem6tixY7W1tWPGjFHgIi0tLSKR6Pr161988cW3336r8pAIIV4Z8BpVwxYhyqutreWOXeIODqFpWiwWK3ZB7iXRn3/+2cvLS6VhEUL8MlA1ym0R8tNPPz148IAdZF/TlEgkjo6OA3HTl9XLmqq5c+eamJgoc/H169fHxsZaWFicO3fOwsJCRZERQryj4hrltgg5fPhwe3s7O8i254IFC2xsbFR4L4X189glJcnlcm9v7xMnTvj6+mZnZ/PhSzdCaCCopkbZLUJoms7Ly+vq6gIAgUAgEonEYrHatgjpkwLHLimprq5u2rRpTU1N8fHxUVFRKr8+QogPlKrRpqamgwcPpqennzx5knvJ/e233w4MDAwICFDs9zMqR3ZNVXZ2tkQiwZdEEdJiitfoo0eP9u7du3jxYnhqi5Dg4GA+bBbHqzVVK1as2LFjh7W1dVlZGb4kipD2UbxGq6ur7ezsjI2Nt27dKpFIFF5lqULPPXaJ+JoqmUzm4uJy587dJUt+jouzV38AhNCAUrZGbWxsuEeNpPB/TVVVVY2v7/CbNy1TUyEoiGwWhJCKkX9fSGEvWlM1EMcuKcnOzvq//gs+/BAiI8HZGfAdUYS0iebVqEasqfqzZcvgl18gNRUWLIDCQsB3RBHSGhpTo9yxS9nZ2R0dHeygmo9dUtIPP8CZM1BaCuvXwzffkE6DEFIRvtfoc49dcnd3J3LskpJMTWHvXpgxA/7xD5gxA2bNIh0IIaQKPK3R5x67xLYnT9ZUKcbDA778Etatg0WLoKICRo0iHQghpDR+1eiNGzeOHDny9LFLQqHQx8dHIpH4+/vz5NglJa1dC6dPw8mTsGgRZGUB6UUECCFl8aJGNe7YJWUIBJCcDPb2kJMD27fDypWkAyGElEOyRtnXNLOyssrKytgRTTl2SUmvvgo7d8K8efDpp+DhAfa4JB8hTUagRtn2TE1NvXTpEjvCbhFCUdTs2bP19fXVH0n9AgJg2TL44QdYuBB+/RU05KwphNBzqLtG6+vrp06dyj73tLCw8Pf3DwwMnDlzJh8ODlGzTZuguBjOn4ePPoKEBNJpEEKKUnd5WVpaUhRlYWERGBjo7u6uq6ur5gD8IRTC/v3g7AyJiTBzJixYQDoQQkgh2vBOvUbbsQNWrABTU5BK4bXXSKdBCL088qfIDXLLl4O/P7S0QHAwPHpEOg1C6OVhjZL3z3/CuHFQWgobNpCOgpR279499oURNHhgjZJnZgZ79oCuLmzYACdOkE6DFFVWVhYeHv7qq6/Onz8/MDCwtbWVdCKkJlijvODpCX//O8jlEB4OTU2k06CXIZPJdu3aZW9v7+TktGfPHoZhjh49eujQIRcXl/Pnz5NOh9QBa5QvYmNhxgyoq4NFi2BgDr1GKnb79u2vv/7aysoqIiLi3Llz5ubmMTEx165dO3/+/LRp0y5fvuzi4hIfH086Jhp4jKKqqqoAwMbGRuEroGfcusWMGMEAMNu3k46CelVQUEBRFLfY2dHRMSEhoaOjg/uATCbjzoJ9//3329raCKZFAw1rlF/S0xkARihkKipIR0F/IpPJkpKS3nzzTbYf9fT0KIoqLCx80ef37NnDvtM8efLkc+fOqTMqUiesUd6JjGQAGDs75uFD0lHQH65evRoTEzNixAi2QC0tLWNiYn7//fc+f/DSpUtTp04FAKFQuHnzZjVEReqHNco7MhkzdSoDwCxbRjrKoCeXy/Py8iiK4l63c3R0TEpK6u7u7v9FcIKv9bBG+ejiRcbAgAFgUlJIRxmsHjx4kJCQYGdnx9afvr4+RVElJSUKXzA5OXnYsGHsBP/8+fMqjIqIwxrlqa1bGQDG1JS5fp10lEHmypUrMTExpqambIGOGTMmNja2oaFB+StXV1ezz1Vxgq9lsEb5a+5cBoDx8GAePyYdZRDo6enJy8sTi8U6fxxI4O7unpaW9ujRIxXeRSaTRUZGstcPCwtrb29X4cURKVij/HX/PmNlxQAwsbGko2i1lpaWzZs3jx8/nm03oVAYFhY2oL9YT0pKYif4NjY2OMHXAlijvHbyJKOryzg5MS/zKw3UX5cuXYqKimIbDQAmTJgQFxfX1NSkhltzE3wDA4PExEQ13BENHKxRvsvJYTo7SYfQLj09PRkZGV5eXs/M3x+r9+lJR0fH3/72N5zgawGsUTSINDc3b968edy4cWx5DR8+PDIy8uLFiwQjcRN8W1vbCxcuEEyCFIY1qhnWr2fMzJ78t3ZtHx/+5JMnn/zssycjf/87Y2rKmJoyYWF93+vWrScfNjVlbt1SNjlPlJWVRUZGGv5x6JW1tXVcXNz9+/dJ52IYhqmqqpoyZQo7wd+5cyfpOOil4dYkmqGzE5qbn/wXFwdVVb19uKPjySdlsv+MtLRASwu0t/d9L7n8yYdbWqCnRwXhCeru7qZp2tvb29HRMTExsbOz08vLKyMj4/LlyzExMWZmZqQDAgDY2tqWlpYuWbJEJpNFRESEh4d3dHSQDoVeAtao5unuhmXLcBeoPty9e3fjxo2TJk0KCgo6fvy4sbFxZGRkZWVlXl6eRCLhnoryhIGBwa5du5KSkgwNDffs2ePk5FRZWUk6FOovrFGNVFAASUmkQ/BVWVnZ0qVLX3vttTVr1ty6dWvy5MmbN2++fft2QkKCjY0N6XS9CQ8PP3v27BtvvFFdXe3q6rp//37SibRfQ0NDYmKikhfBGtUwb7zx5A+ffYYbPP8/7Pzd3d3dyckpMTGxu7tbLBbn5eVVV1dHR0dzq5p4ztbW9syZM6Ghoe3t7QsXLsQJ/sCRSqXsP7dLly4tKipS5lJYoxpm5kyYPRsAoLERYmJIp+GHO3fubNy4ccKECUFBQcXFxaamplFRUdeuXcvMzHx6VZOmMDIy2rdvHzfBd3Z2xgm+Cj169CglJcXd3d3BwYF9XC4Wi7nfPSoGa1TzbNwIAgEAwD//CQUFpNMQxR5/NG7cuDVr1tTV1YlEooSEhLq6uvj4+Nc0/Ljq8PDwX3/91c7Orqqqys3N7cCBA6QTaTz2cfnEiRNDQkKKi4vZx+UXL17MzMwUiUTKXBlrVPPY28OSJQAADAMffjgYj2Xu6upKTk7mjj+Sy+Xs/L28vPzpVU2azs7O7syZMyEhIW1tbaGhoTjBV1gvj8ttbW2Vvz7WqEbasAHYHYgqK+G//5t0GjVijz969dVXFy1axB1/VFtby87fSadTveHDh+/fv5+b4Ht4eFy9epV0KI3x58fl7HI3lT8uxxrVSObmsHbtkz+vXw+//UY0jVoUFhYGBQWNGzdu3bp1jY2N7PFH169fj4uL+8tf/kI63cAKDw8vLCycNGmSVCp1cHBISUkhnYjv6uvrn35cbmJiwj4uH6DlbkNUezmkNitXwu7dUFkJMhmsWgXZ2f36qV9+AReXPj7T3a18OpVpb2/fv3//tm3bLly4AH8cfxQdHe3u7k46mlqJRCL2kUVKSkpISEhOTk5CQoKBgQHpXLxTVlYWHx+fkpLy6NEjABCJRMuWLXv//fcH9lGPwu8/4cug6vTFFwwAA8CsXPmfwdOnGR2dJ+MHD/5n/MMPn/3wRx89GXnZ/whuGq3w8UfaLSkpiW1PkUhUU1NDOg5fdHZ2JiUlTZs2jf3boquryz4uV8/d8duoBvP0hPffhz17AACio2H2bBg+vI8f0dcHE5M+PiOXQ2OjahIqgGGY/Pz8xMTEQ4cO9fT0AICjo2NUVFRISMjQoUOJxeKN8PDwKVOmBAcHS6VS9g3X4OBg0qFIun37dmJi4vbt2xsbGwHA3Nz8r3/96/Lly62srNQXQuECxm+j6vTcb6MMw9TXM6amT/6nNWueDPbybTQgoO973bhB5tvon48/CgsLk0ql6kugOVpbW7n2jIyM7ByUeykWFBRQFDVkyJPvguzj8o6ODvUnwV8xaTYLC/jmmyd//p//gZoaomkUVVNTs2bNGisrq6VLl1ZVVbHHH/3+++/sqibS6fjI2Ng4JSUlISFBX18/MTFx+vTp165dIx1KTTo7O5OTk6dOnerp6UnTtEAgoCiqsLDw7NmzkZGRRJ4X46Re4y1fDnv2QGkpdHfD55/D4cOkA/WbXC4/ceJEfHx8dnY2wzAA4O7uHh0dHRAQwH3FQL2IjIx0cnIKCgoqLy93cHDYuXNnUFAQ6VAD6Nq1azt37ty5c+f9+/cBwNLSctGiRatWrRo7dizZYPiXVeMJBLBlC0yfDnI5HDkCJ0+SDtQPra2t//rXv+Lj42trawFAKBRSFPXZZ5+x52qg/nNwcCgvL4+IiEhLSwsODs7Pz9+6dauenh7pXKrE8P5xOU7qtYGLC/xxGgV8/DHI5UTT9OrSpUvR0dFjx45dvXp1bW3txIkT4+Li6urqkpOTsUMVY2xsnJqampCQoKenx07wf9OWhcRtbW2JiYlTpkzx9vamaXrIkCEURZWUlJw9ezY8PJwnHQr4bVRrfPcdHDoEjY1QUQH37pFO8zwnTpxYv379qVOnAEBHR8fHx2flypW+vr4CAf5brgKRkZGOjo5BQUFlZWUikWjXrl0URZEOpbiamprdu3cnJCS0tLQAwJgxYyIiIlauXDlq1CjS0Z4D/wZriREjIC7uyZ9v3yYa5QVu37596tQp9vijCxcu5Obm+vn5YYeqkKOjo1QqpSjqwYMHwcHB0dHR3bx6laIf5HL58ePHJRLJ5MmTN27c2NLSwp42eOPGja+//pqfHQr4bVSbLF4Myclw+jTpHC9AURS7h+bwPle3IkWxE3x3d/fPP/98y5YtxcXFqampEyZMIJ2rb899XP7pp59OnTqVdLS+4XcB7aGjA9u2AW+eFz1LX19/2bJl2KEDTUdHJzo6uqioaMKECWfPnnVwcDh48CDpUL25fPny04/LJ0yYwD0u14gOBaxRLfPmm7B8OekQiAecnJzKy8vnz5/f2toaFBTEwwm+XC7PzMz09va2tbXdsmXLw4cP2fn7lStXnn4JWCNgjWoGoRDMzMDMDPpcXLx+PVhbP/thQ0MwNQVTUzAy6vteAsGTD5uagq6usskRKSYmJmlpaZs3bx46dOiWLVs8PDzY+TJxLS0t8fHxEyZMeO+9944fP25kZMRun1xYWEhRlK4m/p1T+P0nfBkUIY1QWlo6fvx4ADAxMTn49B42aldWVvb0vtrW1tZxcXH3798nGEklsEb5LieHGZQvTCNVamxs9PPzAwAdHZ2oqKju7m513r2rqystLY3bV1sgELDbJ8vlcnXGGDhYo7x28iSjq8s4OTHq/WuPtJBcLmcn+ADg7Oz822+/qeGm9fX1T++rzR5/VF1drYZbqxPWKH/dv89YWTEATGws6ShIW5SWlrKH/Y0cOZLdymCAsBuFCIVCtkDZ44/a29sH7o4EYY3y19y5DADj4cE8fkw6CtIijY2Nc+bMGaAJPjt/nz59Ojd/Z7dP1pr5+3NhjfLU1q0MAGNqSnL/eaStnp7gv/XWW7W1tcpf8/bt23FxcdxmS6amplFRUSq5Mv9hjfLRxYuMgQEDwKSkkI6CtNe///1vdoI/atSonJwcha9z9uzZsLAwbqMQkUiUkJDw8OFDFUblOaxR3pHJmKlTGQBm2TLSUZC2a2ho8PX1VWyCT/b4I17BGuWdyEgGgLGzYwbTP+eImKcn+G+//XZ/Dg2sq6uLjY3lNgoxNzePiYm5efOmGtLyE9Yov6SnMwCMUMhUVJCOggaTU6dOjRkzhp3g5+bmvuhj/Dn+iFewRnnk1i1mxAgGgNm+nXQUNPg0NDT4+Pg8d4Lf1taWkJDA7autp6fHHn9EMC2vYI3yRU8PM2MGA8D4+TFavTgE8dczE/y6urqrV68+vVGIpaVlTExMfyb+gwruN8oX69bBL7/A2LGQlAQ6OqTToEGJ3WRv2rRpoaGhp0+ftrGxefjwoVwuBwBPT8+VK1cGBATw5+gO/sAa5YWCAvj2WxAIIDkZRo4knQYNbjNmzDh37lxYWNioUaMOHjz43nvvffzxx66urqRz8RfWKHnNzRAWBj09sHYtvPMO6TQIAYwePTonJ6enpyc+Pn4k/sPeF6xR8hYvhhs34K234MsvSUdB6A8CgUAgEGCH9gdu20zYjh1w5AiYmkJqKn/P/0AI9YLAt9EFCxZYWFgEBga6u7tr5E7XqlNZCZ9+CgDwv/8Lr71GOAxCSDHqrtH6+nqapuVy+ZYtWywsLPz9/QMDA2fOnMkt6B08OjshNBRkMoiMhAULSKdBCClK3ZN6S0vL8+fPx8bG2tjY3L17NyEhYfbs2ebm5uHh4ZmZmV1dXWrOQ9BHH8H582BnB5s2kY6CEFKCDsMwiv1kdXW1nZ2djY1NdXW1YleorKykaTorK6usrIwdMTQ0fOeddyiKmjdvnlF/Tl/TWIcPw7x5oK8PZ86AvT3pNAghJZCsUU5tbW1GRgZN08XFxWweAwODWbNmURTl7+9vbGys5PX55vffwd4emppg61ZYuZJ0GoSQcnhRo5wbN24cOXKEpumSkhL23QmhUOjl5SWRSPz9/c3NzVV1I4LkcvDygpMnYc4cyMrCF5YQ0nj8qlFOQ0NDbm4uTdNHjx59/PgxAOjq6rq6ulIUFRwcbGlpqfI7qs3XX8O6dTB2LFRUwB87jSGENBhPa5TT1NSUnZ1N0/SxY8e6u7sBQCAQuLm5URQ1f/587sQCTVFYCDNmAMPAsWMwaxbpNAghVeB7jXKam5szMzOzsrKys7M7OjrYQTs7O4qiFi5caG1trYYMSmppAXt7uHEDvvwSvvmGdBqEkIpoTI1yOjo68vPzaZo+fPhwe3s7O8j26YIFC2xsbNQZ5qUsWACpqfDWW1BYiC8sIaQ9NK9GOTKZ7Pjx4zRN//TTTw8ePGAH7ezsJBKJWCz28PAgkupFfvgBPvwQTExAKoXx40mnQQipjgbXKKerq6ugoCAzMzMlJeXevXvs4Pjx4yUSCUVR7u7uOqR/HV5VVePrO/zmTcvUVAgKIpsFIaRiytaosbHx1q1bJRKJmZmZapMpoKenp6SkhKZpmqbv3LnDDo4bN27u3LkURU2fPl0gILAVi0wmc3FxuXPn7pIlP8fF4VJ7hLSN4jX66NGjPXv2LFmyBPi3GkkulxcXF9M0nZ6eXldXxw6OHj3ax8eHoihfX191vsK/YsWKHTt2WFtbl5WVDR8+XG33RQiph+I1CgBNTU0HDx5MT08/efIkt7rz7bffDgwMDAgIYA8aJI595XT//v01NTXsyMiRI+fMmUNR1Lvvvqunpzegd8/OzpZIJHp6eiUlJSKRaEDvhRAiQqka5bCrkWiazsvLY7cXEQgEIpFILBbzZzUS26dpaWncw1wzMzOxWCyRSObMmTNs2DCV37Gurm7atGlNTU3x8fFRUVEqvz5CiA9UU6McjViNVFlZmZWVlZmZWVRUxI5wW6IEBASoat4tl8u9vb1PnDjh6+ubnZ1N/NdcCKEBouIa5bxoNRJFURKJxNHRcSBu+rJ62RJl7ty5JiYmylx8/fr1sbGxFhYW586ds7CwUFFkhBDvDFSNcvi/GgkAamtrDx06lJ6efubMGfb/EKFQSNO0WCxW7IKlpaUeHh49PT0///yzl5eXSsMihPhlwGuUw61GSktLq6+vZweJr0Z6BrclyrFjx2praxX7LVlLS4tIJLp+/foXX3zx7bffqjwkQohX1FejHF6tRnqR5uZmhVfChoSEpKSkODs7FxUVDcW3PhHSdgRq9GlkVyMNhMTExKVLl5qYmEil0vH41idCgwDhGuX0shrJz8/P0NCQbLx+qqqqcnZ27ujo2LdvX2hoKOk4CCF14EuNcnpZjcTzA5q6urpcXV0rKioWL168e/du0nEQQmrCuxrlaNwBTatWrdq2bdukSZPKy8vxpU+EBg/+1ijn5s2bhw8ffvqAJn19fU9PT7FYHBISwpMDmnJycsRiMb70idAgpAE1yunlgKagoKBXXnmFVLC6ujp7e/vGxsZNmzatXr2aVAyEEBGaVKMcXh3QJJfLZ8+enZ+f7+Pjk5OTw4e3CRBC6qSRNcrhDmjKycl5+PAhO6jmA5o2bNjw1VdfWVhYVFRU8GGTQISQmml2jXK4LVGOHDnS1tbGDrJ9GhwcbGtrO0D35V76PHr0qLe39wDdBSHEZ1pSoxxuS5SMjIzW1lZ2cIAOaGptbRWJRLW1tWvWrPnuu+9UeGWEkAbRthrlqGFLlNDQ0AMHDjg5ORUVFWniC1cIIZXQ2hrlPPeAJisrK39/f2W2RNm1a1dERISRkVFZWdnrr7+u0sgIIU2i/TXK6X1LFB8fn/5vI3L16lUHB4e2tra9e/cuXLhwwCIjhDTAIKrRp7Gv8B84cODKlSvsSP+3ROnq6nJzc5NKpR988MGPP/6olrwIIf4apDXKUeCAptWrV8fHx0+cOFEqleJLnwihwV6jnH4e0JSbm+vn5zdkyJCCggIXFxdyeRFCfIE1+qw/b4kiFAq9vLwoinJzc/P09Lx79+7333//ySefkE6KEOIFrNEXun79OndAE7sliq6ubk9PD770iRB6GtZo37gtUS5evBgWFrZixQo86RMhxMEafQmdnZ1CoZB0CoQQv2CNIoSQUsifaYwQQhoNaxQhhJSCNYoQQkrBGkUIIaX8H4KYl4U0UqYtAAAAvnpUWHRyZGtpdFBLTCByZGtpdCAyMDIwLjA5LjEAAHice79v7T0GIOABYkYGCOAEYg4gbmBkY0gAiTNDaCYmdJrdwQJIMwMVOGSAGMyMeBjcQAsYmRiYmBlYWBNY2TKY2NgT2DkymDhYEkQYWVlY2dg5xGFOYOCcsVLQofMY9z4QJ/p/psP1S8/sQOxn0bNBbLD48bbjIDVgcYPgxfZnz/gsAbErH560Oxcy2x7E/qbYsp+BwQGZvR9JDZgtBgCP0DC2cKORFAAAANF6VFh0TU9MIHJka2l0IDIwMjAuMDkuMQAAeJylkksOAyEIhveeggtIUMTHutNV01l00Tt03/unyhhnmrRJY4kx/CBfwGCg2W25PJ4wzC/GABSADEAfTykF7p6IjL5Hn2INgyV0HKXFCGuW4ATfEMejFMYo0W+1QSjPUQJmkd7LPCWiZ+m9zE/kDv+itTtl/Z1CyCVw5zme68U6JPa81SZ562WGYv+g7BPZ+YlcXTu9m6geD1G9MERlxyGql/RxE0nXu2eyrnvPFC3TzBXgvC7mBQwNe/pcM9inAAAAjnpUWHRTTUlMRVMgcmRraXQgMjAyMC4wOS4xAAB4nE2NOwpCMRREt2KpcDPcbz6kTOMexCrYuoK3eKOQh+UMZ86MMQamzMf7/pxTLsdVoSVapcQQy0HdkCOrEMODqxbqjhreyhfZVYZaKO1NErCptlWUFRnWXJdUoGJh1GW/ONy5tleS8selEzxN6ae6HR+CtyVNWhIbKQAAAABJRU5ErkJggg==\n"
},
"metadata": {}
}
]
},
{
"metadata": {},
"cell_type": "markdown",
"source": "Get the environments"
},
{
"metadata": {
"trusted": true
},
"cell_type": "code",
"source": "bond = m.GetBondWithIdx(3)\nbeginAtomIdx = bond.GetBeginAtomIdx()\nendAtomIdx = bond.GetEndAtomIdx()\nradius = 1\nenv1 = Chem.FindAtomEnvironmentOfRadiusN(fm,rootedAtAtom=beginAtomIdx,radius=1)\nenv2 = Chem.FindAtomEnvironmentOfRadiusN(fm,rootedAtAtom=endAtomIdx,radius=1)",
"execution_count": 25,
"outputs": []
},
{
"metadata": {
"trusted": true
},
"cell_type": "code",
"source": "submol1 = Chem.PathToSubmol(fm,env1)\nsubmol2 = Chem.PathToSubmol(fm,env2)\nprint(Chem.MolToSmiles(submol1),Chem.MolToSmiles(submol2))",
"execution_count": 26,
"outputs": [
{
"output_type": "stream",
"text": "CC c[nH]c\n",
"name": "stdout"
}
]
},
{
"metadata": {
"trusted": true
},
"cell_type": "code",
"source": "",
"execution_count": null,
"outputs": []
}
],
"metadata": {
"kernelspec": {
"name": "python3",
"display_name": "Python 3",
"language": "python"
},
"language_info": {
"name": "python",
"version": "3.8.2",
"mimetype": "text/x-python",
"codemirror_mode": {
"name": "ipython",
"version": 3
},
"pygments_lexer": "ipython3",
"nbconvert_exporter": "python",
"file_extension": ".py"
},
"gist": {
"id": "",
"data": {
"description": "github3669.ipynb",
"public": true
}
}
},
"nbformat": 4,
"nbformat_minor": 4
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment