Cover Feature: Regulation of pyruvate dehydrogenase complex related to lactate switch in CHO cells.

DOI: 10.1002/elsc.202000037 The cover feature visualizes our recent article about the investigation of the regulation of the Pyruvate Dehydrogenase Complex (PDC) during the lactate switch in batch cultures of Chinese Hamster Ovary cells. The relevance of this work to bioprocess engineering is highlighted in the background and the central cellular metabolic regulations are shown symbolically on the right-hand side. The regulation of PDC through phosphorylation was quantified at three regulating sites using a novel indirect flow cytometry protocol, shown as "glowing" antibodies. For details see article DOI 10.1002/elsc.202000037 on page 99.

Regulation of pyruvate dehydrogenase complex related to lactate switch in CHO cells.

The metabolism of Chinese hamster ovary (CHO) cell lines is typically characterized by high rates of aerobic glycolysis with increased lactate formation, known as the "Warburg" effect. Although this metabolic state can switch to lactate consumption, the involved regulations of the central metabolism have only been partially studied so far. An important reaction transferring the lactate precursor, pyruvate, into the tricarboxylic acid cycle is the decarboxylation reaction catalyzed by the pyruvate dehydrogenase enzyme complex (PDC). Among other mechanisms, PDC is mainly regulated by phosphorylation-dephosphorylation at the three sites Ser232, Ser293, and Ser300. In this work, the PDC phosphorylation in antibody-producing CHO DP-12 cell culture is investigated during the lactate switch. Batch cultivations were carried out with frequent sampling (every 6 h) during the transition from lactate formation to lactate uptake, and the PDC phosphorylation levels were quantified using a novel indirect flow cytometry protocol. Contrary to the expected activation of PDC (i.e., reduced PDC phosphorylation) during lactate consumption, Ser293 and Ser300 phosphorylation levels were 33% higher compared to the phase of glucose excess. At the same time, the relative phosphorylation level of Ser232 increased steadily throughout the cultivation (66% increase overall). The intracellular pyruvate was found to accumulate only during the period of high lactate production, while acetyl-CoA showed nearly no accumulation. These results indicate a deactivation of PDC and reduced oxidative metabolism during lactate switch even though the cells undergo a metabolic transition to lactate-based cell growth and metabolism. Overall, this study provides a unique view on the regulation of PDC during the lactate switch, which contributes to an improved understanding of PDC and its interaction with the bioprocess.

Quantification of the dynamics of population heterogeneities in CHO cultures with stably integrated fluorescent markers.

Cell population heterogeneities and their changes in mammalian cell culture processes are still not well characterized. In this study, the formation and dynamics of cell population heterogeneities were investigated with flow cytometry and stably integrated fluorescent markers based on the lentiviral gene ontology (LeGO) vector system. To achieve this, antibody-producing CHO cells were transduced with different LeGO vectors to stably express single or multiple fluorescent proteins. This enables the tracking of the transduced populations and is discussed in two case studies from the field of bioprocess engineering: In case study I, cells were co-transduced to express red, green, and blue fluorescent proteins and the development of sub-populations and expression heterogeneities were investigated in high passage cultivations (total 130 days). The formation of a fast-growing and more productive population was observed with a simultaneous increase in cell density and product titer. In case study II, different preculture growth phases and their influence on the population dynamics were investigated in mixed batch cultures with flow cytometry (offline and automated). Four cell line derivatives, each expressing a different fluorescent protein, were generated and cultivated for different time intervals, corresponding to different growth phases. Mixed cultures were inoculated from them, and changes in the composition of the cell populations were observed during the first 48 h of cultivation with reduced process productivity. In summary, we showed how the dynamics of population heterogeneities can be characterized. This represents a novel approach to investigate the dynamics of cell population heterogeneities under near-physiological conditions with changing productivity in mammalian cell culture processes.

Measurement of length distribution of beta-lactoglobulin fibrils by multiwavelength analytical ultracentrifugation.

