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.

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.

What can a unicellular parasite tell us about mitochondrial biogenesis?

In the unicellular parasite Trypanosoma brucei, as in other eukaryotes, more than 95% of all mitochondrial proteins are imported from the cytosol. The recently characterized multisubunit ATOM complex, the functional analogue of the TOM complex of yeast, mediates import of essentially all proteins across the outer mitochondrial membrane in T. brucei. Moreover, an additional protein termed pATOM36, which is loosely associated with the ATOM complex, has been implicated in the import of only a subset of mitochondrial proteins. Here we have investigated more precisely which role pATOM36 plays in mitochondrial protein import. RNAi mediated ablation of pATOM36 specifically depletes a subset of outer mitochondrial membrane proteins including ATOM complex subunits and as a consequence results in the collapse of the ATOM complex as shown by Blue native PAGE. In addition, a SILAC-based global proteomic analysis of uninduced and induced pATOM36 RNAi cells together with in vitro import experiments suggest that pATOM36 might be a novel protein import factor acting on a subset of alpha-helically anchored mitochondrial outer membrane proteins. Identification of pATOM36 interaction partners by co-immunoprecipitation together with immunofluorescence analysis shows that unexpectedly a fraction of the protein is associated with the tripartite attachment complex (TAC). This complex is essential for proper inheritance of the mitochondrial DNA in T. brucei. It forms a physical connection between the single unit mitochondrial DNA and the basal body of the flagellum that is stable throughout the cell cycle. Thus, pATOM36 simultaneously mediates ATOM assembly, and thus protein import, as well as mitochondrial DNA inheritance since it is an essential component of the TAC.

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.

Probing the Structure, Pseudorotation, and Radial Vibrations of Cyclopentane by Femtosecond Rotational Raman Coherence Spectroscopy

Femtosecond time-resolved Raman rotational coherence spectroscopy (RCS) is employed to determine accurate rotational, vibration–rotation coupling constants, and centrifugal distortion constants of cyclopentane (C⁵H¹⁰). Its lowest-frequency vibration is a pseudorotating ring deformation that interconverts 10 permutationally distinct but energetically degenerate “twist” minima interspersed by 10 “bent” conformers. While the individual twist and bent structures are polar asymmetric tops, the pseudorotation is fast on the time scale of external rotation, rendering cyclopentane a fluxionally nonpolar symmetric top molecule. The pseudorotational level pattern corresponds to a one-dimensional internal rotor with a pseudorotation constant Bps ≈ 2.8 cm⁻¹. The pseudorotational levels are significantly populated up to l = ± 13 at 298 K; <10% of the molecules are in the l = 0 level. The next-higher vibration is the “radial” ν²³ ring deformation mode at 273 cm⁻¹, which is far above the pseudorotational fundamental. Femtosecond Raman RCS measurements were performed in a gas cell at T = 293 K and in a pulsed supersonic jet at T ≈ 90 K. The jet cooling reduces the pseudorotational distribution to l < ±8 and eliminates the population of ν²³, allowing one to determine the rotational constant as A0 = B0 = 6484.930(11) MHz. This value is ∼300 times more precise than the previous value. The fit of the RCS transients reveals that the rotation–pseudorotation coupling constant αe,psB = −0.00070(1) MHz is diminutive, implying that excitation of the pseudorotation has virtually no effect on the B0 rotational constant of cyclopentane. The smallness of αe,psB can be realized when comparing to the vibration–rotation coupling constant of the ν²³ vibration, αe,23B = −9.547(1) MHz, which is about 10⁴ times larger.

Rotational constants and structure of para-difluorobenzene determined by femtosecond Raman coherence spectroscopy: A new transient type

Femtosecond Raman rotational coherence spectroscopy (RCS) detected by degenerate four-wave mixing is a background-free method that allows to determine accurate gas-phase rotational constants of non-polar molecules. Raman RCS has so far mostly been applied to the regular coherence patterns of symmetric-top molecules, while its application to nonpolar asymmetric tops has been hampered by the large number of RCS transient types, the resulting variability of the RCS patterns, and the 10³–10⁴ times larger computational effort to simulate and fit rotational Raman RCS transients. We present the rotational Raman RCS spectra of the nonpolar asymmetric top 1,4-difluorobenzene (para-difluorobenzene, p-DFB) measured in a pulsed Ar supersonic jet and in a gas cell over delay times up to ~2.5 ns. p-DFB exhibits rotational Raman transitions with ΔJ = 0, 1, 2 and ΔK = 0, 2, leading to the observation of J −, K −, A −, and C–type transients, as well as a novel transient (S–type) that has not been characterized so far. The jet and gas cell RCS measurements were fully analyzed and yield the ground-state (v = 0) rotational constants Aₒ = 5637.68(20) MHz, Bₒ = 1428.23(37) MHz, and Cₒ = 1138.90(48) MHz (1σ uncertainties). Combining the Aₒ, Bₒ, and Cₒ constants with coupled-cluster with single-, double- and perturbatively corrected triple-excitation calculations using large basis sets allows to determine the semi-experimental equilibrium bond lengths rₑ(C₁–C₂) = 1.3849(4) Å, rₑ(C₂–C³) = 1.3917(4) Å, rₑ(C–F) = 1.3422(3) Å, and rₑ(C₂–H₂) = 1.0791(5) Å.

