Continuous maturation of proliferating erythroid precursors

Abstract. This study examines published steady state cell kinetic (mitotic and DNA synthesis phase) data from the recognizable proliferating erythroid precursors in humans, rats, and guinea-pigs, and human neutrophilic precursors, for consistency with a continuous maturation-proliferation model of the cell cycle. We find that these data are completely consistent with the hypothesis that maturation between morphological compartments may take place at any point in the cell cycle. A number of compartmental parameters are derived under this assumption.     

Calcium transport and distribution in Ehrlich mouse ascites tumor cells

Data from isotopic uptake experiments were used to measure calcium fluxes and compartment sizes in ascites tumor cells. The data were analyzed with two kinetic models, A and B. In 80% of the experiments model A, which is based on one exchangeable calcium compartment, was rejected in favor of Model B, which predicts two exchangeable compartments. A statistical evaluation of the model's performance, when fit to the experimental data was used to select between the two models. The results show that calcium was distributed between three cellular compartments in the ratio, non-exchangeable (88%): rapidly exchanging (7%): slowly exchanging (5%). The undirectional fluxes suggested that calcium transport could be described as a series system with the temporal sequence: environment ? rapidly exchanging ? slowly exchanging.     

Analysis of Chinese hamster ovary cell metabolism through a combined computational and experimental approach

Optimization of cell culture processes can benefit from the systematic analysis of experimental data and their organization in mathematical models, which can be used to decipher the effect of individual process variables on multiple outputs of interest. Towards this goal, a kinetic model of cytosolic glucose metabolism coupled with a population-level model of Chinese hamster ovary cells was used to analyse metabolic behavior under batch and fed-batch cell culture conditions. The model was parameterized using experimental data for cell growth dynamics, extracellular and intracellular metabolite profiles. The results highlight significant differences between the two culture conditions in terms of metabolic efficiency and motivate the exploration of lactate as a secondary carbon source. Finally, the application of global sensitivity analysis to the model parameters highlights the need for additional experimental information on cell cycle distribution to complement metabolomic analyses with a view to parameterize kinetic models.     

Developmental compartments and planar polarity in Drosophila

BACKGROUND: Planar polarity refers to the asymmetry of a cell within the plane of the epithelium; for example, cells may form hairs that point in a posterior direction, or cilia may beat in one way. This property implies that cells have information about their orientation; we wish to understand the nature of this information. Relevant also is the body plan of insects, which, in the ectoderm and somatic mesoderm, consists of a chain of alternating anterior and posterior compartments - basic units of development with independent cell lineage and subject to independent genetic control.RESULTS: Using the abdomen of adult Drosophila, we have taken genes required for normal polarity and either removed the gene or constitutively expressed it in small clones of cells and observed the effects on polarity. Hitherto, all such studies of polarity genes have not found any difference of behavior between the different compartments. We report here that the three genes, four-jointed, dachsous, and fat, cause opposite effects in anterior and posterior compartments. For example, in anterior compartments, clones ectopically expressing four-jointed reverse the polarity of cells in front of the clone, while, in posterior compartments, they reverse behind the clone. These three genes have been reported by others to be functionally linked.CONCLUSIONS: This discovery impacts on models of how cells read polarity. At the heart of one class of models is the hypothesis that cell polarity is determined by the vector of a morphogen gradient. Here, we present evidence that cell polarity in the abdomen depends on at least two protein gradients (Fj and Ds), each of which is reflected at compartment borders. Consequently, these gradients have opposing slopes in the two compartments. Because all polarized structures made by abdominal cells point posteriorly, we surmise that cells in each compartment are programmed to interpret these protein gradients with opposite signs, pointing up the gradient in one compartment and down the gradient in the other.     

Kinetic model of cell growth and secondary metabolite synthesis in plant cell culture ofThalictrum rugosum

