Construction of a three-dimensional pharmacophore for Bcl-2 inhibitors by flexible docking and the multiple copy simultaneous search method

B-Cell lymphoma-2 (Bcl-2) protein is a new promising target for anticancer drugs. A number of anticancer Bcl-2 inhibitors with diverse chemical structures have been discovered in recent years. In this paper, the flexible docking was performed to determine the binding modes of the representative inhibitors from different structural types. Subsequently, the binding modes of inhibitor were used to construct a primary three- dimensional (3D) pharmacophore model. It proved that this model can effectively disrupt the binding of the BH3 domain of proapoptotic Bcl-2 family members to Bcl-2, and match the structural requirement of a new type of Bcl-2 inhibitors. However, these distances between pharmacophoric points are not optimal due to the fact that not all of individual functional groups are located in the ideal position when inhibitors bind to its receptor. In this paper, we proposed a new idea to improve the quality of the pharmacophore model: the multiple copy simultaneous search (MCSS) method was performed to determine the energetically favorable distribution of functional groups with similar features to these pharmacophoric points in the active site of Bcl-2 first. Then their most energetically favorable minima in the positions near the pharmacophoric points were used to optimize the distances between pharmacophoric points. By examining the binding modes of several inhibitors from the same structural type, it was found that the more potent the inhibitor was, the closer it was to the optimized distances between pharmacophoric points. The optimized 3D pharmacophore model obtained in this paper may provide a good starting point for further rational design of Bcl-2 inhibitors.