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.

Aqueous solubilities of five N-(diethylaminothiocarbonyl)benzimido derivatives at T = 298.15 K

N-(diethylaminothiocarbonyl)benzimido derivatives are polar multifunctional substances. A set of these compounds was synthesised by successive substitution on the enamine side, resulting in similar substances with different polarities, providing a set of model compounds with respect to the study of substituent effects on physico-chemical properties. Experimental aqueous solubility data, at T = 298.15 K, of N-(diethylaminothiocarbonyl)benzamidine, PhCNH2NCSNEt2 (1),N-(diethylaminothiocarbonyl)-N'-phenylbenzamidine, PhCNHPhNCSNEt2 (2), N-(diethylaminothiocarbonyl)-N'-monoethylbenzamidine, PhCNHEtNCSNEt2 (3), N-(diethylaminothiocarbonyl)-N',N'-diethylbenzamidine, PhCNEt2NCSNEt2 (4), and N-(diethylaminothiocarbonyl)benzimido ethylester, PhCOEtNCSNEt2 (5) were measured at T = 298.15 K. The obtained data are supplemented by COSMO-RS aqueous solubility predictions as well as other environmentally important partition coefficients. This information is shown in a two-dimensional chemical space diagram, providing indications about the compartment into which the bulk of the compounds is likely to concentrate. The expected quality of COSMO-RS predictions for this type of screening exercise is illustrated on a set of pesticides with established thermophysical property data.

Operating space diagrams: a tool for designers of wastewater treatment plants

Recent years have seen the introduction of new and varied designs of activated sludge plants. With increasing needs for higher efficiencies and lower costs, the possibility of a plant that operates more effectively has created the need for tools that can be used to evaluate and compare designs at the design stage. One such tool is the operating space diagram. It is the aim of this paper to present this tool and demonstrate its application and relevance to design using a simple case study. In the case study, use of the operating space diagram suggested changes in design that would improve the flexibility of the process. It also was useful for designing suitable control strategies.

Exploring the space of histidine containing dipeptides in search of novel efficient RCS sequestering agents.

The study reports a set of forty proteinogenic histidine-containing dipeptides as potential carbonyl quenchers. The peptides were chosen to cover as exhaustively as possible the accessible chemical space, and their quenching activities toward 4-hydroxy-2-nonenal (HNE) and pyridoxal were evaluated by HPLC analyses. The peptides were capped at the C-terminus as methyl esters or amides to favor their resistance to proteolysis and diastereoisomeric pairs were considered to reveal the influence of configuration on quenching. On average, the examined dipeptides are less active than the parent compound carnosine (βAla + His) thus emphasizing the unfavorable effect of the shortening of the βAla residue as confirmed by the control dipeptide Gly-His. Nevertheless, some peptides show promising activities toward HNE combined with a remarkable selectivity. The results emphasize the beneficial role of aromatic and positively charged residues, while negatively charged and H-bonding side chains show a detrimental effect on quenching. As a trend, ester derivatives are slightly more active than amides while heterochiral peptides are more active than their homochiral diastereoisomer. Overall, the results reveal that quenching activity strongly depends on conformational effects and vicinal residues (as evidenced by the reported QSAR analysis), offering insightful clues for the design of improved carbonyl quenchers and to rationalize the specific reactivity of histidine residues within proteins.

Generation of Aptamers with an Expanded Chemical Repertoire

The enzymatic co-polymerization of modified nucleoside triphosphates (dN*TPs and N*TPs) is a versatile method for the expansion and exploration of expanded chemical space in SELEX and related combinatorial methods of in vitro selection. This strategy can be exploited to generate aptamers with improved or hitherto unknown properties. In this review, we discuss the nature of the functionalities appended to nucleoside triphosphates and their impact on selection experiments. The properties of the resulting modified aptamers will be described, particularly those integrated in the fields of biomolecular diagnostics, therapeutics, and in the expansion of genetic systems (XNAs).

DNA Catalysis: The Chemical Repertoire of DNAzymes

