Whilst shape-based complementarity between a ligand and the protein active site are critical for high affinity binding there are very few high performance methods. This recent publication offers an interesting methodology.
Riemannian geometry and molecular similarity I: spectrum of the Laplacian [DOI], all the code is on GitHub.
https://github.com/RPirie96/RGMolSA