Designing optimal protocols in Bayesian quantum parameter estimation with higher-order operations