A dynamic model linking cell growth to intracellular metabolism and extracellular by-product accumulation

Mathematical modeling of animal cell growth and metabolism is essential for the understanding and improvement of the production of biopharmaceuticals. Models can explain the dynamic behavior of cell growth and product formation, support the identification of the most relevant parameters for process design, and significantly reduce the number of experiments to be performed for process optimization. Few dynamic models have been established that describe both extracellular and intracellular dynamics of growth and metabolism of animal cells. In this study, a model was developed, which comprises a set of 33 ordinary differential equations to describe batch cultivations of suspension AGE1.HN.AAT cells considered for the production of α1-antitrypsin. This model combines a segregated cell growth model with a structured model of intracellular metabolism. Overall, it considers the viable cell concentration, mean cell diameter, viable cell volume, concentration of extracellular substrates, and intracellular concentrations of key metabolites from the central carbon metabolism. Furthermore, the release of metabolic by-products such as lactate and ammonium was estimated directly from the intracellular reactions. Based on the same set of parameters, this model simulates well the dynamics of four independent batch cultivations. Analysis of the simulated intracellular rates revealed at least two distinct cellular physiological states. The first physiological state was characterized by a high glycolytic rate and high lactate production. Whereas the second state was characterized by efficient adenosine triphosphate production, a low glycolytic rate, and reactions of the TCA cycle running in the reverse direction from α-ketoglutarate to citrate. Finally, we show possible applications of the model for cell line engineering and media optimization with two case studies.     

Enzyme kinetics and distinct modulation of the protein kinase N family of kinases by lipid activators and small molecule inhibitors

The PKN (protein kinase N) family of Ser/Thr protein kinases regulates a diverse set of cellular functions, such as cell migration and cytoskeletal organization. Inhibition of tumour PKN activity has been explored as an oncology therapeutic approach, with a PKN3-targeted RNAi (RNA interference)-derived therapeutic agent in Phase I clinical trials. To better understand this important family of kinases, we performed detailed enzymatic characterization, determining the kinetic mechanism and lipid sensitivity of each PKN isoform using full-length enzymes and synthetic peptide substrate. Steady-state kinetic analysis revealed that PKN1–3 follows a sequential ordered Bi–Bi kinetic mechanism, where peptide substrate binding is preceded by ATP binding. This kinetic mechanism was confirmed by additional kinetic studies for product inhibition and affinity of small molecule inhibitors. The known lipid effector, arachidonic acid, increased the catalytic efficiency of each isoform, mainly through an increase in k_cat for PKN1 and PKN2, and a decrease in peptide K_M for PKN3. In addition, a number of PKN inhibitors with various degrees of isoform selectivity, including potent (K_i<10 nM) and selective PKN3 inhibitors, were identified by testing commercial libraries of small molecule kinase inhibitors. This study provides a kinetic framework and useful chemical probes for understanding PKN biology and the discovery of isoform-selective PKN-targeted inhibitors.     

Changes in Triacylglycerol Composition in Maturing Sea Buckthorn (Hippophaë rhamnoides) Seeds

Sea buckthorn (Hippophaë rhamnoides L.) seeds on the 29th, 53rd, 80th, and 107th day after pollination were used for determining, by lipase hydrolysis, the qualitative and quantitative composition of the triacylglycerol (TAG) positional types, groups, and positional species, as well as the factor of selectivity of incorporation of unsaturated fatty acids, octadecenoic, linoleic, and linolenic, into the sn-2-position of TAGs. Until the 80th day after pollination, there was a predominant formation of triunsaturated TAGs, which included linolenic and linoleic acid residues. After the 80th day, the absolute content of these major components of total TAGs markedly decreased, and an increase in total TAG content was mainly accounted for by the rise in the level of those TAG species, which included saturated fatty acids, palmitic and stearic (monosaturated–diunsaturated and disaturated–monounsaturated), as well as in the level of sn-2-octadecenoyl species belonging to the triunsaturated and palmito–diunsaturated types of TAGs. At each maturation stage, the quantitative dynamics of separate TAG species was determined by the content of fatty acid species available for TAG formation and the factor of selectivity of these species. The decrease in the content of a certain group of triunsaturated TAGs found here seems to be caused by their metabolization during seed maturation.     

Properties of spinach chloroplast fructose-1,6-diphosphatase

