Visualizing Protein-Ligand Interactions

Problem

You want to visualize protein-ligand interactions. See example in Figure 1.

Figure 1. Example of visualizing protein-ligand interactions (PDB: 13GS)

Ingredients

  • OEChem TK - cheminformatics toolkit (including OEBio TK)
  • OEDepict TK - molecule depiction toolkit
  • Grapheme TK - molecule and property visualization toolkit

Difficulty level

../_images/chilly7.png ../_images/chilly7.png

Download

Download code

complex2img.py

See also the Usage subsection.

Solution

The depict_complex illustrates how to generate protein-ligand iteration images.

  1. OEInteractionHintContainer object is constructed that stores information about possible interactions between the ligand and the protein.
  2. The interactions are perceived by calling the OEPerceiveInteractionHints function.
  3. The image is then separated into two frames: the ligand and the residues around it is going to rendered into one frame, while the corresponding legend into the other frame.
  4. The active site is then prepared for 2D depiction by invoking the OEPrepareActiveSiteDepiction function.
  5. When the OE2DActiveSiteDisplay object is constructed, residues are positioned around the ligand close to those atoms which they are interacting with. The OERenderActiveSite function generates an image in which each residue cycle are colored based on the types in which they are interact with the ligand. The interactions between the residues and the ligand atoms are marked by using different style of lines. See more details in the Discussion section below.
  6. The legend associated with the active site is rendered by invoking the OEDrawActiveSiteLegend function.
def depict_complex(image, protein, ligand, opts, legend_frame=None):
    """
    :type image: oedepict.OEImageBase
    :type protein: oechem.OEMolBase
    :type ligand: oechem.OEMolBase
    :type opts: oedepict.OE2DMolDisplayOptions
    :type legend_frame: oedepict.OEImageBase
    """

    # perceive interactions

    asite = oechem.OEInteractionHintContainer(protein, ligand)
    if not asite.IsValid():
        oechem.OEThrow.Fatal("Cannot initialize active site!")
    asite.SetTitle(ligand.GetTitle())

    oechem.OEPerceiveInteractionHints(asite)

    # depiction

    oegrapheme.OEPrepareActiveSiteDepiction(asite)
    adisp = oegrapheme.OE2DActiveSiteDisplay(asite, opts)
    oegrapheme.OERenderActiveSite(image, adisp)

    if legend_frame is not None:
        lopts = oegrapheme.OE2DActiveSiteLegendDisplayOptions(12, 1)
        oegrapheme.OEDrawActiveSiteLegend(legend_frame, adisp, lopts)

Usage

Download code

complex2img.py

Usage

The following commands will generate the image shown in Figure 1.

prompt > wget https://files.rcsb.org/download/13gs.pdb
prompt > python3 complex2img.py -complex 13gs.pdb -out 13gs.svg -interactive-legend

Command Line Parameters

Simple parameter list
    -height : Height of output image
    -width : Width of output image

    SplitMolComplex options :
      -covalentligand : Split covalent ligands
      -ligandname : Ligand name

    molecule display options :
      -aromstyle : Aromatic ring display style

    input/output options :
      -complex : Input filename of the protein-ligand complex
      -protein : Input filename of the protein
      -ligand : Input filename of the ligand
      -out : Output filename of the generated image

    active site display options:
      -interactive-legend : Visualize legend on mouse hover (SVG feature)
      -magnify-residue : Scaling factor to magnify residue glyph on mouse hover
                         (SVG feature)

Discussion

Interaction perception

Currently the OEPerceiveInteractionHints function perceives the following interaction types:

Table 1. Interaction types currently available in OEChem TK
name corresponding interaction class corresponding interaction type namespace
cation-pi OECationPiInteractionHint OECationPiInteractionHintType
chelator OEChelatorInteractionHint OEChelatorInteractionHintType
clash OEClashInteractionHint None
contact OEContactInteractionHint None
covalent OECovalentInteractionHint None
halogen bond OEHalogenBondInteractionHint OEHalogenBondInteractionHintType
hydrogen bond OEHBondInteractionHint OEHBondInteractionHintType
salt-bridge OESaltBridgeInteractionHint OESaltBridgeInteractionHintType
stacking (T and Pi) OEStackingInteractionHint OEStackingInteractionHintType

