Adrian Kolodzik

NAOMI: on the almost trivial task of reading molecules from different file formats.

In most cheminformatics workflows, chemical information is stored in files which provide the necessary data for subsequent calculations. The correct interpretation of the file formats is an important prerequisite to obtain meaningful results. Consistent reading of molecules from files, however, is not an easy task. Each file format implicitly represents an underlying chemical model, which has to be taken into consideration when the input data is processed. Additionally, many data sources contain invalid molecules. These have to be identified and either corrected or discarded. We present the chemical file format converter NAOMI, which provides efficient procedures for reliable handling of molecules from the common chemical file formats SDF, MOL2, and SMILES. These procedures are based on a consistent chemical model which has been designed for the appropriate representation of molecules relevant in the context of drug discovery. NAOMIs functionality is tested by round robin file IO exercises with public data sets, which we believe should become a standard test for every cheminformatics tool.

Reading PDB: perception of molecules from 3D atomic coordinates.

The analysis of small molecule crystal structures is a common way to gather valuable information for drug development. The necessary structural data is usually provided in specific file formats containing only element identities and three-dimensional atomic coordinates as reliable chemical information. Consequently, the automated perception of molecular structures from atomic coordinates has become a standard task in cheminformatics. The molecules generated by such methods must be both chemically valid and reasonable to provide a reliable basis for subsequent calculations. This can be a difficult task since the provided coordinates may deviate from ideal molecular geometries due to experimental uncertainties or low resolution. Additionally, the quality of the input data often differs significantly thus making it difficult to distinguish between actual structural features and mere geometric distortions. We present a method for the generation of molecular structures from atomic coordinates based on the recently published NAOMI model. By making use of this consistent chemical description, our method is able to generate reliable results even with input data of low quality. Molecules from 363 Protein Data Bank (PDB) entries could be perceived with a success rate of 98%, a result which could not be achieved with previously described methods. The robustness of our approach has been assessed by processing all small molecules from the PDB and comparing them to reference structures. The complete data set can be processed in less than 3 min, thus showing that our approach is suitable for large scale applications.

The Valence State Combination Model: A Generic Framework for Handling Tautomers and Protonation States

The consistent handling of molecules is probably the most basic and important requirement in the field of cheminformatics. Reliable results can only be obtained if the underlying calculations are independent of the specific way molecules are represented in the input data. However, ensuring consistency is a complex task with many pitfalls, an important one being the fact that the same molecule can be represented by different valence bond structures. In order to achieve reliability, a cheminformatics system needs to solve two fundamental problems. First, different choices of valence bond structures must be identified as the same molecule. Second, for each molecule all valence bond structures relevant to the context must be taken into consideration. The latter is especially important with regard to tautomers and protonation states, as these have considerable influence on physicochemical properties of molecules. We present a comprehensive method for the rapid and consistent generation of reasonable tautomers and protonation states for molecules relevant in the context of drug design. This method is based on a generic scheme, the Valence State Combination Model, which has been designed for the enumeration and scoring of valence bond structures in large data sets. In order to ensure our methods consistency, we have developed procedures which can serve as a general validation scheme for similar approaches. The analysis of both the average number of generated structures and the associated runtimes shows that our method is perfectly suited for typical cheminformatics applications. By comparison with frequently used and curated public data sets, we can demonstrate that the tautomers and protonation state produced by our method are chemically reasonable.

Unique Ring Families: A Chemically Meaningful Description of Molecular Ring Topologies

The perception of a set of rings forms the basis for a number of chemoinformatics applications, e.g. the systematic naming of compounds, the calculation of molecular descriptors, the matching of SMARTS expressions, and the generation of atomic coordinates. We introduce the concept of unique ring families (URFs) as an extension of the concept of relevant cycles (RCs). URFs are consistent for different atom orders and represent an intuitive description of the rings of a molecular graph. Furthermore, in contrast to RCs, URFs are polynomial in number. We provide an algorithm to efficiently calculate URFs in polynomial time and demonstrate their suitability for real-time applications by providing computing time benchmarks for the PubChem Database. URFs combine three important properties of chemical ring descriptions, for the first time, namely being unique, chemically meaningful, and efficient to compute. Therefore, URFs are a valuable alternative to the commonly used concept of the smallest set of smallest rings (SSSR) and would be suited to become the standard measure for ring topologies of small molecules.

In silico Design, Synthesis, and Screening of Novel Deoxyhypusine Synthase Inhibitors Targeting HIV-1 Replication

The human enzyme deoxyhypusine synthase (DHS) is an important host cell factor that participates in the post‐translational hypusine modification of eukaryotic initiation factor 5A (eIF‐5A). Hypusine‐modified eIF‐5A plays a role in a number of diseases, including HIV infection/AIDS. Thus, DHS represents a novel and attractive drug target. So far, four crystal structures are available, and various substances have been tested for inhibition of human DHS. Among these inhibitors, N‐1‐guanyl‐1,7‐diaminoheptane (GC7) has been co‐crystallized in the active site of DHS. However, despite its potency, GC7 is not selective enough to be used in drug applications. Therefore, new compounds that target DHS are needed. Herein we report the in silico design, chemical synthesis, and biological evaluation of new DHS inhibitors. One of these inhibitors showed dose‐dependent inhibition of DHS in vitro, as well as suppression of HIV replication in cell cultures. Furthermore, the compound exhibited no cytotoxic effects at active concentrations. Thus, this designed compound demonstrated proof of principle and represents a promising starting point for the development of new drug candidates to specifically interfere with DHS activity.

Linker-Region Modified Derivatives of the Deoxyhypusine Synthase Inhibitor CNI-1493 Suppress HIV-1 Replication

The inhibition of cellular factors that are involved in viral replication may be an important alternative to the commonly used strategy of targeting viral enzymes. The guanylhydrazone CNI‐1493, a potent inhibitor of the deoxyhypusine synthase (DHS), prevents the activation of the cellular factor eIF‐5A and thereby suppresses HIV replication and a number of other diseases. Here, we report on the design, synthesis and biological evaluation of a series of CNI‐1493 analogues. The sebacoyl linker in CNI‐1493 was replaced by different alkyl or aryl dicarboxylic acids. Most of the tested derivatives suppress HIV‐1 replication efficiently in a dose‐dependent manner without showing toxic side effects. The unexpected antiviral activity of the rigid derivatives point to a second binding mode as previously assumed for CNI‐1493. Moreover, the chemical stability of CNI‐1493 was analysed, showing a successive hydrolysis of the imino bonds. By molecular dynamics simulations, the behaviour of the parent CNI‐1493 in solution and its interactions with DHS were investigated.