The accurate prediction of protein–ligand binding affinity belongs to one of the central goals in computer-based drug design. Molecular dynamics (MD)-based free energy calculations have become increasingly popular in this respect due to their accuracy and solid theoretical basis. Here, we present a combined study which encompasses experimental and computational studies on two series of factor Xa ligands, which enclose a broad chemical space including large modifications of the central scaffold. Using this integrated approach, we identified several new ligands with different heterocyclic scaffolds different from the previously identified indole-2-carboxamides that show superior or similar affinity. Furthermore, the so far underexplored terminal alkyne moiety proved to be a suitable non-classical bioisosteric replacement for the higher halogen−π aryl interactions. With this challenging example, we demonstrated the ability of the MD-based non-equilibrium free energy calculation approach for guiding crucial modifications in the lead optimization process, such as scaffold replacement and single-site modifications at molecular interaction hot spots.