Three-State Single-Molecule Naphthalenediimide Switch: Integration of a Pendant Redox Unit for Conductance Tuning

We studied charge transport through core-substituted naphthalenediimide (NDI) single-molecule junctions using the electrochemical STM-based break-junction technique in combination with DFT calculations. Conductance switching among three well-defined states was demonstrated by electrochemically controlling the redox state of the pendent diimide unit of the molecule in an ionic liquid. The electrical conductances of the dianion and neutral states differ by more than one order of magnitude. The potential-dependence of the charge-transport characteristics of the NDI molecules was confirmed by DFT calculations, which account for electrochemical double-layer effects on the conductance of the NDI junctions. This study suggests that integration of a pendant redox unit with strong coupling to a molecular backbone enables the tuning of charge transport through single-molecule devices by controlling their redox states.

Electrochemical CO2 Reduction - A Critical View on Fundamentals, Materials and Applications

The electrochemical reduction of CO2 has been extensively studied over the past decades. Nevertheless, this topic has been tackled so far only by using a very fundamental approach and mostly by trying to improve kinetics and selectivities toward specific products in half-cell configurations and liquid-based electrolytes. The main drawback of this approach is that, due to the low solubility of CO2 in water, the maximum CO2 reduction current which could be drawn falls in the range of 0.01–0.02 A cm–2. This is at least an order of magnitude lower current density than the requirement to make CO2-electrolysis a technically and economically feasible option for transformation of CO2 into chemical feedstock or fuel thereby closing the CO2 cycle. This work attempts to give a short overview on the status of electrochemical CO2 reduction with respect to challenges at the electrolysis cell as well as at the catalyst level. We will critically discuss possible pathways to increase both operating current density and conversion efficiency in order to close the gap with established energy conversion technologies.

Tandem Mass Spectrometric Analysis of Metallocene Adducts with Short Oligonucleotides

Metallocene dichlorides constitute a remarkable class of antineoplastic agents that are highly effective against several cancer cell lines. They were shown to accumulate in the DNA-rich region, which suggests DNA as the primary target. These compounds exhibit two cyclopentadienyl ligands and two labile halide ligands, resulting in a bent sandwich structure. The cis-dihalide motif is structurally related to the cis-chloro configuration of cisplatin and similar modes of action can thus be assumed. Cisplatin binds to two neighboring guanine nucleobases in DNA and consequently, distorts the double-helix, thereby inhibiting DNA replication and transcription. Platinum is classified as a soft Lewis acid and binds preferentially to the nitrogen atoms within the nucleobases. The metallocene dichlorides investigated in this study comprise the metal centers Ti, V, Nb, Mo, Hf, and W, which are classified as hard or intermediate Lewis acids, and thus, favor binding to the phosphate oxygen. Although several studies reported adduct formation of metallocene dichlorides with nucleic acids, substantial information about the adduct composition, the binding pattern, and the nucleobase selectivity has not been provided yet. ESI-MS analyses gave evidence for the formation of metallocene adducts (M = Ti, V, Mo, and W) with single-stranded DNA homologues at pH 7. No adducts were formed with Nb and Hf at neutral pH, albeit adducts with Nb were observed at a low pH. MS2 data revealed considerable differences of the adduct compositions. The product ion spectra of DNA adducts with hard Lewis acids (Ti, V) gave evidence for the loss of metallocene ligands and only moderate backbone fragmentation was observed. By contrast, adducts with intermediate Lewis acids (Mo, W) retained the hydroxy ligands. Preliminary results are in good agreement with the Pearson concept and DFT calculations. Since the metallodrugs were not lost upon CID, the nucleobase selectivity, stoichiometry, and binding patterns can be elucidated by means of tandem mass spectrometry.

Structural elucidation and antisense properties of a sugar-modified DNA analogue