The whey protein beta-lactoglobulin is the building block of amyloid fibrils which exhibit a great potential in various applications. These include stabilization of gels or emulsions. During biotechnological processing, high shear forces lead to fragmentation of fibrils and therefore to smaller fibril lengths. To provide insight into such processes, pure straight amyloid fibril dispersions (prepared at pH 2) were produced and sheared using the rotor stator setup of an Ultra Turrax. In the first part of this work, the sedimentation properties of fragmented amyloid fibrils sheared at different stress levels were analyzed with mulitwavelength analytical ultracentrifugation (AUC). Sedimentation data analysis was carried out with the boundary condition that fragmented fibrils were of cylindrical shape, for which frictional properties are known. These results were compared with complementary atomic force microscopy (AFM) measurements. We demonstrate how the sedimentation coefficient distribution from AUC experiments is influenced by the underlying length and diameter distribution of amyloid fibrils.In the second part of this work, we show how to correlate the fibril size reduction kinetics with the applied rotor revolution and the resulting energy density, respectively, using modal values of the sedimentation coefficients obtained from AUC. Remarkably, the determined scaling laws for the size reduction are in agreement with the results for other material systems, such as emulsification processes or the size reduction of graphene oxide sheets.

Near-Physiological Cell Cycle Synchronization with Countercurrent Centrifugal Elutriation.

The bioreactor conditions and cell diversity in mammalian cell cultures are often regarded as homogeneous. Recently, the influence of various kinds of heterogeneities on production rates receives increasing attention. Besides spatial gradients within the cultivation system, the variation between cell populations and the progress of the cells through the cell cycle can affect the dynamics of the cultivation process. Strong metabolic up- and down-regulations leading to variable productivities, even in exponentially growing cell cultures, have been identified in CHO cell cultivations. Consequently, scientific studies of cell cycle-related effects and metabolic regulations require experiments utilizing cell cycle-enriched subpopulations. Importantly, the enrichment procedure itself must not strongly interfere with the cell culture under investigation. Such subpopulations can be generated by near-physiological countercurrent centrifugal elutriation, which is described in the following chapter. At first, a brief overview regarding the cell cycle, currently identified effects and commonly used methods, and their applicability is outlined. Then, the experimental setup and the synchronization itself are explained.

Process-induced cell cycle oscillations in CHO cultures: Online monitoring and model-based investigation.

The influence of process strategies on the dynamics of cell population heterogeneities in mammalian cell culture is still not well understood. We recently found that the progression of cells through the cell cycle causes metabolic regulations with variable productivities in antibody-producing Chimese hamster ovary (CHO) cells. On the other hand, it is so far unknown how bulk cultivation conditions, for example, variable nutrient concentrations depending on process strategies, can influence cell cycle-derived population dynamics. In this study, process-induced cell cycle synchronization was assessed in repeated-batch and fed-batch cultures. An automated flow cytometry set-up was developed to measure the cell cycle distribution online, using antibody-producing CHO DP-12 cells transduced with the cell cycle-specific fluorescent ubiquitination-based cell cycle indicator (FUCCI) system. On the basis of the population-resolved model, feeding-induced partial self-synchronization was predicted and the results were evaluated experimentally. In the repeated-batch culture, stable cell cycle oscillations were confirmed with an oscillating G1 phase distribution between 41% and 72%. Furthermore, oscillations of the cell cycle distribution were simulated and determined in a (bolus) fed-batch process with up to 25×106 cells/ml. The cell cycle synchronization arose with pulse feeding only and ceased with continuous feeding. Both simulated and observed oscillations occurred at higher frequencies than those observable based on regular (e.g., daily) sample analysis, thus demonstrating the need for high-frequency online cell cycle analysis. In summary, we showed how experimental methods combined with simulations enable the improved assessment of the effects of process strategies on the dynamics of cell cycle-dependent population heterogeneities. This provides a novel approach to understand cell cycle regulations, control cell population dynamics, avoid inadvertently induced oscillations of cell cycle distributions and thus to improve process stability and efficiency.

