PubMed Health⌕ Search

Biomedical subjects

S K Kearsley

Publications and source records attributed to S K Kearsley.

10 recordsLinked to original sources

Chemical similarity searches using latent semantic structural indexing (LaSSI) and comparison to TOPOSIM.

Similarity searches based on chemical descriptors have proven extremely useful in aiding large-scale drug screening. Here we present results of similarity searching using Latent Semantic Structure Indexing (LaSSI). LaSSI uses a singular value decomposition on chemical descriptors to project molecules into a k-dimensional descriptor space, where k is the number of retained singular values. The effect of the projection is that certain descriptors are emphasized over others and some descriptors may count as partially equivalent to others. We compare LaSSI searches to searches done with TOPOSIM, our standard in-house method, which uses the Dice similarity definition. Standard descriptor-based methods such as TOPOSIM count all descriptors equally and treat all descriptors as independent. For this work we use atom pairs and topological torsions as examples of chemical descriptors. Using objective criteria to determine how effective one similarity method is versus another in selecting active compounds from a large database, we find for a series of 16 drug-like probes that LaSSI is as good as or better than TOPOSIM in selecting active compounds from the MDDR database, if the user is allowed to treat k as an adjustable parameter. Typically, LaSSI selects very different sets of actives than does TOPOSIM, so it can find classes of actives that TOPOSIM would miss.

Algorithms↗

Latent semantic structure indexing (LaSSI) for defining chemical similarity.

A novel method for computing chemical similarity from chemical substructure descriptors is described. This new method, called LaSSI, uses the singular value decomposition (SVD) of a chemical descriptor-molecule matrix to create a low-dimensional representation of the original descriptor space. Ranking molecules by similarity to a probe molecule in the reduced-dimensional space has several advantages over analogous ranking in the original descriptor space: matching latent structures is more robust than matching discrete descriptors, choosing the number of singular values provides a rational way to vary the "fuzziness" of the search, and the reduction in the dimensionality of the chemical space increases searching speed. LaSSI also allows the calculation of the similarity between two descriptors and between a descriptor and a molecule.

Algorithms↗

SQ: a program for rapidly producing pharmacophorically relevent molecular superpositions.

A new method SQ has been developed to provide fast, automatic, and objective pairwise three-dimensional molecular alignments. SQ uses an atom-based clique-matching step followed by an alignment scoring function that has been parametrized to recognize pharmacologically relevant atomic properties. Molecular alignments from SQ are consistent with known drug-receptor interactions. We demonstrate this with six pairs of receptor-ligand complexes from the Brookhaven Protein Data Bank. The SQ-generated alignment of one isolated ligand onto another is shown to approximate the alignment of the ligands when the receptors are superimposed. SQ appears to be better than its predecessor SEAL (Kearsley and Smith, Tetrahedron Comput. Methodol. 1990, 3, 615-633) in this regard. SQ has been tailored so that, given one molecule as a probe, it can be used to routinely scan large chemical databases for which precomputed conformations have been stored. The SQ score, a measure of 3D similarity of each candidate molecule to the probe, can be used to rank compounds for the purposes of chemical screening. We demonstrate this with three probes (a thrombin inhibitor, an HIV protease inhibitor, and a model for angiotensin II). In each case SQ can preferentially select from the MDDR database other compounds with the same activity as the probe. We further show, using the angiotensin example, how SQ can identify topologically diverse compounds with the same activity.

Algorithms↗

FLOG: a system to select 'quasi-flexible' ligands complementary to a receptor of known three-dimensional structure.

We present a system, FLOG (Flexible Ligands Oriented on Grid), that searches a database of 3D coordinates to find molecules complementary to a macromolecular receptor of known 3D structure. The philosophy of FLOG is similar to that reported for DOCK [Shoichet, B.K. et al., J. Comput. Chem., 13 (1992) 380]. In common with that system, we use a match center representation of the volume of the binding cavity and we use a clique-finding algorithm to generate trial orientations of each candidate ligand in the binding site. Also we use a grid representation of the receptor to assess the fit of each orientation. We have introduced a number of novel features within this paradigm. First, we address ligand flexibility by including up to 25 explicit conformations of each structure in our databases. Nonhydrogen atoms in each database entry are assigned one of seven atom types (anion, cation, donor, acceptor, polar, hydrophobic and other) based on their local bonded chemical environments. Second, we have devised a new grid-based scoring function compatible with this 'heavy atom' representation of the ligands. This includes several potentials (electrostatic, hydrogen bonding, hydrophobic and van der Waals) calculated from the location of the receptor atoms. Third, we have improved the fitting stage of the search. Initial dockings are generated with a more efficient clique-finding algorithm. This new algorithm includes the concept of 'essential points', match centers that must be paired with a ligand atom. Also, we introduce the use of a rapid simplex-based rigid-body optimizer to refine the orientations. We demonstrate, using dihydrofolate reductase as a sample receptor, that the FLOG system can select known inhibitors from a large database of drug-like compounds.

