SBMLBioModels: W - Ç

The electrical activity of insulin-secreting pancreatic islets of Langerhans is characterized by bursts of action potentials. Most often this bursting is periodic, but in some cases it is modulated by an underlying slower rhythm. We suggest that the modulatory rhythm for this complex bursting pattern is due to oscillations in glycolysis, while the bursting itself is generated by some other slow process. To demonstrate this hypothesis, we couple a minimal model of glycolytic oscillations to a minimal model for activity-dependent bursting in islets. We show that the combined model can reproduce several complex bursting patterns from mouse islets published in the literature, and we illustrate how these complex oscillations are produced through the use of a fast/slow analysis. link: http://identifiers.org/pubmed/15178199

Bistability underlies basic biological phenomena, such as cell division, differentiation, cancer onset, and apoptosis. So far biologists identified two necessary conditions for bistability: positive feedback and ultrasensitivity.Biological systems are based upon elementary mono- and bimolecular chemical reactions. In order to definitely clarify all necessary conditions for bistability we here present the corresponding minimal system. According to our definition, it contains the minimal number of (i) reactants, (ii) reactions, and (iii) terms in the corresponding ordinary differential equations (decreasing importance from i-iii). The minimal bistable system contains two reactants and four irreversible reactions (three bimolecular, one monomolecular).We discuss the roles of the reactions with respect to the necessary conditions for bistability: two reactions comprise the positive feedback loop, a third reaction filters out small stimuli thus enabling a stable 'off' state, and the fourth reaction prevents explosions. We argue that prevention of explosion is a third general necessary condition for bistability, which is so far lacking discussion in the literature.Moreover, in addition to proving that in two-component systems three steady states are necessary for bistability (five for tristability, etc.), we also present a simple general method to design such systems: one just needs one production and three different degradation mechanisms (one production, five degradations for tristability, etc.). This helps modelling multistable systems and it is important for corresponding synthetic biology projects.The presented minimal bistable system finally clarifies the often discussed question for the necessary conditions for bistability. The three necessary conditions are: positive feedback, a mechanism to filter out small stimuli and a mechanism to prevent explosions. This is important for modelling bistability with simple systems and for synthetically designing new bistable systems. Our simple model system is also well suited for corresponding teaching purposes. link: http://identifiers.org/pubmed/19737387

Cancer dormancy, a state in which cancer cells persist in a host without significant growth, is a natural forestallment of progression to manifest disease and is thus of great clinical interest. Experimental work in mice suggests that in immune-induced dormancy, the longer a cancer remains dormant in a host, the more resistant the cancer cells become to cytotoxic T-cell-mediated killing. In this work, mathematical models are used to analyse the possible causative mechanisms of cancer escape from immune-induced dormancy. Using a data-driven approach, both decaying efficacy in immune predation and immune recruitment are analysed with results suggesting that decline in recruitment is a stronger determinant of escape than increased resistance to predation. Using a mechanistic approach, the existence of an immune-resistant cancer cell subpopulation is considered, and the effects on cancer dormancy and potential immunoediting mechanisms of cancer escape are analysed and discussed. The immunoediting mechanism assumes that the immune system selectively prunes the cancer of immune-sensitive cells, which is shown to cause an initially heterogeneous population to become a more homogeneous, and more resistant, population. The fact that this selection may result in the appearance of decreasing efficacy in T-cell cytotoxic effect with time in dormancy is also demonstrated. This work suggests that through actions that temporarily delay cancer growth through the targeted removal of immune-sensitive subpopulations, the immune response may actually progress the cancer to a more aggressive state. link: http://identifiers.org/pubmed/24511375

Cancer dormancy, a state in which cancer cells persist in a host without significant growth, is a natural forestallment of progression to manifest disease and is thus of great clinical interest. Experimental work in mice suggests that in immune-induced dormancy, the longer a cancer remains dormant in a host, the more resistant the cancer cells become to cytotoxic T-cell-mediated killing. In this work, mathematical models are used to analyse the possible causative mechanisms of cancer escape from immune-induced dormancy. Using a data-driven approach, both decaying efficacy in immune predation and immune recruitment are analysed with results suggesting that decline in recruitment is a stronger determinant of escape than increased resistance to predation. Using a mechanistic approach, the existence of an immune-resistant cancer cell subpopulation is considered, and the effects on cancer dormancy and potential immunoediting mechanisms of cancer escape are analysed and discussed. The immunoediting mechanism assumes that the immune system selectively prunes the cancer of immune-sensitive cells, which is shown to cause an initially heterogeneous population to become a more homogeneous, and more resistant, population. The fact that this selection may result in the appearance of decreasing efficacy in T-cell cytotoxic effect with time in dormancy is also demonstrated. This work suggests that through actions that temporarily delay cancer growth through the targeted removal of immune-sensitive subpopulations, the immune response may actually progress the cancer to a more aggressive state. link: http://identifiers.org/pubmed/24511375

A simulation model for the production of thrombin in plasma is presented. Values of the reaction rate constants as determined in purified systems are used and the model is tested by comparison of simulations of factor Xa, factor Va and thrombin generation curves with experimental data obtained in thromboplastin-activated plasma. Simulations of the effect of hirudin indicate that factor V is predominantly activated by thrombin and not by factor Xa. The model predicts a threshold value for the factor Xa production which, if exceeded, results in explosive and complete activation of prothrombinase. The dependence of this threshold value on different negative feedback reactions, e.g. the inactivation of thrombin and factor Xa by antithrombin III (+ heparin), is investigated. The threshold value, for control plasma in the range of 1-10 pM total factor Xa production, can be raised two orders of magnitude by accelerated inactivation of factor Xa and prothrombinase but is hardly affected by a tenfold increase in the rate of thrombin inactivation or by increased production of activated protein C. This latter effect, however, results in a more gradual input-response relation between factor Xa input and the extent of prothrombinase activation. link: http://identifiers.org/pubmed/1794746

TGF-β is an immunoregulatory protein that contributes to inadequate antitumor immune responses in cancer patients. Recent experimental data suggests that TGF-β inhibition alone, provides few clinical benefits, yet it can significantly amplify the anti-tumor immune response when combined with a tumor vaccine. We develop a mathematical model in order to gain insight into the cooperative interaction between anti-TGF-β and vaccine treatments. The mathematical model follows the dynamics of the tumor size, TGF-β concentration, activated cytotoxic effector cells, and regulatory T cells. Using numerical simulations and stability analysis, we study the following scenarios: a control case of no treatment, anti-TGF-β treatment, vaccine treatment, and combined anti-TGF-β vaccine treatments. We show that our model is capable of capturing the observed experimental results, and hence can be potentially used in designing future experiments involving this approach to immunotherapy. link: http://identifiers.org/pubmed/22438084

Blood oxygen level-dependent functional magnetic resonance imaging (BOLD-fMRI) is a standard clinical tool for the detection of brain activation. In Alzheimer’s disease (AD), task-related and resting state fMRI have been used to detect brain dysfunction. It has been shown that the shape of the BOLD response is affected in early AD. To correctly interpret these changes, the mechanisms responsible for the observed behaviour need to be known. The parameters of the canonical hemodynamic response function (HRF) commonly used in the analysis of fMRI data have no direct biological interpretation and cannot be used to answer this question. We here present a model that allows relating AD-specific changes in the BOLD shape to changes in the underlying energy metabolism. According to our findings, the classic view that differences in the BOLD shape are only attributed to changes in strength and duration of the stimulus does not hold. Instead, peak height, peak timing and full width at half maximum are sensitive to changes in the reaction rate of several metabolic reactions. Our systems-theoretic approach allows the use of patient-specific clinical data to predict dementia- driven changes in the HRF, which can be used to improve the results of fMRI analyses in AD patients. link: http://identifiers.org/doi/10.1177/0271678X17693024

We use mathematical models to study the relationship between HIV and the immune system during the natural course of infection and in the context of different antiviral treatment regimes. The models suggest that an efficient cytotoxic T lymphocyte (CTL) memory response is required to control the virus. We define CTL memory as long-term persistence of CTL precursors in the absence of antigen. Infection and depletion of CD4(+) T helper cells interfere with CTL memory generation, resulting in persistent viral replication and disease progression. We find that antiviral drug therapy during primary infection can enable the development of CTL memory. In chronically infected patients, specific treatment schedules, either including deliberate drug holidays or antigenic boosts of the immune system, can lead to a re-establishment of CTL memory. Whether such treatment regimes would lead to long-term immunologic control deserves investigation under carefully controlled conditions. link: http://identifiers.org/pubmed/10588728

Cytotoxic T lymphocytes (CTL) responses are required to fight many diseases such as viral infections and tumors. At the same time, they can cause disease when induced inappropriately. Which factors regulate CTL and decide whether they should remain silent or react is open to debate. The phenomenon called cross-priming has received attention in this respect. That is, CTL expansion occurs if antigen is recognized on the surface of professional antigen presenting cells (APCs). This is in contrast to direct presentation where antigen is seen on the surface of the target cells (e.g. infected cells or tumor cells). Here we introduce a mathematical model, which takes the phenomenon of cross-priming into account. We propose a new mechanism of regulation which is implicit in the dynamics of the CTL: According to the model, the ability of a CTL response to become established depends on the ratio of cross-presentation to direct presentation of the antigen. If this ratio is relatively high, CTL responses are likely to become established. If this ratio is relatively low, tolerance is the likely outcome. The behavior of the model includes a parameter region where the outcome depends on the initial conditions. We discuss our results with respect to the idea of self/non-self discrimination and the danger signal hypothesis. We apply the model to study the role of CTL in cancer initiation, cancer evolution/progression, and therapeutic vaccination against cancers. link: http://identifiers.org/pubmed/12706524

Memory is a central characteristic of immune responses. It is defined as an elevated number of specific immune cells that remain after resolution of infection and can protect the host against reinfection. The evolution of immunological memory is subject to debate. The advantages of memory discussed so far include protection from reinfection, control of chronic infection, and the transfer of immune function to the next generation. Mathematical models are used to identify a new force that can drive the evolution of immunological memory: the duration of memory can regulate the degree of competition between different pathogens. While a long duration of memory provides lasting protection against reinfection, it may also allow an inferior pathogen species to persist. This can be detrimental for the host if the inferior pathogen is more virulent. On the other hand, a shorter duration of memory ensures that an inferior pathogen species is excluded. This can be beneficial for the host if the inferior pathogen is more virulent. Thus, while in the absence of pathogen diversity memory is always expected to evolve to a long duration, under specific circumstances, memory can evolve toward shorter durations in the presence of pathogen diversity. link: http://identifiers.org/pubmed/13678598

Upon acute viral infection, a typical cytotoxic T lymphocyte (CTL) response is characterized by a phase of expansion and contraction after which it settles at a relatively stable memory level. Recently, experimental data from mice infected with murine cytomegalovirus (MCMV) showed different and unusual dynamics. After acute infection had resolved, some antigen specific CTL started to expand over time despite the fact that no replicative virus was detectable. This phenomenon has been termed as "CTL memory inflation". In order to examine the dynamics of this system further, we developed a mathematical model analysing the impact of innate and adaptive immune responses. According to this model, a potentially important contributor to CTL inflation is competition between the specific CTL response and an innate natural killer (NK) cell response. Inflation occurs most readily if the NK cell response is more efficient than the CTL at reducing virus load during acute infection, but thereafter maintains a chronic virus load which is sufficient to induce CTL proliferation. The model further suggests that weaker NK cell mediated protection can correlate with more pronounced CTL inflation dynamics over time. We present experimental data from mice infected with MCMV which are consistent with the theoretical predictions. This model provides valuable information and may help to explain the inflation of CMV specific CD8+T cells seen in humans as they age. link: http://identifiers.org/pubmed/17251133