Direct and highly sensitive measurement of fluorescent molecules in bulk solutions using flow cytometry.

Utilizing flow cytometry to monitor progress of bulk biochemical reactions and concentration of chemical species normally relies on the utilization of cells carrying intrinsic fluorescence or modified beads. We present a method for a simple measurement of the fluorescent marker molecule fluorescein and GFPuv in bulk solutions with high sensitivity using a CytoFLEX flow cytometer and without the need for modified beads. Polystyrene beads were used to trigger measurements based on their high scatter signal, to detect the fluorescence signal from two different fluorophores present in the sample solution. We report sensitivities of 33 pg/mL for fluorescein and 50 ng/mL for GFPuv. This method is comparable in sensitivity to a typical spectrometric fluorescence assay tested with fluorescein, and approximately ten times more sensitive for the measurement of GFPuv. PEG was added to the sample at a low concentration of 0.001% (w/v) to block unspecific GFPuv binding to the beads. The method was further applied to measure the GFPuv concentration in crude cell lysate samples used for cell free protein expression. An advantage of this method over spectrometric assays is the ability to differentiate signal subpopulations in the sample based on their individual fluorescence intensities.

Toward Multiscale Modeling of Proteins and Bioagglomerates: An Orientation-Sensitive Diffusion Model for the Integration of Molecular Dynamics and the Discrete Element Method.

Most processes involved in biological and biotechnological systems spread over many scales in space and time. For example, the interaction of multiple enzymes in heterogeneous enzymatic agglomerates or clusters, necessary for efficient enzymatic conversion, is of high interest for research and enzyme engineering. In order to understand and predict their overall behavior and performance, it is important to describe these scales as completely as possible, known as multiscale modeling. While many different approaches have been presented in recent years, knowledge about protein formation and bioagglomeration at the micro scale is still very limited. In an attempt to address such systems, we propose a bottom- up multiscale modeling methodology, bridging the gaps between molecular dynamics (MD) with an explicit solvent and the larger scale discrete element method (DEM) using an implicit solvent and abstracting macromolecules (e.g., proteins) as objects with anisotropic properties. We term this approach the molecular discrete element method (MDEM). For this, we present an orientation-sensitive diffusion model for DEM, which describes the dynamics of anisotropic translational and rotational diffusion, while implicitly considering solvent molecules and enforcing a canonical ensemble. A general-purpose model and parametrization approach is presented, which can be used to simulate any process involving diffusion of discrete particles. Effects of temperature and viscosity changes can be considered, and guidance is provided concerning time step selection. This model is generally applicable and serves as a precondition to enforce the proper dynamics (i.e., diffusion characteristics and canonical ensemble, similar to a thermostat in MD) for the proposed multiscale modeling methodology with anisotropic properties. Thereby, it presents a first step toward modeling at the micro scale and is integral to enforcing dynamics of such systems and therefore extensively validated. As a next step, interaction models are to be defined and added to the presented model. In comparison to atomistic and coarse-grained (CG) MD, a speedup of 5-7 orders of magnitude can be achieved. The approach is demonstrated on multiple components of the pyruvate dehydrogenase enzyme complex, a multienzymatic machinery that involves very different types of enzymes and is of high value to further elucidate the mechanisms of bioagglomeration and metabolic channeling.

Model-based identification of cell-cycle-dependent metabolism and putative autocrine effects in antibody producing CHO cell culture.

