Analysis of c-Myc function in mouse epidermis

Abstract The c-myc gene is one of the most frequently mutated oncogenes found in human tumors. c-Myc has been implicated in the regulation of various biological processes including cell cycle progression, cellular growth, differentiation, angiogenesis, immortalization and apoptosis. To assess the normal role of c-Myc in epithelial cell types in vitro and in vivo we have deleted the c-myc gene in keratinocytes and in the adult skin epidermis by conditional Cre/loxP mediated recombination. Similar to what we have previously shown in mouse embryonic fibroblasts acute elimination of c-Myc activity in cultured keratinocytes causes cells to cease proliferation and adapt a flat cell morphology. Mutant cells accumulate in a diploid Ki67neg stage, indicative of a quiescent Go stage. This demonstrates that c-Myc activity is essential to maintain keratinocytes in a productive cell cycle. In addition, mutant keratinocytes showed a defect in Ca2+ induced induction of the differentiation marker Keratin 1 suggesting a role for c-Myc during differentiation. To assess the in vivo role of c-Myc we used a tamoxifen inducible K5::CreERT transgene to delete the c-myc gene in the adult skin epidermis. Unexpectedly, despite strong c-Myc expression in the basal compartment it is not required for maintenance of the skin epidermis in the adult mouse. The epidermis appeared normal with respect to both proliferation and differentiation. In addition, no selection against c-Myc deficient epidermal cells occurred over many months, further confirming that c-Myc is dispensable for normal skin homeostasis. Even more surprising, TPA induced hyperproliferation also occurred in a c-Myc independent manner. Treatment of the skin with the mutagen DMBA prior to TPA is a classical way to induce papillomas by selecting for mutations that lead to dominant activation of the oncogene Ha-Ras. Most interestingly tumor formation was severely inhibited suggesting that tumor progression requires endogenous c-Myc. Further studies are required to address whether the role of c-Myc in the activation of telomerase or the Werner protein, or its role to induce angiogenesis is required for skin tumor progression, In conclusion, this work shows that while c-Myc is not required for maintenance or hyperplasia of mouse epidermis, it is essential for skin tumor progression in collaboration with Ras. Résumé Le gène c-myc est un des oncogènes les plus fréquemment mutés dans les tumeurs humaines. c-Myc est impliqué dans la régulation de processus biologiques variés, comme la progression du cycle cellulaire, la croissance cellulaire, la différenciation, l'angiogenèse, l'immortalisation et l'apoptose. Pour caractériser le rôle physiologique de c-Myc dans les cellules de type épithélial in vitro et in vivo, le gène c-myc a été délété dans des kératinocytes primaires et dans l'épiderme de peau de souris adultes par des recombinaisons conditionnelles (système Cre/loxP). De la même façon que dans les fibroblastes d'embryon de souris, l'élimination aiguë de l'activité de c-Myc dans les kératinocytes en culture primaire provoque l'arrêt de la prolifération des cellules et leur applatissement morphologique. Les cellules mutantes restent dans un stade diploïde Ki67neg, indiquant un stade quiescent Go. Cela démontre que l'activité de c-Myc est essentielle pour maintenir les kératinocytes dans le cycle cellulaire. De plus, les kératinocytes mutants montrent une déficience pour le marqueur de différenciation Kératine 1 au cours de la différenciation induite par le calcium, suggérant un rôle de c-Myc dans la différenciation cellulaire. Pour comprendre le rôle de c-Myc in vivo, le transgène K5::CreERT inductible par le tamoxifen a été utilisé pour déléter le gène c-inyc dans l'épiderme de souris adultes. Etonnemment, malgré une forte expression de c-Myc dans le compartiment basal de l'épiderme, ce gène n'est pas nécessaire pour la maintenance de l'épiderme de la peau chez la souris adulte. L'épiderme apparait normal avec une prolifération et une différenciation physiologique des cellules. De plus, il n'y a pas de sélection contre les cellules épidennales c-Myc déficientes après plusieurs mois, ce qui confirme que c-Myc n'est pas nécessaire pour l'homéostasie normale de la peau. Encore plus surprenant, une hyperprolifération est également induite par du TPA chez les souris mutantes, impliquant une voie de prolifération indépendante de c-Myc. Le traitement de la peau par le mutagène DMBA avant le traitement au TPA est une voie classique d'induction de papillomes, par sélection de mutations conduisant à l'activation de l'oncogène Ha-Ras. La formation des tumeurs est fortement inhibée chez les souris mutantes, suggérant que la progression des tumeurs nécessite la présence endogène de c-Myc. De nouvelles études sont nécessaires pour savoir si c-Myc a un rôle dans l'activation de la télomérase ou de la protéine de Werner, ou encore dans l'angiogénèse, qui sont nécessaires pour la progression tumorale. En conclusion, ce travail montre que même si c-Myc n'est pas nécessaire pour la maintenance ou l'hyperplasie de la peau de souris, il est essentiel pour la progression des tumeurs de la peau en collaboration avec Ras.

