A kinetic model for calcium dynamics in RAW 264.7 cells: 2. Knockdown response and long-term response

Modeling transcriptional regulation of the cell cycle using a novel cybernetic-inspired approach

Quantitative understanding of cellular processes, such as cell cycle and differentiation, is impeded by various forms of complexity ranging from myriad molecular players and their multilevel regulatory interactions, cellular evolution with multiple intermediate stages, lack of elucidation of cause-effect relationships among the many system players, and the computational complexity associated with the profusion of variables and parameters. In this paper, we present a modeling framework based on the cybernetic concept that biological regulation is inspired by objectives embedding rational strategies for dimension reduction, process stage specification through the system dynamics, and innovative causal association of regulatory events with the ability to predict the evolution of the dynamical system. The elementary step of the modeling strategy involves stage-specific objective functions that are computationally determined from experiments, augmented with dynamical network computations involving endpoint objective functions, mutual information, change-point detection, and maximal clique centrality. We demonstrate the power of the method through application to the mammalian cell cycle, which involves thousands of biomolecules engaged in signaling, transcription, and regulation. Starting with a fine-grained transcriptional description obtained from RNA sequencing measurements, we develop an initial model, which is then dynamically modeled using the cybernetic-inspired method, based on the strategies described above. The cybernetic-inspired method is able to distill the most significant interactions from a multitude of possibilities. In addition to capturing the complexity of regulatory processes in a mechanistically causal and stage-specific manner, we identify the functional network modules, including novel cell cycle stages. Our model is able to predict future cell cycles consistent with experimental measurements. We posit that this innovative framework has the promise to extend to the dynamics of other biological processes, with a potential to provide novel mechanistic insights.

Modeling transcriptional regulation of the cell cycle using a novel cybernetic-inspired approach

Quantitative understanding of cellular processes, such as cell cycle and differentiation, is impeded by various forms of complexity ranging from myriad molecular players and their multilevel regulatory interactions, cellular evolution with multiple intermediate stages, lack of elucidation of cause-effect relationships among the many system players, and the computational complexity associated with the profusion of variables and parameters. In this paper, we present an elegant modeling framework based on the cybernetic concept that biological regulation is inspired by objectives embedding entirely novel strategies for dimension reduction, process stage specification through the system dynamics, and innovative causal association of regulatory events with the ability to predict the evolution of the dynamical system. The elementary step of the modeling strategy involves stage-specific objective functions that are computationally-determined from experiments, augmented with dynamical network computations involving end point objective functions, mutual information, change point detection, and maximal clique centrality. We demonstrate the power of the method through application to the mammalian cell cycle, which involves thousands of biomolecules engaged in signaling, transcription, and regulation. Starting with a fine-grained transcriptional description obtained from RNA sequencing measurements, we develop an initial model, which is then dynamically modeled using the cybernetic-inspired method (CIM), utilizing the strategies described above. The CIM is able to distill the most significant interactions from a multitude of possibilities. In addition to capturing the complexity of regulatory processes in a mechanistically causal and stage-specific manner, we identify the functional network modules, including novel cell cycle stages. Our model is able to predict future cell cycles consistent with experimental measurements. We posit that this state-of-the-art framework has the promise to extend to the dynamics of other biological processes, with a potential to provide novel mechanistic insights.

STATEMENT OF SIGNIFICANCE

Cellular processes like cell cycle are overly complex, involving multiple players interacting at multiple levels, and explicit modeling of such systems is challenging. The availability of longitudinal RNA measurements provides an opportunity to "reverse-engineer" for novel regulatory models. We develop a novel framework, inspired using goal-oriented cybernetic model, to implicitly model transcriptional regulation by constraining the system using inferred temporal goals. A preliminary causal network based on information-theory is used as a starting point, and our framework is used to distill the network to temporally-based networks containing essential molecular players. The strength of this approach is its ability to dynamically model the RNA temporal measurements. The approach developed paves the way for inferring regulatory processes in many complex cellular processes.