Upon acute viral infection, a typical cytotoxic T lymphocyte (CTL) response is characterized by a phase of expansion and contraction after which it settles at a relatively stable memory level. Recently, experimental data from mice infected with murine cytomegalovirus (MCMV) showed different and unusual dynamics. After acute infection had resolved, some antigen specific CTL started to expand over time despite the fact that no replicative virus was detectable. This phenomenon has been termed as "CTL memory inflation". In order to examine the dynamics of this system further, we developed a mathematical model analysing the impact of innate and adaptive immune responses. According to this model, a potentially important contributor to CTL inflation is competition between the specific CTL response and an innate natural killer (NK) cell response. Inflation occurs most readily if the NK cell response is more efficient than the CTL at reducing virus load during acute infection, but thereafter maintains a chronic virus load which is sufficient to induce CTL proliferation. The model further suggests that weaker NK cell mediated protection can correlate with more pronounced CTL inflation dynamics over time. We present experimental data from mice infected with MCMV which are consistent with the theoretical predictions. This model provides valuable information and may help to explain the inflation of CMV specific CD8+T cells seen in humans as they age. link: http://identifiers.org/pubmed/17251133

Upon acute viral infection, a typical cytotoxic T lymphocyte (CTL) response is characterized by a phase of expansion and contraction after which it settles at a relatively stable memory level. Recently, experimental data from mice infected with murine cytomegalovirus (MCMV) showed different and unusual dynamics. After acute infection had resolved, some antigen specific CTL started to expand over time despite the fact that no replicative virus was detectable. This phenomenon has been termed as "CTL memory inflation". In order to examine the dynamics of this system further, we developed a mathematical model analysing the impact of innate and adaptive immune responses. According to this model, a potentially important contributor to CTL inflation is competition between the specific CTL response and an innate natural killer (NK) cell response. Inflation occurs most readily if the NK cell response is more efficient than the CTL at reducing virus load during acute infection, but thereafter maintains a chronic virus load which is sufficient to induce CTL proliferation. The model further suggests that weaker NK cell mediated protection can correlate with more pronounced CTL inflation dynamics over time. We present experimental data from mice infected with MCMV which are consistent with the theoretical predictions. This model provides valuable information and may help to explain the inflation of CMV specific CD8+T cells seen in humans as they age. link: http://identifiers.org/pubmed/17251133

Recent experimental data have shown that HIV-specific CD4 T cells provide a very important target for HIV replication. We use mathematical models to explore the effect of specific CD4 T cell infection on the dynamics of virus spread and immune responses. Infected CD4 T cells can provide antigen for their own stimulation. We show that such autocatalytic cell division can significantly enhance virus spread, and can also provide an additional reservoir for virus persistence during anti-viral drug therapy. In addition, the initial number of HIV-specific CD4 T cells is an important determinant of acute infection dynamics. A high initial number of HIV-specific CD4 T cells can lead to a sudden and fast drop of the population of HIV-specific CD4 T cells which results quickly in their extinction. On the other hand, a low initial number of HIV-specific CD4 T cells can lead to a prolonged persistence of HIV-specific CD4 T cell help at higher levels. The model suggests that boosting the population of HIV-specific CD4 T cells can increase the amount of virus-induced immune impairment, lead to less efficient anti-viral effector responses, and thus speed up disease progression, especially if effector responses such as CTL have not been sufficiently boosted at the same time. link: http://identifiers.org/pubmed/17379260

Tissues are maintained by adult stem cells that self-renew and also differentiate into functioning tissue cells. Homeostasis is achieved by a set of complex mechanisms that involve regulatory feedback loops. Similarly, tumors are believed to be maintained by a minority population of cancer stem cells, while the bulk of the tumor is made up of more differentiated cells, and there is indication that some of the feedback loops that operate in tissues continue to be functional in tumors. Mathematical models of such tissue hierarchies, including feedback loops, have been analyzed in a variety of different contexts. Apart from stem cells giving rise to differentiated cells, it has also been observed that more differentiated cells can de-differentiate into stem cells, both in healthy tissue and tumors, aspects of which have also been investigated mathematically. This paper analyses the effect of de-differentiation on the basic and evolutionary dynamics of cells in the context of tissue hierarchy models that include negative feedback regulation of the cell populations. The models predict that in the presence of de-differentiation, the fixation probability of a neutral mutant is lower than in its absence. Therefore, if de-differentiation occurs, a mutant with identical parameters compared to the wild-type cell population behaves like a disadvantageous mutant. Similarly, the process of de-differentiation is found to lower the fixation probability of an advantageous mutant. These results indicate that the presence of de-differentiation can lower the rates of tumor initiation and progression in the context of the models considered here. link: http://identifiers.org/pubmed/29605227

Tissues are maintained by adult stem cells that self-renew and also differentiate into functioning tissue cells. Homeostasis is achieved by a set of complex mechanisms that involve regulatory feedback loops. Similarly, tumors are believed to be maintained by a minority population of cancer stem cells, while the bulk of the tumor is made up of more differentiated cells, and there is indication that some of the feedback loops that operate in tissues continue to be functional in tumors. Mathematical models of such tissue hierarchies, including feedback loops, have been analyzed in a variety of different contexts. Apart from stem cells giving rise to differentiated cells, it has also been observed that more differentiated cells can de-differentiate into stem cells, both in healthy tissue and tumors, aspects of which have also been investigated mathematically. This paper analyses the effect of de-differentiation on the basic and evolutionary dynamics of cells in the context of tissue hierarchy models that include negative feedback regulation of the cell populations. The models predict that in the presence of de-differentiation, the fixation probability of a neutral mutant is lower than in its absence. Therefore, if de-differentiation occurs, a mutant with identical parameters compared to the wild-type cell population behaves like a disadvantageous mutant. Similarly, the process of de-differentiation is found to lower the fixation probability of an advantageous mutant. These results indicate that the presence of de-differentiation can lower the rates of tumor initiation and progression in the context of the models considered here. link: http://identifiers.org/pubmed/29605227

Mycoplasma pneumoniae, a threatening pathogen with a minimal genome, is a model organism for bacterial systems biology for which substantial experimental information is available. With the goal of understanding the complex interactions underlying its metabolism, we analyzed and characterized the metabolic network of M. pneumoniae in great detail, integrating data from different omics analyses under a range of conditions into a constraint-based model backbone. Iterating model predictions, hypothesis generation, experimental testing, and model refinement, we accurately curated the network and quantitatively explored the energy metabolism. In contrast to other bacteria, M. pneumoniae uses most of its energy for maintenance tasks instead of growth. We show that in highly linear networks the prediction of flux distributions for different growth times allows analysis of time-dependent changes, albeit using a static model. By performing an in silico knock-out study as well as analyzing flux distributions in single and double mutant phenotypes, we demonstrated that the model accurately represents the metabolism of M. pneumoniae. The experimentally validated model provides a solid basis for understanding its metabolic regulatory mechanisms. link: http://identifiers.org/pubmed/23549481

On the basis of a detailed model of yeast glycolysis, the effect of intercellular dynamics is analysed theoretically. The model includes the main steps of anaerobic glycolysis, and the production of ethanol and glycerol. Transmembrane diffusion of acetaldehyde is included, since it has been hypothesized that this substance mediates the interaction. Depending on the kinetic parameter, the single-cell model shows both stationary and oscillatory behaviour. This agrees with experimental data with respect to metabolite concentrations and phase shifts. The inclusion of intercellular coupling leads to a variety of dynamical modes, such as synchronous oscillations, and different kinds of asynchronous behavior. These oscillations can co-exist, leading to bi- and tri-rhythmicity. The corresponding parameter regions have been identified by a bifurcation analysis. The oscillatory dynamics of synchronized cell populations are investigated by calculating the phase responses to acetaldehyde pulses. Simulations are performed with respect to the synchronization of two subpopulations that are oscillating out of phase before mixing. The effect of the various process on synchronization is characterized quantitatively. While continuous exchange of acetaldehyde might synchronize the oscillations for appropriate sets of parameter values, the calculated synchronization time is longer than that observed experimentally. It is concluded either that addition to the transmembrane exchange of acetaldehyde, other processes may contribute to intercellular coupling, or that intracellular regulator feedback plays a role in the acceleration of the synchronization. for appropriate sets of parameter values, the calculated synchronization time is longer than that observed experimentally. It is concluded either that addition to the transmembrane exchange of acetaldehyde, other processes may contribute to intercellular coupling, or that intracellular regulator feedback plays a role in the acceleration of the synchronization. link: http://identifiers.org/pubmed/10702114

Under certain well-defined conditions, a population of yeast cells exhibits glycolytic oscillations that synchronize through intercellular acetaldehyde. This implies that the dynamic phenomenon of the oscillation propagates within and between cells. We here develop a method to establish by which route dynamics propagate through a biological reaction network. Application of the method to yeast demonstrates how the oscillations and the synchronization signal can be transduced. That transduction is not so much through the backbone of glycolysis, as via the Gibbs energy and redox coenzyme couples (ATP/ADP, and NADH/NAD), and via both intra- and intercellular acetaldehyde. link: http://identifiers.org/pubmed/10692304

Autonomous metabolic oscillations were observed in aerobic continuous culture of Saccharomyces cerevisiae. Experimental investigation of the underlying mechanism revealed that several pathways and regulatory couplings are involved. Here a hypothetical mechanism including the sulfate assimilation pathway, ethanol degradation and respiration is transformed into a mathematical model. Simulations confirm the ability of the model to produce limit cycle oscillations which reproduce most of the characteristic features of the system. link: http://identifiers.org/pubmed/11423122

BACKGROUND: Recent findings suggest that in pancreatic acinar cells stimulated with bile acid, a pro-apoptotic effect of reactive oxygen species (ROS) dominates their effect on necrosis and spreading of inflammation. The first effect presumably occurs via cytochrome C release from the inner mitochondrial membrane. A pro-necrotic effect - similar to the one of Ca2+ - can be strong opening of mitochondrial pores leading to breakdown of the membrane potential, ATP depletion, sustained Ca2+ increase and premature activation of digestive enzymes. To explain published data and to understand ROS effects during the onset of acute pancreatitis, a model using multi-valued logic is constructed. Formal concept analysis (FCA) is used to validate the model against data as well as to analyze and visualize rules that capture the dynamics. RESULTS: Simulations for two different levels of bile stimulation and for inhibition or addition of antioxidants reproduce the qualitative behaviour shown in the experiments. Based on reported differences of ROS production and of ROS induced pore opening, the model predicts a more uniform apoptosis/necrosis ratio for higher and lower bile stimulation in liver cells than in pancreatic acinar cells. FCA confirms that essential dynamical features of the data are captured by the model. For instance, high necrosis always occurs together with at least a medium level of apoptosis. At the same time, FCA helps to reveal subtle differences between data and simulations. The FCA visualization underlines the protective role of ROS against necrosis. CONCLUSIONS: The analysis of the model demonstrates how ROS and decreased antioxidant levels contribute to apoptosis. Studying the induction of necrosis via a sustained Ca2+ increase, we implemented the commonly accepted hypothesis of ATP depletion after strong bile stimulation. Using an alternative model, we demonstrate that this process is not necessary to generate the dynamics of the measured variables. Opening of plasma membrane channels could also lead to a prolonged increase of Ca2+ and to necrosis. Finally, the analysis of the model suggests a direct experimental testing for the model-based hypothesis of a self-enhancing cycle of cytochrome C release and ROS production by interruption of the mitochondrial electron transport chain. link: http://identifiers.org/pubmed/25315877

Metabolism is one of the attributes of life and supplies energy and building blocks to organisms. Therefore, understanding metabolism is crucial for the understanding of complex biological phenomena. Despite having been in the focus of research for centuries, our picture of metabolism is still incomplete. Metabolomics, the systematic analysis of all small molecules in a biological system, aims to close this gap. In order to facilitate such investigations a blueprint of the metabolic network is required. Recently, several metabolic network reconstructions for the model organism Caenorhabditis elegans have been published, each having unique features. We have established the WormJam Community to merge and reconcile these (and other unpublished models) into a single consensus metabolic reconstruction. In a series of workshops and annotation seminars this model was refined with manual correction of incorrect assignments, metabolite structure and identifier curation as well as addition of new pathways. The WormJam consensus metabolic reconstruction represents a rich data source not only for in silico network-based approaches like flux balance analysis, but also for metabolomics, as it includes a database of metabolites present in C. elegans, which can be used for annotation. Here we present the process of model merging, correction and curation and give a detailed overview of the model. In the future it is intended to expand the model toward different tissues and put special emphasizes on lipid metabolism and secondary metabolism including ascaroside metabolism in accordance to their central role in C. elegans physiology. link: http://identifiers.org/pubmed/30488036

