RESEARCH

 

Preliminary modeling of the perfusion culture of mammalian cells with a spinfilter as retention device

 

Modelación preliminar del cultivo en perfusión de células de mamíferos en tanque agitado con spinfilter como dispositivo de retención

 

 

Luis Y Hernández¹, Abel González¹, Jorge Bouza², Orestes Mayo³, Elena Kulich³, Guido Riera³

¹ Centro de Inmunología Molecular, CIM. Esquina 15 y 216, Atabey, CP 11600, Ciudad de La Habana, Cuba
² Centro Nacional de Investigaciones Científicas, CNIC. Ave. 25 No. 15202 esq. 158, CP 11600, Cubanacán, Playa, Ciudad de La Habana, Cuba
³ Instituto Superior Politécnico José Antonio Echeverría, ISPJAE. Calle 114 No. 11901, entre 119 y 129, Marianao, Ciudad de La Habana, Cuba

 

 

---

ABSTRACT

Specific equations describing the behavior of cell growth and the filtration mechanism of a perfusion culture in stirred tank fermentors at 30 L scale were derived from basic equations for mass balance and mechanical energy in a spinfilter. These equations, when used for modeling the operation process in the Matlab package together with previously reported experimental data, yielded results similar to those of culture kinetics. The operational variables with the highest influence on the process were analyzed with the Matlab module, comparing them to a basal case using the spin rate of the spinfilter and the filtration area usually employed in the production runs as comparison parameters. Two additional comparisons were also performed, using cases with different filtration areas (smaller or higher than that of the basal case) in which the stirring rate was varied to analyze the behavior of the perfusion flow capacity during the run. The influence of the filtration area on the fermentation life was corroborated, with higher values for the latter as filtration area and spin rate of the filter increased.

Keywords: spinfilter, clogging, perfusion flow

---

RESUMEN

Partiendo de las ecuaciones básicas del balance de masa y de energía mecánica aplicado en el spinfilter, se llegó a las ecuaciones particulares que describen el comportamiento del crecimiento celular y al mecanismo de filtración en el cultivo en perfusión en los fermentadores de tanque agitado a escala de banco (30 L). Con las ecuaciones desarrolladas, los datos experimentales reportados en trabajos anteriores y utilizando el programa Matlab, se simuló el proceso de operación y los resultados fueron semejantes a los de la cinética del cultivo. Se analizaron las variables de operación que tienen una marcada influencia durante el proceso, utilizando el módulo del Matlab. Estos se compararon con un caso base, tomando como parámetros: la velocidad de agitación del spinfilter y el área de filtración normalmente empleados en producción. También se compararon con dos casos a dos niveles de áreas de filtración (menor y mayor que la del caso base) en los que se varió la velocidad de agitación para observar el comportamiento de la capacidad de flujo de perfusión durante la corrida. Se comprobó la influencia sobre el tiempo de vida de la fermentación (de la filtración), cuyo valor ascendió para la mayor área de filtración y la mayor velocidad de giro del spinfilter.

Palabras clave: spinfilter, colmatación, flujo de perfusión

---

 

 

INTRODUCTION

The culture of mammalian cells by perfusion in a stirred tank fitted with a spinfilter as a retention device is a methodology aimed at retaining the highest number of cells inside the fermentor, which therefore allows the obtaining of high cell densities and, consequently, high product concentrations in a relatively short period of time, with high volumetric flows in small-scale facilities (1). The filtration time, in the case of the bioreactors of the Center for Molecular Immunology (Havana, Cuba) fitted with spinfilters for the perfusion culture of mammalian cells, is of 18 days. Filtration time, however, can be extended to 90 days according to literature (1-4). Given that no mathematical models for this system have been published, it was decided to model its behavior as a tool to guide further efforts in extending culture time (5) or filtration or perfusion time by varying the spin rate of the spinfilter and its filtration area (2, 3).

MATERIALS AND METHODS

Bioreactor

The fermentor used for the experimental runs (CMF 400, manufactured by Chemap AG) has a total volume of 41 L and an effective volume of 30 L, with a diameter of 0.27 m, a height of 0.7164 m and an effective height of 0.52 m. It has a propeller-type impellent with a diameter of 0.088 m (1, 3).