Photosynthetic fructose-1,6-diphosphatase (FDPase) fractions I and II, earlier purified from spinach leaves, show a similar amino acid composition, with the exception of a higher glutamic acid content in the latter. In both fractions glutamic and aspartic acids are the main amino acids. pH activity profiles of fractions I and II are similar, with optima at 8·65–8·70, both showing a high specificity for fructose- 1,6-diphosphate. These two fractions are Mg²⁺-dependent for activity, with an Optimum Mg²⁺ concentration of 10 mM in standard conditions, which shifts to 5 mM when the MG²⁺/EDTA ratio is increased to 10; Mn²⁺ and Co²⁺ are slightly active. EDTA enhances FDPase activity slightly, with an optimum at 0·4–0·8 mM. Cysteine has no activating effect, and acts as an inhibitor above 10 mM. Both I and II have an optimum substrate concentration of 4 mM, and the substrate inhibits at concns above this value. Kinetic velocity curves are sigmoidal, with the concave zone located in the range of physiological substrate concns. (Hill coefficient 1·75 for both). This suggests a strong regulatory role of fructose-1,6-diphosphate. K_m values are 1·4 × 10⁻³ M (fraction I) and 1·1 × 10⁻³ M (fraction II). The highest activity rate occurs at 60°, in accordance with the high thermostability of both fractions; the activation energies are 14·3 kcal/mol (fraction I) and 13·0 kcal/mol (fraction II).     

A thermodynamic and crystallographic study of complexes of the highly preorganized ligand 8-hydroxyquinoline-2-carboxylic acid

The metal ion complexing properties of the ligand HQC (8-hydroxyquinoline-2-carboxylic acid) are reported. The structures of [Zn(HQCH)₂] · 3H₂O (1) and [Cd(HQCH)₂] · 3H₂O (2) were determined (HQCH = HQC with phenol protonated). Both 1 and 2 are triclinic, space group , with Z = 2. For 1 a = 7.152(3), b = 9.227(4), c = 15.629(7) Å,  = 103.978(7)°, β = 94.896(7)°, γ = 108.033(8)°, R = 0.0499. For 2 a = 7.0897(5), b = 9.1674(7), c = 16.0672(11) Å,  = 105.0240(10)°, β = 93.9910(10)°, γ = 107.1270(10)°, R = 0.0330. In 1 the Zn has a distorted octahedral coordination geometry, with Zn–N of 2.00 and 2.15 Å, and Zn–O to the protonated phenolic oxygens of 2.431 and 2.220 Å. The structure of 2 is similar, with Cd–N bonds of 2.220 and 2.228 Å, with Cd–O bonds to the protonated phenolate oxygens of 2.334 and 2.463 Å. The structures of 1 and 2, and isomorphous Ni(II) and Co(II) HQC complexes reported in the literature, show very interesting short (<2.5 Å) O–O distances in H-bonds involving the protons on the coordinated phenolates and lattice water molecules. These are discussed in relation to the possible role of short low-energy H-bonds in alcohol dehydrogenase in mediating the transfer of the hydroxyl proton of the alcohol to an adjacent serine oxygen.

The formation constants for HQC are determined by UV–Visible spectroscopy at 25 °C in 0.1 M NaClO₄ with Mg(II), Ca(II), Sr(II), Ba(II), La(III), Gd(III), Zn(II), Cd(II), Ni(II), Cu(II), and Pb(II). These show greatest stabilization with metal ions with an ionic radius above 1.0 Å. This is as would be expected from the fact that HQC forms two five-membered chelate rings on complex-formation, which favors larger metal ions. The ligand design concept of using rigid aromatic backbones in ligands to achieve high levels of preorganization, and hence the high log K values (for a tridentate ligand) and strong metal ion selectivities observed for HQC, is discussed.     

杜仲细胞悬浮培养产黄酮及其动力学研究

本文应用正交设计对杜仲细胞悬浮培养的基本培养基和植物生长物质浓度进行了筛选,并对影响杜仲细胞悬浮培养和总黄酮含量的不同因素进行了考察。结果表明,B5培养基+0.5mg/L NAA+0.6mg/L 6-BA、蔗糖30g/L、初始pH 5.0-5.5、接种量20g(FW)/L以及摇床转速110r/min为杜仲细胞悬浮培养的适宜条件。通过对杜仲悬浮细胞生长和代谢动力学的分析表明：杜仲细胞悬浮培养生长符合Logistic生长模型,最大比生长速率( m)为0.417d-1;细胞基于蔗糖的真正比生长得率(YG)与维持系数(m)分别为0.619g/g和0.0206g/(g·d-1);黄酮合成属部分生长耦联型,可用Luedeking-Piret模型进行描述。研究结果为杜仲细胞大规模悬浮培养生产天然活性成分奠定了基础。     

Thermodynamic and kinetic determinants of Thermotoga maritima cold shock protein stability: a structural and dynamic analysis

