SPLICE: A program to assemble partial query solutions from three-dimensional database searches into novel ligands

Summary SPLICE is a program that processes partial query solutions retrieved from 3D, structural databases to generate novel, aggregate ligands. It is designed to interface with the database searching program FOUNDATION, which retrieves fragments containing any combination of a user-specified minimum number of matching query elements. SPLICE eliminates aspects of structures that are physically incapable of binding within the active site. Then, a systematic rule-based procedure is performed upon the remaining fragments to ensure receptor complementarity. All modifications are automated and remain transparent to the user. Ligands are then assembled by linking components into composite structures through overlapping bonds. As a control experiment, FOUNDATION and SPLICE were used to reconstruct a know HIV-1 protease inhibitor after it had been fragmented, reoriented, and added to a sham database of fifty different small molecules. To illustrate the capabilities of this program, a 3D search query containing the pharmacophoric elements of an aspartic proteinase-inhibitor crystal complex was searched using FOUNDATION against a subset of the Cambridge Structural Database. One hundred thirty-one compounds were retrieved, each containing any combination of at least four query elements. Compounds were automatically screened and edited for receptor complementarity. Numerous combinations of fragments were discovered that could be linked to form novel structures, containing a greater number of pharmacophoric elements than any single retrieved fragment.     

Characterising the geometric diversity of functional groups in chemical databases

Summary We have developed a program, HookSpace, which provides a simplistic approach to assessing the diversity of molecular databases. The spatial relationship between pairs of intramolecular functional groups can be analysed in a variety of ways to provide both qualitative and quantitative measures of diversity. Results are described and contrasted for two commercially available databases and a combinatorial library of benzodiazepam derivatives. HookSpace highlights the main differences in molecular content of these data sets.     

The computer program LUDI: A new method for the de novo design of enzyme inhibitors

Summary A new computer program is described, which positions small molecules into clefts of protein structures (e.g. an active site of an enzyme) in such a way that hydrogen bonds can be formed with the enzyme and hydrophobic pockets are filled with hydrophobic groups. The program works in three steps. First it calculates interaction sites, which are discrete positions in space suitable to form hydrogen bonds or to fill a hydrophobic pocket. The interaction sites are derived from distributions of nonbonded contacts generated by a search through the Cambridge Structural Database. An alternative route to generate the interaction sites is the use of rules. The second step is the fit of molecular fragments onto the interaction sites. Currently we use a library of 600 fragments for the fitting. The final step in the present program is the connection of some or all of the fitted fragments to a single molecule. This is done by bridge fragments. Applications are presented for the crystal packing of benzoic acid and the enzymes dihydrofolate reductase and trypsin.     

PRO_LIGAND: An approach to de novo molecular design. 3. A genetic algorithm for structure refinement

Summary Recently, the development of computer programs which permit the de novo design of molecular structures satisfying a set of steric and chemical constraints has become a burgeoning area of research and many operational systems have been reported in the literature. Experience with PRO_LIGAND—the de novo design methodology embodied in our in-house molecular design and simulation system PRO-METHEUS—has suggested that the addition of a genetic algorithm (GA) structure refinement procedure can add value to an already useful tool. Starting with the set of designed molecules as an initial population, the GA can combine features from both high- and low-scoring structures and, over a number of generations, produce individuals of better score than any of the starting structures. This paper describes how we have implemented such a procedure and demonstrates its efficacy in improving two sets of molecules generated by different de novo design projects.     

CAVEAT: A program to facilitate the design of organic molecules

Summary A frequently encountered problem in the design of enzyme inhibitors and other biologically active molecules is the identification of molecular frameworks to serve as templates or linking units that can position functional groups in specific relative orientations. The program CAVEAT was designed to address this problem by searching 3D databases for such molecular fragments. Key innovations introduced in CAVEAT are a focus on relationships between bonds and the provision of automated methods to identify and classify structural frameworks. Performance has been a particular concern in formulating CAVEAT, since it is intended to be used in an interactive manner. The focus in this report is the design and implementation of the principal algorthms and the performance achieved.CAVEAT is available from the Office of Technology Licensing, University of California, Berkeley, CA, and from Molecular Simulations, Inc., Burlington, MA; further information is available from P.A. Bartlett.     

LUDI: rule-based automatic design of new substituents for enzyme inhibitor leads

Summary Recent advances in a new method for the de novo design of enzyme inhibitors are reported. A new set of rules to define the possible nonbonded contacts between protein and ligand is presented. This method was derived from published statistical analyses of nonbonded contacts in crystal packings of organic molecules and has been implemented in the recently described computer program LUDI. Moreover, LUDI can now append a new substituent onto an already existing ligand. Applications are reported for the design of inhibitors of HIV protease and dihydrofolate reductase. The results demonstrate that LUDI is indeed capable of designing new ligands with improved binding when compared to the reference compound.     

Automated site-directed drug design: The generation of a basic set of fragments to be used for automated structure assembly