Data from (31)P-nuclear magnetic resonance spectroscopy of human forearm flexor muscle were analyzed based on a previously developed model of mitochondrial oxidative phosphorylation (PLoS Comp Bio 1: e36, 2005) to test the hypothesis that substrate level (concentrations of ADP and inorganic phosphate) represents the primary signal governing the rate of mitochondrial ATP synthesis and maintaining the cellular ATP hydrolysis potential in skeletal muscle. Model-based predictions of cytoplasmic concentrations of phosphate metabolites (ATP, ADP, and P(i)) matched data obtained from 20 healthy volunteers and indicated that as work rate is varied from rest to submaximal exercise commensurate increases in the rate of mitochondrial ATP synthesis are effected by changes in concentrations of available ADP and P(i). Additional data from patients with a defect of complex I of the respiratory chain and a patient with a deficiency in the mitochondrial adenine nucleotide translocase were also predicted the by the model by making the appropriate adjustments to the activities of the affected proteins associates with the defects, providing both further validation of the biophysical model of the control of oxidative phosphorylation and insight into the impact of these diseases on the ability of the cell to maintain its energetic state. link: http://identifiers.org/pubmed/16837647

The ATP-sensitive K(+) (K(ATP)) channels are composed of sulfonylurea receptor and inwardly rectifying K(+) channel (Kir6.2) subunit. These channels are regulated by intracellular ADP/ATP ratio and play a role in cellular metabolism. Diethyl pyrocarbonate (DEPC), a histidine-specific alkylating reagent, is known to modify the histidine residues of the structure of proteins. The objective of this study was to determine whether DEPC modifies K(ATP)-channel activity in pituitary GH(3) cells. Steady-state fluctuation analyses of macroscopic K(+) current at -120 mV produced power spectra that could be fitted with a single Lorentzian curve in these cells. The time constants in the presence of DEPC were increased. Consistent with fluctuation analyses, the mean open time of K(ATP)-channels was significantly increased during exposure to DEPC. However, DEPC produced no change in single-channel conductance, despite the ability of this compound to enhance K(ATP)-channel activity in a concentration-dependent manner with an EC(50) value of 16 microM. DEPC-stimulated K(ATP)-channel activity was attenuated by pretreatment with glibenclamide. In current-clamp configuration, DEPC decreased the firing of action potentials in GH(3) cells. A further application of glibenclamide reversed DEPC-induced inhibition of spontaneous action potentials. Intracellullar Ca(2+) measurements revealed the ability of DEPC to decrease Ca(2+) oscillations in GH(3) cells. Simulation studies also demonstrated that the increased conductance of K(ATP)-channels used to mimic DEPC actions reduced the frequency of spontaneous action potentials and fluctuation of intracellular Ca(2+). The results indicate that chemical modification with DEPC enhances K(ATP)-channel activity and influences functional activities of pituitary GH(3) cells. link: http://identifiers.org/pubmed/16375866

A computational model of mitochondrial metabolism and electrophysiology is introduced and applied to analysis of data from isolated cardiac mitochondria and data on phosphate metabolites in striated muscle in vivo. This model is constructed based on detailed kinetics and thermodynamically balanced reaction mechanisms and a strict accounting of rapidly equilibrating biochemical species. Since building such a model requires introducing a large number of adjustable kinetic parameters, a correspondingly large amount of independent data from isolated mitochondria respiring on different substrates and subject to a variety of protocols is used to parameterize the model and ensure that it is challenged by a wide range of data corresponding to diverse conditions. The developed model is further validated by both in vitro data on isolated cardiac mitochondria and in vivo experimental measurements on human skeletal muscle. The validated model is used to predict the roles of NAD and ADP in regulating the tricarboxylic acid cycle dehydrogenase fluxes, demonstrating that NAD is the more important regulator. Further model predictions reveal that a decrease of cytosolic pH value results in decreases in mitochondrial membrane potential and a corresponding drop in the ability of the mitochondria to synthesize ATP at the hydrolysis potential required for cellular function. link: http://identifiers.org/pubmed/17591785

Our objective was to determine the pharmacokinetics, bioavailability and lymph node uptake of the monoclonal antibody bevacizumab, labeled with the near-infrared (IR) dye 800CW, after intravenous (IV) and subcutaneous (SC) administration in mice. Fluorescence imaging and enzyme-linked immunosorbent assay (ELISA) assays were developed and validated to measure the concentration of bevacizumab in plasma. The bevacizumab-IRDye conjugate remained predominantly intact in plasma and in lymph node homogenate samples over a 24-h period, as determined by sodium dodecyl sulfate polyacrylamide gel electrophoresis and size exclusion chromatography. The plasma concentration vs. time plots obtained by fluorescence and ELISA measurements were similar; however, unlike ELISA, fluorescent imaging was only able to quantitate concentrations for 24 h after administration. At a low dose of 0.45 mg/kg, the plasma clearance of bevacizumab was 6.96 mL/h/kg after IV administration; this clearance is higher than that reported after higher doses. Half-lives of bevacizumab after SC and IV administration were 4.6 and 3.9 days, respectively. After SC administration, bevacizumab-IRDye800CW was present in the axillary lymph nodes that drain the SC site; lymph node uptake of bevacizumab-IRDye 800CW was negligible after IV administration. Bevacizumab exhibited complete bioavailability after SC administration. Using a compartmental pharmacokinetic model, the fraction absorbed through the lymphatics after SC administration was estimated to be about 1%. This is the first report evaluating the use of fluorescent imaging to determine the pharmacokinetics, lymphatic uptake, and bioavailability of a near-infrared dye-labeled antibody conjugate. link: http://identifiers.org/pubmed/22391791

Circadian rhythms of gene activity, metabolism, physiology and behaviour are observed in all the eukaryotes and some prokaryotes. In this study, we present a model to represent the transcriptional regulatory network essential for the circadian rhythmicity in Drosophila. The model incorporates the transcriptional feedback loops revealed so far in the network of the circadian clock (PER/TIM and VRI/PDP1 loops). Conventional Hill functions are not assumed to describe the regulation of genes, instead of the explicit reactions of binding and unbinding processes of transcription factors to promoters are modelled. The model simulates sustained circadian oscillations in mRNA and protein concentrations in constant darkness in agreement with experimental observations. It also simulates entrainment by light-dark cycles, disappearance of the rhythmicity in constant light and the shape of phase response curves resembling that of the experimental results. The model is robust over a wide range of parameter variations. In addition, the simulated E-box mutation, per(S) and per(L) mutants are similar to that observed in the experiments. The deficiency between the simulated mRNA levels and experimental observations in per(01), tim(01) and clk(Jrk) mutants suggests some difference on the part of the model from reality. link: http://identifiers.org/pubmed/17157878

In order to confirm which process is the most important in the blood coagulation cascade, a dynamic model of the function of platelets in blood coagulation is provided based on biochemical experiments. A series of conclusions based on qualitative analysis and mathematical simulation are drawn about the influence of the activation rate of factor VIII and factor IX on the generation of thrombin (IIa). It is evident that the pro-coagulation stimulus must exceed a threshold value to initiate the coagulation cascade. The value is related to the rate of platelet activation, the binding constant d2. The stability of the fixed value is also related to the pro-coagulation stimulus. This article also evaluates the influence of the stimulus strength and the activated rate parameter of platelets on thrombin. The proportion of platelets activated at any given time is designated c. To each c, we obtain a maximum concentration of thrombin. It is evident that when the level of factor IX is below 1% of normal levels, the rate of thrombin generation reduces dramatically resulting in severe bleeding tendency. link: http://identifiers.org/pubmed/12376045

We studied the bradykinin-induced changes in phosphoinositide composition of N1E-115 neuroblastoma cells using a combination of biochemistry, microscope imaging, and mathematical modeling. Phosphatidylinositol-4,5-bisphosphate (PIP2) decreased over the first 30 s, and then recovered over the following 2-3 min. However, the rate and amount of inositol-1,4,5-trisphosphate (InsP3) production were much greater than the rate or amount of PIP2 decline. A mathematical model of phosphoinositide turnover based on this data predicted that PIP2 synthesis is also stimulated by bradykinin, causing an early transient increase in its concentration. This was subsequently confirmed experimentally. Then, we used single-cell microscopy to further examine phosphoinositide turnover by following the translocation of the pleckstrin homology domain of PLCdelta1 fused to green fluorescent protein (PH-GFP). The observed time course could be simulated by incorporating binding of PIP2 and InsP3 to PH-GFP into the model that had been used to analyze the biochemistry. Furthermore, this analysis could help to resolve a controversy over whether the translocation of PH-GFP from membrane to cytosol is due to a decrease in PIP2 on the membrane or an increase in InsP3 in cytosol; by computationally clamping the concentrations of each of these compounds, the model shows how both contribute to the dynamics of probe translocation. link: http://identifiers.org/pubmed/12771127

Blood coagulation is a complex physiological process. With the development of biochemistry and molecular biology, especially the discovery of the tissue factor pathway inhibitor (TFPI) and further research, the researches of tissue factor (TF) and TFPI have been the hot topics now. Based on this result, we provide a new concept of blood coagulation system and provide a dynamic model to describe the role of TFPI of this pathway. Through stability analysis of the equilibrium point and numerical simulations, we obtain some important conclusions. All these results may have great significance for lesions of cardiac vessels and cerebral vessels, haemophilia, DIC, thrombosis formation and malignancy. link: http://identifiers.org/pubmed/15767116

The cellular and molecular mechanisms of immunomodulatory actions of glucocorticoids (GC) remain to be identified. Using our experimental findings, a mathematical model based on a system of ordinary differential equations for characterization of the regulation of anti-tumor immune activity by the both direct and indirect GC effects was generated to study the effects of GC treatment on effector CD8+ T cells, GC-generated tolerogenic dendritic cells (DC), regulatory T cells and the growth of lymphoma cells. In addition, we compared the data from in vivo and in silico experiments. The mathematical simulations indicated that treatment with GCs may suppress anti-tumor immune response in a dose-dependent manner. Our in silico results were in line with our earlier experimental findings of inhibitory effects of GCs on T, NK and dendritic cells. The model simulations were in line with earlier experimental observations of inhibitory effects of GCs on T and NK cells and DCs. The results of this study might be useful for predicting clinical outcomes in patients receiving GC therapy. link: http://identifiers.org/pubmed/31734335

Suppressor of cytokine signaling-1 (SOCS1) was identified as the negative regulator of Janus kinase (JAK) and signal transducer and activator of transcription (STAT) signal transduction pathway. However, the kinetics and control mechanism of the pathway have not yet been fully understood. We have developed the computer simulation of the JAK/STAT pathway. Without nuclear phosphatase, SOCS1's binding to JAK did not cause the decrease in nuclear phosphorylated STAT1. However, without SH2 domain-containing tyrosine phosphatase 2 (SHP-2) or cytoplasmic phosphatase, it did. So nuclear phosphatase is considered to be the most important in this system. By changing parameters of the model, dynamical characteristics and control mechanism were investigated. link: http://identifiers.org/pubmed/12527385

Suppressor of cytokine signaling-1 (SOCS1) was identified as the negative regulator of Janus kinase (JAK) and signal transducer and activator of transcription (STAT) signal transduction pathway. However, the kinetics and control mechanism of the pathway have not yet been fully understood. We have developed the computer simulation of the JAK/STAT pathway. Without nuclear phosphatase, SOCS1's binding to JAK did not cause the decrease in nuclear phosphorylated STAT1. However, without SH2 domain-containing tyrosine phosphatase 2 (SHP-2) or cytoplasmic phosphatase, it did. So nuclear phosphatase is considered to be the most important in this system. By changing parameters of the model, dynamical characteristics and control mechanism were investigated. link: http://identifiers.org/pubmed/12527385