The understanding of cell-cycle-dependent population heterogeneities in mammalian cell culture and their influence on production rates is still limited. Furthermore, metabolic regulations arising from self-expressed signaling factors (autocrine/autoinhibitory factors) have been postulated in the past, but no determination of such effects have been made so far for fast-growing production Chinese hamster ovary (CHO) cells in chemically defined media. In this study, a novel approach combining near-physiological treatment of cells (including synchronization), population resolved mechanistic modeling and statistical analysis was developed to identify population inhomogeneities. Cell-cycle-dependent population dynamics and metabolic regulations due to a putative autocrine factor were examined and their impact on the metabolic rates and antibody production of near-physiologically synchronized CHO DP-12 cell cultures was determined. To achieve this, a population resolved model was extended to describe putative autocrine-dependentt and cell-cycle-related metabolic rates for glucose, glutamine, lactate, ammonia, and antibody production. The model parameters were estimated based on data of two repeated batch cultivations (three batches each), with main substrates in excess and potentially inhibiting waste products (lactate and ammonium) controlled within narrow ranges. Significant changes, due to a putative autocrine factor, were identified for lactate and ammonia formation and antibody production. The cell growth and the uptake of glucose and glutamine were only partially affected by the putative autocrine under the given conditions. The results indicate the presence of a self-expressed autocrine factor and its strong impact on the metabolism of CHO DP-12 cells. Furthermore, glucose and glutamine uptake, as well as the formation of ammonium and antibody were found to be significantly cell-cycle-dependent. The combined approach has a strong potential to improve the understanding of the interplay of population heterogeneities and signal factors in mammalian cell culture.

Full Enzyme Complex Simulation: Interactions in Human Pyruvate Dehydrogenase Complex.

The pyruvate dehydrogenase complex (PDC) is a large macromolecular machine consisting of dozens of interacting enzymes that are connected and regulated by highly flexible domains, also called swinging arms. The overall structure and function of these domains and how they organize the complex function have not been elucidated in detail to date. This lack of structural and dynamic understanding is frequently observed in multidomain enzymatic complexes. Here we present the first full and dynamic structural model of full human PDC (hPDC), including binding of the linking arms to the surrounding E1 and E3 enzymes via their binding domains with variable stoichiometries. All of the linking domains were modeled at atomistic and coarse-grained levels, and the latter was parametrized to reproduce the same properties of those from the atomistic model. The radii of gyration of the wild-type full complex and functional trimeric subunits were in agreement with available experimental data. Furthermore, the E1 and E3 population effect on the overall structure of the full complex was studied. The results indicated that decreasing the number of E1s increases the flexibility of the now nonoccupied arms. Furthermore, their flexibility depends on the presence of other E1s and E3s in the vicinity, even if they are associated with other arms. As one consequence, the radius of gyration decreases with decreasing number of E1s. This effect also provides an indication of the optimal configuration of E1 and E3 on the basis of the assumption that a certain stability of the enymatic cloud is necessary to avoid free metabolic diffusion of intermediates (metabolic channeling). Our approach and results open a window for future enzyme engineering in a more effective way by evaluating the effect of different linker arm lengths, flexibilities, and combinations of mutations on the activity of other complex enzymes that involve flexible domains, including for example processive enzymes.

CHO cells engineered for fluorescence read out of cell cycle and growth rate in real time.

For efficient production of recombinant proteins by mammalian cells in a bioreactor, optimal growth rates are required and represent the most important process parameter. We present the first successful attempt to monitor the growth behavior and cell cycle state of a mammalian production relevant cell line under bioreactor cultivation conditions up to 1.2 l, utilizing a fluorescent read-out without the need of additional staining or marking. For this purpose, we developed two new production relevant cell line derivatives (CHO-K1 FUCCI CM & CHO-K1 FUCCI CN) and corresponding analytical methods. The approach is easily scalable, applicable to mammalian recombinant protein production cell lines, and it allows for real-time monitoring using appropriate fluorescence probes. It is based on the Ubiquitination-based Cell Cycle Indicator (FUCCI) system developed by Miyawaki et al. CHO-K1 was chosen as a model cell line due to its close relationship to several production cell lines.1 We defined a new process parameter ired , a quantitative and numerically robust representation of the cell cycle distribution, and demonstrate it to be linearly correlated with the cell cycle state and inversely related to the real time growth rate. Detection of growth rate limitations is possible earlier than using cell-count-based approaches. Analytics were compatible with bulk fluorescence methods, using a plate reader as well as a flow cytometer. For future real time applications in industry scale bioreactors we recommend the use of on-line or at-line fluorescence probes. © 2017 American Institute of Chemical Engineers Biotechnol. Prog., 33:1408-1417, 2017.

