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.

Gibbs Energy of Formation of Ca3Ti8Al12O37 and Phase Relations and Chemical Potentials in the System Al2O3-TiO2-CaO

The quaternary oxide in the system Al2O3-CaO-TiO2 is found to have the composition Ca3Ti8Al12O37 rather than CaTi3Al8O19 as reported in the literature. The standard Gibbs energy of formation of Ca3Ti8Al12O37 from component binary oxides is measured in the temperature range from 900 to 1250 K using a solid-state electrochemical cell incorporating single crystal CaF2 as the solid electrolyte. The results can be represented by the equation: delta G(f(0x))(0) (+/- 70)/J mol(-1) = -248474 - 15.706(T/K). Combining this information with thermodynamic data on calcium aluminates and titanates available in the literature, subsolidus phase relations in the pseudo-ternary system Al2O3-CaO-TiO2 are computed and presented as isothermal sections. The evolution of phase relations with temperature is highlighted. Chemical potential diagrams are computed at 1200 K, showing the stability domains of the various phases in the chemical potential-composition space. In each chemical potential diagram, chemical potential of one component is plotted against the cationic fraction of the other two components. The diagrams are valid at relatively high oxygen potentials where Ti is present in its four-valent state in all the oxide phases.

System Ho-Rh-O: Phase Equilibria, Chemical Potentials and Gibbs Energy of Formation of HoRhO3

The thermodynamic properties of the HoRhO3 were determined in the temperature range from 900 to 1300 K by using a solid-state electrochemical cell incorporating calcia-stabilized zirconia as the electrolyte. The standard Gibbs free energy of formation of orthorhombic perovskite HoRhO3, from Ho2O3 with C-rare earth structure and Rh2O3 with orthorhombic structure, can be expressed by the equation; Delta G(f)degrees((ox)) (+/- 78)/(J/mol) = -50535 + 3.85(T/K) Using the thermodynamic data of HoRhO3 and auxiliary data for binary oxides from the literature, the phase relations in the Ho-Rh-O system were computed at 1273 K. Thermodynamic data for intermetallic phases in the binary Ho-Rh were estimated from experimental enthalpy of formation for three compositions from the literature and Miedema's model, consistent with the phase diagram. The oxygen potential-composition diagram and three-dimensional chemical potential diagram at 1273 K, and temperature-composition diagrams at constant oxygen partial pressures were computed for the system Ho-Rh-O. The decomposition temperature of HoRhO3 is 1717(+/- 2) K in pure O-2 and 1610(+/- 2) K in air at a total pressure p(o) = 0.1 MPa.

Thermodynamics of TmRhO3, Phase Equilibria, and Chemical Potentials in the System Tm-Rh-O

Phase equilibria in the system Tm-Rh-O at 1200 K is established by isothermal equilibration of selected compositions and phase identification after quenching to room temperature. Six intermetallic phases (Tm3Rh, Tm7Rh3, Tm5Rh3, Tm3Rh2, TmRh, TmRh2 +/-delta) and a ternary oxide TmRhO3 are identified. Based on experimentally determined phase relations, a solid-state electrochemical cell is devised to measure the standard free energy of formation of orthorhombic perovskite TmRhO3 from cubic Tm2O3 and beta-Rh2O3 in the temperature range from (900 to 1300) K. The results can be summarized as: Delta G(f,ox)(o) +/- 104/J.mol(-1) = -46474 + 3.925(T/K). Invoking the Neumann-Kopp rule, the standard enthalpy of formation of TmRhO3 from its constituent elements at 298.15 K is estimated as -1193.89 (+/- 2.86) kJ.mol(-1). The standard entropy of TmRhO3 at 298.15 K is evaluated as 103.8 (+/- 1.6) J.mol(-1).K-1. The oxygen potential-composition diagram and three-dimensional chemical potential diagram at 1200 K and temperature-composition diagrams at constant partial pressures of oxygen are computed from thermodynamic data. The compound TmRhO3 decomposes at 1688 (+/- 2) K in pure oxygen and at 1583 (+/- 2) K in air at standard pressure.

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.

Calciothermic reduction of TiO2: A diagrammatic assessment of the thermodynamic limit of deoxidation

Calciothermic reduction of TiO2 provides a potentially low-cost route to titanium production. Presented in this article is a suitably designed diagram, useful for assessing the degree of reduction of TiO2 and residual oxygen contamination in metal as a function of reduction temperature and other process parameters. The oxygen chemical potential diagram à la Ellingham-Richardson-Jeffes is useful for visualization of the thermodynamics of reduction reactions at high temperatures. Although traditionally the diagram depicts oxygen potentials corresponding to the oxidation of different metals to their corresponding oxides or of lower oxides to higher oxides, oxygen potentials associated with solution phases at constant composition can be readily superimposed. The usefulness of the diagram for an insightful analysis of calciothermic reduction, either direct or through an electrochemical process, is discussed. Identified are possible process variations, modeling and optimization strategies.

Electrochemical deoxidation of RE–O (RE=Gd, Tb, Dy, Er) solid solutions

The removal of oxygen from rare-earth metals (RE, RE=Gd, Tb, Dy, Er) by an electrochemical deoxidation method was investigated. A titanium basket containing the rare-earth metal sample, submerged in molten CaCl2 electrolyte, formed the cathode of an electrolysis cell. A high-purity graphite anode was used. The calcium metal produced at the cathode effectively deoxidized the rare-earth metal. Carbon monoxide and dioxide were generated at the graphite anode. Rare-earth metals containing more than 2000 mass ppm oxygen were deoxidized to 10–50 mass ppm level by electrolysis at 1189 K for 36 ks (10 h). Cyclic voltammetry was used to characterize the molten salt at different stages of the process. The effectiveness of the process is discussed with the aid of a chemical potential diagram for RE–O solid solutions. The new electrochemical technique is compared with the conventional deoxidation methods reported in the literature. The possibility of nitrogen removal from the rare-earth metals by the electrochemical method is outlined.