We simulated myocardial Ca2+ (Ca) and cross-bridge (CB) kinetics to get insight into the experimentally observed constancy and variability of cardiac contractile efficiency in generating total mechanical energy under various inotropic and pathological conditions. The simulation consisted of a Ca transient, Ca association and dissociation rate constants of troponin C, and CB on and off rate constants. We evaluated sarcomere isometric twitch contractions at a constant muscle length. We assumed that each CB cycle hydrolyzes one ATP and that the force-length area (FLA) quantifies the total mechanical energy generated by CB cycles in a twitch contraction. FLA is a linear version of pressure-volume area, which quantifies the total mechanical energy of cardiac twitch contraction and correlates linearly with cardiac oxygen consumption (H. Suga, Physiol. Rev. 70: 247-277, 1990). The simulation shows that the contractile efficiency varies with changes in the Ca transient and Ca and CB kinetics except when they simultaneously speed up or slow down proportionally. These results point to possible mechanisms underlying the constancy and variability of cardiac contractile efficiency. link: http://identifiers.org/pubmed/8967394

MiRNAs, which are a family of small non-coding RNAs, regulate a broad array of physiological and developmental processes. However, their regulatory roles have remained largely mysterious. E2F is a positive regulator of cell cycle progression and also a potent inducer of apoptosis. Positive feedback loops in the regulation of Rb-E2F pathway are predicted and shown experimentally. Recently, it has been discovered that E2F induce a cluster of miRNAs called miR449. In turn, E2F is inhibited by miR449 through regulating different transcripts, thus forming negative feedback loops in the interaction network. Here, based on the integration of experimental evidence and quantitative data, we studied Rb-E2F pathway coupling the positive feedback loops and negative feedback loops mediated by miR449. Therefore, a mathematical model is constructed based in part on the model proposed in Yao-Lee et al. (2008) and nonlinear dynamical behaviors including the stability and bifurcations of the model are discussed. A comparison is given to reveal the implication of the fundamental differences of Rb-E2F pathway between regulation and deregulation of miR449. Coherent with the experiments it predicts that miR449 plays a critical role in regulating the cell cycle progression and provides a twofold safety mechanism to avoid excessive E2F-induced proliferation by cell cycle arrest and apoptosis. Moreover, numerical simulation and bifurcation analysis shows that the mechanisms of the negative regulation of miR449 to three different transcripts are quite distinctive which needs to be verified experimentally. This study may help us to analyze the whole cell cycle process mediated by other miRNAs more easily. A better knowledge of the dynamical behaviors of miRNAs mediated networks is also of interest for bio-engineering and artificial control. link: http://identifiers.org/pubmed/23028477

The pacemaker activity of the S-A node cell was explained by reconstructing the pacemaker potential using a Hodgkin-Huxley type mathematical model which was based on the reported voltage clamp data. In this model four dynamic currents, the sodium current iNa, the slow inward current, is, the potassium current, iK, and the hyperpolarization-activated current, ih, in addition to a time-dependent leak current, i1 were included. The model simulated the spontaneous action potential the current voltage relationship, and the voltage clamp experiment in normal Tyrode solution of the rabbit S-A node. Furthermore, the changes of activity induced by the potassium current blocker Ba2+, by applying constant current, acetylcholine, and epinephrine were also reconstructed. It was strongly suggested that the pacemaker depolarization in the S-A node cell is mainly due to a gradual increase of iS during diastole, and that the contribution of iK is much less compared to the potassium current iK2 in the Purkinje fiber pacemaker depolarization. The rising phase of the action potential was due to iS and the plateau phase is determined by both the inactivation of iS and activation of iK. link: http://identifiers.org/pubmed/7265560

Limited pharmacokinetic (PK) and pharmacodynamic (PD) data are available to use in methadone dosing recommendations in pediatric patients for either opioid abstinence or analgesia. Considering the extreme inter-individual variability of absorption and metabolism of methadone, population-based PK would be useful to provide insight into the relationship between dose, blood concentrations, and clinical effects of methadone. To address this need, an age-dependent physiologically based pharmacokinetic (PBPK) model has been constructed to systematically study methadone metabolism and PK. The model will facilitate the design of cost-effective studies that will evaluate methadone PK and PD relationships, and may be useful to guide methadone dosing in children. The PBPK model, which includes whole-body multi-organ distribution, plasma protein binding, metabolism, and clearance, is parameterized based on a database of pediatric PK parameters and data collected from clinical experiments. The model is further tailored and verified based on PK data from individual adults, then scaled appropriately to apply to children aged 0-24 months. Based on measured variability in CYP3A enzyme expression levels and plasma orosomucoid (ORM2) concentrations, a Monte-Carlo-based simulation of methadone kinetics in a pediatric population was performed. The simulation predicts extreme variability in plasma concentrations and clearance kinetics for methadone in the pediatric population, based on standard dosing protocols. In addition, it is shown that when doses are designed for individuals based on prior protein expression information, inter-individual variability in methadone kinetics may be greatly reduced. link: http://identifiers.org/pubmed/16758333

Drug molecules not only interact with specific targets, but also alter the state and function of the associated biological network. How to design drugs and evaluate their functions at the systems level becomes a key issue in highly efficient and low-side-effect drug design. The arachidonic acid metabolic network is the network that produces inflammatory mediators, in which several enzymes, including cyclooxygenase-2 (COX-2), have been used as targets for anti-inflammatory drugs. However, neither the century-old nonsteriodal anti-inflammatory drugs nor the recently revocatory Vioxx have provided completely successful anti-inflammatory treatment. To gain more insights into the anti-inflammatory drug design, the authors have studied the dynamic properties of arachidonic acid (AA) metabolic network in human polymorphous leukocytes. Metabolic flux, exogenous AA effects, and drug efficacy have been analyzed using ordinary differential equations. The flux balance in the AA network was found to be important for efficient and safe drug design. When only the 5-lipoxygenase (5-LOX) inhibitor was used, the flux of the COX-2 pathway was increased significantly, showing that a single functional inhibitor cannot effectively control the production of inflammatory mediators. When both COX-2 and 5-LOX were blocked, the production of inflammatory mediators could be completely shut off. The authors have also investigated the differences between a dual-functional COX-2 and 5-LOX inhibitor and a mixture of these two types of inhibitors. Their work provides an example for the integration of systems biology and drug discovery. link: http://identifiers.org/pubmed/17381237

Drugs against multiple targets may overcome the many limitations of single targets and achieve a more effective and safer control of the disease. Numerous high-throughput experiments have been performed in this emerging field. However, systematic identification of multiple drug targets and their best intervention requires knowledge of the underlying disease network and calls for innovative computational methods that exploit the network structure and dynamics. Here, we develop a robust computational algorithm for finding multiple target optimal intervention (MTOI) solutions in a disease network. MTOI identifies potential drug targets and suggests optimal combinations of the target intervention that best restore the network to a normal state, which can be customer designed. We applied MTOI to an inflammation-related network. The well-known side effects of the traditional non-steriodal anti-inflammatory drugs and the recently recalled Vioxx were correctly accounted for in our network model. A number of promising MTOI solutions were found to be both effective and safer. link: http://identifiers.org/pubmed/18985027

Drugs against multiple targets may overcome the many limitations of single targets and achieve a more effective and safer control of the disease. Numerous high-throughput experiments have been performed in this emerging field. However, systematic identification of multiple drug targets and their best intervention requires knowledge of the underlying disease network and calls for innovative computational methods that exploit the network structure and dynamics. Here, we develop a robust computational algorithm for finding multiple target optimal intervention (MTOI) solutions in a disease network. MTOI identifies potential drug targets and suggests optimal combinations of the target intervention that best restore the network to a normal state, which can be customer designed. We applied MTOI to an inflammation-related network. The well-known side effects of the traditional non-steriodal anti-inflammatory drugs and the recently recalled Vioxx were correctly accounted for in our network model. A number of promising MTOI solutions were found to be both effective and safer. link: http://identifiers.org/pubmed/18985027

Drugs against multiple targets may overcome the many limitations of single targets and achieve a more effective and safer control of the disease. Numerous high-throughput experiments have been performed in this emerging field. However, systematic identification of multiple drug targets and their best intervention requires knowledge of the underlying disease network and calls for innovative computational methods that exploit the network structure and dynamics. Here, we develop a robust computational algorithm for finding multiple target optimal intervention (MTOI) solutions in a disease network. MTOI identifies potential drug targets and suggests optimal combinations of the target intervention that best restore the network to a normal state, which can be customer designed. We applied MTOI to an inflammation-related network. The well-known side effects of the traditional non-steriodal anti-inflammatory drugs and the recently recalled Vioxx were correctly accounted for in our network model. A number of promising MTOI solutions were found to be both effective and safer. link: http://identifiers.org/pubmed/18985027

BACKGROUND: Cancer arises when within a single cell multiple malfunctions of control systems occur, which are, broadly, the system that promote cell growth and the system that protect against erratic growth. Additional systems within the cell must be corrupted so that a cancer cell, to form a mass of any real size, produces substances that promote the growth of new blood vessels. Multiple mutations are required before a normal cell can become a cancer cell by corruption of multiple growth-promoting systems. METHODS: We develop a simple mathematical model to describe the solid cancer growth dynamics inducing angiogenesis in the absence of cancer controlling mechanisms. RESULTS: The initial conditions supplied to the dynamical system consist of a perturbation in form of pulse: The origin of cancer cells from normal cells of an organ of human body. Thresholds of interacting parameters were obtained from the steady states analysis. The existence of two equilibrium points determine the strong dependency of dynamical trajectories on the initial conditions. The thresholds can be used to control cancer. CONCLUSIONS: Cancer can be settled in an organ if the following combination matches: better fitness of cancer cells, decrease in the efficiency of the repairing systems, increase in the capacity of sprouting from existing vascularization, and higher capacity of mounting up new vascularization. However, we show that cancer is rarely induced in organs (or tissues) displaying an efficient (numerically and functionally) reparative or regenerative mechanism. link: http://identifiers.org/pubmed/22300422

The restriction point (R-point) marks the critical event when a mammalian cell commits to proliferation and becomes independent of growth stimulation. It is fundamental for normal differentiation and tissue homeostasis, and seems to be dysregulated in virtually all cancers. Although the R-point has been linked to various activities involved in the regulation of G1-S transition of the mammalian cell cycle, the underlying mechanism remains unclear. Using single-cell measurements, we show here that the Rb-E2F pathway functions as a bistable switch to convert graded serum inputs into all-or-none E2F responses. Once turned ON by sufficient serum stimulation, E2F can memorize and maintain this ON state independently of continuous serum stimulation. We further show that, at critical concentrations and duration of serum stimulation, bistable E2F activation correlates directly with the ability of a cell to traverse the R-point. link: http://identifiers.org/pubmed/18364697

The heterogeneity in mammalian cells signaling response is largely a result of preexisting cell to cell variability. It is unknown whether cell to cell variability rises from biochemical stochastic fluctuations or distinct cellular states. Here, we utilize calcium response to adenosine trisphosphate as a model for investigating the structure of heterogeneity within a population of cells and analyze whether distinct cellular response states coexist. We use a functional definition of cellular state that is based on a mechanistic dynamical systems model of calcium signaling. Using Bayesian parameter inference, we obtain high confidence parameter value distributions for several hundred cells, each fitted individually. Clustering the inferred parameter distributions revealed three major distinct cellular states within the population. The existence of distinct cellular states raises the possibility that the observed variability in response is a result of structured heterogeneity between cells. The inferred parameter distribution predicts, and experiments confirm that variability in IP3R response explains the majority of calcium heterogeneity. Our work shows how mechanistic models and single cell parameter fitting can uncover hidden population structure and demonstrate the need for parameter inference at the single cell level. link: http://identifiers.org/doi/10.15252/msb.20167137