Behavioural queueing : cellular automata and a laboratory experiment

Dans cette thèse, nous étudions les aspects comportementaux d'agents qui interagissent dans des systèmes de files d'attente à l'aide de modèles de simulation et de méthodologies expérimentales. Chaque période les clients doivent choisir un prestataire de servivce. L'objectif est d'analyser l'impact des décisions des clients et des prestataires sur la formation des files d'attente. Dans un premier cas nous considérons des clients ayant un certain degré d'aversion au risque. Sur la base de leur perception de l'attente moyenne et de la variabilité de cette attente, ils forment une estimation de la limite supérieure de l'attente chez chacun des prestataires. Chaque période, ils choisissent le prestataire pour lequel cette estimation est la plus basse. Nos résultats indiquent qu'il n'y a pas de relation monotone entre le degré d'aversion au risque et la performance globale. En effet, une population de clients ayant un degré d'aversion au risque intermédiaire encoure généralement une attente moyenne plus élevée qu'une population d'agents indifférents au risque ou très averses au risque. Ensuite, nous incorporons les décisions des prestataires en leur permettant d'ajuster leur capacité de service sur la base de leur perception de la fréquence moyenne d'arrivées. Les résultats montrent que le comportement des clients et les décisions des prestataires présentent une forte "dépendance au sentier". En outre, nous montrons que les décisions des prestataires font converger l'attente moyenne pondérée vers l'attente de référence du marché. Finalement, une expérience de laboratoire dans laquelle des sujets jouent le rôle de prestataire de service nous a permis de conclure que les délais d'installation et de démantèlement de capacité affectent de manière significative la performance et les décisions des sujets. En particulier, les décisions du prestataire, sont influencées par ses commandes en carnet, sa capacité de service actuellement disponible et les décisions d'ajustement de capacité qu'il a prises, mais pas encore implémentées. - Queuing is a fact of life that we witness daily. We all have had the experience of waiting in line for some reason and we also know that it is an annoying situation. As the adage says "time is money"; this is perhaps the best way of stating what queuing problems mean for customers. Human beings are not very tolerant, but they are even less so when having to wait in line for service. Banks, roads, post offices and restaurants are just some examples where people must wait for service. Studies of queuing phenomena have typically addressed the optimisation of performance measures (e.g. average waiting time, queue length and server utilisation rates) and the analysis of equilibrium solutions. The individual behaviour of the agents involved in queueing systems and their decision making process have received little attention. Although this work has been useful to improve the efficiency of many queueing systems, or to design new processes in social and physical systems, it has only provided us with a limited ability to explain the behaviour observed in many real queues. In this dissertation we differ from this traditional research by analysing how the agents involved in the system make decisions instead of focusing on optimising performance measures or analysing an equilibrium solution. This dissertation builds on and extends the framework proposed by van Ackere and Larsen (2004) and van Ackere et al. (2010). We focus on studying behavioural aspects in queueing systems and incorporate this still underdeveloped framework into the operations management field. In the first chapter of this thesis we provide a general introduction to the area, as well as an overview of the results. In Chapters 2 and 3, we use Cellular Automata (CA) to model service systems where captive interacting customers must decide each period which facility to join for service. They base this decision on their expectations of sojourn times. Each period, customers use new information (their most recent experience and that of their best performing neighbour) to form expectations of sojourn time at the different facilities. Customers update their expectations using an adaptive expectations process to combine their memory and their new information. We label "conservative" those customers who give more weight to their memory than to the xiv Summary new information. In contrast, when they give more weight to new information, we call them "reactive". In Chapter 2, we consider customers with different degree of risk-aversion who take into account uncertainty. They choose which facility to join based on an estimated upper-bound of the sojourn time which they compute using their perceptions of the average sojourn time and the level of uncertainty. We assume the same exogenous service capacity for all facilities, which remains constant throughout. We first analyse the collective behaviour generated by the customers' decisions. We show that the system achieves low weighted average sojourn times when the collective behaviour results in neighbourhoods of customers loyal to a facility and the customers are approximately equally split among all facilities. The lowest weighted average sojourn time is achieved when exactly the same number of customers patronises each facility, implying that they do not wish to switch facility. In this case, the system has achieved the Nash equilibrium. We show that there is a non-monotonic relationship between the degree of risk-aversion and system performance. Customers with an intermediate degree of riskaversion typically achieve higher sojourn times; in particular they rarely achieve the Nash equilibrium. Risk-neutral customers have the highest probability of achieving the Nash Equilibrium. Chapter 3 considers a service system similar to the previous one but with risk-neutral customers, and relaxes the assumption of exogenous service rates. In this sense, we model a queueing system with endogenous service rates by enabling managers to adjust the service capacity of the facilities. We assume that managers do so based on their perceptions of the arrival rates and use the same principle of adaptive expectations to model these perceptions. We consider service systems in which the managers' decisions take time to be implemented. Managers are characterised by a profile which is determined by the speed at which they update their perceptions, the speed at which they take decisions, and how coherent they are when accounting for their previous decisions still to be implemented when taking their next decision. We find that the managers' decisions exhibit a strong path-dependence: owing to the initial conditions of the model, the facilities of managers with identical profiles can evolve completely differently. In some cases the system becomes "locked-in" into a monopoly or duopoly situation. The competition between managers causes the weighted average sojourn time of the system to converge to the exogenous benchmark value which they use to estimate their desired capacity. Concerning the managers' profile, we found that the more conservative Summary xv a manager is regarding new information, the larger the market share his facility achieves. Additionally, the faster he takes decisions, the higher the probability that he achieves a monopoly position. In Chapter 4 we consider a one-server queueing system with non-captive customers. We carry out an experiment aimed at analysing the way human subjects, taking on the role of the manager, take decisions in a laboratory regarding the capacity of a service facility. We adapt the model proposed by van Ackere et al (2010). This model relaxes the assumption of a captive market and allows current customers to decide whether or not to use the facility. Additionally the facility also has potential customers who currently do not patronise it, but might consider doing so in the future. We identify three groups of subjects whose decisions cause similar behavioural patterns. These groups are labelled: gradual investors, lumpy investors, and random investor. Using an autocorrelation analysis of the subjects' decisions, we illustrate that these decisions are positively correlated to the decisions taken one period early. Subsequently we formulate a heuristic to model the decision rule considered by subjects in the laboratory. We found that this decision rule fits very well for those subjects who gradually adjust capacity, but it does not capture the behaviour of the subjects of the other two groups. In Chapter 5 we summarise the results and provide suggestions for further work. Our main contribution is the use of simulation and experimental methodologies to explain the collective behaviour generated by customers' and managers' decisions in queueing systems as well as the analysis of the individual behaviour of these agents. In this way, we differ from the typical literature related to queueing systems which focuses on optimising performance measures and the analysis of equilibrium solutions. Our work can be seen as a first step towards understanding the interaction between customer behaviour and the capacity adjustment process in queueing systems. This framework is still in its early stages and accordingly there is a large potential for further work that spans several research topics. Interesting extensions to this work include incorporating other characteristics of queueing systems which affect the customers' experience (e.g. balking, reneging and jockeying); providing customers and managers with additional information to take their decisions (e.g. service price, quality, customers' profile); analysing different decision rules and studying other characteristics which determine the profile of customers and managers.