Investigation of Core Structure and Stability of Human Pyruvate Dehydrogenase Complex: A Coarse-Grained Approach.

The human pyruvate dehydrogenase complex (hPDC) is a large macromolecular machine, and its unique structural and functional properties make it a versatile target for manipulation aiming for the design of new types of artificial multienzyme cascades. However, model-based and hence systematic understanding of the structure-function relationship of this kind of complexes is yet poor. However, with new structure data, modeling techniques, and increasing computation power available, this shortfall is about to cease. Recently, we have built new atomistic models of E2/E3BP, the two subunits of the massive hPDC core. Here, we present developed coarse-grained models of the same, on the basis of which we built and simulated the full core of hPDC, which is so far the first simulation of the catalytic core of any member in the branched-chain α-keto acid dehydrogenase complex family. We explored the stability of two previously proposed substitutional models of hPDC core: 40E2+20E3BP and 48E2+12E3BP. Our molecular dynamics simulations showed a higher stability and sphericity for the second model. Our core's radius of gyration was found to be in strong agreement with previously published experimental dynamic light scattering (DLS) data. Finally, in the direction of our experimental effort to design active minimized complexes, we simulated C-terminal truncated E2/E3BP cores of different lengths, which clearly showed the instability of the core assembly and symmetry due to subunit separations. This correlated very well with the findings of our experimental investigations by analysis of DLS data for variable truncation lengths. The use of polarizable water and an increased total ion concentration did not lead to a substantially different initial stability of the truncated mutants compared to that of standard MARTINI water; however, a different rearrangement behavior of the fragments was observed. The obtained structure models will serve as a basis for full complex simulations in the future, providing the possibility for the model-assisted targeted manipulation of such a complex enzymatic system.

Weak cell cycle dependency but strong distortive effects of transfection with Lipofectamine 2000 in near-physiologically synchronized cell culture.

Previously, we reported a method to generate and validate cell cycle-synchronized cultures of multiple mammalian suspension cell lines under near-physiological conditions. This method was applied to elucidate the putative interdependencies of the cell cycle and recombinant protein expression in the human producer cell line HEK293s using Lipofectamine 2000 and the reporter plasmid pcDNA3.3 enhanced green fluorescent protein, destabilized using PEST sequence. A population-resolved modeling approach was applied to quantitatively assess putative variations of cell cycle dependent expression rates based on the obtained experimental data. We could not confirm results published earlier by other groups, based on nonphysiological synchronization attempts, reporting transfection efficiency being strongly dependent on the cell cycle phase at transfection time point. On the other hand, it is demonstrated that transfection and protein expression distort the progression of the cell cycle.

Reengineering of the human pyruvate dehydrogenase complex: from disintegration to highly active agglomerates.

The pyruvate dehydrogenase complex (PDC) plays a central role in cellular metabolism and regulation. As a metabolite-channeling multi-enzyme complex it acts as a complete nanomachine due to its unique geometry and by coupling a cascade of catalytic reactions using 'swinging arms'. Mammalian and specifically human PDC (hPDC) is assembled from multiple copies of E1 and E3 bound to a large E2/E3BP 60-meric core. A less restrictive and smaller catalytic core, which is still active, is highly desired for both fundamental research on channeling mechanisms and also to create a basis for further modification and engineering of new enzyme cascades. Here, we present the first experimental results of the successful disintegration of the E2/E3BP core while retaining its activity. This was achieved by C-terminal α-helixes double truncations (eight residues from E2 and seven residues from E3BP). Disintegration of the hPDC core via double truncations led to the formation of highly active (approximately 70% of wildtype) apparently unordered clusters or agglomerates and inactive non-agglomerated species (hexamer/trimer). After additional deletion of N-terminal 'swinging arms', the aforementioned C-terminal truncations also caused the formation of agglomerates of minimized E2/E3BP complexes. It is likely that these 'swinging arm' regions are not solely responsible for the formation of the large agglomerates.

