TransferFunctionModel[{{{k (1 + s)}}, ((-1 + s) s) (16 + 4 s + s^2)}, s]

tf[f_, s_] := TransferFunctionModel[f, s];

tf[{{{k (1 + s)}}, ((-1 + s) s) (16 + 4 s + s^2)}, s]