Spinfilters

The 41 L bioreactors used cylindrical stainless steel spinfilters (Chemap AG) with a diameter of 0.088 m and a height of 0.152 m (1, 3), fitted with a 15 μm pore size stainless steel mesh.

Cell line

The study used the NSO/H7 host cell line (1, 3).

Culture medium

The study used the PFHM II protein-free culture medium (1).

RESULTS AND DISCUSSION

Derivation of hydrodynamic equations

At the Center for Molecular Immunology, perfusion cultures are usually performed in a stirred tank, using a spinfilter as a separation device. A spinfilter is a rotatory cylinder spinning on its axis that allows the continuous separation of cultured cells from the culture media and, therefore, the obtaining of a clarified culture supernatant.

Modeling a perfusion culture in a stirred tank with a spinfilter presents a number of challenges, derived from the simultaneous operation of different phenomena such as centrifugal effects, axial and sweeping forces, among others (3, 6, 7). However, the system can still be analyzed by decomposition into individual parts, modeling: 1) Its behavior as a rotatory filter; 2) Its behavior as a filtrating centrifuge and 3) Elements of mechanic energy balance in the interface of the outer and inner surface of the mesh (Figure 1).

By simultaneously using the equations corresponding to each individual part according to figure 1 and applying the filtration mechanism n = 3/2 (8, 9), the equation of Cozzeny-Karman (10) and the equation of Bernoulli (3, 8, 11, 12), the following is obtained:

Fp: perfusion flow (L/h) q₀*: constant (mL/s) Ks: separation power of the spinfilter tp: time of perfusion (s) ρ_SS: Density of the suspension (kg/m³) ϕ: viscosity of the fluid (Pa-s) The asterisk indicates the exit flow.

Therefore:

Fp= f(X(t), n, S, V, t...)

Where:

X(t): cell density in relationship to elapsed time n: spin rate of the spinfilter S: filtration area V: volume of the bioreactor (L) t: time (h)

Biomass balance equations

Previous work on this subject (1) dismissed the influence of the exchange flow through the filter (F) in order to simplify the calculations; however, this variable has a large influence on the process of clarification and the operation time (13). Therefore, it was decided to include this variable in the biomass balance corresponding to this system (Figures 2 and 3).

Modeling the E1 exponential growth phase

Modeling the exponential growth phase (E1) is simple, since it involves the same equations used for discontinuous cultures (1, 3, 14, 15):

Where:

X₀: starting biomass concentration in the bioreactor (cells/mL) μMAX: maximum specific growth rate (h^-1)

For edge conditions: X₀= 0.5 x 10⁶ cell/mL; X(t₁)= 10⁶ cell/mL; t₀= 24h and t= t₁

Modeling the E2 exponential growth phase under continuous flow

Exponential growth can be modeled in this case with the following differential equation:

Where:

Xs: biomass concentration in the spinfilter (cells/L)

Taking into account that feeding and extraction are described by the same function as that of perfusion flow capacity through the membrane, and assuming that the system reaches 90% retention instantaneously upon starting the perfusion, the differential equation can be arranged to depend on biomass concentration in the bioreactor (3):

Where:

C1= m_MAX= 0.003 h_-1
C*2= 8.561 x 10_-4 (g/cell)
C3= 0.02314 (mL/10₆ cells)
C4= r; r: density of the pure liquid (kg/m₃)
C5= 2 ρq* ₀(1-H); H: humidity of the cake

Where: η: filling coefficient, D: spinfilter diameter (m), α: specific resistance of the apparent cake (m/kg).

Where:

X_K: cell density at perfusion start (around 10⁶ cell/mL)

And for the perfusion flow capacity:

For edge conditions: X(t₁)= 10⁶ cells/mL, X(t₂)= 9 x 10⁶ cells/mL

Modeling the phase of limited growth under continuous flow (stationary phase), E3

During the third phase the stationary state (regarding biomass) is finally reached; that is, biomass remains constant. Therefore:

And for the perfusion flow:

In this phase perfusion flow capacity depends only on time, and therefore:

Fp= Fp(t)

For edge conditions:

