Operon dynamics with state dependent transcription and/or translation delays.

Transcription and translation retrieve and operationalize gene encoded information in cells. These processes are not instantaneous and incur significant delays. In this paper we study Goodwin models of both inducible and repressible operons with state-dependent delays. The paper provides justification and derivation of the model, detailed analysis of the appropriate setting of the corresponding dynamical system, and extensive numerical analysis of its dynamics. Comparison with constant delay models shows significant differences in dynamics that include existence of stable periodic orbits in inducible systems and multistability in repressible systems. A combination of parameter space exploration, numerics, analysis of steady state linearization and bifurcation theory indicates the likely presence of Shilnikov-type homoclinic bifurcations in the repressible operon model.

How can we describe density evolution under delayed dynamics?

Although the theory of density evolution in maps and ordinary differential equations is well developed, the situation is far from satisfactory in continuous time systems with delay. This paper reviews some of the work that has been done numerically, the interesting dynamics that have emerged, and the largely unsuccessful attempts that have been made to analytically treat the evolution of densities in differential delay equations. We also present a new approach to the problem and illustrate it with a simple example.

The timing of cyclic cytotoxic chemotherapy can worsen neutropenia and neutrophilia.

Despite recent advances in immunotherapies, cytotoxic chemotherapy continues to be a first-line treatment option for the majority of cancers. Unfortunately, a common side effect in patients undergoing chemotherapy treatment is neutropenia. To mitigate the risk of neutropenia and febrile neutropenia, prophylactic treatment with granulocyte-colony stimulating factor (G-CSF) is administered. Extensive pharmacokinetic/pharmacodynamic modelling of myelosuppression during chemotherapy has suggested avenues for therapy optimization to mitigate this neutropenia. However, the issue of resonance, whereby neutrophil oscillations are induced by the periodic administration of cytotoxic chemotherapy and the coadministration of G-CSF, potentially aggravating a patient's neutropenic/neutrophilic status, is not well-characterized in the clinical literature. Here, through analysis of neutrophil data from young acute lymphoblastic leukaemia patients, we find that resonance is occurring during cyclic chemotherapy treatment in 26% of these patients. Motivated by these data and our previous modelling studies on adult lymphoma patients, we examined resonance during treatment with or without G-CSF. Using our quantitative systems pharmacology model of granulopoiesis, we show that the timing of cyclic chemotherapy can worsen neutropenia or neutrophilia, and suggest clinically-actionable schedules to reduce the resonant effect. We emphasize that delaying supportive G-CSF therapy to 6-7 days after chemotherapy can mitigate myelosuppressive effects. This study therefore highlights the importance of quantitative systems pharmacology for the clinical practice for developing rational therapeutic strategies.

Characterizing Chemotherapy-Induced Neutropenia and Monocytopenia Through Mathematical Modelling.

In spite of the recent focus on the development of novel targeted drugs to treat cancer, cytotoxic chemotherapy remains the standard treatment for the vast majority of patients. Unfortunately, chemotherapy is associated with high hematopoietic toxicity that may limit its efficacy. We have previously established potential strategies to mitigate chemotherapy-induced neutropenia (a lack of circulating neutrophils) using a mechanistic model of granulopoiesis to predict the interactions defining the neutrophil response to chemotherapy and to define optimal strategies for concurrent chemotherapy/prophylactic granulocyte colony-stimulating factor (G-CSF). Here, we extend our analyses to include monocyte production by constructing and parameterizing a model of monocytopoiesis. Using data for neutrophil and monocyte concentrations during chemotherapy in a large cohort of childhood acute lymphoblastic leukemia patients, we leveraged our model to determine the relationship between the monocyte and neutrophil nadirs during cyclic chemotherapy. We show that monocytopenia precedes neutropenia by 3 days, and rationalize the use of G-CSF during chemotherapy by establishing that the onset of monocytopenia can be used as a clinical marker for G-CSF dosing post-chemotherapy. This work therefore has important clinical applications as a comprehensive approach to understanding the relationship between monocyte and neutrophils after cyclic chemotherapy with or without G-CSF support.

Periodic hematological disorders: Quintessential examples of dynamical diseases.

This paper summarizes the evidence supporting the classification of cyclic neutropenia as a dynamical disease and periodic chronic myelogenous leukemia is also considered. The unsatisfactory state of knowledge concerning the genesis of cyclic thrombocytopenia and periodic autoimmune hemolytic anemia is detailed.

