- Open Access
LeView: automatic and interactive generation of 2D diagrams for biomacromolecule/ligand interactions
© Caboche; licensee Chemistry Central Ltd. 2013
- Received: 23 June 2013
- Accepted: 28 August 2013
- Published: 29 August 2013
2D diagrams are widely used in the scientific literature to represent interactions between ligands and biomacromolecules. Such schematic diagrams are very helpful to better understand the chemical interactions and biological processes in which ligands are involved. Here, a new tool for automatic and interactive generation of 2D diagrams for biomacromolecule/ligand interactions is presented. LeView (Ligand-Environment Viewer) produces customised and high-quality figures, with a good compromise between a faithful representation of the 3D data (structures and interactions) and aesthetic criteria. LeView can be freely downloaded at http://www.pegase-biosciences.com/tools/leview/.
- 2D diagram
Non-covalent interactions, such as hydrogen bonds, between ligands and biomacromolecules forming a complex provide important clues in the understanding of biological processes and the design of compounds with desirable binding properties. Schematic 2D diagrams are widely used to visualise the binding interactions of ligands, which can provide important clues about their role and function. A lot of software packages are available for visualising structure in 3D, but only a handful exists for generating 2D protein/ligand interaction diagrams. LIGPLOT  is the most commonly used and its successor LigPlot+  includes a new interface, superposition of related diagrams and links to PyMOL and RasMol. PoseView  automatically generates structure diagrams of molecular complexes. A module of the commercial software MOE - Molecular Operating Environment -  generates schematic diagrams for protein-ligand complexes, and other commercial software vendors offer similar tools. In this paper, LeView (Ligand-Environment Viewer) is introduced. This new tool automatically detects ligands in a PDB file, generates 2D diagrams for biomacromolecule/ligand interactions and allows the user to customise diagrams through an interactive and intuitive graphical interface. LeView is written in Java and distributed as an executable file (source code is also available under GNU Public Licence, Additional file 1) in two versions: a version with an interactive graphical interface and a command-line version for use in pipelines. LeView does not require any special installation steps, can be run on all operating systems.
For each structural element, all the atoms and their neighbours are recursively placed. The recursive algorithm starts with the longest chain or with the larger ring otherwise. In the case of EPO (Figure 2a), algorithm starts with CG1 atom, one of the terminal atom of the longest chain, keeping the initial orientation of the chain. CB, one of the CG1 neighbours, is then processed followed by its neighbours (CG2 and CA). All the chain atoms are recursively processed with their neighbours (C, N, C4, O4, C3, C2, O2, C1, O1, O5, C5, C6 and CG1). The following atom in the call stack is C: 2D coordinates of C atom are previously calculated and are then calculated for its neighbours O and N2 atoms. when N2 atom, which is included in a ring, is processed the next step place all the ring atoms (CD3, CG3, CB1, and CA1). C7 atom is finally placed with their neighbours OXT and O3. Figure 2b shows the simple orthographic projection of the 3D structure of EPO and Figure 2c shows the resulting EPO coordinates obtained with the ligand layout algorithm.
Keeping the initial orientation of structural elements could lead to atom clashes and bond crossings. For example, Figure 2c shows that some conflicts occur. The next step of the ligand layout algorithm tries to resolves conflicts. A conflict is detected when the distance between two atoms is less than 0.4 Å. In this case, all torsion angles which do not imply ring atoms are flipped (0 to 180 degrees) if this change reduces the overlaps in the molecule. An explicit residue, the cystein 29 of chain A (CYS29) is connected to the ligand in the PDB entry 1CSB through a LINK record wich specify connectivity between residues that is not implied by the primary structure. The ligand layout algorithm is applied to the explicit residue CYS29. Once the residue 2D coordinates are obtained, the residue is randomly placed on a circle with a radius of the real 3D distance. The number of bond crossings and atoms overlaps with the ligand are calculated for each position obtained by rotating the residue through 10 degrees each time. The returned position is the position minimising the number of bond crossings and atom overlaps. Figure 2d shows the resulting conflict-free layout of EPO and the explicit residue CYS29.
A list of neighbouring residues that have at least one non-hydrogen atom within a cut-off distance of 4 Å by default from any ligand atom is constructed. Hydrogen-bond donor and acceptor atoms in residues on this list that lie within a cut-off distance to a ligand atom are identified. Donor and acceptor atoms are also identified in the ligand. For each pair of donor(D)/acceptor (A) atoms, a hydrogen bond is created if an angle (D,A,aa) is less than 90 degrees exists (with aa, an acceptor’s neighbour atom). For water-mediated hydrogen bonds, water molecules within 3.3 Å of the donor/acceptor atoms in the ligand are identified as well as all paths between a ligand atom and a residue atom containing between one and four water molecules. Neighbouring residues are defined as biomacromolecule residues which have at least one heavy atom approaching within a cut-off distance to the ligand. Once all interactions are identified, the next step is to place all interactions around the ligand, minimising bond crossings and overlapping graphical elements. The same method is used for placing all interaction types. The interacting element is randomly placed on a circle with a radius of the real 3D distance. The number of bond crossings and nearby atom, i.e. atoms approaching within a distance less than 0.4 Å are calculated for each position obtained by rotating the interacting element through 5 degrees steps. The best position is that which minimises the number of bond crossings. If several positions show the same minimum number of bond crossings, the selection is made based on the position that has minimal close interactions. LeView first places the hydrogen bonds, then the nearby residues and finally any water-mediated hydrogen bonds. LeView starts to place hydrogen bonds because they are represented by lines implying a higher probability of atom overlaps and bond crossings than nearby residues which are represented by ovals and placed in second. Water-mediated hydrogen bonds are not displayed by default and when the user add one, the coordinates are calculed to include it on the existing diagram. Figure 2e shows the layout of hydrogen bonds and nearby residues in the example of EPO ligand.
The finally step of the layout algorithm is to add text on the diagram. The name of each residue component of the ligand and the name of explicit residues are placed on the diagram in order to minimise overlaps with other diagram elements. The residue name is randomly placed on a circle with the centre being the centroid of the residue which is the arithmetic mean position of all the residue atoms and a radius of M a x(d i s t a n c e(c e n t r o i d,i))+0.5 with i the residue atoms, i.e. the distance between the centroid and the further residue atom plus 0.5 Å. The number of nearby elements approaching within a distance less than 2 Å is calculated for each position obtained by rotating the text through 10 degrees each time. The position minimising the number of nearby elements is kept. Figure 2f shows the final diagram generated with default settings by LeView for EPO in the PDB entry 1CSB.
LeView automatically detects biomacromolecules, ligands and metal ions in a PDB file and shows the user a complete list of these. Composite ligands can also be user specified by entering the appropriate residue range. It is possible to create an interactive 2D diagram for each ligand and metal ion and this shows ligands, metals, hydrogen bonds and nearby residues. Nearby residues are defined as biomacromolecule residues (i.e. in proteins, DNA or RNA) which have at least one heavy atom approaching within a cut-off distance to the ligand. LeView allows the user to display several ligands and metal ions from the same or from several different PDB entries, which makes it easy to compare diagrams. Figure 2f shows the 2D diagram produced by LeView using default settings for EPO in PDB entry 1CSB. LeView offers the user a highly functional and intuitive interactive graphical interface to customise the 2D diagram. The user can change the cut-off distance for hydrogen bonds and nearby residues by sliding the corresponding distance cut-off bar. Possible water-mediated interactions involving the ligand or ion, with up to four bridging water molecules, are not displayed by default but can be added to the diagram through the menu. A number of display options are available as well. Atom labels (e.g. C15) can be displayed or hidden and the user can choose between standard colours (e.g. red for oxygen) or plain-colour mode, i.e. ligand atoms appear in the same colour than the ligand chemical bonds while several colour schemes are available for nearby residues. Colour can be used to represent hydrophobicity/hydrophilicity, standard amino acid properties similar to the Shapely scheme available in RasMol , charge or secondary structure type. The colour of every element in the diagram can be changed via the colour menu. Hydrogen bonds can be represented by solid lines (default) or arrows from donor to acceptor, and hydrogen-bond distances additionally displayed. LeView allows the user to reposition all the graphical elements in the diagram at will. If desired, individual hydrogen bonds and nearby residues can be deleted from the diagram, by a right click on it. Finally, the diagram can be exported in a variety of raster and vector graphics formats: PNG, GIF, JPG, PDF, SVG and EPS. For the latter three formats, LeView uses the VectorGraphics package of the FreeHEP Java Library (http://java.freehep.org). The list of the interactions represented in the diagram, with the atoms involved and the associated distance, can also be exported as a flat text.
Main features of non-commercial programs to generate 2D diagrams of protein-ligand interactions
number of entries
PNG, JPG, GIF, PDF
SVG, EPS, flat text
LeView produces customisable and high-quality figures often necessary for scientific publications, in several popular formats. It is able to work with a wide range of complex ligands and environments not possible with existing tools. This tool can be freely downloaded and easily used without installation and offers the user an interactive graphical interface with numerous options to customise the figure such as: varying the cut-off distances, changing the diagram element colours, moving and deleting elements. The diagrams are a good compromise between a faithful representation of the 3D data (structures and interactions) and aesthetic criteria. Ligand, hydrogen bonds and nearby residues are included and the list of possible water-mediated interactions is also available.
Project name: LeViewProject home page: http://www.pegase-biosciences.com/tools/leview/Operating system(s): Platform independentProgramming language: JavaOther requirements: Java 1.5 or higherLicense: GNU General Public Licence Any restrictions to use by non-academics: no restrictions
The author gratefully acknowledges the EBI for funding and Protein Data Bank in Europe (PDBe; http://pdbe.org) staff for help, testing and discussions. The author thanks Gerard Kleywegt and Tom Oldfield for advice and for manuscript proof reading and Roman Laskowski for his help in ligand detection from PDB files. This work was supported by the European Bioinformatics Institute, part of the EMBL. Funding for publication charge: IFR142
- Wallace AC, Laskowski RA, Thornton JM: LIGPLOT: a program to generate schematic diagrams of protein-ligand interactions. Protein Eng. 1995, 8 (2): 127-134. 10.1093/protein/8.2.127.View ArticleGoogle Scholar
- Laskowski RA, Swindells MB: LigPlot+: multiple ligand-protein interaction diagrams for drug discovery. J Chem Inf Model. 2011, 51 (10): 2778-2786. 10.1021/ci200227u.View ArticleGoogle Scholar
- Stierand K, Maass PC, Rarey M: Molecular complexes at a glance: automated generation of two-dimensional complex diagrams. Bioinformatics. 2006, 22 (14): 1710-1716. 10.1093/bioinformatics/btl150.View ArticleGoogle Scholar
- Clark AM, Labute P: 2D depiction of protein-ligand complexes. J Chem Inf Model. 2007, 47 (5): 1933-1944. 10.1021/ci7001473.View ArticleGoogle Scholar
- Laskowski RA, Hutchinson EG, Michie AD, Wallace AC, Jones ML, Thornton JM: PDBsum: a Web-based database of summaries and analyses of all PDB structures. Trends Biochem Sci. 1997, 22 (12): 488-490. 10.1016/S0968-0004(97)01140-7.View ArticleGoogle Scholar
- Turk D, Podobnik M, Popovic T, Katunuma N, Bode W, Huber R, Turk V: Crystal structure of cathepsin B inhibited with CA030 at 2.0-A resolution: a basis for the design of specific epoxysuccinyl inhibitors. Biochem. 1995, 34 (14): 4791-4797. 10.1021/bi00014a037.View ArticleGoogle Scholar
- Sayle RA, Milner-White EJ: RASMOL: biomolecular graphics for all. Trends Biochem Sci. 1995, 20 (9): 374-10.1016/S0968-0004(00)89080-5.View ArticleGoogle Scholar
- Sung BJ, Hwang KY, Jeon YH, Lee JI, Heo YS, Kim JH, Moon J, Yoon JM, Hyun YL, Kim E, Eum SJ, Park SY, Lee JO, Lee TG, Ro S, Cho JM: Structure of the catalytic domain of human phosphodiesterase 5 with bound drug molecules. Nature. 2003, 425 (6953): 98-102. 10.1038/nature01914.View ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.