Optimal and continuous anaemia control in a cohort of dialysis patients in Switzerland.

BACKGROUND: Guidelines for the management of anaemia in patients with chronic kidney disease (CKD) recommend a minimal haemoglobin (Hb) target of 11 g/dL. Recent surveys indicate that this requirement is not met in many patients in Europe. In most studies, Hb is only assessed over a short-term period. The aim of this study was to examine the control of anaemia over a continuous long-term period in Switzerland. METHODS: A prospective multi-centre observational study was conducted in dialysed patients treated with recombinant human epoetin (EPO) beta, over a one-year follow-up period, with monthly assessments of anaemia parameters. RESULTS: Three hundred and fifty patients from 27 centres, representing 14% of the dialysis population in Switzerland, were included. Mean Hb was 11.9 +/- 1.0 g/dL, and remained stable over time. Eighty-five % of the patients achieved mean Hb >or= 11 g/dL. Mean EPO dose was 155 +/- 118 IU/kg/week, being delivered mostly by subcutaneous route (64-71%). Mean serum ferritin and transferrin saturation were 435 +/- 253 microg/L and 30 +/- 11%, respectively. At month 12, adequate iron stores were found in 72.5% of patients, whereas absolute and functional iron deficiencies were observed in only 5.1% and 17.8%, respectively. Multivariate analysis showed that diabetes unexpectedly influenced Hb towards higher levels (12.1 +/- 0.9 g/dL; p = 0.02). One year survival was significantly higher in patients with Hb >or= 11 g/dL than in those with Hb <11 g/dL (19.7% vs 7.3%, p = 0.006). CONCLUSION: In comparison to European studies of reference, this survey shows a remarkable and continuous control of anaemia in Swiss dialysis centres. These results were reached through moderately high EPO doses, mostly given subcutaneously, and careful iron therapy management.

Optimal growth strategies when mortality and production rates are size-dependent

Pontryagin's maximum principle from optimal control theory is used to find the optimal allocation of energy between growth and reproduction when lifespan may be finite and the trade-off between growth and reproduction is linear. Analyses of the optimal allocation problem to date have generally yielded bang-bang solutions, i.e. determinate growth: life-histories in which growth is followed by reproduction, with no intermediate phase of simultaneous reproduction and growth. Here we show that an intermediate strategy (indeterminate growth) can be selected for if the rates of production and mortality either both increase or both decrease with increasing body size, this arises as a singular solution to the problem. Our conclusion is that indeterminate growth is optimal in more cases than was previously realized. The relevance of our results to natural situations is discussed.

Early olfactory involvement in Alzheimer's disease.

BACKGROUND: In Alzheimer's disease (AD) the olfactory system, including the olfactory bulb, a limbic paleocortex is severely damaged. The occurrence of early olfactory deficits and the presence of senile plaques and neurofibrillary tangles in olfactory bulb were reported previously by a few authors. The goal of the present study was to analyze the occurrence of AD-type degenerative changes in the peripheral part of the olfactory system and to answer the question whether the frequency and severity of changes in the olfactory bulb and tract are associated with those of the cerebral cortex in AD. MATERIAL AND METHODS: In 110 autopsy cases several cortical areas and the olfactory bulb and tract were analyzed using histo- and immunohistochemical techniques. Based on a semiquantitative analysis of cortical senile plaques, neurofibrillary tangles and curly fibers, the 110 cases were divided into four groups: 19 cases with severe (definite AD), 14 cases with moderate, 58 cases with discrete and 19 control cases without AD-type cortical changes. RESULTS: The number of cases with olfactory involvement was very high, more than 84% in the three groups with cortical AD-type lesions. Degenerative olfactory changes were present in all 19 definite AD cases, and in two of the 19 controls. The statistical analysis showed a significant association between the peripheral olfactory and cortical degenerative changes with respect to their frequency and severity (P < 0.001). Neurofibrillary tangles and neuropil threads appear in the olfactory system as early as in entorhinal cortex. CONCLUSION: The results indicate a close relationship between the olfactory and cortical degenerative changes and indicate that the involvement of the olfactory bulb and tract is one of the earliest events in the degenerative process of the central nervous system in AD.

A discrete forest in a continuous landscape : investigating interactions between natural resource management and sustainable development

