Skip to content

Instantly share code, notes, and snippets.

@greglandrum
Created November 14, 2019 07:33
Show Gist options
  • Save greglandrum/5d45b56afe75603b955103cdd0d8e038 to your computer and use it in GitHub Desktop.
Save greglandrum/5d45b56afe75603b955103cdd0d8e038 to your computer and use it in GitHub Desktop.
HellerSubstructureHighlighting.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\nimport rdkit\nrdkit.__version__",
"execution_count": 2,
"outputs": [
{
"output_type": "execute_result",
"execution_count": 2,
"data": {
"text/plain": "'2019.09.1'"
},
"metadata": {}
}
]
},
{
"metadata": {
"trusted": true
},
"cell_type": "code",
"source": "m = Chem.MolFromSmiles('c1cc(C(=O)O)c(OC(=O)C)cc1')\nm",
"execution_count": 5,
"outputs": [
{
"output_type": "execute_result",
"execution_count": 5,
"data": {
"text/plain": "<rdkit.Chem.rdchem.Mol at 0x7eff54f83940>",
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAcIAAACWCAIAAADCEh9HAAAABmJLR0QA/wD/AP+gvaeTAAAeKUlEQVR4nO3deVQT594H8G/YV1HZBBEBW1GU4lVxQ7GWqKAEERq3F9Sj9fVctWit1taeFnpt+2pbrmC97e0mlYtWDW4EEaUoGoRbi1APKlosKBVFVBbZwpLM+8dYZVWWhMmE3+f4R3kIk2+O+O3MPDPPCBiGASGEkO7S4ToAIYTwG9UoIYT0CNUoUYPTp+HtDVNTmJtDKER6+rNvTZ2Kjz9u8WKBAGlpvZuPEFWiGiWqdugQAgIQGIiCAly/Dm9vCIVISeE6FiHqosd1AKJdmpqwfj02bcLGjU9GPvwQxcVYvx7XrnGajBB1ob1RolKXL6OkBKGhLQaXLkVeHoqKOMpEiHpRjRKVevgQAAYPbjFobw8ADx48+fL//g9WVs/+EMJzdFBPVIqtxbt3MXz4s8G7dwHA2vrJl2vXIizs2XeHDOm1dISoA+2NEpXy8ICtLeLiWgzGxcHVFY6OT77s1w8ODs/+EMJztDdKVEpPDzt3YsUKDByIJUugUGDPHsTE4PhxrpMRoi5Uo0TVFi9GXR2++gpbt0JHBxMmPLmMlBAtJaB76omKZWdjzhyYmeHCBdjacp2GELWjc6NEpc6exYwZuH8fQ4fC2JjrNIT0BqpRojpHjmDOHDx+jP/5H5w8iX79uA5ESG+gGiUq8uWXEIshlyMsDLGxMDDgOhAhvYRqlPQYwyAiAmFhYBiEhyM6Gjr0e0X6EPp154OOFkxSKpGSgsZGLrMpFPjf/8VHH0FPD999h4gILsMQwgWqUY33nAWTLlzArFmwtcWCBYiNRU1Nb2erqcG8efj+e5iaIiEBK1f2dgBCNABd8KTZmpowZAhWrmyxRufq1ZDJcO0akpOxaROuXn0ybm6OOXMQHAw/P5iZqT1bWRlEImRkYOBAJCZi8mS1vyMhGolqVLNduoTx43H9Olxdnw1euICpU3H79pPbKwsLkZAAiQQZGWD/No2MIBRCJML8+c/uZFet27fh64vr1+HkhOTkFvEI6WOoRjXbqVPw9UVVVYu9y8JCuLggKwvjxrV4cVERjh5FYiLS0tDUBAC6upg0CWIxFi7EoEGqCpWbmzt8+XLD7GyMGYOkJNjZqWrLhPARnRvVbE8XTGqu1YJJTzk6Yv16pKTg9m3s3o3XXgOACxewYQMcHPDaa4UxMcXFxT1MdO7cOW9v76CKioaAAKSlUYcSQjWq2TqzYFJb9vZYuxapqSgtxd698PeHnh7Onv337t0ODg6jRo2KiIj4/fffuxHn+PHjfn5+FRUVhh4eyoMHYWHRjY0Qom0YouH272eMjJidO5n795m7d5mPP2YMDZkjRxi5vAsbqahQxsUtCg42bnaD5tixYz/55JO8vLxObuOHH37Q09MDsG7dOoVC0Z3PQog2ohrlg6QkxsuLMTZmTE2ZGTOYc+eYDz9kTEwYf39m717m8ePOb6m2tjYhISE0NLRfszs1XVxcwsLCZDKZUqns6Ae3b98OQCAQhIeHq+ATEaJFaIqJnxYtwsGDT/7bxAR+fggOxty5nb+Nvb6+XiaTSaXSAwcOlJaWsoNOTk4BAQFisdjLy0sgELCDCoVi7dq133zzja6u7tdff71q1SpVfxhC+I1qlLdu38axYy2uczI0xLRp8PfHokWdX6GusbHx7Nmzhw8fPnbs2NM+HTJkSHBwcGRkZGNjY0hISHx8vImJycGDB/39/dX0aQjhL6pR/rtzB0lJkEqRnNz6Oiex+Mnj5DpBqVRmZGQkJibGx8f/8ccfU6ZMSUxMnDdvnkwmGzBggFQq9fLyUuOnIIS3qEa1yP37OHYMhw8jLe3JjfY6OvDyQnBwRXBw/04/9YhhmIsXL965c2fbtm2XL1+2t7dPTk52d3dXY3JC+IxqVBuVl0MqRWIikpJQUwMTE1uBwGroULFYvGjRohEjRrxwA3l5eb6+vkVFRW5ubsnJyUPo4Z2EdIxqVKtVVeHEiaKMDLc9e2r+WrjEw8MjKCgoKCho9OjR7f6QUqkcM2ZMbm6ut7f38ePH+/fv34uJCeEfqtE+QS6Xp6SkSCQSqVRaUVHBDjo7O4tEolbz8gB++eWXZcuW1dfX5+XlGRkZcRSZEN6gGu1bFApFZmamRCI5dOhQSUkJO+jo6BgYGCgSiV599VU9Pb38/Pzhw4fb2dndbXUTKiGkPVSjfZRCoTh37hx7ndPTurS3tw8MDAwNDfX393/06NGff/7p0OmJKUL6LKrRvk6pVObk5Eil0p9++om90f6DDz64ePHiqVOn4uPjg4ODuQ5IiKajpUn6Oh0dnXHjxkVERNy4cSMrK+u9995buHDhxIkTAVy8eJHrdITwgB7XAYgGGTdu3Lhx4wCwNfrLL79wnYgQHqCDetKOR48eWVtbm5iYVFRUsKs6EUI6Qgf1pB2WlpbOzs41NTV5eXlcZyFE01GNkvbRcT0hnUQ1StpHNappIiIidu3axXUK0g467UXaN2HCBNBkvWZQKBRr1qz59ttv9fX1AwICnJycuE5EWqApJtI+uVxuYWGhUCjKy8vNzc25jtN31dfXh4aGSiQSQ0PDuLi4119/netEpDU6qCftMzIyeuWVVxQKRXZ2NtdZ+q6KiopZs2ZJJJIBAwakpKRQh2omqlHSITo9yq179+69+uqr58+ft7e3T0tLmzZtGteJSPuoRkmH6PQoh/Ly8iZNmnT58uWRI0dmZma+8sorXCciHaIaJR2ivVGuXLx4cfr06UVFRRMnTjx//ryjoyPXicjzUI2SDg0fPnzAgAF37twpLi7mOksfkpiYOGPGjAcPHohEojNnzlhZWXGdiLwA1SjpkEAg8PT0BB3X96LY2NigoKDa2tply5YdOXLExMSE60TkxahGyfPQ6dHeFB0dvXz58sbGxi1btsTExNBqBnxBf0/keej0aO9gGGbz5s2RkZG6urq7du1as2YN14lIF9DeKHketkazsrIUCgXXWXrm9Gl4e8PUFObmEAqRns51oGcaGhqWLFkSGRlpaGi4f/9+6lDeoRolz2Ntbe3k5FRVVcXvpZ4OHUJAAAIDUVCA69fh7Q2hECkpXMcCgOrqapFIdODAATMzM6lUumDBAq4TkS6jGiUvwPuV8JuasH49Nm3Cxo2wtcXgwfjwQyxbhvXruU6GkpISb2/v06dPDxo0SCaTzZw5k+tEpDuoRskL8P706OXLKClBaGiLwaVLkZeHoiKOMgFAQUHBtGnTcnJyhg0bJpPJxowZw2EY0hNUo+QFeD9Z//AhAAwe3GLQ3h4Azp3DggU4eBBVVb0cKisra/LkyTdv3vT09MzMzHzppZd6OQBRIapR8gJjx47V19fPzc2tqanhOku3sJev//UQ6SfYL9PTIZFg0SLY2EAkwrfforS0FxKlpqb6+PiUlpYKhcLU1FRra+teeFOiPlSj5AWMjY3d3d15vNSThwdsbREX12IwLg6urnj/ffzzn5g6FQ0NSEzE6tVwcMDs2Wrt03379vn5+T1+/DgkJCQpKYkWIdQGDCEv8ve//x3A559/znWQLvrxRyYri2EYZv9+xsiI2bmTuX+fuXuX+fhjxtCQSU5+9soHD5i9exl/f0ZfnwEYgNHRYby8mO3bmZs3VZgoKipKR0cHQFhYmEKhUOGWCYeoRsmLxcTEABCLxVwH6YpPP2UEAsbGhikrYxiGSUpivLwYY2PG1JSZMYM5d679nyore9KnhoZP+hRg3NyY8HDm+vWexFEqlVu2bAEgEAg+++yznmyKaBqqUfJi165dA+Do6Mh1kM5RKplNmxiAEQiYyMhubqSyktm3jwkOZkxMnvbphtdf37Zt29WrV7u6scbGxhUrVgAwMDDYv39/NyMRTUUPESEvxjDMwIEDKyoqiouL7dk5bo3V0IBly3DgAAwMEBuLhQt7usG6Ovz8MySSyrNnbUpLGxoaALi4uPj7+4vFYi8vL4FA8PwN1NTUiMXikydPmpmZSSQSX1/fnkYimobrHif8IBQKARw7dozrIM9VVcXMns0AjJkZc+qUarfdIJcnJiauWLHC0tLy6T8fZ2fnTZs2ZWZmKpXKdn/q4cOHkydPBmBra5vFnqglWodqlHTK1q1bAWzdupXrIB27d48ZO5YBmEGDmOxs9b1PU1OTTCYLCwtrvmNubW0dGhqakJDQ0NDw9JWFhYXDhw9n2/bGjRvqi0S4RQf1pFMSEhLmzZvn4+Pz888/c52lPQUF8PVFfj5cXHDqFHrlanalUpmTkyOVSvft23fz5k120NLScs6cOWKx2MHBQSQSFRcXu7u7Jycna/rJENIDVKOkU0pLS21tbc3NzSsqKthLdjRIVhbmzkVpKTw9kZgIG5tefn+GYbKysg4fPnz48OGnfWpmZlZdXS0UCkNCQiQSyapVq+bNm9fLwUjvoBolneXk5HT79u2rV6+6ublxnaWZM2cwfz4eP4aPD44eBadXszc0NJw4ceL333+XSqWTJk2qra2Nior64osv3n///fXr10dFRXGYjagPLdtMOout0bfeeuuDDz6YMmWKRuyT7tuHFSvQ0ICQEOzZA319DrM0NTVZWVnV1tZWVFRs2bJFoVDo6upCCxYlIC+iAf8SCB/8+OOP6enplpaWp0+fnjZtmq2t7dKlS6VSaWNjI1eRHv373wgNRUMDNm9GbCy3HQpAT09vxIgRCoUiKysLANuhADw9PXV0dLKzs+vr6zkNSNSFapS8WHR09IoVKxQKxcyZM99++21nZ+eHDx/+5z//CQgIsLOzW7lyZVJSUm92BMMwERERHu+8UztiBHbswGef4UUXb/YOdsez1aKCFhYWrq6u9fX1ubm5HOUiasbpdQJE0ymVyk2bNgEQCARffPHF0/ErV66Eh4c3P0lqYmLi7++/d+/eqqoqtUZqaGgICQkBYGBgcPinn9T6Xl0VGxsLICgoqNX48uXLAezevZuTVETdqEZJh+rr6xcvXswW1k8dFNYff/wRFRXV/GYeY2Njtk8rKipUHqm6utrPzw+AqanpyZMnVb79Hrpx4wYAe3v7VuNfffUVgKVLl3KSiqgb1ShpX1VV1ezZswGYmZmd6sQdQYWFha361NDQUCgURkVF3b9/XyWRHj16xN4RZGlpmZmZqZJtqpZSqWTvcbpz507z8UuXLgFwdXXlKhhRK6pR0o6SkpKxY8cCGDRoUHYX7wgqKir65ptv/P39nz5mXVdX18vLKyoqqri4uNuRCgsLXV1dofF3BM2aNQvA4cOHmw82NDSYmJgIBIIydrkpol2oRklrBQUFL7/8MgAXF5f8/Pxub+fhw4d79+719/c3MDBg+1RHR8fLy2v79u1d3Wxubu7gwYMBuLu7t9rR0zQffPABgC1btrQa9/LyAtCZ/XrCO1SjpIWsrCwbGxsA48ePV9XBeFlZ2d69e8Visamp6dMpKTc3t/Dw8OudWMTz7NmzFhYWAGbMmFFZWamSSOqTmJgI4NVXX201vnHjRgD/+Mc/OElF1IpqlDyTmprar18/AD4+Ps8vrF9++WXp0qXHjx+vq6vr/PZramoSEhJCQ0ObPzmD7dOOVj86cuSIkZERO/3dpffiSmlpKXtCuampqfn4gQMHAPj7+3MVjKgP1Sh5Ij4+ni2skJCQ5ssUtWvDhg2t5uW7tJ9YV1fH9im7m8lydnYOCwuTyWRPF53bvXs3e6/UunXrePTIDWdnZwC5ubnNBwsLCwFYW1tzlYqoD9UoYRiGiY6O7tIzgvLz8z/99NNx48Y9LUFjY+PAwMDY2Njy8vLOv69cLj9x4kTbRTzffvvtVatWsderhoeHd/+DcWHRokUAvv/++1bjtra2AAoKCjhJRdSHarSvUyqV4eHhbGFt3769qz9++/btqKgooVDYdl7+3r17nd9O20U8raysdHV1Y2JiuhqJczt37gSwatWqVuP+/v4AOroCl/AX1Wif1tjYuHLlSgB6eno//PBDTzb14MEDdl5e/69729l5+aioqC7NrSsUCplMxl42NH369J5E4sqFCxcAeHh4tBrftm0bgHfefpuTVER9qEb7rurq6jlz5rB3BCUlJalqs48ePWL71NDQsNU8Uuev97x69SqAoUOHqipVb6qrqzMwMNDV1W11X2zJmTN1L72knDqVq2BETWi90T6qrKxMJBJlZGRYWlpKpVL27iDVqq2tTU1NlUgkR48era6uZgfd3NzEYvGCBQuev2gp89dD9O7evWtnZ6fybOq2fNYs96qq0MhImylTno1WVmLgQBgY4PFjztejIqrEdY8TDvTyHUG1tbXsvDx7NRXLxcWl1bx8Kz4+PgCOHz+u7nhqsWYNAzBtn0c/YgQDMJcucZGJqAstlNfnXLlyZdq0aTdu3Bg9erRMJmOfuaZWxsbGIpEoNjb2/v37CQkJy5cvHzhwYEFBwa5du6ZNm+bu7q5UKtv+1MSJE9Fm0TnemDABANou1cyO8/RDkQ5QjfYtaWlpU6dOvXPnzowZM9LT09k7LHuNkZGRSCSKiYkpLS1l5+Xt7OyGDRvW7kL6/F40fuJEoL26ZMd5+qFIB+jcaB9y7NixxYsXy+Xy+fPn79+/n73YnlsKhaK8vNzKyqrtt0pKSuzs7CwsLMrKyjTigSVdwjCwtER5OYqL0fyZoFlZ8PTEyJG4do27cETF+PbbSbrrX//6V3BwsFwuX7du3dMbljinq6vbbocCGDRokKOjY2VlJbuIJ88IBBg/Hmiz4+nhAWNjXL+OigpOchF1oBrtE3bs2LFu3TqGYcLDw7/88ku+7Nzx+/Rou8fv+voYMwYMg6wsTkIRdeDHPyfSbQqFYvXq1e+++66uru4333wTERHBdaIu4Pfp0Y5mk+j0qNahByxrs/r6+pCQkPj4eBMTk4MHD7I3I/KINuyN/vorFAr89ZRQgCbrtRBNMWmt8vLygICA9PT0AQMGSKVSdtlgfqmtrbWwsBAIBBUVFSYmJlzH6TpnZ9y6hStXMGrUs8GCAgwbBhsb3L/PXTKiSnRQr53Ky8unTp2anp4+dOjQjIwMPnYoABMTk9GjRzc2Nubk5HCdpVva3fF0cYGNDUpLUVTESSiiclSj2unPP/80MjJydHSUyWQjRozgOk738fv0aEenQdlJfDqu1xZUo9opIyMjOzt7ypQpQ4YM4TpLj2jD6VGaZdJ2VKPaiW2f7OxsroP0FLs3ytcaHTsW+vrIzUVNTYtxmmXSLlSj2snd3d3U1DQ/P//Ro0dcZ+kRNze3fv363bp16z4fJ2SMjTF6NBQKtDq3O2ECdHRadyvhLapR7aSnpzdmzBiGYbLaXOa9Z8+eN954o66ujpNgXaWjo8M+quTXX3/lOku3tHtcP3Agystx6RIniYjKUY1qLfa4vu3kzK5du3744QcezX3z+/RoQADeegttl3NttmYg4TuqUa3VUfvwrpX4PVnv54d//hPNF28mWodqVGt11D68a6VJkyYBuHjxYrvLkvLJ6dPw9oapKczNIRQiPZ3rQEQ1qEa1lpOTk62t7YMHD9gnpD/Fu71ROzs7BweHioqK/Px8rrP0wKFDCAhAYCAKCnD9Ory9IRQiJYXrWEQFqEa1WbtXC7Fz34WFhTya++Zd9bfW1IT167FpEzZuhK0tBg/Ghx9i2TKsX891MqICVKParN324ePcN+9ORLR2+TJKShAa2mJw6VLk5dEtoVqAalSbddQ+HU3iayze740+fAgArR7Zwq6K/+ABB3mISlGNarMJEybo6OhkZ2c3NDS0GgevWmn8+PF6enqXL1/my+WurbEr/N+922KQ/dLamoM8RKWoRrWZhYXF8OHD5XJ5bm5u83HezX2bmpq6ubk1Njb+9ttvXGfpFg8P2NoiLq7FYFwcXF3h6MhRJqIyVKNart3DYT7OffP79KieHnbuxOefIyoKpaW4dw+ffIKYGERHc52MqADVqJZ7/tWjPDqu5/3p0cWLceQI4uPh5ISXX0ZqKk6fxuzZXMciKkA1quWefy8Tj3bueF+jAPz8kJ6O2lpUV+PMGXh7cx2IqAbVqJZ75ZVXjI2Nb9y4UV5e3nycd600atQoc3PzgoKC0tJSrrMQ0gLVqJbT19f/29/+1napJ97Nfevo6IwdOxZA2zWrCOEW1aj2a/f0qKmp6ciRI/k19827PWjSR1CNaj+tOT3Ku2kx0kdQjWo/fq2Y9/jx4+3bt7d7QSsb+Ndff6WnghONQjWq/ZydnW1sbEpLS2/dutV8XAN37kpKSqZPn/7ee+999NFHbb/r4OAwePDgsrKymzdv9n42QjpCNdoneHp6os3xu6bNfRcUFHh7e//222/Dhg0LbbWKBwBALpfr6+u7urqWlZX1fjxCOkI12ie0e/yuq6urOXPfWVlZkydPzs/P9/T0zMzMfOmll1q9oLy8fNasWbdu3dLX12f/r0CIhqAa7RM6On7XkNOjqampPj4+paWlQqEwNTXVus1qHffu3XvttddkMpmdnV1cXJyODv3eEg1Cv459wsSJEwUCwaVLlxobG5uPa8Lp0X379vn5+T1+/DgkJCQpKcnc3LzVC/Ly8iZNmvTbb7+NHDnyv//9r4eHByc5CekI1Wif0L9//5dfflkul1+5cqX5OOdz39HR0UuXLm1sbAwLC9u7d6++vn6rF1y8eHH69OlFRUUTJ048f/68I62HRDQP1Whf0e7xO4dz3wzDvPvuuxs2bGAY5rPPPouOjm57qJ6SkiIUCh88eCASic6cOWPFrtpJiIahGu0rNGqpp6ampjfeeGPHjh0GBgZxcXGbN29u+5rY2Ni5c+dWVVUtW7bsyJEjJiYmvZmQkC5gSN/AFqibm1ur8RMnTnz99de3b9/utSTV1dV+fn4AzMzMTp482e5roqKiBAIBgC1btiiVyl7LRkg3CBi6IaRvaGxstLCwkMvlZWVl/fv35ypGWVmZv79/Zmamra3tiRMn2IfrNccwzObNmyMjIwUCwRdffLFx40ZOchLSeXRQ31fo6+t7eHgwDHPp0iWuMty6dWvKlCmZmZnOzs7nz59v26ENDQ1LliyJjIw0NDT86aefqEMJL1CN9iHcXt505cqVqVOn3rhxw93dPT09ffjw4a1eUF1dLRKJDhw4YGZmJpVKFy5cyElOQrpKj+sApPdwuKRTWlpaYGBgZWXljBkzjh071q9fv1YvKCkpmTNnTk5OzqBBg06ePDlmzJjeD0lI99DeaB/C1T1LR48e9fPzq6ysDAoKSkpKatuh7N30OTk5w4YNk8lk1KGEX6hG+xAXFxcrK6uSkhKZTNZrb7p79+7XX39dLpe/+eabEonEyMio1QteeDc9IRqOarQPEQgEa9euHTlypLe396hRoyIiItQ63cQwTERExJtvvskwTHh4+K5du9peYP/Cu+kJ4QFOL7ciHFi1apWFhcXTX4ARI0Zs3br10qVLqn2XpqamVatWAdDV1f3uu+/afU1cXBx792dISEhDQ4NqAxDSa6hG+6KmpiaZTBYWFmZra/u0Tx0dHcPCwmQymUKh6Plb1NfXz5w509TUNDExsd0XREVFsTunYWFhKnlHQrhCl9/3aQqFIjMzUyKRxMfH3717lx20trb29fUVi8W+vr5t1wrpvMrKyj/++INd0rQ5hmHee++9HTt2CASCHTt2tHsnKCE8QjVKAECpVObk5Eil0v379+fn57ODAwcOnDt3rlgsnj17toGBgUreqKmpafXq1Xv27DEwMIiJiVmyZIlKNksIh6hGSWtXr16VSCSHDh3Ky8tjR/r37z9z5kx/f/+goCAzM7Nub7mmpkYsFp88edLMzEwikfj6+qooMiFcoholHbp69WpiYqJUKr1w4QI7Ymxs7OPjIxaLAwMD217++XyPHj0SiUTPuZueEJ6iGiUvduvWrePHj0skkoyMDPYXxsjISCgUikSiwMBAGxubzmxh9uzZv//+u7Ozc3Jycts7QQnhL6pR0gV//vnnkSNHEhMT09LSmpqaAOjq6k6aNEksFi9YsMDOzq7dn7py5Yqvr29xcbG7u3tycrK9vX3vpiZEvahGSXc8fPgwKSlJIpGcOnWKfb6Tjo7O5MmTxWJxcHCwg4PD01eePXt2/vz5z7mbnhC+oxolPVJeXi6VSiUSSUpKSn19PTvo5uYmFosXL1587dq1JUuWyOXyoKCgffv2tb0TlBAtQDVKVKOysvLEiROHDx9OTk6ura1lB3V1dRUKxYYNGyIjI+mpyERbUY0SFaurq/v5558lEklCQsKaNWv69+//zjvvcB2KEDWiGiXqUl9fb2BgwD5SiRAtRjVKCCE9QqerCCGkR6hGCSGkR6hGCSGkR6hGCSGkR/4fopAJphO25WkAAAAASUVORK5CYII=\n"
},
"metadata": {}
}
]
},
{
"metadata": {
"trusted": true
},
"cell_type": "code",
"source": "print(m.GetSubstructMatches(Chem.MolFromSmarts('C(=O)O')))\nm",
"execution_count": 6,
"outputs": [
{
"output_type": "stream",
"text": "((3, 4, 5), (8, 9, 7))\n",
"name": "stdout"
},
{
"output_type": "execute_result",
"execution_count": 6,
"data": {
"text/plain": "<rdkit.Chem.rdchem.Mol at 0x7eff54f83940>",
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAcIAAACWCAIAAADCEh9HAAAABmJLR0QA/wD/AP+gvaeTAAAgAElEQVR4nO3dd1xUZ9o38N8UBpCu0qQpKDCDYqIIcZVYiIqIYCMSZHFNsqu7mxDdrDGbLdHPvu8b3TVZ8UnWlCfGGkkmWAALi6sIRgQsaCgK0puASBWY/v5xCGUYmpSZYa7vX3jmzJlrEvx57nM3lkKhACGEkOfFVncBhBCi3ShGCSFkWChGCSFkWChGCSFkWChGCSFkWLjqLoDoBpkMRUUoKUFNDaRSGBvD3h7Tp8PCQt2VETJcLBrwREaXTIbUVKSkgMWCWIzO3zcuFwBsbREQABsbNRZIyDBRjJLR1NiIY8fQ0gKxWPUJLBY4HCxahIULx7YyQkYMxSgZNY2N+PJLtLdDLh/gTD09eHvjlVfGpCxCRhh1MZHRIZfj+HG0tQ2coQAkEqSnIzt79MsiZORRjJLRkZaG5mYMvq0jkSA+HhLJaNZEyKigGCWjQC5HcnKfz0P7IpPhzp3RKYiQUUQDnsgoKC3tuA+9exc3buDpUxgZYdYsLFnS0UH/6ad45RW4u3ecn5yMykqEhuLOHfj4qK1sQp4LxSgZBbm5EIk6xjkFBWHaNNTVIS4OdXUIDe3vjXV1kMvBpkYS0Sb0+0pGWkUFbt+GWIykJKxaBXd36OtjyhSEhiI/HyUl/b2Xw0Fz81gVSsjIoBglI6qwEMeOQSZDZSUkEvD5XS+ZmcHODkVFA1xBJhvVAgkZcdSoJyPn3j3ExnaMcGptxYQJys1zY2O0tnb8fOYMOJyOnyUSuLhAoYBMhgkTxrBiQkYAxSgZIWlpuHSp64+GhmhtVX7Q2dICK6uOn1eswPTpXe+tqwOLBR4PBgZjVTEhI4NiVAu1tuL+feTmoq4OYjG4XEycCFdXeHrC3FwN9SgUuHwZN270OGhnBy4Xubnw8Og40tiIigr4+XX8ccIEmJp2/KyvDwAsVlffPSHag2JUq0iluHoVGRkAukaqSySoqMDjx0hJAZ+PlSthaDh2JclkOHsWWVnKx3k8LF6M8+fB5WLqVDx9ithYuLjAyanPS3E4eOmlUS2WkNFAMao9Wlpw5AiamlRP9WF6ZnJzUVCAzZu72s6jSiTCd9/12Ws0fTp+/BEXL6K5GRMmYNYsLF3a56WYW9GxKZuQEUVLk2iJ9nZ88QWamgaeos48Ydy6ddSX8nz2DCdPoqpK9atFRYiOhkiEmTOxYcMAl2KzYWqKbds6WveEaBUa8KQlzpxBc/OglvlQKCAS4fBhtLWNYj319Th8uM8Mzc3FyZMQiSAQYM2aAS7F5cLUFK+/ThlKtBTdjWqDkhKcPDnkZTvYbEydCldXzJwJI6ORrKeyEt9+i2fPVL+ano6LF6FQwMcH/v5gsfqrkM2GhwcCAsDjjWSFhIwhilFtcPw4CguVD/Y1Xb03Fgv29vDwgEAAE5PhFlNUhO++g0ik4iWFAteuISkJABYvxuLFfV6ExYKxMdzd4e2NyZOHWxIhakUxqvGkUuzdqzy3p/d0dTOzjunq//kPAPD5sLdXvhNk8tTNDXw+Jk58nmJ++gnnzqmeaCSX4/x53L4NNhuBgZgzR/UVjI2xdi0cHfsMfUK0DcWoxqusxLFjPe7+xGJ8/DGCgnoMyTx4EBERsLPDP//ZcbKREaZPh0CA6dO75gt1srSEQICZM4dwM5iWhoQE1UuISiQQCpGXBx4PISGYMUP1FSwtER7eNVyUkHGB7gg0XkuL8pF+pqs7OmLzZuTl4aefUFeHe/dw7x4MDeHqqpyntbW4dg3XrnXkqZsbbG37rIFprV+7pvrVtjacOoXSUhgaIiwMDg6qT7OzQ1gYzfUk4w/FqDZQapv3M12dxcKUKZgyBYsXo6YGOTnIysKTJx15amAAFxe4uoLP7+rS6cxTc3O4uUEggINDj09sa8OlS7h/X3Vtzc04cQLV1TA3R3h4n/e2bm7YsIEa8mRcol9rjWdiotyOHnC6OsPKClZWWLwY1dXIzUVuLqqrkZ2N7GycP48ZMxQCwbMZM4w787ShAWlpSEuDmRlcXaGnh9JSPH4MmayvvUCya2r+dObM6dZWrpVVf631F17A6tW0iigZryhGNZ6VlXKXzoDT1QEoFF13lNbWsLbG4sWor8fDh8jJQWkpsrOb6upszpx5xdl5tZvbGnd3q85BUY2NHfNN+5VcUhJ86lRDe/vfZ87cExjY55IiCxbQlp9kfKMuJm3w7bd49KjHLeGNG7h+HcHBXdPVTUywfv1gR7A3NCA3N76mJjgzU65QANDjcJZOm7bO3X0Nn281iEGm5x48eC0mpk0iWePu/u369YZ6eipOYrEQEAAvr8F9SUK0FcWoNqiowNGjysPvb9/GzZt4+rRjuvqSJYiKgqEhPDzg4QFLyx43pH140tp6IT9fmJ2dUFAgkckAsFms+Q4OIQLBeoHAvo9G+jd37/4mLk4ql//e2/vgypVslZ/C4WDt2q77ZULGL4pRLfHDD3j4EFJpnyfU1uKrr7o247SxAZ8PPn+Qi308bWuLz8sTZmcnFhaKfv4UgaVliIfHazNnunXrONp3/fr7ly8D2LVw4d6+WusGBggN7W8xJ0LGEYpRLSGR4Isv0NDQ3x4bUikKCpCTg4cP0d7ecdDCAq6u8PBQ7n/vQ2N7e3xe3unc3EuPHrX+fP/7oq1tlL//Lxwc3rpw4fNbtzhs9mcBAVv7aq0bG2PTJtjYDPEbEqKtKEa1R1sbTpxAbe3Ak+vlcpSXd3TKdw47NTODuzvc3DB1ap+d5t2eA7RJJJcLC4U5ObEPHza2t9/ZuvWj69eF2dn6XO6Jdes2CASqr2BhgV/+ctQXlyJEk1CMahW5HDduICUFQFf7ncFiqRiWJJejpAQ5OXjwoGvHTROT+hdeuOvsvMjJiTOIQUgiqfRCfv6BmzeTS0osDA3PhYb69tVanzIFmzbRAHuiayhGtZBYjNxcPHiA6uqOTUQmT4arKxwdUVKCvDwUFysvqadQoKqqc3bTlalT/YqLJxoarnJ1DREIlru46Pc9ML6quXnlyZP3Hj+2NTG5FB7uaW2t+jxnZ2zcSAs1ER1EMToetbYiPx85OXj0SMUSpZWVRx492nv//sMnT5gD5gYGga6u6wWCFS4uvYcuXS4sDDx5cpqFxX9++UsHMzPVn+jpieBgGmBPdBPF6LjW1oa8POTkoKCgd99Udk1NfF5eXF7ej6WlzBFDPT2/adNCPDzWuLub/jwEdf7//u/N8vIzGzeu6T6LvzsfH6xYMZj+K0LGJYpR3SCRoKgI2dl48ED5oSpQ3NBw7sEDYU7OjdJS5rfBgMtd6OgY6Or62qxZ/y8lJermzf/r5/eBr6+KK+vr4/33R71+QjQYxaiOkUpRWIjsbDx82GPxPYUCLFZpY+Pp3NzTubk/lpZ2zm76P0uX7kpMDHJzO/faayou6OEx8FZLhIxrFKO6Sibr6I/Kyuq9HUjn7Kb71dWXIyJc/+d/rIyMqnfu7HGSQgF9fbz6Klxcxq5sQjQPxajOUyhQVoacHNy+3XuWlEgq1edyrf/5z5pnz0p27HBU6mKaOBFvvUVPRYmOo65VncdiwdER/v54553ew5WYgVBeU6YASCsv7/Eal4t16yhDCaEYJT8zNkZ4OFSt1eRtZwcgraKi65CeHvz9YWc3ZtURorEoRkk3Dg6IiICBgdIy9T729gDSmRhls6Gnh8BAzJ2rlhoJ0TQUo6Qne3tERsLTE1xu5+ql3nZ2LOB2ZaWUzca0adi2DZ6e6i1TFxw5cuTtt99WdxVkYNTFRPrQ3t4xr7S+HnK565//nF9VdTcl5YWFC9VdmU7Ys2fP7t27ASQnJ/uqHLFLNAZtIkL6YGAAT8/Ou06fK1fyT5xIy86mGB1tMpns7bffPnToEIfD+fTTTylDNR816smgeHt7A0hPT1d3Ic9FIkFpKX76Cffu4dGj3uNkNYdIJAoLCzt06JC+vn50dPS2bdvUXREZGN2NkkHx8fEBkJaWpu5Chqi8HElJKC4Gl9uxkCCbDYkEFhZYuBCenho1YKuhoSE4ODg5Odnc3Dw2NpbuQ7UFPRslgyIWi01NTSUSSX19vWlfGylrFKkUcXHIze1zlWseD2ZmeO01DVlkuqqqKiAgIDMz09bW9uLFi7Nnz1Z3RWSwqFFPBoXH473wwgtyufzWrVvqrmUQxGJ8/TVycvrbKUAkwpMn+OILVFePYWWqFRQU+Pr6ZmZm8vn8mzdvUoZqF4pRMlha83hUoUB0NJ48GWC3FWa/AJEIR4+q92lpRkbG/PnzCwoKvL29k5OTHR0d1VgMeQ4Uo2SwtObx6N27qKiAVDrY555tbfj8c9y712PJq7GSmJjo5+dXW1u7evXqq1evTu62CSvRFvRslAzWo0ePZsyYYWtrW1lZqe5a+iaX4+OP0dr6PO/lcuHsDFdX8Pljs6PUsWPH3nzzTYlEsnnz5q+++kpP1UxcovkoRslgKRQKS0vLurq60tJSBwcHdZfTh4ICCIUq7ivv3sWNG3j6FEZGmDULS5ag7+2nwGLB3h4eHvDwgLHxKFUaFRW1Y8cOhUIRGRl54MABliaNGSBDQgOeyGCxWKx58+ZdunQpPT1dc2M0P7/38v5ITUVKCoKCMG0a6uoQF4e6OoSG9nkRZvHAsjIkJHTkqUAAE5ORqlGhULz33nv79+9nsVj79+9/9913R+rKRC3o2SgZAi14PFpVpbzRtFiMpCSsWgV3d+jrY8oUhIYiPx8lJcrv7d0yY/L00iV88gn+/W8kJaGubpgFisXisLCw/fv383i8U6dOUYaOA3Q3SoZAC2K0d3O+shISCbrvx2dmBjs7FBXByanHmf03q2trce0arl2DpSUEAnh4wNJyqNW1tLRs2LAhISHB2Ng4JiZm+fLlQ70C0UAUo2QIfHx8WCzWrVu3pFIpt59ni2pkYKB8pLUVEyYob/5sbIzaWnz2WccDUCYQFYpBde535qmFBVxdIRDAwQFSKUpKUFOD1lbo6WHiRDg5oec8herq6oCAgDt37tjY2Fy4cOHFF18c1jclGkMj/yYQTTVx4kRnZ+eCgoKcnBxPzVwrz9oapaU9mueGhmhthVzeI0lbWsBiobYWSUlISoKVFfh88PmwsRnCZ9XXIy0NaWng8SCRgMuFXA6ZDCwWeDzIZLC0xLJlmDYNQGFhob+/f35+vrOzc0JCwvTp00fqGxO1oxglQ+Pj41NQUJCWlqahMTp9etcIUJkMbDbs7MDlIjcXHh4d5zQ2oqICmzYBwMOHyMpCTQ1qanDtGszN4eYGDw84OAxhuj3TqdU52p8Z1Q+gqgrR0XBzu21vH7B6dU1NjZeX1/nz562srEbuCxP14zBrGhIySBUVFZcuXbKysgoKClJ3LapYWODGDUilEIvx/fcoKwOfDw4H//0vJk+GiQlqaxETAxsbvPwyLCwwYwbmz4eLC3g8NDWhqQkVFbh7F5mZePIEACZOHNbyJTLZlVu3Av74x6eNjX5+fhcvXpw4ceJIfVeiIWjcKBmamzdvzp8/f9asWffv31d3LX24dQtnzuDoUVRUwMgI27bBxAS3b+PmTTx9igkTMGsWli5VMW5ULkdpKXJzkZuLpqaOg8bGcHcvnTlziqMjlz3kkS0n799//dw5sUwW/vLLhy9fpgH24xLFKBkakUhkZmam0Us9FRXhF7/A48ewsEB4OCZNep6L1NQgJwdZWcw96TxLy4KWlkBX19WurgEzZhj12kJVpYNpaTsuXZIrFJE+Pv8KCmKHhsLZ+XmKIZqNYpQMmbe3d0ZGxtWrVxcvXqzuWnrJyoK/PyoqMGUKwsNHYE5nVdWzvDzvrKyc2lrmgKm+/ipX13V8fsCMGRP6uLtUKBR7rl3bk5TEAvYuW/beggUAYG0NWoZ5PKJno2TI7t+/n56ezufzF2rahiJJSVi+HLW1WLIEV66gvR1NTZDJlE8b5MAmhokJb+rU33t7R8ye7Wxh0SwWFzx9+lNNjTA7+5PU1LTycqlcPtXc3KDbIwKpXL41Pv5fqak8DufYunXbvLw6XhCJ4OmpYkgW0XJ0N0qG7Pjx4xEREevWrYuJiVF3Ld2cOYOwMLS3Y906nDzZkVbFxbh6FRUVHavfKxRgsyGTQU8PhoZ4+vQ5PqekoeHsgwfCnJzUsjK5QgGAw2a/ZG8fIhCEzpxpzOOFCIUX8/ONeDxhSMjKGTO63snjYeVKvPDCyHxfojEoRsmQ5eXlubm5TZkypYLZuV4TfPYZIiMhl+OttxAVpTzYvr0dFRVoaoJcDiMj2NrCzAwAGhvx4AGys1FW9hyfWd7UdCY393RubkppqUwuB8BlsycZGlY/e2ZlZLR32bK61taI2bOtjIw63sBiYeFCLF06vK9KNA7FKBmyzqWeysvL7ezs1F4N9uzBnj1gsfC3v+H5HlI1NSE/H3l5ePQIcvlQ313X2hqflyfMzk4sLFzm4lJYX382NPStCxcSCwrOhIaucXfvOtXbGytXPk+FRIPR8HsyZCwWy9XVNTU19S9/+cvu3budlGamjyWZDL/7Hb78EhwODh3Cr3/9nNcxNcXcuZg7F21tyMtDTg4KClQ8VO3DP2/c+HdGxpE1a46vW/dMIrE1Nuaw2d52dokFBekVFV0xymaPzTKmZIzRCk9kyGJjY+/cuWNtbX3kyJGpU6d6eHjs3r07Nzd3jMtobW1tef11fPkljIwQG/v8GdqdoSFmz8Zrr+GPf8SaNRAIMIiRnjwOp1kkSq+osDA0tDc15bDZAHzs7ACklZd3naenB5q/NB5RjJKhOXLkyPr160UikaenZ0hIiLGxcU5Ozp49ewQCwaxZs3bv3j02w/Lr6+tXrFgRePWqdMYMJCQgIGCEP8DAALNnIyQEO3fitdfg6Ym+x4qqSEzAx94eQEZlpazzKYFMxsyvJ+MMNerJEOzbt+/9998HsGvXrr179wJob29PTEyMj48/e/ZsVlZWVlbWnj17pk6dGhQUFBISsmDBgtFY1L2srMzf3z8nJ2fq1KlFcXEz3NxG/CO66OnB1RWurpBKUViI7Gzk5aG9HegaOOVtZwfgVmWlTC7n/Ny7ZWVkNNXcvLih4cGTJx7MTaidHY12Gpeoi4kMikwme/vttw8dOsThcD799NNtvYaRy2Sy1NRUoVAoFAqrqqqYgw4ODmvXrl29evXixYtHamG9nJwcf3//srIyDw+PS5cu2dvbj8hlh0ChQGkpvvsObW2dx1yiogrr6+//9rezrK07D24UCr/Pzj4cHLyFWRNv0SJo4IQFMmzUqCcDE4lEYWFhhw4d0tfXj46O7p2hADgczsKFC6OiosrLy1NSUiIjI+3t7cvKyg4ePLhs2TJbW9uIiIi4uDhJ/5seDyQtLW3RokVlZWWLFi26fv26GjIUAIsFJye88Ub3Zj5zQ5recwQY065P6zz4XMNUieajGCUDaGhoWL58+ffff29ubp6YmLhhw4b+z2ez2UyelpWVZWVlffjhh25ubk+ePDl+/HhQUJC1tXVERIRQKGwd+uadsbGxS5YsefLkyZo1ay5evGhubv6832kkTJqEsDBwOACgUCgnJoDez0w1Z5gtGVEUo6Q/VVVVS5YsSU5OtrW1TUpK8vX1HdLbmU78Bw8eMHkqEAjq6+uPHz/+6quvWlparl69+tixY83NzYO5FNO11dbW9vrrrwuFQkNDw+f6QiPKyQmrVgEAi6WcmAoFgDm2tnoczk81Nc+YBUmfPu3+HICMGxSjpE8FBQW+vr6ZmZl8Pv/mzZuzZ89WeVpxcfETZmnOvjF5mp2dXVBQcODAgQULFrS1tcXHx2/evNnKyorJ08bGxr7evm/fvi1btkil0l27dn399dcatH+JQMD0Mr1oa8vjcLJra1uYxGSxABjq6c20spLJ5XcfP+44n25IxyOKUaJaRkbG/PnzCwoKvL29k5OTHR0d+zrzgw8+sLGxYRryA04PdXZ2fuedd65fv15UVMTkqUgkYvJ00qRJzEWqq6s7z5fJZL/73e/ef/99Dofz+eefM8MDNIi+PiZOBGDA5c6ytpbJ5Xd+7l5jULteF1CMEhUSExP9/Pxqa2tXr1599erVyZMn93OyXC7ncDg//vjj9u3bHR0dFyxY8MknnxQXF/f/EU5OTkyeVldXHz16NDAwkM1mMxexs7Nj8rS4uLh719bWrVtH8kuOlJ+7uVSOHmW6nrqemVKMjkc04IkoO3bs2JtvvimRSDZv3vzVV18NZsH2hoaGxMTEuLi4M2fOtLS0MAcFAkFISMjGjRv53Tc37ltNTc3Zs2djYmKuXr3KdOizWCyFQjFx4sTY2NgFzJKdGigjAxcuADh2797mM2fWCwQ/vPpq54s5tbUen33mZG5evH07ABgaYufOYe1KQjQPxSjpISoqaseOHQqFYteuXR999NFQB8+3tbVdvnxZKBSeO3eu6ed9OJydnQMDAwc/Gr8zlM+ePdva2nr16tWhdm2NqcpKfPUVgAdPnvA//dTBzKx0x47OF+UKhcXevU0iUfUf/2hlbAwAkZGwsFBXsWQ0UKOedFAoFDt37ty+fTuA/fv379279zkmIBkaGjL9RdXV1bGxsb/5zW+srKwKCwsPHjzo6+vb+WC0/3+8zc3NQ0JCjh075uXlJZPJ+ul60gg2Nsy2Tm6TJlkYGpY1NlZ2G3vAZrGy+HyZkZFVZ3O+Z6ufjAMUowQAxGJxWFjY/v37eTzeqVOn3n333WFe0MDAYPXq1V988UVlZSUzGt/W1ra4uJjJUycnp61bt8bFxUml0n4u4uPjAyAtLW2YxYwuNpvZ3Z7FYm2fMeOUkxOnZy+Tg4kJ+9kz0OPR8YtilKClpSUoKCg6OtrY2DguLm7jxo0jeHGl2U27du1ycXEpKyv78ssvg4KCbGxsmNlNYmacUE/aEaPo6mX6m7l5aEmJdWlpj1eZJVmps378ohjVddXV1YsWLUpISLCxsUlOTl6+fPkofRAzu2nv3r35+fmpqak7d+50dnauq6tjZjc5Ozv3vjNlYjQ9PV0+9KWUx1Tn2tXMD0pByYRsZSUzJh9VVYNfyZRoBYpRnVZYWOjr63vnzh1nZ+eUlJQXmRU0RhmLxXrppZf+8Y9/FBQUMLOb5s6dO2/evN6D6m1tbR0cHBobG/Pz88egsOfXGaOdidk9942NYWaG9nZmr2bIZOg2MJaMAxSjuuv27dvz58/Pz8/38vJKTU2dPn362NfAzG66devWyZMnVZ7g7e0NzW/XW1iA2XDJyAjm5hCJoDSti4lXejw6TlGM6qgrV64sXbq0pqbGz8/vypUrVupelX1CH7trdLbrx7acoZsypeMHpcRkKDX2qbN+fKEY1UUnT55cuXJlU1NTeHj4xYsXTUxM1F1Rn7TjbhQDPR6lXqZxjWJU5xw8eDAiIkIsFkdGRh49enQwk5TUyMvLi8vl3rt3r03D10bqXPlUKTEZU6aAzUZNDZjlVuvqOtbPJ+MCxagOUSgUu3fvfueddxQKxb59+6KiothsTf8FMDIyEggEEokkMzNT3bX0q7NRb2sLDgc1Neg+hIvZzE4mQ+eQUrohHUc0/W8RGSlSqfTXv/71nj17eDzeiRMn3nvvPXVXNFja8XjU0JBZ6qkjMeVy9ByEr9zYpxgdRyhGdUVFRUVsbKyRkVF8fHxYWJi6yxkCrXk8qtSuVzl6lGJ0PKIY1RUHDhyora194403li1bpu5ahkZr5jIpjR7tv7O+uBi0KtB4QTGqK+bMmQOgXAuH2nh4eJiamhYWFtbU1Ki7ln4pddYr/ae2tIS+PurrwSwkKBbj5Elo+OwsMjgUo7qCaRpr+hNGVdhsNvNvwK1bt9RdS79sbDp2uJs8GQYGaGxE922mWCzw+fD0ROec1+JinD6thjrJSKMY1RWurq4WFhbl5eVK+3xkZmbu3btXw+NVOx6PcjhgNitlsTo67isre5ywZg3WrUPnhqYyGfLzce/e2FZJRh7FqK5gsVjz5s1DrxtSoVD4pz/96dy5c2qqa1C0o7NeocDPK1VjxQpERsLNbYC3iMVISKCVSrQdxagOURlGWtHY7yxSozdrePSoq9fI2rpj/BOAu3fx2Wf4+9/xySdITITSQlYyGR48GNM6yUijGNUhKpvGTLZmZGRo8mJ09vb2dnZ2T58+ffTokbpr6VtOjor7ytRUJCbCzw/vvYfQUBQW4ocfepwgFiMra8xqJKOBYlSHMIl569YtWbe/7TY2No6Ojo2NjQ8fPlRfaQPTgsej5eXKY5jEYiQlYdUquLtDXx9TpiA0FPn5KCnpcZrSQH2ibShGdYilpeW0adOam5tzc3O7H9eKgZla8Hi096z/ykpIJOi+MaqZGezsUFQ0wBuJVqEY1S0q7+m06PGoRmd97wUKWlsxYYLycWNjtLYO8EaiVej/n25ReU+nFXejXl5eHA4nMzNTJBKpu5Y+mJoqHzE0RGur8hj7lhYoLa6qwQsVksGgGNUtKhNz7ty5XC73p59+alW6S9IkJiYmfD5fLBZr7lJPLi7K95V2duBy0f0RSmMjKiowbVrXETYbzs5jVCEZHRSjumXOnDl6enpZWVnPnj3rPDhhwoSZM2dKJJK7d++qsbYBafrj0VmzlGOUx8PixTh/Hg8fQiRCVRWio+HiAienrnM4HHh6jnGlZGRRjOoWAwODWbNmyWSyO3fudD9Oj0dHwOTJcHBQTtJf/AJ+frh8Gf/4B779FtOm4dVXu15lsWBl1bVWKdFOFKM6R2W7Xisej2r63SiAoKCOafXdzZ2L3/8ef/0r3n0Xy5ej+waoXC7WrBnLAslooBjVOf0MwtfwGJ05c6axsfGjR4/q6urUXUsfzM0REoJB7svC5SI4GJMnj3JNZNRRjOoclfd0fD7f1NS0uLi4WoO3UOdwOHPmzFEoFBkZGequpW8zZmDjRvB4Km5LO7HZ0NPD+vXw8BjDyshooRjVOe7u7ubm5rr8M9EAAAaaSURBVKWlpZXd1h9is9lz584FoNEJpfmPRxkuLnj7bfD54HLB4/V4SU8PXC5cXfHWW3B3V1N9ZIRRjOocFovl5eWFXompFe16LXg8yjA2xvr12L4d/v6YPRtOTnB0hKcn/P3xzjvYuFHFIFOitbgDn0LGHR8fn8uXL6enpwcHB3ce1KgbvW+++YbH423atEnpeGeRCoWCxWKpo7ShMDLCiy/ixRfVXQcZXXQ3qotUJuZLL70EID09Xe1LPX300UdvvPHGli1beq/n5OjoaGtrW1dXV1hYqJbaCOmNYlQXqUxMW1tbe3v7xsbG/Px8dRWmUCh27tz5wQcfANi7d+/06dOVTsjJyWltbQ0ODtYbZG84IaOPYlQXWVlZOTk5NTc3P+i5YLB6H4+KxeKwsLD9+/fzeLxTp0794Q9/UDrh+vXrCxcubGxsFIlEjo6OaimSkN4oRnVUP0s9qSVGW1pagoKCoqOjjY2N4+LiNm7cqHRCbGzs8uXL6+vrg4ODT9NOcESTUIzqqH6Wehr7fvDq6upFixYlJCTY2NgkJycvX75c6YQjR46sX7++ra1ty5YtP/zwg6Gh4RhXSEg/KEZ1lMr2u5eXF5fLvXfvXtsYLiRcWFjo6+t7584dZ2fnlJSUF3v1a+/bt2/Lli1SqXTXrl2HDx/mcml4CdEwCqKTnj17xuVyuVxuS0tL9+Oenp4Abty4MTZlZGRkWFlZAZg3b151dbXSq1Kp9Le//S0ADodz6NChsSmJkKGiu1EdxSyOJ5VKlZbvXLVq1dq1a8emH/zKlSt+fn41NTV+fn7//e9/mTztJBKJwsLCDh06pK+vf+rUqW3bto1BSYQ8D3XnOFGbrVu3Avj444/V8uknTpzg8XgAwsPDxWKx0qv19fUvv/wyAHNz8+TkZLVUSMgg0d2o7lJjv3xUVFRERIRYLI6MjDx69KjSzW9VVdWSJUuSk5NtbW2TkpJ8fX3HvkJCBo9iVHeppV9eoVDs3r17+/btCoVi3759UVFR7J7rHBcUFPj6+mZmZvL5/Js3b86ePXssyyPkObAUSjtrE50hl8stLCyampoeP35sbW09Bp8olUq3bdv29ddf83i8b775JiwsTOmEjIyMVatW1dbWent7nz9/fjKtxUm0Ad2N6i42m80s9TQ2N6TPnj0LCgr6+uuvjYyMzp071ztDL1++7OfnV1tbGxgYePXqVcpQoi0oRnUa064/ePBgUVHRqH7Q06dPly1bdvHixUmTJl2+fNnf31/phOPHjwcEBDQ3N0dERJw+fXqC0hbEhGgyNXdxEbW6e/fuG2+8wfwmCASCDz/8MDs7e8Q/paioyM3NDcC0adMePnzY+4QDBw4wq95FRkbK5fIRL4CQUUUxqutSUlJCQ0NNTEw6/2UVCAR//etf7969OyLXr62ttbW1BTBnzpzHjx8rvSqXy3fu3AmAxWLt379/RD6RkDFGXUwEAEQiUUpKSlxcXHR0dE1NDXPQyckpODg4JCRkwYIFw1kj+c9//nNqaurZs2dNey75LhaLf/WrX506dYrH4x09ejQ0NHRY34EQNaEYJT3IZLLU1FShUCgUCquqqpiD9vb2AQEBgYGBK1eufI4p7QqFQiqVKg0ObWlp2bBhQ0JCgrGxcUxMTO/lSAjRFhSjRDW5XH7jxo34+PiYmJjOVegnTZoUEBAQEhKyYsUKntJmbUNRXV0dEBBw584dGxubCxcu9F6OhBAtQjFKBpadnS0UCqOjox8+fMgcsbCwCAwMXL16dUBAgJGR0ZCuVlRUtGLFivz8fGdn54SEhN5L3BOiXShGyRAweRofH3/79m3myIQJE5YuXRoSErJ27dru/VR9uX37dkBAQE1NjZeX1/nz55WWIyFEG1GMkudRVFQUGxsrFAqZJfUAGBgYvPLKKyEhIUFBQebm5irfdeXKlbVr1zY1Nfn5+Z0+fdqUNhkm4wLFKBmWkpKSs2fPCoXC1NRUZoM8Dofz0ksvhYSEbNy40cbGpvPMmJiY8PDw9vb28PDww4cP0550ZNygGCUjo7a29uLFi0KhMCEhQSKRoFuebtiwISYmZseOHXK5PDIy8l//+pfSciSEaDWKUTLCamtrz549e/r06StXrojFYgBsNlsulzMD7Hvv90mItqMYJaOloaEhMTExLi7u3Llzb7755pw5czZt2qTuoggZeRSjZNS1t7cbGBiouwpCRgvFKCGEDAs96SeEkGGhGCWEkGGhGCWEkGGhGCWEkGH5/7kuvZTlGXrJAAAAAElFTkSuQmCC\n"
},
"metadata": {}
}
]
},
{
"metadata": {
"trusted": true
},
"cell_type": "code",
"source": "m.__sssAtoms = None\nm",
"execution_count": 7,
"outputs": [
{
"output_type": "execute_result",
"execution_count": 7,
"data": {
"text/plain": "<rdkit.Chem.rdchem.Mol at 0x7eff54f83940>",
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAcIAAACWCAIAAADCEh9HAAAABmJLR0QA/wD/AP+gvaeTAAAeKUlEQVR4nO3deVQT594H8G/YV1HZBBEBW1GU4lVxQ7GWqKAEERq3F9Sj9fVctWit1taeFnpt+2pbrmC97e0mlYtWDW4EEaUoGoRbi1APKlosKBVFVBbZwpLM+8dYZVWWhMmE3+f4R3kIk2+O+O3MPDPPCBiGASGEkO7S4ToAIYTwG9UoIYT0CNUoUYPTp+HtDVNTmJtDKER6+rNvTZ2Kjz9u8WKBAGlpvZuPEFWiGiWqdugQAgIQGIiCAly/Dm9vCIVISeE6FiHqosd1AKJdmpqwfj02bcLGjU9GPvwQxcVYvx7XrnGajBB1ob1RolKXL6OkBKGhLQaXLkVeHoqKOMpEiHpRjRKVevgQAAYPbjFobw8ADx48+fL//g9WVs/+EMJzdFBPVIqtxbt3MXz4s8G7dwHA2vrJl2vXIizs2XeHDOm1dISoA+2NEpXy8ICtLeLiWgzGxcHVFY6OT77s1w8ODs/+EMJztDdKVEpPDzt3YsUKDByIJUugUGDPHsTE4PhxrpMRoi5Uo0TVFi9GXR2++gpbt0JHBxMmPLmMlBAtJaB76omKZWdjzhyYmeHCBdjacp2GELWjc6NEpc6exYwZuH8fQ4fC2JjrNIT0BqpRojpHjmDOHDx+jP/5H5w8iX79uA5ESG+gGiUq8uWXEIshlyMsDLGxMDDgOhAhvYRqlPQYwyAiAmFhYBiEhyM6Gjr0e0X6EPp154OOFkxSKpGSgsZGLrMpFPjf/8VHH0FPD999h4gILsMQwgWqUY33nAWTLlzArFmwtcWCBYiNRU1Nb2erqcG8efj+e5iaIiEBK1f2dgBCNABd8KTZmpowZAhWrmyxRufq1ZDJcO0akpOxaROuXn0ybm6OOXMQHAw/P5iZqT1bWRlEImRkYOBAJCZi8mS1vyMhGolqVLNduoTx43H9Olxdnw1euICpU3H79pPbKwsLkZAAiQQZGWD/No2MIBRCJML8+c/uZFet27fh64vr1+HkhOTkFvEI6WOoRjXbqVPw9UVVVYu9y8JCuLggKwvjxrV4cVERjh5FYiLS0tDUBAC6upg0CWIxFi7EoEGqCpWbmzt8+XLD7GyMGYOkJNjZqWrLhPARnRvVbE8XTGqu1YJJTzk6Yv16pKTg9m3s3o3XXgOACxewYQMcHPDaa4UxMcXFxT1MdO7cOW9v76CKioaAAKSlUYcSQjWq2TqzYFJb9vZYuxapqSgtxd698PeHnh7Onv337t0ODg6jRo2KiIj4/fffuxHn+PHjfn5+FRUVhh4eyoMHYWHRjY0Qom0YouH272eMjJidO5n795m7d5mPP2YMDZkjRxi5vAsbqahQxsUtCg42bnaD5tixYz/55JO8vLxObuOHH37Q09MDsG7dOoVC0Z3PQog2ohrlg6QkxsuLMTZmTE2ZGTOYc+eYDz9kTEwYf39m717m8ePOb6m2tjYhISE0NLRfszs1XVxcwsLCZDKZUqns6Ae3b98OQCAQhIeHq+ATEaJFaIqJnxYtwsGDT/7bxAR+fggOxty5nb+Nvb6+XiaTSaXSAwcOlJaWsoNOTk4BAQFisdjLy0sgELCDCoVi7dq133zzja6u7tdff71q1SpVfxhC+I1qlLdu38axYy2uczI0xLRp8PfHokWdX6GusbHx7Nmzhw8fPnbs2NM+HTJkSHBwcGRkZGNjY0hISHx8vImJycGDB/39/dX0aQjhL6pR/rtzB0lJkEqRnNz6Oiex+Mnj5DpBqVRmZGQkJibGx8f/8ccfU6ZMSUxMnDdvnkwmGzBggFQq9fLyUuOnIIS3qEa1yP37OHYMhw8jLe3JjfY6OvDyQnBwRXBw/04/9YhhmIsXL965c2fbtm2XL1+2t7dPTk52d3dXY3JC+IxqVBuVl0MqRWIikpJQUwMTE1uBwGroULFYvGjRohEjRrxwA3l5eb6+vkVFRW5ubsnJyUPo4Z2EdIxqVKtVVeHEiaKMDLc9e2r+WrjEw8MjKCgoKCho9OjR7f6QUqkcM2ZMbm6ut7f38ePH+/fv34uJCeEfqtE+QS6Xp6SkSCQSqVRaUVHBDjo7O4tEolbz8gB++eWXZcuW1dfX5+XlGRkZcRSZEN6gGu1bFApFZmamRCI5dOhQSUkJO+jo6BgYGCgSiV599VU9Pb38/Pzhw4fb2dndbXUTKiGkPVSjfZRCoTh37hx7ndPTurS3tw8MDAwNDfX393/06NGff/7p0OmJKUL6LKrRvk6pVObk5Eil0p9++om90f6DDz64ePHiqVOn4uPjg4ODuQ5IiKajpUn6Oh0dnXHjxkVERNy4cSMrK+u9995buHDhxIkTAVy8eJHrdITwgB7XAYgGGTdu3Lhx4wCwNfrLL79wnYgQHqCDetKOR48eWVtbm5iYVFRUsKs6EUI6Qgf1pB2WlpbOzs41NTV5eXlcZyFE01GNkvbRcT0hnUQ1StpHNappIiIidu3axXUK0g467UXaN2HCBNBkvWZQKBRr1qz59ttv9fX1AwICnJycuE5EWqApJtI+uVxuYWGhUCjKy8vNzc25jtN31dfXh4aGSiQSQ0PDuLi4119/netEpDU6qCftMzIyeuWVVxQKRXZ2NtdZ+q6KiopZs2ZJJJIBAwakpKRQh2omqlHSITo9yq179+69+uqr58+ft7e3T0tLmzZtGteJSPuoRkmH6PQoh/Ly8iZNmnT58uWRI0dmZma+8sorXCciHaIaJR2ivVGuXLx4cfr06UVFRRMnTjx//ryjoyPXicjzUI2SDg0fPnzAgAF37twpLi7mOksfkpiYOGPGjAcPHohEojNnzlhZWXGdiLwA1SjpkEAg8PT0BB3X96LY2NigoKDa2tply5YdOXLExMSE60TkxahGyfPQ6dHeFB0dvXz58sbGxi1btsTExNBqBnxBf0/keej0aO9gGGbz5s2RkZG6urq7du1as2YN14lIF9DeKHketkazsrIUCgXXWXrm9Gl4e8PUFObmEAqRns51oGcaGhqWLFkSGRlpaGi4f/9+6lDeoRolz2Ntbe3k5FRVVcXvpZ4OHUJAAAIDUVCA69fh7Q2hECkpXMcCgOrqapFIdODAATMzM6lUumDBAq4TkS6jGiUvwPuV8JuasH49Nm3Cxo2wtcXgwfjwQyxbhvXruU6GkpISb2/v06dPDxo0SCaTzZw5k+tEpDuoRskL8P706OXLKClBaGiLwaVLkZeHoiKOMgFAQUHBtGnTcnJyhg0bJpPJxowZw2EY0hNUo+QFeD9Z//AhAAwe3GLQ3h4Azp3DggU4eBBVVb0cKisra/LkyTdv3vT09MzMzHzppZd6OQBRIapR8gJjx47V19fPzc2tqanhOku3sJev//UQ6SfYL9PTIZFg0SLY2EAkwrfforS0FxKlpqb6+PiUlpYKhcLU1FRra+teeFOiPlSj5AWMjY3d3d15vNSThwdsbREX12IwLg6urnj/ffzzn5g6FQ0NSEzE6tVwcMDs2Wrt03379vn5+T1+/DgkJCQpKYkWIdQGDCEv8ve//x3A559/znWQLvrxRyYri2EYZv9+xsiI2bmTuX+fuXuX+fhjxtCQSU5+9soHD5i9exl/f0ZfnwEYgNHRYby8mO3bmZs3VZgoKipKR0cHQFhYmEKhUOGWCYeoRsmLxcTEABCLxVwH6YpPP2UEAsbGhikrYxiGSUpivLwYY2PG1JSZMYM5d679nyore9KnhoZP+hRg3NyY8HDm+vWexFEqlVu2bAEgEAg+++yznmyKaBqqUfJi165dA+Do6Mh1kM5RKplNmxiAEQiYyMhubqSyktm3jwkOZkxMnvbphtdf37Zt29WrV7u6scbGxhUrVgAwMDDYv39/NyMRTUUPESEvxjDMwIEDKyoqiouL7dk5bo3V0IBly3DgAAwMEBuLhQt7usG6Ovz8MySSyrNnbUpLGxoaALi4uPj7+4vFYi8vL4FA8PwN1NTUiMXikydPmpmZSSQSX1/fnkYimobrHif8IBQKARw7dozrIM9VVcXMns0AjJkZc+qUarfdIJcnJiauWLHC0tLy6T8fZ2fnTZs2ZWZmKpXKdn/q4cOHkydPBmBra5vFnqglWodqlHTK1q1bAWzdupXrIB27d48ZO5YBmEGDmOxs9b1PU1OTTCYLCwtrvmNubW0dGhqakJDQ0NDw9JWFhYXDhw9n2/bGjRvqi0S4RQf1pFMSEhLmzZvn4+Pz888/c52lPQUF8PVFfj5cXHDqFHrlanalUpmTkyOVSvft23fz5k120NLScs6cOWKx2MHBQSQSFRcXu7u7Jycna/rJENIDVKOkU0pLS21tbc3NzSsqKthLdjRIVhbmzkVpKTw9kZgIG5tefn+GYbKysg4fPnz48OGnfWpmZlZdXS0UCkNCQiQSyapVq+bNm9fLwUjvoBolneXk5HT79u2rV6+6ublxnaWZM2cwfz4eP4aPD44eBadXszc0NJw4ceL333+XSqWTJk2qra2Nior64osv3n///fXr10dFRXGYjagPLdtMOout0bfeeuuDDz6YMmWKRuyT7tuHFSvQ0ICQEOzZA319DrM0NTVZWVnV1tZWVFRs2bJFoVDo6upCCxYlIC+iAf8SCB/8+OOP6enplpaWp0+fnjZtmq2t7dKlS6VSaWNjI1eRHv373wgNRUMDNm9GbCy3HQpAT09vxIgRCoUiKysLANuhADw9PXV0dLKzs+vr6zkNSNSFapS8WHR09IoVKxQKxcyZM99++21nZ+eHDx/+5z//CQgIsLOzW7lyZVJSUm92BMMwERERHu+8UztiBHbswGef4UUXb/YOdsez1aKCFhYWrq6u9fX1ubm5HOUiasbpdQJE0ymVyk2bNgEQCARffPHF0/ErV66Eh4c3P0lqYmLi7++/d+/eqqoqtUZqaGgICQkBYGBgcPinn9T6Xl0VGxsLICgoqNX48uXLAezevZuTVETdqEZJh+rr6xcvXswW1k8dFNYff/wRFRXV/GYeY2Njtk8rKipUHqm6utrPzw+AqanpyZMnVb79Hrpx4wYAe3v7VuNfffUVgKVLl3KSiqgb1ShpX1VV1ezZswGYmZmd6sQdQYWFha361NDQUCgURkVF3b9/XyWRHj16xN4RZGlpmZmZqZJtqpZSqWTvcbpz507z8UuXLgFwdXXlKhhRK6pR0o6SkpKxY8cCGDRoUHYX7wgqKir65ptv/P39nz5mXVdX18vLKyoqqri4uNuRCgsLXV1dofF3BM2aNQvA4cOHmw82NDSYmJgIBIIydrkpol2oRklrBQUFL7/8MgAXF5f8/Pxub+fhw4d79+719/c3MDBg+1RHR8fLy2v79u1d3Wxubu7gwYMBuLu7t9rR0zQffPABgC1btrQa9/LyAtCZ/XrCO1SjpIWsrCwbGxsA48ePV9XBeFlZ2d69e8Visamp6dMpKTc3t/Dw8OudWMTz7NmzFhYWAGbMmFFZWamSSOqTmJgI4NVXX201vnHjRgD/+Mc/OElF1IpqlDyTmprar18/AD4+Ps8vrF9++WXp0qXHjx+vq6vr/PZramoSEhJCQ0ObPzmD7dOOVj86cuSIkZERO/3dpffiSmlpKXtCuampqfn4gQMHAPj7+3MVjKgP1Sh5Ij4+ni2skJCQ5ssUtWvDhg2t5uW7tJ9YV1fH9im7m8lydnYOCwuTyWRPF53bvXs3e6/UunXrePTIDWdnZwC5ubnNBwsLCwFYW1tzlYqoD9UoYRiGiY6O7tIzgvLz8z/99NNx48Y9LUFjY+PAwMDY2Njy8vLOv69cLj9x4kTbRTzffvvtVatWsderhoeHd/+DcWHRokUAvv/++1bjtra2AAoKCjhJRdSHarSvUyqV4eHhbGFt3769qz9++/btqKgooVDYdl7+3r17nd9O20U8raysdHV1Y2JiuhqJczt37gSwatWqVuP+/v4AOroCl/AX1Wif1tjYuHLlSgB6eno//PBDTzb14MEDdl5e/69729l5+aioqC7NrSsUCplMxl42NH369J5E4sqFCxcAeHh4tBrftm0bgHfefpuTVER9qEb7rurq6jlz5rB3BCUlJalqs48ePWL71NDQsNU8Uuev97x69SqAoUOHqipVb6qrqzMwMNDV1W11X2zJmTN1L72knDqVq2BETWi90T6qrKxMJBJlZGRYWlpKpVL27iDVqq2tTU1NlUgkR48era6uZgfd3NzEYvGCBQuev2gp89dD9O7evWtnZ6fybOq2fNYs96qq0MhImylTno1WVmLgQBgY4PFjztejIqrEdY8TDvTyHUG1tbXsvDx7NRXLxcWl1bx8Kz4+PgCOHz+u7nhqsWYNAzBtn0c/YgQDMJcucZGJqAstlNfnXLlyZdq0aTdu3Bg9erRMJmOfuaZWxsbGIpEoNjb2/v37CQkJy5cvHzhwYEFBwa5du6ZNm+bu7q5UKtv+1MSJE9Fm0TnemDABANou1cyO8/RDkQ5QjfYtaWlpU6dOvXPnzowZM9LT09k7LHuNkZGRSCSKiYkpLS1l5+Xt7OyGDRvW7kL6/F40fuJEoL26ZMd5+qFIB+jcaB9y7NixxYsXy+Xy+fPn79+/n73YnlsKhaK8vNzKyqrtt0pKSuzs7CwsLMrKyjTigSVdwjCwtER5OYqL0fyZoFlZ8PTEyJG4do27cETF+PbbSbrrX//6V3BwsFwuX7du3dMbljinq6vbbocCGDRokKOjY2VlJbuIJ88IBBg/Hmiz4+nhAWNjXL+OigpOchF1oBrtE3bs2LFu3TqGYcLDw7/88ku+7Nzx+/Rou8fv+voYMwYMg6wsTkIRdeDHPyfSbQqFYvXq1e+++66uru4333wTERHBdaIu4Pfp0Y5mk+j0qNahByxrs/r6+pCQkPj4eBMTk4MHD7I3I/KINuyN/vorFAr89ZRQgCbrtRBNMWmt8vLygICA9PT0AQMGSKVSdtlgfqmtrbWwsBAIBBUVFSYmJlzH6TpnZ9y6hStXMGrUs8GCAgwbBhsb3L/PXTKiSnRQr53Ky8unTp2anp4+dOjQjIwMPnYoABMTk9GjRzc2Nubk5HCdpVva3fF0cYGNDUpLUVTESSiiclSj2unPP/80MjJydHSUyWQjRozgOk738fv0aEenQdlJfDqu1xZUo9opIyMjOzt7ypQpQ4YM4TpLj2jD6VGaZdJ2VKPaiW2f7OxsroP0FLs3ytcaHTsW+vrIzUVNTYtxmmXSLlSj2snd3d3U1DQ/P//Ro0dcZ+kRNze3fv363bp16z4fJ2SMjTF6NBQKtDq3O2ECdHRadyvhLapR7aSnpzdmzBiGYbLaXOa9Z8+eN954o66ujpNgXaWjo8M+quTXX3/lOku3tHtcP3Agystx6RIniYjKUY1qLfa4vu3kzK5du3744QcezX3z+/RoQADeegttl3NttmYg4TuqUa3VUfvwrpX4PVnv54d//hPNF28mWodqVGt11D68a6VJkyYBuHjxYrvLkvLJ6dPw9oapKczNIRQiPZ3rQEQ1qEa1lpOTk62t7YMHD9gnpD/Fu71ROzs7BweHioqK/Px8rrP0wKFDCAhAYCAKCnD9Ory9IRQiJYXrWEQFqEa1WbtXC7Fz34WFhTya++Zd9bfW1IT167FpEzZuhK0tBg/Ghx9i2TKsX891MqICVKParN324ePcN+9ORLR2+TJKShAa2mJw6VLk5dEtoVqAalSbddQ+HU3iayze740+fAgArR7Zwq6K/+ABB3mISlGNarMJEybo6OhkZ2c3NDS0GgevWmn8+PF6enqXL1/my+WurbEr/N+922KQ/dLamoM8RKWoRrWZhYXF8OHD5XJ5bm5u83HezX2bmpq6ubk1Njb+9ttvXGfpFg8P2NoiLq7FYFwcXF3h6MhRJqIyVKNart3DYT7OffP79KieHnbuxOefIyoKpaW4dw+ffIKYGERHc52MqADVqJZ7/tWjPDqu5/3p0cWLceQI4uPh5ISXX0ZqKk6fxuzZXMciKkA1quWefy8Tj3bueF+jAPz8kJ6O2lpUV+PMGXh7cx2IqAbVqJZ75ZVXjI2Nb9y4UV5e3nycd600atQoc3PzgoKC0tJSrrMQ0gLVqJbT19f/29/+1napJ97Nfevo6IwdOxZA2zWrCOEW1aj2a/f0qKmp6ciRI/k19827PWjSR1CNaj+tOT3Ku2kx0kdQjWo/fq2Y9/jx4+3bt7d7QSsb+Ndff6WnghONQjWq/ZydnW1sbEpLS2/dutV8XAN37kpKSqZPn/7ee+999NFHbb/r4OAwePDgsrKymzdv9n42QjpCNdoneHp6os3xu6bNfRcUFHh7e//222/Dhg0LbbWKBwBALpfr6+u7urqWlZX1fjxCOkI12ie0e/yuq6urOXPfWVlZkydPzs/P9/T0zMzMfOmll1q9oLy8fNasWbdu3dLX12f/r0CIhqAa7RM6On7XkNOjqampPj4+paWlQqEwNTXVus1qHffu3XvttddkMpmdnV1cXJyODv3eEg1Cv459wsSJEwUCwaVLlxobG5uPa8Lp0X379vn5+T1+/DgkJCQpKcnc3LzVC/Ly8iZNmvTbb7+NHDnyv//9r4eHByc5CekI1Wif0L9//5dfflkul1+5cqX5OOdz39HR0UuXLm1sbAwLC9u7d6++vn6rF1y8eHH69OlFRUUTJ048f/68I62HRDQP1Whf0e7xO4dz3wzDvPvuuxs2bGAY5rPPPouOjm57qJ6SkiIUCh88eCASic6cOWPFrtpJiIahGu0rNGqpp6ampjfeeGPHjh0GBgZxcXGbN29u+5rY2Ni5c+dWVVUtW7bsyJEjJiYmvZmQkC5gSN/AFqibm1ur8RMnTnz99de3b9/utSTV1dV+fn4AzMzMTp482e5roqKiBAIBgC1btiiVyl7LRkg3CBi6IaRvaGxstLCwkMvlZWVl/fv35ypGWVmZv79/Zmamra3tiRMn2IfrNccwzObNmyMjIwUCwRdffLFx40ZOchLSeXRQ31fo6+t7eHgwDHPp0iWuMty6dWvKlCmZmZnOzs7nz59v26ENDQ1LliyJjIw0NDT86aefqEMJL1CN9iHcXt505cqVqVOn3rhxw93dPT09ffjw4a1eUF1dLRKJDhw4YGZmJpVKFy5cyElOQrpKj+sApPdwuKRTWlpaYGBgZWXljBkzjh071q9fv1YvKCkpmTNnTk5OzqBBg06ePDlmzJjeD0lI99DeaB/C1T1LR48e9fPzq6ysDAoKSkpKatuh7N30OTk5w4YNk8lk1KGEX6hG+xAXFxcrK6uSkhKZTNZrb7p79+7XX39dLpe/+eabEonEyMio1QteeDc9IRqOarQPEQgEa9euHTlypLe396hRoyIiItQ63cQwTERExJtvvskwTHh4+K5du9peYP/Cu+kJ4QFOL7ciHFi1apWFhcXTX4ARI0Zs3br10qVLqn2XpqamVatWAdDV1f3uu+/afU1cXBx792dISEhDQ4NqAxDSa6hG+6KmpiaZTBYWFmZra/u0Tx0dHcPCwmQymUKh6Plb1NfXz5w509TUNDExsd0XREVFsTunYWFhKnlHQrhCl9/3aQqFIjMzUyKRxMfH3717lx20trb29fUVi8W+vr5t1wrpvMrKyj/++INd0rQ5hmHee++9HTt2CASCHTt2tHsnKCE8QjVKAECpVObk5Eil0v379+fn57ODAwcOnDt3rlgsnj17toGBgUreqKmpafXq1Xv27DEwMIiJiVmyZIlKNksIh6hGSWtXr16VSCSHDh3Ky8tjR/r37z9z5kx/f/+goCAzM7Nub7mmpkYsFp88edLMzEwikfj6+qooMiFcoholHbp69WpiYqJUKr1w4QI7Ymxs7OPjIxaLAwMD217++XyPHj0SiUTPuZueEJ6iGiUvduvWrePHj0skkoyMDPYXxsjISCgUikSiwMBAGxubzmxh9uzZv//+u7Ozc3Jycts7QQnhL6pR0gV//vnnkSNHEhMT09LSmpqaAOjq6k6aNEksFi9YsMDOzq7dn7py5Yqvr29xcbG7u3tycrK9vX3vpiZEvahGSXc8fPgwKSlJIpGcOnWKfb6Tjo7O5MmTxWJxcHCwg4PD01eePXt2/vz5z7mbnhC+oxolPVJeXi6VSiUSSUpKSn19PTvo5uYmFosXL1587dq1JUuWyOXyoKCgffv2tb0TlBAtQDVKVKOysvLEiROHDx9OTk6ura1lB3V1dRUKxYYNGyIjI+mpyERbUY0SFaurq/v5558lEklCQsKaNWv69+//zjvvcB2KEDWiGiXqUl9fb2BgwD5SiRAtRjVKCCE9QqerCCGkR6hGCSGkR6hGCSGkR6hGCSGkR/4fopAJphO25WkAAAAASUVORK5CYII=\n"
},
"metadata": {}
}
]
},
{
"metadata": {
"trusted": true
},
"cell_type": "code",
"source": "#disable the monkey-patched highlighting:\nChem.Mol.GetSubstructMatch = Chem.Mol.__GetSubstructMatch\nChem.Mol.GetSubstructMatches = Chem.Mol.__GetSubstructMatches\nprint(m.GetSubstructMatches(Chem.MolFromSmarts('C(=O)O')))\nm",
"execution_count": 8,
"outputs": [
{
"output_type": "stream",
"text": "((3, 4, 5), (8, 9, 7))\n",
"name": "stdout"
},
{
"output_type": "execute_result",
"execution_count": 8,
"data": {
"text/plain": "<rdkit.Chem.rdchem.Mol at 0x7eff54f83940>",
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAcIAAACWCAIAAADCEh9HAAAABmJLR0QA/wD/AP+gvaeTAAAeKUlEQVR4nO3deVQT594H8G/YV1HZBBEBW1GU4lVxQ7GWqKAEERq3F9Sj9fVctWit1taeFnpt+2pbrmC97e0mlYtWDW4EEaUoGoRbi1APKlosKBVFVBbZwpLM+8dYZVWWhMmE3+f4R3kIk2+O+O3MPDPPCBiGASGEkO7S4ToAIYTwG9UoIYT0CNUoUYPTp+HtDVNTmJtDKER6+rNvTZ2Kjz9u8WKBAGlpvZuPEFWiGiWqdugQAgIQGIiCAly/Dm9vCIVISeE6FiHqosd1AKJdmpqwfj02bcLGjU9GPvwQxcVYvx7XrnGajBB1ob1RolKXL6OkBKGhLQaXLkVeHoqKOMpEiHpRjRKVevgQAAYPbjFobw8ADx48+fL//g9WVs/+EMJzdFBPVIqtxbt3MXz4s8G7dwHA2vrJl2vXIizs2XeHDOm1dISoA+2NEpXy8ICtLeLiWgzGxcHVFY6OT77s1w8ODs/+EMJztDdKVEpPDzt3YsUKDByIJUugUGDPHsTE4PhxrpMRoi5Uo0TVFi9GXR2++gpbt0JHBxMmPLmMlBAtJaB76omKZWdjzhyYmeHCBdjacp2GELWjc6NEpc6exYwZuH8fQ4fC2JjrNIT0BqpRojpHjmDOHDx+jP/5H5w8iX79uA5ESG+gGiUq8uWXEIshlyMsDLGxMDDgOhAhvYRqlPQYwyAiAmFhYBiEhyM6Gjr0e0X6EPp154OOFkxSKpGSgsZGLrMpFPjf/8VHH0FPD999h4gILsMQwgWqUY33nAWTLlzArFmwtcWCBYiNRU1Nb2erqcG8efj+e5iaIiEBK1f2dgBCNABd8KTZmpowZAhWrmyxRufq1ZDJcO0akpOxaROuXn0ybm6OOXMQHAw/P5iZqT1bWRlEImRkYOBAJCZi8mS1vyMhGolqVLNduoTx43H9Olxdnw1euICpU3H79pPbKwsLkZAAiQQZGWD/No2MIBRCJML8+c/uZFet27fh64vr1+HkhOTkFvEI6WOoRjXbqVPw9UVVVYu9y8JCuLggKwvjxrV4cVERjh5FYiLS0tDUBAC6upg0CWIxFi7EoEGqCpWbmzt8+XLD7GyMGYOkJNjZqWrLhPARnRvVbE8XTGqu1YJJTzk6Yv16pKTg9m3s3o3XXgOACxewYQMcHPDaa4UxMcXFxT1MdO7cOW9v76CKioaAAKSlUYcSQjWq2TqzYFJb9vZYuxapqSgtxd698PeHnh7Onv337t0ODg6jRo2KiIj4/fffuxHn+PHjfn5+FRUVhh4eyoMHYWHRjY0Qom0YouH272eMjJidO5n795m7d5mPP2YMDZkjRxi5vAsbqahQxsUtCg42bnaD5tixYz/55JO8vLxObuOHH37Q09MDsG7dOoVC0Z3PQog2ohrlg6QkxsuLMTZmTE2ZGTOYc+eYDz9kTEwYf39m717m8ePOb6m2tjYhISE0NLRfszs1XVxcwsLCZDKZUqns6Ae3b98OQCAQhIeHq+ATEaJFaIqJnxYtwsGDT/7bxAR+fggOxty5nb+Nvb6+XiaTSaXSAwcOlJaWsoNOTk4BAQFisdjLy0sgELCDCoVi7dq133zzja6u7tdff71q1SpVfxhC+I1qlLdu38axYy2uczI0xLRp8PfHokWdX6GusbHx7Nmzhw8fPnbs2NM+HTJkSHBwcGRkZGNjY0hISHx8vImJycGDB/39/dX0aQjhL6pR/rtzB0lJkEqRnNz6Oiex+Mnj5DpBqVRmZGQkJibGx8f/8ccfU6ZMSUxMnDdvnkwmGzBggFQq9fLyUuOnIIS3qEa1yP37OHYMhw8jLe3JjfY6OvDyQnBwRXBw/04/9YhhmIsXL965c2fbtm2XL1+2t7dPTk52d3dXY3JC+IxqVBuVl0MqRWIikpJQUwMTE1uBwGroULFYvGjRohEjRrxwA3l5eb6+vkVFRW5ubsnJyUPo4Z2EdIxqVKtVVeHEiaKMDLc9e2r+WrjEw8MjKCgoKCho9OjR7f6QUqkcM2ZMbm6ut7f38ePH+/fv34uJCeEfqtE+QS6Xp6SkSCQSqVRaUVHBDjo7O4tEolbz8gB++eWXZcuW1dfX5+XlGRkZcRSZEN6gGu1bFApFZmamRCI5dOhQSUkJO+jo6BgYGCgSiV599VU9Pb38/Pzhw4fb2dndbXUTKiGkPVSjfZRCoTh37hx7ndPTurS3tw8MDAwNDfX393/06NGff/7p0OmJKUL6LKrRvk6pVObk5Eil0p9++om90f6DDz64ePHiqVOn4uPjg4ODuQ5IiKajpUn6Oh0dnXHjxkVERNy4cSMrK+u9995buHDhxIkTAVy8eJHrdITwgB7XAYgGGTdu3Lhx4wCwNfrLL79wnYgQHqCDetKOR48eWVtbm5iYVFRUsKs6EUI6Qgf1pB2WlpbOzs41NTV5eXlcZyFE01GNkvbRcT0hnUQ1StpHNappIiIidu3axXUK0g467UXaN2HCBNBkvWZQKBRr1qz59ttv9fX1AwICnJycuE5EWqApJtI+uVxuYWGhUCjKy8vNzc25jtN31dfXh4aGSiQSQ0PDuLi4119/netEpDU6qCftMzIyeuWVVxQKRXZ2NtdZ+q6KiopZs2ZJJJIBAwakpKRQh2omqlHSITo9yq179+69+uqr58+ft7e3T0tLmzZtGteJSPuoRkmH6PQoh/Ly8iZNmnT58uWRI0dmZma+8sorXCciHaIaJR2ivVGuXLx4cfr06UVFRRMnTjx//ryjoyPXicjzUI2SDg0fPnzAgAF37twpLi7mOksfkpiYOGPGjAcPHohEojNnzlhZWXGdiLwA1SjpkEAg8PT0BB3X96LY2NigoKDa2tply5YdOXLExMSE60TkxahGyfPQ6dHeFB0dvXz58sbGxi1btsTExNBqBnxBf0/keej0aO9gGGbz5s2RkZG6urq7du1as2YN14lIF9DeKHketkazsrIUCgXXWXrm9Gl4e8PUFObmEAqRns51oGcaGhqWLFkSGRlpaGi4f/9+6lDeoRolz2Ntbe3k5FRVVcXvpZ4OHUJAAAIDUVCA69fh7Q2hECkpXMcCgOrqapFIdODAATMzM6lUumDBAq4TkS6jGiUvwPuV8JuasH49Nm3Cxo2wtcXgwfjwQyxbhvXruU6GkpISb2/v06dPDxo0SCaTzZw5k+tEpDuoRskL8P706OXLKClBaGiLwaVLkZeHoiKOMgFAQUHBtGnTcnJyhg0bJpPJxowZw2EY0hNUo+QFeD9Z//AhAAwe3GLQ3h4Azp3DggU4eBBVVb0cKisra/LkyTdv3vT09MzMzHzppZd6OQBRIapR8gJjx47V19fPzc2tqanhOku3sJev//UQ6SfYL9PTIZFg0SLY2EAkwrfforS0FxKlpqb6+PiUlpYKhcLU1FRra+teeFOiPlSj5AWMjY3d3d15vNSThwdsbREX12IwLg6urnj/ffzzn5g6FQ0NSEzE6tVwcMDs2Wrt03379vn5+T1+/DgkJCQpKYkWIdQGDCEv8ve//x3A559/znWQLvrxRyYri2EYZv9+xsiI2bmTuX+fuXuX+fhjxtCQSU5+9soHD5i9exl/f0ZfnwEYgNHRYby8mO3bmZs3VZgoKipKR0cHQFhYmEKhUOGWCYeoRsmLxcTEABCLxVwH6YpPP2UEAsbGhikrYxiGSUpivLwYY2PG1JSZMYM5d679nyore9KnhoZP+hRg3NyY8HDm+vWexFEqlVu2bAEgEAg+++yznmyKaBqqUfJi165dA+Do6Mh1kM5RKplNmxiAEQiYyMhubqSyktm3jwkOZkxMnvbphtdf37Zt29WrV7u6scbGxhUrVgAwMDDYv39/NyMRTUUPESEvxjDMwIEDKyoqiouL7dk5bo3V0IBly3DgAAwMEBuLhQt7usG6Ovz8MySSyrNnbUpLGxoaALi4uPj7+4vFYi8vL4FA8PwN1NTUiMXikydPmpmZSSQSX1/fnkYimobrHif8IBQKARw7dozrIM9VVcXMns0AjJkZc+qUarfdIJcnJiauWLHC0tLy6T8fZ2fnTZs2ZWZmKpXKdn/q4cOHkydPBmBra5vFnqglWodqlHTK1q1bAWzdupXrIB27d48ZO5YBmEGDmOxs9b1PU1OTTCYLCwtrvmNubW0dGhqakJDQ0NDw9JWFhYXDhw9n2/bGjRvqi0S4RQf1pFMSEhLmzZvn4+Pz888/c52lPQUF8PVFfj5cXHDqFHrlanalUpmTkyOVSvft23fz5k120NLScs6cOWKx2MHBQSQSFRcXu7u7Jycna/rJENIDVKOkU0pLS21tbc3NzSsqKthLdjRIVhbmzkVpKTw9kZgIG5tefn+GYbKysg4fPnz48OGnfWpmZlZdXS0UCkNCQiQSyapVq+bNm9fLwUjvoBolneXk5HT79u2rV6+6ublxnaWZM2cwfz4eP4aPD44eBadXszc0NJw4ceL333+XSqWTJk2qra2Nior64osv3n///fXr10dFRXGYjagPLdtMOout0bfeeuuDDz6YMmWKRuyT7tuHFSvQ0ICQEOzZA319DrM0NTVZWVnV1tZWVFRs2bJFoVDo6upCCxYlIC+iAf8SCB/8+OOP6enplpaWp0+fnjZtmq2t7dKlS6VSaWNjI1eRHv373wgNRUMDNm9GbCy3HQpAT09vxIgRCoUiKysLANuhADw9PXV0dLKzs+vr6zkNSNSFapS8WHR09IoVKxQKxcyZM99++21nZ+eHDx/+5z//CQgIsLOzW7lyZVJSUm92BMMwERERHu+8UztiBHbswGef4UUXb/YOdsez1aKCFhYWrq6u9fX1ubm5HOUiasbpdQJE0ymVyk2bNgEQCARffPHF0/ErV66Eh4c3P0lqYmLi7++/d+/eqqoqtUZqaGgICQkBYGBgcPinn9T6Xl0VGxsLICgoqNX48uXLAezevZuTVETdqEZJh+rr6xcvXswW1k8dFNYff/wRFRXV/GYeY2Njtk8rKipUHqm6utrPzw+AqanpyZMnVb79Hrpx4wYAe3v7VuNfffUVgKVLl3KSiqgb1ShpX1VV1ezZswGYmZmd6sQdQYWFha361NDQUCgURkVF3b9/XyWRHj16xN4RZGlpmZmZqZJtqpZSqWTvcbpz507z8UuXLgFwdXXlKhhRK6pR0o6SkpKxY8cCGDRoUHYX7wgqKir65ptv/P39nz5mXVdX18vLKyoqqri4uNuRCgsLXV1dofF3BM2aNQvA4cOHmw82NDSYmJgIBIIydrkpol2oRklrBQUFL7/8MgAXF5f8/Pxub+fhw4d79+719/c3MDBg+1RHR8fLy2v79u1d3Wxubu7gwYMBuLu7t9rR0zQffPABgC1btrQa9/LyAtCZ/XrCO1SjpIWsrCwbGxsA48ePV9XBeFlZ2d69e8Visamp6dMpKTc3t/Dw8OudWMTz7NmzFhYWAGbMmFFZWamSSOqTmJgI4NVXX201vnHjRgD/+Mc/OElF1IpqlDyTmprar18/AD4+Ps8vrF9++WXp0qXHjx+vq6vr/PZramoSEhJCQ0ObPzmD7dOOVj86cuSIkZERO/3dpffiSmlpKXtCuampqfn4gQMHAPj7+3MVjKgP1Sh5Ij4+ni2skJCQ5ssUtWvDhg2t5uW7tJ9YV1fH9im7m8lydnYOCwuTyWRPF53bvXs3e6/UunXrePTIDWdnZwC5ubnNBwsLCwFYW1tzlYqoD9UoYRiGiY6O7tIzgvLz8z/99NNx48Y9LUFjY+PAwMDY2Njy8vLOv69cLj9x4kTbRTzffvvtVatWsderhoeHd/+DcWHRokUAvv/++1bjtra2AAoKCjhJRdSHarSvUyqV4eHhbGFt3769qz9++/btqKgooVDYdl7+3r17nd9O20U8raysdHV1Y2JiuhqJczt37gSwatWqVuP+/v4AOroCl/AX1Wif1tjYuHLlSgB6eno//PBDTzb14MEDdl5e/69729l5+aioqC7NrSsUCplMxl42NH369J5E4sqFCxcAeHh4tBrftm0bgHfefpuTVER9qEb7rurq6jlz5rB3BCUlJalqs48ePWL71NDQsNU8Uuev97x69SqAoUOHqipVb6qrqzMwMNDV1W11X2zJmTN1L72knDqVq2BETWi90T6qrKxMJBJlZGRYWlpKpVL27iDVqq2tTU1NlUgkR48era6uZgfd3NzEYvGCBQuev2gp89dD9O7evWtnZ6fybOq2fNYs96qq0MhImylTno1WVmLgQBgY4PFjztejIqrEdY8TDvTyHUG1tbXsvDx7NRXLxcWl1bx8Kz4+PgCOHz+u7nhqsWYNAzBtn0c/YgQDMJcucZGJqAstlNfnXLlyZdq0aTdu3Bg9erRMJmOfuaZWxsbGIpEoNjb2/v37CQkJy5cvHzhwYEFBwa5du6ZNm+bu7q5UKtv+1MSJE9Fm0TnemDABANou1cyO8/RDkQ5QjfYtaWlpU6dOvXPnzowZM9LT09k7LHuNkZGRSCSKiYkpLS1l5+Xt7OyGDRvW7kL6/F40fuJEoL26ZMd5+qFIB+jcaB9y7NixxYsXy+Xy+fPn79+/n73YnlsKhaK8vNzKyqrtt0pKSuzs7CwsLMrKyjTigSVdwjCwtER5OYqL0fyZoFlZ8PTEyJG4do27cETF+PbbSbrrX//6V3BwsFwuX7du3dMbljinq6vbbocCGDRokKOjY2VlJbuIJ88IBBg/Hmiz4+nhAWNjXL+OigpOchF1oBrtE3bs2LFu3TqGYcLDw7/88ku+7Nzx+/Rou8fv+voYMwYMg6wsTkIRdeDHPyfSbQqFYvXq1e+++66uru4333wTERHBdaIu4Pfp0Y5mk+j0qNahByxrs/r6+pCQkPj4eBMTk4MHD7I3I/KINuyN/vorFAr89ZRQgCbrtRBNMWmt8vLygICA9PT0AQMGSKVSdtlgfqmtrbWwsBAIBBUVFSYmJlzH6TpnZ9y6hStXMGrUs8GCAgwbBhsb3L/PXTKiSnRQr53Ky8unTp2anp4+dOjQjIwMPnYoABMTk9GjRzc2Nubk5HCdpVva3fF0cYGNDUpLUVTESSiiclSj2unPP/80MjJydHSUyWQjRozgOk738fv0aEenQdlJfDqu1xZUo9opIyMjOzt7ypQpQ4YM4TpLj2jD6VGaZdJ2VKPaiW2f7OxsroP0FLs3ytcaHTsW+vrIzUVNTYtxmmXSLlSj2snd3d3U1DQ/P//Ro0dcZ+kRNze3fv363bp16z4fJ2SMjTF6NBQKtDq3O2ECdHRadyvhLapR7aSnpzdmzBiGYbLaXOa9Z8+eN954o66ujpNgXaWjo8M+quTXX3/lOku3tHtcP3Agystx6RIniYjKUY1qLfa4vu3kzK5du3744QcezX3z+/RoQADeegttl3NttmYg4TuqUa3VUfvwrpX4PVnv54d//hPNF28mWodqVGt11D68a6VJkyYBuHjxYrvLkvLJ6dPw9oapKczNIRQiPZ3rQEQ1qEa1lpOTk62t7YMHD9gnpD/Fu71ROzs7BweHioqK/Px8rrP0wKFDCAhAYCAKCnD9Ory9IRQiJYXrWEQFqEa1WbtXC7Fz34WFhTya++Zd9bfW1IT167FpEzZuhK0tBg/Ghx9i2TKsX891MqICVKParN324ePcN+9ORLR2+TJKShAa2mJw6VLk5dEtoVqAalSbddQ+HU3iayze740+fAgArR7Zwq6K/+ABB3mISlGNarMJEybo6OhkZ2c3NDS0GgevWmn8+PF6enqXL1/my+WurbEr/N+922KQ/dLamoM8RKWoRrWZhYXF8OHD5XJ5bm5u83HezX2bmpq6ubk1Njb+9ttvXGfpFg8P2NoiLq7FYFwcXF3h6MhRJqIyVKNart3DYT7OffP79KieHnbuxOefIyoKpaW4dw+ffIKYGERHc52MqADVqJZ7/tWjPDqu5/3p0cWLceQI4uPh5ISXX0ZqKk6fxuzZXMciKkA1quWefy8Tj3bueF+jAPz8kJ6O2lpUV+PMGXh7cx2IqAbVqJZ75ZVXjI2Nb9y4UV5e3nycd600atQoc3PzgoKC0tJSrrMQ0gLVqJbT19f/29/+1napJ97Nfevo6IwdOxZA2zWrCOEW1aj2a/f0qKmp6ciRI/k19827PWjSR1CNaj+tOT3Ku2kx0kdQjWo/fq2Y9/jx4+3bt7d7QSsb+Ndff6WnghONQjWq/ZydnW1sbEpLS2/dutV8XAN37kpKSqZPn/7ee+999NFHbb/r4OAwePDgsrKymzdv9n42QjpCNdoneHp6os3xu6bNfRcUFHh7e//222/Dhg0LbbWKBwBALpfr6+u7urqWlZX1fjxCOkI12ie0e/yuq6urOXPfWVlZkydPzs/P9/T0zMzMfOmll1q9oLy8fNasWbdu3dLX12f/r0CIhqAa7RM6On7XkNOjqampPj4+paWlQqEwNTXVus1qHffu3XvttddkMpmdnV1cXJyODv3eEg1Cv459wsSJEwUCwaVLlxobG5uPa8Lp0X379vn5+T1+/DgkJCQpKcnc3LzVC/Ly8iZNmvTbb7+NHDnyv//9r4eHByc5CekI1Wif0L9//5dfflkul1+5cqX5OOdz39HR0UuXLm1sbAwLC9u7d6++vn6rF1y8eHH69OlFRUUTJ048f/68I62HRDQP1Whf0e7xO4dz3wzDvPvuuxs2bGAY5rPPPouOjm57qJ6SkiIUCh88eCASic6cOWPFrtpJiIahGu0rNGqpp6ampjfeeGPHjh0GBgZxcXGbN29u+5rY2Ni5c+dWVVUtW7bsyJEjJiYmvZmQkC5gSN/AFqibm1ur8RMnTnz99de3b9/utSTV1dV+fn4AzMzMTp482e5roqKiBAIBgC1btiiVyl7LRkg3CBi6IaRvaGxstLCwkMvlZWVl/fv35ypGWVmZv79/Zmamra3tiRMn2IfrNccwzObNmyMjIwUCwRdffLFx40ZOchLSeXRQ31fo6+t7eHgwDHPp0iWuMty6dWvKlCmZmZnOzs7nz59v26ENDQ1LliyJjIw0NDT86aefqEMJL1CN9iHcXt505cqVqVOn3rhxw93dPT09ffjw4a1eUF1dLRKJDhw4YGZmJpVKFy5cyElOQrpKj+sApPdwuKRTWlpaYGBgZWXljBkzjh071q9fv1YvKCkpmTNnTk5OzqBBg06ePDlmzJjeD0lI99DeaB/C1T1LR48e9fPzq6ysDAoKSkpKatuh7N30OTk5w4YNk8lk1KGEX6hG+xAXFxcrK6uSkhKZTNZrb7p79+7XX39dLpe/+eabEonEyMio1QteeDc9IRqOarQPEQgEa9euHTlypLe396hRoyIiItQ63cQwTERExJtvvskwTHh4+K5du9peYP/Cu+kJ4QFOL7ciHFi1apWFhcXTX4ARI0Zs3br10qVLqn2XpqamVatWAdDV1f3uu+/afU1cXBx792dISEhDQ4NqAxDSa6hG+6KmpiaZTBYWFmZra/u0Tx0dHcPCwmQymUKh6Plb1NfXz5w509TUNDExsd0XREVFsTunYWFhKnlHQrhCl9/3aQqFIjMzUyKRxMfH3717lx20trb29fUVi8W+vr5t1wrpvMrKyj/++INd0rQ5hmHee++9HTt2CASCHTt2tHsnKCE8QjVKAECpVObk5Eil0v379+fn57ODAwcOnDt3rlgsnj17toGBgUreqKmpafXq1Xv27DEwMIiJiVmyZIlKNksIh6hGSWtXr16VSCSHDh3Ky8tjR/r37z9z5kx/f/+goCAzM7Nub7mmpkYsFp88edLMzEwikfj6+qooMiFcoholHbp69WpiYqJUKr1w4QI7Ymxs7OPjIxaLAwMD217++XyPHj0SiUTPuZueEJ6iGiUvduvWrePHj0skkoyMDPYXxsjISCgUikSiwMBAGxubzmxh9uzZv//+u7Ozc3Jycts7QQnhL6pR0gV//vnnkSNHEhMT09LSmpqaAOjq6k6aNEksFi9YsMDOzq7dn7py5Yqvr29xcbG7u3tycrK9vX3vpiZEvahGSXc8fPgwKSlJIpGcOnWKfb6Tjo7O5MmTxWJxcHCwg4PD01eePXt2/vz5z7mbnhC+oxolPVJeXi6VSiUSSUpKSn19PTvo5uYmFosXL1587dq1JUuWyOXyoKCgffv2tb0TlBAtQDVKVKOysvLEiROHDx9OTk6ura1lB3V1dRUKxYYNGyIjI+mpyERbUY0SFaurq/v5558lEklCQsKaNWv69+//zjvvcB2KEDWiGiXqUl9fb2BgwD5SiRAtRjVKCCE9QqerCCGkR6hGCSGkR6hGCSGkR6hGCSGkR/4fopAJphO25WkAAAAASUVORK5CYII=\n"
},
"metadata": {}
}
]
},
{
"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.7.3",
"mimetype": "text/x-python",
"codemirror_mode": {
"name": "ipython",
"version": 3
},
"pygments_lexer": "ipython3",
"nbconvert_exporter": "python",
"file_extension": ".py"
},
"gist": {
"id": "",
"data": {
"description": "HellerSubstructureHighlighting.ipynb",
"public": true
}
}
},
"nbformat": 4,
"nbformat_minor": 2
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment