PDB-Explorer: a web-based interactive map of the protein data bank in shape space

Background 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. Methods 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. Results The 3DP property space groups proteins and protein assemblies according to their 3D-shape similarity, yet shows exquisite ability to distinguish between closely related structures. An interactive website called PDB-Explorer is presented featuring a color-coded interactive map of PDB in 3DP-space. Each pixel of the map contains one or more PDB-entries which are directly visualized as ribbon diagrams when the pixel is selected. The PDB-Explorer website allows performing 3DP-nearest neighbor searches of any PDB-entry or of any structure uploaded as protein-type PDB file. All functionalities on the website are implemented in JavaScript in a platform-independent manner and draw data from a server that is updated daily with the latest PDB additions, ensuring complete and up-to-date coverage. The essentially instantaneous 3DP-similarity search with the PDB-Explorer provides results comparable to those of much slower 3D-alignment algorithms, and automatically clusters proteins from the same superfamilies in tight groups. Conclusion 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.

Similarity Mapplet: Interactive Visualization of the Directory of Useful Decoys and ChEMBL in High Dimensional Chemical Spaces

An Internet portal accessible at www.gdb.unibe.ch has been set up to automatically generate color-coded similarity maps of the ChEMBL database in relation to up to two sets of active compounds taken from the enhanced Directory of Useful Decoys (eDUD), a random set of molecules, or up to two sets of user-defined reference molecules. These maps visualize the relationships between the selected compounds and ChEMBL in six different high dimensional chemical spaces, namely MQN (42-D molecular quantum numbers), SMIfp (34-D SMILES fingerprint), APfp (20-D shape fingerprint), Xfp (55-D pharmacophore fingerprint), Sfp (1024-bit substructure fingerprint), and ECfp4 (1024-bit extended connectivity fingerprint). The maps are supplied in form of Java based desktop applications called “similarity mapplets” allowing interactive content browsing and linked to a “Multifingerprint Browser for ChEMBL” (also accessible directly at www.gdb.unibe.ch) to perform nearest neighbor searches. One can obtain six similarity mapplets of ChEMBL relative to random reference compounds, 606 similarity mapplets relative to single eDUD active sets, 30 300 similarity mapplets relative to pairs of eDUD active sets, and any number of similarity mapplets relative to user-defined reference sets to help visualize the structural diversity of compound series in drug optimization projects and their relationship to other known bioactive compounds.

The Chemical Space Project

One of the simplest questions that can be asked about molecular diversity is how many organic molecules are possible in total? To answer this question, my research group has computationally enumerated all possible organic molecules up to a certain size to gain an unbiased insight into the entire chemical space. Our latest database, GDB-17, contains 166.4 billion molecules of up to 17 atoms of C, N, O, S, and halogens, by far the largest small molecule database reported to date. Molecules allowed by valency rules but unstable or nonsynthesizable due to strained topologies or reactive functional groups were not considered, which reduced the enumeration by at least 10 orders of magnitude and was essential to arrive at a manageable database size. Despite these restrictions, GDB-17 is highly relevant with respect to known molecules. Beyond enumeration, understanding and exploiting GDBs (generated databases) led us to develop methods for virtual screening and visualization of very large databases in the form of a “periodic system of molecules” comprising six different fingerprint spaces, with web-browsers for nearest neighbor searches, and the MQN- and SMIfp-Mapplet application for exploring color-coded principal component maps of GDB and other large databases. Proof-of-concept applications of GDB for drug discovery were realized by combining virtual screening with chemical synthesis and activity testing for neurotransmitter receptor and transporter ligands. One surprising lesson from using GDB for drug analog searches is the incredible depth of chemical space, that is, the fact that millions of very close analogs of any molecule can be readily identified by nearest-neighbor searches in the MQN-space of the various GDBs. The chemical space project has opened an unprecedented door on chemical diversity. Ongoing and yet unmet challenges concern enumerating molecules beyond 17 atoms and synthesizing GDB molecules with innovative scaffolds and pharmacophores.

Stereoselective virtual screening of the ZINC database using atom pair 3D-fingerprints

Background Tools to explore large compound databases in search for analogs of query molecules provide a strategically important support in drug discovery to help identify available analogs of any given reference or hit compound by ligand based virtual screening (LBVS). We recently showed that large databases can be formatted for very fast searching with various 2D-fingerprints using the city-block distance as similarity measure, in particular a 2D-atom pair fingerprint (APfp) and the related category extended atom pair fingerprint (Xfp) which efficiently encode molecular shape and pharmacophores, but do not perceive stereochemistry. Here we investigated related 3D-atom pair fingerprints to enable rapid stereoselective searches in the ZINC database (23.2 million 3D structures). Results Molecular fingerprints counting atom pairs at increasing through-space distance intervals were designed using either all atoms (16-bit 3DAPfp) or different atom categories (80-bit 3DXfp). These 3D-fingerprints retrieved molecular shape and pharmacophore analogs (defined by OpenEye ROCS scoring functions) of 110,000 compounds from the Cambridge Structural Database with equal or better accuracy than the 2D-fingerprints APfp and Xfp, and showed comparable performance in recovering actives from decoys in the DUD database. LBVS by 3DXfp or 3DAPfp similarity was stereoselective and gave very different analogs when starting from different diastereomers of the same chiral drug. Results were also different from LBVS with the parent 2D-fingerprints Xfp or APfp. 3D- and 2D-fingerprints also gave very different results in LBVS of folded molecules where through-space distances between atom pairs are much shorter than topological distances. Conclusions 3DAPfp and 3DXfp are suitable for stereoselective searches for shape and pharmacophore analogs of query molecules in large databases. Web-browsers for searching ZINC by 3DAPfp and 3DXfp similarity are accessible at www.gdb.unibe.ch webcite and should provide useful assistance to drug discovery projects.

TimX, a novel player in protein import across the inner mitochondrial membrane of T. brucei

Mitochondrial protein import is an essential function of the unique mitochondrion in T. brucei as roughly 1000 different nuclear encoded proteins need to be correctly localized to their mitochondrial subcompartment. For this reason the responsible import machinery is expected to be similarly complex as in other Eukaryotes. This was recently demonstrated for the translocation machinery in the outer mitochondrial membrane. In contrast, the composition of the inner membrane import machinery and the exact molecular pathway(s) taken by various substrates are still ill-defined. To elucidate this further, we performed a pulldown analysis of epitope tagged TbTim17 in combination with quantitative mass spectrometry. By this we identified novel components of the mitochondrial import machinery in trypanosomes. One of these, TimX, is an essential mitochondrial membrane protein of 42 kDa that is unique to kinetoplastids. This protein migrates on Blue Native PAGE in a high molecular weight complex similar to TbTim17. Ablation of either of the two proteins leads to a destabilization of the complex containing the other protein. Furthermore, its involvement in protein import could be demonstrated by in vivo and in vitro protein import assays. This corroborates that TimX together with TbTim17 forms a protein import complex in the inner mitochondrial membrane. As TbTim17 the TimX protein was subjected to pulldown analysis in combination with quantitative mass spectrometry. The overlap of candidates defined by these two sets of IPs likely defines further components of the inner membrane translocase which are presently being analyzed. In summary our study on novel components of the trypanosome mitochondrial protein import system gives us fascinating new insights into evolution of the mitochondrion.

The TIM translocase in T. brucei contains Rhomboid-like proteins that are essential for mitochondrial protein import

The parasitic protozoon Trypanosoma brucei is one of the earliest branching eukaryotes that have mitochondria capable of oxidative phosphorylation. Their protein import systems are of similar complexity yet different composition than those in other eukaryotes. To elucidate the composition of the trypanosomal translocase of the inner mitochondrial membrane (TIM) we performed CoIPs of epitope-tagged TbTim17 and two other candidates in combination with SILAC-based quantitative mass spectrometry. This led to the identification of ten candidates for core TIM subunits. Eight of them were present in the previously determined inner membrane proteome and four show homology to small Tim chaperones. Three candidates, a trypanosomatid-specific 42 kDa protein (Tim42) and two putative orthologues of inactive rhomboid proteases were analyzed further. All three proteins are essential in both life cycle stages and their ablation results in a strong protein import defect in vivo and in vitro. Blue native PAGE revealed their presence in a high molecular weight complex. Unlike anticipated, trypanosomes have a highly complex TIM translocase that has extensively been redesigned. None of the three novel TIM subunits has ever been associated with mitochondrial protein import. Two of them belong to the rhomboid protease family, a member of which recently has been implicated in the ERAD translocation system. This suggests an exciting analogy between protein translocases of mitochondria and the ER.

The TIM translocase in T. brucei contains Rhomboid-like proteins that are essential for mitochondrial protein import

The parasitic protozoon Trypanosoma brucei is often considered as one of the earliest branching eukaryotes that have mitochondria capable of oxidative phosphorylation. Its protein import systems are therefore of great interest. Recently, it was shown that the outer mitochondrial membrane protein translocase is of similar complexity yet different composition than in other eukaryotes (1). In the inner membrane however, only a single orthologue of the pore forming Tim17/22/23 protein family was identified and termed TbTim17. Based on this finding it has been suggested that, instead of separate TIM22 and TIM23 complexes as in other eukaryotes, trypanosomes may have a single multifunctional translocase of the inner mitochondrial membrane (TIM) of reduced complexity. To elucidate the composition of the trypanosomal TIM complex we performed co-immunoprecipitations (CoIP) of epitope-tagged TbTim17 in combination with SILAC-based quantitative mass spectrometry. This led to the identification of 22 highly enriched TbTim17-interacting proteins. We tagged two of the top-scoring proteins for reciprocal CoIP analyses and recovered a set of ten proteins that are highly enriched in all three CoIPs. These proteins are excellent candidates for core subunits of the trypanosomal TIM complex. Eight of them were present in the previously determined inner membrane proteome and four show homology to small Tim chaperones. Three candidates, a novel trypanosomatid-specific 42 kDa protein, termed Tim42, and two putative orthologues of probably inactive rhomboid proteases were chosen for further analysis. All three proteins are essential in both life cycle stages and in a cell line that can grow in the absence of mitochondrial DNA. Additionally, their ablation by RNAi results in a strong protein import defect both in vivo and in vitro. Blue native PAGE reveals that Tim42, like TbTim17 is present in a high molecular weight complex. Moreover, ablation of either Tim42 or TbTim17 leads to a destabilization of the complex containing the other protein, suggesting a tight interaction of the two proteins. In summary our study shows that unlike anticipated trypanosomes have a highly complex TIM translocase that has extensively been redesigned. We have characterized three novel TIM subunits that have never been associated with mitochondrial protein import before. Two of them belong to the rhomboid protease family, a member of which recently has been implicated in the ERAD translocation system. Our study provides insight into mitochondrial evolution over large phylogenetic distances and suggests an exciting analogy between protein translocation systems of mitochondria and the ER.

An intergenic non-coding RNA targets and regulates the ribosome in H. volcanii

As translation is the final step in gene expression it is particularly important to understand the processes involved in translation regulation. It was shown in the last years that a class of RNA, the non-protein-coding RNAs (ncRNAs), is involved in regulation of gene expression via various mechanisms (e.g. gene silencing by microRNAs). Almost all of these ncRNA discovered so far target the mRNA in order to modulate protein biosynthesis, this is rather unexpected considering the crucial role of the ribosome during gene expression. However, recent data from our laboratory showed that there is a new class of ncRNAs, which target the ribosome itself [Gebetsberger et al., 2012/ Pircher et al, 2014]. These so called ribosome-associated ncRNAs (rancRNAs) have an impact on translation regulation, mainly by interfering / modulating the rate of protein biosynthesis. The main goal of this project is to identify and describe novel potential regulatory rancRNAs in H. volcanii with the focus on intergenic candidates. Northern blot analyses already revealed interactions with the ribosome and showed differential expression of rancRNAs during different growth phases or under specific stress conditions. To investigate the biological relevance of these rancRNAs, knock-outs were generated in H. volcanii which were used for phenotypic characterization studies. The rancRNA s194 showed association with the 50S ribosomal subunit in vitro and in vivo and was capable of inhibiting peptide bond formation and seems to inhibit translation in vitro. These preliminary data for the rancRNA s194 make it an interesting candidate for further functional studies to identify the molecular mechanisms by which rancRNAs can modulate protein biosynthesis. Characterization of further rancRNA candidates are also underway.

Mechanistic insights into tRNA translocation on the ribosome

The ribosome is a highly conserved cellular complex and constitutes the center of protein biosynthesis. As the ribosome consists to about 2/3 of ribosomal RNA (rRNA), the rRNA is involved in most steps of translation. In order to investigate the role of some defined rRNA residues in different aspects of translation we use the atomic mutagenesis approach. This method allows the site-specific incorporation of unnatural nucleosides into the rRNA in the context of the complete 70S from Thermus aquaticus, and thereby exceeds the possibilities of conventional mutagenesis. We first studied ribosome-stimulated EF-G GTP hydrolysis. Here, we could show that the non-bridging phosphate oxygen of A2662, which is part of the Sarcin-Ricin-Loop, is required for EF-G GTPase activation by the ribosome. EF-G GTPase is a crucial step for tRNA translocation from the A- to the P-site, and from the P- to the E-site, respectively. We furthermore used the atomic mutagenesis approach to more precisely characterize the 23S rRNA functional groups involved in E-site tRNA binding. While the ribosomal A- and P-sites have been functionally well characterized in the past, the contribution of the E-site to protein biosynthesis is still poorly understood in molecular terms. Our data disclose the importance of the highly conserved E-site base pair G2421-C2395 for effective translation. Ribosomes with a disrupted G2421-C2395 base pair are defective in tRNA binding to the E-site. This results in an impaired translation of genuine mRNAs, while homo-polymeric templates are not affected. Cumulatively our data emphasize the importance of E-site tRNA occupancy and in particular the intactness of the 23S rRNA base pair G2421-C2395 for productive protein biosynthesis.

Identification and characterisation of ncRNAs in dormant bacteria

Recent years have led to increasing interest and appreciation of the possible importance of single cell heterogeneity in various biological processes. One of the examples of phenotypic heterogeneity in bacterial populations is antibiotic tolerant persister cells. Such an antibiotic tolerance phenotype is of considerable clinical relevance since dormant bacteria can re-establish infections rapidly after the antibiotic treatment has been terminated. Up to now mechanisms for establishing the persistence phenomenon in bacteria have remained largely enigmatic. Persisters are cells considered to be in a dormant state with down regulated gene expression. Only recently small regulatory RNAs (sRNAs) have been appreciated as important regulators of gene expression in response to environmental stimuli and several theoretical studies have suggested a possible involvement of sRNAs in the mechanisms of regulated heterogeneity in bacteria. We have experimentally addressed this potential link between sRNAs and persistence/dormancy in E. coli as an example of heterogeneity. Beside classical sRNAs we are focusing also on sRNAs directly associating with and possibly regulating the ribosome, the central enzyme of gene expression. The persister and dormant cell specific sRNA profile is studied by the comparative analysis of sRNA profile changes of the whole bacterial population after antibiotic killing. From RNA-Seq data ~ 25 000 potentially stable RNA fragments were identified and initial analysis predicted ~300 of them to be dormant/persister cell specific. After further evaluation the most prominent dormant/persister cell specific sRNAs are functionally characterized and their potential role in the persistence/dormancy will be evaluated by applying genetic, molecular and biochemical tools. The potential results of this project will provide a better understanding on the molecular mechanism of bacterial persistence/dormancy and on the role of ribosome-bound sRNA molecules in fine-tuning gene expression.

Identification and characterisation of ncRNAs in dormant bacteria

Recent years have led to increasing interest and appreciation of the possible importance of single cell heterogeneity in various biological processes. One of the examples of phenotypic heterogeneity in bacterial populations is antibiotic tolerant persister cells. Such an antibiotic tolerance phenotype is of considerable clinical relevance since dormant bacteria can re-establish infections rapidly after the antibiotic treatment has been terminated. Up to now mechanisms for establishing the persistence phenomenon in bacteria have remained largely enigmatic. Persisters are cells considered to be in a dormant state with down regulated gene expression. Only recently small regulatory RNAs (sRNAs) have been appreciated as important regulators of gene expression in response to environmental stimuli and several theoretical studies have suggested a possible involvement of sRNAs in the mechanisms of regulated heterogeneity in bacteria. We have experimentally addressed this potential link between sRNAs and persistence/dormancy in E. coli as an example of heterogeneity. Beside classical sRNAs we are focusing also on sRNAs directly associating with and possibly regulating the ribosome, the central enzyme of gene expression. The persister and dormant cell specific sRNA profile is studied by the comparative analysis of sRNA profile changes of the whole bacterial population after antibiotic killing. From RNA-Seq data ~ 25 000 potentially stable RNA fragments were identified and initial analysis predicted ~300 of them to be dormant/persister cell specific. After further evaluation the most prominent dormant/persister cell specific sRNAs are functionally characterized and their potential role in the persistence/dormancy will be evaluated by applying genetic, molecular and biochemical tools. The potential results of this project will provide a better understanding on the molecular mechanism of bacterial persistence/dormancy and on the role of ribosome-bound sRNA molecules in fine-tuning gene expression.

Identification and characterisation of ncRNAs in dormant bacteria

Recent years have led to increasing interest and appreciation of the possible importance of single cell heterogeneity in various biological processes. One of the examples of phenotypic heterogeneity in bacterial populations is antibiotic tolerant persister cells. Such an antibiotic tolerance phenotype is of considerable clinical relevance since dormant bacteria can re-establish infections rapidly after the antibiotic treatment has been terminated. Up to now mechanisms for establishing the persistence phenomenon in bacteria have remained largely enigmatic. Persisters are cells considered to be in a dormant state with down regulated gene expression. Only recently small regulatory RNAs (sRNAs) have been appreciated as important regulators of gene expression in response to environmental stimuli and several theoretical studies have suggested a possible involvement of sRNAs in the mechanisms of regulated heterogeneity in bacteria. We have experimentally addressed this potential link between sRNAs and persistence/dormancy in E. coli as an example of heterogeneity. Beside classical sRNAs we are focusing also on sRNAs directly associating with and possibly regulating the ribosome, the central enzyme of gene expression. The persister and dormant cell specific sRNA profile is studied by the comparative analysis of sRNA profile changes of the whole bacterial population after antibiotic killing. From RNA-Seq data ~ 25 000 potentially stable RNA fragments were identified and initial analysis predicted ~300 of them to be dormant/persister cell specific. After further evaluation the most prominent dormant/persister cell specific sRNAs are functionally characterized and their potential role in the persistence/dormancy will be evaluated by applying genetic, molecular and biochemical tools. The potential results of this project will provide a better understanding on the molecular mechanism of bacterial persistence/dormancy and on the role of ribosome-bound sRNA molecules in fine-tuning gene expression.

Effect of lipid bilayer properties on the photocycle of green proteorhodopsin

The significance of specific lipids for proton pumping by the bacterial rhodopsin proteorhodopsin (pR) was studied. To this end, it was examined whether pR preferentially binds certain lipids and whether molecular properties of the lipid environment affect the photocycle. pR's photocycle was followed by microsecond flash-photolysis in the visible spectral range. It was fastest in phosphatidylcholine liposomes (soy bean lipid), intermediate in 3-[(3-cholamidopropyl) dimethylammonio] propanesulfonate (CHAPS): 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC) bicelles and in Triton X-100, and slowest when pR was solubilized in CHAPS. In bicelles with different lipid compositions, the nature of the head groups, the unsaturation level and the fatty acid chain length had small effects on the photocycle. The specific affinity of pR for lipids of the expression host Escherichia coli was investigated by an optimized method of lipid isolation from purified membrane protein using two different concentrations of the detergent N-dodecyl-β-d-maltoside (DDM). We found that 11 lipids were copurified per pR molecule at 0.1% DDM, whereas essentially all lipids were stripped off from pR by 1% DDM. The relative amounts of copurified phosphatidylethanolamine, phosphatidylglycerol, and cardiolipin did not correlate with the molar percentages normally present in E. coli cells. The results indicate a predominance of phosphatidylethanolamine species in the lipid annulus around recombinant pR that are less polar than the dominant species in the cell membrane of the expression host E. coli.