Response of an oscillatory differential delay equation to a periodic stimulus.

Periodic hematological diseases such as cyclical neutropenia or cyclical thrombocytopenia, with their characteristic oscillations of circulating neutrophils or platelets, may pose grave problems for patients. Likewise, periodically administered chemotherapy has the unintended side effect of establishing periodic fluctuations in circulating white cells, red cell precursors and/or platelets. These fluctuations, either spontaneous or induced, often have serious consequences for the patient (e.g. neutropenia, anemia, or thrombocytopenia respectively) which exogenously administered cytokines can partially correct. The question of when and how to administer these drugs is a difficult one for clinicians and not easily answered. In this paper we use a simple model consisting of a delay differential equation with a piecewise linear nonlinearity, that has a periodic solution, to model the effect of a periodic disease or periodic chemotherapy. We then examine the response of this toy model to both single and periodic perturbations, meant to mimic the drug administration, as a function of the drug dose and the duration and frequency of its administration to best determine how to avoid side effects.

Neutropenia in Barth syndrome: characteristics, risks, and management.

PURPOSE OF REVIEW: Barth syndrome (BTHS) is an X-linked disease characterized by defective remodeling of phospholipid side chains in mitochondrial membranes. Major features include neutropenia, dilated cardiomyopathy, motor delay and proximal myopathy, feeding problems, and constitutional growth delay. We conducted this review of neutropenia in BTHS to aid in the diagnosis of this disease, and to improve understanding of both the consequences of neutropenia and the benefits of treatment with granulocyte colony-stimulating factor (G-CSF). RECENT FINDINGS: In 88 patients with BTHS, neutropenia, that is, at least one count below 1.5 × 10/l, was detected in 74 (84%) and 44% had severe chronic neutropenia, with multiple counts below 0.5 × 10/l. The pattern of neutropenia varied between intermittent and unpredictable, chronic and severe, or cyclical with mathematically regular oscillations. Monocytosis, that is, monocytes more than 1.0 × 10/l, was observed at least once in 64 of 85 (75%) patients. G-CSF was administered to 39 of 88 patients (44%). Weekly average G-CSF doses ranged from 0.12 to 10.92 µg/kg/day (mean 1.16 µg/kg/day, median 1.16 µg/kg/day). Antibiotic prophylaxis was additionally employed in 21 of 26 neutropenic patients. Pretreatment bone marrow evaluations predominantly showed reduced myeloid maturation which normalized on G-CSF therapy in seven of 13 examined. Consistent clinical improvement, with reduced signs and symptoms of infections, was observed in response to prophylactic G-CSF ±â€Šprophylactic antibiotics. However, despite G-CSF and antibiotics, one adult patient died with multiple infections related to indwelling medical devices and gastrostomy site infection after 15.5 years on G-CSF and a pediatric patient required gastrostomy removal for recurrent abdominal wall cellulitis. SUMMARY: BTHS should be considered in any men with neutropenia accompanied by any of the characteristic features of this syndrome. Prophylaxis with G-CSF ±â€Šantibiotics prevents serious bacterial infections in the more severe neutropenic patients although infections remain a threat even in patients who are very compliant with therapy, especially in those with indwelling devices.

Origins of oscillation patterns in cyclical thrombocytopenia.

Cyclical thrombocytopenia (CT) is a rare hematological disease characterized by periodic oscillations in circulating platelet counts. In almost all CT patients, other cell lines show no sign of oscillation, but recently a CT patient was reported with significant oscillations in circulating neutrophils (in the same period as the platelets). In this paper, we attempt to understand this phenomenon through a previously published model of human hematopoiesis. We have investigated a variety of possible oscillation patterns that may appear when alterations occur in the control parameters in the platelet regulatory dynamics. Our results indicate that the platelet maturation time and the differentiation rate from hematopoietic stem cells (HSCs) into the platelet cell line play important roles in the emergence of various types of CT like oscillations. Moreover, we find different oscillation patterns, including CT and cyclical neutropenia like oscillations, with certain parameter values in the platelet compartment. A bifurcation analysis revealed the different origins of these oscillation patterns. We also identified bistable dynamics which indicate the potential importance of system history in the treatment of these diseases. Together, these results demonstrate the possible origins for various oscillation patterns dependent on alterations in the platelet cell line control mechanisms. One of the important origins may be related to the regulation of apoptosis in platelet precursors.

