MODELLING β2AR REGULATION

The β2 adrenergic receptor (β2AR) regulates smooth muscle relaxation in the vasculature and airways. Long- and Short-acting β-agonists (LABAs/SABAs) are widely used in treatment of chronic obstructive pulmonary disorder (COPD) and asthma. Despite their widespread clinical use we do not understand well the dominant β2AR regulatory pathways that are stimulated during therapy and bring about tachyphylaxis, which is the loss of drug effects. Thus, an understanding of how the β2AR responds to various β-agonists is crucial to their rational use. Towards that end we have developed deterministic models that explore the mechanism of drug- induced β2AR regulation. These mathematical models can be classified into three classes; (i) Six quantitative models of SABA-induced G protein coupled receptor kinase (GRK)-mediated β2AR regulation; (ii) Three phenomenological models of salmeterol (a LABA)-induced GRK-mediated β2AR regulation; and (iii) One semi-quantitative, unified model of SABA-induced GRK-, protein kinase A (PKA)-, and phosphodiesterase (PDE)-mediated regulation of β2AR signalling. The various models were constrained with all or some of the following experimental data; (i) GRK-mediated β2AR phosphorylation in response to various LABAs/SABAs; (ii) dephosphorylation of the GRK site on the β2AR; (iii) β2AR internalisation; (iv) β2AR recycling; (v) β2AR desensitisation; (vi) β2AR resensitisation; (vii) PKA-mediated β2AR phosphorylation in response to a SABA; and (viii) LABA/SABA induced cAMP profile ± PDE inhibitors. The models of GRK-mediated β2AR regulation show that plasma membrane dephosphorylation and recycling of the phosphorylated β2AR are required to reconcile with the measured dephosphorylation kinetics. We further used a consensus model to predict the consequences of rapid pulsatile agonist stimulation and found that although resensitisation was rapid, the β2AR system retained the memory of prior stimuli and desensitised much more rapidly and strongly in response to subsequent stimuli. This could explain tachyphylaxis of SABAs over repeated use in rescue therapy of asthma patients. The LABA models show that the long action of salmeterol can be explained due to decreased stability of the arrestin/β2AR/salmeterol complex. This could explain long action of β-agonists used in maintenance therapy of asthma patients. Our consensus model of PKA/PDE/GRK-mediated β2AR regulation is being used to identify the dominant β2AR desensitisation pathways under different therapeutic regimens in human airway cells. In summary our models represent a significant advance towards understanding agonist-specific β2AR regulation that will aid in a more rational use of the β2AR agonists in the treatment of asthma

Electronic data sources for kinetic models of cell signaling

Functional understanding of signaling pathways requires detailed information about the constituent molecules and their interactions. Simulations of signaling pathways therefore build upon a great deal of data from various sources. We first survey electronic data resources for cell signaling modeling and then based on the type of data representation the data sources are broadly classified into five groups. None of the data sources surveyed provide all required data in a ready-to-be-modeled fashion. We then put forward a wish list for the desired attributes for an ideal modeling centric database. Finally, we close with perspectives on how electronic data sources for cell signaling modeling have developed. We suggest that future directions in such data sources are largely model-driven and are hinged on interoperability of data sources

Fluvoxamine: A review of its mechanism of action and its role in COVID-19

Fluvoxamine is a well-tolerated, widely available, inexpensive selective serotonin reuptake inhibitor that has been shown in a small, double-blind, placebo-controlled, randomized study to prevent clinical deterioration of patients with mild coronavirus disease 2019 (COVID-19). Fluvoxamine is also an agonist for the sigma-1 receptor, through which it controls inflammation. We review here a body of literature that shows important mechanisms of action of fluvoxamine and other SSRIs that could play a role in COVID-19 treatment. These effects include: reduction in platelet aggregation, decreased mast cell degranulation, interference with endolysosomal viral trafficking, regulation of inositol-requiring enzyme 1α-driven inflammation and increased melatonin levels, which collectively have a direct antiviral effect, regulate coagulopathy or mitigate cytokine storm, which are known hallmarks of severe COVID-19

When kinases meet mathematics: the systems biology of MAPK signalling

The mitogen activated protein kinase/extracellular signal regulated kinase pathway regulates fundamental cellular function such as cell proliferation, survival, differentiation and motility, raising the question how these diverse functions are specified and coordinated. They are encoded through the activation kinetics of the pathway, a multitude of feedback loops, scaffold proteins, subcellular compartmentalisation, and crosstalk with other pathways. These regulatory motifs alone or in combination can generate a multitude of complex behaviour. Systems biology tries to decode this complexity through mathematical modelling and prediction in order to gain a deeper insight into the inner works of signalling networks

Quantitative Modeling of GRK-Mediated β2AR Regulation