Brief dopamine events are critical actors of reward-mediated learning in the striatum; the intracellular cAMP-protein kinase A (PKA) response of striatal medium spiny neurons to such events was studied dynamically using a combination of biosensor imaging in mouse brain slices and in silico simulations. Both D1 and D2 medium spiny neurons can sense brief dopamine transients in the sub-micromolar range. While dopamine transients profoundly change cAMP levels in both types of medium spiny neurons, the PKA-dependent phosphorylation level remains unaffected in D2 neurons. At the level of PKA-dependent phosphorylation, D2 unresponsiveness depends on protein phosphatase-1 (PP1) inhibition by DARPP-32. Simulations suggest that D2 medium spiny neurons could detect transient dips in dopamine level.The phasic release of dopamine in the striatum determines various aspects of reward and action selection, but the dynamics of the dopamine effect on intracellular signalling remains poorly understood. We used genetically encoded FRET biosensors in striatal brain slices to quantify the effect of transient dopamine on cAMP or PKA-dependent phosphorylation levels, and computational modelling to further explore the dynamics of this signalling pathway. Medium-sized spiny neurons (MSNs), which express either D1 or D2 dopamine receptors, responded to dopamine by an increase or a decrease in cAMP, respectively. Transient dopamine showed similar sub-micromolar efficacies on cAMP in both D1 and D2 MSNs, thus challenging the commonly accepted notion that dopamine efficacy is much higher on D2 than on D1 receptors. However, in D2 MSNs, the large decrease in cAMP level triggered by transient dopamine did not translate to a decrease in PKA-dependent phosphorylation level, owing to the efficient inhibition of protein phosphatase 1 by DARPP-32. Simulations further suggested that D2 MSNs can also operate in a 'tone-sensing' mode, allowing them to detect transient dips in basal dopamine. Overall, our results show that D2 MSNs may sense much more complex patterns of dopamine than previously thought. link: http://identifiers.org/pubmed/28782235

Brief dopamine events are critical actors of reward-mediated learning in the striatum; the intracellular cAMP-protein kinase A (PKA) response of striatal medium spiny neurons to such events was studied dynamically using a combination of biosensor imaging in mouse brain slices and in silico simulations. Both D1 and D2 medium spiny neurons can sense brief dopamine transients in the sub-micromolar range. While dopamine transients profoundly change cAMP levels in both types of medium spiny neurons, the PKA-dependent phosphorylation level remains unaffected in D2 neurons. At the level of PKA-dependent phosphorylation, D2 unresponsiveness depends on protein phosphatase-1 (PP1) inhibition by DARPP-32. Simulations suggest that D2 medium spiny neurons could detect transient dips in dopamine level.The phasic release of dopamine in the striatum determines various aspects of reward and action selection, but the dynamics of the dopamine effect on intracellular signalling remains poorly understood. We used genetically encoded FRET biosensors in striatal brain slices to quantify the effect of transient dopamine on cAMP or PKA-dependent phosphorylation levels, and computational modelling to further explore the dynamics of this signalling pathway. Medium-sized spiny neurons (MSNs), which express either D1 or D2 dopamine receptors, responded to dopamine by an increase or a decrease in cAMP, respectively. Transient dopamine showed similar sub-micromolar efficacies on cAMP in both D1 and D2 MSNs, thus challenging the commonly accepted notion that dopamine efficacy is much higher on D2 than on D1 receptors. However, in D2 MSNs, the large decrease in cAMP level triggered by transient dopamine did not translate to a decrease in PKA-dependent phosphorylation level, owing to the efficient inhibition of protein phosphatase 1 by DARPP-32. Simulations further suggested that D2 MSNs can also operate in a 'tone-sensing' mode, allowing them to detect transient dips in basal dopamine. Overall, our results show that D2 MSNs may sense much more complex patterns of dopamine than previously thought. link: http://identifiers.org/pubmed/28782235

The asymptomatic phase of HIV infection is characterised by a slow decline of peripheral blood CD4(+) T cells. Why this decline is slow is not understood. One potential explanation is that the low average rate of homeostatic proliferation or immune activation dictates the pace of a "runaway" decline of memory CD4(+) T cells, in which activation drives infection, higher viral loads, more recruitment of cells into an activated state, and further infection events. We explore this hypothesis using mathematical models.Using simple mathematical models of the dynamics of T cell homeostasis and proliferation, we find that this mechanism fails to explain the time scale of CD4(+) memory T cell loss. Instead it predicts the rapid attainment of a stable set point, so other mechanisms must be invoked to explain the slow decline in CD4(+) cells.A runaway cycle in which elevated CD4(+) T cell activation and proliferation drive HIV production and vice versa cannot explain the pace of depletion during chronic HIV infection. We summarize some alternative mechanisms by which the CD4(+) memory T cell homeostatic set point might slowly diminish. While none are mutually exclusive, the phenomenon of viral rebound, in which interruption of antiretroviral therapy causes a rapid return to pretreatment viral load and T cell counts, supports the model of virus adaptation as a major force driving depletion. link: http://identifiers.org/pubmed/17518516

BACKGROUND AND OBJECTIVES:In recent decades, cancer has become one of the most fatal and destructive diseases which is threatening humans life. Accordingly, different types of cancer treatment are studied with the main aim to have the best treatment with minimum side effects. Anti-angiogenic is a molecular targeted therapy which can be coupled with chemotherapy and radiotherapy. Although this method does not eliminate the whole tumor, but it can keep the tumor size in a given state by preventing the formation of new blood vessels. In this paper, a novel model-free method based on reinforcement learning (RL) framework is used to design a closed-loop control of anti-angiogenic drug dosing administration. METHODS:A Q-learning algorithm is developed for the drug dosing closed-loop control. This controller is designed using two different values of the maximum drug dosage to reduce the tumor volume up to a desired value. The mathematical model of tumor growth under anti-angiogenic inhibitor is used to simulate a real patient. RESULTS:The effectiveness of the proposed method is shown through in silico simulation and its robustness to patient parameters variation is demonstrated. It is demonstrated that the tumor reaches its minimal volume in 84 days with maximum drug inlet of 30 mg/kg/day. Also, it is shown that the designed controller is robust with respect to  ± 20% of tumor growth parameters changes. CONCLUSION:The proposed closed-loop reinforcement learning-based controller for cancer treatment using anti-angiogenic inhibitor provides an effective and novel result such that with a clinically valid and safe dosage of drug, the volume reduces up to 1mm3 in a reasonable short period compared to the literature. link: http://identifiers.org/pubmed/31046990

Chinese hamster ovary (CHO) cells are most prevalently used for producing recombinant therapeutics in biomanufacturing. Recently, more rational and systems approaches have been increasingly exploited to identify key metabolic bottlenecks and engineering targets for cell line engineering and process development based on the CHO genome-scale metabolic model which mechanistically characterizes cell culture behaviours. However, it is still challenging to quantify plausible intracellular fluxes and discern metabolic pathway usages considering various clonal traits and bioprocessing conditions. Thus, we newly incorporated enzyme kinetic information into the updated CHO genome-scale model (iCHO2291) and added enzyme capacity constraints within the flux balance analysis framework (ecFBA) to significantly reduce the flux variability in biologically meaningful manner, as such improving the accuracy of intracellular flux prediction. Interestingly, ecFBA could capture the overflow metabolism under the glucose excess condition where the usage of oxidative phosphorylation is limited by the enzyme capacity. In addition, its applicability was successfully demonstrated via a case study where the clone- and media-specific lactate metabolism was deciphered, suggesting that the lactate-pyruvate cycling could be beneficial for CHO cells to efficiently utilize the mitochondrial redox capacity. In summary, iCHO2296 with ecFBA can be used to confidently elucidate cell cultures and effectively identify key engineering targets, thus guiding bioprocess optimization and cell engineering efforts as a part of digital twin model for advanced biomanufacturing in future. link: http://identifiers.org/pubmed/32330653

The yeast mating response is one of the best understood heterotrimeric G protein signaling pathways. Yet, most descriptions of this system have been qualitative. We have quantitatively characterized the heterotrimeric G protein cycle in yeast based on direct in vivo measurements. We used fluorescence resonance energy transfer to monitor the association state of cyan fluorescent protein (CFP)-Galpha and Gbetagamma-yellow fluorescent protein (YFP), and we found that receptor-mediated G protein activation produced a loss of fluorescence resonance energy transfer. Quantitative time course and dose-response data were obtained for both wild-type and mutant cells possessing an altered pheromone response. These results paint a quantitative portrait of how regulators such as Sst2p and the C-terminal tail of alpha-factor receptor modulate the kinetics and sensitivity of G protein signaling. We have explored critical features of the dynamics including the rapid rise and subsequent decline of active G proteins during the early response, and the relationship between the G protein activation dose-response curve and the downstream dose-response curves for cell-cycle arrest and transcriptional induction. Fitting the data to a mathematical model produced estimates of the in vivo rates of heterotrimeric G protein activation and deactivation in yeast. link: http://identifiers.org/pubmed/12960402

A mathematical model for the regulation of induction in the lac operon in Escherichia coli is presented. This model takes into account the dynamics of the permease facilitating the internalization of external lactose; internal lactose; beta-galactosidase, which is involved in the conversion of lactose to allolactose, glucose and galactose; the allolactose interactions with the lac repressor; and mRNA. The final model consists of five nonlinear differential delay equations with delays due to the transcription and translation process. We have paid particular attention to the estimation of the parameters in the model. We have tested our model against two sets of beta-galactosidase activity versus time data, as well as a set of data on beta-galactosidase activity during periodic phosphate feeding. In all three cases we find excellent agreement between the data and the model predictions. Analytical and numerical studies also indicate that for physiologically realistic values of the external lactose and the bacterial growth rate, a regime exists where there may be bistable steady-state behavior, and that this corresponds to a cusp bifurcation in the model dynamics. link: http://identifiers.org/pubmed/12719218

Caenorhabditis elegans is a powerful model to study metabolism and how it relates to nutrition, gene expression, and life history traits. However, while numerous experimental techniques that enable perturbation of its diet and gene function are available, a high-quality metabolic network model has been lacking. Here, we reconstruct an initial version of the C. elegans metabolic network. This network model contains 1,273 genes, 623 enzymes, and 1,985 metabolic reactions and is referred to as iCEL1273. Using flux balance analysis, we show that iCEL1273 is capable of representing the conversion of bacterial biomass into C. elegans biomass during growth and enables the predictions of gene essentiality and other phenotypes. In addition, we demonstrate that gene expression data can be integrated with the model by comparing metabolic rewiring in dauer animals versus growing larvae. iCEL1273 is available at a dedicated website (wormflux.umassmed.edu) and will enable the unraveling of the mechanisms by which different macro- and micronutrients contribute to the animal's physiology. link: http://identifiers.org/pubmed/27211857

Metabolism is a highly compartmentalized process that provides building blocks for biomass generation during development, homeostasis and wound healing, and energy to support cellular and organismal processes. In metazoans, different cells and tissues specialize in different aspects of metabolism. However, studying the compartmentalization of metabolism in different cell types in a whole animal and for a particular stage of life is difficult. Here, we present MERGE (MEtabolic models Reconciled with Gene Expression), a computational pipeline that we used to predict tissue-relevant metabolic function at the network-, pathway-, reaction-, and metabolite-level based on single-cell RNA-sequencing (scRNA-seq) data from the nematode Caenorhabditis elegans. Our analysis recapitulated known tissue functions in C. elegans, captured metabolic properties that are shared with similar tissues in human, and also provided predictions for novel metabolic functions. MERGE is versatile and applicable to other systems. We envision this work as a starting point for the development of metabolic network models for individual cells as scRNA-seq continues to provide higher resolution gene expression data. link: http://identifiers.org/doi/10.15252/msb.20209649