Human Pyruvate Dehydrogenase Complex E2 and E3BP Core Subunits: New Models and Insights from Molecular Dynamics Simulations.

Targeted manipulation and exploitation of beneficial properties of multienzyme complexes, especially for the design of novel and efficiently structured enzymatic reaction cascades, require a solid model understanding of mechanistic principles governing the structure and functionality of the complexes. This type of system-level and quantitative knowledge has been very scarce thus far. We utilize the human pyruvate dehydrogenase complex (hPDC) as a versatile template to conduct corresponding studies. Here we present new homology models of the core subunits of the hPDC, namely E2 and E3BP, as the first time effort to elucidate the assembly of hPDC core based on molecular dynamic simulation. New models of E2 and E3BP were generated and validated at atomistic level for different properties of the proteins. The results of the wild type dimer simulations showed a strong hydrophobic interaction between the C-terminal and the hydrophobic pocket which is the main driving force in the intertrimer binding and the core self-assembly. On the contrary, the C-terminal truncated versions exhibited a drastic loss of hydrophobic interaction leading to a dimeric separation. This study represents a significant step toward a model-based understanding of structure and function of large multienzyme systems like PDC for developing highly efficient biocatalyst or bioreaction cascades.

Synchronized mammalian cell culture: part I--a physical strategy for synchronized cultivation under physiological conditions.

Conventional analysis and optimization procedures of mammalian cell culture processes mostly treat the culture as a homogeneous population. Hence, the focus is on cell physiology and metabolism, cell line development, and process control strategy. Impact on cultivations caused by potential variations in cellular properties between different subpopulations, however, has not yet been evaluated systematically. One main cause for the formation of such subpopulations is the progress of all cells through the cell cycle. The interaction of potential cell cycle specific variations in the cell behavior with large-scale process conditions can be optimally determined by means of (partially) synchronized cultivations, with subsequent population resolved model analysis. Therefore, it is desirable to synchronize a culture with minimal perturbation, which is possible with different yield and quality using physical selection methods, but not with frequently used chemical or whole-culture methods. Conventional nonsynchronizing methods with subsequent cell-specific, for example, flow cytometric analysis, can only resolve cell-limited effects of the cell cycle. In this work, we demonstrate countercurrent-flow centrifugal elutriation as a useful physical method to enrich mammalian cell populations within different phases of a cell cycle, which can be further cultivated for synchronized growth in bioreactors under physiological conditions. The presented combined approach contrasts with other physical selection methods especially with respect to the achievable yield, which makes it suitable for bioreactor scale cultivations. As shown with two industrial cell lines (CHO-K1 and human AGE1.HN), synchronous inocula can be obtained with overall synchrony degrees of up to 82% in the G1 phase, 53% in the S phase and 60% in the G2/M phase, with enrichment factors (Ysync) of 1.71, 1.79, and 4.24 respectively. Cells are able to grow with synchrony in bioreactors over several cell cycles. This strategy, combined with population-resolved model analysis and parameter extraction as described in the accompanying paper, offers new possibilities for studies of cell lines and processes at levels of cell cycle and population under physiological conditions.

Synchronized mammalian cell culture: part II--population ensemble modeling and analysis for development of reproducible processes.

The consideration of inherent population inhomogeneities of mammalian cell cultures becomes increasingly important for systems biology study and for developing more stable and efficient processes. However, variations of cellular properties belonging to different sub-populations and their potential effects on cellular physiology and kinetics of culture productivity under bioproduction conditions have not yet been much in the focus of research. Culture heterogeneity is strongly determined by the advance of the cell cycle. The assignment of cell-cycle specific cellular variations to large-scale process conditions can be optimally determined based on the combination of (partially) synchronized cultivation under otherwise physiological conditions and subsequent population-resolved model adaptation. The first step has been achieved using the physical selection method of countercurrent flow centrifugal elutriation, recently established in our group for different mammalian cell lines which is presented in Part I of this paper series. In this second part, we demonstrate the successful adaptation and application of a cell-cycle dependent population balance ensemble model to describe and understand synchronized bioreactor cultivations performed with two model mammalian cell lines, AGE1.HNAAT and CHO-K1. Numerical adaptation of the model to experimental data allows for detection of phase-specific parameters and for determination of significant variations between different phases and different cell lines. It shows that special care must be taken with regard to the sampling frequency in such oscillation cultures to minimize phase shift (jitter) artifacts. Based on predictions of long-term oscillation behavior of a culture depending on its start conditions, optimal elutriation setup trade-offs between high cell yields and high synchronization efficiency are proposed.