We developed a unified model of the GRK-mediated β2 adrenergic receptor (β2AR) regulation that simultaneously accounts for six different biochemical measurements of the system obtained over a wide range of agonist concentrations. Using a single deterministic model we accounted for (1) GRK phosphorylation in response to various full and partial agonists; (2) dephosphorylation of the GRK site on the β2AR; (3) β2AR internalization; (4) recycling of the β2AR post isoproterenol treatment; (5) β2AR desensitization; and (6) β2AR resensitization. Simulations of our model show that plasma membrane dephosphorylation and recycling of the phosphorylated receptor are necessary to adequately account for the measured dephosphorylation kinetics. We further used the model to predict the consequences of (1) modifying rates such as GRK phosphorylation of the receptor, arrestin binding and dissociation from the receptor, and receptor dephosphorylation that should reflect effects of knockdowns and overexpressions of these components; and (2) varying concentration and frequency of agonist stimulation “seen” by the β2AR to better mimic hormonal, neurophysiological and pharmacological stimulations of the β2AR. Exploring the consequences of rapid pulsatile agonist stimulation, we found that although resensitization was rapid, the β2AR system retained the memory of the previous stimuli and desensitized faster and much more strongly in response to subsequent stimuli. The latent memory that we predict is due to slower membrane dephosphorylation, which allows for progressive accumulation of phosphorylated receptor on the surface. This primes the receptor for faster arrestin binding on subsequent agonist activation leading to a greater extent of desensitization. In summary, the model is unique in accounting for the behavior of the β2AR system across multiple types of biochemical measurements using a single set of experimentally constrained parameters. It also provides insight into how the signaling machinery can retain memory of prior stimulation long after near complete resensitization has been achieved

Wnt Signalling Pathway Parameters for Mammalian Cells

Wnt/β-catenin signalling regulates cell fate, survival, proliferation and differentiation at many stages of mammalian development and pathology. Mutations of two key proteins in the pathway, APC and β-catenin, have been implicated in a range of cancers, including colorectal cancer. Activation of Wnt signalling has been associated with the stabilization and nuclear accumulation of β-catenin and consequential up-regulation of β-catenin/TCF gene transcription. In 2003, Lee et al. constructed a computational model of Wnt signalling supported by experimental data from analysis of time-dependent concentration of Wnt signalling proteins in Xenopus egg extracts. Subsequent studies have used the Xenopus quantitative data to infer Wnt pathway dynamics in other systems. As a basis for understanding Wnt signalling in mammalian cells, a confocal live cell imaging measurement technique is developed to measure the cell and nuclear volumes of MDCK, HEK293T cells and 3 human colorectal cancer cell lines and the concentrations of Wnt signalling proteins β-catenin, Axin, APC, GSK3β and E-cadherin. These parameters provide the basis for formulating Wnt signalling models for kidney/intestinal epithelial mammalian cells. There are significant differences in concentrations of key proteins between Xenopus extracts and mammalian whole cell lysates. Higher concentrations of Axin and lower concentrations of APC are present in mammalian cells. Axin concentrations are greater than APC in kidney epithelial cells, whereas in intestinal epithelial cells the APC concentration is higher than Axin. Computational simulations based on Lee's model, with this new data, suggest a need for a recalibration of the model

Why Are Computational Neuroscience and Systems Biology So Separate?

Despite similar computational approaches, there is surprisingly little interaction between the computational neuroscience and the systems biology research communities. In this review I reconstruct the history of the two disciplines and show that this may explain why they grew up apart. The separation is a pity, as both fields can learn quite a bit from each other. Several examples are given, covering sociological, software technical, and methodological aspects. Systems biology is a better organized community which is very effective at sharing resources, while computational neuroscience has more experience in multiscale modeling and the analysis of information processing by biological systems. Finally, I speculate about how the relationship between the two fields may evolve in the near future

Developing complex signaling models using GENESIS/Kinetikit

The development of biologically realistic models of signaling pathways is a demanding process, involving computational challenges as well as those arising from the complexity of detailed pathway models. We have developed the General Neural Simulation System (GENESIS) and Kinetikit (GENESIS/Kinetikit), a graphical simulation environment for modeling biochemical signaling pathways using deterministic and stochastic methods. A library of models of several common signaling pathways complements the software. This combination of numerical computation engines, graphical modeling tools, and library of models is designed to build on the cumulative development of models and techniques from many sources. The complete simulation environment and demonstration models are available from (http://stke.sciencemag.org/cgi/content/full/sigtrans;2004/219/pl4/DC1; also at http://www.ncbs.res.in/~bhalla/kkit/download.html). The associated library of signaling pathways is based on published experimental and simulation studies and is curated to ensure that the simulation outcomes match published results. Models in the library are maintained in a database (http://doqcs.ncbs.res.in). Individual pathway models can be combined to build complex signaling network simulations. The overall goal of this process is to attain sufficient biological realism in models to directly compare their outcomes with experiments and to improve our understanding of complex signaling