The default geometric parameters used by the OEPerceiveInteractionHints function have been set based on literature data ([Kumar-2002], [Cavallo-2016], [Bissantz-2010], and [Marcou-2007] ). The interaction parameters can be customized by using the OEPerceiveInteractionOptions class.

Interaction depiction

Figure 2. Example of visualizing protein-ligand interactions (PDB: 1GKC)

Protein-ligand interaction images (such as Figure 2) are generated in the following steps:

  1. Determining the 2D layout of the ligand

    The 2D coordinates of the ligand are generated using the OEPrepareDepictionFrom3D function. The 2D layout is driven by the 3D coordinates of the bound ligand.

    ../_images/complex2img-ligand-layout.png

    Step 1. Determining the 2D layout of the ligand

  2. Visualizing the shape of the active site

    The shape of active site is represented with a continuous gray line around the ligand. First, 2D molecule surface is constructed from adjoining arcs. The radius of each arc reflects the distance between the molecule surface of the ligand and the protein in 3D. The absence of an arc indicates regions where the ligand is exposed to the solvent. This 2D molecule surface is then smoothed to get the final aesthetic pleasing representation.

    ../_images/complex2img-pocket.png

    Step 2. Visualizing the shape of the active site

  3. Positioning the residues

    The residues are positioned close to the ligand atoms they interacts with. The distances of residue glyphs for the 2D ligand do not indicate distances in the 3D complex. Residues that only have contact interaction with the ligand are positioned further away, while those that have other type of interactions such as hydrogen bonding, salt bridge etc. are positioned closer to the atoms they are interacting with

    Step 3. Positioning the residues (click on any residues in the image above)

  4. Visualizing the residues and interactions

    The type of the residues and other components around the ligand are represented by different shapes. For example, nearby water is depicted with a glyph resembling a water drop. Basic amino acid residues are represented with a circle and further annotated whether their atoms that are interact with ligand are in the backbone, side-chain or in both.

    Legend of residue and interaction styles

    Apart from the color each interaction type is also associated with a linker type. These linkers indicates which ligand atom(s) are interacting with which residue(s). Some of these linkers also imply the direction of the iterations. For example in case of the hydrogen bond interaction (see the first image in the table below), the chevron arrow shows the direction of the proton.

    Examples of various interaction types
    ../_images/complex2img-hbond.png ../_images/complex2img-clash.png ../_images/complex2img-chelator.png ../_images/complex2img-covalent.png
    hydrogen bond clash chelator covalent
    ../_images/complex2img-salt-bridge-ligand-positive.png ../_images/complex2img-salt-bridge-ligand-negative.png ../_images/complex2img-cation-pi-ligand-pi.png ../_images/complex2img-cation-pi-protein-pi.png
    salt-bridge ligand (+) salt-bridge ligand (-) cation-pi (ligand aromatic ring) cation-pi (protein aromatic ring)
    ../_images/complex2img-stacking-pi.png ../_images/complex2img-stacking-t.png ../_images/complex2img-halogen.png  
    pi-stacking t-stacking halogen bond  

    If a residue participates in more than one type of interactions than the residue glyph has multiple colors (see examples below).

    Example of of residues contributing to more than one types of interactions
    ../_images/complex2img-multi-amino.png ../_images/complex2img-multi-other.png

Hydrogen position optimization

Since the interaction perception dependents on the position of hydrogens, it is highly recommended to optimize those positions prior to perceiving the interactions. The two images below reveal the effect of optimizing the hydrogen bond network in a protein-ligand complex: less atom clash(es) and more hydrogen bond interaction(s).

Table 2. Example of visualizing interactions before and after optimizing hydrogen positions
original complex downloaded from PDB Database (PDB 1D3H) complex after optimizing hydrogen positions

See also

See also in OEChem TK manual

Theory

API

See also in OEDepict TK manual

Theory

API