Hi, good question. So I think in this case because the operators inside and outside the exponential are all commuting, you can just write the exponential of the sum as a product of exponentials and then just locally exponentiate each Sz_k operator (this can just be done on paper even, ahead of time). So then to measure such a string operator which is now a product of local operators, you can either make an MPO which just has a bond dimension of 1 (product MPO) or just apply each operator to one of the MPS in the expectation value and finally compute the overlap with the other MPS.
Hope that helps you find the answer you need - if not please comment below.
Miles