Skip to content

Instantly share code, notes, and snippets.

@iwatobipen
Created April 4, 2024 12:32
Show Gist options
  • Save iwatobipen/81a1c28e376fd5fc144c61c5101769b3 to your computer and use it in GitHub Desktop.
Save iwatobipen/81a1c28e376fd5fc144c61c5101769b3 to your computer and use it in GitHub Desktop.
Display the source blob
Display the rendered blob
Raw
{
"cells": [
{
"cell_type": "code",
"execution_count": 1,
"id": "48b30764-f427-4d15-af95-402592bbc77f",
"metadata": {},
"outputs": [],
"source": [
"from rdkit import Chem\n",
"from rdkit.Chem import rdmolops\n",
"from rdkit.Chem import Draw\n",
"from rdkit.Chem.Draw import IPythonConsole\n",
"IPythonConsole.drawOptions.addAtomIndices = True"
]
},
{
"cell_type": "code",
"execution_count": 2,
"id": "407c3511-de71-456f-b387-a7ec275dde8f",
"metadata": {},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAcIAAACWCAIAAADCEh9HAAAABmJLR0QA/wD/AP+gvaeTAAAPN0lEQVR4nO3deVCU5wEG8GdhEQQ5FSU1wQPxjKLEOwkEB9MQCKZaaqcTyMSMzphYnEbbxaCSocaucTrZVmNG0+AsGVNDMh0HJTgxEu/ERPAigDqeUYFyCAoo7LJv//gIUl0M+O5+nyzPb/xjd7+Dhx19fL9bJ4QAERE9LDetAxAR9WysUSIiKaxRIiIprFEiIimsUSIiKaxRIiIpeq0DkKu5du3ajh07LBbL888/P3bsWK3jEDkdR6PkSMeOHZs6deq1a9fq6uqefvrpgoICrRMROR1Ho+RI/v7+n3/++cyZMwE0NzdnZWXNmjVL61BEzsUaJUcKDw8PDw9XXnt7e7u7u2ubh0gF3Kgnp2hubt62bdvcuXO1DkLkdDpeU08OZ7PZXn311dbW1k8//VTrLEROx9EoOVhdXd2cOXMsFovZbNY6C5EaWKPkSKdOnZo8ebKXl9fy5ctPnTpVXFysdSIip+NGPTnSli1bvv766/a3gwYN2rBhg4Z5iFTAGiUiksKNeiIiKaxRIiIprFEiIimsUSIiKaxRIiIpvKaeHOrSJRw8iOvX0dKCgQMxZQomToQb/7cmV8YaJQc5fRpvvYUOJ422CQ/H+vWYM0eLTERq4Hmj5AgHDyIhATdvYtAgvPIKJkyAXo+LF7F9O4qLodNh/XosW6Z1SiKnYI2StJoajB+P8nLExeGzz+Dre3dSaytWroTRCDc3FBQgOlq7lETOwp1WJO2f/0R5OUJD7+1QAO7u+NvfEB8Pmw3p6RrlI3Iu1ihJy84GgDfeuLdD2ykFevgwzp9XLxWRWlijJOf6dVy6BACxsZ3OM306/P0B4MgRdUIRqYk1SnIuXmx7MXJkp/PodG1TlcIlci2sUZJTXw8Abm7o1+9Bsymj0bo6NSIRqYs1SnI8PQHAZoPF8qDZ7twBgL591YhEpC7WKMkZMKDtxX//+6DZKisBoH9/p+chUh1rlOSMGQMvLwAoKup0nrq6tmP0kZEqpSJSEWuU5PTpg+eeA4Dt2zud59//hs2GwEBMnapWLCL1sEZJ2pIlAPDZZzhwwM7UykpkZgLAwoXcN0ouiTVK0uLjkZQEmw0vvQSzGVbr3UkHDyI6GhUVCAvDqlXaRSRyIl5TT45w+zZ+/3vk5gJAUBDGjYOHB86fx+XLADB6NPLzMXSophGJnIU1Sg4iBP7zH2zZggMH2k5vcnNDZCT+8Ae88UbbeVFErog1So4mBKqrYbEgOBgeHlqnIXI61igRkRQeYiJHKCtDnz6IiLAz6U9/Qp8+2LRJ9UxEKmGNkiMIAYsFLS12JlmtsFjQ2qp6JiKVsEaJiKSwRomIpLBGiYiksEaJiKSwRomIpLBGiYiksEaJiKSwRomIpLBGiYiksEaJiKSwRomIpLBGiYiksEaJiKSwRomIpLBGiYiksEaJiKSwRomIpLBGiYiksEaJiKSwRomIpLBGiYiksEaJiKTotQ5ArqDJza3pqaesv/pVyH2TKoKC9E89Jby9gzXIRaQGjkbJAS7bbMGFhTHnzt0/6d3a2uDCwu1NTeqnIlIHa5SISAprlIhICmuUiEgKa5SISAprlIhICmuUiEgKa5SISAprlIhICmuUiEgKa5SISAprlIhICmuUiEgKa5SISAprlIhICmuUiEgKa5SISAprlIhICmuUiEgKa5SISAprlIhICmuUHGzt2rX+/v6+vr6vv/66xWLROs6jq6io6MUXX6yoqNA6CMlijZLDtLS0DB8+PD09/ebNmw0NDVlZWREREVeuXNE61yOnsrJy4cKFU6ZMyc/PX7t2rdZxSBafU08O0NTUBODixYtCCJ1OFxsbm5iY+OGHH5aUlJSWlgIoLy/XOuMjwWKxbNq0KSMjo76+3sPDY8mSJZmZmVqHImmCSEJra6vZbB44cCAAnU6XnJx89uxZZVJLS4vJZPLx8QGg1+tTU1Pr6uq0TautPXv2jB07Vvl3Fxsb++OPP2qdiByDNUoP77vvvps+fbrSC1OnTv3222/vn6e6ujo1NdXd3R1A//79TSaT1WpVP6q2zpw5Ex8fr3xRI0eO3LVrl9aJyJFYo/Qwfvrpp+TkZJ1OB2Dw4MFms9lmsz1g/qKioqioKKVHJk2atH//ftWiaqu2ttZgMPTp0wdAQECA0Whsbm7ubGabzTZv3rywsLATJ06oGZIksUapexobG41GY79+/QB4e3sbDIZbt251nKG1tTU3NzctLe3+ZXNzc4cOHaqUaUJCwoULF9TJbLVat27dunjx4tWrV1++fFmdH9pxd4ebm1tycnJlZeWDF8nKyoqOjh4+fPjRo0fVCUkOwRqlbsjNzR0yZEh7DyrHlDq6cOFCRERETEyMXq+3u4ampiaj0ejr6wugb9++BoPh5s2bzo6dkpIyf/78bdu2rVixIiQkRIVdtAUFBRMmTFC+qJiYmK6MLisqKkJCQkpKSsLCwlijPQtrlLqksLDwmWeeUXohMjLywIEDdmezWCwNDQ0VFRWd1aji6tWr3donIOn06dOtra3K66CgoCNHjjjvZ125ciU5OVn5op544gmz2dzFBZOSkt555x0hBGu0x2GN0i+oqqpqP0Y0YMCArhwj+sUaVXz//fczZsxQGmfKlClObTeLxVJRUfHBBx+Eh4c3NjY640c0NDRkZGR4eXkpuzsyMjJu377dxWV37doVHh6uzM8a7XFYo9Qp5YwlPz8/AB4eHqmpqfX19V1ZsIs1KoSw2WxmszkkJKT9fKny8nK51PadPn06MjLSz8/PZDI5fOU2my0nJyc0NFT5LZKSkrq1B7alpWXw4MGzZ882GAwGgyEgICAlJeX69esOz0lOwhol+3Jzc8PCwtpPciwpKen6sl2vUUXHcZyPj0+3xnHdUlNT89hjj+3bt8+B6/zhhx9mzpypfFGTJ08+fPhwd9fQ3Ny8uYPg4OAVK1ZUV1c7MCQ5FWuU7lVWVhYXF6f0wqhRo/Ly8rq7hu7WqOLcuXOJiYkA9HqvmJhbO3Z0dwWdamhoaH/97LPPbt682SGrvXpVvPKKmDz5BWUPb3Z2tkP28HKjvsdhjdJdtbW1qamper0eQGBg4INPcuxMXl7eRx995O7unpOTU1RU1N3Fv/rqq9/9Lh8QgJg9WxQXd3cFdkRFRS1fvnz37t3r16/38/M7f/685AqbmsSaNcLHRwBi7Ngf09PTOza1pEOHDnVx5wk9InRCCIlLSclFWK3WrKyslStXVlVV6fX6BQsWrFmzJjg4+CFW9d5779XW1iqvZ8yYMWfOnO6HQVYWVq5EVRX0eixYgDVr8FBZ2tTW1m7cuLGkpCQ4OHjRokXjx49/+HUBO3di6VJcvAgACQn4xz8wfLjM+qjn07rHSXt79+5tb5ZZs2adPHlS60RCCFFbK1JThV4vABEYKIxG0f2RsYOVlIhf/1ooI+UxY8Tu3RrnoUcEa7RXO3fuXFJSklKgI0aMyMnJ0TrRvcrKRFxcW3ONGiW6v5/WMWpqRGqqcHcXgAgKEiaT6H03BqBOsUZ7kubm5lWrVs2YMSM+Pl7ycLNqB8cdIjdXhIW1lWlsrOjOWQOyWlqEySQCAgQgPDzEokWiqkq9n049Amu0J1mxYsW0adMKCwu3bt3q6+v7cJeHq3aqpmMpdebn11ZnqalChcMwe/aIcePu1rdDjneR62GN9hhWqzUoKKj9tMT2awe7Rc0Lh5yhquruxvWAAU7cuD5zRiQktBVoeLh49PZ20COENdpjXLt2DcCdO3eUt0ajcf78+V1fXOXL2J2qsFA880xbx0VGik6u739IN24Ig0F4egpABAQIo1H8/JUT2cdnMfUYt27d0uv1np6eylsfH5/6+vquLHj79u1169aNGTPmk08+8fLyMhgMpaWlKSkpSqX2RJGROHgQubkYMgRFRYiKwksv4dKl/5vn0CGkpSEtDaWldtZw5gzS0rB+/f99aLMhOxujR2PdOlgsSE5GWRkMBvz8lRN1Qusep666ceMGgNraWuXt6tWrX3vttV9cSqtbfKqjsVEYjaJfPwEIb29hMIj2e5++/37bcHXaNPHz3Z3uyssTgAgLu/vJvn0iIqJtkehocfy4Sr8CuQCORnuMgICA8ePHf/nllwCEEPn5+e33k7fr+PHj0dHRiYmJly5dUm44v3PnzmHDhqmVVw3e3jAYUFqK5GTcvo116zB6NLKz0fGakqNH8fHHD1rJ1atISUFMDE6exOOPw2zGN99g4kRnZycXonWPUzd88cUXQUFB6enpiYmJ48aN6+wUpd75+KP9+8WkSQIQoaGiqaltNDpihNDpRFCQuOfG8x1Ho5mZAhA+PuKvfxVNTZpkp56ND1juSZQH9RQUFLz88stJSUnKWZ8d3fP83jfffDMzM9Pf31+TtCqLisKxY/j4YwwahL592z6cOBGRkcjJwZ//DLPZ/oLLl6O6Gn/5CwYPVi0suRate5wchs/v7UgZjf72t+LyZeHtLQCxd+/dqffvGyV6aNw36grOnj2bkJAwe/bskpIS5fm9HSu1lwsNxVtvAcDixbhzR+s05IpYoz3e22+//eSTT+bl5QUGBppMpuLi4vZHopMiLQ2hoTh79t4znIgcgjXa4/n4+LS2tiYnJ5eVlS1dutTDw0PrRI8cHx+8/z4AvPsuzpzROg25HNZoj7ds2bITJ05kZ2crj0Qnu+bORVwcmpuxbJnWUcjlsEZ7PC8vL8n7EPcSGzbAywt5ecjP1zoKuRbWKPUWYWFISwOAZctgtWqdhlwIa5R6kbQ0jBqF0lL8619aRyEXwhqlXsTTEyYTAOzcqXUUciGsUepdXngB8+ZpHYJcC2uUeh2TCf36aR2CXAivqSfXNHIkkpIwfbqdSY8/jo0bkZeHQYNUj0WuiM+pJyKSwo16ck1//zvi4rB5s/2p5eWIi0NcHJqa1I1Frogb9eSaTp3C7t0YOdL+1MZG7N4NgCeQkgNwNEpEJIU1SkQkhTVKRCSFNUpEJIU1SkQkhUfqyZVt24ZvvrHzeUuL6lHIdbFGyZXV1KCmRusQ5Oq4UU+u7I9/hNVq5w8fJUIOxNEouTKdDu7udj534/iBHId/m4iIpLBGiYiksEaJiKSwRomIpLBGiYik8Eg9uabf/AbDhtm/+z2AoCBkZACAp6eaocg18e73RERSuFFPRCSFNUpEJIU1SkQkhTVKRCSFNUpEJIU1SkQk5X/yoDvUNtY1jAAAAIV6VFh0cmRraXRQS0wgcmRraXQgMjAyNC4wMy4xAAB4nHu/b+09BiAQAGJGBghgBWIWIG5gZGNIAIkzszFoAGlmFg4IzcTOkAGm4Qq4gZoZmTSYGJkVmFkYnECmiLshGcnAuujki/0MDA77QZxGVlYQrYpgH7AHsZ063O2R1DgAxZeC2GIAQOoS7q0w1coAAADNelRYdE1PTCByZGtpdCAyMDI0LjAzLjEAAHicfVDbCsMgDH33K/IDlcRLWx/bWsYYVdi6/cPe9/8sUpyW0eYoxOTkchSQ7O5v7w/8THkhAPDkOOfgpRFRLJAcGOfLNcC0DmOOTPEZ1gdYMFzB2DOHNS45QjBBQ7KnvrUaGpSd5c5cIhE3JzNVYqK05FqyVf6PqCFWRJL2iGggcLTrnFHnoy2PVhK73rSUiIej5+B36ja9Ywy+6E1QRVSCLqtzDkxZkPjaunvdK73zr7MvvrR1V++brqz7AAAAgXpUWHRTTUlMRVMgcmRraXQgMjAyNC4wMy4xAAB4nGXJMQ6AIAyF4as4alKbFihgjRO7egUHR4Mxjh5enBic3v/lpTm1qZuW5mkNUojOA8FIGMLgjIW+lMBYRnjwLN/ZM0aOXv4no0AH252P9cqnEn4553tHp1xh1VSw2gpSV2GeF2oaKdIxJSA6AAAAAElFTkSuQmCC",
"text/plain": [
"<rdkit.Chem.rdchem.Mol at 0x7437ceb0e880>"
]
},
"execution_count": 2,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"mol = Chem.MolFromSmiles('CC(=O)N(C)')\n",
"Chem.MolToSmiles(mol)\n",
"mol"
]
},
{
"cell_type": "code",
"execution_count": 3,
"id": "9ffae000-45eb-413d-afee-1ee4a6089b08",
"metadata": {},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAcIAAACWCAIAAADCEh9HAAAABmJLR0QA/wD/AP+gvaeTAAAPYUlEQVR4nO3de0xUVwIG8G9GUEAEEfGBillHFEK1WB+sD6oVdn1Ak24t7MYWam2FWO3o9iHV1mXbBJVsY7BhG0fd7CI2GmzTBrNoFkjrWvEF1EoF0eq24mIRkQKCMDBz9o9LsSqy6Jm5l7l+v7+GuWfufGL8PPfeM3cMQggQEdHDMmodgIjItbFGiYiksEaJiKSwRomIpLBGiYiksEaJiKS4aR2AXN5//4udO2E04k9/6n7Anj347jtERSEyUt1kRKowcN0oSTp5EhERMBphs3U/YMEC/OtfSEvDhg3qJiNSBQ/qiYiksEaJiKSwRomIpLBGiYiksEaJiKRwwRM5hhBYsaL7TeXl6kYhUhdrlBxDCOzapXUIIi2wRskxjEYcOdL9pj/+ESdPqpuGSEWsUXKYWbO6f37wYHVzEKmLl5iIiKSwRomIpLBGiYiksEaJiKSwRomIpPBKPckaPBgLF6Jfv/sOmD4dRiNMJhUzEamI9xslIpLCg3qSZbMhPh7x8fjzn7sfsHo14uNx+rSqqYhUw9koyerogLt75+P8fERH3z1g3Dj85z84eBALF6ocjUgNnI2SI61cidZWrUMQqYs1Sg7zxBP47jts3qx1DiJ1sUbJYTZtgtGI9HRUVmodhUhFrFFymPBwLF+OtjYkJ4On3OnRwRolR0pLw+DBOHwYH3+sdRQitbBGyZGGDetc9vTGG7hxQ+Mwfd+pU6e0jkAOwBolB1u1CuHhuHYNb7+tdZQ+rKqqKjExccaMGZ9++qnWWUgWa5QczM0NFguMRvztb7zpfTdu3ry5YcOGCRMmZGdne3t719XVaZ2IZLFGyfFmzMCyZbDbsXYtrzXdJoTYv39/WFjY5s2b29ra4uLizp49m5SUpHUuksUaJadIT4e/P44dQ06O1lH6huLi4jlz5sTHx1++fHnatGlfffVVTk5OUFCQ1rnIAVij9JDOnsW6dffdOnQoNm0CgA0bYLXeft5iQW6u07P1KdXV1cnJyREREUVFRYGBgRaL5cSJE7Pu88VVp06d2r9//w1ennMtgugB1dWJ1auFm5sARE6OaG8XgADEjz/eMcxmEzNndm4CxMGDoqpKDBwoAPGb34hvv9UovYra2toyMjIGDRoEoH///mazubGxsYfxV65cGTVq1JAhQ06cOKFaSJLH2Sg9gI4O7NiBkBBkZsJuR0IC5s2772CjEZmZd9yHdMQIbN2KgADk5yM8HMnJqK11fmiNHDhwIDQ0dO3atU1NTbGxsRUVFdu2bVMq9X5Wr1795ptv+vn5qRaSHEPrHieXUVgoJk/unFo+9ZT45pvO5+83G1WsWnV7Nqq4cUOYzZ2TWT8/kZEh2ttV+iOoo7y8fMGCBcq/r9DQ0EOHDvXmVfv27Zs8eXJ7e7vJZOJs1LWwRun/u3BBxMV1tmFQkMjKumOrzSZSUkRKimhq6ua19fXi7bdFSoqorLzj+YoKsWhR5z4nThT//KcT86umrq7ObDb369cPwJAhQzIyMjo6Onr5wpEjRx47dkwIwRp1OaxR6snNmyI1VXh4CEAMHChSU8WtW47cf26uMJk6yzQ6WpSXO3LnarJarRaLZejQoQDc3d2TkpJqa2t7//Jly5aZzWblMWvU5bBGqXt2u8jKEiNGCEAYDCIhQVy96pQ3slpFRobw8RGAcHcXZrNoaHDKGzlPfn5+WFiYchQfHR1dVlb2QC9vbm42mUzjfubu7j5q1Khy1/0v5dHDGqVunDx5+yL79OmiqMjp71hdLZKSRL9+AhBDh4qMDNG7o2GNnT9/Pi4uTinQ4ODgnJwc+X1yNupyWKN0hytXREKCMBgEIAIDRVaWsNvVe/eSEjFnTmd9P/GEOHJEvbd+UPX19SkpKQMGDADg7e2dmpra2trqkD0nJiaeO3fOIbsidbBGqVNLi9iyRQwaJADh6SlSUkSPaxydKDdXjB3bWaaxseL777WJcT82my0rK2v48OEAjEZjQkLCj92uUaBHBmuUhBAiN1f86le3m+vSJY3zNDeLLVuEt7cAhJfXfZcBqO/LL78MDw9XjuIjIiKOHz+udSLSHmv0UVdaKp58srNAp0wRhw/L7vDChQvvv//+O++8U1hYKLmrqqrbZxhGj5Y6w9DY2JiZmSkTqaqqKiEhwWAwABg9enRWVpZdzfMd1IexRl1MfX19QUFBcXGx/L/hmpqaV181+/vfAMSwYWLHDmGzycY7ceJEcHDwpk2bPvroo8DAwOzsbNk9CnH4sAgP7yz6F174a2lp6YPuYfv27SaTKTg4eM2aNQ8RoLm5OTU11dPTE4CXl1dqampLS8tD7If0ijXqSoqKivz9/WNiYkJCQmJiYnq5tPtebW1tH3zwga+vL4Cnnlrzxhvip58ck7C2traiokJ5vHHjxueff94hu+3oEDt2iOjoH5TTkStWrKipqen9y2tqajo6OjZu3PigNWq323NycsaOHQvAYDDExcV939fO1FIfwBp1JXPnzk1PTxdCNDc3jx8//rPPPnuIndy1yPFbJ9wjxG63nzt3btKkSXv27HHgbpuamlJTUx/64viD1mhxcfHs2bOVX9TUqVOP9OV1A6Qp1qjLaGpqAlBdXa38uG7duhUrVjzQHiorK2NiYpRemDBhwoEDB5wQUwghli5dGhAQEBkZ+ZOjZrm/8NBLNXtfo9XV1UlJSUajEcDIkSMtFstDT/zpUcAadRnnz593c3PrOiX64YcfxsTE9PK1yiLH/v37Axg8ePCWLVva2tqcllQIIWw228qVK5cuXeqk/RcUFDz22GNKmUZFRfXmg0O9qVGr1ZqRkeHj46N8ptNsNje43GeqSHW8UZ7L8PDwsNlsHR0dyo+tra3KRY+e2e323bt3T5w4MT09vaOjIyEhobKysqtSHU4I0d7eDsBoNM6bN6+srMwZ7wIgKirq66+/Vj7GXlhYOGXKlOTk5OvXr8vss+vWdo2NjV23tlMqlagnWvc49ZbVavX09Oyadr388svr1q3r+SVffPHF448/rvxFz5s37/Tp084OuXfv3jlz5uTl5RUWFk6dOvWtt95y9jvee1Ol9nvuu3f16tWcnJwlS5YsXrw4Jyfn3nsnV1RULFq0SPlFhYSE5OXlOTs26Qlr1JUsX748Li7u1q1bpaWlvr6+PRzJXr58OSEhQemFMWPGZN11bzunsdvtWVlZS5Ysefrpp7dt23ZvozlJeXn5woULu3rwYNfNTYUQQpSVlSX9wtVf3GRFaWE3NzcAfn5+3bYwUc9Yo66ktrb2mWee8fX1DQoK2rlzZ7djlEWOHh4e+HmR4y3H3tuuD8vNzR03bpxSprGxsRcvXuxhcHt7u8ViCQgIAODm5paUlHTt2jXVopKesEb1Q1nkqHzZpLLI8YcfftA6lNqUrz9STmgqX3/U7TWigoKCSZMmKYU7f/78M2fOqB+VdIM1qhOnTp3q+rLJadOmHT16VOtEWvrliqWAgACLxWL7+RNapaWlXeulxo8f75Bb29EjziCEcOIFLHK+6urq9957b9euXXa7PTAwMDU19ZVXXlEa5BEXFxf3ySefKI+nTp364osv7tu3r6ioqGtAS0tLb1Y7EPXMTesAJGvZsmX5+fkeHh6vv/76+vXrvb29tU7UVwwbNgxAYmJiQUFBSUlJSUmJ8vysWbOUbz1SPhBFJIlzFpeXlpb23HPPlZeXp6WlsUPvNWPGjMrKyujoaC8vryFDhmRnZx89elTrUKQrnI26vOnTp+/fv1/rFH2at7d3fn6+1ilItzgbJSKSwholIpLCGiUiksIaJSKSwholIpLCGiUiksIaJSKSwholIpLCGiUiksIaJSKSwholIpLCGiUiksIaJSKSwholIpLCGiUiksIaJSKSwholIpLCGiUiksIaJSKSwholIpLCr7Qj3Xpt2LBXIyMHdvdtqWVPPgm73WgwqJ+K9Ic1SroVcu0ajhzB739/76awf/8bQkAIsElJGg/qiYiksEaJiKSwRomIpLBGiYiksEaJiKSwRomIpLBGiYiksEaJiKSwRomIpLBGiYiksEaJiKSwRomIpLBGiYiksEaJiKSwRomIpLBGiYiksEaJiKSwRomIpLBGiYiksEaJiKSwRomIpLBGiYiksEaJiKQYhBBaZyByjkuXcPUqTCaMGHH3pqIiCIHZs7WIRXrDGiUikuKmdQAiZxICeXn4/HOUleH6dXh5ISgIUVFITIS/v9bhSCc4GyX9qq7GkiU4frzzxwED0NbW+djXF5mZeOEFraKRnvASE+lUfT0iI3H8OEaPxq5dqK1FayusVhw5gpgYNDQgMRF79midkvSAs1HSqZdewj/+gbFjUVSEwMA7NgmBlSthscDHB5WV3VyAInoQnI2SHlVX4+OPAWDbtrs7FIDBgIwMjBmDxkZs365+OtIZ1ijp0aFDaG9HQABiY7sf4OGBxEQAOHBAzVykS6xR0qPiYgCYPh39+t13zMyZAHDmDKxWlVKRTrFGSY9qagBg1KiexoweDQAdHairUyMS6RdrlPTo1i0A8PDoaYynZ+eDlhan5yFdY42SHvn4AMDNmz2NaWrqfODr6/Q8pGusUdIj5YD90qWexly8CAADB8LPT41IpF+sUdKjiAgAKCnp6YD98GHg/12GIuoF1ijp0cKF8PHBzZvYvbv7AXV12LsXAP7wBzVzkS6xRkmPBg2C2QwA69fj9Om7t1qteOkl1Ndj7Fh+rJ7k8cOgpFNWK+bPx9GjGDgQa9bgd7/DmDFoasLx4/jLX3DmDLy8cOgQIiO1DkoujzVK+tXcjFWrkJ0Nu/3uTSEh+Pvf8etfaxGL9IY1SnpXXo7PP8c33+DGDXh6IigIv/0tFi+GG2+2S47BGiUiksJLTEREUlijRERSWKNERFJYo0REUlijRERSuOaD9K+lpSUzM7O0tHT48OFms9lkMmmdiHSFs1HSv2efffbcuXPJyck+Pj6RkZFNXbfII3IErhsl/btw4YLJZDIajQACAwP37t07d+5crUORfvCgnvQvODhYedDe3t7S0uLv769tHtIZHtTTI2Tnzp2hoaFhYWFaByFd4WyUHhV5eXnp6en5+fkGg0HrLKQrrFHSPyHE1q1bd+zYcfDgwQkTJmgdh/SGl5hI54QQ8fHxeXl57777rnJWdObMmZMmTdI6F+kHZ6Okc3a73WQyvfbaaw0NDQ0NDQBCQ0O1DkW6wtkoEZEUXqknIpLCGiUiksIaJSKSwholIpLCGiUiksIaJSKS8j/6chZethfNCAAAAIx6VFh0cmRraXRQS0wgcmRraXQgMjAyNC4wMy4xAAB4nHu/b+09BiAQAGJGBghgBWIWIG5gZGNIAIkzszFoAGlmFg4IzcTOkAGm4Qq4WRgZGJk0mBiZFZgZGJxApoi7IRnJwNrIynogLU1NDcQBsu3PnjmzBMZmYDiwH8R26nAH0g72IPaiky/sYWwxAFQ/FIVT5CEKAAAAyXpUWHRNT0wgcmRraXQgMjAyNC4wMy4xAAB4nH2QzQ7CIAzH7zxFX0DS8jHGcRuLMWaQ6PQdvPv+sWgQdthaSNryo/SPgGy3cH294W8qCAGAB8t7D0+NiGKBHMA4ny8RpnUYS2VKj7jewYLhG+xbcljTUioEE5yURNebjgAlfq0JCqgYRGnJd2ThdABqSA1I0u6BBmJu5Jw3SvO5szug5adJ9tR39pCbY9ho+6kdUwxVLXE3VTVl13Vy4tTU+SjvtnvbK+flzzkWHwL4V5X4iK5NAAAAgHpUWHRTTUlMRVMgcmRraXQgMjAyNC4wMy4xAAB4nF3JMQqAMAxA0as4KqQhaZu2Rpy6q1dwcJSKOHp4dSq4/cfPU25zN87N3RqLFJMPQDAYwhh7bx28ITAQCveB5XuMiVOQ/zGMAh2sV9mXsxxK+OVUrg1JucKprWB1FV59hb0fLnIpeHuLtKUAAAAASUVORK5CYII=",
"text/plain": [
"<rdkit.Chem.rdchem.Mol at 0x7437ceb0e810>"
]
},
"execution_count": 3,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"newmol1 = rdmolops.RenumberAtoms(mol, [4, 1, 2, 3, 0])\n",
"Chem.MolToSmiles(newmol1)\n",
"newmol1"
]
},
{
"cell_type": "code",
"execution_count": 4,
"id": "f0be7d1a-45f4-4fc6-b239-351cd19c63df",
"metadata": {},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAcIAAACWCAIAAADCEh9HAAAABmJLR0QA/wD/AP+gvaeTAAAPKUlEQVR4nO3dbVBU1wEG4HeXRQERBAQVBaeuWgjR0ChSP1CjtGrhR5oK7Y9AiImQRLOxiQORyXSbH8SPZhzM2KQrdlpqMjpoJhkc0Sl0jDXiFzBGI4pfNWrXD4oUEZAF9vTHJVgVCHB272Uv7/Nr2T1795UZXs+999y7BiEEiIhooIxaByAi8mysUSIiKaxRIiIprFEiIimsUSIiKaxRIiIpJq0DkMf7979RUACjEb/7XfcDPv0Uly5h8WIkJKibjEgVBq4bJUknTiA+HkYjOjq6H7BkCf7+d+TlITdX3WREquBOPRGRFNYoEZEU1igRkRTWKBGRFNYoEZEULngi1xACK1d2/1J1tbpRiNTFGiXXEALbt2sdgkgLrFFyDaMRhw93/9Jvf4sTJ9RNQ6Qi1ii5zJw53T8/apS6OYjUxVNMRERSWKNERFJYo0REUlijRERSWKNERFJ4pp5kjRqFpUvh5dXjgLg4GI0wm1XMRKQi3m+UiEgKd+pJVkcHUlORmorf/777AatXIzUVp06pmopINZyNkqz2dnh7dz4uLUVi4uMDJk3Cv/6F/fuxdKnK0YjUwNkoudLrr+PBA61DEKmLNUou8+yzuHQJ69drnYNIXaxRcpkPPoDRiI0bUVOjdRQiFbFGyWViY7FiBVpbkZUFHnKnoYM1Sq6Ul4dRo3DoED77TOsoRGphjZIrhYV1Lnt65x3cvatxmMHv5MmTWkcgF2CNkoutWoXYWNy5g3ff1TrKIHb9+vX09PRZs2Z9/vnnWmchWaxRcjGTCTYbjEb8+c+86X037t+/n5ubO3Xq1B07dvj7+9fV1WmdiGSxRsn1Zs1CRgacTqxZw3NNDwkhdu/eHRMTs379+tbW1pSUlLNnz2ZmZmqdi2SxRsktNm5ESAiOHkVRkdZRBoeKiop58+alpqZeu3Zt5syZX3/9dVFRUWRkpNa5yAVYozRAZ88iO7vHV0ePxgcfAEBuLhyOh8/bbCgudnu2QcVut2dlZcXHx5eXl4eHh9tstuPHj8/p4YurTp48uXv37rs8PedZBFE/1dWJ1auFySQAUVQk2toEIABx69Yjwzo6xOzZnS8BYv9+cf26GDFCAOJnPxPffqtRehW1trbm5+ePHDkSwLBhwywWy71793oZf+PGjfHjxwcHBx8/fly1kCSPs1Hqh/Z2bNuGqChs3QqnE2lpWLiwx8FGI7ZufeQ+pGPHYvNmhIaitBSxscjKQm2t+0NrZO/evdHR0WvWrGlsbExOTj537tyWLVuUSu3J6tWr165dGxQUpFpIcg2te5w8xj/+IaZP75xaPvec+Oabzud7mo0qVq16OBtV3L0rLJbOyWxQkMjPF21tKv0T1FFdXb1kyRLl7ys6OvrAgQN9edeuXbumT5/e1tZmNps5G/UsrFH6YRcvipSUzjaMjBSFhY+82tEhcnJETo5obOzmvfX14t13RU6OqKl55Plz58SyZZ3b/PGPxb59bsyvmrq6OovF4uXlBSA4ODg/P7+9vb2Pbxw3btzRo0eFEKxRj8Mapd7cvy+sVuHjIwAxYoSwWkVLiyu3X1wszObOMk1MFNXVrty4mhwOh81mGz16NABvb+/MzMza2tq+vz0jI8NisSiPWaMehzVK3XM6RWGhGDtWAMJgEGlp4uZNt3yQwyHy80VAgACEt7ewWERDg1s+yH1KS0tjYmKUvfjExMQzZ8706+1NTU1ms3nS97y9vcePH1/tuf+lDD2sUerGiRMPT7LHxYnycrd/ot0uMjOFl5cAxOjRIj9f9G1vWGMXLlxISUlRCnTKlClFRUXy2+Rs1OOwRukRN26ItDRhMAhAhIeLwkLhdKr36ZWVYt68zvp+9llx+LB6H91f9fX1OTk5w4cPB+Dv72+1Wh88eOCSLaenp58/f94lmyJ1sEapU3Oz2LBBjBwpAOHrK3JyRK9rHN2ouFhMnNhZpsnJ4upVbWL0pKOjo7CwcMyYMQCMRmNaWtqtbtco0JDBGiUhhCguFj/60cPmunJF4zxNTWLDBuHvLwDh59fjMgD1ffXVV7GxscpefHx8/LFjx7RORNpjjQ51VVVi/vzOAv3JT8ShQy7Y5sGDB/Py8uS3c/36wyMMEyYM/AiD0+ncs2dPdnb2pk2bbt++PdAw19PS0gwGA4AJEyYUFhY61TzeQYMYa9TD1NfXl5WVVVRUyP8N3759+403LCEhdwERFia2bRMdHbLxrl69GhcXFxcXZzKZZLf1vUOHRGxsZ9G/+OIfq6qq+ruF1157LTk52WazrVq1KiIiovcrMp/U1NRktVp9fX0B+Pn5Wa3W5ubm/mYgHWONepLy8vKQkJCkpKSoqKikpKQ+Lu1+Umtr64cffhgYGAjguefeeucd8d//uiahw+Gora29deuWC2tUCNHeLrZtE4mJ3ymHI1euXNmvSeXRo0e7flfBwcHKKve+cDqdRUVFEydOBGAwGFJSUq4OtiO1NAiwRj3JggULNm7cKIRoamqaPHnyF198MYCNPLbI8Vs33CPE5TWqaGxstFqtAz453tLSUlhYGBkZ2cfZaEVFxdy5c5Vf1IwZMw4P5nUDpCnWqMdobGwEYLfblR+zs7NXrlzZry3U1NQkJSUpvTB16tS9e/e6IaYQbqtRxcCWap45c+app57y8fH5+OOPf3Cw3W7PzMw0Go0Axo0bZ7PZBjzxp6GANeoxLly4YDKZug6JfvTRR0lJSX18r7LIcdiwYQBGjRq1YcOG1tZWtyV1b40qysrKnn76aaVMFy9e3McLh+x2e1hY2JEjR3oa4HA48vPzAwIClGs6LRZLg8ddU0WqY416jGvXrhkMBofDofy4adOm5cuX/+C7lEWOYWFhXYscB3yquu9UqFEhRFtbW9dl7CaTqZfL2P///4z58+fbbLZuhxUXF5vNZqWak5OTL1265JbcpDusUY/hcDh8fX27pl2vvPJKdnZ27285ePDgM888o/TCwoULT5065f6YYt++fQUFBV5eXkVFRQM4q95fT95Uqe2J++4tXLjQarWWl5d/8skn/v7+NY/dbEqIc+fOLVu2TPlFRUVFlZSUuDs26Qlr1JOsWLEiJSWlpaWlqqoqMDCwlz3Za9eupaWlKb0QERFR+Ni97dwpNzc383s7d+5U50Orq6uXLl3a1YP7u25uKoQQ4ubNm2vXrl22bFlGRsZjC+aVFjaZTACCgoK6bWGi3rFGPUltbe3zzz8fGBgYGRlZUFDQ7RhlkaOPj0/XIscW197bbhArLi6eNGlS11755cuXexmsHBMIDQ3tOiZw584d1aKSnrBG9UNZ5Kh82aSyyPG7777TOpTalK8/Us4RKV9/1O05orKysmnTpimFu2jRotOnT6sflXSDNaoTJ0+e7PqyyZkzZ/ZyMnoo+P8VS6GhoTabreP7K7Sqqqq61ktNnjzZJbe2oyHOIITowzc20eBlt9vff//97du3O53O8PBwq9X66quvKg0yxKWkpOzZs0d5PGPGjJdeemnXrl3l5eVdA5qbm5VLPIlkmLQOQLIyMjJKS0t9fHzefvvtdevW+fv7a51osFCWeaWnp5eVlVVWVlZWVirPz5kzR7keVLkgikgS5yweLy8vb/ny5dXV1Xl5eezQJ82aNaumpiYxMdHPzy84OHjHjh1HjhzROhTpCmejHi8uLm737t1apxjU/P39S0tLtU5BusXZKBGRFNYoEZEU1igRkRTWKBGRFNYoEZEU1igRkRTWKBGRFNYoEZEU1igRkRTWKBGRFNYoEZEU1igRkRTWKBGRFNYoEZEU1igRkRTWKBGRFNYoEZEU1igRkRTWKBGRFNYoEZEUfqUd6dabYWFvJCSM6O7bUs/Mnw+n02gwqJ+K9Ic1SroVdecODh/Gr3/95Esx//wnhIAQYJOSNO7UExFJYY0SEUlhjRIRSWGNEhFJYY0SEUlhjRIRSWGNEhFJYY0SEUlhjRIRSWGNEhFJYY0SEUlhjRIRSWGNEhFJYY0SEUlhjRIRSWGNEhFJYY0SEUlhjRIRSWGNEhFJYY0SEUlhjRIRSWGNEhFJYY0SEUkxCCG0zkDkHleu4OZNmM0YO/bxl8rLIQTmztUiFukNa5SISIpJ6wBE7iQESkrw5Zc4cwb/+Q/8/BAZicWLkZ6OkBCtw5FOcDZK+mW341e/wrFjnT8OH47W1s7HgYHYuhUvvqhVNNITnmIinaqvR0ICjh3DhAnYvh21tXjwAA4HDh9GUhIaGpCejk8/1Tol6QFno6RTL7+Mv/4VEyeivBzh4Y+8JARefx02GwICUFPTzQkoov7gbJT0yG7HZ58BwJYtj3coAIMB+fmIiMC9e/jTn9RPRzrDGiU9OnAAbW0IDUVycvcDfHyQng4Ae/eqmYt0iTVKelRRAQBxcfDy6nHM7NkAcPo0HA6VUpFOsUZJj27fBoDx43sbM2ECALS3o65OjUikX6xR0qOWFgDw8eltjK9v54PmZrfnIV1jjZIeBQQAwP37vY1pbOx8EBjo9jyka6xR0iNlh/3Kld7GXL4MACNGIChIjUikX6xR0qP4eACorOxth/3QIeCHTkMR9QFrlPRo6VIEBOD+ffztb90PqKvDzp0A8JvfqJmLdIk1Sno0ciQsFgBYtw6nTj3+qsOBl19GfT0mTuRl9SSPF4OSTjkcWLQIR45gxAi89RZ++UtERKCxEceO4Q9/wOnT8PPDgQNISNA6KHk81ijpV1MTVq3Cjh1wOh9/KSoKf/kLfvpTLWKR3rBGSe+qq/Hll/jmG9y9C19fREbi5z/HL34BE2+2S67BGiUiksJTTEREUlijRERSWKNERFJYo0REUlijRERSuOaD9K+5uXnr1q1VVVVjxoyxWCxms1nrRKQrnI2S/r3wwgvnz5/PysoKCAhISEho7LpFHpErcN0o6d/FixfNZrPRaAQQHh6+c+fOBQsWaB2K9IM79aR/U6ZMUR60tbU1NzeHhIRom4d0hjv1NIQUFBRER0fHxMRoHYR0hbNRGipKSko2btxYWlpqMBi0zkK6whol/RNCbN68edu2bfv37586darWcUhveIqJdE4IkZqaWlJS8t577ylHRWfPnj1t2jStc5F+cDZKOud0Os1m85tvvtnQ0NDQ0AAgOjpa61CkK5yNEhFJ4Zl6IiIprFEiIimsUSIiKaxRIiIprFEiIimsUSIiKf8DYbnsz7BxsiAAAACMelRYdHJka2l0UEtMIHJka2l0IDIwMjQuMDMuMQAAeJx7v2/tPQYgEABiRgYIYAViFiBuYGRjSACJM7MzZABpZiZGDgYNMIMNQrPAFHCzMANFNZiYGRWApjiBTBF3QzKSgbWRlfVAWpqaGojj1OG+n4HBwR7EBooD6QP7YeyzZ84sAbEXnXxhD1MjBgA0PxSH0s1g/gAAAMx6VFh0TU9MIHJka2l0IDIwMjQuMDMuMQAAeJx9UEkOwjAMvOcV84FGcZamOdKmQgg1kaDwB+78X7igkPbQOjmM7fEyFljsFq+vN/6moxCAOvghBDyNUkpMWAD68XxJGOZTXyJDfqT5DgfLFfy2zNOcpxIhDGi0VL6zLUFJ9bUVKESNhEZJ74PVhvPe7RANMkcdhZYcGpJur6Pl0ZV4MNoxkWRHXesOJ48pbrT91PY5xaqWlqNUTZq3MHVzYlfX/YjztO6+7rX45eaMxQcM3FeXkFuu/gAAAIB6VFh0U01JTEVTIHJka2l0IDIwMjQuMDMuMQAAeJxdyTEKgDAMQNGrOCqkIWmbtkacuqtXcHCUijh6eHUquP3Hz1NuczfOzd0aixSTD0AwGMIYe28dvCEwEAr3geV7jIlTkP8xjAIdrFfZl7McSvjlVK4NSbmC1VY4dRVefYW9Hy5YKXiDEOREAAAAAElFTkSuQmCC",
"text/plain": [
"<rdkit.Chem.rdchem.Mol at 0x7437ceb0e8f0>"
]
},
"execution_count": 4,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"newmol2 = rdmolops.RenumberAtoms(mol, [4, 3, 2, 1, 0])\n",
"Chem.MolToSmiles(newmol2)\n",
"newmol2"
]
},
{
"cell_type": "code",
"execution_count": 5,
"id": "cfcbdbd4-0158-4a17-9d86-2a3a55c48455",
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"[4,3,1,0,2,] CNC(C)=O\n",
"[0,3,1,4,2,] CNC(C)=O\n",
"[0,1,3,4,2,] CNC(C)=O\n"
]
}
],
"source": [
"print(mol.GetProp('_smilesAtomOutputOrder'), Chem.MolToSmiles(mol)) \n",
"print(newmol1.GetProp('_smilesAtomOutputOrder'), Chem.MolToSmiles(newmol1))\n",
"print(newmol2.GetProp('_smilesAtomOutputOrder'), Chem.MolToSmiles(newmol2))"
]
},
{
"cell_type": "code",
"execution_count": null,
"id": "7e7a2970-88a4-47f7-a259-ccf2d1d0abf5",
"metadata": {},
"outputs": [],
"source": []
}
],
"metadata": {
"kernelspec": {
"display_name": "Python 3 (ipykernel)",
"language": "python",
"name": "python3"
},
"language_info": {
"codemirror_mode": {
"name": "ipython",
"version": 3
},
"file_extension": ".py",
"mimetype": "text/x-python",
"name": "python",
"nbconvert_exporter": "python",
"pygments_lexer": "ipython3",
"version": "3.10.13"
}
},
"nbformat": 4,
"nbformat_minor": 5
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment