I like to use ludi, I read the paper Automated denovo design workflow with physics-based scoring function for fast lead identification and optimization.
But i can't find the protocol