1
|
Schindler DJ, Clarke J, Barahona M. Multiscale mobility patterns and the restriction of human movement. ROYAL SOCIETY OPEN SCIENCE 2023; 10:230405. [PMID: 37830024 PMCID: PMC10565406 DOI: 10.1098/rsos.230405] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 04/11/2023] [Accepted: 09/18/2023] [Indexed: 10/14/2023]
Abstract
From the perspective of human mobility, the COVID-19 pandemic constituted a natural experiment of enormous reach in space and time. Here, we analyse the inherent multiple scales of human mobility using Facebook Movement maps collected before and during the first UK lockdown. Firstly, we obtain the pre-lockdown UK mobility graph and employ multiscale community detection to extract, in an unsupervised manner, a set of robust partitions into flow communities at different levels of coarseness. The partitions so obtained capture intrinsic mobility scales with better coverage than nomenclature of territorial units for statistics (NUTS) regions, which suffer from mismatches between human mobility and administrative divisions. Furthermore, the flow communities in the fine-scale partition not only match well the UK travel to work areas but also capture mobility patterns beyond commuting to work. We also examine the evolution of mobility under lockdown and show that mobility first reverted towards fine-scale flow communities already found in the pre-lockdown data, and then expanded back towards coarser flow communities as restrictions were lifted. The improved coverage induced by lockdown is well captured by a linear decay shock model, which allows us to quantify regional differences in both the strength of the effect and the recovery time from the lockdown shock.
Collapse
Affiliation(s)
| | - Jonathan Clarke
- Department of Mathematics, Imperial College London, London SW7 2BX, UK
| | - Mauricio Barahona
- Department of Mathematics, Imperial College London, London SW7 2BX, UK
| |
Collapse
|
2
|
Rombouts J, Verplaetse S, Gelens L. The ups and downs of biological oscillators: a comparison of time-delayed negative feedback mechanisms. J R Soc Interface 2023; 20:20230123. [PMID: 37376871 PMCID: PMC10300510 DOI: 10.1098/rsif.2023.0123] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/06/2023] [Accepted: 05/31/2023] [Indexed: 06/29/2023] Open
Abstract
Many biochemical oscillators are driven by the periodic rise and fall of protein concentrations or activities. A negative feedback loop underlies such oscillations. The feedback can act on different parts of the biochemical network. Here, we mathematically compare time-delay models where the feedback affects production and degradation. We show a mathematical connection between the linear stability of the two models, and derive how both mechanisms impose different constraints on the production and degradation rates that allow oscillations. We show how oscillations are affected by the inclusion of a distributed delay, of double regulation (acting on production and degradation) and of enzymatic degradation.
Collapse
Affiliation(s)
- Jan Rombouts
- Cell Biology and Biophysics Unit, European Molecular Biology Laboratory (EMBL), Heidelberg, Germany
- Developmental Biology Unit, European Molecular Biology Laboratory (EMBL), Heidelberg, Germany
- Department of Cellular and Molecular Medicine, KU Leuven, Belgium
| | - Sarah Verplaetse
- Department of Cellular and Molecular Medicine, KU Leuven, Belgium
| | - Lendert Gelens
- Department of Cellular and Molecular Medicine, KU Leuven, Belgium
| |
Collapse
|
3
|
Kim DW, Hong H, Kim JK. Systematic inference identifies a major source of heterogeneity in cell signaling dynamics: The rate-limiting step number. SCIENCE ADVANCES 2022; 8:eabl4598. [PMID: 35302852 PMCID: PMC8932658 DOI: 10.1126/sciadv.abl4598] [Citation(s) in RCA: 11] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 07/14/2021] [Accepted: 01/26/2022] [Indexed: 06/14/2023]
Abstract
Identifying the sources of cell-to-cell variability in signaling dynamics is essential to understand drug response variability and develop effective therapeutics. However, it is challenging because not all signaling intermediate reactions can be experimentally measured simultaneously. This can be overcome by replacing them with a single random time delay, but the resulting process is non-Markovian, making it difficult to infer cell-to-cell heterogeneity in reaction rates and time delays. To address this, we developed an efficient and scalable moment-based Bayesian inference method (MBI) with a user-friendly computational package that infers cell-to-cell heterogeneity in the non-Markovian signaling process. We applied MBI to single-cell expression profiles from promoters responding to antibiotics and discovered a major source of cell-to-cell variability in antibiotic stress response: the number of rate-limiting steps in signaling cascades. This knowledge can help identify effective therapies that destroy all pathogenic or cancer cells, and the approach can be applied to precision medicine.
Collapse
Affiliation(s)
- Dae Wook Kim
- Department of Mathematical Sciences, Korea Advanced Institute of Science and Technology, Daejeon 34141, Republic of Korea
- Biomedical Mathematics Group, Institute for Basic Science, Daejeon 34126, Republic of Korea
| | - Hyukpyo Hong
- Department of Mathematical Sciences, Korea Advanced Institute of Science and Technology, Daejeon 34141, Republic of Korea
- Biomedical Mathematics Group, Institute for Basic Science, Daejeon 34126, Republic of Korea
| | - Jae Kyoung Kim
- Department of Mathematical Sciences, Korea Advanced Institute of Science and Technology, Daejeon 34141, Republic of Korea
- Biomedical Mathematics Group, Institute for Basic Science, Daejeon 34126, Republic of Korea
| |
Collapse
|
4
|
Yesilkanal AE, Yang D, Valdespino A, Tiwari P, Sabino AU, Nguyen LC, Lee J, Xie XH, Sun S, Dann C, Robinson-Mailman L, Steinberg E, Stuhlmiller T, Frankenberger C, Goldsmith E, Johnson GL, Ramos AF, Rosner MR. Limited inhibition of multiple nodes in a driver network blocks metastasis. eLife 2021; 10:e59696. [PMID: 33973518 PMCID: PMC8128439 DOI: 10.7554/elife.59696] [Citation(s) in RCA: 17] [Impact Index Per Article: 5.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/05/2020] [Accepted: 04/29/2021] [Indexed: 12/14/2022] Open
Abstract
Metastasis suppression by high-dose, multi-drug targeting is unsuccessful due to network heterogeneity and compensatory network activation. Here, we show that targeting driver network signaling capacity by limited inhibition of core pathways is a more effective anti-metastatic strategy. This principle underlies the action of a physiological metastasis suppressor, Raf Kinase Inhibitory Protein (RKIP), that moderately decreases stress-regulated MAP kinase network activity, reducing output to transcription factors such as pro-metastastic BACH1 and motility-related target genes. We developed a low-dose four-drug mimic that blocks metastatic colonization in mouse breast cancer models and increases survival. Experiments and network flow modeling show limited inhibition of multiple pathways is required to overcome variation in MAPK network topology and suppress signaling output across heterogeneous tumor cells. Restricting inhibition of individual kinases dissipates surplus signal, preventing threshold activation of compensatory kinase networks. This low-dose multi-drug approach to decrease signaling capacity of driver networks represents a transformative, clinically relevant strategy for anti-metastatic treatment.
Collapse
Affiliation(s)
- Ali Ekrem Yesilkanal
- Ben May Department for Cancer Research, University of ChicagoChicagoUnited States
| | - Dongbo Yang
- Ben May Department for Cancer Research, University of ChicagoChicagoUnited States
| | - Andrea Valdespino
- Ben May Department for Cancer Research, University of ChicagoChicagoUnited States
| | - Payal Tiwari
- Ben May Department for Cancer Research, University of ChicagoChicagoUnited States
| | - Alan U Sabino
- Instituto do Câncer do Estado de São Paulo, Faculdade de Medicina and Escola de Artes, Ciências e Humanidades; University of São PauloSão PauloBrazil
| | - Long Chi Nguyen
- Ben May Department for Cancer Research, University of ChicagoChicagoUnited States
| | - Jiyoung Lee
- Ben May Department for Cancer Research, University of ChicagoChicagoUnited States
| | - Xiao-He Xie
- Ben May Department for Cancer Research, University of ChicagoChicagoUnited States
| | - Siqi Sun
- Ben May Department for Cancer Research, University of ChicagoChicagoUnited States
| | - Christopher Dann
- Ben May Department for Cancer Research, University of ChicagoChicagoUnited States
| | | | - Ethan Steinberg
- Ben May Department for Cancer Research, University of ChicagoChicagoUnited States
| | - Timothy Stuhlmiller
- Department of Pharmacology, University of North Carolina at Chapel HillChapel HillUnited States
| | - Casey Frankenberger
- Ben May Department for Cancer Research, University of ChicagoChicagoUnited States
| | | | - Gary L Johnson
- Department of Pharmacology, University of North Carolina at Chapel HillChapel HillUnited States
| | - Alexandre F Ramos
- Instituto do Câncer do Estado de São Paulo, Faculdade de Medicina and Escola de Artes, Ciências e Humanidades; University of São PauloSão PauloBrazil
| | - Marsha R Rosner
- Ben May Department for Cancer Research, University of ChicagoChicagoUnited States
| |
Collapse
|
5
|
Jiang Y, Hao N. Memorizing environmental signals through feedback and feedforward loops. Curr Opin Cell Biol 2021; 69:96-102. [PMID: 33549848 PMCID: PMC8058236 DOI: 10.1016/j.ceb.2020.11.008] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/30/2020] [Revised: 11/25/2020] [Accepted: 11/26/2020] [Indexed: 12/12/2022]
Abstract
Cells in diverse organisms can store the information of previous environmental conditions for long periods of time. This form of cellular memory adjusts the cell's responses to future challenges, providing fitness advantages in fluctuating environments. Many biological functions, including cellular memory, are mediated by specific recurring patterns of interactions among proteins and genes, known as 'network motifs.' In this review, we focus on three well-characterized network motifs - negative feedback loops, positive feedback loops, and feedforward loops, which underlie different types of cellular memories. We describe the latest studies identifying these motifs in various molecular processes and discuss how the topologies and dynamics of these motifs can enable memory encoding and storage.
Collapse
Affiliation(s)
- Yanfei Jiang
- Section of Molecular Biology, Division of Biological Sciences, University of California San Diego, 9500 Gilman Drive, La Jolla, CA, 92093, USA
| | - Nan Hao
- Section of Molecular Biology, Division of Biological Sciences, University of California San Diego, 9500 Gilman Drive, La Jolla, CA, 92093, USA.
| |
Collapse
|
6
|
Molecular switch architecture determines response properties of signaling pathways. Proc Natl Acad Sci U S A 2021; 118:2013401118. [PMID: 33688042 DOI: 10.1073/pnas.2013401118] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022] Open
Abstract
Many intracellular signaling pathways are composed of molecular switches, proteins that transition between two states-on and off Typically, signaling is initiated when an external stimulus activates its cognate receptor that, in turn, causes downstream switches to transition from off to on using one of the following mechanisms: activation, in which the transition rate from the off state to the on state increases; derepression, in which the transition rate from the on state to the off state decreases; and concerted, in which activation and derepression operate simultaneously. We use mathematical modeling to compare these signaling mechanisms in terms of their dose-response curves, response times, and abilities to process upstream fluctuations. Our analysis elucidates several operating principles for molecular switches. First, activation increases the sensitivity of the pathway, whereas derepression decreases sensitivity. Second, activation generates response times that decrease with signal strength, whereas derepression causes response times to increase with signal strength. These opposing features allow the concerted mechanism to not only show dose-response alignment, but also to decouple the response time from stimulus strength. However, these potentially beneficial properties come at the expense of increased susceptibility to upstream fluctuations. We demonstrate that these operating principles also hold when the models are extended to include additional features, such as receptor removal, kinetic proofreading, and cascades of switches. In total, we show how the architecture of molecular switches govern their response properties. We also discuss the biological implications of our findings.
Collapse
|
7
|
Exact solutions and equi-dosing regimen regions for multi-dose pharmacokinetics models with transit compartments. J Pharmacokinet Pharmacodyn 2020; 48:99-131. [PMID: 33040255 PMCID: PMC7979639 DOI: 10.1007/s10928-020-09719-8] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/08/2020] [Accepted: 09/17/2020] [Indexed: 12/02/2022]
Abstract
Compartmental models which yield linear ordinary differential equations (ODEs) provide common tools for pharmacokinetics (PK) analysis, with exact solutions for drug levels or concentrations readily obtainable for low-dimensional compartment models. Exact solutions enable valuable insights and further analysis of these systems. Transit compartment models are a popular semi-mechanistic approach for generalising simple PK models to allow for delayed kinetics, but computing exact solutions for multi-dosing inputs to transit compartment systems leading to different final compartments is nontrivial. Here, we find exact solutions for drug levels as functions of time throughout a linear transit compartment cascade followed by an absorption compartment and a central blood compartment, for the general case of n transit compartments and M equi-bolus doses to the first compartment. We further show the utility of exact solutions to PK ODE models in finding constraints on equi-dosing regimen parameters imposed by a prescribed therapeutic range. This leads to the construction of equi-dosing regimen regions (EDRRs), providing new, novel visualisations which summarise the safe and effective dosing parameter space. EDRRs are computed for classical and transit compartment models with two- and three-dimensional parameter spaces, and are proposed as useful graphical tools for informing drug dosing regimen design.
Collapse
|
8
|
It's about time: Analysing simplifying assumptions for modelling multi-step pathways in systems biology. PLoS Comput Biol 2020; 16:e1007982. [PMID: 32598362 PMCID: PMC7351226 DOI: 10.1371/journal.pcbi.1007982] [Citation(s) in RCA: 14] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/26/2019] [Revised: 07/10/2020] [Accepted: 05/27/2020] [Indexed: 11/19/2022] Open
Abstract
Thoughtful use of simplifying assumptions is crucial to make systems biology models tractable while still representative of the underlying biology. A useful simplification can elucidate the core dynamics of a system. A poorly chosen assumption can, however, either render a model too complicated for making conclusions or it can prevent an otherwise accurate model from describing experimentally observed dynamics. Here, we perform a computational investigation of sequential multi-step pathway models that contain fewer pathway steps than the system they are designed to emulate. We demonstrate when such models will fail to reproduce data and how detrimental truncation of a pathway leads to detectable signatures in model dynamics and its optimised parameters. An alternative assumption is suggested for simplifying such pathways. Rather than assuming a truncated number of pathway steps, we propose to use the assumption that the rates of information propagation along the pathway is homogeneous and, instead, letting the length of the pathway be a free parameter. We first focus on linear pathways that are sequential and have first-order kinetics, and we show how this assumption results in a three-parameter model that consistently outperforms its truncated rival and a delay differential equation alternative in recapitulating observed dynamics. We then show how the proposed assumption allows for similarly terse and effective models of non-linear pathways. Our results provide a foundation for well-informed decision making during model simplifications.
Collapse
|
9
|
Li J, Xie Y, Zhang C, Wang J, Wu Y, Yang Y, Xie Y, Lv Z. A network-based analysis for mining the risk pathways in glioblastoma. Oncol Lett 2019; 18:2712-2717. [PMID: 31402957 PMCID: PMC6676740 DOI: 10.3892/ol.2019.10598] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/01/2018] [Accepted: 06/13/2019] [Indexed: 11/05/2022] Open
Abstract
The most malignant type of brain tumour is glioblastoma multiforme (GBM). Patients with GBM often have a poor prognosis, as a result of incomplete or inaccurate diagnoses. Regulatory pathways have been demonstrated to serve important roles in complex human diseases. Therefore, deciphering these risk pathways may shed light on the molecular mechanisms underlying GBM progression. In the present study, differentially expressed genes and microRNAs (miRNAs) in a publicly available database were identified between normal and tumour samples. To determine the pathophysiology and molecular mechanisms underlying GBM, integrated network analysis was performed to mine GBM-specific risk pathways. Specifically, a GBM-specific regulatory network was constructed that integrated manually curated GBM-associated transcription and post-transcriptional data resources, including transcription factors and miRNAs. A total of 1,827 differentially expressed genes and 30 miRNAs were identified. The differentially expressed genes were significantly enriched in a number of immune response-associated functions. Based on the GBM-specific regulatory network, 15 risk regulatory pathways containing not only known regulators, but also potential novel targets that might be involved in tumourigenesis were identified. Network analysis provides a strategy for leveraging genomic data to identify potential oncogenic pathways and molecular targets for GBM.
Collapse
Affiliation(s)
- Jing Li
- Department of Hepatobiliary Surgery, Affiliated Hospital of Southwest Medical University, Luzhou, Sichuan 646000, P.R. China
| | - Yujie Xie
- Department of Rehabilitation Medicine, Affiliated Hospital of Southwest Medical University, Luzhou, Sichuan 646000, P.R. China
| | - Chi Zhang
- Department of Rehabilitation Medicine, Affiliated Hospital of Southwest Medical University, Luzhou, Sichuan 646000, P.R. China
| | - Jianxiong Wang
- Department of Rehabilitation Medicine, Affiliated Hospital of Southwest Medical University, Luzhou, Sichuan 646000, P.R. China
| | - Yong Wu
- Department of Neurology, Affiliated Hospital of Southwest Medical University, Luzhou, Sichuan 646000, P.R. China
| | - Yuan Yang
- Department of Neurology, Affiliated Hospital of Southwest Medical University, Luzhou, Sichuan 646000, P.R. China
| | - Yang Xie
- Department of Neurology, Affiliated Hospital of Southwest Medical University, Luzhou, Sichuan 646000, P.R. China
| | - Zhiyu Lv
- Department of Neurology, Affiliated Hospital of Southwest Medical University, Luzhou, Sichuan 646000, P.R. China
| |
Collapse
|
10
|
Seigal A, Beguerisse-Díaz M, Schoeberl B, Niepel M, Harrington HA. Tensor clustering with algebraic constraints gives interpretable groups of crosstalk mechanisms in breast cancer. J R Soc Interface 2019; 16:20180661. [PMID: 30958184 PMCID: PMC6408352 DOI: 10.1098/rsif.2018.0661] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022] Open
Abstract
We introduce a tensor-based clustering method to extract sparse, low-dimensional structure from high-dimensional, multi-indexed datasets. This framework is designed to enable detection of clusters of data in the presence of structural requirements which we encode as algebraic constraints in a linear program. Our clustering method is general and can be tailored to a variety of applications in science and industry. We illustrate our method on a collection of experiments measuring the response of genetically diverse breast cancer cell lines to an array of ligands. Each experiment consists of a cell line–ligand combination, and contains time-course measurements of the early signalling kinases MAPK and AKT at two different ligand dose levels. By imposing appropriate structural constraints and respecting the multi-indexed structure of the data, the analysis of clusters can be optimized for biological interpretation and therapeutic understanding. We then perform a systematic, large-scale exploration of mechanistic models of MAPK–AKT crosstalk for each cluster. This analysis allows us to quantify the heterogeneity of breast cancer cell subtypes, and leads to hypotheses about the signalling mechanisms that mediate the response of the cell lines to ligands.
Collapse
Affiliation(s)
- Anna Seigal
- 1 Department of Mathematics, University of California , Berkeley, CA 94702 , USA
| | | | - Birgit Schoeberl
- 3 Novartis Institutes for BioMedical Research , Cambridge, MA 02139 , USA
| | - Mario Niepel
- 4 Ribon Therapeutics , Lexington, MA 02421 , USA
| | | |
Collapse
|
11
|
Giese W, Milicic G, Schröder A, Klipp E. Spatial modeling of the membrane-cytosolic interface in protein kinase signal transduction. PLoS Comput Biol 2018; 14:e1006075. [PMID: 29630597 PMCID: PMC5908195 DOI: 10.1371/journal.pcbi.1006075] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/10/2017] [Revised: 04/19/2018] [Accepted: 03/07/2018] [Indexed: 02/06/2023] Open
Abstract
The spatial architecture of signaling pathways and the interaction with cell size and morphology are complex, but little understood. With the advances of single cell imaging and single cell biology, it becomes crucial to understand intracellular processes in time and space. Activation of cell surface receptors often triggers a signaling cascade including the activation of membrane-attached and cytosolic signaling components, which eventually transmit the signal to the cell nucleus. Signaling proteins can form steep gradients in the cytosol, which cause strong cell size dependence. We show that the kinetics at the membrane-cytosolic interface and the ratio of cell membrane area to the enclosed cytosolic volume change the behavior of signaling cascades significantly. We suggest an estimate of average concentration for arbitrary cell shapes depending on the cell volume and cell surface area. The normalized variance, known from image analysis, is suggested as an alternative measure to quantify the deviation from the average concentration. A mathematical analysis of signal transduction in time and space is presented, providing analytical solutions for different spatial arrangements of linear signaling cascades. Quantification of signaling time scales reveals that signal propagation is faster at the membrane than at the nucleus, while this time difference decreases with the number of signaling components in the cytosol. Our investigations are complemented by numerical simulations of non-linear cascades with feedback and asymmetric cell shapes. We conclude that intracellular signal propagation is highly dependent on cell geometry and, thereby, conveys information on cell size and shape to the nucleus.
Collapse
Affiliation(s)
- Wolfgang Giese
- Mathematical Cell Physiology, Max Delbrück Center for Molecular Medicine, Berlin, Germany
| | - Gregor Milicic
- Department of Mathematics, University of Salzburg, Salzburg, Austria
| | - Andreas Schröder
- Department of Mathematics, University of Salzburg, Salzburg, Austria
| | - Edda Klipp
- Theoretische Biophysik, Humboldt-Universität zu Berlin, Berlin, Germany
- * E-mail:
| |
Collapse
|
12
|
Hass H, Kipkeew F, Gauhar A, Bouché E, May P, Timmer J, Bock HH. Mathematical model of early Reelin-induced Src family kinase-mediated signaling. PLoS One 2017; 12:e0186927. [PMID: 29049379 PMCID: PMC5648249 DOI: 10.1371/journal.pone.0186927] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/07/2017] [Accepted: 10/10/2017] [Indexed: 12/23/2022] Open
Abstract
Reelin is a large glycoprotein with a dual role in the mammalian brain. It regulates the positioning and differentiation of postmitotic neurons during brain development and modulates neurotransmission and memory formation in the adult brain. Alterations in the Reelin signaling pathway have been described in different psychiatric disorders. Reelin mainly signals by binding to the lipoprotein receptors Vldlr and ApoER2, which induces tyrosine phosphorylation of the adaptor protein Dab1 mediated by Src family kinases (SFKs). In turn, phosphorylated Dab1 activates downstream signaling cascades, including PI3-kinase-dependent signaling. In this work, a mechanistic model based on ordinary differential equations was built to model early dynamics of the Reelin-mediated signaling cascade. Mechanistic models are frequently used to disentangle the highly complex mechanisms underlying cellular processes and obtain new biological insights. The model was calibrated on time-resolved data and a dose-response measurement of protein concentrations measured in cortical neurons treated with Reelin. It focusses on the interplay between Dab1 and SFKs with a special emphasis on the tyrosine phosphorylation of Dab1, and their role for the regulation of Reelin-induced signaling. Model selection was performed on different model structures and a comprehensive mechanistic model of the early Reelin signaling cascade is provided in this work. It emphasizes the importance of Reelin-induced lipoprotein receptor clustering for SFK-mediated Dab1 trans-phosphorylation and does not require co-receptors to describe the measured data. The model is freely available within the open-source framework Data2Dynamics (www.data2dynamics.org). It can be used to generate predictions that can be validated experimentally, and provides a platform for model extensions both to downstream targets such as transcription factors and interactions with other transmembrane proteins and neuronal signaling pathways.
Collapse
Affiliation(s)
- Helge Hass
- Institute of Physics, University of Freiburg, Freiburg, Germany
- * E-mail: (HH); (JT); (HHB)
| | - Friederike Kipkeew
- Clinic of Gastroenterology, Hepatology and Infectiology, Heinrich-Heine University Düsseldorf, Düsseldorf, Germany
| | - Aziz Gauhar
- Clinic of Gastroenterology, Hepatology and Infectiology, Heinrich-Heine University Düsseldorf, Düsseldorf, Germany
| | - Elisabeth Bouché
- Centre for Neuroscience, University of Freiburg, Freiburg, Germany
| | - Petra May
- Clinic of Gastroenterology, Hepatology and Infectiology, Heinrich-Heine University Düsseldorf, Düsseldorf, Germany
| | - Jens Timmer
- Institute of Physics, University of Freiburg, Freiburg, Germany
- BIOSS Centre for Biological Signaling Studies, University of Freiburg, Freiburg, Germany
- * E-mail: (HH); (JT); (HHB)
| | - Hans H. Bock
- Clinic of Gastroenterology, Hepatology and Infectiology, Heinrich-Heine University Düsseldorf, Düsseldorf, Germany
- * E-mail: (HH); (JT); (HHB)
| |
Collapse
|
13
|
Bravi B, Sollich P. Statistical physics approaches to subnetwork dynamics in biochemical systems. Phys Biol 2017; 14:045010. [PMID: 28510539 DOI: 10.1088/1478-3975/aa7363] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Abstract
We apply a Gaussian variational approximation to model reduction in large biochemical networks of unary and binary reactions. We focus on a small subset of variables (subnetwork) of interest, e.g. because they are accessible experimentally, embedded in a larger network (bulk). The key goal is to write dynamical equations reduced to the subnetwork but still retaining the effects of the bulk. As a result, the subnetwork-reduced dynamics contains a memory term and an extrinsic noise term with non-trivial temporal correlations. We first derive expressions for this memory and noise in the linearized (Gaussian) dynamics and then use a perturbative power expansion to obtain first order nonlinear corrections. For the case of vanishing intrinsic noise, our description is explicitly shown to be equivalent to projection methods up to quadratic terms, but it is applicable also in the presence of stochastic fluctuations in the original dynamics. An example from the epidermal growth factor receptor signalling pathway is provided to probe the increased prediction accuracy and computational efficiency of our method.
Collapse
Affiliation(s)
- B Bravi
- Current affiliation: Institute of Theoretical Physics, Ecole Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| | | |
Collapse
|