- Open Access
iPBAvizu: a PyMOL plugin for an efficient 3D protein structure superimposition approach
Source Code for Biology and Medicine volume 14, Article number: 5 (2019)
Protein 3D structure is the support of its function. Comparison of 3D protein structures provides insight on their evolution and their functional specificities and can be done efficiently via protein structure superimposition analysis. Multiple approaches have been developed to perform such task and are often based on structural superimposition deduced from sequence alignment, which does not take into account structural features. Our methodology is based on the use of a Structural Alphabet (SA), i.e. a library of 3D local protein prototypes able to approximate protein backbone. The interest of a SA is to translate into 1D sequences into the 3D structures.
We used Protein blocks (PB), a widely used SA consisting of 16 prototypes, each representing a conformation of the pentapeptide skeleton defined in terms of dihedral angles. Proteins are described using PB from which we have previously developed a sequence alignment procedure based on dynamic programming with a dedicated PB Substitution Matrix. We improved the procedure with a specific two-step search: (i) very similar regions are selected using very high weights and aligned, and (ii) the alignment is completed (if possible) with less stringent parameters. Our approach, iPBA, has shown to perform better than other available tools in benchmark tests. To facilitate the usage of iPBA, we designed and implemented iPBAvizu, a plugin for PyMOL that allows users to run iPBA in an easy way and analyse protein superimpositions.
iPBAvizu is an implementation of iPBA within the well-known and widely used PyMOL software. iPBAvizu enables to generate iPBA alignments, create and interactively explore structural superimposition, and assess the quality of the protein alignments.
The detection of structural analogy between protein folds requires development of methods and tools to compare and classify them. This is extremely helpful for studying evolutionary relationships between proteins especially in the low sequence identity ranges . However, an optimal superposition is far from being a trivial task. Popular methods such as DALI  and CE , use a reduced representation of backbone conformation in terms of distance matrices.
Protein backbone conformation can be characterized by a set of local structure prototypes, namely Structural Alphabets (SAs), which enables the transformation of 3D information into a 1D sequence of alphabets . Hence a 3D structure comparison can be obtained by aligning sequences of SAs (protein structures encoded in terms of SA). A SA consisting of 16 pentapeptide conformations, called Protein Blocks (PBs), was developed in our group . Based on this library, a protein superimposition approach was developed. A substitution matrix for PBs  was generated based on all PB substitutions observed in pairwise structure alignments in PALI dataset . The superimposition was carried out with simple dynamic programming approaches . We recently improved the efficiency of our structural alignment algorithm by (i) refining the substitution matrix and (ii) designing an improved dynamic programming algorithm based on preference for well-aligned regions as anchors. This improvement (improved Protein Block Alignment, iPBA) resulted in a better performance over other established methods like MUSTANG  for 89% of the alignments and DALI for 79% . Benchmarks on difficult cases of alignment also show similar results [11, 12]. Protein Blocks were also recently used to analyse Molecular Dynamic simulations [13, 14] underlining their abilities to apprehend protein flexibility .
We present here a plugin, iPBAvizu, which integrates the efficient protein structure alignment approach iPBA with the very popular molecular graphics viewer PyMOL (The PyMOL Molecular Graphics System, Version 1.7, Schrödinger, LLC) from which several plugins like PyKnoT  or PyETV  have been integrated in. iPBAvizu enables interactive visualization and analysis of protein structure superposition and the resulting sequence alignment. Different scores to assess the quality of the alignment are also given.
After installing all the dependencies, iPBAvizu can be easily integrated within PyMOL using the ‘Plugin’ menu on the PyMOL console, choosing ‘Install’ under the ‘Manage Plugins’ and then locating and selecting the iPBAvizu.py file. The installation procedures as well as few examples of structural alignments are illustrated in a series of videos (see http://www.dsimb.inserm.fr/dsimb_tools/iPBAVizu/). The plugin is easy to use and does not require any command line or programming skills. It is fully controlled by the PyMOL GUI.
To launch iPBAvizu from the PyMOL Wizard menu, at least two protein structures must be loaded and made available in the PyMOL session. iPBAvizu menu appears in PyMOL GUI, like the Measurement or Fit native functions. Users can select two chains among the available loaded structures, and then select ‘Align!’ to run iPBA program. Once the alignment process is over, results are displayed as two new protein objects in PyMOL. The two new objects correspond to the two aligned structures. A new window containing different alignment scores (e.g., GDT-TS, RMSD, see Methods) and an interactive sequence alignment manager is also displayed. Both residue and Protein Block sequences of aligned structures are given. Users can highlight any residue or PB of one or both sequences. Highlighting selects the residues directly in the 2 new aligned protein objects created in PyMOL 3D window. This interactive functionality provides an efficient way to explore sequence and structural alignment.
Figure 1 shows an example of structural superposition of two proteins of the monooxygenase protein family using iPBAvizu plugin: Cyclohexanone Monooxygenase (CHMO, PDB code 3GWD) and Phenylacetone Monooxygenase (PAMO, PDB code 1W4X) . The obtained results were also compared with other popular superimposition tools (e.g., cealign  and TM-align ). The alignment generated by iPBA based on PBs was compared to alignment generated with cealign and TM-align and the iPBA alignment show a better Cα RMSD score (1.5 Å versus values between 1.9–2.7 Å for the 2 other approaches). The values are provided for the aligned residues that are on average larger than with other superimposition tools.
Discussion & Conclusion
A structural alphabet is a library of protein fragments able to approximate every part of protein structures (for a review ). These libraries yielded prototypes that are representative of local folds found in proteins. The structural alphabet allows the translation of three-dimensional protein structures into a series of letters. As a result, it is possible to use classical sequence alignment methodologies to perform structural alignments. The main difficulty lies in obtaining a pertinent substitution matrix that gives the similarity score between alphabets, which guides the alignments. Few teams have used this approach to perform structural comparisons and/or PDB mining:
Guyon and co-workers had used a structural alphabet based on Hidden Markov Model and proposed an approach named SA-search (http://bioserv.rpbs.jussieu.fr/cgi-bin/SA-Search, ). Their substitution matrix is generated from a transition matrix, however the details of the method are uncleared. The webserver gives only C-alpha coordinates for superimposition and does not provide a fully interactive interface to explore structural alignment. Finally, SA-Search webserver has not been updated since 2006 and miss modern web-technology interactivity based.
3D-BLAST was developed late 2006 and is based on the BLAST methods . The structural alphabet proposed is based on optimization of nearest-neighbor clustering (NNC). Interestingly the substitution matrix was generated based on SCOP classification. Since 3D-BLAST was initially developed to search for structural similarity and not to specifically compare two protein structures of interest, it was not benchmark. The webserver (http://3d-blast.life.nctu.edu.tw/) needs Chime applet, and users do not have a direct access to simple alignment results.
SA-FAST was developed for the same purpose  but was based on FASTA algorithm. Structural alphabet was generated using a Self-Organizing Map, taking into account the most frequent clusters. The final benchmark was done using 50 proteins. The webserver (http://bioinfo.cis.nctu.edu.tw/safast/) is very fast. However, it is not possible to do simple pairwise alignments and the output needs Chime applet which is not very easy to install. The major drawback is that users do not have access to the alignment by itself for further analysis.
CLePAPS  is based on the use of a dedicated structural alphabet built only to perform database search. In the first step, aligned fragment pairs (AFP) are found, which correspond to fragments that involve exact matches of similar letters. CLePAPS then joins consistent AFPs guided by their similarity scores to extend the alignment by several “zoom-in” iteration steps; it does not use dynamic programming. CLePAPS was tested on a limited number of protein structure pairs. A stand-alone program is reported to be available, but not found.
Hence, iPBAvizu is quite interesting approach. Indeed, it is an easy-to-use plugin for PyMOL that allows users to superimpose protein structures using iPBA methodology, an efficient way to superimpose protein 3D structures  and explore the structural alignment results. Its total integration as a plugin into PyMOL molecular viewer offers an easy but powerful way to process and study structural alignment with quantitative measurements.
Materials and methods
iPBA program is fully written in Python (2.7+). It depends on ProFit program stand-alone version (Martin, A.C.R., http://www.bioinf.org.uk/software/profit) for generating the final structural alignment. iPBA provides an efficient way to align two protein structures using anchor-based alignment methodology [11, 12].
iPBAvizu package has an installer to configure iPBA and manage its dependencies on the local machine before integrating it into PyMOL. Due to ProFit requirements, iPBAvizu is only available on Unix-based operating systems. iPBAvizu is embedded into PyMOL as a wizard plugin, and all iPBA functionalities are totally integrated into the graphic interface of PyMOL. iPBAvizu can be launched with the current PyMOL internal GUI. Users can easily align structures with a few clicks and access both scores and the alignment results that are displayed in PyMOL itself, as a Tkinter GUI. The alignment window is interactive; it is linked to 3D PyMOL interface for the best interpretation and exploration of results.
iPBA and iPBAvizu can estimate the quality of the superimposition via a score. The GDT score (GDT_TS) is widely used for the assessment of structural models generated in CASP structure prediction trials , it is supposed to be less sensible to large deviation as seen with Root Mean Square Deviation (RMSD). The GDT_TS is the combination of set of superimposed residues for fixed thresholds at 1, 2, 4 and 8 Å. GDT_PB scores (calculated in a similar way as that of GDT_TS, but using PB substitution scores [11, 12] instead of distances) are also provided for the hits obtained (see for [11, 12] more details).
Protein Blocks (PB) and amino acid sequences are provided. PB is the most widely used structural alphabet and is composed of 16 local prototypes  of five residue length, it is dedicated to analyse local conformations of protein structures from the Protein DataBank (PDB) . Each PB is characterized by the φ and ψ dihedral angles of five consecutive residues. PBs give a reasonable approximation of all local protein 3D structures [14, 27, 28]. PBs are labelled from a to p. PBs m and d can be roughly described as prototypes for α-helix and central β-strand, respectively. PBs a to c primarily represent β-strand N-caps and PBs e and f representing β-strand C-caps; PBs g to j are specific to coils; PBs k and l to α-helix N-caps while PBs n to p to α-helix C-caps. For each PB is associated 5 residues, its assignment is done on the central residue. As PBs are overlapping, a structure of length N is translated in N-4 PBs, the two first and two last residues are associated to letter Z (see Fig. 1). Missing residues are also associated to the letter Z.
Availability of data and materials
iPBAvizu is a PyMOL plugin freely available to the academic scientific community, i.e. the data is only informatics codes. It is composed of the PyMOL script code and the iPBA code. This last used python and some C codes. The downloadable archive can be freely accessed at our academic website: http://www.dsimb.inserm.fr/dsimb_tools/iPBAVizu/. As it is a PyMOL plugin, user needs to install independently PyMOL software: https://pymol.org. There is no restriction for use or modifications of iPBAvizu by any academic scientists. For commercial usage, please contact the authors.
Agarwal G, Rajavel M, Gopal B, Srinivasan N. Structure-based phylogeny as a diagnostic for functional characterization of proteins with a cupin fold. PLoS One. 2009;4(5):e5736.
Holm L, Sander C. Protein structure comparison by alignment of distance matrices. J Mol Biol. 1993;233(1):123–38.
Shindyalov IN, Bourne PE. Protein structure alignment by incremental combinatorial extension (CE) of the optimal path. Protein Eng. 1998;11(9):739–47.
de Brevern AG, Etchebest C, Hazout S. Bayesian probabilistic approach for predicting backbone structures in terms of protein blocks. Proteins. 2000;41(3):271–87.
Joseph AP, Agarwal G, Mahajan S, Gelly JC, Swapna LS, Offmann B, Cadet F, Bornot A, Tyagi M, Valadié H. A short survey on protein blocks. Biophys Rev. 2010;2:137–45.
Tyagi M, Gowri VS, Srinivasan N, de Brevern AG, Offmann B. A substitution matrix for structural alphabet based on structural alignment of homologous proteins and its applications. Proteins. 2006;65(1):32–9.
Balaji S, Sujatha S, Kumar SS, Srinivasan N. PALI-a database of phylogeny and ALIgnment of homologous protein structures. Nucleic Acids Res. 2001;29(1):61–5.
Tyagi M, de Brevern AG, Srinivasan N, Offmann B. Protein structure mining using a structural alphabet. Proteins. 2008;71(2):920–37.
Konagurthu AS, Whisstock JC, Stuckey PJ, Lesk AM. MUSTANG: a multiple structural alignment algorithm. Proteins. 2006;64(3):559–74.
Holm L, Park J. DaliLite workbench for protein structure comparison. Bioinformatics. 2000;16(6):566–7.
Joseph AP, Srinivasan N, de Brevern AG. Improvement of protein structure comparison using a structural alphabet. Biochimie. 2011;93(9):1434–45.
Gelly JC, Joseph AP, Srinivasan N, de Brevern AG. iPBA: a tool for protein structure comparison using sequence alignment strategies. Nucleic Acids Res. 2011;39(Web Server issue):W18–23.
Barnoud J, Santuz H, Craveur P, Joseph AP, Jallu V, de Brevern AG, Poulain P. PBxplore: a tool to analyze local protein structure and deformability with protein blocks. PeerJ. 2017;5:e4013.
Goguet M, Narwani TJ, Petermann R, Jallu V, de Brevern AG. In silico analysis of Glanzmann variants of Calf-1 domain of alphaIIbbeta3 integrin revealed dynamic allosteric effect. Sci Rep. 2017;7(1):8001.
Craveur P, Joseph AP, Esque J, Narwani TJ, Noel F, Shinada N, Goguet M, Leonard S, Poulain P, Bertrand O, et al. Protein flexibility in the light of structural alphabets. Front Mol Biosci. 2015;2:20.
Lua RC. PyKnot: a PyMOL tool for the discovery and analysis of knots in proteins. Bioinformatics. 2012;28(15):2069–71.
Lua RC, Lichtarge O. PyETV: a PyMOL evolutionary trace viewer to analyze functional site predictions in protein complexes. Bioinformatics. 2010;26(23):2981–2.
Rebehmed J, Alphand V, de Berardinis V, de Brevern AG. Evolution study of the Baeyer-Villiger monooxygenases enzyme family: functional importance of the highly conserved residues. Biochimie. 2013;95(7):1394–402.
Zhang Y, Skolnick J. TM-align: a protein structure alignment algorithm based on the TM-score. Nucleic Acids Res. 2005;33(7):2302–9.
Offmann B, Tyagi M, de Brevern AG. Local Protein Structures. Curr Bioinforma. 2007;3:165–202.
Guyon F, Camproux AC, Hochez J, Tuffery P. SA-Search: a web tool for protein structure mining based on a Structural Alphabet. Nucleic Acids Res. 2004;32(Web Server issue):W545–8.
Yang JM, Tung CH. Protein structure database search and evolutionary classification. Nucleic Acids Res. 2006;34(13):3646–59.
Ku SY, Hu YJ. Protein structure search and local structure characterization. BMC Bioinformatics. 2008;9:349.
Wang S, Zheng WM. CLePAPS: fast pair alignment of protein structures based on conformational letters. J Bioinforma Comput Biol. 2008;6(2):347–66.
Zemla A. LGA: a method for finding 3D similarities in protein structures. Nucleic Acids Res. 2003;31(13):3370–4.
Berman HM, Westbrook J, Feng Z, Gilliland G, Bhat TN, Weissig H, Shindyalov IN, Bourne PE. The protein data Bank. Nucleic Acids Res. 2000;28(1):235–42.
Joseph AP, Agarwal G, Mahajan S, Gelly J-C, Swapna LS, Offmann B, Cadet F, Bornot A, Tyagi M, Valadié H, et al. A short survey on protein blocks. Biophys Rev. 2010;2(3):137–45.
Narwani TJ, Craveur P, Shinada NK, Floch A, Santuz H, Melarkode Vattekatte A, Srinivasan N, Rebehmed J, Gelly JC, Etchebest C, et al. Discrete analyses of protein dynamics. J Biomol Struct Dyn. 2019:1–23.
We would like to thank Nicolas Shinada and Akhila Melarkode Vattekatte for fruitful discussions.
This work was supported by grants from the French Ministry of Research, University of Paris Diderot – Sorbonne Paris Cité, University de la Réunion, University des Antilles, French National Institute for Blood Transfusion (INTS), French Institute for Health and Medical Research (INSERM). AGdB, APJ, NS and TJN acknowledge the Indo-French Centre for the Promotion of Advanced Research / CEFIPRA for collaborative grants (number 3903-E and 5203–2). AdB, and JR acknowledge ANR NaturaDyRe (France, ANR-2010-CD2I-014-04). This study was supported by grants from the Laboratory of Excellence GR-Ex (reference ANR-11-LABX-0051). The labex GR-Ex is funded by the programme “Investissements d’avenir” of the French National Research Agency, reference ANR-11-IDEX-0005-02. Calculations were performed on an SGI cluster granted by Conseil Régional Ile de France and INTS (SESAME Grant). The authors were granted access to high performance computing (HPC) resources at the French National Computing Centre CINES under grants no. c2013037147, c2016077621 and A0010707621 funded by the GENCI (Grand Equipement National de Calcul Intensif).
Research in NS group is supported by Mathematical Biology program and FIST program sponsored by the Department of Science and Technology and also by the Department of Biotechnology, Government of India in the form of IISc-DBT partnership programme. Support from UGC, India – Centre for Advanced Studies and Ministry of Human Resource Development, India is gratefully acknowledged. NS is a J. C. Bose National Fellow.
The funding bodies have no roles in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Faure, G., Joseph, A.P., Craveur, P. et al. iPBAvizu: a PyMOL plugin for an efficient 3D protein structure superimposition approach. Source Code Biol Med 14, 5 (2019) doi:10.1186/s13029-019-0075-3
- Protein superimposition
- Structural alphabet
- Structural alignment
- Structural bioinformatics