A computational protocol to fragment-based drug design at PDB scale

Fragment- based drug design has emerged in the last decade. This exciting field has been recently reviewed [1]. Obtaining structural information on the fragment complexed to the protein target is a key factor and also a major limitation. Therefore, computational methods are needed to mine efficiently all the available 3D structures of ligands complexed to proteins, both treated as a whole and as smaller fragment to increase the likelihood of fragment hopping from one target to another. MED-SuMo [2] is a target based drug design tool. This tool is based on the identification of local shape and 3D Surface Chemical Features similarities in the target binding site with other proteins (with their co-crystallized ligands). MED-SuMo uses the binding site of the target as a query to search either the whole PDB (or any corporate protein structure databank) for all the binding sites that display a local match with the query. This valuable information can then be used to identify which residues of the binding site are potentially important for affinity and selectivity. As the similar binding sites are overlaid, the co-crystallized ligands are aligned and are therefore a starting material for ligand hybridization. Among the hundreds of overlaid binding sites generated by MED-SuMo, we found the protein-ligand complexes overlaid by Pierce et al. [3] for ligand hybridization in the BREED method. Interestingly, they found relevant combinations of fragments starting from only a few protein-ligand complexes structures and we believe that the output of MED-SuMo is a very promising input for automatic methods like BREED. In this work, we used MED-SuMo on a fragment database derived from the PDB: each pdb file is converted into one or more pdb files containing a single fragment as ligand. MED-SuMo is used to query and mine the Surface Chemical Functions surrounding the fragments, seeking similarities with the protein of interest. A library of fragments positioned in the active site and annotated with the counts of contacts and h-bonds is generated. This is the basis of a protocol which ranks and hybridizes fragments. Fragmentation and hybridization steps are coupled with substructure searches in a database of chemical compounds from suppliers to increase the likelihood of fragments and hybridized molecules as synthetically feasible molecules.