A computational model of PKD and CERT interactions at the trans-Golgi network of mammalian cells

Benchmarking methods for computing local sensitivities in ordinary differential equation models at dynamic and steady states

Estimating parameters of dynamic models from experimental data is a challenging, and often computationally-demanding task. It requires a large number of model simulations and objective function gradient computations, if gradient-based optimization is used. In many cases, steady-state computation is a part of model simulation, either due to steady-state data or an assumption that the system is at steady state at the initial time point. Various methods are available for steady-state and gradient computation. Yet, the most efficient pair of methods (one for steady states, one for gradients) for a particular model is often not clear. In order to facilitate the selection of methods, we explore six method pairs for computing the steady state and sensitivities at steady state using six real-world problems. The method pairs involve numerical integration or Newton's method to compute the steady-state, and-for both forward and adjoint sensitivity analysis-numerical integration or a tailored method to compute the sensitivities at steady-state. Our evaluation shows that all method pairs provide accurate steady-state and gradient values, and that the two method pairs that combine numerical integration for the steady-state with a tailored method for the sensitivities at steady-state were the most robust, and amongst the most computationally-efficient. We also observed that while Newton's method for steady-state computation yields a substantial speedup compared to numerical integration, it may lead to a large number of simulation failures. Overall, our study provides a concise overview across current methods for computing sensitivities at steady state. While our study shows that there is no universally-best method pair, it also provides guidance to modelers in choosing the right methods for a problem at hand.

How has the evolution of our understanding of the compartmentalization of sphingolipid biosynthesis over the past 30 years altered our view of the evolution of the pathway?

Sphingolipids are unique among cellular lipids inasmuch as their biosynthesis is compartmentalized between the endoplasmic reticulum (ER) and the Golgi apparatus. This compartmentalization was first recognized about thirty years ago, and the current review not only updates studies on the compartmentalization of sphingolipid biosynthesis, but also discusses the ramifications of this feature for our understanding of how the pathway could have evolved. Thus, we augment some of our recent studies by inclusion of two further molecular pathways that need to be considered when analyzing the evolutionary requirements for generation of sphingolipids, namely contact sites between the ER and the Golgi apparatus, and the mechanism(s) of vesicular transport between these two organelles. Along with evolution of the individual enzymes of the pathway, their subcellular localization, and the supply of essential metabolites via the anteome, it becomes apparent that current models to describe evolution of the sphingolipid biosynthetic pathway may need substantial refinement.

Ceramide Transfer Protein (CERT): An Overlooked Molecular Player in Cancer

Sphingolipids are a class of essential lipids implicated in constructing cellular membranes and regulating nearly all cellular functions. Sphingolipid metabolic network is centered with the ceramide-sphingomyelin axis. Ceramide is well-recognized as a pro-apoptotic signal; while sphingomyelin, as the most abundant type of sphingolipids, is required for cell growth. Therefore, the balance between these two sphingolipids can be critical for cancer cell survival and functioning. Ceramide transfer protein (CERT) dictates the ratio of ceramide to sphingomyelin within the cell. It is the only lipid transfer protein that specifically delivers ceramide from the endoplasmic reticulum to the Golgi apparatus, where ceramide serves as the substrate for sphingomyelin synthesis. In the past two decades, an increasing body of evidence has suggested a critical role of CERT in cancer, but much more intensive efforts are required to draw a definite conclusion. Herein, we review all research findings of CERT, focusing on its molecular structure, cellular functions and implications in cancer. This comprehensive review of CERT will help to better understand the molecular mechanism of cancer and inspire to identify novel druggable targets.

Benchmark problems for dynamic modeling of intracellular processes

Motivation

Dynamic models are used in systems biology to study and understand cellular processes like gene regulation or signal transduction. Frequently, ordinary differential equation (ODE) models are used to model the time and dose dependency of the abundances of molecular compounds as well as interactions and translocations. A multitude of computational approaches, e.g. for parameter estimation or uncertainty analysis have been developed within recent years. However, many of these approaches lack proper testing in application settings because a comprehensive set of benchmark problems is yet missing.

