Freeze-Drying of Pharmaceutical Proteins in Vials: Modeling of Freezing and Sublimation Steps

The commercial finite element code FEMLAB was used to perform two-dimensional axisymmetric simulations of the temperature profiles and the moving front velocities of standard BSA (bovine serum albumin)-based formulations used to stabilize pharmaceutical proteins during the freeze-drying process. The simulations were validated with both experimental and numerical approaches. In an initial step, the heat transfer phenomena taking place during the cooling of liquid solutions was studied in commercial size glass vials without freezing or sublimation. Then, this model was extended and validated for the freezing process of aqueous BSA-based solutions encountered in the industrial freeze-drying processes with the same vials in order to confirm the identified values of the different thermal conductances between the product and the shelf and between the product and the surroundings. Finally, the conductances between the product and the shelf and between the vial and the surroundings thus determined were used in a dynamic sublimation model with two zones and a moving sublimation front similar to the ones previously proposed in the literature. The simulations showed a satisfactory agreement between experimental and simulated data. The results of this study demonstrated that the freeze-drying process of pharmaceutical proteins in glass vials for standard industrial operating conditions was mainly controlled by the heat transfer from the shelf and the surroundings to the product sublimation front.

[1]  A. I. Liapis,et al.  An analysis of the lyophilization process using a sorption‐sublimation model and various operational policies , 1985 .

[2]  J. Andrieu,et al.  Determination of mass and heat transfer parameters during freeze-drying cycles of pharmaceutical products. , 2005, PDA journal of pharmaceutical science and technology.

[3]  Marco P. Schoen,et al.  A simulation model for the primary drying phase of the freeze-drying cycle , 1995 .

[4]  A. I. Liapis,et al.  Optimal control of a freeze dryer—I Theoretical development and quasi steady state analysis , 1979 .

[5]  A. I. Liapis,et al.  OPTIMAL CONTROL OF THE PRIMARY AND SECONDARY DRYING STAGES OF BULK SOLUTION FREEZE DRYING IN TRAYS , 1998 .

[6]  A. Rasmuson,et al.  Heat transfer in vial lyophilization. , 2002, International journal of pharmaceutics.

[7]  Herman Vromans,et al.  Vial Lyophilization: Calculations on Rate Limitation During Primary Drying , 2004, Pharmaceutical Research.

[8]  K. Kobe.  The properties of gases and liquids , 1959 .

[9]  A. Liapis,et al.  Advances in the Modeling and Control of Freeze-Drying , 1984 .

[10]  Michael J. Pikal,et al.  The secondary drying stage of freeze drying: drying kinetics as a function of temperature and chamber pressure☆ , 1990 .

[11]  A. I. Liapis,et al.  An adsorption-sublimation model for a freeze dryer , 1979 .

[12]  Steven L Nail,et al.  Fundamentals of freeze-drying. , 2002, Pharmaceutical biotechnology.

[13]  A. I. Liapis,et al.  ESTIMATION OF THE EFFECT OF PRODUCT SHRINKAGE ON THE DRYING TIMES, HEAT INPUT AND CONDENSER LOAD OF THE PRIMARY AND SECONDARY DRYING STAGES OF THE LYOPHILIZATION PROCESS IN VIALS , 1999 .

[14]  Hasan U. Akay,et al.  A computational model for finite element analysis of the freeze-drying process , 1997 .

[15]  M. Pikal,et al.  Mass and heat transfer in vial freeze-drying of pharmaceuticals: role of the vial. , 1984, Journal of pharmaceutical sciences.

[16]  P Sheehan,et al.  Modeling of the primary and secondary drying stages of the freeze drying of pharmaceutical products in vials: numerical results obtained from the solution of a dynamic and spatially multi-dimensional lyophilization model for different operational policies. , 1998, Biotechnology and bioengineering.

[17]  A. I. Liapis,et al.  A theory for the primary and secondary drying stages of the freeze-drying of pharmaceutical crystalline and amorphous solutes: comparison between experimental data and theory , 1994 .

[18]  Hasan Sadikoglu,et al.  Mathematical Modelling of the Primary and Secondary Drying Stages of Bulk Solution Freeze-Drying in Trays: Parameter Estimation and Model Discrimination by Comparison of Theoretical Results With Experimental Data , 1997 .