Modeling of Endothelial Calcium Responses within a Microfluidic Generator of Spatio-Temporal ATP and Shear Stress Signals

Intracellular calcium dynamics play essential roles in the proper functioning of cellular activities. It is a well known important chemosensing and mechanosensing process regulated by the spatio-temporal microenvironment. Nevertheless, how spatio-temporal biochemical and biomechanical stimuli affect calcium dynamics is not fully understood and the underlying regulation mechanism remains missing. Herein, based on a developed microfluidic generator of biochemical and biomechanical signals, we theoretically analyzed the generation of spatio-temporal ATP and shear stress signals within the microfluidic platform and investigated the effect of spatial combination of ATP and shear stress stimuli on the intracellular calcium dynamics. The simulation results demonstrate the capacity and flexibility of the microfluidic system in generating spatio-temporal ATP and shear stress. Along the transverse direction of the microchannel, dynamic ATP signals of distinct amplitudes coupled with identical shear stress are created, which induce the spatio-temporal diversity in calcium responses. Interestingly, to the multiple combinations of stimuli, the intracellular calcium dynamics reveal two main modes: unimodal and oscillatory modes, showing significant dependence on the features of the spatio-temporal ATP and shear stress stimuli. The present study provides essential information for controlling calcium dynamics by regulating spatio-temporal biochemical and biomechanical stimuli, which shows the potential in directing cellular activities and understanding the occurrence and development of disease.

GPCR mediated control of calcium dynamics: A systems perspective

G-protein coupled receptor (GPCR) mediated calcium (Ca²⁺)-signaling transduction remains crucial in designing drugs for various complex diseases including neurodegeneration, chronic heart failure as well as respiratory diseases. Although there are several reviews detailing various aspects of Ca²⁺-signaling such as the role of IP₃ receptors and Ca²⁺-induced-Ca²⁺-release, none of them provide an integrated view of the mathematical descriptions of GPCR signal transduction and investigations on dose-response curves. This article is the first study in reviewing the network structures underlying GPCR signal transduction that control downstream [Ca_c²⁺]-oscillations. The central theme of this paper is to present the biochemical pathways, as well as molecular mechanisms underlying the GPCR-mediated Ca²⁺-dynamics in order to facilitate a better understanding of how agonist concentration is encoded in Ca²⁺-signals for G_αq, G_αs, and G_αi/o signaling pathways. Moreover, we present the GPCR targeting drugs that are relevant for treating cardiac, respiratory, and neuro-diseases. The current paper presents the ODE formulation for various models along with the detailed schematics of signaling networks. To provide a systems perspective, we present the network motifs that can provide readers an insight into the complex and intriguing science of agonist-mediated Ca²⁺-dynamics. One of the features of this review is to pinpoint the interplay between positive and negative feedback loops that are involved in controlling intracellular [Ca_c²⁺]-oscillations. Furthermore, we review several examples of dose-response curves obtained from [Ca_c²⁺]-spiking for various GPCR pathways. This paper is expected to be useful for pharmacologists and computational biologists for designing clinical applications of GPCR targeting drugs through modulation of Ca²⁺-dynamics.

Inhibition of osteoclast activity by complement regulation with DF3016A, a novel small-molecular-weight C5aR inhibitor

Recent insights have indicated an active role of the complex complement system not only in immunity, but also in bone remodeling. Evidence from knockout mice and observations from skeletal diseases have drawn attention to the C5a/C5aR axis of the complement cascade in the modulation of osteoclast functions and as potential therapeutic targets for treatment of bone pathologies. With the aim to identify novel C5aR regulators, a medicinal chemistry program was initiated, driven by structural information on a minor pocket of C5aR that has been proposed to be a key motif for C5aR intracellular activation. The impact of the peptidomimetic orthosteric C5aR antagonist (PMX-53), of two newly synthesized allosteric C5aR antagonists (DF2593A, DF3016A), and of C5aR down-regulation by specific siRNAs, were examined for regulation of osteoclastogenesis, using a well-validated in-vitro model starting from RAW264.7 precursor cells. Both pharmacological and molecular approaches reduced osteoclast maturation of RAW264.7 cells induced by receptor-activator of nuclear factor kappa-B ligand (RANKL), which limited the transcription of several differentiation markers evaluated by real-time PCR, including nuclear factor of activated T-cell 1, matrix metalloproteinase-9, cathepsin-K, and tartrate-resistant acid phosphatase. These treatments were ineffective on the subsequent step of osteoclast syncytium formation, apparently as a consequence of reduction of C5aR mRNA levels in the course of osteoclastogenesis, as monitored by real-time PCR. Among the C5aR antagonists analyzed, DF3016A inhibited osteoclast degradation activity through inhibition of C5aR signal transduction and transcription. These data confirm the preclinical relevance of this novel therapeutic candidate.

