In this paper, we present a systematic robust dynamic optimization framework applied to the benzaldehyde lyase-catalyzed carboligation of propanal and benzaldehyde to produce (R)-2-hydroxy-1-phenylbutan-1-one (BA). First, the elementary process functions approach was used to screen between different dosing concepts, and it was found that simultaneously dosing propanal and benzaldehyde leads to the highest final concentration of BA. Next, we applied global sensitivity analysis and found that 10 out of 13 kinetic parameters are relevant. Time-varying back-offs were then used to handle parametric uncertainties due to these 10 parameters. A major contribution in our work is the use of the point estimate method instead of Monte Carlo simulations to calculate the back-offs in an efficient and reproducible manner. We show that this new approach is at least 10 times faster than the conventional Monte Carlo approach while achieving low approximation errors. (C) 2018 Elsevier Ltd. All rights reserved.