In this study, we developed a novel scoring program, HotLig, which applies the Connolly surface of a protein to calculate hydrophobic interaction and paired pharmacophore interactions with ligands. In addition to molecular surface distance, ligand-contacting areas and hydrogen-bond angles were also introduced to the scoring functions in HotLig