Calcium Dynamics of Ex Vivo Long-Term Cultured CD8+ T Cells Are Regulated by Changes in Redox Metabolism

T cells reach a state of replicative senescence characterized by a decreased ability to proliferate and respond to foreign antigens. Calcium release associated with TCR engagement is widely used as a surrogate measure of T cell response. Using an ex vivo culture model that partially replicates features of organismal aging, we observe that while the amplitude of Ca2+ signaling does not change with time in culture, older T cells exhibit faster Ca2+ rise and a faster decay. Gene expression analysis of Ca2+ channels and pumps expressed in T cells by RT-qPCR identified overexpression of the plasma membrane CRAC channel subunit ORAI1 and PMCA in older T cells. To test whether overexpression of the plasma membrane Ca2+ channel is sufficient to explain the kinetic information, we adapted a previously published computational model by Maurya and Subramaniam to include additional details on the store-operated calcium entry (SOCE) process to recapitulate Ca2+ dynamics after T cell receptor stimulation. Simulations demonstrated that upregulation of ORAI1 and PMCA channels is not sufficient to explain the observed alterations in Ca2+ signaling. Instead, modeling analysis identified kinetic parameters associated with the IP3R and STIM1 channels as potential causes for alterations in Ca2+ dynamics associated with the long term ex vivo culturing protocol. Due to these proteins having known cysteine residues susceptible to oxidation, we subsequently investigated and observed transcriptional remodeling of metabolic enzymes, a shift to more oxidized redox couples, and post-translational thiol oxidation of STIM1. The model-directed findings from this study highlight changes in the cellular redox environment that may ultimately lead to altered T cell calcium dynamics during immunosenescence or organismal aging.

The macrophage and its role in inflammation and tissue repair: mathematical and systems biology approaches

Macrophages are central to the inflammatory response and its ability to resolve effectively. They are complex cells that adopt a range of subtypes depending on the tissue type and stimulus that they find themselves under. This flexibility allows them to play multiple, sometimes opposing, roles in inflammation and tissue repair. Their central role in the inflammatory process is reflected in macrophage dysfunction being implicated in chronic inflammation and poorly healing wounds. In this study, we discuss recent attempts to model mathematically and computationally the macrophage and how it partakes in the complex processes of inflammation and tissue repair. There are increasing data describing the variety of macrophage phenotypes and their underlying transcriptional programs. Dynamic mathematical and computational models are an ideal way to test biological hypotheses against experimental data and could aid in understanding this multi-functional cell and its potential role as an attractive therapeutic target for inflammatory conditions and tissue repair.

Microfluidic interrogation and mathematical modeling of multi-regime calcium signaling dynamics

Through microfluidic interrogation we analyzed real-time calcium responses of HEK293 cells stimulated with short pulses of the M3 muscarinic receptor ligand carbachol in two different concentration regimes. Lower ligand concentrations elicit oscillatory calcium signals while higher concentrations trigger a rapid rise that eventually settles down at a steady-state slightly above pre-stimulus levels, referred to as an acute signal. Cells were periodically pulsed with carbachol at these two concentration regimes using a custom-made microfluidic platform, and the resulting calcium signals were measured with a single fluorescent readout. Pulsed stimulations at these two concentration regimes resulted in multiple types of response patterns that each delivered complementary information about the M3 muscarinic receptor signaling pathway. These multiple types of calcium response patterns enabled development of a comprehensive mathematical model of multi-regime calcium signaling. The resulting model suggests that dephosphorylation of deactivated receptors is rate limiting for recovery of calcium signals in the acute regime (high ligand concentration), while calcium replenishment and IP3 production determine signal recovery in the oscillatory regime (low ligand concentration). This study not only provides mechanistic insight into multi-regime signaling of the M3 muscarinic receptor pathway, but also provides a general strategy for analyzing multi-regime pathways using only one fluorescent readout.

