I posted new function about scaffold tree which is implemented in rdkit 2020 03 before.
In previous my post, I showed example to draw scaffold tree with networkx. It could draw the scaffold tree as a network but molecular structures are not shown on the node. For chemist, structure image is important so I tried to draw network with structure image on node.
Fortunately cyjupyter is very useful package to do it, it can call cytoscape.js internally so user can draw network with many options.
OK let’s try. At first import required packages. And set example strcutures.
from rdkit import Chem
from rdkit.Chem import rdBase
from rdkit.Chem import Draw
from rdkit.Chem.Draw import rdMolDraw2D
from rdkit.Chem.Draw import IPythonConsole
import networkx as nx
from networkx.readwrite import cytoscape_data
import cyjupyter
from cyjupyter import Cytoscape
from rdkit.Chem import AllChem
from rdkit.Chem.Scaffolds import rdScaffoldNetwork
from urllib import parse
smiles_lsit = [
"CC1=C(C=C(C=C1)NC(=O)C2=CC=C(C=C2)CN3CCN(CC3)C)NC4=NC=CC(=N4)C5=CN=CC=C5",
"CC1CCN(CC1N(C)C2=NC=NC3=C2C=CN3)C(=O)CC#N"
]
Then create scaffold tree following code is almost as same as previous post.
params = rdScaffoldNetwork.ScaffoldNetworkParams() netwks = rdScaffoldNetwork.CreateScaffoldNetwork([mols[1]], params)
Then define smiles to image function by using MolDraw2DSVG which is implemented in RDKit. The function generates svg from SMIELS. It is used for rendering image on the network node.
def smi2svg(smi):
mol = Chem.MolFromSmiles(smi)
try:
Chem.rdmolops.Kekulize(mol)
except:
pass
drawer = rdMolDraw2D.MolDraw2DSVG(690, 400)
AllChem.Compute2DCoords(mol)
drawer.DrawMolecule(mol)
drawer.FinishDrawing()
svg = drawer.GetDrawingText().replace("svg:", "")
return svg
def smi2image(smi):
svg_string = smi2svg(smi)
impath = 'data:image/svg+xml;charset=utf-8,' + parse.quote(svg_string, safe="")
return impath
Finally draw network with cytoscape.js!
g = nx.graph.Graph()
for idx, node in enumerate(netwks.nodes):
g.add_node(idx, smiles=node, img=smi2image(node), hac=Chem.MolFromSmiles(node).GetNumAtoms())
g.add_edges_from([(e.beginIdx,e.endIdx) for e in netwks.edges])
cy_g = cytoscape_data(g)
stobj=[
{'style': [{'css': {
'background-color': 'blue',
'shape' : 'rectangle',
'width':600,
'height':400,
'border-color': 'rgb(0,0,0)',
'border-opacity': 1.0,
'border-width': 0.0,
'color': '#4579e8',
'background-image':'data(img)',
'background-fit':'contain'
},
'selector': 'node'},
{'css': {
'width': 20.0,
},
'selector': 'edge'}
],
}]
cyobj=Cytoscape(data=cy_g, visual_style=stobj[0]['style'], layout_name='breadthfirst')
cyobj
Now I could get scaffold tree with molecular image. ;)
I tried to change atom font size but failed. I’m trying to do that now.
I uploaded today’s code on my gist. RDKit has many functions for chemoinformatics. I would like to buy a beer to the developers! Thanks!