The transferrin (Tf) trafficking pathway is a promising mechanism for use in targeted cancer therapy due to the overexpression of transferrin receptors (TfRs) on cancerous cells. We have previously developed a mathematical model of the Tf/TfR trafficking pathway to improve the efficiency of Tf as a drug carrier. By using diphtheria toxin (DT) as a model toxin, we found that mutating the Tf protein to change its iron release rate improves cellular association and efficacy of the drug. Though this is an improvement upon using wild-type Tf as the targeting ligand, conjugated toxins like DT are unfortunately still highly cytotoxic at off-target sites. In this work, we address this hurdle in cancer research by developing a mathematical model to predict the efficacy and selectivity of Tf conjugates that use an alternative toxin. For this purpose, we have chosen to study a mutant of DT, cross-reacting material 107 (CRM107). First, we developed a mathematical model of the Tf-DT trafficking pathway by extending our Tf/TfR model to include intracellular trafficking via DT and DT receptors. Using this mathematical model, we subsequently investigated the efficacy of several conjugates in cancer cells: DT and CRM107 conjugated to wild-type Tf, as well as to our engineered mutant Tf proteins (K206E/R632A Tf and K206E/R534A Tf). We also investigated the selectivity of mutant Tf-CRM107 against non-neoplastic cells. Through the use of our mathematical model, we predicted that (i) mutant Tf-CRM107 exhibits a greater cytotoxicity than wild-type Tf-CRM107 against cancerous cells, (ii) this improvement was more drastic with CRM107 conjugates than with DT conjugates, and (iii) mutant Tf-CRM107 conjugates were selective against non-neoplastic cells. These predictions were validated with in vitro cytotoxicity experiments, demonstrating that mutant Tf-CRM107 conjugates is indeed a more suitable therapeutic agent. Validation from in vitro experiments also confirmed that such whole-cell kinetic models can be useful in cancer therapeutic design. link: http://identifiers.org/pubmed/28065783

Messenger RNA (mRNA) translation is an essential step in eukaryotic gene expression that contributes to the regulation of this process. We describe a deterministic model based on ordinary differential equations that describe mRNA translation in Saccharomyces cerevisiae. This model, which was parameterized using published data, was developed to examine the kinetic behaviour of translation initiation factors in response to amino acid availability. The model predicts that the abundance of the eIF1-eIF3-eIF5 complex increases under amino acid starvation conditions, suggesting a possible auxiliary role for these factors in modulating translation initiation in addition to the known mechanisms involving eIF2. Our analyses of the robustness of the mRNA translation model suggest that individual cells within a randomly generated population are sensitive to external perturbations (such as changes in amino acid availability) through Gcn2 signalling. However, the model predicts that individual cells exhibit robustness against internal perturbations (such as changes in the abundance of translation initiation factors and kinetic parameters). Gcn2 appears to enhance this robustness within the system. These findings suggest a trade-off between the robustness and performance of this biological network. The model also predicts that individual cells exhibit considerable heterogeneity with respect to their absolute translation rates, due to random internal perturbations. Therefore, averaging the kinetic behaviour of cell populations probably obscures the dynamic robustness of individual cells. This highlights the importance of single-cell measurements for evaluating network properties. link: http://identifiers.org/pubmed/20306461

Blockade of immune checkpoints has recently been shown as a revolutionary strategy in the fight against cancers. Based on recent mouse experiments and clinical trials, large tumors can be completely suppressed with an additional blockade of immune checkpoints. We construct mathematical models capturing key interactions among malignant tumor cells, CD4+ T cells, anti-tumor cytokines, and immune checkpoint inhibitor of CTLA-4 to explore the importance of immune checkpoints on regression of tumor. Our study shows that blockade of immune checkpoints plays essential roles in immune responses. Continuous and one day pulse immune therapies by either T cells, anti-tumor cytokines, anti-CTLA-4 or a joint therapy are administered to exam the effectiveness of immune therapies. Our investigation indicates anti-tumor cytokine is potentially a key factor in determining the future of the malignant tumor. The malignant tumor can be suppressed thoroughly with reasonable dosages of anti-tumor cytokines if pre-radiation along with anti-CTLA-4 therapy are implemented. link: http://identifiers.org/doi/10.1016/j.amc.2019.06.037

Cellular homeostasis is regulated by signals through multiple molecular networks that include protein phosphorylation and metabolites. However, where and when the signal flows through a network and regulates homeostasis has not been explored. We have developed a reconstruction method for the signal flow based on time-course phosphoproteome and metabolome data, using multiple databases, and have applied it to acute action of insulin, an important hormone for metabolic homeostasis. An insulin signal flows through a network, through signaling pathways that involve 13 protein kinases, 26 phosphorylated metabolic enzymes, and 35 allosteric effectors, resulting in quantitative changes in 44 metabolites. Analysis of the network reveals that insulin induces phosphorylation and activation of liver-type phosphofructokinase 1, thereby controlling a key reaction in glycolysis. We thus provide a versatile method of reconstruction of signal flow through the network using phosphoproteome and metabolome data. link: http://identifiers.org/pubmed/25131207

Cellular homeostasis is regulated by signals through multiple molecular networks that include protein phosphorylation and metabolites. However, where and when the signal flows through a network and regulates homeostasis has not been explored. We have developed a reconstruction method for the signal flow based on time-course phosphoproteome and metabolome data, using multiple databases, and have applied it to acute action of insulin, an important hormone for metabolic homeostasis. An insulin signal flows through a network, through signaling pathways that involve 13 protein kinases, 26 phosphorylated metabolic enzymes, and 35 allosteric effectors, resulting in quantitative changes in 44 metabolites. Analysis of the network reveals that insulin induces phosphorylation and activation of liver-type phosphofructokinase 1, thereby controlling a key reaction in glycolysis. We thus provide a versatile method of reconstruction of signal flow through the network using phosphoproteome and metabolome data. link: http://identifiers.org/pubmed/25131207

The multifactorial nature of disease motivates the use of systems-level analyses to understand their pathology. We used a systems biology approach to study tau aggregation, one of the hallmark features of Alzheimer's disease. A mathematical model was constructed to capture the current state of knowledge concerning tau's behavior and interactions in cells. The model was implemented in silico in the form of ordinary differential equations. The identifiability of the model was assessed and parameters were estimated to generate two cellular states: a population of solutions that corresponds to normal tau homeostasis and a population of solutions that displays aggregation-prone behavior. The model of normal tau homeostasis was robust to perturbations, and disturbances in multiple processes were required to achieve an aggregation-prone state. The aggregation-prone state was ultrasensitive to perturbations in diverse subsets of networks. Tau aggregation requires that multiple cellular parameters are set coordinately to a set of values that drive pathological assembly of tau. This model provides a foundation on which to build and increase our understanding of the series of events that lead to tau aggregation and may ultimately be used to identify critical intervention points that can direct the cell away from tau aggregation to aid in the treatment of tau-mediated (or related) aggregation diseases including Alzheimer's. link: http://identifiers.org/pubmed/21085645

Genistein is an endocrine-active compound (EAC) found in soy products. It has been linked to beneficial effects such as mammary tumor growth suppression and adverse endocrine-related effects such as reduced birth weight in rats and humans. In its conjugated form, genistein is excreted in the bile, which is a significant factor in its pharmacokinetics. Experimental data suggest that genistein induces a concentration-dependent suppression of biliary excretion. In this article, we describe a physiologically based pharmacokinetic (PBPK) model that focuses on biliary excretion with the goal of accurately simulating the observed suppression. The mathematical model is a system of nonlinear differential equations with state-dependent delay to describe biliary excretion. The model was analyzed to examine local existence and uniqueness of a solution to the equations. Furthermore, unknown parameters were estimated, and the mathematical model was compared against published experimental data. link: http://identifiers.org/pubmed/17024552

We developed and analyzed the mathematical model of the intrinsic pathway based on the current biochemical data on the kinetics of blood coagulation individual stages. The model includes eight differential equations describing the spatio-temporal dynamics of activation of factors XI, IX, X, II, I, VIII, V, and protein C. The assembly of tenase and prothrombinase complexes is considered as a function of calcium concentration. The spatial dynamics of coagulation was analyzed for the one-dimensional case. We examined the formation of active factors, their spreading, and growth of the clot from the site of injury in the direction perpendicular to the vessel wall, into the blood thickness. We assumed that the site of injury (in the model one boundary of the space segment under examination) becomes a source of the continuous influx of factor XIa. In the first part, we described the model, selected the parameters, etc. In the second part, we compared the model with experimental data obtained in the homogeneous system and analyzed the spatial dynamics of the clot growth. link: http://identifiers.org/pubmed/8948047

Understanding the dynamics and variability of protein circuitry requires accurate measurements in living cells as well as theoretical models. To address this, we employed one of the best-studied protein circuits in human cells, the negative feedback loop between the tumor suppressor p53 and the oncogene Mdm2. We measured the dynamics of fluorescently tagged p53 and Mdm2 over several days in individual living cells. We found that isogenic cells in the same environment behaved in highly variable ways following DNA-damaging gamma irradiation: some cells showed undamped oscillations for at least 3 days (more than 10 peaks). The amplitude of the oscillations was much more variable than the period. Sister cells continued to oscillate in a correlated way after cell division, but lost correlation after about 11 h on average. Other cells showed low-frequency fluctuations that did not resemble oscillations. We also analyzed different families of mathematical models of the system, including a novel checkpoint mechanism. The models point to the possible source of the variability in the oscillations: low-frequency noise in protein production rates, rather than noise in other parameters such as degradation rates. This study provides a view of the extensive variability of the behavior of a protein circuit in living human cells, both from cell to cell and in the same cell over time. link: http://identifiers.org/pubmed/16773083

Understanding the dynamics and variability of protein circuitry requires accurate measurements in living cells as well as theoretical models. To address this, we employed one of the best-studied protein circuits in human cells, the negative feedback loop between the tumor suppressor p53 and the oncogene Mdm2. We measured the dynamics of fluorescently tagged p53 and Mdm2 over several days in individual living cells. We found that isogenic cells in the same environment behaved in highly variable ways following DNA-damaging gamma irradiation: some cells showed undamped oscillations for at least 3 days (more than 10 peaks). The amplitude of the oscillations was much more variable than the period. Sister cells continued to oscillate in a correlated way after cell division, but lost correlation after about 11 h on average. Other cells showed low-frequency fluctuations that did not resemble oscillations. We also analyzed different families of mathematical models of the system, including a novel checkpoint mechanism. The models point to the possible source of the variability in the oscillations: low-frequency noise in protein production rates, rather than noise in other parameters such as degradation rates. This study provides a view of the extensive variability of the behavior of a protein circuit in living human cells, both from cell to cell and in the same cell over time. link: http://identifiers.org/pubmed/16773083

Understanding the dynamics and variability of protein circuitry requires accurate measurements in living cells as well as theoretical models. To address this, we employed one of the best-studied protein circuits in human cells, the negative feedback loop between the tumor suppressor p53 and the oncogene Mdm2. We measured the dynamics of fluorescently tagged p53 and Mdm2 over several days in individual living cells. We found that isogenic cells in the same environment behaved in highly variable ways following DNA-damaging gamma irradiation: some cells showed undamped oscillations for at least 3 days (more than 10 peaks). The amplitude of the oscillations was much more variable than the period. Sister cells continued to oscillate in a correlated way after cell division, but lost correlation after about 11 h on average. Other cells showed low-frequency fluctuations that did not resemble oscillations. We also analyzed different families of mathematical models of the system, including a novel checkpoint mechanism. The models point to the possible source of the variability in the oscillations: low-frequency noise in protein production rates, rather than noise in other parameters such as degradation rates. This study provides a view of the extensive variability of the behavior of a protein circuit in living human cells, both from cell to cell and in the same cell over time. link: http://identifiers.org/pubmed/16773083

