In drug design, it is crucial to have reliable information on how a chemical entity behaves in the presence of metabolizing enzymes. This requires substantial experimental efforts. Consequently, being able to predict the likely site/s of metabolism in any compound, synthesized or virtual, would be highly beneficial and time efficient. In this work, six different methodologies for predictions of the site of metabolism (SOM) have been compared and validated using structurally diverse data sets of drug-like molecules with well-established metabolic pattern in CYP3A4, CYP2C9, or both. Three of the methods predict the SOM based on the ligand's chemical structure, two additional methods use structural information of the enzymes, and the sixth method combines structure and ligand similarity and reactivity. The SOM is correctly predicted in 50 to 90% of the cases, depending on method and enzyme, which is an encouraging rate. We also discuss the underlying mechanisms of cytochrome P450 metabolism in the light of the results from this comparison.

Molecular docking plays an important role in drug discovery as a tool for the structure-based design of small organic ligands for macromolecules. Possible applications of docking are identification of the bioactive conformation of a protein-ligand complex and the ranking of different ligands with respect to their strength of binding to a particular target. We have investigated the effect of implicit water on the postprocessing of binding poses generated by molecular docking using MM-PB/GB-SA (molecular mechanics Poisson-Boltzmann and generalized Born surface area) methodology. The investigation was divided into three parts: geometry optimization, pose selection, and estimation of the relative binding energies of docked protein-ligand complexes. Appropriate geometry optimization afforded more accurate binding poses for 20% of the complexes investigated. The time required for this step was greatly reduced by minimizing the energy of the binding site using GB solvation models rather than minimizing the entire complex using the PB model. By optimizing the geometries of docking poses using the GB(HCT+SA) model then calculating their free energies of binding using the PB implicit solvent model, binding poses similar to those observed in crystal structures were obtained. Rescoring of these poses according to their calculated binding energies resulted in improved correlations with experimental binding data. These correlations could be further improved by applying the postprocessing to several of the most highly ranked poses rather than focusing exclusively on the top-scored pose. The postprocessing protocol was successfully applied to the analysis of a set of Factor Xa inhibitors and a set of glycopeptide ligands for the class II major histocompatibility complex (MHC) A(q) protein. These results indicate that the protocol for the postprocessing of docked protein-ligand complexes developed in this paper may be generally useful for structure-based design in drug discovery.

This work describes a data driven method for scaffold hopping by fragment replacement. A search database of scaffolds is created by cutting bonds of existing compounds in a combinatorial fashion. Three-dimensional structures of the scaffolds are then generated and made searchable based on the relative orientation of the broken bonds using an auxiliary index file. The retrieved scaffolds are ranked using volume overlap and electrostatic similarity scores. A similar approach has been used before in the program CAVEAT and others. The present work introduces a novel indexing scheme for the attachment vector geometry, which allows for fast searching. A scaffold shape descriptor is defined, which allows for queries with a single attachment vector (R-groups) and improves the shape similarity between the query and the suggested replacement fragments. The program, called Scaffold Hopping, is shown to retrieve relevant bioisosteric replacement scaffolds for a set of example queries in a reasonable time frame, making the program suitable to be used in drug design work.

The vinylic polymerization of 5-norbornene-2-carboxylic esters with different palladium catalysts to polymers of high molecular weights is reported. Computer simulations based on a force field approach show that these polymers display a rigid statistical chain conformation and are thus a further example of rotationally strongly constrained polymers. The polymers are soluble in a variety of solvents despite their rigidity, they are amorphous, possess glass transition temperatures well above 250°C, and have a high packing density. The dipoles located in the lateral groups perform a secondary relaxation process similar to the case of flexible or rigid rod-like polymers containing ester groups. Acta Polymer., 48,385-391 Q VCH Verlagsgesellschaft mbH, D-69451 Weinheim 1997 0323-7648/97/0808-0385$17.50+.50/0 385