Deoxyribozymes or DNAzymes are single-stranded catalytic DNA molecules that are obtained by combinatorial in vitro selection methods. Initially conceived to function as gene silencing agents, the scope of DNAzymes has rapidly expanded into diverse fields, including biosensing, diagnostics, logic gate operations, and the development of novel synthetic and biological tools. In this review, an overview of all the different chemical reactions catalyzed by DNAzymes is given with an emphasis on RNA cleavage and the use of non-nucleosidic substrates. The use of modified nucleoside triphosphates (dN*TPs) to expand the chemical space to be explored in selection experiments and ultimately to generate DNAzymes with an expanded chemical repertoire is also highlighted.

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.

New chemical modalities for leishmaniasis drug discovery

Leishmaniasis is one of the major parasitic diseases among neglected tropical diseases with a high rate of morbidity and mortality. Human migration and climate change have spread the disease from limited endemic areas all over the world, also reaching regions in Southern Europe, and causing significant health and economic burden. The currently available treatments are far from ideal due to host toxicity, elevated cost, and increasing rates of drug resistance. Safer and more effective drugs are thus urgently required. Nevertheless, the identification of new chemical entities for leishmaniasis has proven to be incredibly hard and exacerbated by the scarcity of well-validated targets. Trypanothione reductase (TR) represents one robustly validated target in Leishmania that fulfils most of the requirements for a good drug target. However, due to the large and featureless active site, TR is considered extremely challenging and almost undruggable by small molecules. This scenario advocates the development of new chemical entities by unlocking new modalities for leishmaniasis drug discovery. The classical toolbox for drug discovery has enormously expanded in the last decade, and medicinal chemists can now strategize across a variety of new chemical modalities and a vast chemical space, to efficiently modulate challenging targets and provide effective treatments. Beyond others, Targeted p Protein Degradation (TPD) is an emerging strategy that uses small molecules to hijack endogenous proteolysis systems to degrade disease-relevant proteins and thus reduce their abundance in the cell. Based on these considerations, this thesis aimed to develop new strategies for leishmaniasis drug discovery while embracing novel chemical modalities and navigating the chemical space by chasing unprecedented chemotypes. This has been achieved by four complementary projects. We believe that these next-generation chemical modalities for leishmaniasis will play an important role in what was previously thought to be a drug discovery landscape dominated by small molecules.

A hedging point strategy - balancing effluent quality, economy and robustness in the control of wastewater treatment plants

An operational space map is an efficient tool to compare a large number of operational strategies to find an optimal choice of setpoints based on a multicriterion. Typically, such a multicriterion includes a weighted sum of cost of operation and effluent quality. Due to the relative high cost of aeration such a definition of optimality result in a relatively high fraction of the effluent total nitrogen in the form of ammonium. Such a strategy may however introduce a risk into operation because a low degree of ammonium removal leads to a low amount of nitrifiers. This in turn leads to a reduced ability to reject event disturbances, such as large variations in the ammonium load, drop in temperature, the presence of toxic/inhibitory compounds in the influent etc. Hedging is a risk minimisation tool, with the aim to "reduce one's risk of loss on a bet or speculation by compensating transactions on the other side" (The Concise Oxford Dictionary (1995)). In wastewater treatment plant operation hedging can be applied by choosing a higher level of ammonium removal to increase the amount of nitrifiers. This is a sensible way to introduce disturbance rejection ability into the multi criterion. In practice, this is done by deciding upon an internal effluent ammonium criterion. In some countries such as Germany, a separate criterion already applies to the level of ammonium in the effluent. However, in most countries the effluent criterion applies to total nitrogen only. In these cases, an internal effluent ammonium criterion should be selected in order to secure proper disturbance rejection ability.

Development of CO-releasing molecular for the treatment of inflammatory diseases

Dissertation presented to obtain a Ph.D. degree (Doutoramento) in Chemistry at the Instituto de Tecnologia Quimica e Biol6gica da Universidade Nova de Lisboa

A platform for anti-biofilm assays combining viability, biomass and matrix quantifications in susceptibility assessments of antimicrobials against Staphylococcus aureus biofilms

