FLAME: A Program to Flexibly Align Molecules

Structure-based design of novel inhibitors of the MDM2-p53 interaction

Structure-based rational design led to the discovery of novel inhibitors of the MDM2-p53 protein-protein interaction. The affinity of these compounds for MDM2 was improved through conformational control of both the piperidinone ring and the appended N-alkyl substituent. Optimization afforded 29 (AM-8553), a potent and selective MDM2 inhibitor with excellent pharmacokinetic properties and in vivo efficacy.

Development and validation of an improved algorithm for overlaying flexible molecules

A program for overlaying multiple flexible molecules has been developed. Candidate overlays are generated by a novel fingerprint algorithm, scored on three objective functions (union volume, hydrogen-bond match, and hydrophobic match), and ranked by constrained Pareto ranking. A diverse subset of the best ranked solutions is chosen using an overlay-dissimilarity metric. If necessary, the solutions can be optimised. A multi-objective genetic algorithm can be used to find additional overlays with a given mapping of chemical features but different ligand conformations. The fingerprint algorithm may also be used to produce constrained overlays, in which user-specified chemical groups are forced to be superimposed. The program has been tested on several sets of ligands, for each of which the true overlay is known from protein–ligand crystal structures. Both objective and subjective success criteria indicate that good results are obtained on the majority of these sets.

Graph-based similarity concepts in virtual screening

Applying similarity for finding new promising compounds is a key issue in drug design. Conversely, quantifying similarity between molecules has remained a difficult task despite the numerous approaches. Here, some general aspects along with recent developments regarding similarity criteria are collected. For the purpose of virtual screening, the compounds have to be encoded into a computer-readable format that permits a comparison, according to given similarity criteria, comprising the use of the 3D structure, fingerprints, graph-based and alignment-based approaches. Whereas finding the most common substructures is the most obvious method, more recent approaches take into account chemical modifications that appear throughout existing drugs, from various therapeutic categories and targets.

The flexible alignment of molecular structures using simulated annealing with weighted Lagrangian multipliers

A framework for superimposing small molecules is presented. The proposed method consists of a simple atom-based, flexible alignment. The optimization procedure used in the alignment is based on a recently published variant of the simulated annealing whereby nonlinear constraints are accommodated using Lagrangian multipliers. It differs from other published superposition algorithms in that any number of nonlinear constraints can be readily imposed on the structural alignment directly through the objective function without assuming an a priori trade-off between competing conditions. These can include equality and equality constraints on distances, angles, and energy states. Examples illustrating the use of the proposed approach are also provided.

pharmACOphore: multiple flexible ligand alignment based on ant colony optimization

The flexible superimposition of biologically active ligands is a crucial step in ligand-based drug design. Here we present pharmACOphore, a new approach for pairwise as well as multiple flexible alignment of ligands based on ant colony optimization (ACO; Dorigo, M.; Stützle, T. Ant Colony Optimization; MIT Press: Cambridge, MA, USA, 2004). An empirical scoring function is used, which describes ligand similarity by minimizing the distance of pharmacophoric features. The scoring function was parametrized on pairwise alignments of ligand sets for four proteins from diverse protein families (cyclooxygenase-2, cyclin-dependent kinase 2, factor Xa and peroxisome proliferator-activated receptor γ). The derived parameters were assessed with respect to pose prediction performance on the independent FlexS data set ( Lemmen, C.; Lengauer, T.; Klebe, G. J. Med. Chem. 1998, 41, 4502 - 4520) in exhausting pairwise alignments. Additionally, multiple flexible alignment experiments were carried out for the pharmacologically relevant targets trypsin and poly (ADP-ribose) polymerase (PARP). The results obtained show that the new procedure provides a robust and efficient way for the pairwise as well as multiple flexible alignment of small molecules.

Training a scoring function for the alignment of small molecules

A comprehensive data set of aligned ligands with highly similar binding pockets from the Protein Data Bank has been built. Based on this data set, a scoring function for recognizing good alignment poses for small molecules has been developed. This function is based on atoms and hydrogen-bond projected features. The concept is simply that atoms and features of a similar type (hydrogen-bond acceptors/donors and hydrophobic) tend to occupy the same space in a binding pocket and atoms of incompatible types often tend to avoid the same space. Comparison with some recently published results of small molecule alignments shows that the current scoring function can lead to performance better than those of several existing methods.

Technique for generating three-dimensional alignments of multiple ligands from one-dimensional alignments

We describe and demonstrate a method for the simultaneous, fully flexible alignment of multiple molecules with a common biological activity. The key aspect of the algorithm is that the alignment problem is first solved in a lower dimensional space, in this case using the one-dimensional representations of the molecules. The three-dimensional alignment is then guided by constraints derived from the one-dimensional alignment. We demonstrate using 10 hERG channel blockers, with a total of 72 rotatable bonds, that the one-dimensional alignment is able to effectively isolate key conserved pharmacophoric features and that these conserved features can effectively guide the three-dimensional alignment. Further using 10 estrogen receptor agonists and 5 estrogen receptor antagonists with publicly available cocrystal structures we show that the method is able to produce superpositions comparable to those derived from crystal structures. Finally, we demonstrate, using examples from peptidic CXCR3 agonists, that the method is able to generate reasonable binding hypotheses.