A structured kinetic model was proposed to describe cell growth and synthesis of a secondary metabolite, berberine, in batch suspension culture ofThalictrum rugosum. The model was developed by representing the physiological state of the cell in terms of the activity and the viability, which can be estimated using the culture fluorescence measurement. In the proposed model, the cells were divided into three types; active-viable, nonactive-viable, and dead cells. The model was formulated in terms of cell growth (dry/fresh weight, activity, and viability), carbon source utilization (sucrose, glucose and fructose), and product formation (intracellular and extracellular berberine). The concept of cell expansion and the death phase were also included in this model to describe the sugar accumulation and the release of intracellular berberine into medium by cell lysis, respectively. The parameters used in this model were estimated based on the experimental results in conjunction with numerical optimization techniques. Satisfactory agreement between the model and experimental data was obtained. The proposed model could accurately predict cell growth and product synthesis as well as the distribution of the secondary metabolite between the cell and the medium. It is suggested that the proposed model could be extended as a useful framework for quantitative analysis of physiological characteristics in the other plant cell culture systems.     

Cell proliferation in the erythroid compartment of guinea-pig bone marrow: studies with 3H-thymidine.

Single and multiple injections of 3H-TdR have been used for measuring the rate of proliferation in morphologically defined cell populations of guinea-pig bone marrow that are committed to erythroid differentiation. The conclusions are based on the analysis of absolute cell numbers in the maturational compartments, the labeling and mitotic indices, labeled mitotic curves, pulse and chase grain counts over dividing and interphase cells, and on the rate or labeling during multiple, repeated injections of 3H-TdR. The average duration of S and the rate of cycling is similar in all maturational compartments of the erythrom. The majority of cells progress to the next maturational compartment by the time they divide for the second time. All proerythroblasts and basophilic erythroblasts are in cycle. Polychromatic erythroblasts incapable of incorporating 3H-TdR reach the orthochromatic population in the span of 5-6 hr. The orthochromatic population is renewed every 20-24 hr. The number of divisions between the proerythroblast and orthochromatic erythroblast does not exceed four and some cells may undergo only two divisions during the maturation pathway. Cell input from a progenitor cell population contributes to the maintenance of the erythron. The kinetic behavior of progenitor cells is similar to that of proerythroblasts. By the time of their second division, progenitor cells may reach either the proerythroblast or basophilic erythroblast compartments. The kinetic behavior of basophilic transitional cells corresponds to the predicted behavior of the erythroblast progenitor cell pool. Several of the conclusions are based on the assumption that grain count halving is the result of cell division. In view of the evidence discussed, this assumption in the present studies seems justified.     

DMSO及RA对GPI-80分子表达的不同影响

以往的研究表明GPI-80的表达可能与髓系细胞的分化相关。DMSO及RA是两种不同的中性粒细胞的诱导分化剂,均可刺激HL-60白血病细胞向中性粒细胞分化。GPI-80是人糖基化磷脂酰肌醇锚糖蛋白,被认为是潜在的β2-黏合素分子依赖的白细胞黏附的调节剂,主要在人中性粒细胞上表达。本研究通过RT—PCR、流式细胞仪及Western—blot分析,检测分化细胞的GPI-80表达,并分析GPI-80的表达与CD11b及CD71表达之间的关系。结果表明GPI-80在RA诱导的类中性粒细胞上只有mRNA水平上的微弱表达,用流式细胞仪和Western—blot分析均检测不到,且RA可抑制GPI-80的表达;相反GPI-80在DMSO诱导的类中性粒细胞上有明显的表达,且随DMSO的浓度增加及诱导时间的延长而增强。GPI-80的表达出现在CD11b上调表达及CD71下调表达之后,提示GPI-80表达与DMSO诱导分化的类中性粒细胞的成熟密切相关。RA不能明确诱导GPI-80的表达,反而抑制GPI-80的表达,提示可能两者诱导HL-60细胞分化时所激活的信号传递通路不同。     

Kinetics of phosphate limited algal growth

The kinetics of phosphate limited growth of two green algae Chorella pyrenoidosa and Selenastrum capricornutum have been studied in chemostats. Several kinetic models which express the specific growth rate as a function of the intracellular phosphours content have been examined, and one of the models was found to be significantly better than the other models. The principles of this model were described in a recent paper by Nyholm. The kinetics of phosphate uptake have been investigated by adding pulses of phosphate to the chemostats. The uptake by phosphours deficient cells could be described by Michaelis–Menten kinetics for phosphate concentrations below approximately 500 μg P/liter. Further, with the assumption of a discontinuous adjustment of the uptake rate at the onset of phosphours deficiency, a complete kinetic model for growth and phosphate removal is proposed. The mean cell size and the contents of chlorophyll and RNA per unit dry weight have been measured for C. pyrenoidosa as a function of the dilution rate.     

