Quantitative risk assessment via uncertainty analysis in combination with error propagation for the determination of the dynamic Design Space of the primary drying step during freeze-drying.

Traditional pharmaceutical freeze-drying is an inefficient batch process often applied to improve the stability of biopharmaceutical drug products. The freeze-drying process is regulated by the (dynamic) settings of the adaptable process parameters shelf temperature Ts and chamber pressure Pc. Mechanistic modelling of the primary drying step allows the computation of the optimal combination of Ts and Pc in function of the primary drying time. In this study, an uncertainty analysis was performed on the mechanistic primary drying model to construct the dynamic Design Space for the primary drying step of a freeze-drying process, allowing to quantitatively estimate and control the risk of cake collapse (i.e., the Risk of Failure (RoF)). The propagation of the error on the estimation of the thickness of the dried layer Ldried as function of primary drying time was included in the uncertainty analysis. The constructed dynamic Design Space and the predicted primary drying endpoint were experimentally verified for different RoF acceptance levels (1%, 25%, 50% and 99% RoF), defined as the chance of macroscopic cake collapse in one or more vials. An acceptable cake structure was only obtained for the verification runs with a preset RoF of 1% and 25%. The run with the nominal values for the input variables (RoF of 50%) led to collapse in a significant number of vials. For each RoF acceptance level, the experimentally determined primary drying endpoint was situated below the computed endpoint, with a certainty of 99%, ensuring sublimation was finished before secondary drying was started. The uncertainty on the model input parameters demonstrates the need of the uncertainty analysis for the determination of the dynamic Design Space to quantitatively estimate the risk of batch rejection due to cake collapse.