Ab initio modeling and molecular dynamics simulation of the alpha 1b-adrenergic receptor activation.

This work describes the ab initio procedure employed to build an activation model for the alpha 1b-adrenergic receptor (alpha 1b-AR). The first version of the model was progressively modified and complicated by means of a many-step iterative procedure characterized by the employment of experimental validations of the model in each upgrading step. A combined simulated (molecular dynamics) and experimental mutagenesis approach was used to determine the structural and dynamic features characterizing the inactive and active states of alpha 1b-AR. The latest version of the model has been successfully challenged with respect to its ability to interpret and predict the functional properties of a large number of mutants. The iterative approach employed to describe alpha 1b-AR activation in terms of molecular structure and dynamics allows further complications of the model to allow prediction and interpretation of an ever-increasing number of experimental data.

Hydrogeophysical data integration at larger scales: Application of Bayesian sequential simulation for the characterization of heterogeneous alluvial aquifers

Knowledge of the spatial distribution of hydraulic conductivity (K) within an aquifer is critical for reliable predictions of solute transport and the development of effective groundwater management and/or remediation strategies. While core analyses and hydraulic logging can provide highly detailed information, such information is inherently localized around boreholes that tend to be sparsely distributed throughout the aquifer volume. Conversely, larger-scale hydraulic experiments like pumping and tracer tests provide relatively low-resolution estimates of K in the investigated subsurface region. As a result, traditional hydrogeological measurement techniques contain a gap in terms of spatial resolution and coverage, and they are often alone inadequate for characterizing heterogeneous aquifers. Geophysical methods have the potential to bridge this gap. The recent increased interest in the application of geophysical methods to hydrogeological problems is clearly evidenced by the formation and rapid growth of the domain of hydrogeophysics over the past decade (e.g., Rubin and Hubbard, 2005).