Cyclic thrombocytopenia with statistically significant neutrophil oscillations.

Cyclic thrombocytopenia is often misdiagnosed as immune thrombocytopenia due to similar clinical features, a fact of significance because cyclic thrombocytopenia generally responds poorly to treatments used successfully in immune thrombocytopenia. A precise diagnosis must establish the statistical significance of periodicity of the platelet counts using statistical methods (eg, Lomb-Scargle periodogram).

Normal and pathological dynamics of platelets in humans.

We develop a mathematical model of platelet, megakaryocyte, and thrombopoietin dynamics in humans. We show that there is a single stationary solution that can undergo a Hopf bifurcation, and use this information to investigate both normal and pathological platelet production, specifically cyclic thrombocytopenia. Carefully estimating model parameters from laboratory and clinical data, we then argue that a subset of parameters are involved in the genesis of cyclic thrombocytopenia based on clinical information. We provide model fits to the existing data for both platelet counts and thrombopoietin levels by changing four parameters that have physiological correlates. Our results indicate that the primary change in cyclic thrombocytopenia is an interference with, or destruction of, the thrombopoietin receptor with secondary changes in other processes, including immune-mediated destruction of platelets and megakaryocyte deficiency and failure in platelet production. This study contributes to the understanding of the origin of cyclic thrombocytopenia as well as extending the modeling of thrombopoiesis.

Response of an oscillatory differential delay equation to a single stimulus.

Here we analytically examine the response of a limit cycle solution to a simple differential delay equation to a single pulse perturbation of the piecewise linear nonlinearity. We construct the unperturbed limit cycle analytically, and are able to completely characterize the perturbed response to a pulse of positive amplitude and duration with onset at different points in the limit cycle. We determine the perturbed minima and maxima and period of the limit cycle and show how the pulse modifies these from the unperturbed case.

Mathematical model of galactose regulation and metabolic consumption in yeast.

The galactose network has been extensively studied at the unicellular level to broaden our understanding of the regulatory mechanisms governing galactose metabolism in multicellular organisms. Although the key molecular players involved in the metabolic and regulatory processes of this system have been known for decades, their interactions and chemical kinetics remain incompletely understood. Mathematical models can provide an alternative method to study the dynamics of this network from a quantitative and a qualitative perspective. Here, we employ this approach to unravel the main properties of the galactose network, including equilibrium binary and temporal responses, as a way to decipher its adaptation to actively-changing inputs. We combine its two main components: the genetic branch, which allows for bistable responses, and a metabolic branch, encompassing the relevant metabolic processes that can be repressed by glucose. We use both computational tools to estimate model parameters based on published experimental data, as well as bifurcation analysis to decipher the properties of the system in various parameter regimes. Our model analysis reveals that the interplay between the inducer (galactose) and the repressor (glucose) creates a bistable regime which dictates the temporal responses of the system. Based on the same bifurcation techniques, we explain why the system is robust to genetic mutations and molecular instabilities. These findings may provide experimentalists with a theoretical framework with which they can determine how the galactose network functions under various conditions.

The limiting dynamics of a bistable molecular switch with and without noise.

We consider the dynamics of a population of organisms containing two mutually inhibitory gene regulatory networks, that can result in a bistable switch-like behaviour. We completely characterize their local and global dynamics in the absence of any noise, and then go on to consider the effects of either noise coming from bursting (transcription or translation), or Gaussian noise in molecular degradation rates when there is a dominant slow variable in the system. We show analytically how the steady state distribution in the population can range from a single unimodal distribution through a bimodal distribution and give the explicit analytic form for the invariant stationary density which is globally asymptotically stable. Rather remarkably, the behaviour of the stationary density with respect to the parameters characterizing the molecular behaviour of the bistable switch is qualitatively identical in the presence of noise coming from bursting as well as in the presence of Gaussian noise in the degradation rate. This implies that one cannot distinguish between either the dominant source or nature of noise based on the stationary molecular distribution in a population of cells. We finally show that the switch model with bursting but two dominant slow genes has an asymptotically stable stationary density.

Neutrophil dynamics during concurrent chemotherapy and G-CSF administration: Mathematical modelling guides dose optimisation to minimise neutropenia.

