MotViz: a tool for sequence motif prediction in parallel to structural visualization and analyses.

Linking similar proteins structurally is a challenging task that may help in finding the novel members of a protein family. In this respect, identification of conserved sequence can facilitate understanding and classifying the exact role of proteins. However, the exact role of these conserved elements cannot be elucidated without structural and physiochemical information. In this work, we present a novel desktop application MotViz designed for searching and analyzing the conserved sequence segments within protein structure. With MotViz, the user can extract a complete list of sequence motifs from loaded 3D structures, annotate the motifs structurally and analyze their physiochemical properties. The conservation value calculated for an individual motif can be visualized graphically. To check the efficiency, predicted motifs from the data sets of 9 protein families were analyzed and MotViz algorithm was more efficient in comparison to other online motif prediction tools. Furthermore, a database was also integrated for storing, retrieving and performing the detailed functional annotation studies. In summary, MotViz effectively predicts motifs with high sensitivity and simultaneously visualizes them into 3D strucures. Moreover, MotViz is user-friendly with optimized graphical parameters and better processing speed due to the inclusion of a database at the back end. MotViz is available at http://www.fi-pk.com/motviz.html.

Introduction

Computationally understanding and annotating the central dogma of molecular biology encompasses the success of bioinformatics. From sequence to its annotation as well as from protein sequence to structure prediction and visualization, all processes are being executed with high efficiency through in silico tools and databases. It is noteworthy that number of protein sequences grows rapidly and is far beyond the number of experimental structures. According to a recent description, SWISS-PROT ( database contains 525, 997 protein sequence entries, while Protein Data Bank (PDB) ( comprises only 71,794 protein structures. The tremendous amount of sequence data is in dire need of analysis in order to understand the biological meanings and to establish appropriate relationships between sequences and structures. Considering that protein structures are relatively conserved and adopt only a limited number of folds, it is possible to predict their 3D structures by knowledge-based information (. Devising and implementing novel algorithms will largely contribute to determination of protein structure or function.

Generally, to achieve insight into structures through Cartesian coordinates, it is essential to have a visualization tool. During recent years, a variety of data visualization tools have been developed to meet increasing data complexity. These include Rasmol, VMD, Cn3D, Swiss PDB Viewer, Chimera, Jmol and PyMol 4, 5, 6, 7, 8, 9, 10. Jmol ( is a popular Java-based freeware and standalone high-end tool for protein structural monitoring. The browser applet in Jmol supports the loading of multiple molecules with independent movement, surfaces and molecular orbitals, biological units and crystal symmetry, cavity visualization, translucency, high-quality rendering, arbitrary objects such as arrows/planes and true slabbing properties.

In the current study, a plug-in MotViz has been built and added into Jmol version 11.2.1 in order to extend its applicability. Specialized features in MotViz include sequence motif retrieval from protein 3D structure in parallel to structural analysis based on the physiochemical properties. Additionally, a local database has been created to facilitate user in re-evaluation of results with better performance. Finally, for detailed annotation of selected motifs, links of web-based tools like Prosite, MEME and Hits 11, 12, 13 have been incorporated into MotViz.

Methods

We updated Jmol protein visualization tool version 11.2.1 by incorporating new plug-in and modules. These new features include sequence annotation, motif searching, motif visualization and graphical output of their physiochemical chart, motifs conservation score calculation and their comparisons through bar graph representation. MotViz algorithm created in this study predicts motifs from the loaded 3D protein structure. This algorithm uses a dynamic approach and latest-updated data from UniProt and InterPro 14, 15 databases. MotViz algorithm works in five phases, which are described as follows.

Phase 1 Retrieving amino acid sequence from 3D coordinates

Protein sequence is dissected from 3D coordinates of determined amino acids at Phase 1. MotViz retrieves protein sequence from the structural coordinates, because the target of MotViz is to hunt structural motifs in parallel to sequence motifs, as some amino acids might be missing in 3D coordinates but likely available in the protein Fasta sequence.

Phase 2 Launching BLASTP and fetching Fasta sequences

At the second phase MotViz retrieves similar protein sequences related to the visualized protein sequence to hunt for conserved motifs. UniProt BlastP is launched by calling the EBI web service to obtain IDs of the 50 most conserved sequences (with a maximum difference of 10−50 E-value). If this E-value difference is found in about topmost 30 BLASTP sequences, then only those sequences are picked. Fasta sequences of selected proteins are then fetched from the UniProt server and saved for further analysis.

Phase 3 Performing multiple sequence alignment

At the third phase multiple sequence alignment (MSA) of closely-related protein sequences (from phase 2) is performed to locate the conservation points using ClustalW.

Phase 4 Locating conservations and determining motifs

This step identifies conservation points by reading ClustalW “*, :, .”, where “*”, “:” and “.” represents identical, conserved and semi-conserved residue, respectively. MotViz algorithms will identify ---*--:--.--- as the worst case and ---***--- as the best case. Furthermore MotViz also locates the positions of the determined motifs in the protein sequence and stores the information in the database.

MotViz Algorithm Steps

Input: PDB ID of protein

Search: motifs of PDB ID in database

If(motifs)

store: motifs in array m of length k

For i 1 to k

Calculate_physiochemical_properties(m(i))

Input: loc of m(i) location

Output: m(i) and loc

Output: m(i) C value graph

Else

Input Amino Acid sequence from 3D coordinates

Input UniProt Psiblast parameters

Launch UniProt_Psiblast()

getResult in Array results

diff Calculate E-value difference (e−50)location from results

For j 1 to diff do sequence(j) results(j)

Store sequence in MotViz database

Launch ClustalW MSA(sequence)

OutputRes Read ClustalW .output file

Input array new_motifs of dynamic length

Input cont_gap, cons_leng, motif_No, start to 0

for k 1 to OutputRes length

if (cont_gap< 3)

If (OutputRes(k) == * || OutputRes(k) ==: || OutputRes(k) == .)

If (cons_leng = 0)

start k

cons_leng++

Else

cont_gap++

Else if (cont_gap > 2 && cons_leng >2)

Calculate C_Value (start to k)

new_motifs(motif_No) = OutputRes(start to k)

start =0, cons-leng =0, cont_gap=0

Store new_motifs array in MotViz database

Phase 5 Calculating conservation value

At the final phase, conservation of each predicted motifs is calculated using the MotViz algorithm according to the following formula.

Where Cv is conservation value; ‘i’ is from 0 to 3, x is number of identical residues; w is identity score (in this case 1 is optimal score); x is number of completely-conserved residues; w is score for completely-conserved residues (maximum value is 0.85); x is number of semi-conserved residues; w is score for semi-conserved residues (maximum value is 0.70) and is the total length of motif.

PSI-BLAST was used during the algorithmic development of MotViz (. Dependency of PSIBLAST ( is on activation.jar, commons-cli-1.1.jar, commons-logging-1.0.2.jar, mail.jar, servlet.jar, xercesimpl.jar, axis.jar, commons-discovery-0.2.jar, jaxrpx. jar, saaj.jar and wsdl4j-1.4.1.jar, which are available at http://www.ebi.ac.uk/Tools/webservices/services/dbfetch (. ClustalW ( was used for performing MSA of protein sequences retrieved from UniProt.

The conserved segments across all aligned sequences are marked by numerical values and subsequently positioned into the structural coordinates. Figure 1 shows the schematics flow of the MotViz algorithm. Predicted conserved areas are considered to be of functional relevance and are stored in a separate database for easy retrieval.

Figure 1

Schematic overview of motif searching and visualization algorithm (MSVA). A. AA sequence is retrieved from the PDB, visualized and stored into local database. User may examine the respective sequence properties. B. Calculation and displaying the AA physiochemical properties of selected sequence region or motif(s). The corresponding structure is rendered with distinct color. C. Searching in database and displaying the list of motifs present in the loaded PDB file. D. If motifs are not found in the MotViz database, UniProt PSI-BLAST ( is launched and results are stored in MotViz database for MSA. E. MSA is performed using ClustalW ( and conserved regions are identified and stored in MotViz database to carry out step В and step C. The circle at the lower right corner represents inspection to approve the quality of product.

Results

Sequence retrieval and analysis

Protein structure files were obtained from PDB (. Here, we took 3ML6, the structure of complex between dishevlled2 and clathrin adaptor AP-2, as example (Figure 2A). Upon clicking the sequence button, there displays a new panel at the left bottom corner bearing the protein sequence (Figure 2B). The sequence segments are highlighted within protein structure with same color representation (Figure 2B1, 2B2). In the same panel, there are sub-panels for bar and pie charts, button to display the physiochemical properties (hydrophobic, acidic, basic and polar residues) together with slider and input spinner for selecting the protein sequence or to visualize the selected amino acids properties. In addition, two radio-buttons are enabled for selecting the visualization option of either whole structure or selected sequence only. For selected sequence segments, pie and bar chart features are also enabled together with slider selection (Figure 2B3). Multiple visualization aspects that have already been embedded in Jmol (, including balls-sticks, cartoon, space-filled (Vander Waals), ribbon model and wireframe-trace, are enabled for detailed analysis.

Figure 2

MotViz working interface encompassing feature’s track individually. A. Loaded 3ML6 (complex between Dishevlled2 and Clathrin adaptor AP-2) structure. B. Sequence retrieved from PDB structure. B1. Selected amino acids in the sequence are visualized in structure. B2. Relative amino acids highlighted in structure only. B3. Individual amino acid physiochemical properties charts (pie and bar), buttons and slider panel to select the region(s) of choice are indicated. Detailed view of amino acid physiochemical properties charts. Green color represents the polar amino acids (G, S, T, Y, C, Q and N); orange color represents the hydrophobic amino acids (A, V, L, I, P, W, F and M); red color indicates acidic amino acids (D and E); blue color is for basic residues (K, R and H); while other amino acids are represented by Tan color. C1. List of predicted motifs in multi-color panel. Pie chart indicates the AA physiochemical properties of predicted motifs in parallel to their structure representation. C2. Predicted motifs in structure. C3. Selected motifs in motif panel as well as in structure. C4. Bar chart representing the predicted motif conservation rate analysis encompassing the whole protein. Moreover, conservation value (Cv) along with individual motif location for each predicted motif is also indicated.

Motif prediction and visualization

Sequence retrieved from the 3D model of loaded protein structure is stored in a local database (MotViz-Database). Individual motifs are retrieved by clicking the search motif button. Motifs are predicted from conserved sites generated using the MotViz algorithm against the query sequence and stored in the local database (MotViz-Database). Later, a motif panel is displayed after retrieving motifs in the right upper corner (Figure 2C1). The motif panel contains a comprehensive list of all predicted motifs. To facilitate the user, motif positions in retrieved sequence are made available. Also, the individual motifs can be selected by clicking check box (Figure 2C2). Radio-buttons help in selecting either the visualization option or select motifs only within structure. In order to select all motifs, a select all check box is added to display all motifs simultaneously. By selecting this checkbox, all motifs can be visualized in 3D structure accordingly containing the same color scheme to that of corresponding motif sequences. By clicking the selected only button, selected motifs are displayed in the structure (Figure 2C3). The conservation value calculated for each predicted motif along with its location in the respective sequence is displayed. Moreover, a comprehensive bar chart helps in categorizing related predicted motif(s) on the basis of conservation rate (Figure 2C4).

Performance results

Evaluation of the motif search result was carried out by querying MotViz predicted motifs using multiple online motif finding tools. To validate MotViz predicted motifs, the performance of these predicted motifs was tested by comparing the MotViz predicted motifs with other online motif prediction tools including Prosite, MEME, PATTERN, BLOCKS, Pfam, ProDom, PRINTS 11, 12, 19, 20, 21, 22, 23. Table 1 summarized the verified sensitivity, specificity, precision and accuracy of results predicted by MotViz and other motif prediction tools 11, 12, 19, 20, 21, 22, 23 for 9 protein families randomly selected. These protein families include Wnt protein family, Gli protein family, histone protein family, ribosomal protein family, HIV protein family, PHD-finger protein family, Hox protein family, HSP100 protein family and Frizzled protein family. These selected protein families were analyzed (four proteins in average for each family) for the presence of motifs using MotViz. Highly conserved motifs which were predicted only by MotViz were listed in Table 2 along with their conservation value score (Cv). MotViz specificity was calculated by averaging Cv greater than 0.60 and an E-value greater than 1009 for MotViz predicted motifs 11, 12, 19, 20, 21, 22, 23. MotViz sensitivity was measured by taking the averages of Cv and E-values of the predicted motifs.

Table 1

Comparative performance of motif searching tools

Database/Method	Time (s)	Sensitivity	Precision	Specificity	Accuracy	No. of predicted motifs
Prosite	230.12	84.43	94.31	88.92	92.65	37
MEME	348.48	92.11	93.52	95.362	92.02	139
PATTERN	165.63	80.78	96.24	98.43	95.13	13
BLOCKS	1330.73	90.59	90.10	93.62	88.73	99
Pfam	219.12	88.34	89.37	95.10	87.15	78
ProDom	198.34	81.23	94.83	98.32	92.45	15
PRINTS	165.63	82.93	94.67	98.10	91.89	27
MotViz	0.0699	96.18	95.82	97.43	93.45	233

Note: Comparative performance of MotViz with online databases like Prosite, MEME, PATTERN, BLOCKS, Pfam, ProDom, PRINTS 11, 12, 19, 20, 21, 22, 23. Precision was calculated by TP/(TP+FP). TP represents true positive and FP indicates false positive.

Table 2

List of highly conserved motifs predicted by MotViz only

PDB ID	Motif sequence	Cv	Location
3ML6	VFGDLS	0.733	434-440
3ML6	RNELFLDV	0.9875	456-464
3ML6	KVVIK	0.85	558-563
3ML6	GMKESQISAEIE	0.99104	964-976
2XQL	AGNAARDNK	0.883	232-241
2XQL	IIPRHLQLA	0.9013	244-253
2XQL	LGKVTIAQG	0.677	263-272
2WOM	KIEELRQH	0.9893	200-208
2WOM	GIRK	0.925	554-558
1KHY	PSQD	0.812	203-207
1KHY	AGQLRTDIN	0.772	315-324
1KHY	ISS	0.85	359-362
1KHY	LSALLN	0.733	426-432
1KHY	GSVS	0.675	435-439
1KHY	NITQAIEQ	0.656	526-534
3MDF	KHRGFAF	0.871	127-134
3MDF	RTIRV	0.679	157-162
3MEY	TEALRFPPVM	0.655	44-54
3MEY	YLKSFPNL	0.9832	63-71
3MEY	CFRREPSKHL	0.825	139-149
3MEY	VDYASDPFF	0.95	192-201
3MEY	LKFELLIPL	0.838	216-225

Note: Highly-conserved motifs predicted from 9 protein families by MotViz were verified by Cv.

A comprehensive comparison of MotViz predicted motifs was made by combining the motif-prediction data of listed protein families. Individual motifs were predicted by afore-mentioned tools, respectively and compared to those predicted by MotViz. Subsequently, the representative sequence motifs were plotted graphically (Figure S1). The pattern clearly indicate that majority of motifs predicted by MotViz are missed by the listed tools. Furthermore, we plotted the number of motifs predicted and found that more motifs were predicted by MotViz than any listed tools for all 9 protein families. For example, MotViz predicted 9 motifs in Hox protein family, and 5-6 motifs were also predicted by MEME, Prosite, BLOCKS, PRINTS and Pfam, while PATTERN predicted only 3 motifs. In addition, MotViz predicted 36 motifs in histone protein family, out of which MEME, BLOCKS and Pfam predicted 21, 16 and 12 motifs, respectively, while the remaining tools predicted less than 10 motifs (5 by Prosite, 4 by PRINTS, 3 by ProDom and 1 by PATTERN). Furthermore, MotViz predicted 16 motifs in ribosomal family proteins, while 7 motifs were predicted by MEME. Prosite and BLOCKS only predicted 1 motif while PATTERN, Pfam, ProDom and PRINTS did not predict any motif (Figure 3). These data illustrate the MotViz algorithm efficiency is significantly high, compared to other tools.

Figure 3

Number of motifs predicted by MotViz for protein families selected randomly. X-axis indicates protein families examined, Y-axis represents number of motifs predicted using different tools while Z-axis represents the motif databases. The performance of MotViz is superior to other tools examined, although the performance of MEME ( is very close to that of MotViz.

Discussion

Combination of conserved sequence segments in protein with visualization in a single click offers a remarkable addition that may result in gathering the in-depth details about functional elaboration of protein structure elements. These details would largely help in dissecting functional information of short segments in parallel to their structural localization. Motifs are important for predicting protein structure or behavior and to categorize anonymous proteins into appropriate families (. Presence of similar motifs can be annotated by analyzing multiple structures, which may provide a unique mode for determining the novel links and characterizing the functional relevance of unknown proteins.

MotViz plug-in preferably isolates the conserved motifs by aligning the query sequence to a similar dataset by performing the PSI-BLAST followed by MSA using ClustalW in order to dissect the narrow regions of sequence similarities. MotViz motif prediction efficiency is much higher than other tools due to the implementation of increased conservation rate and accuracy in picking by applying filters in two consecutive steps. In the first step, PSI-BLAST at UniProt ( is launched resulting in retrieving the 50 most related sequences and later in the second step, these sequences are further narrowed down by performing MSA using the ClustalW ( algorithm. The subsequent identification of conserved sites was further enhanced by selecting the alignment regions of consecutive amino acids with less than three gaps.

Conclusion

MotViz tool provides better usability due to many reasons. Firstly, it generates the query sequence from 3D coordinates, searches similar proteins (with maximum E-value of 10−50) and then aligns them to find the conserved residues in a single effort. Secondly, these residues are further annotated by online motif prediction tools and are verified before appearing for visualization in protein structure. Thirdly, MotViz keeps features that have already been embedded in Jmol intact (. Lastly, addition of database, incorporation of conservation score calculation, graphical representations and individual sequence to structural compatibility differentiate MotViz. The user may visualize each motif individually within the corresponding structure and examine structural details of the highlighted motif. The graphical output of physiochemical properties and conservation score in parallel to motif prediction make it a unique platform for user.

Authors’ contributions

MSN designed the tool and carried out the data analysis and comparison assays. SR conceived and supervised the project, developed algorithm and methods, and prepared the manuscript. All authors read and approved the final manuscript.

Competing interests

The authors have no competing interests to declare.