Effort, reward and self-reported mental health: a simulation study on negative affectivity bias.

Background : In the present article, we propose an alternative method for dealing with negative affectivity (NA) biases in research, while investigating the association between a deleterious psychosocial environment at work and poor mental health. First, we investigated how strong NA must be to cause an observed correlation between the independent and dependent variables. Second, we subjectively assessed whether NA can have a large enough impact on a large enough number of subjects to invalidate the observed correlations between dependent and independent variables.Methods : We simulated 10,000 populations of 300 subjects each, using the marginal distribution of workers in an actual population that had answered the Siegrist's questionnaire on effort and reward imbalance (ERI) and the General Health Questionnaire (GHQ).Results : The results of the present study suggested that simulated NA has a minimal effect on the mean scores for effort and reward. However, the correlations between the effort and reward imbalance (ERI) ratio and the GHQ score might be important, even in simulated populations with a limited NA.Conclusions : When investigating the relationship between the ERI ratio and the GHQ score, we suggest the following rules for the interpretation of the results: correlations with an explained variance of 5% and below should be considered with caution; correlations with an explained variance between 5% and 10% may result from NA, although this effect does not seem likely; and correlations with an explained variance of 10% and above are not likely to be the result of NA biases. [Authors]

Microsatellites can be misleading: an empirical and simulation study.

It has been long recognized that highly polymorphic genetic markers can lead to underestimation of divergence between populations when migration is low. Microsatellite loci, which are characterized by extremely high mutation rates, are particularly likely to be affected. Here, we report genetic differentiation estimates in a contact zone between two chromosome races of the common shrew (Sorex araneus), based on 10 autosomal microsatellites, a newly developed Y-chromosome microsatellite, and mitochondrial DNA. These results are compared to previous data on proteins and karyotypes. Estimates of genetic differentiation based on F- and R-statistics are much lower for autosomal microsatellites than for all other genetic markers. We show by simulations that this discrepancy stems mainly from the high mutation rate of microsatellite markers for F-statistics and from deviations from a single-step mutation model for R-statistics. The sex-linked genetic markers show that all gene exchange between races is mediated by females. The absence of male-mediated gene flow most likely results from male hybrid sterility.

Treatment of essential hypertension with calcium channel blockers: what is the place of lercanidipine?

In all actual clinical guidelines, dihydropyridine calcium channel blockers (CCBs) belong to the recommended first line antihypertensive drugs to treat essential hypertension. Several recent large clinical trials have confirmed their efficacy not only in lowering blood pressure but also in reducing cardiovascular morbidity and mortality in hypertensive patients with a normal or high cardiovascular risk profile. In clinical trials such as ALLHAT, VALUE or ASCOT, an amlodipine-based therapy was at least as effective, when not slightly superior, in lowering blood pressure and sometimes more effective in preventing target organ damages than blood pressure lowering strategies based on the use of diuretics, beta-blockers and blockers of the renin-angiotensin system. One of the main clinical side effects of the first and second generation CCBs including amlodipine is the development of peripheral edema. The incidence of leg edema can be markedly reduced by combining the CCB with a blocker of the renin-angiotensin system. This strategy has now led to the development of several fixed-dose combinations of amlodipine and angiotensin II receptor antagonists. Another alternative to lower the incidence of edema is to use CCBs of the third generation such as lercanidipine. Indeed, although no major clinical trials have been conducted with this compound, clinical studies have shown that lercanidipine and amlodipine have a comparable antihypertensive efficacy but with significantly less peripheral edema in patients receiving lercanidipine. In some countries, lercanidipine is now available in a single-pill association with an ACE inhibitor thereby further improving its efficacy and tolerability profile.