Results

We present a collection of 20 benchmark problems in order to evaluate new and existing methodologies, where an ODE model with corresponding experimental data is referred to as problem. In addition to the equations of the dynamical system, the benchmark collection provides observation functions as well as assumptions about measurement noise distributions and parameters. The presented benchmark models comprise problems of different size, complexity and numerical demands. Important characteristics of the models and methodological requirements are summarized, estimated parameters are provided, and some example studies were performed for illustrating the capabilities of the presented benchmark collection.

Availability and implementation

The models are provided in several standardized formats, including an easy-to-use human readable form and machine-readable SBML files. The data is provided as Excel sheets. All files are available at https://github.com/Benchmarking-Initiative/Benchmark-Models, including step-by-step explanations and MATLAB code to process and simulate the models.

Supplementary information

Supplementary data are available at Bioinformatics online.

The Great Escape: how phosphatidylinositol 4-kinases and PI4P promote vesicle exit from the Golgi (and drive cancer)

Phosphatidylinositol 4-phosphate (PI4P) is a membrane glycerophospholipid and a major regulator of the characteristic appearance of the Golgi complex as well as its vesicular trafficking, signalling and metabolic functions. Phosphatidylinositol 4-kinases, and in particular the PI4KIIIβ isoform, act in concert with PI4P to recruit macromolecular complexes to initiate the biogenesis of trafficking vesicles for several Golgi exit routes. Dysregulation of Golgi PI4P metabolism and the PI4P protein interactome features in many cancers and is often associated with tumour progression and a poor prognosis. Increased expression of PI4P-binding proteins, such as GOLPH3 or PITPNC1, induces a malignant secretory phenotype and the release of proteins that can remodel the extracellular matrix, promote angiogenesis and enhance cell motility. Aberrant Golgi PI4P metabolism can also result in the impaired post-translational modification of proteins required for focal adhesion formation and cell–matrix interactions, thereby potentiating the development of aggressive metastatic and invasive tumours. Altered expression of the Golgi-targeted PI 4-kinases, PI4KIIIβ, PI4KIIα and PI4KIIβ, or the PI4P phosphate Sac1, can also modulate oncogenic signalling through effects on TGN-endosomal trafficking. A Golgi trafficking role for a PIP 5-kinase has been recently described, which indicates that PI4P is not the only functionally important phosphoinositide at this subcellular location. This review charts new developments in our understanding of phosphatidylinositol 4-kinase function at the Golgi and how PI4P-dependent trafficking can be deregulated in malignant disease.

Lipid-dependent coupling of secretory cargo sorting and trafficking at the trans-Golgi network

In eukaryotic cells, the trans-Golgi network (TGN) serves as a platform for secretory cargo sorting and trafficking. In recent years, it has become evident that a complex network of lipid–lipid and lipid–protein interactions contributes to these key functions. This review addresses the role of lipids at the TGN with a particular emphasis on sphingolipids and diacylglycerol. We further highlight how these lipids couple secretory cargo sorting and trafficking for spatiotemporal coordination of protein transport to the plasma membrane.

Structure, functions and regulation of CERT, a lipid-transfer protein for the delivery of ceramide at the ER-Golgi membrane contact sites

