PDB-Explorer: a web-based interactive map of the protein data bank in shape space
BMC Bioinformatics volume 16, Article number: 339 (2015)
The RCSB Protein Data Bank (PDB) provides public access to experimentally determined 3D-structures of biological macromolecules (proteins, peptides and nucleic acids). While various tools are available to explore the PDB, options to access the global structural diversity of the entire PDB and to perceive relationships between PDB structures remain very limited.
A 136-dimensional atom pair 3D-fingerprint for proteins (3DP) counting categorized atom pairs at increasing through-space distances was designed to represent the molecular shape of PDB-entries. Nearest neighbor searches examples were reported exemplifying the ability of 3DP-similarity to identify closely related biomolecules from small peptides to enzyme and large multiprotein complexes such as virus particles. The principle component analysis was used to obtain the visualization of PDB in 3DP-space.
A chemical space classification of PDB based on molecular shape was obtained using a new atom-pair 3D-fingerprint for proteins and implemented in a web-based database exploration tool comprising an interactive color-coded map of the PDB chemical space and a nearest neighbor search tool. The PDB-Explorer website is freely available at www.cheminfo.org/pdbexplorer and represents an unprecedented opportunity to interactively visualize and explore the structural diversity of the PDB.
One of the striking features of biomolecules is their extremely large diversity spanning from small organic molecules such as metabolites and drugs to large supramolecular complexes such as the ribosome or viral particles. A vast amount of knowledge about these biomolecules has been collected in various public databases, in particular the Protein Data Bank (PDB) which collects over 100,000 different 3-dimensional (3D) structures of biological macromolecules determined by X-ray crystallography, NMR spectroscopy and electron microscopy [1–4]. Despite of this vast amount of information, the overall structural diversity available in the PDB is difficult to perceive. Indeed various tools are available to search the PDB for analogs of specific proteins according to similarities in evolutionary history, sequences, secondary structures and subdomains [5–12]. In the case of 3D-SURFER [13, 14] the PDB is classified according to similarities in protein surface allowing to search for shape analogs among PDB-entries. The CATH [7, 8] contains a subset of the PDB which is directly visible via an overview interface using a hierarchical classification by structural domains. However none of these interfaces provides a direct, global yet comprehensive overview of the PDB irrespective of a specific query, which would be desirable to understand its overall contents.
Herein we report a new exploration tool for the PDB called PDB-Explorer which addresses the need for a global perception of the database by giving direct access to all PDB-entries via an interactive color-coded map representing its entire contents in molecular shape space. This application follows the principle of our recently reported MQN-mapplet and SMIfp-mapplet applications designed to visualize the chemical space of small organic molecules [15–18]. Each individual PDB-entry is placed on the map of the PDB-Explorer according to its 3D-shape as encoded by a new fingerprint called 3DP featuring a generalized version of our recently reported 3D-atom pair fingerprints for small molecules [19, 20]. The PDB-Explorer provides an unprecedented global view of the PDB allowing a detailed exploration of its entire content in a curiosity-driven manner with or without specific queries. This tool is freely available at www.cheminfo.org/pdbexplorer and should greatly facilitate the perception and understanding of the overall diversity of proteins and biological assemblies available in the PDB.
The X-ray structures in PDB-Explorer were downloaded from http://www.rcsb.org. The water molecules and hydrogen atoms of each PDB molecule were removed at the beginning. If the number of atoms assigned as “HETATM” occupied more than 20 % of total heavy atom count, this molecule was not included in the database.
3D protein atom-pair fingerprint (3DP)
City-Block Distance (CBD): The CBD between two molecules, A and B, with 136 dimensions of 3DP is calculated as shown in Eq. 3:
The 3DP-fingerprint was transformed into a 200-dimensional 3DP-similarity fingerprint (3DPs) by calculating the similarity value of each of the PDB-entry in the database relative to 200 randomly selected PDB molecules. Similarity values S3DP were calculated from the city-block distance CBD3DP as shown in Eq. 4, with X = median city-block distance observed for distribution of CBDs computed for 1 M random pairs of molecules within the database:
Principal component analysis (PCA)
The PCA calculation used the source code from Java program developed based on the tutorial of Lindsay I. Smith (http://www.cs.otago.ac.nz/cosc453/student_tutorials/principal_components.pdf). The Java source code uses mathematical functions from the JSci (A science API for Java: http://jsci.sourceforge.net/) library.
PDB-maps generation and color-coding
PCA of the 3DP-similarity fingerprint space of PDB database was carried out, and the PC-1 and PC-2 values were computed for each molecule in database. The PC-1 and PC-2 values were binned onto the 2D-grid of size 300 × 300 using the same absolute bin size on the PC-1 and PC-2. The largest (PCmax) and smallest (PCmin) PC values appearing in the PC-1 or PC-2 values were used to define the value range ΔPC = PCmax − PCmin and set the binning scale as ΔPC/300. Afterwards, each molecule was assigned to a point (or bin) on this 2D-grid.
The following molecular properties were computed at each grid point: HAC, percentage of positive, negative, and hydrophobic atoms, and molecular volume occupancy (mvo). The molecular volume occupancy (mvo) describing the compactness of a molecule was computed with the following formula (Eq. 5):
Each of the grid point was color coded according to the average and standard deviation of the molecular property at that grid point using the Hue-Saturation-Lightness (HSL) color space, with the Hue value (blue-cyan-green-yellow-red-magenta) representing the average value and the Saturation (fading to grey) representing the standard deviation.
The calculation of normalized principal moment of inertia (nPMI1, nPMI2) was implemented by an in-house Java program written based on the work from Sauer and Schwarz . The position in the (nPMI1, nPMI2) triangle was color coded using the RGB color space assigning the distance to each triangle summit as the relative R, G, and B values.
Results and discussion
Design of the protein shape fingerprint 3DP
In view of a global analysis of the PDB we set out to identify a fingerprint encoding the 3D-structure of biomacromolecules since it is known to play an important role in their biological function, their interaction with other molecules, [22–25] and their evolution [26–29]. While 3D-SURFER uses a 121-dimensional scalar fingerprint to encode the shape of the molecular surface of proteins using 3D Zernike descriptors, [13, 14] we searched for a simpler yet more detailed encoding considering not only the shape of the molecular surface, but also the atom types and the internal structure of the protein. Inspired by the concept of atom-pair fingerprints proposed by Carhart,  Sheridan  and Schneider  to encode pharmacophores in small organic molecules, we recently reported a detailed analysis establishing the suitability of atom-pair fingerprints for 3D-shape and pharmacophore similarity searches in very large databases such as ChEMBL  and ZINC  using both topological distances read from the 2D-structures  and through-space distance read from the 3D-structures . While topological distances cannot be extracted easily from the structures of macromolecules or even do not exist within non-covalent assemblies, an atom-pair analysis of macromolecules should be possible using through-space distances which can be directly computed from the atomic coordinates available in each pdb file. An atom pair 3D-fingerprint, here called 3DP, was therefore envisioned to encode 3D-structures in the PDB.
The 3DP fingerprint was designed to encode the 3D structure of proteins and other biomacromolecules in PDB considering the biological assembly as defined by the authors in each entry. Following our previous approach for small molecules, each atom pair distance was converted to a Gaussian centered on the atom pair distance with a width of 18 % of the distance itself, and the Gaussian was sampled at 34 values between 1.45 Å and 342.53 Å covering all atom pair distances present in PDB (Fig. 1a). For each of the 34 values increments were summed across all atom pairs. The 34 resulting sums were normalized to the heavy atom count to the power of 1.5 (HAC1.5) to reduce sensitivity to size (Fig. 1b). The 34 values were computed separately for four different atom categories, using the corresponding sum of category atoms as HAC for normalization. These four categories considered as important for macromolecular properties comprised 1) all atoms; 2) positive charges: lysine and arginine side chains; 3) negative charges: aspartate and glutamate side chains, phosphate groups; 4) hydrophobic atoms, defined here as carbon atoms with covalent bonds only to C or H atoms. Due to the smaller number of charged and hydrophobic atoms compared to all atoms the bit values in the corresponding categories were on average much smaller and were therefore multiplied by 5 for charged atom categories and by 2 for the hydrophobic atom category. The final bit values were expressed in percent and rounded to the integer value. The four combined sets of normalized, rounded values formed the 136-dimensional 3DP fingerprint.
Because most PDB-entries contain thousands of atoms which would require explicit calculation of many millions of atom pair distances (Fig. 1d), the 3DP calculation was simplified by computing the fingerprint using formal “sum atoms”. A maximum of 1728 sum atoms (resulting in a maximum of 1.5 million sum atom pairs) were defined for each PDB-entry by fitting its biological unit into a 12×12×12 grid, using the average coordinates of each atom category within each box as the sum atom coordinates and the number of category atoms within this box as its relative weight. 3DP bit values computed with sum atoms were essentially indistinguishable from those computed with actual atoms.
The 3DP calculation was performed 91,223 X-ray structures downloaded from the PDB in September 2014 (Additional file 1: Supplement 1), considering in each case the biological assembly as defined by the authors . The 3DP fingerprint had a high resolution, with 99.99 % of PDB molecules having unique 3DP fingerprints. The average bit values peaked at 33.76 Å (B20, B122) for all atoms and hydrophobic atom pairs and at 39.83 Å (B55, B89) for positive and negative charge atom pairs, while the corresponding standard deviation covered almost the entire bit value range (Fig. 1c).
3DP encodes protein conformations
The 3DP fingerprint had a remarkable ability to precisely encode the shape of proteins, as evidenced by investigating the correlation between 3DP similarity and the root mean square deviation (RMSD) [37, 38] in different conformers of the same molecule, as presented with the following three test cases. As a first test case a random coil 24-mer peptide was obtained from PDB-entry 4GOF and simulated using the ff99SB force field from the AMBER12 package,  with 1 ns simulated annealing followed by 50 ns free simulation in water solvent, generating a large variety of conformers. The simulation was repeated 50 times, and in each case the last structure was selected as reference. Its conformer analogs were defined as all conformers with RMSD lower than 2 Å to that reference, which comprised up to 136 conformers that were always the last series of frames in the MD run. Retrieval of these conformers from the 2000 conformers of the corresponding trajectory was then performed by sorting using the city-block distance CBD3DP as similarity measure. The recovery was excellent, with an area under the curve (AUC) of 95.85 % for average receiver operator characteristic (ROC) curve, showing that 3DP indeed provided a good encoding of molecular shape in terms of similarity searching (Fig. 2a). All RMSD analogs were found within CBD3DP < 519, while the average distance of these RMSD analogs from the query conformer was 264 and the largest distance was 2164 (Fig. 2b).
The correlation between 3DP similarity and RMSD was tested in a second case for a larger protein by considering ten domain movement frames for glutamine binding protein (1762 atoms) from the Protein Motion Database [40, 41]. These ten conformers represent different conformations of the binding domain spanning between open (PDB-entry 1GGG, purple structure in Fig. 3a)  and closed state (PDB-entry 1WDN, red structure in Fig. 3a) . A good correlation was observed between 3DP similarity and RMSD (calculated with all-atom alignment by Maestro 8.5 ), in particular when considering conformer pairs from one extreme of the movement range (purple and red lines in Fig. 3b). Representation of the bit value changes showed that 3DP perceived the conformational change at the level of each of the four different atom type categories (Fig. 3c).
The sensitivity of 3DP similarity to slight changes in protein conformation was evaluated in a third test case considering CDK2 (cyclin-dependent kinase 2), an important player in cell cycle regulation . The activity of CDK2 is induced by conformational changes occurring upon ligand binding, in particular at the level of the large T-loop . A total of 245 X-ray structures of CDK2 monomers bound with various small molecule inhibitors are available in the PDB (Additional file 1: Supplement 2). These CDK2 structures are nearly identical in sequence (>97 % identity), but show different conformation at the level of the T-loop due to an induced-fit adaptation of the protein around the various ligands. Although these 245 CDK2 structures were almost identical in terms of their global molecular shape, the 3DP analysis readily distinguished between the different T-loop conformers, with pairwise CBD3DP distances between 75 and 1415 CBD3DP units and maximum population around 300 ≤ CBD3DP ≤ 500 (Fig. 4a). Three different pairs were analyzed closer in terms of differences in their 3DP values (Fig. 4b) and their overall structures superimposed by all-atom alignment (Fig. 4c) . In terms of the 3DP values differences were visible in all four atom pair categories. In terms of the protein structures, the closest pair 3QZH/3ROY at CBD3DP = 89 and the intermediate pair 3QRT/2C5Y at CBD3DP = 424 only significantly differed at the level of the T-loop. For the 3QQH/4EZ7 pair at CBD3DP = 1021, deviations were not only visible at the level of the T-loop but also in the α-helix and β-sheet at the top of the structures. The CDK2 analysis clearly showed that 3DP was able to perceive small conformational differences residing in a loop orientation between otherwise highly shape similar proteins.
3DP analysis of the RSCB protein databank
To facilitate access to PDB-entries a graphical user interface was created based on an interactive color-coded map representing the 3DP chemical space using the Mapplet principle previously reported for small molecule databases, by which database entries are displayed in a visualization window as the mouse cursor moves on the map, and connects to a fingerprint similarity search window to allow nearest neighbor searches [15–18]. Although 78 % of the data variability was represented in (PC1, PC2)-plane obtained by principal component analysis (PCA) of the 3DP dataset, this direct PCA map contained many scattered pixels with uneven occupancies and was not suitable as an interface (data not shown). We therefore generated an alternative representation of 3DP by similarity mapping, which produces more compact and evenly populated maps with a fairly good rendering for various chemical spaces [18, 48, 49].
To create the similarity map, 3DP similarity values relative to 200 randomly selected reference molecules from the PDB were calculated, and the (PC1, PC2)-plane obtained by PCA of the resulting 200-dimensional 3DP-similarity fingerprint was represented. This (PC1,PC2)-plane, called similarity map, covered 98 % of the similarity fingerprint data variability and had a compact, comet-like shape without peripheral pixels. Furthermore the map had a relatively even occupancy of pixels which was highly suitable as a representation of the PDB in 3DP space (Fig. 5a). PDB-entries were distributed on the map according to their size along the edge of the comet (Fig. 5b). The position on the map was influenced by the molecular volume occupancy (mvo, Å3/atom), a property defined here as the volume of the sphere with a radius corresponding to the average distance of all atoms to the center of gravity, divided by the total number of atoms, with the more compact, globular structure present at the comet edge and less compact structures at the center (Fig. 5c). The map also separated PDB-entries with different fractions of positively charged atoms, negatively charged atoms, and hydrophobic atoms, whereby the distribution pattern for the fraction of hydrophobic and positively charged atoms was somewhat similar (Fig. 5d–f). Color-coding by the normalized principal moments of inertia vector (nPMI1, nPMI2) distinguishing between rod-like, disc-like and sphere-like shapes,  showed that rod-like PDB-entries were located at the comet center, while the more spherical entries were distributed around the rest of the map (Fig. 5g).
The color-coded similarity maps were combined with the 3DP-similarity search tool into a web-based application called “PDB-Explorer” for interactive visualization and 3DP-similarity search through the entire PDB (Fig. 6). The website uses the same principles as our previously reported MQN-mapplet, [15–17] and is based on the open-source project visualizer (http://github.com/NPellet/visualizer) that was already successfully used for another cheminformatic project . The PDB-Explorer consists of a main window to browse through various color-coded rendering of the 3DP-similarity map. The average PDB molecule in each pixel is shown in the 3D-viewer of PV,  and the “Show Bin” table details the contents of any selected pixel on the map from which each PDB-entry, shown as ribbon image, can be inspected closer as 3D-model by opening a secondary JSmol window. The “Locate Molecule” function allows locating any given PDB-entries on the map and the “Similarity Search” option offers nearest neighbors of any PDB-entry in the 3DP space within seconds. The application also contains uploading function to locate and search nearest neighbors of any structure represented as pdb file. All computations are performed locally by the browser of the user machine. The color-coded similarity maps and the database for similarity search are updated every day by adding new entries in the PDB. The detailed PDB-Explorer functions are described in the HELP page.
Exploring PDB using the PDB-Explorer
The PDB-Explorer allows the rapid analysis and overview of the entire PDB as well as detailed searches around selected PDB-entries using 3DP similarity as a guiding principle. Its use is exemplified with three case studies detailed below which further demonstrate the remarkable ability of 3DP to classify proteins according to their 3D structure.
The case of the conotoxins, a group of 10 to 30 residues neurotoxic peptides containing multiple intramolecular disulfide bridges and a variety of secondary structures, [52–54] provides an example of the smallest molecules listed in the PDB. Alpha-conotoxin PnIA (PDB-entry 1PEN) with 17 residues serves as reference molecule (Fig. 7a). The closest 3DP analog of 1PEN retrieved by 3DP similarity is 1AKG, which belongs to the similar α-conotoxin family and is also retrieved as nearest neighbor by BLAST search with 83 % sequence similarity (Fig. 7b) . By contrast the second analog found by 3DP similarity is a short fibril peptide complex (Fig. 7c), which has a completely different secondary structure and sequence, but a similar overall size and shape. The 3rd (Fig. 7d) and 4th analog (Fig. 7e) are both from the enterotoxin family featuring a conotoxin-like shape containing double disulfide bridges, again without significant sequence similarity to the reference. Three further conotoxins missed by sequence similarity appear in the 3DP nearest neighbor search at rank 7 (1NOT, Fig. 7f), rank 9 (1HJE, Fig. 7g) and rank 33 (4TTL, Fig. 7h). In the latter case the bit value profile shows that 4TTL is the only sequence containing negatively charged residues, while all other cases do not contain any charged residues (Fig. 7i).
3DP similarity also retrieves shape analogs of molecular assemblies, as exemplified here with triose phosphate isomerase (TIM), a homodimeric enzyme which contains around 250 amino acids in the monomer . The wild type TIM dimer from yeast (PDB-entry 1YPI) is selected as the reference. The top 100 neighbors of 1YPI are also protein dimers, with essentially all neighbors of 1YPI up to CBD3DP = 800 being TIM enzymes. All TIM structures in PDB (Additional file 1: Supplement 3) in fact occur within a CBD3DP distance of 1200 except 13 TIM enzymes ranked beyond rank 300 (Fig. 8a). The differences in 3DP values are illustrated for the nearest neighbor 4FF7, the rank 10 analog 8TIM at distance CBD3DP = 493 representative of the bulk of TIM structures, and 4GNJ found at rank 312 and CBD3DP = 1181 representative of the more distant group around CBD3DP = 1100 (Fig. 8b). The nearest neighbor of 1YPI, 4FF7, is also a yeast TIM with 99 % sequence similarity and nearly identical structure. PDB-entry 8TIM at rank 10 is a TIM dimer from a different organism (gallus gallus) with only 52 % of sequence similarity. In both cases the shape and fold differences to the reference are quite small, and the larger distance of 8TIM stems from differences in positively charged atoms that are not directly visible in the protein shape (Fig. 8b/c). The difference in the number and position of charged atoms also explains the further distance from the reference of 4GNJ, a TIM dimer from leishmania siamensis.
The PDB-Explorer furthermore provides interesting insight at the level of very large protein assemblies, as exemplified here for the case of a virus capsid. Starting with the subviral particle of the bursal disease virus capsid (PDB-entry 2GSY)  containing 27,120 residues (Fig. 9a), four analogs are readily identified within the distance range CBD3DP < 10,000. The nearest neighbor is 3FBM (Fig. 9b) which is a mutant protein of the query 2GSY,  and the second and third closest structures 1WCD and 2DF7 are capsids from the same bursal disease virus (Fig. 9c, d). The fourth analog (PDB-entry 3IDE) is the protein coat of Infectious Pancreatic Necrosis Virus (IPNV)  which has similar spikes and forms a similar icosahedral capsid organization as the capsid of infectious bursal disease virus (Fig. 9e) . Interestingly the parvovirus capsid protein 1DNV (Fig. 9f), which is very close in size to the reference 2GSY, only appears much further away in CBD3DP space because its spherical shape does not feature spikes, which significantly impact its 3DP fingerprint profile (Fig. 9g).
3DP space automatically clusters proteins from the same superfamily in tight groups provided that they are in a similar size range. This property is illustrated here for the recovery of each family from a benchmark dataset of 150 CATH superfamilies, each containing between 12 and 1533 PDB entries [7, 8], considering proteins in the majority size range (Additional file 1: Supplement 4). The ROC (receiver operator characteristic) curves for recovering superfamily members from the PDB-entry closest to all same superfamily members in 3DP-space give very high AUC (area under the curve, all above 90 %) and generally high EF0.1% values (enrichment factor at 0.1 % database coverage in the range 29–1000, Fig. 10a). Most superfamilies also appear as tight groups on the PDB-map (Fig. 10b). The clustering of protein superfamilies in 3DP-space reflects the fact that the definition of these families considers similarities in folds in addition to function.
Comparing 3DP with protein structure alignment tools
The performance of 3DP was compared with three protein structure alignment tools, Fr-TMalign [60, 61], SPalign  and MATT . Fr-TMalign is applied to pairwise structure alignment based on fragment similarity, while SPalign is designed for detecting proteins with similar fold and similar function of DNA or RNA binding. MATT (Multiple Alignment with Translations and Twists) is a program to align multiple protein structures allowing certain flexibility between fragments. These alignment methods are computationally intensive and therefore only applicable to a limited number of comparisons. They are size independent and focus on backbone alignment resulting in a focus on secondary structures. On the other hand 3DP is size dependent and considers all protein atoms indiscriminately, resulting in a sensitivity to the overall shape rather than to secondary structures. Remarkably, 3DP allows an essentially instantaneous comparison with the entire PDB when using the PDB-Explorer website.
In a first comparative study, all pairwise alignment scores were computed for the ten domain movement frames for glutamine binding protein and compared with CBD3DP. The data showed that CBD3DP, which performs comparably to overall structure RMSD as dicussed above (Fig. 3), had similar trend but higher sensitivity to conformer differences than Fr-TMalign or SPalign. MATT was highly sensitive to small conformational changes and classified all non-identical conformer pairs as low scoring (Fig. 10c).
A further comparison of 3DP with structure alignment tools was carried out for 50 homologous CDK2 proteins (Fig. 4) and 50 non-CDK2 decoy proteins. The 100 structures were in similar size range from 2100–2600 heavy atoms, and decoys consisted of non-homologous proteins with pairwise sequence identity lower than 30 %. Alignment scores and 3DP distances were computed for the 1225 CDK2 pairs and the 2500 CDK2-decoy pairs (Fig. 10d). 3DP made a relatively clear cut between CDK2 pairs and CDK2-decoy cross-pairs at CBD3DP = 750, with all CDK2 pairs found within the range CBD3DP < 1000. However the 3DP comparison recognized some decoys such as 1A8E (serum transferrin) as CDK2-like due to an overall similar shape, although this decoy had a clearly different fold as correctly analyzed by each of the three alignment tools (Fig. 10e, left). The three alignment methods correctly assigned a high score to all CDK2 pairs, but also returned a high core with part of the CDK2-decoys which were not recognized as CDK2 like by 3DP. For example decoy 4W9X, which is a non-CDK2 kinase, clearly showed a partly homologous fold to CDK2 leading to a high alignment score, but also showed substantial differences with the presence of a central helix and an extended terminal absent from CDK2 which resulted in a relatively high 3DP-distance to CDK2s such as 3PY1 (Fig. 10e, right).
The 3D-structure of biomolecules in the PDB were encoded in the 136-dimensional 3D atom-pair fingerprint 3DP counting the number of atom pairs at increasing distance intervals for all atoms, positively charged, negatively charged, and hydrophobic atoms. The 3DP fingerprint perceives the spatial distribution of shape, hydrophobicity and charges in molecular objects across a very broad size range. 3DP nearest neighbors are shown in various examples to be closely related shape and fold analogs.
Availability of supporting data
The data sets supporting the results of this article are included within the article and its additional file.
Rose PW, Bi C, Bluhm WF, Christie CH, Dimitropoulos D, Dutta S, et al. The RCSB Protein Data Bank: new resources for research and education. Nucleic Acids Res. 2013;41:D475–82.
Gutmanas A, Alhroub Y, Battle GM, Berrisford JM, Bochet E, Conroy MJ, et al. PDBe: Protein Data Bank in Europe. Nucleic Acids Res. 2014;42:D285–91.
Rose PW, Prlic A, Bi C, Bluhm WF, Christie CH, Dutta S, et al. The RCSB Protein Data Bank: views of structural biology for basic and applied research and education. Nucleic Acids Res. 2015;43:D345–56.
Touw WG, Baakman C, Black J, Te Beek TA, Krieger E, Joosten RP, et al. A series of PDB-related databanks for everyday needs. Nucleic Acids Res. 2015;43:D364–8.
Andreeva A, Howorth D, Chandonia JM, Brenner SE, Hubbard TJ, Chothia C, et al. Data growth and its impact on the SCOP database: new developments. Nucleic Acids Res. 2008;36:D419–25.
Andreeva A, Howorth D, Chothia C, Kulesha E, Murzin AG. SCOP2 prototype: a new approach to protein structure mining. Nucleic Acids Res. 2014;42:D310–4.
Orengo CA, Michie AD, Jones S, Jones DT, Swindells MB, Thornton JM. CATH–a hierarchic classification of protein domain structures. Structure (London, England : 1993). 1997;5:1093–108.
Cuff A, Redfern OC, Greene L, Sillitoe I, Lewis T, Dibley M, et al. The CATH hierarchy revisited-structural divergence in domain superfamilies and the continuity of fold space. Structure (London, England : 1993). 2009;17:1051–62.
Holm L, Rosenstrom P. Dali server: conservation mapping in 3D. Nucleic Acids Res. 2010;38:W545–9.
Shindyalov IN, Bourne PE. Protein structure alignment by incremental combinatorial extension (CE) of the optimal path. Protein Eng. 1998;11:739–47.
Prlic A, Bliven S, Rose PW, Bluhm WF, Bizon C, Godzik A, et al. Pre-calculated protein structure alignments at the RCSB PDB website. Bioinformatics (Oxford, England). 2010;26:2983–5.
Madej T, Lanczycki CJ, Zhang D, Thiessen PA, Geer RC, Marchler-Bauer A, et al. MMDB and VAST+: tracking structural similarities between macromolecular complexes. Nucleic Acids Res. 2014;42:D297–303.
Sael L, Li B, La D, Fang Y, Ramani K, Rustamov R, et al. Fast protein tertiary structure retrieval based on global surface shape similarity. Proteins. 2008;72:1259–73.
La D, Esquivel-Rodriguez J, Venkatraman V, Li B, Sael L, Ueng S, et al. 3D-SURFER: software for high-throughput protein surface comparison and analysis. Bioinformatics (Oxford, England). 2009;25:2843–4.
Awale M, van Deursen R, Reymond JL. MQN-mapplet: visualization of chemical space with interactive maps of DrugBank, ChEMBL, PubChem, GDB-11, and GDB-13. J Chem Inf Model. 2013;53:509–18.
Schwartz J, Awale M, Reymond JL. SMIfp (SMILES fingerprint) chemical space for virtual screening and visualization of large databases of organic molecules. J Chem Inf Model. 2013;53:1979–89.
Reymond JL. The chemical space project. Acc Chem Res. 2015;48:722–30.
Awale M, Reymond JL. Similarity mapplet: interactive visualization of the directory of useful decoys and chembl in high dimensional chemical spaces. J Chem Inf Model. 2015, Epub ahead of print.
Hopkins AL, Keseru GM, Leeson PD, Rees DC, Reynolds CH. The role of ligand efficiency metrics in drug discovery. Nat Rev Drug Discov. 2014;13:105–21.
Awale M, Jin X, Reymond JL. Stereoselective virtual screening of the ZINC database using atom pair 3D-fingerprints. J Cheminform. 2015, 7: doi:10.1186/s13321-13014-10051-13325.
Wirth M, Volkamer A, Zoete V, Rippmann F, Michielin O, Rarey M, et al. Protein pocket and ligand shape comparison and its application in virtual screening. J Comput Aided Mol Des. 2013;27:511–24.
Chen K, Kurgan L. Investigation of atomic level patterns in protein--small ligand interactions. PLoS One. 2009;4, e4473.
Hammes GG. Multiple conformational changes in enzyme catalysis. Biochemistry. 2002;41:8221–8.
James LC, Tawfik DS. Conformational diversity and protein evolution–a 60-year-old hypothesis revisited. Trends Biochem Sci. 2003;28:361–8.
Flores TP, Orengo CA, Moss DS, Thornton JM. Comparison of conformational characteristics in structurally similar protein pairs. Protein Sci. 1993;2:1811–26.
Wood TC, Pearson WR. Evolution of protein sequences and structures. J Mol Biol. 1999;291:977–95.
Orengo CA, Thornton JM. Protein families and their evolution-a structural perspective. Annu Rev Biochem. 2005;74:867–900.
Lee D, Redfern O, Orengo C. Predicting protein function from sequence and structure. Nat Rev Mol Cell Biol. 2007;8:995–1005.
Carhart RE, Smith DH, Venkataraghavan R. Atom pairs as molecular features in structure-activity studies: definition and applications. J Chem Inf Comput Sci. 1985;25:64–73.
Sheridan RP, Miller MD, Underwood DJ, Kearsley SK. Chemical similarity using geometric atom pair descriptors. J Chem Inf Comput Sci. 1996;36:128–36.
Schneider G, Neidhart W, Giller T, Schmid G. “Scaffold-hopping” by topological pharmacophore search: a contribution to virtual screening. Angew Chem Int Ed Engl. 1999;38:2894–6.
Gaulton A, Bellis LJ, Bento AP, Chambers J, Davies M, Hersey A, et al. ChEMBL: a large-scale bioactivity database for drug discovery. Nucleic Acids Res. 2012;40:D1100–7.
Irwin JJ, Sterling T, Mysinger MM, Bolstad ES, Coleman RG. ZINC: a free tool to discover chemistry for biology. J Chem Inf Model. 2012;52:1757–68.
Awale M, Reymond JL. Atom pair 2D-fingerprints perceive 3D-molecular shape and pharmacophores for very fast virtual screening of ZINC and GDB-17. J Chem Inf Model. 2014;54:1892–7.
Krissinel E, Henrick K. Inference of macromolecular assemblies from crystalline state. J Mol Biol. 2007;372:774–97.
Sutcliffe MJ, Haneef I, Carney D, Blundell TL. Knowledge based modelling of homologous proteins, part I: three-dimensional frameworks derived from the simultaneous superposition of multiple structures. Protein Eng. 1987;1:377–84.
Mizuguchi K, Go N. Seeking significance in three-dimensional protein structure comparisons. Curr Opin Struct Biol. 1995;5:377–82.
Case DA, Cheatham 3rd TE, Darden T, Gohlke H, Luo R, Merz Jr KM, et al. The Amber biomolecular simulation programs. J Comput Chem. 2005;26:1668–88.
Gerstein M, Krebs W. A database of macromolecular motions. Nucleic Acids Res. 1998;26:4280–90.
Flores S, Echols N, Milburn D, Hespenheide B, Keating K, Lu J, et al. The database of macromolecular motions: new features added at the decade mark. Nucleic Acids Res. 2006;34:D296–301.
Hsiao CD, Sun YJ, Rose J, Wang BC. The crystal structure of glutamine-binding protein from Escherichia coli. J Mol Biol. 1996;262:225–42.
Sun YJ, Rose J, Wang BC, Hsiao CD. The structure of glutamine-binding protein complexed with glutamine at 1.94 A resolution: comparisons with other amino acid binding proteins. J Mol Biol. 1998;278:219–29.
Bowers KJ, Chow E, Huageng X, Dror RO, Eastwood MP, Gregersen BA, et al. Scalable algorithms for molecular dynamics simulations on commodity clusters. In: SC 2006 Conference, Proceedings of the ACM/IEEE; 11–17 Nov. 2006. 2006. p. 43.
Johnson LN, Noble ME, Owen DJ. Active and inactive protein kinases: structural basis for regulation. Cell. 1996;85:149–58.
Kontopidis G, McInnes C, Pandalaneni SR, McNae I, Gibson D, Mezna M, et al. Differential binding of inhibitors to active and inactive CDK2 provides insights for drug design. Chem Biol. 2006;13:201–11.
Schrodinger LLC. The PyMOL molecular graphics system, version 1.3r1. 2010.
Garrity GM, Lilburn TG. Mapping taxonomic space: an overview of the road map to the second edition of Bergey’s Manual of Systematic Bacteriology. WFCC Newsl. 2002;35:5–15.
Medina-Franco JL, Maggiora GM, Giulianotti MA, Pinilla C, Houghten RA. A similarity-based data-fusion approach to the visual characterization and comparison of compound databases. Chem Biol Drug Des. 2007;70:393–412.
Sauer WH, Schwarz MK. Molecular shape diversity of combinatorial libraries: a prerequisite for broad bioactivity. J Chem Inf Comput Sci. 2003;43:987–1003.
Ertl P, Patiny L, Sander T, Rufener C, Zasso M. Wikipedia Chemical Structure Explorer: substructure and similarity searching of molecules from Wikipedia. J Cheminform. 2015;7:10.
Martinez A, Ortiz Mellet C, Garcia Fernandez JM. Cyclodextrin-based multivalent glycodisplays: covalent and supramolecular conjugates to assess carbohydrate-protein interactions. Chem Soc Rev. 2013;42:4746–73.
Craik DJ, Adams DJ. Chemical modification of conotoxins to improve stability and activity. ACS Chem Biol. 2007;2:457–68.
Hu SH, Gehrmann J, Guddat LW, Alewood PF, Craik DJ, Martin JL. The 1.1 A crystal structure of the neuronal acetylcholine receptor antagonist, alpha-conotoxin PnIA from Conus pennaceus. Structure (London, England : 1993). 1996;4:417–23.
Farr-Jones S, Miljanich GP, Nadasdi L, Ramachandran J, Basus VJ. Solution structure of omega-conotoxin MVIIC, a high affinity ligand of P-type calcium channels, using 1H NMR spectroscopy and complete relaxation matrix analysis. J Mol Biol. 1995;248:106–24.
McGinnis S, Madden TL. BLAST: at the core of a powerful and diverse set of sequence analysis tools. Nucleic Acids Res. 2004;32:W20–5.
Brylinski M, Skolnick J. Comparison of structure-based and threading-based approaches to protein functional annotation. Proteins. 2010;78:118–34.
Garriga D, Querol-Audi J, Abaitua F, Saugar I, Pous J, Verdaguer N, et al. The 2.6-Angstrom structure of infectious bursal disease virus-derived T = 1 particles reveals new stabilizing elements of the virus capsid. J Virol. 2006;80:6895–905.
Irigoyen N, Garriga D, Navarro A, Verdaguer N, Rodriguez JF, Caston JR. Autoproteolytic activity derived from the infectious bursal disease virus capsid protein. J Biol Chem. 2009;284:8064–72.
Coulibaly F, Chevalier C, Delmas B, Rey FA. Crystal structure of an Aquabirnavirus particle: insights into antigenic diversity and virulence determinism. J Virol. 2010;84:1792–9.
Zhang Y, Skolnick J. TM-align: a protein structure alignment algorithm based on the TM-score. Nucleic Acids Res. 2005;33:2302–9.
Pandit SB, Skolnick J. Fr-TM-align: a new protein structural alignment method based on fragment alignments and the TM-score. BMC bioinformatics. 2008;9:531.
Yang Y, Zhan J, Zhao H, Zhou Y. A new size-independent score for pairwise protein structure alignment and its application to structure classification and nucleic-acid binding prediction. Proteins. 2012;80:2080–8.
Menke M, Berger B, Cowen L. Matt: local flexibility aids protein multiple structure alignment. PLoS Comput Biol. 2008;4, e10.
This work was supported financially by the University of Berne, the Swiss National Science Foundation, and the NCCR TransCure.
The authors declare no competing financial interests.
XJ designed and programmed 3DP, carried out the PDB study, designed and programmed the Protein Explorer website, and wrote the paper. MA designed and programmed 3DP and wrote the paper. MZ, DK and LP designed and programmed the Protein Explorer website. JLR designed and supervised the study and wrote the paper. All authors have read and approved the final manuscript.
Supplement 1. PDB IDs of 91,223 X-ray structures downloaded from the PDB in September 2014. Supplement 2. PDB IDs of X-ray structures of cyclin-dependent kinase 2 (CDK2). Supplement 3. PDB IDs of X-ray structures of triose phosphate isomerase (TIM) dimer. Supplement 4. IDs of CATH superfamilies and PDB IDs in each CATH superfamily. (DOCX 297 kb)
About this article
Cite this article
Jin, X., Awale, M., Zasso, M. et al. PDB-Explorer: a web-based interactive map of the protein data bank in shape space. BMC Bioinformatics 16, 339 (2015). https://doi.org/10.1186/s12859-015-0776-9
- Protein data bank
- Chemical space
- Molecular shape