This article has Open Peer Review reports available.
SequenceCEROSENE: a computational method and web server to visualize spatial residue neighborhoods at the sequence level
© Heinke et al. 2016
Received: 23 September 2015
Accepted: 17 January 2016
Published: 27 January 2016
To understand the molecular function of biopolymers, studying their structural characteristics is of central importance. Graphics programs are often utilized to conceive these properties, but with the increasing number of available structures in databases or structure models produced by automated modeling frameworks this process requires assistance from tools that allow automated structure visualization. In this paper a web server and its underlying method for generating graphical sequence representations of molecular structures is presented.
The method, called SequenceCEROSENE (color encoding of residues obtained by spatial neighborhood embedding), retrieves the sequence of each amino acid or nucleotide chain in a given structure and produces a color coding for each residue based on three-dimensional structure information. From this, color-highlighted
sequences are obtained, where residue coloring represent three-dimensional residue locations in the structure. This color encoding thus provides a one-dimensional representation, from which spatial interactions, proximity and relations between residues or entire chains can be deduced quickly and solely from color similarity. Furthermore, additional heteroatoms and chemical compounds bound to the structure, like ligands or coenzymes, are processed and reported as well.
To provide free access to SequenceCEROSENE, a web server has been implemented that allows generating color codings for structures deposited in the Protein Data Bank or structure models uploaded by the user. Besides retrieving visualizations in popular graphic formats, underlying raw data can be downloaded as well. In addition, the server provides user interactivity with generated visualizations and the three-dimensional structure in question.
Color encoded sequences generated by SequenceCEROSENE can aid to quickly perceive the general characteristics of a structure of interest (or entire sets of complexes), thus supporting the researcher in the initial phase of structure-based studies. In this respect, the web server can be a valuable tool, as users are allowed to process multiple structures, quickly switch between results, and interact with generated visualizations in an intuitive manner.
The SequenceCEROSENE web server is available at https://biosciences.hs-mittweida.de/seqcerosene.
To understand the role governed by a protein or biopolymer of interest, the study of its structure is often of great importance. In the process of investigation, schematic representations of interactions between structural elements provide simplified but informative visualizations that help the researcher to understand and illustrate but also to communicate essential molecular characteristics . Such characteristics range from the general topological arrangement of entire domains and secondary structure elements to the level of residue interaction networks.
One can distinguish between two categories of representations. The first category includes schematics which we here refer to as 2.5-dimensional representations. These diagrammatic representations illustrate spatial neighborhood and arrangements of structural elements. Given atomic coordinate data as input, locations of the elements in two-dimensional space are computed that resemble the three-dimensional topology as closely as possible. For example, TopDraw , HERA , Pro-origami  and TOPS  are tools for generating such representations. Furthermore, precomputed topology diagrams of solved structures can be found at the PDBsum database . PROTTER  is a recent addition to this class of tools and is specifically tailored for visualizing α-helical membrane protein topology. In contrast to aforementioned tools, PROTTER derives topology information only from predictions made from sequence, whereas available structure information is not considered. On the level of residue interactions, intuitive and clear representations are restricted to subsets of residues and their associated interactions, due to the vast amount of interactions present in a structure. Approaches such as RING  propose 2.5-dimensional graph-based representations of residue/interaction sets. In this respect, the Protein Graph Repository provides access to nearly 190,000 graphs generated from about 94,000 protein structures .
Although diagrammatic representations feature visual clarity, the process of low-dimensional mapping achieved by reducing dimensionality yields a drawback not obvious to the user: dimensionality reduction is generally accompanied by information loss and unavoidable morphing effects, both depending on structural complexity. Eventually, information loss and morphing effects reduce mapping quality. Thus, resemblance of structure topology between the actual three-dimensional structure and its corresponding 2.5-dimensional visualization is not necessarily of the same quality for any given protein structure.
The second category of techniques avoid this problem by employing bijective projections. Here, color maps illustrating residue-residue adjacencies or distances are straightforward visualizations and can be produced using a number of available tools, such as CMView  and CMA . However, compared to 2.5-dimensional representations such maps cannot achieve the same degree of clarity, which makes them difficult to interpret intuitively. Thus, these visualizations are of greater use if interactivity with 3D structure viewers is implemented (such as realized by CMView).
In this paper, an effective and straightforward approach for generating intuitive representations of spatial residue neighborhood is proposed, which we refer to as SequenceCEROSENE (color encoding of residues obtained by spatial neighborhood embedding). SequenceCEROSENE produces color-highlighted sequences of amino acid, DNA, and RNA chains present in the query structure, including bound chemical compounds and heteroatoms. Using a straightforward transformation of the atomic coordinates into RGB color space, color-highlighting of individual residues and compounds corresponds to their location in the three-dimensional space. Thus, similarity of color present in color encoded sequences between residues or compounds illustrates spatial proximity. The overall presentation of the visualization resembles a colorized set of FASTA formatted sequences, where bound compounds/heteroatoms are reported individually. SequenceCEROSENE aims at producing representations containing as much structural information as possible while allowing intuitive interpretation. In fact, by decoding the RGB content of individual residues, the three-dimensional structure of the complex can be well approximated.
In the following sections, the SequenceCEROSENE method is introduced and demonstrated, followed by a brief presentation of the web server.
Method and discussion
Finally, SequenceCEROSENE provides color-encoded illustrations of corresponding amino acid and nucleotide sequences. Bound compounds and heteroatoms are visualized separately. This representation is both intuitive and powerful, since residue neighborhood can be directly perceived by color similarity and the structural information between representative coordinates is maintained. In theory, RGB color contents for each residue in the generated representation can be decoded in order to restore the actual three-dimensional coordinates. Information loss only occurs during the representative coordinate determination process, as centroid coordinates and single atoms are considered.
As an example, the structure of transcription promoter CAP (catabolite activator protein, PDB-Id: 1cgp ) and the corresponding sequence representation generated by SequenceCEROSENE are shown in Fig. 1 b and c, respectively.
A web server has been implemented to provide free access to SequenceCEROSENE for the scientific community. Upon entering a valid Protein Data Bank (PDB) identifier  or uploading single or multiple structures in PDB format, the server computes color encodings and presents these to the user interactively. As discussed in the Method and discussion section, the set of residue representative coordinates needs to be determined prior to processing. By default, C β atoms (C α atoms for glycines) are selected for amino acids. However, as an advanced option, the user can also choose between using C α atoms, terminal heavy side chain atoms, or side chain centroid coordinates. For nucleotide residues and compounds, the centroid coordinate is always considered.
Furthermore, multiple structures can be processed in one submission. Therefore, the user has to provide an archive containing the PDB files in question. All commonly used archive file formats are supported. If a PDB-Id is provided, the web server retrieves the data from a local weekly-updated snapshot of the PDB. If no local PDB file for a user query is present, a fallback option is implemented to automatically retrieve the data from the PDB using RESTful web services.
The utilization of dedicated 3D structure visualization programs is the method of choice for studying and understanding general spatial characteristics of biopolymer structural data. Techniques for generating "condensed" visualizations portraying such characteristics are widely used in the process, as they support the researcher in constructing a "mental image" of the structure in question.
Generated visualizations can be classified into two general categories: 2.5-dimensional diagrammatic representations of protein topology (the general spatial arrangements of structural elements) and visualizations on residue-residue interaction level (such as residue distance or interaction matrices). As discussed, both visualization categories are limited by trade-offs between visual clarity and amount of information presented.
SequenceCEROSENE aims at filling the gap between techniques of both categories by providing an intuitive visualization on sequence level while keeping information loss minimal. Furthermore, in contrast to most methods, structure data of protein, DNA, RNA and complexes thereof can be processed, making this technique applicable to a variety of biopolymers.
Analogous to available automated structure visualization techniques, the main intention of the method is to provide visual guidance to researchers and students in the initial phase of studying structural data, especially if the number of structures is large. The implemented web server allows processing, downloading, and interacting with generated visualizations, however future implementations of SequenceCEROSENE as plug-ins or add-ins for common visualization programs, such as PyMOL  and VMD , could improve accessibility. Considering the simplicity of the methodological idea, we hope that researchers are inspired to adapt this technique, thus giving rise to implementations that are specifically tailored toward their own research and analysis pipelines.
Availability and Requirements
The server can be accessed using common web browsers. For displaying protein structures, PV requires WebGL to be supported, which, however, is the case for the majority of available browsers. To enable these WebGL-based features for unsupported browsers, manual installation of browser-specific WebGL add-ons and libraries is required. No further requirements are necessary. The Sequence- CEROSENE web server is freely available to the scientific community.
The authors thank the Free State of Saxony and the Saxon Ministry of Science and the Fine Arts for funding.
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Hutchinson EG, Thornton JM. HERA–a program to draw schematic diagrams of protein secondary structures. Proteins. 1990; 8(3):203–12. doi:http://dx.doi.org/10.1002/prot.340080303.View ArticlePubMedGoogle Scholar
- Bond CS. TopDraw: a sketchpad for protein structure topology cartoons. Bioinformatics. 2003; 19(2):311–2.View ArticlePubMedGoogle Scholar
- Stivala A, Wybrow M, Wirth A, Whisstock JC, Stuckey PJ. Automatic generation of protein structure cartoons with Pro-origami. Bioinformatics. 2011; 27(23):3315–6. doi:http://dx.doi.org/10.1093/bioinformatics/btr575.View ArticlePubMedGoogle Scholar
- Westhead DR, Slidel TW, Flores TP, Thornton JM. Protein structural topology: Automated analysis and diagrammatic representation. Protein Sci. 1999; 8(4):897–904. doi:http://dx.doi.org/10.1110/ps.8.4.897.PubMed CentralView ArticlePubMedGoogle Scholar
- de Beer TAP, Berka K, Thornton JM, Laskowski RA. PDBsum additions. Nucleic Acids Res. 2014; 42(Database issue):292–6. doi:http://dx.doi.org/10.1093/nar/gkt940.View ArticleGoogle Scholar
- Omasits U, Ahrens CH, Müiller S, Wollscheid B. Protter: interactive protein feature visualization and integration with experimental proteomic data. Bioinformatics. 2014; 30(6):884–6. doi:http://dx.doi.org/10.1093/bioinformatics/btt607.
- Martin AJM, Vidotto M, Boscariol F, Di Domenico T, Walsh I, Tosatto SCE. RING: networking interacting residues, evolutionary information and energetics in protein structures. Bioinformatics. 2011; 27(14):2003–5. doi:http://dx.doi.org/10.1093/bioinformatics/btr191.View ArticlePubMedGoogle Scholar
- Dhifli W, Diallo B. PGR: A Novel Graph Repository of Protein 3D-Structures. J Data Mining Genomics Proteomics. 2015; 6(2):1–4.Google Scholar
- Vehlow C, Stehr H, Winkelmann M, Duarte JM, Petzold L, Dinse J, et al.CMView: interactive contact map visualization and analysis. Bioinformatics. 2011; 27(11):1573–4. doi:http://dx.doi.org/10.1093/bioinformatics/btr163.View ArticlePubMedGoogle Scholar
- Sobolev V, Eyal E, Gerzon S, Potapov V, Babor M, Prilusky J, et al.SPACE: a suite of tools for protein structure prediction and analysis based on complementarity and environment. Nucleic Acids Res. 2005; 33(Web Server issue):39–43. doi:http://dx.doi.org/10.1093/nar/gki398.View ArticleGoogle Scholar
- Berman HM, Kleywegt GJ, Nakamura H, Markley JL. The Protein Data Bank archive as an open data resource. J Comput Aided Mol Des. 2014; 28(10):1009–14. doi:http://dx.doi.org/10.1007/s10822-014-9770-y.PubMed CentralView ArticlePubMedGoogle Scholar
- Schultz SC, Shields GC, Steitz TA. Crystal structure of a CAP-DNA complex: the DNA is bent by 90 degrees. Science. 1991; 253(5023):1001–7.View ArticlePubMedGoogle Scholar
- The PyMOL Molecular Graphics System. CA: DeLano Scientific. 2002. Available online at https://www.pymol.org/. Accessed 22 Sept 2015.
- Humphrey W, Dalke A, Schulten K. VMD: visual molecular dynamics. J Mol Graph. 1996; 14(1):8233–78.View ArticleGoogle Scholar
- Biasini M. PV - WebGL-based protein viewer. 2014. doi:http://dx.doi.org/10.5281/zenodo.12620.