View on GitHub

Download Notebook

Run in Colab

Visualizing Ligand Networks

When constructing a ligand network, it is very helpful to visualize it, both to quickly check that it makes sense and to allow rapid iteration of settings.

Setup

%matplotlib inline

import rdkit.Chem

import openfe

This cookbook assumes you’ve already loaded a LigandNetwork. For more information, see Generate a Ligand Network Automatically:

ligand_network = openfe.ligand_network_planning.generate_minimal_spanning_network(
    ligands=[
        openfe.SmallMoleculeComponent(mol)
        for mol in rdkit.Chem.SDMolSupplier(
            "assets/somebenzenes.sdf",
            removeHs=False,
        )
    ],
    mappers=[openfe.setup.LomapAtomMapper()],
    scorer=openfe.lomap_scorers.default_lomap_score,
)

Visualize Network Topology

The topology of the network describes how the ligands are transformed into each other. The network itself forms a graph that can be visualized with the plot_atommapping_network function:

from openfe.utils.atommapping_network_plotting import plot_atommapping_network

plot_atommapping_network(ligand_network)

Visualizing Edges

The edges of the network each represent a single transformation. They include not only the identities of the two molecules being mutated between, but also the complete mapping between atoms in those molecules.

Displaying an edge in the notebook will provide a 2D description of the atom mapping:

benzene_to_toluene = [
    edge for edge in ligand_network.edges
    if {edge.componentA.name, edge.componentB.name} == {"benzene", "toluene"}
][0]

benzene_to_toluene

Atoms created or destroyed over the course of the transformation are highlighted in red, while atoms that are transformed are highlighted in blue.

An edge can also be visualized in 3D with the view_mapping_3d function:

from openfe.utils.visualization_3D import view_mapping_3d

view_mapping_3d(benzene_to_toluene)

You appear to be running in JupyterLab (or JavaScript failed to load for some other reason). You need to install the 3dmol extension:
jupyter labextension install jupyterlab_3dmol

<py3Dmol.view at 0x7f35fb4ff4d0>

In the middle, both molecules are overlaid to demonstrate their 3D alignment. Each molecule is also depicted on its own on either side, with colored spheres indicating atoms that are mutated into each other. For instance, one of the benzene’s hydrogen atoms is highlighted in the same color as the toluene’s methyl carbon, indicating that this hydrogen mutates into the methyl carbon. Atoms without spheres vanish over the transformation, rather than being mutated into another atom.

Visualizing Ligand Overlap

Sometimes it is helpful to visualize all the ligands together, for instance to inspect their alignment. This can be done with the view_components_3d function.

from openfe.utils.visualization_3D import view_components_3d

view_components_3d(ligand_network.nodes)

You appear to be running in JupyterLab (or JavaScript failed to load for some other reason). You need to install the 3dmol extension:
jupyter labextension install jupyterlab_3dmol

<py3Dmol.view at 0x7f35f5eca510>