Dear Sebastian,
one possibility might be to ask the Vina extension to compute many modes (e.g. 1000) and then use a script to filter out the conformations that do not satisfy the distance constraint. Another possibility might be to directly modify the source code of Vina to add a term to the scoring function that would penalize one or more unsatisfied constraint(s) (similar to a NMR energy function). We're going to look into this. Thanks for the suggestion.
Best,
Stephane