Understanding the dynamics and variability of protein circuitry requires accurate measurements in living cells as well as theoretical models. To address this, we employed one of the best-studied protein circuits in human cells, the negative feedback loop between the tumor suppressor p53 and the oncogene Mdm2. We measured the dynamics of fluorescently tagged p53 and Mdm2 over several days in individual living cells. We found that isogenic cells in the same environment behaved in highly variable ways following DNA-damaging gamma irradiation: some cells showed undamped oscillations for at least 3 days (more than 10 peaks). The amplitude of the oscillations was much more variable than the period. Sister cells continued to oscillate in a correlated way after cell division, but lost correlation after about 11 h on average. Other cells showed low-frequency fluctuations that did not resemble oscillations. We also analyzed different families of mathematical models of the system, including a novel checkpoint mechanism. The models point to the possible source of the variability in the oscillations: low-frequency noise in protein production rates, rather than noise in other parameters such as degradation rates. This study provides a view of the extensive variability of the behavior of a protein circuit in living human cells, both from cell to cell and in the same cell over time. link: http://identifiers.org/pubmed/16773083

Understanding the dynamics and variability of protein circuitry requires accurate measurements in living cells as well as theoretical models. To address this, we employed one of the best-studied protein circuits in human cells, the negative feedback loop between the tumor suppressor p53 and the oncogene Mdm2. We measured the dynamics of fluorescently tagged p53 and Mdm2 over several days in individual living cells. We found that isogenic cells in the same environment behaved in highly variable ways following DNA-damaging gamma irradiation: some cells showed undamped oscillations for at least 3 days (more than 10 peaks). The amplitude of the oscillations was much more variable than the period. Sister cells continued to oscillate in a correlated way after cell division, but lost correlation after about 11 h on average. Other cells showed low-frequency fluctuations that did not resemble oscillations. We also analyzed different families of mathematical models of the system, including a novel checkpoint mechanism. The models point to the possible source of the variability in the oscillations: low-frequency noise in protein production rates, rather than noise in other parameters such as degradation rates. This study provides a view of the extensive variability of the behavior of a protein circuit in living human cells, both from cell to cell and in the same cell over time. link: http://identifiers.org/pubmed/16773083

Understanding the dynamics and variability of protein circuitry requires accurate measurements in living cells as well as theoretical models. To address this, we employed one of the best-studied protein circuits in human cells, the negative feedback loop between the tumor suppressor p53 and the oncogene Mdm2. We measured the dynamics of fluorescently tagged p53 and Mdm2 over several days in individual living cells. We found that isogenic cells in the same environment behaved in highly variable ways following DNA-damaging gamma irradiation: some cells showed undamped oscillations for at least 3 days (more than 10 peaks). The amplitude of the oscillations was much more variable than the period. Sister cells continued to oscillate in a correlated way after cell division, but lost correlation after about 11 h on average. Other cells showed low-frequency fluctuations that did not resemble oscillations. We also analyzed different families of mathematical models of the system, including a novel checkpoint mechanism. The models point to the possible source of the variability in the oscillations: low-frequency noise in protein production rates, rather than noise in other parameters such as degradation rates. This study provides a view of the extensive variability of the behavior of a protein circuit in living human cells, both from cell to cell and in the same cell over time. link: http://identifiers.org/pubmed/16773083

In plants, as in animals, the core mechanism to retain rhythmic gene expression relies on the interaction of multiple feedback loops. In recent years, molecular genetic techniques have revealed a complex network of clock components in Arabidopsis. To gain insight into the dynamics of these interactions, new components need to be integrated into the mathematical model of the plant clock. Our approach accelerates the iterative process of model identification, to incorporate new components, and to systematically test different proposed structural hypotheses. Recent studies indicate that the pseudo-response regulators PRR7 and PRR9 play a key role in the core clock of Arabidopsis. We incorporate PRR7 and PRR9 into an existing model involving the transcription factors TIMING OF CAB (TOC1), LATE ELONGATED HYPOCOTYL (LHY) and CIRCADIAN CLOCK ASSOCIATED (CCA1). We propose candidate models based on experimental hypotheses and identify the computational models with the application of an optimization routine. Validation is accomplished through systematic analysis of various mutant phenotypes. We introduce and apply sensitivity analysis as a novel tool for analyzing and distinguishing the characteristics of proposed architectures, which also allows for further validation of the hypothesized structures. link: http://identifiers.org/pubmed/17102803

In plants, as in animals, the core mechanism to retain rhythmic gene expression relies on the interaction of multiple feedback loops. In recent years, molecular genetic techniques have revealed a complex network of clock components in Arabidopsis. To gain insight into the dynamics of these interactions, new components need to be integrated into the mathematical model of the plant clock. Our approach accelerates the iterative process of model identification, to incorporate new components, and to systematically test different proposed structural hypotheses. Recent studies indicate that the pseudo-response regulators PRR7 and PRR9 play a key role in the core clock of Arabidopsis. We incorporate PRR7 and PRR9 into an existing model involving the transcription factors TIMING OF CAB (TOC1), LATE ELONGATED HYPOCOTYL (LHY) and CIRCADIAN CLOCK ASSOCIATED (CCA1). We propose candidate models based on experimental hypotheses and identify the computational models with the application of an optimization routine. Validation is accomplished through systematic analysis of various mutant phenotypes. We introduce and apply sensitivity analysis as a novel tool for analyzing and distinguishing the characteristics of proposed architectures, which also allows for further validation of the hypothesized structures. link: http://identifiers.org/pubmed/17102803

In plants, as in animals, the core mechanism to retain rhythmic gene expression relies on the interaction of multiple feedback loops. In recent years, molecular genetic techniques have revealed a complex network of clock components in Arabidopsis. To gain insight into the dynamics of these interactions, new components need to be integrated into the mathematical model of the plant clock. Our approach accelerates the iterative process of model identification, to incorporate new components, and to systematically test different proposed structural hypotheses. Recent studies indicate that the pseudo-response regulators PRR7 and PRR9 play a key role in the core clock of Arabidopsis. We incorporate PRR7 and PRR9 into an existing model involving the transcription factors TIMING OF CAB (TOC1), LATE ELONGATED HYPOCOTYL (LHY) and CIRCADIAN CLOCK ASSOCIATED (CCA1). We propose candidate models based on experimental hypotheses and identify the computational models with the application of an optimization routine. Validation is accomplished through systematic analysis of various mutant phenotypes. We introduce and apply sensitivity analysis as a novel tool for analyzing and distinguishing the characteristics of proposed architectures, which also allows for further validation of the hypothesized structures. link: http://identifiers.org/pubmed/17102803

Mathematical models of the action potential in the periphery and center of the rabbit sinoatrial (SA) node have been developed on the basis of published experimental data. Simulated action potentials are consistent with those recorded experimentally: the model-generated peripheral action potential has a more negative takeoff potential, faster upstroke, more positive peak value, prominent phase 1 repolarization, greater amplitude, shorter duration, and more negative maximum diastolic potential than the model-generated central action potential. In addition, the model peripheral cell shows faster pacemaking. The models behave qualitatively the same as tissue from the periphery and center of the SA node in response to block of tetrodotoxin-sensitive Na(+) current, L- and T-type Ca(2+) currents, 4-aminopyridine-sensitive transient outward current, rapid and slow delayed rectifying K(+) currents, and hyperpolarization-activated current. A one-dimensional model of a string of SA node tissue, incorporating regional heterogeneity, coupled to a string of atrial tissue has been constructed to simulate the behavior of the intact SA node. In the one-dimensional model, the spontaneous action potential initiated in the center propagates to the periphery at approximately 0.06 m/s and then into the atrial muscle at 0.62 m/s. link: http://identifiers.org/pubmed/10899081

The transcription factor p53 plays a central role in maintaining genomic integrity. Recent experiments in MCF7 cells have shown that p53 protein level rises and falls in distinct pulses in response to DNA damage. The amplitudes of and intervals between pulses seem to be independent of the extent of damage, and some cells generate regular pulses of p53 over many days. Identifying the molecular mechanisms responsible for such interesting behavior is an important and challenging problem. This paper describes four dual-feedback mechanisms that combine both positive and negative feedback loops, which have been identified in the signaling network responsible for p53 regulation. Mathematical models of all four mechanisms are analyzed to determine if they are consistent with experimental observations and to characterize subtle differences among the possible mechanisms. In addition, a novel molecular mechanism is proposed whereby p53 pulses may induce, at first, cell cycle arrest and, if sustained, cell death. The proposal accounts for basic features of p53-mediated responses to DNA damage and suggests new experiments to probe the dynamics of p53 signaling. link: http://identifiers.org/pubmed/17245126

The transcription factor p53 plays a central role in maintaining genomic integrity. Recent experiments in MCF7 cells have shown that p53 protein level rises and falls in distinct pulses in response to DNA damage. The amplitudes of and intervals between pulses seem to be independent of the extent of damage, and some cells generate regular pulses of p53 over many days. Identifying the molecular mechanisms responsible for such interesting behavior is an important and challenging problem. This paper describes four dual-feedback mechanisms that combine both positive and negative feedback loops, which have been identified in the signaling network responsible for p53 regulation. Mathematical models of all four mechanisms are analyzed to determine if they are consistent with experimental observations and to characterize subtle differences among the possible mechanisms. In addition, a novel molecular mechanism is proposed whereby p53 pulses may induce, at first, cell cycle arrest and, if sustained, cell death. The proposal accounts for basic features of p53-mediated responses to DNA damage and suggests new experiments to probe the dynamics of p53 signaling. link: http://identifiers.org/pubmed/17245126

The transcription factor p53 plays a central role in maintaining genomic integrity. Recent experiments in MCF7 cells have shown that p53 protein level rises and falls in distinct pulses in response to DNA damage. The amplitudes of and intervals between pulses seem to be independent of the extent of damage, and some cells generate regular pulses of p53 over many days. Identifying the molecular mechanisms responsible for such interesting behavior is an important and challenging problem. This paper describes four dual-feedback mechanisms that combine both positive and negative feedback loops, which have been identified in the signaling network responsible for p53 regulation. Mathematical models of all four mechanisms are analyzed to determine if they are consistent with experimental observations and to characterize subtle differences among the possible mechanisms. In addition, a novel molecular mechanism is proposed whereby p53 pulses may induce, at first, cell cycle arrest and, if sustained, cell death. The proposal accounts for basic features of p53-mediated responses to DNA damage and suggests new experiments to probe the dynamics of p53 signaling. link: http://identifiers.org/pubmed/17245126

Metabolic pathways have traditionally been described in terms of biochemical reactions and metabolites. With the use of structural genomics and systems biology, we generated a three-dimensional reconstruction of the central metabolic network of the bacterium Thermotoga maritima. The network encompassed 478 proteins, of which 120 were determined by experiment and 358 were modeled. Structural analysis revealed that proteins forming the network are dominated by a small number (only 182) of basic shapes (folds) performing diverse but mostly related functions. Most of these folds are already present in the essential core (approximately 30%) of the network, and its expansion by nonessential proteins is achieved with relatively few additional folds. Thus, integration of structural data with networks analysis generates insight into the function, mechanism, and evolution of biological networks. link: http://identifiers.org/pubmed/19762644

Existing experimental data have shown hypoxia to be an important factor affecting the proliferation of mesenchymal stromal cells (MSCs), but the contrasting observations made at various hypoxic levels raise the questions of whether hypoxia accelerates proliferation, and how. On the other hand, in order to meet the increasing demand of MSCs, an optimised bioreactor control strategy is needed to enhance in vitro production.A comprehensive, single-cell mathematical model has been constructed in this work, which combines cellular oxygen sensing with hypoxia-mediated cell cycle progression to predict cell cycle commitment as a proxy to proliferation rate. With oxygen levels defined for in vitro cell culture, the model predicts enhanced proliferation under intermediate (2-8%) and mild (8-15%) hypoxia and cell quiescence under severe (< 2%) hypoxia. Global sensitivity analysis and quasi-Monte Carlo simulation revealed that within a certain range (+/- 100%), model parameters affect (with varying significance) the minimum commitment time, but the existence of a range of optimal oxygen tension could be preserved with the hypothesized effects of Hif2α and reactive oxygen species (ROS). It appears that Hif2α counteracts Hif1α and ROS-mediated protein deactivation under intermediate hypoxia and normoxia (20%), respectively, to regulate the response of cell cycle commitment to oxygen tension.Overall, this modelling study offered an integrative framework to capture several interacting mechanisms and allowed in silico analysis of their individual and collective roles in shaping the hypoxia-mediated commitment to cell cycle. The model offers a starting point to the establishment of a suitable mechanism that can satisfactorily explain the different existing experimental observations from different studies, and warrants future extension and dedicated experimental validation to eventually support bioreactor optimisation. link: http://identifiers.org/pubmed/29606139

