We report a novel method, ChooseLD (CHOOse biological information Semi-Empirically on the Ligand Docking), which uses simulated annealing (SA) based on bioinformatics for protein–ligand flexible docking. The fingerprint alignment score (FPAScore) value is used to determine the docking conformation of the ligand. This method includes the matching of chemical descriptors such as fingerprints (FPs) and the root mean square deviation (rmsd) calculation of the coordinates of atoms of the chemical descriptors. Here, the FPAScore optimization for the translation and rotation of a rigid body is performed using the Metropolis Monte Carlo method. Our ChooseLD method will find wide application in the field of biochemistry and medicine to improve the search for new drugs targeting various proteins implicated in diseases.