Kinetics and cellularity of erythroblasts in the growing rat

Abstracts. Parameters of recognizable erythroid cell proliferation were measured in four groups of normal rats weighing 50, 100, 150 and 300 g in order to provide a comprehensive set of data suitable for a re-investigation of erythropoietic models. Total erythroblast cellularity was measured by the ⁵⁹Fe technique, DNA synthesis time by quantitative ¹⁴C-autoradiography, and the erythrocyte production rate was derived from the increase with time of the erythrocyte labelling index after repeated injections of ³H-leucine. Furthermore, the relative erythroblast density was determined in the various morphological compartments. From the total erythroid cell mass in DNA synthesis and the absolute erythrocyte production rate, figures were derived for the mean DNA synthesis time of erythroid cells and compared with the directly measured ones. The discrepancies in all weight groups between direct and indirect determination of DNA synthesis time were considerable. In a previous study re-evaluation of comparable data in literature had revealed comparable inconsistencies. Since a critical discussion of possible errors in the experimental techniques does not indicate data acquisition to be the principal source of disagreement it is concluded that the type of model applied to describe how cells pass the boundaries of morphological cell compartments is of high significance. Models based on a sequential flux of cells through the individual compartments are inadequate for evaluation of the presented set of data.     

Binding of cell-penetrating penetratin peptides to plasma membrane vesicles correlates directly with cellular uptake

Cell-penetrating peptides (CPPs) gain access to intracellular compartments mainly via endocytosis and have capacity to deliver macromolecular cargo into cells. Although the involvement of various endocytic routes has been described it is still unclear which interactions are involved in eliciting an uptake response and to what extent affinity for particular cell surface components may determine the efficiency of a particular CPP. Previous biophysical studies of the interaction between CPPs and either lipid vesicles or soluble sugar-mimics of cell surface proteoglycans, the two most commonly suggested CPP binding targets, have not allowed quantitative correlations to be established. We here explore the use of plasma membrane vesicles (PMVs) derived from cultured mammalian cells as cell surface models in biophysical experiments. Further, we examine the relationship between affinity for PMVs and uptake into live cells using the CPP penetratin and two analogs enriched in arginines and lysines respectively. We show, using centrifugation to sediment PMVs, that the amount of peptide in the pellet fraction correlates linearly with the degree of cell internalization and that the relative efficiency of all-arginine and all-lysine variants of penetratin can be ascribed to their respective cell surface affinities. Our data show differences between arginine- and lysine-rich variants of penetratin that has not been previously accounted for in studies using lipid vesicles. Our data also indicate greater differences in binding affinity to PMVs than to heparin, a commonly used cell surface proteoglycan mimic. Taken together, this suggests that the cell surface interactions of CPPs are dependent on several cell surface moieties and their molecular organization on the plasma membrane.     

CELL PROLIFERATION IN THE ERYTHROID COMPARTMENT OF GUINEA-PIG BONE MARROW: STUDIES WITH 3H-THYMIDINE

Single and multiple injections of ³H-TdR have been used for measuring the rate of proliferation in morphologically defined cell populations of guinea-pig bone marrow that are committed to erythroid differentiation. The conclusions are based on the analysis of absolute cell numbers in the maturational compartments, the labeling and mitotic indices, labeled mitotic curves, pulse and chase grain counts over dividing and interphase cells, and on the rate of labeling during multiple, repeated injections of ³H-TdR. The average duration of S and the rate of cycling is similar in all maturational compartments of the erythron. The majority of cells progress to the next maturational compartment by the time they divide for the second time. All proerythroblasts and basophilic erythroblasts are in cycle. Polychromatic erythroblasts incapable of incorporating ³H-TdR reach the orthochromatic population in the span of 5–6 hr. The orthochromatic population is renewed every 20–24 hr. The number of divisions between the proerythroblast and orthochromatic erythroblast does not exceed four and some cells may undergo only two divisions during the maturation pathway. Cell input from a progenitor cell population contributes to the maintenance of the erythron. The kinetic behavior of progenitor cells is similar to that of proerythroblasts. By the time of their second division, progenitor cells may reach either the proerythroblast or basophilic erythroblast compartments. The kinetic behavior of basophilic transitional cells corresponds to the predicted behavior of the erythroblast progenitor cell pool. Several of the conclusions are based on the assumption that grain count halving is the result of cell division. In view of the evidence discussed, this assumption in the present studies seems justified.     