Summary Forests are key ecosystems of the earth and associated with a large range of functions. Many of these functions are beneficial to humans and are referred to as ecosystem services. Sustainable development requires that all relevant ecosystem services are quantified, managed and monitored equally. Natural resource management therefore targets the services associated with ecosystems. The main hypothesis of this thesis is that the spatial and temporal domains of relevant services do not correspond to a discrete forest ecosystem. As a consequence, the services are not quantified, managed and monitored in an equal and sustainable manner. The thesis aims were therefore to test this hypothesis, establish an improved conceptual approach and provide spatial applications for the relevant land cover and structure variables. The study was carried out in western Switzerland and based primarily on data from a countrywide landscape inventory. This inventory is part of the third Swiss national forest inventory and assesses continuous landscape variables based on a regular sampling of true colour aerial imagery. In addition, land cover variables were derived from Landsat 5 TM passive sensor data and land structure variables from active sensor data from a small footprint laserscanning system. The results confirmed the main hypothesis, as relevant services did not scale well with the forest ecosystem. Instead, a new conceptual approach for sustainable management of natural resources was described. This concept quantifies the services as a continuous function of the landscape, rather than a discrete function of the forest ecosystem. The explanatory landscape variables are therefore called continuous fields and the forest becomes a dependent and function-driven management unit. Continuous field mapping methods were established for land cover and structure variables. In conclusion, the discrete forest ecosystem is an adequate planning and management unit. However, monitoring the state of and trends in sustainability of services requires them to be quantified as a continuous function of the landscape. Sustainable natural resource management iteratively combines the ecosystem and gradient approaches. Résumé Les forêts sont des écosystèmes-clés de la terre et on leur attribue un grand nombre de fonctions. Beaucoup de ces fonctions sont bénéfiques pour l'homme et sont nommées services écosystémiques. Le développement durable exige que ces services écosystémiques soient tous quantifiés, gérés et surveillés de façon égale. La gestion des ressources naturelles a donc pour cible les services attribués aux écosystèmes. L'hypothèse principale de cette thèse est que les domaines spatiaux et temporels des services attribués à la forêt ne correspondent pas à un écosystème discret. Par conséquent, les services ne sont pas quantifiés, aménagés et surveillés d'une manière équivalente et durable. Les buts de la thèse étaient de tester cette hypothèse, d'établir une nouvelle approche conceptuelle de la gestion des ressources naturelles et de préparer des applications spatiales pour les variables paysagères et structurelles appropriées. L'étude a été menée en Suisse occidentale principalement sur la base d'un inventaire de paysage à l'échelon national. Cet inventaire fait partie du troisième inventaire forestier national suisse et mesure de façon continue des variables paysagères sur la base d'un échantillonnage régulier sur des photos aériennes couleur. En outre, des variables de couverture ? terrestre ont été dérivées des données d'un senseur passif Landsat 5 TM, ainsi que des variables structurelles, dérivées du laserscanning, un senseur actif. Les résultats confirment l'hypothèse principale, car l'échelle des services ne correspond pas à celle de l'écosystème forestier. Au lieu de cela, une nouvelle approche a été élaborée pour la gestion durable des ressources naturelles. Ce concept représente les services comme une fonction continue du paysage, plutôt qu'une fonction discrète de l'écosystème forestier. En conséquence, les variables explicatives de paysage sont dénommées continuous fields et la forêt devient une entité dépendante, définie par la fonction principale du paysage. Des méthodes correspondantes pour la couverture terrestre et la structure ont été élaborées. En conclusion, l'écosystème forestier discret est une unité adéquate pour la planification et la gestion. En revanche, la surveillance de la durabilité de l'état et de son évolution exige que les services soient quantifiés comme fonction continue du paysage. La gestion durable des ressources naturelles joint donc l'approche écosystémique avec celle du gradient de manière itérative.

Effects of brood manipulation costs on optimal sex allocation in social hymenoptera.

In eusocial Hymenoptera, queens and workers are in conflict over optimal sex allocation. Sex ratio theory, while generating predictions on the extent of this conflict under a wide range of conditions, has largely neglected the fact that worker control of investment almost certainly requires the manipulation of brood sex ratio. This manipulation is likely to incur costs, for example, if workers eliminate male larvae or rear more females as sexuals rather than workers. In this article, we present a model of sex ratio evolution under worker control that incorporates costs of brood manipulation. We assume cost to be a continuous, increasing function of the magnitude of sex ratio manipulation. We demonstrate that costs counterselect sex ratio biasing, which leads to less female-biased population sex ratios than expected on the basis of relatedness asymmetry. Furthermore, differently shaped cost functions lead to different equilibria of manipulation at the colony level. While linear and accelerating cost functions generate monomorphic equilibria, decelerating costs lead to a process of evolutionary branching and hence split sex ratios.