Functional classification of skeletal muscle networks. I. Normal physiology

Extensive measurements of the parts list of human skeletal muscle through transcriptomics and other phenotypic assays offer the opportunity to reconstruct detailed functional models. Through integration of vast amounts of data present in databases and extant knowledge of muscle function combined with robust analyses that include a clustering approach, we present both a protein parts list and network models for skeletal muscle function. The model comprises the four key functional family networks that coexist within a functional space; namely, excitation-activation family (forward pathways that transmit a motoneuronal command signal into the spatial volume of the cell and then use Ca(2+) fluxes to bind Ca(2+) to troponin C sites on F-actin filaments, plus transmembrane pumps that maintain transmission capacity); mechanical transmission family (a sophisticated three-dimensional mechanical apparatus that bidirectionally couples the millions of actin-myosin nanomotors with external axial tensile forces at insertion sites); metabolic and bioenergetics family (pathways that supply energy for the skeletal muscle function under widely varying demands and provide for other cellular processes); and signaling-production family (which represents various sensing, signal transduction, and nuclear infrastructure that controls the turn over and structural integrity and regulates the maintenance, regeneration, and remodeling of the muscle). Within each family, we identify subfamilies that function as a unit through analysis of large-scale transcription profiles of muscle and other tissues. This comprehensive network model provides a framework for exploring functional mechanisms of the skeletal muscle in normal and pathophysiology, as well as for quantitative modeling.

Exploring phospholipase C-coupled Ca(2+) signalling networks using Boolean modelling

In this study, the authors explored the utility of a descriptive and predictive bionetwork model for phospholipase C-coupled calcium signalling pathways, built with non-kinetic experimental information. Boolean models generated from these data yield oscillatory activity patterns for both the endoplasmic reticulum resident inositol-1,4,5-trisphosphate receptor (IP(3)R) and the plasma-membrane resident canonical transient receptor potential channel 3 (TRPC3). These results are specific as randomisation of the Boolean operators ablates oscillatory pattern formation. Furthermore, knock-out simulations of the IP(3)R, TRPC3 and multiple other proteins recapitulate experimentally derived results. The potential of this approach can be observed by its ability to predict previously undescribed cellular phenotypes using in vitro experimental data. Indeed, our cellular analysis of the developmental and calcium-regulatory protein, DANGER1a, confirms the counter-intuitive predictions from our Boolean models in two highly relevant cellular models. Based on these results, the authors theorise that with sufficient legacy knowledge and/or computational biology predictions, Boolean networks can provide a robust method for predictive modelling of any biological system. [Includes supplementary material].

Stochastic hybrid modeling of intracellular calcium dynamics

