First, molecules were docked with constraints to the core pattern of the Vioxx bound ligand, where the RMSD between the maximum common substructure between the ligand and Vioxx was limited to 2.5 Å. Then, the constraints were lifted, and the molecules docked again. The two sets of results were combined, and the pose with a lower GlideSP Score for each molecule was selected as the final pose.