Summary If a method is to be developed to assemble putative ligands structures in site-directed drug design, from molecular graphs generated in the site, then basic building blocks are needed. Structure assembly is a combinatoric process that needs to be optimised if it is to be tractable. What has to be determined is whether small molecular fragments can have transferable properties from one molecule to another. In this paper we determine all possible combinations of 3-, 4- and 5-atom aliphatic fragments from a small set of atoms H, C, N, O, F or Cl. The frequency of occurrence of these candidate fragments is searched for in the Cambridge Structural Database. A similar analysis is performed on charged fragments. A more restricted search is carried out for P and S and aromatic structures. A basic set of fragments can be derived that have a significant frequency in known crystal structures. The transferability of fragment properties is discussed in subsequent papers.     

Similarity searching in files of three-dimensional chemical structures: Representation and searching of molecular electrostatic potentials using field-graphs

This paper reports a method for the identification of those molecules in a database of rigid 3D structures with molecular electrostatic potential (MEP) grids that are most similar to that of a user-defined target molecule. The most important features of an MEP grid are encoded in field-graphs, and a target molecule is matched against a database molecule by a comparison of the corresponding field-graphs. The matching is effected using a maximal common subgraph isomorphism algorithm, which provides an alignment of the target molecule's field- graph with those of each of the database molecules in turn. These alignments are used in the second stage of the search algorithm to calculate the intermolecular MEP similarities. Several different ways of generating field-graphs are evaluated, in terms of the effectiveness of the resulting similarity measures and of the associated computational costs. The most appropriate procedure has been implemented in an operational system that searches a corporate database, containing ca. 173,000 3D structures.     

The development of a simple empirical scoring function to estimate the binding constant for a protein-ligand complex of known three-dimensional structure

Summary A new simple empirical function has been developed that estimates the free energy of binding for a given protein-ligand complex of known 3D structure. The function takes into account hydrogen bonds, ionic interactions, the lipophilic protein-ligand contact surface and the number of rotatable bonds in the ligand. The dataset for the calibration of the function consists of 45 protein-ligand complexes. The new energy function reproduces the binding constants (ranging from 2.5·10^-2 to 4·10^-14 M, corresponding to binding energies between -9 and -76 kJ/mol) of the dataset with a standard deviation of 7.9 kJ/mol, corresponding to 1.4 orders of magnitude in binding affinity. The individual contributions to protein-ligand binding obtained from the scoring function are: ideal neutral hydrogen bond: -4.7 kJ/mol; ideal ionic interaction: -8.3 kJ/mol; lipophilic contact: -0.17 kJ/mol Ų; one rotatable bond in the ligand: +1.4 kJ/mol. The function also contains a constant contribution (+5.4 kJ/mol) which may be rationalized as loss of translational and rotational entropy. The function can be evaluated very fast and is therefore also suitable for application in a 3D database search or de novo ligand design program such as LUDI.     

A shape- and chemistry-based docking method and its use in the design of HIV-1 protease inhibitors

Summary The program DOCK [1,2] has been used successfully to identify molecules which will bind to a specified receptor [3]. The original method ranks molecules based on their shape complementarity to the receptor site and relies on the chemist to bring the appropriate electrostatic or hydrogen bond properties into the molecular skeletons obtained in the search. This is useful when screening a small database of compounds, where it is not likely that molecules with both the correct shape and electrostatic properties will be found. As large databases are more likely to have redundant molecular shapes with a variety of functionality (e.g., members of a congeneric series), it would be useful to have a method which identifies molecules with both the correct shape and functionality. To this end we have modified the DOCK 1.0 method to target user-specified atom types to selected positions in the receptor site. The target sites can be chosen based on structural evidence, calculation or inspection. Targeted-DOCK improves the ability of the DOCK method to find the crystallographically determined binding mode of a ligand. Additionally, targeted-DOCK searches a database of small molecules at 100–1000 times the rate of DOCK 1.0, allowing more molecules to be screened and more sophisticated scoring schemes to be employed. Targeted-DOCK has been used successfully in the design of a novel non-peptide inhibitor of HIV-1 protease.     

On the use of LUDI to search the Fine Chemicals Directory for ligands of proteins of known three-dimensional structure

Summary It is shown that the computer program LUDI can be used to search large databases of three-dimensional structures for putative ligands of proteins with known 3D structure. As an example, a subset of 30 000 small molecules (with less than 40 atoms and 0–2 rotatable bonds) from the Fine Chemicals Directory has been used in the search for possible novel ligands for four different proteins (trypsin, streptavidin, purine nucleoside phosphorylase and HIV protease). For trypsin and streptavidin, known ligands or substructures of known ligands are retrieved as top-scoring hits. In addition, a number of new interesting structures are found in all considered cases. Therefore, the method holds promise to retrieve automatically protein ligands from a 3D database if the 3D structure of the target protein is known.     

A genetic algorithm for the automated generation of molecules within constraints

Summary A genetic algorithm has been designed which generates molecular structures within constraints. The constraints may be any useful function, for example an enzyme active site, a pharmacophore or molecular properties from pattern recognition or rule-induction analyses. The starting point may be random or may utilise known molecules. These are modified to grow into families of structures which, using the evolutionary operators of selection, crossover and mutation evolve to better fit the constraints. The basis of the algorithm is described together with some applications in lead generation, 3D database construction and drug design. Genetic algorithms of this type may have wider applications in chemistry, for example in the design and optimisation of new polymers, materials (e.g. superconducting materials) or synthetic enzymes.     