Algorithms↗

Flexibases: a way to enhance the use of molecular docking methods.

Specially expanded databases containing three-dimensional structures are created to enhance the utility of docking methods to find new leads, i.e., active compounds of pharmacological interest. The expansion is based on the automatic generation of a set of maximally dissimilar conformations. The ligand receptor system of methotrexate and dihydrofolate reductase is used to demonstrate the feasibility of creating flexibases and their utility in docking studies.

Algorithms↗

Three-dimensional structure of acyl carrier protein determined by NMR pseudoenergy and distance geometry calculations.

Distance constraints from two-dimensional NMR cross-relaxation data are used to derive a three-dimensional structure for acyl carrier protein from Escherichia coli. Several approaches to structure determination are explored. The most successful proves to be an approach that combines the early stages of a distance geometry program with energy minimization in the presence of NMR constraints represented as pseudopotentials. Approximately 450 proton to proton distance constraints including 50 long-range constraints were included in these programs. Starting structures were generated at random by the distance geometry program and energies minimized by a molecular mechanics module to give final structures. Seven of the structures were deemed acceptable on the basis of agreement with experimentally determined distances. Root-mean-square deviations from the mean of these structures for backbone atoms range from 2 to 3 A. All structures show three roughly parallel helices with hydrophobic residues facing inward and hydrophilic residues facing outward. A hydrophobic cleft is recognizable and is identified as a likely site for acyl chain binding.

Acyl Carrier Protein↗

Designing targeted libraries with genetic algorithms.

In combinatorial synthesis, molecules are assembled by linking chemically similar fragments. Because the number of available chemical fragments often greatly exceeds the number that can be used in one synthetic experiment, one needs a rational method for choosing a subset of desirable fragments. If a combinatorial library is to be targeted against a particular biological activity, virtual screening methods can be used to predict which molecules in a virtual library are most likely to be active. When the number of possible molecules in a virtual library is very large, genetic algorithms (GAs) or simulated annealing can be used to quickly find high-scoring molecules by sampling a small subset of the total combinatorial space. We previously demonstrated how a GA can be used to select a subset of fragments for a combinatorial library, and we used topology-based methods of scoring. Here we extend that earlier work in three ways. (1) We demonstrate use of the GA with 3D scoring methods developed in our laboratory. (2) We show that the approach of assembling libraries from fragments in high-scoring molecules is a reasonable one. (3) We compare results from a library-based GA to those from a molecule-based GA.

Algorithms↗

Comparison of knowledge-based and distance geometry approaches for generation of molecular conformations.

A knowledge-based approach for generating conformations of molecules has been developed. The method described here provides a good sampling of the molecule's conformational space by restricting the generated conformations to those consistent with the reference database. The present approach, internally named et for enumerate torsions, differs from previous database-mining approaches by employing a library of much larger substructures while treating open chains, rings, and combinations of chains and rings in the same manner. In addition to knowledge in the form of observed torsion angles, some knowledge from the medicinal chemist is captured in the form of which substructures are identified. The knowledge-based approach is compared to Blaney et al.'s distance geometry (DG) algorithm for sampling the conformational space of molecules. The structures of 113 protein-bound molecules, determined by X-ray crystallography, were used to compare the methods. The present knowledge-based approach (i) generates conformations closer to the experimentally determined conformation, (ii) generates them sooner, and (iii) is significantly faster than the DG method.

Algorithms↗

Protocols for bridging the peptide to nonpeptide gap in topological similarity searches.

Similarity searches based on chemical descriptors have proven extremely useful in aiding large-scale drug screening. Typically an investigator starts with a "probe", a drug-like molecule with an interesting biological activity, and searches a database to find similar compounds. In some projects, however, the only known actives are peptides, and the investigator needs to identify drug-like actives. 3D similarity methods are able to help in this endeavor but suffer from the necessity of having to specify the active conformation of the probe, something that is not always possible at the beginning of a project. Also, 3D methods are slow and are complicated by the need to generate low-energy conformations. In contrast, topological methods are relatively rapid and do not depend on conformation. However, unmodified topological similarity methods, given a peptide probe, will preferentially select other peptides from a database. In this paper we show some simple protocols that, if used with a standard topological similarity search method, are sufficient to select nonpeptide actives given a peptide probe. We demonstrate these protocols by using 10 peptide-like probes to select appropriate nonpeptide actives from the MDDR database.

Computer Simulation↗