Bacteria can exist as planktonic, the lifestyle in which single cells exist in suspension, and as biofilms, which are surface-attached bacterial communities embedded in a selfproduced matrix. Most of the antibiotics and the methods for antimicrobial work have been developed for planktonic bacteria. However, the majority of the bacteria in natural habitats live as biofilms. Biofilms develop dauntingly fast high resistance towards conventional antibacterial treatments and thus, there is a great need to meet the demands of effective anti-biofilm therapy. In this thesis project it was attempted to fill the void of anti-biofilm screening methods by developing a platform of assays that evaluate the effect that screened compounds have on the total biomass, viability and the extracellular polysaccharide (EPS) layer of the biofilms. Additionally, a new method for studying biofilms and their interactions with compounds in a continuous flow system was developed using capillary electrochromatography (CEC). The screening platform was utilized with a screening campaign using a small library of cinchona alkaloids. The assays were optimized to be statistically robust enough for screening. The first assay, based on crystal violet staining, measures total biofilm biomass, and it was automated using a liquid handling workstation to decrease the manual workload and signal variation. The second assay, based on resazurin staining, measures viability of the biofilm, and it was thoroughly optimized for the strain used, but was then a very simple and fast method to be used for primary screening. The fluorescent resazurin probe is not toxic to the biofilms. In fact, it was also shown in this project that staining the biofilms with resazurin prior to staining with crystal violet had no effect on the latter and they can be used in sequence on the same screening plate. This sequential addition step was indeed a major improvement on the use of reagents and consumables and also shortened the work time. As a third assay in the platform a wheat germ agglutinin based assay was added to evaluate the effect a compound has on the EPS layer. Using this assay it was found that even if compounds might have clear effect on both biomass and viability, the EPS layer can be left untouched or even be increased. This is a clear implication of the importance of using several assays to be able to find “true hits” in a screening setting. In the pilot study of screening for antimicrobial and anti-biofilm effects using a cinchona alkaloid library, one compound was found to have antimicrobial effect against planktonic bacteria and prevent biofilm formation at low micromolar concentration. To eradicate biofilms, a higher concentration was needed. It was also shown that the chemical space occupied by the active compound was slightly different than the rest of the cinchona alkaloids as well as the rest of the compounds used for validatory screening during the optimization processes of the separate assays.

Biomolecular screening for inhibitors of butyrylcholinesterase : identification and characterization using in vitro and in silico tools

Drug discovery is a continuous process where researchers are constantly trying to find new and better drugs for the treatment of various conditions. Alzheimer’s disease, a neurodegenerative disease mostly affecting the elderly, has a complex etiology with several possible drug targets. Some of these targets have been known for years while other new targets and theories have emerged more recently. Cholinesterase inhibitors are the major class of drugs currently used for the symptomatic treatment of Alzheimer’s disease. In the Alzheimer’s disease brain there is a deficit of acetylcholine and an impairment in signal transmission. Acetylcholinesterase has therefore been the main target as this is the main enzyme hydrolysing acetylcholine and ending neurotransmission. It is believed that by inhibiting acetylcholinesterase the cholinergic signalling can be enhanced and the cognitive symptoms that arise in Alzheimer’s disease can be improved. Butyrylcholinesterase, the second enzyme of the cholinesterase family, has more recently attracted interest among researchers. Its function is still not fully known, but it is believed to play a role in several diseases, one of them being Alzheimer’s disease. In this contribution the aim has primarily been to identify butyrylcholinesterase inhibitors to be used as drug molecules or molecular probes in the future. Both synthetic and natural compounds in diverse and targeted screening libraries have been used for this purpose. The active compounds have been further characterized regarding their potencies, cytotoxicity, and furthermore, in two of the publications, the inhibitors ability to also inhibit Aβ aggregation in an attempt to discover bifunctional compounds. Further, in silico methods were used to evaluate the binding position of the active compounds with the enzyme targets. Mostly to differentiate between the selectivity towards acetylcholinesterase and butyrylcholinesterase, but also to assess the structural features required for enzyme inhibition. We also evaluated the compounds, active and non-active, in chemical space using the web-based tool ChemGPS-NP to try and determine the relevant chemical space occupied by cholinesterase inhibitors. In this study, we have succeeded in finding potent butyrylcholinesterase inhibitors with a diverse set of structures, nine chemical classes in total. In addition, some of the compounds are bifunctional as they also inhibit Aβ aggregation. The data gathered from all publications regarding the chemical space occupied by butyrylcholinesterase inhibitors we believe will give an insight into the chemically active space occupied by this type of inhibitors and will hopefully facilitate future screening and result in an even deeper knowledge of butyrylcholinesterase inhibitors.

Calculs ab initio de structures électroniques pour un meilleur design de polymères photovoltaïques