Cell-based computational modeling of vascular morphogenesis using Tissue Simulation Toolkit.

Computational modeling has become a widely used tool for unraveling the mechanisms of higher level cooperative cell behavior during vascular morphogenesis. However, experimenting with published simulation models or adding new assumptions to those models can be daunting for novice and even for experienced computational scientists. Here, we present a step-by-step, practical tutorial for building cell-based simulations of vascular morphogenesis using the Tissue Simulation Toolkit (TST). The TST is a freely available, open-source C++ library for developing simulations with the two-dimensional cellular Potts model, a stochastic, agent-based framework to simulate collective cell behavior. We will show the basic use of the TST to simulate and experiment with published simulations of vascular network formation. Then, we will present step-by-step instructions and explanations for building a recent simulation model of tumor angiogenesis. Demonstrated mechanisms include cell-cell adhesion, chemotaxis, cell elongation, haptotaxis, and haptokinesis.

Angiotensin-converting enzyme inhibitor versus calcium antagonist in the treatment of hypertension.

Sixteen patients with essential hypertension were treated for 2 consecutive 6-week periods with either the angiotensin-converting enzyme (ACE) inhibitor enalapril (20 mg once daily) or the calcium antagonist diltiazem (120 mg twice daily). The sequence of the treatment phases was randomly allocated. Blood pressure decreased from 154/102 +/- 5/2 mm Hg (mean +/- SEM) to 135/96 +/- 4/2 and 140/98 +/- 3/2 mm Hg during treatment with enalapril and diltiazem, respectively. It was impossible in the individual hypertensive patient to predict the long-term blood pressure response to one of the agents studied based on the long-term blood pressure response to the other agent.

Interaction between calcium intake and menarcheal age on bone mass gain: an eight-year follow-up study from prepuberty to postmenarche.

Both late menarcheal age and low calcium intake (Ca intake) during growth are risk factors for osteoporosis, probably by impairing peak bone mass. We investigated whether lasting gain in areal bone mineral density (aBMD) in response to increased Ca intake varies according to menarcheal age and, conversely, whether Ca intake could influence menarcheal age. In an initial study, 144 prepubertal girls were randomized in a double-blind controlled trial to receive either a Ca supplement (Ca-suppl.) of 850 mg/d or placebo from age 7.9-8.9 yr. Mean aBMD gain determined by dual energy x-ray absorptiometry at six sites (radius metaphysis, radius diaphysis, femoral neck, trochanter, femoral diaphysis, and L2-L4) was significantly (P = 0.004) greater in the Ca-suppl. than in the placebo group (27 vs. 21 mg/cm(2)). In 122 girls followed up, menarcheal age was recorded, and aBMD was determined at 16.4 yr of age. Menarcheal age was lower in the Ca-suppl. than in the placebo group (P = 0.048). Menarcheal age and Ca intake were negatively correlated (r = -0.35; P < 0.001), as were aBMD gains from age 7.9-16.4 yr and menarcheal age at all skeletal sites (range: r = -0.41 to r = -0.22; P < 0.001 to P = 0.016). The positive effect of Ca-suppl. on the mean aBMD gain from baseline remained significantly greater in girls below, but not in those above, the median of menarcheal age (13.0 yr). Early menarcheal age (12.1 +/- 0.5 yr): placebo, 286 +/- 36 mg/cm(2); Ca-suppl., 317 +/- 46 (P = 0.009); late menarcheal age (13.9 +/- 0.5 yr): placebo, 284 +/- 58; Ca-suppl., 276 +/- 50 (P > 0.05). The level of Ca intake during prepuberty may influence the timing of menarche, which, in turn, could influence long-term bone mass gain in response to Ca supplementation. Thus, both determinants of early menarcheal age and high Ca intake may positively interact on bone mineral mass accrual.

Simulation-based systematic review of imatinib population pharmacokinetics and PK-PD relationships in chronic myeloid leukemia (CML) patients

