This might help. If we assume z>0
then we get some simplification or at least a convenient rewriting in terms of lower level functions.
Simplify[
FunctionExpand[HypergeometricU[1, 1, z*k],
Assumptions -> {0 < k < 1, z > 0}],
Assumptions -> {0 < k < 1, z > 0}]
(* Out[1204]= -E^(k z) ExpIntegralEi[-k z] *)
We integrate that part.
Integrate[
E^(k z) ExpIntegralEi[-k z] k^(a - 1) (1 - k)^(-a), {k, 0, 1},
Assumptions -> z > 0]
(* Out[1207]= ConditionalExpression[-\[Pi] Csc[a \[Pi]] Gamma[
a] HypergeometricU[a, 1, z], 0 < Re[a] < 1] *)
It might still be difficult to verify this but it's probably easier than working with the original integrand.