X(t³)= X(t²)= 9x10⁶ cells/mL t= t² , Fp= Fp (t²), t= t³ , Fp= Fp (t³)

Limitations of the model

These models have some limitations, since upon starting perfusion 1) the spinfilter is assumed to reach instantaneously a retention of 90%; 2) extraction flow is described by the same function as perfusion capacity through the mesh of the filter; 3) The study was performed only for 3 discrete spin rates of the spinfilter and the same number of different filtration areas; 4) The small-scale study was performed with only one spinfilter mesh; 5) The influence of the spin rate of the impellent over perfusion flow capacity is not taken into account; 6) The influence of the pressure inside the bioreactor on perfusion flow capacity is not taken into account; 7) The washing phase of the fermentor was not modeled and 8) The determination of limiting substrate was not performed, and neither μ = f(S, μ^MAX) nor O² balance nor product formation were adjusted.

Biomass simulation

The simulation of cell density with time in the three growth phases was performed with the real filtration area (S = 4.2 x 10^-2 m²) and the spin rate of the spinfilter (200 rpm), using time (from t = 0 to X(t ) = 9 x 10⁶ cells/mL) as input parameter (16, 17).

Figure 4 shows the stages of exponential growth and exponential growth under constant flow, which for this model fall midrange between the actual run data (Figure 4C), revealing practically the same maximum specific growth rate. In other words, the curve of the model has the same slope as those of the actual runs, with the exception of run 3231TA-0207, which differs in this aspect due to a longer adaptation phase in comparison to the others. Something similar is observed for the stationary phase; although the actual run data follow a curve characterized by a fall that is not present in the model. While this difference is expected, given that the model did not take into account the washing phase that takes place upon complete clogging of the spinfilter, it still constitutes a weakness of the model.

Simulation of perfusion flow capacity

This simulation used as input variables the variation of cell density, the spin rate of the spinfilter, the area of the spinfilter and the operation time, setting the perfusion flow capacity of the system as output variable. Three discrete values for the filtration area were analyzed: that corresponds to the actual filter, and both larger or smaller values, varying on each case the spin rate of the spinfilter (100 rpm, 200 rpm and 300 rpm) (Figure 5) (18, 19).

A similar behavior was observed at the beginning (Figure 6A, B and C) if the fixed area is taken into account, since an increase of the spin rate of the rotatory filter results in an increase of both perfusion flow capacity and filtration time, with a non-linear dependence according to equation 4 (3). The perfusion flow from the pump (Fp_B), ranging from 7.5 to 30 L/day, would be added to each Fp value corresponding to the spin rate of the spinfilter. Although literature mentions this behavior (2, 13, 20-22), this phenomenon is not described analytically, and therefore the results reveal a direct influence of the spin rate on perfusion flow capacity. The effect probably derives from an enlarged zone of laminarity (Figure 7) due to the larger spin rate of the filter (8, 10, 13, 22), which allows a sweep of the cellular profile that approaches the membrane to occlude it, due to the drag of the exchange flow, that appears due to the centrifugal action of the filter itself on the cellular profile.

If the filtration area increases while the spin rate is kept constant, perfusion flow capacity increases almost exponentially (Figure 8). This effect is also described in literature (2, 13, 20-22) although it is not specified analytically, and constitutes another direct result of the influence of the spin rate of the spinfilter and the area on perfusion flow capacity. This must also be a result of the increase in the zone of laminarity (Figure 7) as the spin rate of the filter increases (8, 10, 13, 22) and of the filtration area, allowing a larger sweep of the cellular profile in spite of the drag due to the exchange flow, appearing due to the centrifugal action of the filter itself on the cellular profile.

CONCLUSIONS

This work obtained mathematical models that predict the behavior of perfusion cultures of mammalian cells in stirred tanks for each of their growth stages, coupling kinetic and hydrodynamic equations (1-3, 9). These novel results, which have no precedents in literature, allowed the analysis of the influence of the spin rate of the filter and the filtration area on the perfusion flow capacity of the system. In spite of the limitations of the models, it was possible to determine that an increase in spin rate at a constant filtration area results in an increase in the perfusion capacity of the system, which is even larger if the filtration area also increases, with a concomitant increase in filtration or perfusion times (2, 3, 13, 20).

 