Antisense oligonucleotides deserve great attention as potential drug candidates for the treatment of genetic disorders. For example, muscle dystrophy can be treated successfully in mice by antisense-induced exon skipping in the pre-mRNA coding for the structural protein dystrophin in muscle cells. For this purpose a sugar- and backbone-modified DNA analogue was designed, in which a tricyclic ring system substitutes the deoxyribose. These chemical modifications stabilize the dimers formed with the targeted RNA relative to native nucleic acid duplexes and increase the biostability of the antisense oligonucleotide. While evading enzymatic degradation constitutes an essential property of antisense oligonucleotides for therapeutic application, it renders the oligonucleotide inaccessible to biochemical sequencing techniques and requires the development of alternative methods based on mass spectrometry. The set of sequences studied includes tcDNA oligonucleotides ranging from 10 to 15 nucleotides in length as well as their hybrid duplexes with DNA and RNA complements. All samples were analyzed on a LTQ Orbitrap XL instrument equipped with a nano-electrospray source. For tandem mass spectrometric experiments collision-induced dissociation was performed, using helium as collision gas. Mass spectrometric sequencing of tcDNA oligomers manifests the applicability of the technique to substrates beyond the scope of enzyme-based methods. Sequencing requires the formation of characteristic backbone fragments, which take the form of a-B- and w-ions in the product ion spectra of tcDNA. These types of product ions are typically associated with unmodified DNA, which suggests a DNA-like fragmentation mechanism in tcDNA. The loss of nucleobases constitutes the second prevalent dissociation pathway observed in tcDNA. Comparison of partially and fully modified oligonucleotides indicates a pronounced impact of the sugar-moiety on the base loss. As this event initiates cleavage of the backbone, the presented results provide new mechanistic insights into the fragmentation of DNA in the gas-phase. The influence of the sugar-moiety on the dissociation extends to tcDNA:DNA and tcDNA:RNA hybrid duplexes, where base loss was found to be much more prominent from sugar-modified oligonucleotides than from their natural complements. Further prominent dissociation channels are strand separation and backbone cleavage of the single strands, as well as the ejection of backbone fragments from the intact duplex. The latter pathway depends noticeably on the base sequence. Moreover, it gives evidence of the high stability of the hybrid dimers, and thus directly reflects the affinity of tcDNA for its target in the cell. As the cellular target of tcDNA is a pre-mRNA, the structure was designed to discriminate RNA from DNA complements, which could be demonstrated by mass spectrometric experiments.

Polar Nature of Biomimetic Fluorapatite/Gelatin Composites: A Comparison of Bipolar Objects and the Polar State of Natural Tissue

The correspondence of the state of alignment of macromolecules in biomimetic materials and natural tissues is demonstrated by investigating a mechanism of electrical polarity formation: An in vitro grown biomimetic FAp/gelatin composite is investigated for its polar properties by second harmonic (SHGM) and scanning pyroelectric microscopy (SPEM). Hexagonal prismatic seed crystals formed in gelatin gels represent a monodomain polar state, due to aligned mineralized gelatin molecules. Later growth stages, showing dumbbell morphologies, develop into a bipolar state because of surface recognition by gelatin functionality: A reversal of the polar alignment of macromolecules, thus, takes place close to that basal plane of the seed. In natural hard tissues (teeth and bone investigated by SPEM) and the biomimetic FAp/gelatin composite, we find a surprising analogy in view of growth-induced states of polarity: The development of polarity in vivo and in vitro can be explained by a Markov-type mechanism of molecular recognition during the attachment of macromolecules.

pATOM36 mediates mitochondrial protein import and mitochondrial DNA inheritance

The multisubunit ATOM complex mediates import of essentially all proteins across the outer mitochondrial membrane in T. brucei. Moreover, an additional protein termed pATOM36, which is loosely associated with the ATOM complex, has been implicated in the import of only a subset of mitochondrial matrix proteins. Here we have investigated more precisely which role pATOM36 plays in mitochondrial protein import. RNAi mediated ablation of pATOM36 specifically depletes a subset of ATOM complex subunits and as a consequence results in the collapse of the ATOM complex as shown by Blue native PAGE. In addition, a SILAC-based global proteomic analysis of uninduced and induced pATOM36 RNAi cells together with in vitro import experiments suggest that pATOM36 might be a novel protein insertase acting on a subset of alpha-helically anchored mitochondrial outer membrane proteins. Identification of pATOM36 interaction partners by co-immunoprecipitation together with immunofluorescence analysis furthermore shows that unexpectedly a fraction of the protein is associated with the tripartite attachment complex (TAC). This complex is essential for proper inheritance of the mtDNA; also called kinetoplast or kDNA; as it forms a physical connection between the kDNA and the basal body of the single flagellum throughout the cell cycle. Thus, the presence of pATOM36 in the TAC provides an exciting link between mitochondrial protein import and kDNA inheritance.