Exploring the Altered Dynamics of Mammalian Central Carbon Metabolic Pathway in Cancer Cells: A Classical Control Theoretic Approach

Utility of constraints reflecting system stability on analyses for biological models

Simulating complex biological models consisting of multiple ordinary differential equations can aid in the prediction of the pharmacological/biological responses; however, they are often hampered by the availability of reliable kinetic parameters. In the present study, we aimed to discover the properties of behaviors without determining an optimal combination of kinetic parameter values (parameter set). The key idea was to collect as many parameter sets as possible. Given that many systems are biologically stable and resilient (BSR), we focused on the dynamics around the steady state and formulated objective functions for BSR by partial linear approximation of the focused region. Using the objective functions and modified global cluster Newton method, we developed an algorithm for a thorough exploration of the allowable parameter space for biological systems (TEAPS). We first applied TEAPS to the NF-κB signaling model. This system shows a damped oscillation after stimulation and seems to fit the BSR constraint. By applying TEAPS, we found several directions in parameter space which stringently determines the BSR property. In such directions, the experimentally fitted parameter values were included in the range of the obtained parameter sets. The arachidonic acid metabolic pathway model was used as a model related to pharmacological responses. The pharmacological effects of nonsteroidal anti-inflammatory drugs were simulated using the parameter sets obtained by TEAPS. The structural properties of the system were partly extracted by analyzing the distribution of the obtained parameter sets. In addition, the simulations showed inter-drug differences in prostacyclin to thromboxane A2 ratio such that aspirin treatment tends to increase the ratio, while rofecoxib treatment tends to decrease it. These trends are comparable to the clinical observations. These results on real biological models suggest that the parameter sets satisfying the BSR condition can help in finding biologically plausible parameter sets and understanding the properties of biological systems.

We propose a new method to analyze the properties of biological dynamic models, which we named TEAPS (Thorough Exploration of Allowable Parameter Space). TEAPS can thoroughly determine combinations of parameter values for ordinary differential equations with which an initial state in a certain range converges to a particular fixed point. This stable and resilient behavior is a characteristic shared with many biological systems, including metabolic systems and intracellular signaling systems. Therefore, this thorough search outlined the possible parameter space as biological systems for target models, which helps to understand the system constraints when the target systems behave dynamically. The obtained parameter space can be used as an initial space for parameter tuning. For models that include a large number of parameters, the parameter space to be searched in the parameter tuning process is too large; therefore, narrowing down the space by TEAPS potentially contributes to the analysis of the dynamics of complicated biological models. Thus, our approach can partly overcome the current problem in parameter tuning and can advance the computational dynamic analyses of biological systems.

A control theoretic three timescale model for analyzing energy management in mammalian cancer cells

Interaction among different pathways, such as metabolic, signaling and gene regulatory networks, of cellular system is responsible to maintain homeostasis in a mammalian cell. Malfunctioning of this cooperation may lead to many complex diseases, such as cancer and type 2 diabetes. Timescale differences among these pathways make their integration a daunting task. Metabolic, signaling and gene regulatory networks have three different timescales, such as, ultrafast, fast and slow respectively. The article deals with this problem by developing a support vector regression (SVR) based three timescale model with the application of genetic algorithm based nonlinear controller. The proposed model can successfully capture the nonlinear transient dynamics and regulations of such integrated biochemical pathway under consideration. Besides, the model is quite capable of predicting the effects of certain drug targets for many types of complex diseases. Here, energy and cell proliferation management of mammalian cancer cells have been explored and analyzed with the help of the proposed novel approach. Previous investigations including in silico/in vivo/in vitro experiments have validated the results (the regulations of glucose transporter 1 (glut1), hexokinase (HK), and hypoxia-inducible factor-1 α (HIF-1 α ) among others, and the switching of pyruvate kinase (M2 isoform) between dimer and tetramer) generated by this model proving its effectiveness. Subsequently, the model predicts the effects of six selected drug targets, such as, the deactivation of transketolase and glucose-6-phosphate isomerase among others, in the case of mammalian malignant cells in terms of growth, proliferation, fermentation, and energy supply in the form of adenosine triphosphate (ATP).

Metabolic pathway engineering: Perspectives and applications

BACKGROUND

Metabolic engineering aims at contriving microbes as biocatalysts for enhanced and cost-effective production of countless secondary metabolites. These secondary metabolites can be treated as the resources of industrial chemicals, pharmaceuticals and fuels. Plants are also crucial targets for metabolic engineers to produce necessary secondary metabolites. Metabolic engineering of both microorganism and plants also contributes towards drug discovery. In order to implement advanced metabolic engineering techniques efficiently, metabolic engineers should have detailed knowledge about cell physiology and metabolism. Principle behind methodologies: Genome-scale mathematical models of integrated metabolic, signal transduction, gene regulatory and protein-protein interaction networks along with experimental validation can provide such knowledge in this context. Incorporation of omics data into these models is crucial in the case of drug discovery. Inverse metabolic engineering and metabolic control analysis (MCA) can help in developing such models. Artificial intelligence methodology can also be applied for efficient and accurate metabolic engineering.

CONCLUSION

In this review, we discuss, at the beginning, the perspectives of metabolic engineering and its application on microorganism and plant leading to drug discovery. At the end, we elaborate why inverse metabolic engineering and MCA are closely related to modern metabolic engineering. In addition, some crucial steps ensuring efficient and optimal metabolic engineering strategies have been discussed. Moreover, we explore the use of genomics data for the activation of silent metabolic clusters and how it can be integrated with metabolic engineering. Finally, we exhibit a few applications of artificial intelligence to metabolic engineering.

