simplifying indefinite sums containing Kronecker Deltas

Posting as an answer at the request of the OP

There is probably no built-in way to achieve what you want because the result you are after is in general not correct: in order to get it you need to rewrite and reorder the summations, which is only possible under special convergence conditions.

That being said, if you are sure beforehand that all your manipulations will be legal you can automatically simplify expressions like the ones you present using custom symbols and TagSetDelayed. For example

MyDelta /: f_[x_ + y_] MyDelta[x_ + z_] := f[y - z]

will make an expression like f[q]f[k+k1+q]MyDelta[k1+q] transform to f[q]f[k], and something like

MySum /: MySum[expr_ KroneckerDelta[x_ + y_], left___, x_, right___] := MySum[expr /. x -> -y, left, right]

will take expressions like MySum[f[i + j] KroneckerDelta[j + k], i, j, k] into MySum[f[i - k], i, k].

Depending on your intended application, you may want to add stuff (e.g. make MySum linear, regularize MyDelta[0] under the summation sign, etc.), or make some of these replacements manually instead of automatically... but you should always keep in mind that this symbolic transformations can go wildly wrong if you are not very careful about how you implement things, and to which problems you apply them...