Therapeutic strategies to improve control of hypertension.

Blood pressure is poorly controlled in most European countries and the control rate is even lower in high-risk patients such as patients with chronic kidney disease, diabetic patients or previous coronary heart disease. Several factors have been associated with poor control, some of which involve the characteristic of the patients themselves, such as socioeconomic factors, or unsuitable life-styles, other factors related to hypertension or to associated comorbidity, but there are also factors directly associated with antihypertensive therapy, mainly involving adherence problems, therapeutic inertia and therapeutic strategies unsuited to difficult-to-control hypertensive patients.It is common knowledge that only 30% of hypertensive patients can be controlled using monotherapy; all the rest require a combination of two or more antihypertensive drugs, and this can be a barrier to good adherence and log-term persistence in patients who also often need to use other drugs, such as antidiabetic agents, statins or antiplatelet agents. The fixed combinations of three antihypertensive agents currently available can facilitate long-term control of these patients in clinical practice. If well tolerated, a long-term therapeutic regimen that includes a diuretic, an ACE inhibitor or an angiotensin receptor blocker, and a calcium channel blocker is the recommended optimal triple therapy.

Continuous determination of optimal cerebral perfusion pressure in traumatic brain injury.

OBJECTIVES: We have sought to develop an automated methodology for the continuous updating of optimal cerebral perfusion pressure (CPPopt) for patients after severe traumatic head injury, using continuous monitoring of cerebrovascular pressure reactivity. We then validated the CPPopt algorithm by determining the association between outcome and the deviation of actual CPP from CPPopt. DESIGN: Retrospective analysis of prospectively collected data. SETTING: Neurosciences critical care unit of a university hospital. PATIENTS: A total of 327 traumatic head-injury patients admitted between 2003 and 2009 with continuous monitoring of arterial blood pressure and intracranial pressure. MEASUREMENTS AND MAIN RESULTS: Arterial blood pressure, intracranial pressure, and CPP were continuously recorded, and pressure reactivity index was calculated online. Outcome was assessed at 6 months. An automated curve fitting method was applied to determine CPP at the minimum value for pressure reactivity index (CPPopt). A time trend of CPPopt was created using a moving 4-hr window, updated every minute. Identification of CPPopt was, on average, feasible during 55% of the whole recording period. Patient outcome correlated with the continuously updated difference between median CPP and CPPopt (chi-square=45, p<.001; outcome dichotomized into fatal and nonfatal). Mortality was associated with relative "hypoperfusion" (CPP<CPPopt), severe disability with "hyperperfusion" (CPP>CPPopt), and favorable outcome was associated with smaller deviations of CPP from the individualized CPPopt. While deviations from global target CPP values of 60 mm Hg and 70 mm Hg were also related to outcome, these relationships were less robust. CONCLUSIONS: Real-time CPPopt could be identified during the recording time of majority of the patients. Patients with a median CPP close to CPPopt were more likely to have a favorable outcome than those in whom median CPP was widely different from CPPopt. Deviations from individualized CPPopt were more predictive of outcome than deviations from a common target CPP. CPP management to optimize cerebrovascular pressure reactivity should be the subject of future clinical trial in severe traumatic head-injury patients.

Evolutionarily stable learning schedules and cumulative culture in discrete generation models.

Individual learning (e.g., trial-and-error) and social learning (e.g., imitation) are alternative ways of acquiring and expressing the appropriate phenotype in an environment. The optimal choice between using individual learning and/or social learning may be dictated by the life-stage or age of an organism. Of special interest is a learning schedule in which social learning precedes individual learning, because such a schedule is apparently a necessary condition for cumulative culture. Assuming two obligatory learning stages per discrete generation, we obtain the evolutionarily stable learning schedules for the three situations where the environment is constant, fluctuates between generations, or fluctuates within generations. During each learning stage, we assume that an organism may target the optimal phenotype in the current environment by individual learning, and/or the mature phenotype of the previous generation by oblique social learning. In the absence of exogenous costs to learning, the evolutionarily stable learning schedules are predicted to be either pure social learning followed by pure individual learning ("bang-bang" control) or pure individual learning at both stages ("flat" control). Moreover, we find for each situation that the evolutionarily stable learning schedule is also the one that optimizes the learned phenotype at equilibrium.