REFERENCES

1. Hernández LY, Castro D, Vitón P, Pérez O, Rodríguez M. Modelation of growth kinetics of mammalian cells in perfusión culture. Biotecnol Apl 2009;26:232-6.

2. Yabannavar VM, Singh V, Connelly N. Mammalian cell retention in a spinfilter perfusion Bioreactor for Mammalian Cell. Biotechnol Bioeng 1992;43:159-64.

3. Hernández LY. Estudio hidrodinámico del cultivo de células de mamífero en perfusión con spinfilter en tanque agitado. Tesis de Maestría. Instituto Superior Pedagógico José Antonio Echeverría. Cuba; 2006.

4. Chu L, Robinson DK. Industrial choices for protein production by large-scale cell cultural. Biotechnology 2001;12:180-7.

5. Varecka R, Scheirer W. Use of rotating wire cage for retention of animal cells in a perfusion fermentor. Dev Biol Stand 1987;66:269-72.

6. Monbouquette. H. G. Modeling high biomass-density cell recycle fermentors. Biotechnol Bioeng 1992;39:498-503.

7. Bailey JE. Mathematical Modeling and analysis in biochemical engineering past accomplishments and future opportunities. Biotechnol Prog 1998;14:8-20.

8. Rushton A, Ward AS, Holdich RG. Solid liquid filtration and separation technology. Second ed. Wiley; 2000.

9. González VA. Estudio del Proceso de formación de tortas de células de mamíferos en la superficie de un spinfilter. Tesis. Instituto Superior Pedagógico José Antonio Echeverría. Cuba; 2004.

10. Bird RB, Stewart WE, Lighfood EN. Transport phenomena. Edición Revolucionaria. La Habana; 1960.

11. Rosabal V. J. Hidrodinámica y Separaciones Mecánicas. Ed. ENPES, La Habana; 1990.

12. Foust SA, Wenzel AL, Clump WC, Maus L, B.Andersen L. Principle of Unit Operation, Edición Revolucionaria. La Habana; 1971.

13. Yashwant MD, Mina DM, Renato F. Practical consideration in operation and scale-up of spinfilter based bioreactors for monoclonal antibody production. Biotechnol Prog 1996;12:57-64.

14. Robinson DK, Memmert KW. Kinetics of recombinant Inmunoglobulin production by mammalian cells in continuous culture. Biotechnol Bioeng 1991;35:972-6.

15. Aiba S, Arthur EH, Nancy FM. Biochemical Engineering. Edición Revolucionaria, La Habana; 1970.

16. Zeng AP, Deckewer WD. Model simulation and analysis f perfusion culture of mammalian cell density. Biotechnol Prog 1999;15:373-82.

17. Ryu DY, Shi Y. Monoclonal antibody productivity and the metabolic pattern of perfusion cultures under varying oxygen tensions. Biotechnol Bioeng 1993;42:430-9.

18. Reuss M, Baltes M, Schneider R, Sturm C. Optimal experimental design for parameter estimation in unstructured growth models. Biotechnol Prog 1994;10:480-8.

19. Ryu DY. Cell cycle kinetics and monoclonal antibody productivity of hybridoma cells during perfusion culture. Biotechnol Bioeng 1994;44:361-7.

20. Yabannavar VM, Singh V, Connelly NV. Scale up of spinfilter perfusion bioreactor for mammalian cell retention. Biotechnol Bioeng 2004;40:925-33.

21. Jan DC-H, Emery AN, Al-Rubeai M. Use of a spin-filter can reduce disruption of hybridoma cells in a bioreactor. Biotechnol Tech 1993;7(5):351-6.

22. Favre E, Thaler T. An engineering analysis of rotating sieves for hybridoma cell retention in stirred tank bioreactors. Cytotechnology 1992;12:180-7.

 

Received in August, 2008.
Accepted for publication in March, 2010.

Luis Y Hernández. Centro de Inmunología Molecular, CIM. Esquina 15 y 216, Atabey, CP 11600, Ciudad de La Habana, Cuba. E-mail: [e-mail]