Automated site-directed drug design: Searches of the Cambridge Structural Database for bond lengths in molecular fragments to be used for automated structure assembly

Summary In this paper a database of small frequently occurring molecular fragments is used for the determination of fragment bond lengths from the Cambridge Structural Database. A large number of bond types are described that have not been reported previously.     

PRO_LIGAND: An approach to de novo molecular design. 6. Flexible fitting in the design of peptides

Summary This paper describes the further development of the functionality of our in-house de novo design program, PRO_LIGAND. In particular, attention is focussed on the implementation and validation of the directed tweak method for the construction of conformationally flexible molecules, such as peptides, from molecular fragments. This flexible fitting method is compared to the original method based on libraries of prestored conformations for each fragment. It is shown that the directed tweak method produces results of comparable quality, with significant time savings. By removing the need to generate a set of representative conformers for any new library fragment, the flexible fitting method increases the speed and simplicity with which new fragments can be included in a fragment library and also reduces the disk space required for library storage. A further improvement to the molecular construction process within PRO_LIGAND is the inclusion of a constrained minimisation procedure which relaxes fragments onto the design model and can be used to reject highly strained structures during the structure generation phase. This relaxation is shown to be very useful in simple test cases, but restricts diversity for more realistic examples. The advantages and disadvantages of these additions to the PRO_LIGAND methodology are illustrated by three examples: similar design to an alpha helix region of dihydrofolate reductase, complementary design to the active site of HIV-1 protease and similar design to an epitope region of lysozyme.     

Empirical scoring functions: I. The development of a fast empirical scoring function to estimate the binding affinity of ligands in receptor complexes

This paper describes the development of a simple empirical scoringfunction designed to estimate the free energy of binding for aprotein–ligand complex when the 3D structure of the complex is knownor can be approximated. The function uses simple contact terms to estimatelipophilic and metal–ligand binding contributions, a simple explicitform for hydrogen bonds and a term which penalises flexibility. Thecoefficients of each term are obtained using a regression based on 82ligand–receptor complexes for which the binding affinity is known. Thefunction reproduces the binding affinity of the complexes with across-validated error of 8.68 kJ/mol. Tests on internal consistency indicatethat the coefficients obtained are stable to changes in the composition ofthe training set. The function is also tested on two test sets containing afurther 20 and 10 complexes, respectively. The deficiencies of this type offunction are discussed and it is compared to approaches by other workers.     

Prediction of binding constants of protein ligands: A fast method for the prioritization of hits obtained from de novo design or 3D database search programs

A dataset of 82 protein–ligand complexes of known 3D structure and binding constant K_i was analysed to elucidate the important factors that determine the strength of protein–ligand interactions. The following parameters were investigated: the number and geometry of hydrogen bonds and ionic interactions between the protein and the ligand, the size of the lipophilic contact surface, the flexibility of the ligand, the electrostatic potential in the binding site, water molecules in the binding site, cavities along the protein–ligand interface and specific interactions between aromatic rings. Based on these parameters, a new empirical scoring function is presented that estimates the free energy of binding for a protein–ligand complex of known 3D structure. The function distinguishes between buried and solvent accessible hydrogen bonds. It tolerates deviations in the hydrogen bond geometry of up to 0.25 Å in the length and up to 30 °Cs in the hydrogen bond angle without penalizing the score. The new energy function reproduces the binding constants (ranging from 3.7 × 10^-2 M to 1 × 10^-14 M, corresponding to binding energies between -8 and -80 kJ/mol) of the dataset with a standard deviation of 7.3 kJ/mol corresponding to 1.3 orders of magnitude in binding affinity. The function can be evaluated very fast and is therefore also suitable for the application in a 3D database search or de novo ligand design program such as LUDI. The physical significance of the individual contributions is discussed.     

Peptide mimetics as enzyme inhibitors: Use of free energy perturbation calculations to evaluate isosteric replacement for amide bonds in a potent HIV protease inhibitor

Summary We present the application of free energy perturbation theory/molecular dynamics to predict the consequence of replacing each of the seven peptide bonds in the potent HIV protease inhibitor JG365: ACE (acetyl)-Ser-Leu-Asn-HEA (hydroxyethylamine analog of Phe-Pro)-Ile-Val-NME (N-methyl) by ethylene or fluoroethylene isosteres. Replacing two of these bonds may well lead to significantly tighter binding; replacing two others is predicted to significantly diminish the binding affinity. Also, for three of the peptide bonds fluoroethylene replacements could lead to increased binding of free energies of the inhibitors. Our results should be considered as predictive since there are, as yet, no experimental results on such peptide replacements as enzyme inhibitors.     

Naphthalimides and analogues as antitumor agents: A review on molecular design,bioactivity and mechanism of action

Our research improves the structure diversity of naphthalimide antitumor agents and distinct variances of antitumor targets and mechanism of action.