We have developed a fast grid-based algorithm, BRUTUS, for rigid-body molecular superposition and similarity searching. BRUTUS aligns molecules using field information derived from charge distributions and van der Waals shapes of the compounds. Molecules can have similar biological properties if their charge distributions and shapes are similar, even though they have different chemical structures; that is, BRUTUS can identify compounds possessing similar properties, regardless of their structures. In this paper, we present two applications of BRUTUS. First, BRUTUS was used to superimpose five sets of inhibitors. Second, two sets of known inhibitors were searched from a database, and the results were analyzed using self-organizing maps. We demonstrate that BRUTUS is successful in superimposing compounds using molecular fields and, importantly, is fast and accurate enough for virtual screening of chemical databases using a standard personal computer. This fast and efficient molecular-field-based algorithm is applicable for virtual screening of structurally diverse, active molecules.

The human constitutive androstane receptor (CAR, NR1I3) is an important regulator of xenobiotic metabolism and other physiological processes. So far, only few CAR agonists are known and no explicit mechanism has been proposed for their action. Thus, we aimed to generate a 3D QSAR model that could explain the molecular determinants of CAR agonist action. To obtain a sufficient number of agonists that cover a wide range of activity, we applied a virtual screening approach using both structure- and ligand-based methods. We identified 27 novel human CAR agonists on which a 3D QSAR model was generated. The model, complemented by coregulator recruitment and mutagenesis results, suggests a potential activation mechanism for human CAR and may serve to predict potential activation of CAR for compounds emerging from drug development projects or for chemicals undergoing toxicological risk assessment.

A panel of 92 catechol-O-methyltransferase (COMT) inhibitors was used to examine the molecular interactions affecting their biological activity. COMT inhibitors are used as therapeutic agents in the treatment of Parkinson's disease, but there are limitations in the currently marketed compounds due to adverse side effects. This study combined molecular docking methods with three-dimensional structure-activity relationships (3D QSAR) to analyse possible interactions between COMT and its inhibitors, and to incite the design of new inhibitors. Comparative molecular field analysis (CoMFA) and GRID/GOLPE models were made by using bioactive conformations from docking experiments, which yielded q2 values of 0.594 and 0.636, respectively. The docking results, the COMT X-ray structure, and the 3D QSAR models are in agreement with each other. The models suggest that an interaction between the inhibitor's catechol oxygens and the Mg2+ ion in the COMT active site is important. Both hydrogen bonding with Lys144, Asn170 and Glu199, and hydrophobic contacts with Trp38, Pro174 and Leu198 influence inhibitor binding. Docking suggests that a large R1 substituent of the catechol ring can form hydrophobic contacts with side chains of Val173, Leu198, Met201 and Val203 on the COMT surface. Our models propose that increasing steric volume of e.g. the diethylamine tail of entacapone is favourable for COMT inhibitory activity.

A set of 113 flexible cyclic urea inhibitors of human immunodeficiency virus protease (HIV-1 PR) was used to compare the quality and predictive power of CoMFA and CoMSIA models for manually or automatically aligned inhibitor set. Inhibitors that were aligned automatically with molecular docking were in agreement with information obtained from existing X-ray structures. Both alignment methods produced statistically significant CoMFA and CoMSIA models, with the best q(2) value being 0.649 and the best predictive r(2) being 0.754. The manual alignment gave statistically higher values, whereas the automated alignment gave more robust models for predicting the activities of an external inhibitor set. Both models utilized similar amino acids in the HIV-1 PR active site, supporting the idea that hydrogen bonds form between an inhibitor and the backbone carbonyl oxygens of Gly48 and Gly48' and also the backbone NH group of Asp30, Gly48, Asp29', and Gly48' of the enzyme. These results suggest that an automated inhibitor alignment can yield predictive 3D QSAR models that are well comparable to manual methods. Thus, an automated alignment method in creating 3D QSAR models is encouragable when a well-characterized structure of the target protein is available.

Finding novel lead molecules is one of the primary goals in early phases of drug discovery projects. However, structurally dissimilar compounds may exhibit similar biological activity, and finding new and structurally diverse lead compounds is difficult for computer algorithms. Molecular energy fields are appropriate for finding structurally novel molecules, but they are demanding to calculate and this limits their usefulness in virtual screening of large chemical databases. In our approach, energy fields are computed only once per superposition and a simple interpolation scheme is devised to allow coarse energy field lattices having fewer grid points to be used without any significant loss of accuracy. The resulting processing speed of about 0.25 s per conformation on a 2.4 GHz Intel Pentium processor allows the method to be used for virtual screening on commonly available desktop machines. Moreover, the results indicate that grid-based superposition methods could be efficiently used for the virtual screening of compound libraries.