La présente thèse porte sur l'utilité de la théorie de la fonctionnelle de la densité dans le design de polymères pour applications photovoltaïques. L'étude porte d'abord sur le rôle des calculs théoriques pour la caractérisation des polymères dans le cadre de collaborations entre la théorie et l'expérience. La stabilité et les niveaux énergétiques de certaines molécules organiques sont étudiés avant et après la sulfuration de leurs groupements carbonyles, un procédé destiné à diminuer le band gap. Les propriétés de dynamique électronique, de séparation des porteurs de charges et de spectres de vibrations Raman sont également explorées dans un polymère à base de polycarbazole. Par la suite, l'utilité des calculs théoriques dans le design de polymères avant leurs synthèses est considérée. La théorie de la fonctionnelle de la densité est étudiée dans le cadre du modèle de Scharber afin de prédire l'efficacité des cellules solaires organiques. Une nouvelle méthode de design de polymères à faible band gaps, basée sur la forme structurale aromatique ou quinoide est également présentée, dont l'efficacité surpasse l'approche actuelle de donneur-accepteur. Ces études sont mises à profit dans l'exploration de l'espace moléculaire et plusieurs candidats de polymères aux propriétés électroniques intéressantes sont présentés.

Self-Organizing Maps of Molecular Descriptors for Sesquiterpene Lactones and Their Application to the Chemotaxonomy of the Asteraceae Family

The Asteraceae, one of the largest families among angiosperms, is chemically characterised by the production of sesquiterpene lactones (SLs). A total of 1,111 SLs, which were extracted from 658 species, 161 genera, 63 subtribes and 15 tribes of Asteraceae, were represented and registered in two dimensions in the SISTEMATX, an in-house software system, and were associated with their botanical sources. The respective 11 block of descriptors: Constitutional, Functional groups, BCUT, Atom-centred, 2D autocorrelations, Topological, Geometrical, RDF, 3D-MoRSE, GETAWAY and WHIM were used as input data to separate the botanical occurrences through self-organising maps. Maps that were generated with each descriptor divided the Asteraceae tribes, with total index values between 66.7% and 83.6%. The analysis of the results shows evident similarities among the Heliantheae, Helenieae and Eupatorieae tribes as well as between the Anthemideae and Inuleae tribes. Those observations are in agreement with systematic classifications that were proposed by Bremer, which use mainly morphological and molecular data, therefore chemical markers partially corroborate with these classifications. The results demonstrate that the atom-centred and RDF descriptors can be used as a tool for taxonomic classification in low hierarchical levels, such as tribes. Descriptors obtained through fragments or by the two-dimensional representation of the SL structures were sufficient to obtain significant results, and better results were not achieved by using descriptors derived from three-dimensional representations of SLs. Such models based on physico-chemical properties can project new design SLs, similar structures from literature or even unreported structures in two-dimensional chemical space. Therefore, the generated SOMs can predict the most probable tribe where a biologically active molecule can be found according Bremer classification.

Identifizierung chemischer Transportinhibitoren zur Modulation (patho)biologischer Genprodukte