Deterministic models of biochemical processes at the subcellular level might become inadequate when a cascade of chemical reactions is induced by a few molecules. Inherent randomness of such phenomena calls for the use of stochastic simulations. However, being computationally intensive, such simulations become infeasible for large and complex reaction networks. To improve their computational efficiency in handling these networks, we present a hybrid approach, in which slow reactions and fluxes are handled through exact stochastic simulation and their fast counterparts are treated partially deterministically through chemical Langevin equation. The classification of reactions as fast or slow is accompanied by the assumption that in the time-scale of fast reactions, slow reactions do not occur and hence do not affect the probability of the state. Our new approach also handles reactions with complex rate expressions such as Michaelis-Menten kinetics. Fluxes which cannot be modeled explicitly through reactions, such as flux of Ca(2+) from endoplasmic reticulum to the cytosol through inositol 1,4,5-trisphosphate receptor channels, are handled deterministically. The proposed hybrid algorithm is used to model the regulation of the dynamics of cytosolic calcium ions in mouse macrophage RAW 264.7 cells. At relatively large number of molecules, the response characteristics obtained with the stochastic and deterministic simulations coincide, which validates our approach in the limit of large numbers. At low doses, the response characteristics of some key chemical species, such as levels of cytosolic calcium, predicted with stochastic simulations, differ quantitatively from their deterministic counterparts. These observations are ubiquitous throughout dose response, sensitivity, and gene-knockdown response analyses. While the relative differences between the peak-heights of the cytosolic [Ca(2+)] time-courses obtained from stochastic (mean of 16 realizations) and deterministic simulations are merely 1%-4% for most perturbations, it is specially sensitive to levels of G(βγ) (relative difference as large as 90% at very low G(βγ)).

PKC and PLA2: probing the complexities of the calcium network

Lipid signaling and phosphorylation cascades are fundamental to calcium signaling networks. In this review, we will discuss the recent laboratory findings for the phospholipase A(2) (PLA(2))/protein kinase C (PKC) pathway within cellular calcium networks. The complexity and connectivity of these ubiquitous cellular signals make interpretation of experimental results extremely challenging. We present here computational methods which have been developed to conquer such complex data, and how they can be used to make models capable of accurately predicting cellular responses within multiple calcium signaling pathways. We propose that information obtained from network analysis and computational techniques provides a rich source of knowledge which can be directly translated to the laboratory benchtop.

A dual receptor crosstalk model of G-protein-coupled signal transduction

Macrophage cells that are stimulated by two different ligands that bind to G-protein-coupled receptors (GPCRs) usually respond as if the stimulus effects are additive, but for a minority of ligand combinations the response is synergistic. The G-protein-coupled receptor system integrates signaling cues from the environment to actuate cell morphology, gene expression, ion homeostasis, and other physiological states. We analyze the effects of the two signaling molecules complement factors 5a (C5a) and uridine diphosphate (UDP) on the intracellular second messenger calcium to elucidate the principles that govern the processing of multiple signals by GPCRs. We have developed a formal hypothesis, in the form of a kinetic model, for the mechanism of action of this GPCR signal transduction system using data obtained from RAW264.7 macrophage cells. Bayesian statistical methods are employed to represent uncertainty in both data and model parameters and formally tie the model to experimental data. When the model is also used as a tool in the design of experiments, it predicts a synergistic region in the calcium peak height dose response that results when cells are simultaneously stimulated by C5a and UDP. An analysis of the model reveals a potential mechanism for crosstalk between the Galphai-coupled C5a receptor and the Galphaq-coupled UDP receptor signaling systems that results in synergistic calcium release.

Systems biology of macrophages

Cells and tissues function in context. Under a given growth or survival medium they perform tasks, replicate and die. Given a stimulus they respond by invoking myriad biomolecular networks that result in a specified cellular outcome. At any given instant it can be argued that the cell is in a "state" defined by its components--their concentrations and locations, the interactions between components--that are modulated in space and time, and the complex circuitry--that involves a large number of interacting networks and a snapshot of the dynamical processes--such as gene expression, cell cycle, transport of components, etc. At present, we can measure, using high and low throughput methods, several cellular components in a context-dependent manner and obtain a partial picture of cellular networks and dynamical processes. Are these measurements sufficient to answer important biological questions and help reconstruct a systems-level understanding of a mammalian cell? This chapter will address systems biology strategies developed to address this question and demonstrate the power of integration of diverse cellular data for answering interesting biological questions in macrophages. We will use this systems biology approach to address the following questions: (1) How good are macrophage cell lines in addressing phenotypic biology of primary macrophages? (2) How do signals associated with inflammatory molecules regulate gene transcription in macrophages? (3) How can we combine proteomic and other cellular measurements to characterize the repertoire of upstream signaling networks invoked by macrophages? (4) How do designed knockdowns of proteins influence cellular phenotypes?