Objectives: Several population pharmacokinetic (PPK) and pharmacokinetic-pharmacodynamic (PK-PD) analyses have been performed with the anticancer drug imatinib. Inspired by the approach of meta-analysis, we aimed to compare and combine results from published studies in a useful way - in particular for improving the clinical interpretation of imatinib concentration measurements in the scope of therapeutic drug monitoring (TDM). Methods: Original PPK analyses and PK-PD studies (PK surrogate: trough concentration Cmin; PD outcomes: optimal early response and specific adverse events) were searched systematically on MEDLINE. From each identified PPK model, a predicted concentration distribution under standard dosage was derived through 1000 simulations (NONMEM), after standardizing model parameters to common covariates. A "reference range" was calculated from pooled simulated concentrations in a semi-quantitative approach (without specific weighting) over the whole dosing interval. Meta-regression summarized relationships between Cmin and optimal/suboptimal early treatment response. Results: 9 PPK models and 6 relevant PK-PD reports in CML patients were identified. Model-based predicted median Cmin ranged from 555 to 1388 ng/ml (grand median: 870 ng/ml and inter-quartile range: 520-1390 ng/ml). The probability to achieve optimal early response was predicted to increase from 60 to 85% from 520 to 1390 ng/ml across PK-PD studies (odds ratio for doubling Cmin: 2.7). Reporting of specific adverse events was too heterogeneous to perform a regression analysis. The general frequency of anemia, rash and fluid retention increased however consistently with Cmin, but less than response probability. Conclusions: Predicted drug exposure may differ substantially between various PPK analyses. In this review, heterogeneity was mainly attributed to 2 "outlying" models. The established reference range seems to cover the range where both good efficacy and acceptable tolerance are expected for most patients. TDM guided dose adjustment appears therefore justified for imatinib in CML patients. Its usefulness remains now to be prospectively validated in a randomized trial.

Experimental calcium oxalate crystal production and dissolution by selected wood rot fungi.

Twenty-six species of white-rotting Agaricomycotina fungi (Basidiomycota) were screened for their ability to produce calcium-oxalate (CaOx) crystals in vitro. Most were able to produce CaOx crystals in malt agar medium in the absence of additional calcium. In the same medium enriched with Ca2+, all the species produced CaOx crystals (weddellite or whewellite). Hyphae of four species (Ganoderma lucidum, Polyporus ciliatus, Pycnoporus cinnabarinus, and Trametes versicolor) were found coated with crystals (weddellite/whewellite). The production of CaOx crystals during the growth phase was confirmed by an investigation of the production kinetics for six of the species considered in the initial screening (Pleurotus citrinopileatus, Pleurotus eryngii, Pleurotus ostreatus, P. cinnabarinus, Trametes suaveolens, and T. versicolor). However, the crystals produced during the growth phase disappeared from the medium over time in four of the six species (P. citrinopileatus, P. eryngii, P. cinnabarinus, and T. suaveolens). For P. cinnabarinus, the disappearance of the crystals was correlated with a decrease in the total oxalate concentration measured in the medium from 0.65 μg mm−2 (at the maximum accumulation rate) to 0.30 μg mm−2. The decrease in the CaOx concentration was correlated with a change in mycelia morphology. The oxalate dissolution capability of all the species was also tested in a medium containing calcium oxalate as the sole source of carbon (modified Schlegel medium). Three species (Agaricus blazei, Pleurotus tuberregium, and P. ciliatus) presented a dissolution halo around the growth zone. This study shows that CaOx crystal production is a widespread phenomenon in white-rot fungi, and that an excess of Ca2+ can enhance CaOx crystal production. In addition, it shows that some white-rot fungal species are capable of dissolving CaOx crystals after growth has ceased. These results highlight a diversity of responses around the production or dissolution of calcium oxalate in white-rot fungi and reveal an unexpected potential importance of fungi on the oxalate cycle in the environment.

Artificial Snow Optimization in Winter Sport Destinations Using a Multi-agent Simulation

This paper presents the Juste-Neige system for predicting the snow height on the ski runs of a resort using a multi-agent simulation software. Its aim is to facilitate snow cover management in order to i) reduce the production cost of artificial snow and to improve the profit margin for the companies managing the ski resorts; and ii) to reduce the water and energy consumption, and thus to reduce the environmental impact, by producing only the snow needed for a good skiing experience. The software provides maps with the predicted snow heights for up to 13 days. On these maps, the areas most exposed to snow erosion are highlighted. The software proceeds in three steps: i) interpolation of snow height measurements with a neural network; ii) local meteorological forecasts for every ski resort; iii) simulation of the impact caused by skiers using a multi-agent system. The software has been evaluated in the Swiss ski resort of Verbier and provides useful predictions.