Automated site-directed drug design: the formation of molecular templates in primary structure generation

In this paper the spacer skeleton concept is used to produce molecular graphs of putative ligands for binding sites. The skeletons are transformed into molecular templates within the constraints of the accessible surface of the ligand-binding site. A distance-matrix method is used to compare ligand points with vertices of the spacer skeleton through a permutation of all possible correspondences. A tolerance parameter is used to screen for poor matches. As a result, a small number of matched vertices and ligand points are produced. These are fitted into the site by a constrained optimization routine using an analytical function. Ligand points fall within the site and are optimally positioned adjacent to the corresponding site points; other vertices of the spacer skeleton lying beneath the accessible surface of the site are clipped off. A molecular template is thereby formed with its vertices linked to the ligand points. The final step is to verify that the bonding integrity of the skeleton remains. The computational methods outlined in this paper have been tested at two binding sites: the pteridine binding site in dihydrofolate reductase and the amidinophenylpyruvate site of trypsin. Molecular graphs for both sites were generated automatically; they showed strong similarity to those of the natural ligands.