The inter-organelle transport of lipids must be regulated to ensure appropriate lipid composition of each organelle. In mammalian cells, ceramide synthesised in the endoplasmic reticulum (ER) is transported to the trans-Golgi regions, where ceramide is converted to sphingomyelin (SM) with the concomitant production of diacylglycerol. Ceramide transport protein (CERT) transports ceramide from the ER to the trans-Golgi regions at the ER-Golgi membrane contact sites (MCS). The function of CERT is down-regulated by multisite phosphorylation of a serine-repeat motif (SRM) and up-regulated by phosphorylation of serine 315 in CERT. Multisite phosphorylation of the SRM is primed by protein kinase D, which is activated by diacylglycerol. The function of CERT is regulated by a phosphorylation-dependent feedback mechanism in response to cellular requirements of SM. CERT-dependent ceramide transport is also affected by the pool of phosphatidylinositol (PtdIns)-4-phosphate (PtdIns(4)P) in the trans-Golgi regions, while the PtdIns(4)P pool is regulated by PtdIns-4-kinases and oxysterol-binding protein. The ER-Golgi MCS may serve as inter-organelle communication zones, in which many factors work in concert to serve as an extensive rheostat of SM, diacylglycerol, cholesterol and PtdIns(4)P.

The tumor suppressor protein DLC1 maintains protein kinase D activity and Golgi secretory function

Many newly synthesized cellular proteins pass through the Golgi complex from where secretory transport carriers sort them to the plasma membrane and the extracellular environment. The formation of these secretory carriers at the trans-Golgi network is promoted by the protein kinase D (PKD) family of serine/threonine kinases. Here, using mathematical modeling and experimental validation of the PKD activation and substrate phosphorylation kinetics, we reveal that the expression level of the PKD substrate deleted in liver cancer 1 (DLC1), a Rho GTPase-activating protein that is inhibited by PKD-mediated phosphorylation, determines PKD activity at the Golgi membranes. RNAi-mediated depletion of DLC1 reduced PKD activity in a Rho-Rho-associated protein kinase (ROCK)-dependent manner, impaired the exocytosis of the cargo protein horseradish peroxidase, and was associated with the accumulation of the small GTPase RAB6 on Golgi membranes, indicating a protein-trafficking defect. In summary, our findings reveal that DLC1 maintains basal activation of PKD at the Golgi and Golgi secretory activity, in part by down-regulating Rho-ROCK signaling. We propose that PKD senses cytoskeletal changes downstream of DLC1 to coordinate Rho signaling with Golgi secretory function.

Ceramide-transfer protein-mediated ceramide transfer is a structurally tunable flow-inducing mechanism with structural feed-forward loops

This paper considers two models of ceramide-transfer protein (CERT)-mediated ceramide transfer at the trans-Golgi network proposed in the literature, short distance shuttle and neck swinging, and seeks structural (parameter-free) features of the two models, which rely exclusively on the peculiar interaction network and not on specific parameter values. In particular, it is shown that both models can be seen as flow-inducing systems, where the flows between pairs of species are tuned by the concentrations of other species, and suitable external inputs can structurally regulate ceramide transfer. In the short distance shuttle model, the amount of transferred ceramide is structurally tuned by active protein kinase D (PKD), both directly and indirectly, in a coherent feed-forward loop motif. In the neck-swinging model, the amount of transferred ceramide is structurally tuned by active PI4KIIIβ, while active PKD has an ambivalent effect, due to the presence of an incoherent feed-forward loop motif that directly inhibits ceramide transfer and indirectly promotes it; the structural role of active PKD is to favour CERT mobility in the cytosol. It is also shown that the influences among key variables often have structurally determined steady-state signs, which can help falsify the models against experimental traces.

Role and Function of Sphingomyelin Biosynthesis in the Development of Cancer

Sphingomyelin (SM) biosynthesis represents a complex, finely regulated process, mostly occurring in vertebrates. It is intimately linked to lipid transport and it is ultimately carried out by two enzymes, SM synthase 1 and 2, selectively localized in the Golgi and plasma membrane. In the course of the SM biosynthetic reaction, various lipids are metabolized. Because these lipids have both structural and signaling functions, the SM biosynthetic process has the potential to affect diverse important cellular processes (such as cell proliferation, cell survival, and migration). Thus defects in SM biosynthesis might directly or indirectly impact the normal physiology of the cell and eventually of the organism. In this chapter, we will focus on evidence supporting a role for SM biosynthesis in specific cellular functions and how its dysregulation can affect neoplastic transformation.