You want to visualize protein-ligand interactions. See example in Figure 1.
Figure 1. Example of visualizing protein-ligand interactions (PDB: 13GS)
The depict_complex illustrates how to generate protein-ligand iteration images.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20
def depict_complex(image, protein, ligand, opts, legendframe=None): # 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 legendframe is not None: lopts = oegrapheme.OE2DActiveSiteLegendDisplayOptions(12, 1) oegrapheme.OEDrawActiveSiteLegend(legendframe, adisp, lopts)
prompt > wget https://files.rcsb.org/download/13gs.pdb prompt > python3 complex2img.py -complex 13gs.pdb -out 13gs.svg -popuplegend
Currently the OEPerceiveInteractionHints function perceives the following interaction types:
¶ 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.
Protein-ligand interaction images (such as Figure 2) are generated in the following steps:
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.
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.
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)
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.
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.
|salt-bridge ligand (+)||salt-bridge ligand (-)||cation-pi (ligand aromatic ring)||cation-pi (protein aromatic ring)|
If a residue participates in more than one type of interactions than the residue glyph has multiple colors (see examples below).
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).
|original complex downloaded from PDB Database (PDB 1D3H)||complex after optimizing hydrogen positions|