Translating biochemical network models between different kinetic formats

Mechanistic biochemical network models describe the dynamics of intracellular metabolite pools in terms of substance concentrations, stoichiometry and reaction kinetics. Data from stimulus response experiments are currently the most informative source for in-vivo parameter estimation in such models. However, only a part of the parameters of classical enzyme kinetic models can usually be estimated from typical stimulus response data. For this reason, several alternative kinetic formats using different “languages” (e.g. linear, power laws, linlog, generic and convenience) have been proposed to reduce the model complexity. The present contribution takes a rigorous “multi-lingual” approach to data evaluation by translating biochemical network models from one kinetic format into another. For this purpose, a new high-performance algorithm has been developed and tested. Starting with a given model, it replaces as many kinetic terms as possible by alternative expressions while still reproducing the experimental data. Application of the algorithm to a published model for Escherichia coli's sugar metabolism demonstrates the power of the new method. It is shown that model translation is a powerful tool to investigate the information content of stimulus response data and the predictive power of models. Moreover, the local and global approximation capabilities of the models are elucidated and some pitfalls of traditional single model approaches to data evaluation are revealed.     

A simple model for cell volume and developmental compartments in nutrient limited cyclostat cultures of algae

Daily light-dark cycles can entrain cell growth and division cycles in populations of algae growing in nutrient limited continuous cultures, or cyclostats. In this study a simple model for the flux of cells between discrete developmental stages is formulated for periodic cyclostat cultures of algae. Cell growth, in terms of volume, was set as being constant within a given developmental compartment, but variable between compartments. Growth within a given compartment or transition between compartments was restricted to specific intervals of the subjective day. The model was calibrated to phosphate limited cyclostat growth of Euglena gracilis, with the intervals for transition between compartments fixed at the times relative to the subjective dawn corresponding to critical transition points in the phased cell cycle of this organism. The model output for mean population volume per cell agreed well with experimental data. Although greatly simplified, the periodic behavior of the model volume frequency distributions for the discrete compartments provide reasonable approximation of experimentally determined distributions.     

Reassessment of fluid-phase endocytosis and diacytosis in monolayer cultures of human fibroblasts

We have investigated the kinetics of fluid-phase endocytosis and diacytosis in confluent monolayers of human fibroblasts by comparing the behavior of three markers that have been previously used to study this process: [14C]sucrose, 125I-labeled polyvinylpyrrolidone ([125I]PVP), and Lucifer Yellow. Three distinct kinetic compartments were observed with all markers. The first was relatively large (10-60 fl/cell), reached steady state within 15 min at 37 degrees C, and was rapidly lost from monolayers after removing the markers at 37 degrees C but not at 0 degree C. These properties indicate that this compartment is the same as that previously proposed to be the major intracellular compartment involved in diacytosis. However, this compartment is probably extracellular fluid trapped between cells since it is rapidly lost into the medium when the cells are either scraped or enzymatically removed from the culture dishes at 0 degree C. In addition, it very slowly undergoes both filling and emptying at 0 degree C. However, we did observe a second, much smaller, kinetic compartment (approximately 2 fl/cell) undergoing rapid diacytosis that does seem to be intracellular. A third compartment that we observed accumulates markers at a linear rate (10-20 fl cell-1 hr-1) and is not lost from cells even after incubation periods greater than 6 hr. The markers [14C]sucrose and [125I]PVP displayed very similar behavior with respect to all three compartments and yielded nearly linear long-term uptake rates, thus indicating that there is little if any absorbed component in their uptake. However, Lucifer Yellow displayed significantly higher incorporation rates and its uptake rate was strongly nonlinear, indicating its uptake in fibroblasts is predominantly adsorptive. Our observations indicate that the rate of fluid-phase endocytosis in fibroblasts is significantly less than previously reported and that any compartment involved in diacytosis is very small and turns over very rapidly. Significantly, we estimate that the constitutive internalization of clathrin-coated pits is sufficient to account for the majority of fluid-phase endocytosis and thus represents a major mechanism of membrane retrieval in these cells.     