Production of niobium powder by electronically mediated reaction (EMR) using calcium as a reductant

Explored in this study is an electronically mediated reaction (EMR) route for the production of niobium powder using calcium as a reductant for niobium oxide (Nb2O5). Feed material, Nb2O5, and reductant calcium alloy containing aluminum and nickel were charged into electronically isolated locations in a molten salt (e.g. CaCl2) at 1173 K. The current flow through an external path between the feed and reductant locations was monitored. A current approximately 0.4 A was measured during the reaction in the external circuit connecting cathode and anode location. Niobium powder with low aluminum and nickel content was obtained although liquid Ca–Al–Ni alloy was used as the reductant. This clearly demonstrates that niobium metal powder can be produced by an electronically mediated reaction (EMR), without direct physical contact between feed (Nb2O5) and reductant (calcium). Mechanism of calciothermic reduction of Nb2O5 in the molten salt is discussed using an isothermal chemical potential diagram.

Solute-solute and solute-solvent interactions in transition metal alloys: Pt-Ti system

An isothermal section of the phase diagram for (silver + rhodium + oxygen) at T = 1173 K has been established by equilibration of samples representing twelve different compositions, and phase identification after quenching by optical and scanning electron microscopy (s.e.m.), X-ray diffraction (x.r.d.), and energy dispersive analysis of X-rays (e.d.x.), Only one ternary oxide, AgRhO2, was found to be stable and a three phase region involving Ag, AgRhO2 and Rh2O3 was identified. The thermodynamic properties of AgRhO2 were measured using a galvanic cell in the temperature range 980 K to 1320 K. Yttria-stabilized zirconia was used as the solid electrolyte and pure oxygen gas at a pressure of 0.1 MPa was used as the reference electrode. The Gibbs free energy of formation of the ternary oxide from the elements, ΔfGo (AgRhO2), can be represented by two linear equations that join at the melting temperature of silver. In the temperature range 980 K to 1235 K, ΔfGo(AgRhO2)/(J . mol-1) = -249080 + 179.08 T/K (±120). Above the melting temperature of silver, in the temperature range 1235 K to 1320 K, ΔfGo(AgRhO2)/(J . mol-1) = -260400 + 188.24 T/K (±95). The thermodynamic properties of AgRhO2 at T = 298.15 K were evaluated from the high temperature data. The chemical potential diagram for (silver + rhodium + oxygen) at T = 1200 K was also computed on the basis of the results of this study.

Electrochemical determination of thermodynamic properties of DyRhO3 and phase relations in the system Dy-Rh-O

Thermodynamic properties of Dysprosium rhodite (DyRhO3) are measured in the temperature range from 900 to 1,300 K using a solid-state electrochemical cell incorporating yttria-stabilized zirconia as the electrolyte. The standard Gibbs free energy of formation of DyRhO3 with O-type perovskite structure from its components binary oxides, Dysprosia with C-rare earth structure and beta-Rh2O3 with orthorhombic structure, can be represented by the equation: Delta G(f(OX))(O) (+/- 182)/J mol(-1) = -52710+3.821(T/K). By using the thermodynamic data for DyRhO3 from experiment and auxiliary data for other phases from the literature, the phase relations in the system Dy-Rh-O are computed. Thermodynamic data for intermetallic phases in the binary system Dy-Rh, required for constructing the chemical potential diagrams, are evaluated using calorimetric data available in the literature for three intermetallics and Miedema's model, consistent with the phase diagram. The results are presented in the form of Gibbs triangle, oxygen potential-composition diagram, and three-dimensional chemical potential diagram at 1,273 K. Temperature-composition diagrams at constant oxygen partial pressures are also developed. The decomposition temperature of DyRhO3 is 1,732 (+/- 2.5) K in pure oxygen and 1,624 (+/- 2.5) K and in air at standard pressure.

System Gd-Rh-O: Thermodynamics and phase relations

Thermodynamic properties of GdRhO3 are investigated in the temperature range from 900 to 1300 K by employing a solid-state electrochemical cell, incorporating calcia-stabilized zirconia as the electrolyte. The standard Gibbs free energy of formation of GdRhO3 from component binary oxide Gd2O3 with C-rare earth structure and Rh2O3 with orthorhombic structure can be expressed as; Delta G(f(ox))(o)(+/- 60)/J mol(-1) = -56603 + 3.78(T/K) Based on the thermodynamic information on GdRhO3 from experiment and auxiliary data for binary oxides from the literature and estimated properties of Gd-Rh alloys, phase relations are computed for the system Gd-Rh-O at 1273 K. Gibbs free energies for intermetallic phases in the binary Gd-Rh are evaluated using calorimetric data available in the literature for two compositions and Miedema's model, consistent with the binary phase diagram. Isothermal section of the ternary phase diagram, oxygen potential-composition diagram and a 3-D chemical potential diagram for the system Gd-Rh-O at 1273 K are developed. Phase relations in the ternary Gd-Rh-O are also computed as a function of temperature at constant oxygen partial pressures. The ternary oxide, GdRhO3 decomposes to Gd2O3 with B-rare earth structure, metallic Rh and O-2 at 1759(+/- 2) K in pure O-2 and 1649(+/- 2) K in air at a total pressure P-0 -0.1 MPa. (c) 2012 Elsevier B.V. All rights reserved.