Modeling the activation of the alternative complement pathway and its effects on hemolysis in health and disease.

The complement system is a powerful mechanism of innate immunity poised to eliminate foreign cells and pathogens. It is an intricate network of >35 proteins, which, once activated, leads to the tagging of the surface to be eliminated, produces potent chemoattractants to recruit immune cells, and inserts cytotoxic pores into nearby lipid surfaces. Although it can be triggered via different pathways, its net output is largely based on the direct or indirect activation of the alternative pathway. Complement dysregulation or deficiencies may cause severe pathologies, such as paroxysmal nocturnal hemoglobinuria (PNH), where a lack of complement control proteins leads to hemolysis and life-threatening anemia. The complexity of the system poses a challenge for the interpretation of experimental data and the design of effective pharmacological therapies. To address this issue, we developed a mathematical model of the alternative complement pathway building on previous modelling efforts. The model links complement activation to the hemolytic activity of the terminal alternative pathway, providing an accurate description of pathway activity as observed in vitro and in vivo, in health and disease. Through adjustment of the parameters describing experimental conditions, the model was capable of reproducing the results of an array of standard assays used in complement research. To demonstrate its clinical applicability, we compared model predictions with clinical observations of the recovery of hematological biomarkers in PNH patients treated with the complement inhibiting anti-C5 antibody eculizumab. In conclusion, the model can enhance the understanding of complement biology and its role in disease pathogenesis, help identifying promising targets for pharmacological intervention, and predict the outcome of complement-targeting pharmacological interventions.

A Shh/Gli-driven three-node timer motif controls temporal identity and fate of neural stem cells.

How time is measured by neural stem cells during temporal neurogenesis has remained unresolved. By combining experiments and computational modeling, we define a Shh/Gli-driven three-node timer underlying the sequential generation of motor neurons (MNs) and serotonergic neurons in the brainstem. The timer is founded on temporal decline of Gli-activator and Gli-repressor activities established through down-regulation of Gli transcription. The circuitry conforms an incoherent feed-forward loop, whereby Gli proteins not only promote expression of Phox2b and thereby MN-fate but also account for a delayed activation of a self-promoting transforming growth factor-ß (Tgfß) node triggering a fate switch by repressing Phox2b. Hysteresis and spatial averaging by diffusion of Tgfß counteract noise and increase temporal accuracy at the population level, providing a functional rationale for the intrinsically programmed activation of extrinsic switch signals in temporal patterning. Our study defines how time is reliably encoded during the sequential specification of neurons.

Growth and size control during development.

The size and shape of organs are characteristic for each species. Even when organisms develop to different sizes due to varying environmental conditions, such as nutrition, organ size follows species-specific rules of proportionality to the rest of the body, a phenomenon referred to as allometry. Therefore, for a given environment, organs stop growth at a predictable size set by the species's genotype. How do organs stop growth? How can related species give rise to organs of strikingly different size? No definitive answer has been given to date. One of the major models for the studies of growth termination is the vinegar fly Drosophila melanogaster. Therefore, this review will focus mostly on work carried out in Drosophila to try to tease apart potential mechanisms and identify routes for further investigation. One general rule, found across the animal kingdom, is that the rate of growth declines with developmental time. Therefore, answers to the problem of growth termination should explain this seemingly universal fact. In addition, growth termination is intimately related to the problems of robustness (i.e. precision) and plasticity in organ size, symmetric and asymmetric organ development, and of how the 'target' size depends on extrinsic, environmental factors.

Growth control in the Drosophila eye disc by the cytokine Unpaired.

A fundamental question in developmental biology is how organ size is controlled. We have previously shown that the area growth rate in the Drosophila eye primordium declines inversely proportionally to the increase in its area. How the observed reduction in the growth rate is achieved is unknown. Here, we explore the dilution of the cytokine Unpaired (Upd) as a possible candidate mechanism. In the developing eye, upd expression is transient, ceasing at the time when the morphogenetic furrow first emerges. We confirm experimentally that the diffusion and stability of the JAK/STAT ligand Upd are sufficient to control eye disc growth via a dilution mechanism. We further show that sequestration of Upd by ectopic expression of an inactive form of the receptor Domeless (Dome) results in a substantially lower growth rate, but the area growth rate still declines inversely proportionally to the area increase. This growth rate-to-area relationship is no longer observed when Upd dilution is prevented by the continuous, ectopic expression of Upd. We conclude that a mechanism based on the dilution of the growth modulator Upd can explain how growth termination is controlled in the eye disc.