How does cellular senescence prevent cancer?

It is widely believed that cellular senescence is a tumor suppressor mechanism; however, it has not been understood why it is advantageous for organisms to retain mutant cells is a postmitotic state rather than simply eliminating them by apoptosis. It has recently been proposed that the primary role of cellular senescence is in mitotic compartments of fixed size in which spatial considerations dictate that a deleted cell is replaced by a neighboring cell. In these situations, rather than eliminating the neoplastic clone, deletion of mutant cells can paradoxically lead to their increased turnover. If mutant cells become senescent, then the compartment is instead progressively filled by senescent cells until the mutant clone is eliminated. Since most of the genetic alterations responsible for malignancy arise in stem cells, this mechanism may have particular relevance to the stem cell niche. In this article the implications of this hypothesis are examined in detail and related to experimental results. It is further proposed here that blockage of stem cell niches by senescent stem cells may account for some of the functional alterations observed in stem cell compartments at old age. Clearly, the existence of senescent stem cells is central to the proposed hypothesis, and although there is preliminary evidence for this assertion it has yet to be proven in vivo. An experimental strategy involving double labeling of stem cells with a nucleotide label is described that can address this question.     

十三碳二元酸发酵过程菌体生长期动力学模型及其应用

介绍了由十三碳烷烃生产十三碳二元酸的发酵过程,对其中的菌体生长期的代谢过程进行了分析。提出了以CO₂释放率判断菌体生长状况的方法,据此可确定进入产酸期的最佳时间．建立了菌体生长期底物消耗及菌体生长的动力学模型,对模型参数进行了回归估值。并对菌体生长期进行了拟合。结果表明,模型的计算值和实测值吻合得较好,平均相对偏差为2．4％。利用所建模型对菌体生长期进行多种操作条件下的模拟计算,结果表明,提高蔗糖浓度及初始菌体浓度均能显著地提高菌体生长期结束时的菌体浓度。     

Modeling Facilitative Sugar Transporters: Transitions Between Single and Double Ligand Occupancy of Multiconformational Channel Models Explain Anomalous Kinetics

The four-state simple carrier model (SCM) is employed to describe ligand translocation by diverse passive membrane transporters. However, its application to systems like facilitative sugar transporters (GLUTs) is controversial: unidirectional fluxes under zero-trans and equilibrium-exchange experimental conditions fit a SCM, but flux data from infinite-cis and infinite-trans experiments appear not to fit the same SCM. More complex kinetic models have been proposed to explain this ``anomalous' behavior of GLUTs, but none of them accounts for all the experimental findings. We propose an alternative model in which GLUTs are channels subject to conformational transitions, and further assume that the results from zero-trans and equilibrium-exchange experiments as well as trans-effects corresponds to a single-occupancy channel regime, whereas the results from the infinite-cis and infinite-trans experiments correspond to a regime including higher channel occupancies. We test the plausibility of this hypothesis by studying a kinetic model of a two-site channel with two conformational states. In each state, the channel can bind the ligand from only one of the compartments. Under single-occupancy, for conditions corresponding to zero-trans and equilibrium-exchange experiments, the model behaves as a SCM capable of exhibiting trans-stimulations. For a regime including higher degrees of occupancy and infinite-cis and infinite-trans conditions, the same channel model can exhibit a behavior qualitatively similar to a SCM, albeit with kinetic parameters different from those for the single-occupancy regime. Numerical results obtained with our model are consistent with available experimental data on facilitative glucose transport across erythrocyte membranes. Hence, if GLUTs are multiconformational channels, their particular kinetic properties can result from transitions between single and double channel occupancies. Received: 12 April 1995/Revised: 28 August 1995     