Artificial evolution on network structures : how time and space influence dynamics

Abstract The main objective of this work is to show how the choice of the temporal dimension and of the spatial structure of the population influences an artificial evolutionary process. In the field of Artificial Evolution we can observe a common trend in synchronously evolv¬ing panmictic populations, i.e., populations in which any individual can be recombined with any other individual. Already in the '90s, the works of Spiessens and Manderick, Sarma and De Jong, and Gorges-Schleuter have pointed out that, if a population is struc¬tured according to a mono- or bi-dimensional regular lattice, the evolutionary process shows a different dynamic with respect to the panmictic case. In particular, Sarma and De Jong have studied the selection pressure (i.e., the diffusion of a best individual when the only selection operator is active) induced by a regular bi-dimensional structure of the population, proposing a logistic modeling of the selection pressure curves. This model supposes that the diffusion of a best individual in a population follows an exponential law. We show that such a model is inadequate to describe the process, since the growth speed must be quadratic or sub-quadratic in the case of a bi-dimensional regular lattice. New linear and sub-quadratic models are proposed for modeling the selection pressure curves in, respectively, mono- and bi-dimensional regu¬lar structures. These models are extended to describe the process when asynchronous evolutions are employed. Different dynamics of the populations imply different search strategies of the resulting algorithm, when the evolutionary process is used to solve optimisation problems. A benchmark of both discrete and continuous test problems is used to study the search characteristics of the different topologies and updates of the populations. In the last decade, the pioneering studies of Watts and Strogatz have shown that most real networks, both in the biological and sociological worlds as well as in man-made structures, have mathematical properties that set them apart from regular and random structures. In particular, they introduced the concepts of small-world graphs, and they showed that this new family of structures has interesting computing capabilities. Populations structured according to these new topologies are proposed, and their evolutionary dynamics are studied and modeled. We also propose asynchronous evolutions for these structures, and the resulting evolutionary behaviors are investigated. Many man-made networks have grown, and are still growing incrementally, and explanations have been proposed for their actual shape, such as Albert and Barabasi's preferential attachment growth rule. However, many actual networks seem to have undergone some kind of Darwinian variation and selection. Thus, how these networks might have come to be selected is an interesting yet unanswered question. In the last part of this work, we show how a simple evolutionary algorithm can enable the emrgence o these kinds of structures for two prototypical problems of the automata networks world, the majority classification and the synchronisation problems. Synopsis L'objectif principal de ce travail est de montrer l'influence du choix de la dimension temporelle et de la structure spatiale d'une population sur un processus évolutionnaire artificiel. Dans le domaine de l'Evolution Artificielle on peut observer une tendence à évoluer d'une façon synchrone des populations panmictiques, où chaque individu peut être récombiné avec tout autre individu dans la population. Déjà dans les année '90, Spiessens et Manderick, Sarma et De Jong, et Gorges-Schleuter ont observé que, si une population possède une structure régulière mono- ou bi-dimensionnelle, le processus évolutionnaire montre une dynamique différente de celle d'une population panmictique. En particulier, Sarma et De Jong ont étudié la pression de sélection (c-à-d la diffusion d'un individu optimal quand seul l'opérateur de sélection est actif) induite par une structure régulière bi-dimensionnelle de la population, proposant une modélisation logistique des courbes de pression de sélection. Ce modèle suppose que la diffusion d'un individu optimal suit une loi exponentielle. On montre que ce modèle est inadéquat pour décrire ce phénomène, étant donné que la vitesse de croissance doit obéir à une loi quadratique ou sous-quadratique dans le cas d'une structure régulière bi-dimensionnelle. De nouveaux modèles linéaires et sous-quadratique sont proposés pour des structures mono- et bi-dimensionnelles. Ces modèles sont étendus pour décrire des processus évolutionnaires asynchrones. Différentes dynamiques de la population impliquent strategies différentes de recherche de l'algorithme résultant lorsque le processus évolutionnaire est utilisé pour résoudre des problèmes d'optimisation. Un ensemble de problèmes discrets et continus est utilisé pour étudier les charactéristiques de recherche des différentes topologies et mises à jour des populations. Ces dernières années, les études de Watts et Strogatz ont montré que beaucoup de réseaux, aussi bien dans les mondes biologiques et sociologiques que dans les structures produites par l'homme, ont des propriétés mathématiques qui les séparent à la fois des structures régulières et des structures aléatoires. En particulier, ils ont introduit la notion de graphe sm,all-world et ont montré que cette nouvelle famille de structures possède des intéressantes propriétés dynamiques. Des populations ayant ces nouvelles topologies sont proposés, et leurs dynamiques évolutionnaires sont étudiées et modélisées. Pour des populations ayant ces structures, des méthodes d'évolution asynchrone sont proposées, et la dynamique résultante est étudiée. Beaucoup de réseaux produits par l'homme se sont formés d'une façon incrémentale, et des explications pour leur forme actuelle ont été proposées, comme le preferential attachment de Albert et Barabàsi. Toutefois, beaucoup de réseaux existants doivent être le produit d'un processus de variation et sélection darwiniennes. Ainsi, la façon dont ces structures ont pu être sélectionnées est une question intéressante restée sans réponse. Dans la dernière partie de ce travail, on montre comment un simple processus évolutif artificiel permet à ce type de topologies d'émerger dans le cas de deux problèmes prototypiques des réseaux d'automates, les tâches de densité et de synchronisation.