The choice of chemotherapy regimens is often constrained by the patient's tolerance to the side effects of chemotherapeutic agents. This dose-limiting issue is a major concern in dose regimen design, which is typically focused on maximising drug benefits. Chemotherapy-induced neutropenia is one of the most prevalent toxic effects patients experience and frequently threatens the efficient use of chemotherapy. In response, granulocyte colony-stimulating factor (G-CSF) is co-administered during chemotherapy to stimulate neutrophil production, increase neutrophil counts, and hopefully avoid neutropenia. Its clinical use is, however, largely dictated by trial and error processes. Based on up-to-date knowledge and rational considerations, we develop a physiologically realistic model to mathematically characterise the neutrophil production in the bone marrow which we then integrate with pharmacokinetic and pharmacodynamic (PKPD) models of a chemotherapeutic agent and an exogenous form of G-CSF (recombinant human G-CSF, or rhG-CSF). In this work, model parameters represent the average values for a general patient and are extracted from the literature or estimated from available data. The dose effect predicted by the model is confirmed through previously published data. Using our model, we were able to determine clinically relevant dosing regimens that advantageously reduce the number of rhG-CSF administrations compared to original studies while significantly improving the neutropenia status. More particularly, we determine that it could be beneficial to delay the first administration of rhG-CSF to day seven post-chemotherapy and reduce the number of administrations from ten to three or four for a patient undergoing 14-day periodic chemotherapy.

Phosphate and ADP differently inhibit coordinated smooth muscle myosin groups.

Actin filaments propelled in vitro by groups of skeletal muscle myosin motors exhibit distinct phases of active sliding or arrest, whose occurrence depends on actin length (L) within a range of up to 1.0 µm. Smooth muscle myosin filaments are exponentially distributed with ≈150 nm average length in vivo--suggesting relevance of the L-dependence of myosin group kinetics. Here, we found L-dependent actin arrest and sliding in in vitro motility assays of smooth muscle myosin. We perturbed individual myosin kinetics with varying, physiological concentrations of phosphate (Pi, release associated with main power stroke) and adenosine diphosphate (ADP, release associated with minor mechanical step). Adenosine triphosphate was kept constant at physiological concentration. Increasing [Pi] lowered the fraction of time for which actin was actively sliding, reflected in reduced average sliding velocity (ν) and motile fraction (fmot, fraction of time that filaments are moving); increasing [ADP] increased the fraction of time actively sliding and reduced the velocity while sliding, reflected in reduced ν and increased fmot. We introduced specific Pi and ADP effects on individual myosin kinetics into our recently developed mathematical model of actin propulsion by myosin groups. Simulations matched our experimental observations and described the inhibition of myosin group kinetics. At low [Pi] and [ADP], actin arrest and sliding were reflected by two distinct chemical states of the myosin group. Upon [Pi] increase, the probability of the active state decreased; upon [ADP] increase, the probability of the active state increased, but the active state became increasingly similar to the arrested state.

Understanding, treating and avoiding hematological disease: better medicine through mathematics?

This paper traces the experimental, clinical and mathematical modeling efforts to understand a periodic hematological disease-cyclical neutropenia. It is primarily a highly personal account by two scientists from quite different backgrounds of their interactions over almost 40 years and their attempts to understand this intriguing disease. It's also a story of their efforts to offer effective treatments for the patients who suffer from cyclic neutropenia and other conditions causing neutropenia and infections.

The utility of simple mathematical models in understanding gene regulatory dynamics.

In this review, we survey work that has been carried out in the attempts of biomathematicians to understand the dynamic behaviour of simple bacterial operons starting with the initial work of the 1960's. We concentrate on the simplest of situations, discussing both repressible and inducible systems and then turning to concrete examples related to the biology of the lactose and tryptophan operons. We conclude with a brief discussion of the role of both extrinsic noise and so-called intrinsic noise in the form of translational and/or transcriptional bursting.

Understanding and treating cytopenia through mathematical modeling.

Here, we briefly review how the study of dynamic hematological diseases with mathematical modeling tools has led to a better understanding of the origin of some types of neutropenia and thrombocytopenia and to improved treatment strategies. In addition, we have briefly discussed how these models suggest improved ways to minimize and/or treat cytopenia induced by chemotherapy.