Structure-based ligand design for flexible proteins: application of new F-DycoBlock

A method of structure-based ligand design - DycoBlock - has been proposed and tested by Liu et al. It was further improved by Zhu et al. and applied to design new selective inhibitors of cyclooxygenase 2. In the current work, we present a new methodology - F-DycoBlock that allows for the incorporation of receptor flexibility. During the designing procedure, both the receptor and molecular building blocks are subjected to the multiple-copy stochastic molecular dynamics (MCSMD) simulation, while the protein moves in the mean field of all copies. It is tested for two enzymes studied previously - cyclooxygenase 2 (COX-2) and human immunodeficiency type 1 (HIV-1) protease. To identify the applicability of F-DycoBlock, the binding protein structure was used as starting point to explore the conformational space around the bound state. This method can be easily extended to accommodate the flexibility in different degree. Four types of treatment of the receptor flexibility - all-atom restrained, backbone restrained, intramolecular hydrogen-bond restrained and active-site flexible - were tested with or without the grid approximation. Two inhibitors, SC-558 for COX-2 and L700417 for HIV-1 protease, are used in this testing study for comparison with previous results. The accuracy of recovery, binding energy, solvent accessible surface area (SASA) and positional root-mean-square (RMS) deviation are used as criteria. The results indicate that F-DycoBlock is a robust methodology for flexible drug design. It is particularly notable that the protein flexibility has been perfectly associated with each stage of drug design - search for the binding sites, dynamic assembly and optimization of candidate compounds. When all protein atoms were restrained, F-DycoBlock yielded higher accuracy of recovery than DycoBlock (100%). If backbone atoms were restrained, the same ratio of accuracy was achieved. Moreover, with the intramolecular hydrogen bonds restrained, reasonable conformational changes were observed for HIV- 1 protease during the long-time MCSMD simulation and L700417 was reassembled at the active site. It makes it possible to study the receptor motion in the binding process.