Die intrazelluläre Lokalisation von Proteinen und Makromolekülen unterliegt in Eukaryoten einer strengen Regulation. Insbesondere erlaubt die Kompartimentierung eukaryotischer Zellen in Zellkern und Zytoplasma den simultanen Ablauf räumlich getrennter biochemischer Reaktionen, und damit die unabhängige Regulation zellulärer Programme. Da trotz intensiver Forschungsbemühungen bis dato die molekularen Details sowie die (patho)biologische Bedeutung von Kern-Zytoplasma-Transportprozessen noch immer nicht vollkommen verstanden sind, wurde im Rahmen der vorliegenden Arbeit ein Fokus auf die Identifizierung von chemischen Transportinhibitoren gelegt. Das zu diesem Zweck entwickelte Translokations-Biosensor-System basiert auf der Kombination von autofluoreszierenden Proteinen, sowie spezifisch ausgewählten Kernexport- und Kernimportsignalen. Nach Etablierung geeigneter Zellmodelle, die effizient und stabil die Translokations-Biosensoren exprimieren, wurde die 17 000 Substanzen umfassende Bibliothek der ChemBioNet-Initiative nach Kernexportinhibitoren mittels einer Fluoreszenzmikroskopie-basierten Hochdurchsatzanalyse-Plattform durchmustert. Zunächst wurden Translokations-Algorithmen, welche eine zuverlässige automatisierte Erkennung von Zellkern und Zytoplasma erlauben, optimiert. Im Folgenden konnten acht neue niedermolekulare Kernexport-Inhibitoren identifiziert werden, die sich in der Stärke, der Geschwindigkeit, sowie in der Beständigkeit der vermittelten Inhibition unterscheiden. Die Aktivität der Inhibitoren konnte auf den isolierten nukleären Exportsignalen (NES) von HIV-1 Rev und Survivin als auch auf den entsprechenden Volllängeproteinen mittels Mikroinjektionsexperimenten sowie durch umfassende in vitro und biochemische Methoden bestätigt werden. Zur Untersuchung der funktionellen Einheiten der Inhibitoren wurden homologe Substanzen auf Ihre Aktivität hin getestet. Dabei konnten für die Aktivität wichtige chemische Gruppen definiert werden. Alle Substanzen stellen neue Inhibitoren des Crm1-abhängigen Exports dar und zeigen keine nachweisbare NES-Selektivität. Interessanterweise konnte jedoch eine zytotoxische und Apoptose-induzierende Wirkung auf verschiedene Krebszellarten festgestellt werden. Da diese Wirkung unabhängig vom p53-Status der Tumorzellen ist und die Inhibitoren C3 und C5 die Vitalität nicht-maligner humaner Zellen signifikant weniger beeinträchtigen, wurden diese Substanzen zum internationalen Patent angemeldet. Da der nukleäre Export besonders für Tumorzellen einen wichtigen Überlebenssignalweg darstellt, könnte dessen reversible Hemmung ausgenutzt werden, um besonders in Kombination mit gängigen Krebstherapien eine therapeutisch relevante Tumorinhibition zu erzeugen. Eine weitere Anwendungsmöglichkeit der neuen Exportinhibitoren ist auf dem Gebiet der Infektionskrankheiten zu sehen, da auch die Aktivität des essentiellen HIV-1 Rev-Proteins inhibiert wird. Zusätzlich konnte in der Arbeit gezeigt werden, dass der zelluläre Kofaktor des Crm1-abhängigen Exports des HIV-1 Rev-Proteins, die RNA-Helikase DDX3, ein eigenes NES enthält. Der Nachweis einer direkten Interaktion des HIV-1 Rev- mit dem DDX3-Protein impliziert, dass multiple Angriffstellen für chemische Modulatoren hinsichtlich einer antiviralen Therapie gegeben sind. Da die Vielfalt des chemischen Strukturraums es unmöglich macht diesen experimentell vollständig zu durchmustern, wurden im Rahmen dieser Arbeit auch Naturstoffe als vielversprechende Wirkstoffquelle untersucht. Um zukünftig umfassend bioaktive Substanzen aus diesen hochkomplexen Stoffgemischen experimentell identifizieren zu können, wurde eine Fluoreszenzmikroskopie-basierte Hochdurchsatzanalyse-Plattform am Mainz Screening Center (MSC) etabliert. Damit konnte bereits ein weiterer, bisher unbekannter Exportinhibitor aus Cyphellopsis anomala identifiziert werden. Neben einer Anwendung dieser Substanz als chemisches Werkzeug zur Aufklärung der Regulation von Transportvorgängen, stellt sich auch die evolutionsbiologisch relevante Frage, wie es dem Pilzproduzenten gelingt die Blockierung des eigenen Kernexports zu umgehen. Weiterführende Projekte müssen sich neben der Aufklärung der molekularen Wirkmechanismen der gefundenen Substanzen mit der Identifizierung spezifischer chemischer „Funktionseinheiten“ beschäftigen. Neben einem verbesserten mechanistischen Verständnis von Transportvorgängen stellen die erarbeiteten Transportinhibitoren Vorstufen zur Weiterentwicklung möglicher Wirkstoffe dar. Die im Rahmen dieser Arbeit etablierte Technologie-Plattform und molekularen Werkzeuge stellen darüber hinaus eine wichtige Voraussetzung dar, um eine systematische Suche nach möglichen Wirkstoffen im Forschungsfeld der „Chemischen Biomedizin“ voranzutreiben.