Modeling species-specific diacylglycerol dynamics in the RAW 264.7 macrophage

A mathematical model for the determination of steady-state cardiolipin remodeling mechanisms using lipidomic data

Technical advances in lipidomic analysis have generated tremendous amounts of quantitative lipid molecular species data, whose value has not been fully explored. We describe a novel computational method to infer mechanisms of de novo lipid synthesis and remodeling from lipidomic data. We focus on the mitochondrial-specific lipid cardiolipin (CL), a polyglycerol phospholipid with four acyl chains. The lengths and degree of unsaturation of these acyl chains vary across CL molecules, and regulation of these differences is important for mitochondrial energy metabolism. We developed a novel mathematical approach to determine mechanisms controlling the steady-state distribution of acyl chain combinations in CL . We analyzed mitochondrial lipids from 18 types of steady-state samples, each with at least 3 replicates, from mouse brain, heart, lung, liver, tumor cells, and tumors grown in vitro. Using a mathematical model for the CL remodeling mechanisms and a maximum likelihood approach to infer parameters, we found that for most samples the four chain positions have an independent and identical distribution, indicating they are remodeled by the same processes. Furthermore, for most brain samples and liver, the distribution of acyl chains is well-fit by a simple linear combination of the pools of acyl chains in phosphatidylcholine (PC), phosphatidylethanolamine (PE), and phosphatidylglycerol (PG). This suggests that headgroup chemistry is the key determinant of acyl donation into CL, with chain length/saturation less important. This canonical remodeling behavior appears damaged in some tumor samples, which display a consistent excess of CL molecules having particular masses. For heart and lung, the “proportional incorporation” assumption is not adequate to explain the CL distribution, suggesting additional acyl CoA-dependent remodeling that is chain-type specific. Our findings indicate that CL remodeling processes can be described by a small set of quantitative relationships, and that bioinformatic approaches can help determine these processes from high-throughput lipidomic data.

Mathematical modelling and analysis of cellular signalling in macrophages

Cell signalling pathways play a crucial role in proper cell development and behaviour, with implications to survival, chemotaxis, proliferation, and even programmed cell death known as apoptosis. In this article, we outline a mathematical model of the G-protein signalling pathway in a particular cell line of macrophages, focusing on activation of a particular G-protein-coupled receptor, P2Y(6). The model is based on the kinetics of P2Y(6) surface receptors, inositol trisphosphate, cytosolic calcium, and differential dynamics of multiple species of diacylglycerol. Insight into the dynamics of the system is given through recently available experimental results and incorporated into the model. Mathematical analysis of the model, including establishment of global existence, uniqueness, positivity, and boundedness of solutions, and global stability of a unique steady-state solution is discussed.