The availability of genome sequences, annotations and knowledge of the biochemistry underlying metabolic transformations has led to the generation of metabolic network reconstructions for a wide range of organisms in bacteria, archaea, and eukaryotes. When modeled using mathematical representations, a reconstruction can simulate underlying genotype-phenotype relationships. Accordingly, genome-scale metabolic models (GEMs) can be used to predict the response of organisms to genetic and environmental variations. A bottom-up reconstruction procedure typically starts by generating a draft model from existing annotation data on a target organism. For model species, this part of the process can be straightforward, due to the abundant organism-specific biochemical data. However, the process becomes complicated for non-model less-annotated species. In this paper, we present a draft liver reconstruction, ReCodLiver0.9, of Atlantic cod (Gadus morhua), a non-model teleost fish, as a practicable guide for cases with comparably few resources. Although the reconstruction is considered a draft version, we show that it already has utility in elucidating metabolic response mechanisms to environmental toxicants by mapping gene expression data of exposure experiments to the resulting model. link: http://identifiers.org/doi/10.3389/fmolb.2020.591406

The availability of genome sequences, annotations and knowledge of the biochemistry underlying metabolic transformations has led to the generation of metabolic network reconstructions for a wide range of organisms in bacteria, archaea, and eukaryotes. When modeled using mathematical representations, a reconstruction can simulate underlying genotype-phenotype relationships. Accordingly, genome-scale metabolic models (GEMs) can be used to predict the response of organisms to genetic and environmental variations. A bottom-up reconstruction procedure typically starts by generating a draft model from existing annotation data on a target organism. For model species, this part of the process can be straightforward, due to the abundant organism-specific biochemical data. However, the process becomes complicated for non-model less-annotated species. In this paper, we present a draft liver reconstruction, ReCodLiver0.9, of Atlantic cod (Gadus morhua), a non-model teleost fish, as a practicable guide for cases with comparably few resources. Although the reconstruction is considered a draft version, we show that it already has utility in elucidating metabolic response mechanisms to environmental toxicants by mapping gene expression data of exposure experiments to the resulting model. link: http://identifiers.org/doi/10.3389/fmolb.2020.591406

The availability of genome sequences, annotations and knowledge of the biochemistry underlying metabolic transformations has led to the generation of metabolic network reconstructions for a wide range of organisms in bacteria, archaea, and eukaryotes. When modeled using mathematical representations, a reconstruction can simulate underlying genotype-phenotype relationships. Accordingly, genome-scale metabolic models (GEMs) can be used to predict the response of organisms to genetic and environmental variations. A bottom-up reconstruction procedure typically starts by generating a draft model from existing annotation data on a target organism. For model species, this part of the process can be straightforward, due to the abundant organism-specific biochemical data. However, the process becomes complicated for non-model less-annotated species. In this paper, we present a draft liver reconstruction, ReCodLiver0.9, of Atlantic cod (Gadus morhua), a non-model teleost fish, as a practicable guide for cases with comparably few resources. Although the reconstruction is considered a draft version, we show that it already has utility in elucidating metabolic response mechanisms to environmental toxicants by mapping gene expression data of exposure experiments to the resulting model. link: http://identifiers.org/doi/10.3389/fmolb.2020.591406

Background:The coronavirus disease 2019 (COVID-19) is rapidly spreading in China and more than 30 countries over last two months. COVID-19 has multiple characteristics distinct from other infectious diseases, including high infectivity during incubation, time delay between real dynamics and daily observed number of confirmed cases, and the intervention effects of implemented quarantine and control measures. Methods:We develop a Susceptible, Un-quanrantined infected, Quarantined infected, Confirmed infected (SUQC) model to characterize the dynamics of COVID-19 and explicitly parameterize the intervention effects of control measures, which is more suitable for analysis than other existing epidemic models. Results:The SUQC model is applied to the daily released data of the confirmed infections to analyze the outbreak of COVID-19 in Wuhan, Hubei (excluding Wuhan), China (excluding Hubei) and four first-tier cities of China. We found that, before January 30, 2020, all these regions except Beijing had a reproductive number R &gt; 1, and after January 30, all regions had a reproductive number R &lt; 1, indicating that the quarantine and control measures are effective in preventing the spread of COVID-19. The confirmation rate of Wuhan estimated by our model is 0.0643, substantially lower than that of Hubei excluding Wuhan (0.1914), and that of China excluding Hubei (0.2189), but it jumps to 0.3229 after February 12 when clinical evidence was adopted in new diagnosis guidelines. The number of unquarantined infected cases in Wuhan on February 12, 2020 is estimated to be 3,509 and declines to 334 on February 21, 2020. After fitting the model with data as of February 21, 2020, we predict that the end time of COVID-19 in Wuhan and Hubei is around late March, around mid March for China excluding Hubei, and before early March 2020 for the four tier-one cities. A total of 80,511 individuals are estimated to be infected in China, among which 49,510 are from Wuhan, 17,679 from Hubei (excluding Wuhan), and the rest 13,322 from other regions of China (excluding Hubei). Note that the estimates are from a deterministic ODE model and should be interpreted with some uncertainty. Conclusions:We suggest that rigorous quarantine and control measures should be kept before early March in Beijing, Shanghai, Guangzhou and Shenzhen, and before late March in Hubei. The model can also be useful to predict the trend of epidemic and provide quantitative guide for other countries at high risk of outbreak, such as South Korea, Japan, Italy and Iran. Supplementary Materials:The supplementary materials can be found online with this article at 10.1007/s40484-020-0199-0. link: http://identifiers.org/pubmed/32219006

Factor Xa (FXa) emerged as a promising target for effective anticoagulation and several FXa inhibitors are now available for the prevention of venous thromboembolism. However, in previously reported pharmacokinetic/pharmacodynamic (PK/PD) models, the complex coagulation processes and detailed information of drug action are usually unclear, which makes it difficult to predict clinical outcome at the drug discovery stage. In this study, a large-scale systems pharmacology model was developed based on several published models and clinical data. It takes into account all pathways of the coagulation network, and captures drug-specific features: plasma pharmacokinetics and drug-target binding kinetics (BKs). We aimed to predict the anticoagulation effects of FXa inhibitors in healthy subjects, and to use this model to compare the effects of compounds with different binding properties. Our model predicts the clotting time and anti-FXa effects and could thus serve as a predictive tool for the anticoagulant potential of a new compound. link: http://identifiers.org/pubmed/26783501

Recent studies have shown that in addition to the transcriptional circadian clock, many organisms, including Arabidopsis, have a circadian redox rhythm driven by the organism's metabolic activities. It has been hypothesized that the redox rhythm is linked to the circadian clock, but the mechanism and the biological significance of this link have only begun to be investigated. Here we report that the master immune regulator NPR1 (non-expressor of pathogenesis-related gene 1) of Arabidopsis is a sensor of the plant's redox state and regulates transcription of core circadian clock genes even in the absence of pathogen challenge. Surprisingly, acute perturbation in the redox status triggered by the immune signal salicylic acid does not compromise the circadian clock but rather leads to its reinforcement. Mathematical modelling and subsequent experiments show that NPR1 reinforces the circadian clock without changing the period by regulating both the morning and the evening clock genes. This balanced network architecture helps plants gate their immune responses towards the morning and minimize costs on growth at night. Our study demonstrates how a sensitive redox rhythm interacts with a robust circadian clock to ensure proper responsiveness to environmental stimuli without compromising fitness of the organism. link: http://identifiers.org/pubmed/26098366

This work presents both deterministic and stochastic models of genetic expression modulated by intracellular calcium (Ca2+) oscillations, based on macroscopic differential equations and chemical Langevin equations, respectively. In deterministic case, the oscillations of intracellular Ca2+ decrease the effective Ca2+ threshold for the activation of transcriptional activator (TF-A). The average activation of TF-A increases with the increase of the average amplitude of intracellular Ca2+ oscillations, but decreases with the increase of the period of intracellular Ca2+ oscillations, which are qualitatively consistent with the experimental results on the gene expression in lymphocytes. In stochastic case, it is found that a large internal fluctuation of the biochemical reaction can enhance gene expression efficiency specifically at a low level of external stimulations or at a small rate of TF-A dimer phosphorylation activated by Ca2+, which reduces the threshold of the average intracellular Ca2+ concentration for gene expression. link: http://identifiers.org/pubmed/17560007

Combination chemotherapy is standard treatment for pancreatic cancer. However, current drugs lack efficacy for most patients, and selection and evaluation of new combination regimens is empirical and time-consuming. The efficacy of gemcitabine, a standard-of-care agent, combined with birinapant, a pro-apoptotic antagonist of Inhibitor of Apoptosis Proteins (IAPs), was investigated in pancreatic cancer cells. PANC-1 cells were treated with vehicle, gemcitabine (6, 10, 20 nM), birinapant (50, 200, 500 nM), and combinations of the two drugs. Temporal changes in cell numbers, cell cycle distribution, and apoptosis were measured. A basic pharmacodynamic (PD) model based on cell numbers, and a mechanism-based PD model integrating all measurements, were developed. The basic PD model indicated that synergistic effects occurred in both cell proliferation and death processes. The mechanism-based model captured key features of drug action: temporary cell cycle arrest in S phase induced by gemcitabine alone, apoptosis induced by birinapant alone, and prolonged cell cycle arrest and enhanced apoptosis induced by the combination. A drug interaction term Ψ was employed in the models to signify interactions of the combination when data were limited. When more experimental information was utilized, Ψ values approaching 1 indicated that specific mechanisms of interactions were captured better. PD modeling identified the potential benefit of combining gemcitabine and birinapant, and characterized the key interaction pathways. An optimal treatment schedule of pretreatment with gemcitabine for 24-48 h was suggested based on model predictions and was verified experimentally. This approach provides a generalizable modeling platform for exploring combinations of cytostatic and cytotoxic agents in cancer cell culture studies. link: http://identifiers.org/pubmed/26252969

Combination chemotherapy is standard treatment for pancreatic cancer. However, current drugs lack efficacy for most patients, and selection and evaluation of new combination regimens is empirical and time-consuming. The efficacy of gemcitabine, a standard-of-care agent, combined with birinapant, a pro-apoptotic antagonist of Inhibitor of Apoptosis Proteins (IAPs), was investigated in pancreatic cancer cells. PANC-1 cells were treated with vehicle, gemcitabine (6, 10, 20 nM), birinapant (50, 200, 500 nM), and combinations of the two drugs. Temporal changes in cell numbers, cell cycle distribution, and apoptosis were measured. A basic pharmacodynamic (PD) model based on cell numbers, and a mechanism-based PD model integrating all measurements, were developed. The basic PD model indicated that synergistic effects occurred in both cell proliferation and death processes. The mechanism-based model captured key features of drug action: temporary cell cycle arrest in S phase induced by gemcitabine alone, apoptosis induced by birinapant alone, and prolonged cell cycle arrest and enhanced apoptosis induced by the combination. A drug interaction term Ψ was employed in the models to signify interactions of the combination when data were limited. When more experimental information was utilized, Ψ values approaching 1 indicated that specific mechanisms of interactions were captured better. PD modeling identified the potential benefit of combining gemcitabine and birinapant, and characterized the key interaction pathways. An optimal treatment schedule of pretreatment with gemcitabine for 24-48 h was suggested based on model predictions and was verified experimentally. This approach provides a generalizable modeling platform for exploring combinations of cytostatic and cytotoxic agents in cancer cell culture studies. link: http://identifiers.org/pubmed/26252969