Catestatin improves insulin sensitivity by attenuating endoplasmic reticulum stress: In vivo and in silico validation

•

An endogenous peptide catestatin alleviates obesity-induced ER stress.

•

Alleviation of ER stress by catestatin improves insulin sensitivity.

•

PID controller based model of ER stress is supported by experimental findings.

•

It predicts AKT phosphorylation achieves insulin sensitivity overcoming ER stress.

Obesity is characterized by a state of chronic, unresolved inflammation in insulin-targeted tissues. Obesity-induced inflammation causes accumulation of proinflammatory macrophages in adipose tissue and liver. Proinflammatory cytokines released from tissue macrophages inhibits insulin sensitivity. Obesity also leads to inflammation-induced endoplasmic reticulum (ER) stress and insulin resistance. In this scenario, based on the data (specifically patterns) generated by our in vivo experiments on both diet-induced obese (DIO) and normal chow diet (NCD) mice, we developed an in silico state space model to integrate ER stress and insulin signaling pathways. Computational results successfully followed the experimental results for both DIO and NCD conditions. Chromogranin A (CgA) peptide catestatin (CST: hCgA352-372) improves obesity-induced hepatic insulin resistance by reducing inflammation and inhibiting proinflammatory macrophage infiltration. We reasoned that the anti-inflammatory effects of CST would alleviate ER stress. CST decreased obesity-induced ER dilation in hepatocytes and macrophages. On application of Proportional-Integral-Derivative (PID) controllers on the in silico model, we checked whether the reduction of phosphorylated PERK resulting in attenuation of ER stress, resembling CST effect, could enhance insulin sensitivity. The simulation results clearly pointed out that CST not only decreased ER stress but also enhanced insulin sensitivity in mammalian cells. In vivo experiment validated the simulation results by depicting that CST caused decrease in phosphorylation of UPR signaling molecules and increased phosphorylation of insulin signaling molecules. Besides simulation results predicted that enhancement of AKT phosphorylation helps in both overcoming ER stress and achieving insulin sensitivity. These effects of CST were verified in hepatocyte culture model.

Succinate aggravates NAFLD progression to liver cancer on the onset of obesity: An in silico model

The incidence and prevalence of nonalcoholic fatty liver disease (NAFLD) have been increasing to epidemic proportions around the world. NAFLD, a chronic liver disease that affects the nondrinkers, is mainly associated with steatohepatitis and cirrhosis. The progression of NAFLD associated with obesity increases the risk of liver cancer, a disease with poor outcomes and limited therapeutic options. In order to investigate the underlying cellular dynamics leading to NAFLD progression towards cancer on the onset of obesity, we have integrated human hepatocyte pathway with hypoxia-inducible factor1- <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"><mml:mi>α</mml:mi></mml:math> (HIF1- <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"><mml:mi>α</mml:mi></mml:math> ) signaling pathway using state space model based on classical control theory. Modified Michaelis-Menten equation and mass action law have been used to define flux vectors of the proposed model. We have incorporated feedback inhibition/activation and allosteric effects into the simulink-based model. The values of kinetic constants have been taken from the literature. It is found that on the onset of obesity, HIF1- <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"><mml:mi>α</mml:mi></mml:math> -induced proteins stabilize approximately 62 times that in the case of a normal cell. Consequently, the HIF1- <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"><mml:mi>α</mml:mi></mml:math> -induced proteins enhance the enzymatic activities of hexokinase (HK), phosphofructo kinase (PFK), lactate dehydrogenase (LDH), and pyruvate dehydrogenase (PDH), which induce Warburg effect promoting an environment suitable for cancer cells.

Computational properties of mitochondria in T cell activation and fate

In this article, we review how mitochondrial Ca²⁺ transport (mitochondrial Ca²⁺ uptake and Na⁺/Ca²⁺ exchange) is involved in T cell biology, including activation and differentiation through shaping cellular Ca²⁺ signals. Based on recent observations, we propose that the Ca²⁺ crosstalk between mitochondria, endoplasmic reticulum and cytoplasm may form a proportional–integral–derivative (PID) controller. This PID mechanism (which is well known in engineering) could be responsible for computing cellular decisions. In addition, we point out the importance of analogue and digital signal processing in T cell life and implication of mitochondrial Ca²⁺ transport in this process.

A fuzzy logic controller based approach to model the switching mechanism of the mammalian central carbon metabolic pathway in normal and cancer cells

Dynamics of large nonlinear complex systems, like metabolic networks, depend on several parameters. A metabolic pathway may switch to another pathway in accordance with the current state of parameters in both normal and cancer cells. Here, most of the parameter values are unknown to us. A fuzzy logic controller (FLC) has been developed here for the purpose of modeling metabolic networks by approximating the reasons for the behaviour of a system and applying expert knowledge to track switching between metabolic pathways. The simulation results can track the switching between glycolysis and gluconeogenesis, as well as glycolysis and pentose phosphate pathways (PPP) in normal cells. Unlike normal cells, pyruvate kinase (M2 isoform) (PKM2) switches alternatively between its two oligomeric forms, i.e. an active tetramer and a relatively low activity dimer, in cancer cells. Besides, there is a coordination among PKM2 switching and enzymes catalyzing PPP. These phenomena help cancer cells to maintain their high energy demand and macromolecular synthesis. However, the reduction of initial adenosine triphosphate (ATP) to a very low concentration, decreasing initial glucose uptake, destroying coordination between glycolysis and PPP, and replacement of PKM2 by its relatively inactive oligomeric form (dimer) or inhibition of the translation of PKM2 may destabilize the mutated control mechanism of the mammalian central carbon metabolic (CCM) pathway in cancer cells. The performance of the model is compared appropriately with some existing ones.