Kinetics of protein and DNA synthesis studied by mathematical modelling of flow cytometric protein and DNA histograms

Abstract. Mathematical models for histograms of cellular protein content as measured by flow cytometry were developed, based on theoretical protein distributions. These were derived from the age distribution of cells and the accumulation function for cellular protein content as a function of age within the cell cycle. A model assuming an exponential age distribution and an exponential protein. accumulation function was found to give the best representation of protein histograms of exponentially growing NHIK 3025 cells. This is in good agreement with the known kinetic behaviour of such cells. By the combined use of the protein histogram model and a similar model for DNA content, and assuming linear DNA accumulation during S, the fraction of cells in S, as a function of cellular protein content, was simulated. This function showed good agreement with values of the [³H]TdR labelling index scored in cells sorted by flow cytometry from 5-channel intervals of the protein histogram. The protein and DNA histogram models were combined into a two-dimensional model for correlated protein/DNA measurements. Comparison between simulated data and experimentally derived two-dimensional protein/DNA histograms gave further support to the cell kinetic assumptions underlying the models, but also identified some minor deviations which could not be recognized in the analysis of the one-dimensional histograms.     

Mathematical model analysis of mouse epidermal cell kinetics measured by bivariate DNA/anti-bromodeoxyuridine flow cytometry and continuous [3H]-thymidine labelling

In a previous study the epidermal cell kinetics of hairless mice were investigated with bivariate DNA/anti-bromodeoxyuridine (BrdU) flow cytometry of isolated basal cells after BrdU pulse labelling. The results confirmed our previous observations of two kinetically distinct sub-populations in the G2 phase. However, the results also showed that almost all BrdU-positive cells had left S phase 6-12 h after pulse labelling, contradicting our previous assumption of a distinct, slowly cycling, major sub-population in S phase. The latter study was based on an experiment combining continuous tritiated thymidine [( 3H]TdR) labelling and cell sorting. The purpose of the present study was to use a mathematical model to analyse epidermal cell kinetics by simulating bivariate DNA/BrdU data in order to get more details about the kinetic organization and cell cycle parameter values. We also wanted to re-evaluate our assumption of slowly cycling cells in S phase. The mathematical model shows a good fit to the experimental BrdU data initiated either at 08.00 hours or 20.00 hours. Simultaneously, it was also possible to obtain a good fit to our previous continuous labelling data without including a sub-population of slowly cycling cells in S phase. This was achieved by improving the way in which the continuous [3H]TdR labelling was simulated. The presence of two distinct subpopulations in G2 phase was confirmed and a similar kinetic organization with rapidly and slowly cycling cells in G1 phase is suggested. The sizes of the slowly cycling fractions in G1 and G2 showed the same distinct circadian dependency. The model analysis indicates that a small fraction of BrdU labelled cells (3-5%) was arrested in G2 phase due to BrdU toxicity. This is insignificant compared with the total number of labelled cells and has a negligible effect on the average cell cycle data. However, it comprises 1/3 to 1/2 of the BrdU positive G2 cells after the pulse labelled cells have been distributed among the cell cycle compartments.     

Experimental analysis and modelling of in vitro HUVECs proliferation in the presence of various types of drugs

Objectives: This study focuses on experimental analysis and corresponding mathematical simulation of in vitro HUVECs (human umbilical vein endothelial cells) proliferation in the presence of various types of drugs. Materials and methods: HUVECs, once seeded in Petri dishes, were expanded to confluence. Temporal profiles of total count obtained by classic haemocytometry and cell size distribution measured using an electronic Coulter counter, are quantitatively simulated by a suitable model based on the population balance approach. Influence of drugs on cell proliferation is also properly simulated by accounting for suitable kinetic equations. Results and discussion: The models’ parameters have been determined by comparison with experimental data related to cell population expansion and cell size distribution in the absence of drugs. Inhibition constant for each type of drug has been estimated by comparing the experimental data with model results concerning temporal profiles of total cell count. The reliability of the model and its predictive capability have been tested by simulating cell size distribution for experiments performed in the presence of drugs. The proposed model will be useful in interpreting effects of selected drugs on expansion of readily available human cells.