The cold shock protein (CSP) from hyperthermophile Thermotoga maritima (TmCSP) is only marginally stable (DeltaG(T(opt)) = 0.3 kcal/mol) at 353 K, the optimum environmental temperature (T(opt)) for T. maritima. In comparison, homologous CSPs from E. coli (DeltaG(T(opt)) = 2.2 kcal/mol) and B. subtilis (DeltaG(T(opt)) = 1.5 kcal/mol) are at least five times more stable at 310 K, the T(opt) for the mesophiles. Yet at the room temperature, TmCSP is more stable (DeltaG(T(R)) = 4.7 kcal/mol) than its homologues (DeltaG(T(R)) = 3.0 kcal/mol for E. coli CSP and DeltaG(T(R)) = 2.1 kcal/mol for B. subtilis CSP). This unique observation suggests that kinetic, rather than thermodynamic, barriers toward unfolding might help TmCSP native structure at high temperatures. Consistently, the unfolding rate of TmCSP is considerably slower than its homologues. High temperature (600 K) complete unfolding molecular dynamics (MD) simulations of TmCSP support our hypothesis and reveal an unfolding scheme unique to TmCSP. For all the studied homologues of TmCSP, the unfolding process first starts at the C-terminal region and N-terminal region unfolds in the end. But for TmCSP, both the terminals resist unfolding for consistently longer simulation times and, in the end, unfold simultaneously. In TmCSP, the C-terminal region is better fortified and has better interactions with the N-terminal region due to the charged residues, R2, E47, E49, H61, K63, and E66, being in spatial vicinity. The electrostatic interactions among these residues are unique to TmCSP. Consistently, the room temperature MD simulations show that TmCSP is more rigid at its N- and C-termini as compared to its homologues from E. coli, B. subtilis, and B. caldolyticus.     

Kinetics of adsorption of extended ligands on DNA at small fillings

In the present work, the adsorption kinetics of extended ligands on DNA duplexes at small fillings when molecules of DNA duplexes are on the underlayer within diffusion layer has been investigated. Both diffusion of ligands in solution (diffusion stage) and adsorption of ligands (kinetic stage) are taken into consideration at adsorption of ligands on DNA duplexes. Nonlinear system of differential equations describing adsorption of ligands where not only diffusion stage but also kinetic stage is taken into account, is obtained, moreover the equations allow localizing duplexes in arbitrary place within diffusion layer. Numeric solution of the equations makes possible to investigate the filling kinetics of DNA duplexes by ligands depending on parameters controlling adsorption process. It has been shown that depending on relation between adsorption parameters different kinetic regimes of adsorption – kinetic, complex, and diffusion regimes may be realized.     

Kinetic locking-on and auxiliary tactics for bioaffinity purification of NADP+-dependent dehydrogenases using N6-linked immobilized NADP+ derivatives: Studies with mammalian and microbial glutamate dehydrogenases

This study is concerned with the development and application of kinetic locking-on and auxiliary tactics for bioaffinity purification of NADP(+)-dependent dehydrogenases, specifically (1) the synthesis and characterization of highly substituted N(6)-linked immobilized NADP(+) derivatives using a rapid solid-phase modular approach; (2) the evaluation of the N(6)-linked immobilized NADP(+) derivatives for use with the kinetic locking-on strategy for bioaffinity purification of NADP(+)-dependent dehydrogenases: Model bioaffinity chromatographic studies with glutamate dehydrogenase from bovine liver (GDH with dual cofactor specificity, EC 1.4.1.3) and glutamate dehydrogenase from Candida utilis (GDH which is NADP(+)-specific, EC 1.4.1.4); (3) the selection of an effective "stripping ligand" for NADP(+)-dehydrogenase bioaffinity purifications using N(6)-linked immobilized NADP(+) derivatives in the locking-on mode; and (4) the application of the developed bioaffinity chromatographic system to the purification of C. utilis GDH from a crude cellular extract.Results confirm that the newly developed N(6)-linked immobilized NADP(+) derivatives are suitable for the one-step bioaffinity purification of NADP(+)-dependent GDH provided that they are used in the locking-on mode, steps are taken to inhibit alkaline phosphatase activity in crude cellular extracts, and 2',5'-ADP is used as the stripping ligand during chromatography. The general principles described here are supported by a specific sample enzyme purification; the purification of C. utilis GDH to electrophoretic homogeneity in a single bioaffinity chromatographic step (specific activity, 9.12 micromol/min/mg; purification factor, 83.7; yield 88%). The potential for development of analogous bioaffinity systems for other NADP(+)-dependent dehydrogenases is also discussed.