Mammalian cell culture synchronization under physiological conditions and population dynamic simulation.

The overall behavior of cell cultures is determined by the actions and regulations of all cells and their interaction in a mixed population. However, the dynamics caused by diversity and heterogeneity within cultures is often neglected in the study of cell culture processes. Usually, a bulk behavior is assumed, although heterogeneity prevails in most cases. It is, however, not valid to conclude from the bulk behavior to the single cell behavior. Instead, it is necessary to include the behavior and kinetics of subpopulations and their interactions into models in order to elucidate the dynamic effects occurring in typical cell cultures. Heterogeneity in cell cultures is largely caused by the progress of the cell cycle. Cell cycle-dependent dynamics resulting for example in variable transfection efficiencies or expression bistability have recently attracted attention. In order to elucidate cell cycle-dependent regulations in cell cultures, it is desirable to synchronize a culture with minimal perturbation, which is possible with different yield and quality using physical methods, but not possible for frequently used chemical, or whole-culture methods. Then, the culture is cultivated again under physiological conditions and subpopulation-resolved analysis and modeling approaches are applied. This should allow to account for the variable contributions of subpopulations to the whole behavior and also for obtaining hereto unaccessible dynamic information of cellular regulation. In this short review, we summarize techniques and key issues to be considered for successful synchronization, cultivation, and modeling in order to achieve the goal of better understanding cell culture at a population level.

Modeling of intracellular transport and compartmentation.

The complexity and internal organization of mammalian cells as well as the regulation of intracellular transport processes has increasingly moved into the focus of investigation during the past two decades. Advanced staining and microscopy techniques help to shed light onto spatial cellular compartmentation and regulation, increasing the demand for improved modeling techniques. In this chapter, we summarize recent developments in the field of quantitative simulation approaches and frameworks for the description of intracellular transport processes. Special focus is therefore laid on compartmented and spatiotemporally resolved simulation approaches. The processes considered include free and facilitated diffusion of molecules, active transport via the microtubule and actin filament network, vesicle distribution, membrane transport, cell cycle-dependent cell growth and morphology variation, and protein production. Commercially and freely available simulation packages are summarized as well as model data exchange and harmonization issues.

Spatiotemporal modeling and analysis of transient gene delivery.

A quantitative and mechanistic understanding of intracellular transport processes in eukaryotic cells during transient transfection is an important prerequisite for the systematic and specific optimization of transient gene expression procedures for pharmaceutic and industrial protein production. There is evidence that intracellular transport processes during gene delivery and their regulation may have significant influence on the transfection efficiency. This contribution describes a compartmented, spatiotemporally resolved and stochastic modeling approach that describes intracellular transport processes responsible for gene delivery during transient transfection. It enables a detailed prediction and analysis and identification of potential bottlenecks. This model is currently being adapted to a model cell line, HEK293s. The simulated results are compared with experimental quantitative polymerase chain reaction (qPCR) data and confocal imaging data obtained with transfected and stained HEK293 cells. Global parameter estimation is performed to qPCR data based on two different novel plasmid constructs in order to identify candidates for plasmid-specific transport parameter variations. The influence of the specific property of HEK293 cells to grow in clusters is investigated and the impact of active microtubule transport depending on cell morphology and clustering is examined. A general sensitivity analysis allows for the identification of the sensitive parameters.