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. NAOMI's 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.

De novo design by pharmacophore-based searches in fragment spaces.

De novo ligand design supports the search for novel molecular scaffolds in medicinal chemistry projects. This search can either be based on structural information of the targeted active site (structure-based approach) or on similarity to known binders (ligand-based approach). In the absence of structural information on the target, pharmacophores provide a way to find topologically novel scaffolds. Fragment spaces have proven to be a valuable source for molecular structures in de novo design that are both diverse and synthetically accessible. They also offer a simple way to formulate custom chemical spaces. We have implemented a new method which stochastically constructs new molecules from fragment spaces under consideration of a three dimensional pharmacophore. The program has been tested on several published pharmacophores and is shown to be able to reproduce scaffold hops from the literature, which resulted in new chemical entities.

Fast automated placement of polar hydrogen atoms in protein-ligand complexes.

BACKGROUND: Hydrogen bonds play a major role in the stabilization of protein-ligand complexes. The ability of a functional group to form them depends on the position of its hydrogen atoms. An accurate knowledge of the positions of hydrogen atoms in proteins is therefore important to correctly identify hydrogen bonds and their properties. The high mobility of hydrogen atoms introduces several degrees of freedom: Tautomeric states, where a hydrogen atom alters its binding partner, torsional changes where the position of the hydrogen atom is rotated around the last heavy-atom bond in a residue, and protonation states, where the number of hydrogen atoms at a functional group may change. Also, side-chain flips in glutamine and asparagine and histidine residues, which are common crystallographic ambiguities must be identified before structure-based calculations can be conducted. RESULTS: We have implemented a method to determine the most probable hydrogen atom positions in a given protein-ligand complex. Optimality of hydrogen bond geometries is determined by an empirical scoring function which is used in molecular docking. This allows to evaluate protein-ligand interactions with an established model. Also, our method allows to resolve common crystallographic ambiguities such as as flipped amide groups and histidine residues. To ensure high speed, we make use of a dynamic programming approach. CONCLUSION: Our results were checked against selected high-resolution structures from an external dataset, for which the positions of the hydrogen atoms have been validated manually. The quality of our results is comparable to that of other programs, with the advantage of being fast enough to be applied on-the-fly for interactive usage or during score evaluation.