An Unbiased Analysis of Candidate Mechanisms for the Regulation of Drosophila Wing Disc Growth.

The control of organ size presents a fundamental open problem in biology. A declining growth rate is observed in all studied higher animals, and the growth limiting mechanism may therefore be evolutionary conserved. Most studies of organ growth control have been carried out in Drosophila imaginal discs. We have previously shown that the area growth rate in the Drosophila eye primordium declines inversely proportional to the increase in its area, which is consistent with a dilution mechanism for growth control. Here, we show that a dilution mechanism cannot explain growth control in the Drosophila wing disc. We computationally evaluate a range of alternative candidate mechanisms and show that the experimental data can be best explained by a biphasic growth law. However, also logistic growth and an exponentially declining growth rate fit the data very well. The three growth laws correspond to fundamentally different growth mechanisms that we discuss. Since, as we show, a fit to the available experimental growth kinetics is insufficient to define the underlying mechanism of growth control, future experimental studies must focus on the molecular mechanisms to define the mechanism of growth control.

A quantitative analysis of growth control in the Drosophila eye disc.

The size and shape of organs is species specific, and even in species in which organ size is strongly influenced by environmental cues, such as nutrition or temperature, it follows defined rules. Therefore, mechanisms must exist to ensure a tight control of organ size within a given species, while being flexible enough to allow for the evolution of different organ sizes in different species. We combined computational modeling and quantitative measurements to analyze growth control in the Drosophila eye disc. We find that the area growth rate declines inversely proportional to the increasing total eye disc area. We identify two growth laws that are consistent with the growth data and that would explain the extraordinary robustness and evolutionary plasticity of the growth process and thus of the final adult eye size. These two growth laws correspond to very different control mechanisms and we discuss how each of these laws constrains the set of candidate biological mechanisms for growth control in the Drosophila eye disc.

Deciphering the mechanism behind Fibroblast Growth Factor (FGF) induced biphasic signal-response profiles.

BACKGROUND: The Fibroblast Growth Factor (FGF) pathway is driving various aspects of cellular responses in both normal and malignant cells. One interesting characteristic of this pathway is the biphasic nature of the cellular response to some FGF ligands like FGF2. Specifically, it has been shown that phenotypic behaviors controlled by FGF signaling, like migration and growth, reach maximal levels in response to intermediate concentrations, while high levels of FGF2 elicit weak responses. The mechanisms leading to the observed biphasic response remains unexplained. RESULTS: A combination of experiments and computational modeling was used to understand the mechanism behind the observed biphasic signaling responses. FGF signaling involves a tertiary surface interaction that we captured with a computational model based on Ordinary Differential Equations (ODEs). It accounts for FGF2 binding to FGF receptors (FGFRs) and heparan sulfate glycosaminoglycans (HSGAGs), followed by receptor-phosphorylation, activation of the FRS2 adapter protein and the Ras-Raf signaling cascade. Quantitative protein assays were used to measure the dynamics of phosphorylated ERK (pERK) in response to a wide range of FGF2 ligand concentrations on a fine-grained time scale for the squamous cell lung cancer cell line H1703. We developed a novel approach combining Particle Swarm Optimization (PSO) and feature-based constraints in the objective function to calibrate the computational model to the experimental data. The model is validated using a series of extracellular and intracellular perturbation experiments. We demonstrate that in silico model predictions are in accordance with the observed in vitro results. CONCLUSIONS: Using a combined approach of computational modeling and experiments we found that competition between binding of the ligand FGF2 to HSGAG and FGF receptor leads to the biphasic response. At low to intermediate concentrations of FGF2 there are sufficient free FGF receptors available for the FGF2-HSGAG complex to enable the formation of the trimeric signaling unit. At high ligand concentrations the ligand binding sites of the receptor become saturated and the trimeric signaling unit cannot be formed. This insight into the pathway is an important consideration for the pharmacological inhibition of this pathway.