Spatiotemporal neuromodulation therapies engaging muscle synergies improve motor control after spinal cord injury.

Electrical neuromodulation of lumbar segments improves motor control after spinal cord injury in animal models and humans. However, the physiological principles underlying the effect of this intervention remain poorly understood, which has limited the therapeutic approach to continuous stimulation applied to restricted spinal cord locations. Here we developed stimulation protocols that reproduce the natural dynamics of motoneuron activation during locomotion. For this, we computed the spatiotemporal activation pattern of muscle synergies during locomotion in healthy rats. Computer simulations identified optimal electrode locations to target each synergy through the recruitment of proprioceptive feedback circuits. This framework steered the design of spatially selective spinal implants and real-time control software that modulate extensor and flexor synergies with precise temporal resolution. Spatiotemporal neuromodulation therapies improved gait quality, weight-bearing capacity, endurance and skilled locomotion in several rodent models of spinal cord injury. These new concepts are directly translatable to strategies to improve motor control in humans.

Periodontitis and gingivitis in inflammatory bowel disease: a case-control study.

BACKGROUND: The oral cavity is frequently affected in patients with inflammatory bowel disease (IBD), especially in patients with Crohn's disease (CD). Periodontitis is thought to influence systemic autoimmune or inflammatory diseases. We aimed to analyze the relationship of periodontitis and gingivitis markers with specific disease characteristics in patients with IBD and to compare these data with healthy controls. METHODS: In a prospective 8-month study, systematic oral examinations were performed in 113 patients with IBD, including 69 patients with CD and 44 patients with ulcerative colitis. For all patients, a structured personal history was taken. One hundred thirteen healthy volunteers served as a control group. Oral examination focussed on established oral health markers for periodontitis (bleeding on probing, loss of attachment, and periodontal pocket depth) and gingivitis (papilla bleeding index). Additionally, visible oral lesions were documented. RESULTS: Both gingivitis and periodontitis markers were higher in patients with IBD than in healthy control. In univariate analysis and logistic regression analysis, perianal disease was a risk factor for periodontitis. Nonsmoking decreased the risk of having periodontitis. No clear association was found between clinical activity and periodontitis in IBD. In only the CD subgroup, high clinical activity (Harvey-Bradshaw index > 10) was associated with 1 periodontitis marker, the loss of attachment at sites of maximal periodontal pocket depth. Oral lesions besides periodontitis and gingivitis were not common, but nevertheless observed in about 10% of patients with IBD. CONCLUSIONS: IBD, and especially perianal disease in CD, is associated with periodontitis. Optimal therapeutic strategies should probably focus on treating both local oral and systemic inflammation.

A prospective randomised multi-centre controlled trial on tight glucose control by intensive insulin therapy in adult intensive care units: the Glucontrol study.

PURPOSE: An optimal target for glucose control in ICU patients remains unclear. This prospective randomized controlled trial compared the effects on ICU mortality of intensive insulin therapy (IIT) with an intermediate glucose control. METHODS: Adult patients admitted to the 21 participating medico-surgical ICUs were randomized to group 1 (target BG 7.8-10.0 mmol/L) or to group 2 (target BG 4.4-6.1 mmol/L). RESULTS: While the required sample size was 1,750 per group, the trial was stopped early due to a high rate of unintended protocol violations. From 1,101 admissions, the outcomes of 542 patients assigned to group 1 and 536 of group 2 were analysed. The groups were well balanced. BG levels averaged in group 1 8.0 mmol/L (IQR 7.1-9.0) (median of all values) and 7.7 mmol/L (IQR 6.7-8.8) (median of morning BG) versus 6.5 mmol/L (IQR 6.0-7.2) and 6.1 mmol/L (IQR 5.5-6.8) for group 2 (p < 0.0001 for both comparisons). The percentage of patients treated with insulin averaged 66.2 and 96.3%, respectively. Proportion of time spent in target BG was similar, averaging 39.5% and 45.1% (median (IQR) 34.3 (18.5-50.0) and 39.3 (26.2-53.6)%) in the groups 1 and 2, respectively. The rate of hypoglycaemia was higher in the group 2 (8.7%) than in group 1 (2.7%, p < 0.0001). ICU mortality was similar in the two groups (15.3 vs. 17.2%). CONCLUSIONS: In this prematurely stopped and therefore underpowered study, there was a lack of clinical benefit of intensive insulin therapy (target 4.4-6.1 mmol/L), associated with an increased incidence of hypoglycaemia, as compared to a 7.8-10.0 mmol/L target. (ClinicalTrials.gov # NCT00107601, EUDRA-CT Number: 200400391440).

Zerebrale Probleme in der pädiatrischen Intensivmedizin [Cerebral problems in pediatric intensive care]

This paper presents a review of different methods enabling the monitoring of cerebral function in neonatal and paediatric intensive care. EEG, evoked potentials, conventional radiological studies, computerized tomography, ultrasound, intracranial pressure measurements, nuclear magnetic resonance, Doppler ultrasound, radioisotope studies, angiography, infra-red spectral analysis and last, but not least, clinical examination produce information regarding the neurological state of the patient which must be critically analysed in order to ensure optimal management of the case. Unfortunately, and in spite of impressive progress in non-invasive monitoring of the cerebral function, we are still forced to make important medical and ethical decisions without precise information about the neurological state of our patients.

Dynamic simulation of regulatory networks using SQUAD.

BACKGROUND: The ambition of most molecular biologists is the understanding of the intricate network of molecular interactions that control biological systems. As scientists uncover the components and the connectivity of these networks, it becomes possible to study their dynamical behavior as a whole and discover what is the specific role of each of their components. Since the behavior of a network is by no means intuitive, it becomes necessary to use computational models to understand its behavior and to be able to make predictions about it. Unfortunately, most current computational models describe small networks due to the scarcity of kinetic data available. To overcome this problem, we previously published a methodology to convert a signaling network into a dynamical system, even in the total absence of kinetic information. In this paper we present a software implementation of such methodology. RESULTS: We developed SQUAD, a software for the dynamic simulation of signaling networks using the standardized qualitative dynamical systems approach. SQUAD converts the network into a discrete dynamical system, and it uses a binary decision diagram algorithm to identify all the steady states of the system. Then, the software creates a continuous dynamical system and localizes its steady states which are located near the steady states of the discrete system. The software permits to make simulations on the continuous system, allowing for the modification of several parameters. Importantly, SQUAD includes a framework for perturbing networks in a manner similar to what is performed in experimental laboratory protocols, for example by activating receptors or knocking out molecular components. Using this software we have been able to successfully reproduce the behavior of the regulatory network implicated in T-helper cell differentiation. CONCLUSION: The simulation of regulatory networks aims at predicting the behavior of a whole system when subject to stimuli, such as drugs, or determine the role of specific components within the network. The predictions can then be used to interpret and/or drive laboratory experiments. SQUAD provides a user-friendly graphical interface, accessible to both computational and experimental biologists for the fast qualitative simulation of large regulatory networks for which kinetic data is not necessarily available.