Stochastic time-concentration activity models for cytotoxicity in 3D brain cell cultures.

BACKGROUND: In vitro aggregating brain cell cultures containing all types of brain cells have been shown to be useful for neurotoxicological investigations. The cultures are used for the detection of nervous system-specific effects of compounds by measuring multiple endpoints, including changes in enzyme activities. Concentration-dependent neurotoxicity is determined at several time points. METHODS: A Markov model was set up to describe the dynamics of brain cell populations exposed to potentially neurotoxic compounds. Brain cells were assumed to be either in a healthy or stressed state, with only stressed cells being susceptible to cell death. Cells may have switched between these states or died with concentration-dependent transition rates. Since cell numbers were not directly measurable, intracellular lactate dehydrogenase (LDH) activity was used as a surrogate. Assuming that changes in cell numbers are proportional to changes in intracellular LDH activity, stochastic enzyme activity models were derived. Maximum likelihood and least squares regression techniques were applied for estimation of the transition rates. Likelihood ratio tests were performed to test hypotheses about the transition rates. Simulation studies were used to investigate the performance of the transition rate estimators and to analyze the error rates of the likelihood ratio tests. The stochastic time-concentration activity model was applied to intracellular LDH activity measurements after 7 and 14 days of continuous exposure to propofol. The model describes transitions from healthy to stressed cells and from stressed cells to death. RESULTS: The model predicted that propofol would affect stressed cells more than healthy cells. Increasing propofol concentration from 10 to 100 μM reduced the mean waiting time for transition to the stressed state by 50%, from 14 to 7 days, whereas the mean duration to cellular death reduced more dramatically from 2.7 days to 6.5 hours. CONCLUSION: The proposed stochastic modeling approach can be used to discriminate between different biological hypotheses regarding the effect of a compound on the transition rates. The effects of different compounds on the transition rate estimates can be quantitatively compared. Data can be extrapolated at late measurement time points to investigate whether costs and time-consuming long-term experiments could possibly be eliminated.

Population pharmacokinetics and clinical response for artemether-lumefantrine in pregnant and nonpregnant women with uncomplicated Plasmodium falciparum malaria in Tanzania.

Artemether-lumefantrine (AL) is the first-line treatment for uncomplicated malaria in the second and third trimesters of pregnancy. Its efficacy during pregnancy has recently been challenged due to altered pharmacokinetic (PK) properties in this vulnerable group. The aim of this study was to determine the PK profile of AL in pregnant and nonpregnant women and assess their therapeutic outcome. Thirty-three pregnant women and 22 nonpregnant women with malaria were treated with AL (80/480 mg) twice daily for 3 days. All patients provided five venous plasma samples for drug quantification at random times over 7 days. Inter- and intraindividual variability was assessed, and the effects of covariates were quantified using a nonlinear mixed-effects modeling approach (NONMEM). A one-compartment model with first-order absorption and elimination with linear metabolism from drug to metabolite fitted the data best for both arthemether (AM) and lumefantrine (LF) and their metabolites. Pregnancy status and diarrhea showed a significant influence on LF PK. The relative bioavailability of lumefantrine and its metabolism rate into desmethyl-lumefantrine were, respectively, 34% lower and 78% higher in pregnant women than in nonpregnant patients. The overall PCR-uncorrected treatment failure rates were 18% in pregnant women and 5% in nonpregnant women (odds ratio [OR] = 4.04; P value of 0.22). A high median day 7 lumefantrine concentration was significantly associated with adequate clinical and parasitological response (P = 0.03). The observed reduction in the relative bioavailability of lumefantrine in pregnant women may explain the higher treatment failure in this group, mostly due to lower posttreatment prophylaxis. Hence, a modified treatment regimen of malaria in pregnancy should be considered.

Use of Two-Part Regression Calibration Model to Correct for Measurement Error in Episodically Consumed Foods in a Single-Replicate Study Design: EPIC Case Study.

In epidemiologic studies, measurement error in dietary variables often attenuates association between dietary intake and disease occurrence. To adjust for the attenuation caused by error in dietary intake, regression calibration is commonly used. To apply regression calibration, unbiased reference measurements are required. Short-term reference measurements for foods that are not consumed daily contain excess zeroes that pose challenges in the calibration model. We adapted two-part regression calibration model, initially developed for multiple replicates of reference measurements per individual to a single-replicate setting. We showed how to handle excess zero reference measurements by two-step modeling approach, how to explore heteroscedasticity in the consumed amount with variance-mean graph, how to explore nonlinearity with the generalized additive modeling (GAM) and the empirical logit approaches, and how to select covariates in the calibration model. The performance of two-part calibration model was compared with the one-part counterpart. We used vegetable intake and mortality data from European Prospective Investigation on Cancer and Nutrition (EPIC) study. In the EPIC, reference measurements were taken with 24-hour recalls. For each of the three vegetable subgroups assessed separately, correcting for error with an appropriately specified two-part calibration model resulted in about three fold increase in the strength of association with all-cause mortality, as measured by the log hazard ratio. Further found is that the standard way of including covariates in the calibration model can lead to over fitting the two-part calibration model. Moreover, the extent of adjusting for error is influenced by the number and forms of covariates in the calibration model. For episodically consumed foods, we advise researchers to pay special attention to response distribution, nonlinearity, and covariate inclusion in specifying the calibration model.

Calibration and Validation of SWAT for the Upper Maquoketa River Watershed, June 2005

A validation study has been performed using the Soil and Water Assessment Tool (SWAT) model with data collected for the Upper Maquoketa River Watershed (UMRW), which drains over 16,000 ha in northeast Iowa. This validation assessment builds on a previous study with nested modeling for the UMRW that required both the Agricultural Policy EXtender (APEX) model and SWAT. In the nested modeling approach, edge-offield flows and pollutant load estimates were generated for manure application fields with APEX and were then subsequently routed to the watershed outlet in SWAT, along with flows and pollutant loadings estimated for the rest of the watershed routed to the watershed outlet. In the current study, the entire UMRW cropland area was simulated in SWAT, which required translating the APEX subareas into SWAT hydrologic response units (HRUs). Calibration and validation of the SWAT output was performed by comparing predicted flow and NO3-N loadings with corresponding in-stream measurements at the watershed outlet from 1999 to 2001. Annual stream flows measured at the watershed outlet were greatly under-predicted when precipitation data collected within the watershed during the 1999-2001 period were used to drive SWAT. Selection of alternative climate data resulted in greatly improved average annual stream predictions, and also relatively strong r2 values of 0.73 and 0.72 for the predicted average monthly flows and NO3-N loads, respectively. The impact of alternative precipitation data shows that as average annual precipitation increases 19%, the relative change in average annual streamflow is about 55%. In summary, the results of this study show that SWAT can replicate measured trends for this watershed and that climate inputs are very important for validating SWAT and other water quality models.

Monoubiquitination and Activity of the Paracaspase MALT1 Requires Glutamate 549 in the Dimerization Interface.

The mucosa-associated lymphoid tissue protein-1 (MALT1, also known as paracaspase) is a protease whose activity is essential for the activation of lymphocytes and the growth of cells derived from human diffuse large B-cell lymphomas of the activated B-cell subtype (ABC DLBCL). Crystallographic approaches have shown that MALT1 can form dimers via its protease domain, but why dimerization is relevant for the biological activity of MALT1 remains largely unknown. Using a molecular modeling approach, we predicted Glu 549 (E549) to be localized within the MALT1 dimer interface and thus potentially relevant. Experimental mutation of this residue into alanine (E549A) led to a complete impairment of MALT1 proteolytic activity. This correlated with an impaired capacity of the mutant to form dimers of the protease domain in vitro, and a reduced capacity to promote NF-κB activation and transcription of the growth-promoting cytokine interleukin-2 in antigen receptor-stimulated lymphocytes. Moreover, this mutant could not rescue the growth of ABC DLBCL cell lines upon MALT1 silencing. Interestingly, the MALT1 mutant E549A was unable to undergo monoubiquitination, which we identified previously as a critical step in MALT1 activation. Collectively, these findings suggest a model in which E549 at the dimerization interface is required for the formation of the enzymatically active, monoubiquitinated form of MALT1.

Dispersal strategies, few dominating or many coexisting: the effect of environmental spatial structure and multiple sources of mortality.

Interspecific competition, life history traits, environmental heterogeneity and spatial structure as well as disturbance are known to impact the successful dispersal strategies in metacommunities. However, studies on the direction of impact of those factors on dispersal have yielded contradictory results and often considered only few competing dispersal strategies at the same time. We used a unifying modeling approach to contrast the combined effects of species traits (adult survival, specialization), environmental heterogeneity and structure (spatial autocorrelation, habitat availability) and disturbance on the selected, maintained and coexisting dispersal strategies in heterogeneous metacommunities. Using a negative exponential dispersal kernel, we allowed for variation of both species dispersal distance and dispersal rate. We showed that strong disturbance promotes species with high dispersal abilities, while low local adult survival and habitat availability select against them. Spatial autocorrelation favors species with higher dispersal ability when adult survival and disturbance rate are low, and selects against them in the opposite situation. Interestingly, several dispersal strategies coexist when disturbance and adult survival act in opposition, as for example when strong disturbance regime favors species with high dispersal abilities while low adult survival selects species with low dispersal. Our results unify apparently contradictory previous results and demonstrate that spatial structure, disturbance and adult survival determine the success and diversity of coexisting dispersal strategies in competing metacommunities.

The three phases of the ensemble forecasting of niche models: geographic range and shifts in climatically suitable areas of Utetheisa ornatrix (Lepidoptera, Arctiidae)

Species' geographic ranges are usually considered as basic units in macroecology and biogeography, yet it is still difficult to measure them accurately for many reasons. About 20 years ago, researchers started using local data on species' occurrences to estimate broad scale ranges, thereby establishing the niche modeling approach. However, there are still many problems in model evaluation and application, and one of the solutions is to find a consensus solution among models derived from different mathematical and statistical models for niche modeling, climatic projections and variable combination, all of which are sources of uncertainty during niche modeling. In this paper, we discuss this approach of ensemble forecasting and propose that it can be divided into three phases with increasing levels of complexity. Phase I is the simple combination of maps to achieve a consensual and hopefully conservative solution. In Phase II, differences among the maps used are described by multivariate analyses, and Phase III consists of the quantitative evaluation of the relative magnitude of uncertainties from different sources and their mapping. To illustrate these developments, we analyzed the occurrence data of the tiger moth, Utetheisa ornatrix (Lepidoptera, Arctiidae), a Neotropical moth species, and modeled its geographic range in current and future climates.

Location models for ceding market share and shrinking services

New location models are presented here for exploring the reduction of facilities in aregion. The first of these models considers firms ceding market share to competitorsunder situations of financial exigency. The goal of this model is to cede the leastmarket share, i.e., retain as much of the customer base as possible while sheddingcostly outlets. The second model considers a firm essentially without competition thatmust shrink it services for economic reasons. This firm is assumed to close outlets sothat the degradation of service is limited. An example is offered within a competitiveenvironment to demonstrate the usefulness of this modeling approach.

In vitro biotransformation of imatinib by the tumor expressed CYP1A1 and CYP1B1.

The main objective of the study was to examine the biotransformation of the anticancer drug imatinib in target cells by incubating it with oxidoreductases expressed in tumor cells. The second objective was to obtain an in silico prediction of the potential activity of imatinib metabolites. An in vitro enzyme kinetic study was performed with cDNA expressed human oxidoreductases and LC-MS/MS analysis. The kinetic parameters (Km and Vmax) were determined for six metabolites. A molecular modeling approach was used to dock these metabolites to the target Abl or Bcr-Abl kinases. CYP3A4 isozyme showed the broadest metabolic capacity, whereas CYP1A1, CYP1B1 and FMO3 isozymes biotransformed imatinib with a high intrinsic clearance. The predicted binding modes for the metabolites to Abl were comparable to that of the parent drug, suggesting potential activity. These findings indicate that CYP1A1 and CYP1B1, which are known to be overexpressed in a wide range of tumors, are involved in the biotransformation of imatinib. They could play a role in imatinib disposition in the targeted stem, progenitor and differentiated cancer cells, with a possible contribution of the metabolites toward the activity of the drug.

A coherent neurobiological framework for functional neuroimaging provided by a model integrating compartmentalized energy metabolism.

Functional neuroimaging has undergone spectacular developments in recent years. Paradoxically, its neurobiological bases have remained elusive, resulting in an intense debate around the cellular mechanisms taking place upon activation that could contribute to the signals measured. Taking advantage of a modeling approach, we propose here a coherent neurobiological framework that not only explains several in vitro and in vivo observations but also provides a physiological basis to interpret imaging signals. First, based on a model of compartmentalized energy metabolism, we show that complex kinetics of NADH changes observed in vitro can be accounted for by distinct metabolic responses in two cell populations reminiscent of neurons and astrocytes. Second, extended application of the model to an in vivo situation allowed us to reproduce the evolution of intraparenchymal oxygen levels upon activation as measured experimentally without substantially altering the initial parameter values. Finally, applying the same model to functional neuroimaging in humans, we were able to determine that the early negative component of the blood oxygenation level-dependent response recorded with functional MRI, known as the initial dip, critically depends on the oxidative response of neurons, whereas the late aspects of the signal correspond to a combination of responses from cell types with two distinct metabolic profiles that could be neurons and astrocytes. In summary, our results, obtained with such a modeling approach, support the concept that both neuronal and glial metabolic responses form essential components of neuroimaging signals.

Light intensity modulates the regulatory network of the shade avoidance response in Arabidopsis.

Plants such as Arabidopsis thaliana respond to foliar shade and neighbors who may become competitors for light resources by elongation growth to secure access to unfiltered sunlight. Challenges faced during this shade avoidance response (SAR) are different under a light-absorbing canopy and during neighbor detection where light remains abundant. In both situations, elongation growth depends on auxin and transcription factors of the phytochrome interacting factor (PIF) class. Using a computational modeling approach to study the SAR regulatory network, we identify and experimentally validate a previously unidentified role for long hypocotyl in far red 1, a negative regulator of the PIFs. Moreover, we find that during neighbor detection, growth is promoted primarily by the production of auxin. In contrast, in true shade, the system operates with less auxin but with an increased sensitivity to the hormonal signal. Our data suggest that this latter signal is less robust, which may reflect a cost-to-robustness tradeoff, a system trait long recognized by engineers and forming the basis of information theory.

Integrating species distribution models (SDMs) and phylogeography for two species of Alpine Primula.

The major intention of the present study was to investigate whether an approach combining the use of niche-based palaeodistribution modeling and phylo-geography would support or modify hypotheses about the Quaternary distributional history derived from phylogeographic methods alone. Our study system comprised two closely related species of Alpine Primula. We used species distribution models based on the extant distribution of the species and last glacial maximum (LGM) climate models to predict the distribution of the two species during the LGM. Phylogeographic data were generated using amplified fragment length polymorphisms (AFLPs). In Primula hirsuta, models of past distribution and phylogeographic data are partly congruent and support the hypothesis of widespread nunatak survival in the Central Alps. Species distribution models (SDMs) allowed us to differentiate between alpine regions that harbor potential nunatak areas and regions that have been colonized from other areas. SDMs revealed that diversity is a good indicator for nunataks, while rarity is a good indicator for peripheral relict populations that were not source for the recolonization of the inner Alps. In P. daonensis, palaeo-distribution models and phylogeographic data are incongruent. Besides the uncertainty inherent to this type of modeling approach (e.g., relatively coarse 1-km grain size), disagreement of models and data may partly be caused by shifts of ecological niche in both species. Nevertheless, we demonstrate that the combination of palaeo-distribution modeling with phylogeographical approaches provides a more differentiated picture of the distributional history of species and partly supports (P. hirsuta) and partly modifies (P. daonensis and P. hirsuta) hypotheses of Quaternary distributional history. Some of the refugial area indicated by palaeodistribution models could not have been identified with phylogeographic data.

Logging residues from regeneration fellings for biofuel production - a GIS-based availability and supply cost analysis

Finland has large forest fuel resources. However, the use of forest fuels for energy production has been low, except for small-scale use in heating. According to national action plans and programs related to wood energy promotion, the utilization of such resources will be multiplied over the next few years. The most significant part of this growth will be based on the utilization of forest fuels, produced from logging residues of regeneration fellings, in industrial and municipal power and heating plants. Availability of logging residues was analyzed by means of resource and demand approaches in order to identify the most suitable regions with focus on increasing the forest fuel usage. The analysis included availability and supply cost comparisons between power plant sites and resource allocation in a least cost manner, and between a predefined power plant structure under demand and supply constraints. Spatial analysis of worksite factors and regional geographies were carried out using the GIS-model environment via geoprocessing and cartographic modeling tools. According to the results of analyses, the cost competitiveness of forest fuel supply should be improved in order to achieve the designed objectives in the near future. Availability and supply costs of forest fuels varied spatially and were very sensitive to worksite factors and transport distances. According to the site-specific analysis the supply potential between differentlocations can be multifold. However, due to technical and economical reasons ofthe fuel supply and dense power plant infrastructure, the supply potential is limited at plant level. Therefore, the potential and supply cost calculations aredepending on site-specific matters, where regional characteristics of resourcesand infrastructure should be taken into consideration, for example by using a GIS-modeling approach constructed in this study.

Modelling the consequences of a reduction in alcohol consumption among patients with alcohol dependence based on real-life observational data.

BACKGROUND: Most available pharmacotherapies for alcohol-dependent patients target abstinence; however, reduced alcohol consumption may be a more realistic goal. Using randomized clinical trial (RCT) data, a previous microsimulation model evaluated the clinical relevance of reduced consumption in terms of avoided alcohol-attributable events. Using real-life observational data, the current analysis aimed to adapt the model and confirm previous findings about the clinical relevance of reduced alcohol consumption. METHODS: Based on the prospective observational CONTROL study, evaluating daily alcohol consumption among alcohol-dependent patients, the model predicted the probability of drinking any alcohol during a given day. Predicted daily alcohol consumption was simulated in a hypothetical sample of 200,000 patients observed over a year. Individual total alcohol consumption (TAC) and number of heavy drinking days (HDD) were derived. Using published risk equations, probabilities of alcohol-attributable adverse health events (e.g., hospitalizations or death) corresponding to simulated consumptions were computed, and aggregated for categories of patients defined by HDDs and TAC (expressed per 100,000 patient-years). Sensitivity analyses tested model robustness. RESULTS: Shifting from >220 HDDs per year to 120-140 HDDs and shifting from 36,000-39,000 g TAC per year (120-130 g/day) to 15,000-18,000 g TAC per year (50-60 g/day) impacted substantially on the incidence of events (14,588 and 6148 events avoided per 100,000 patient-years, respectively). Results were robust to sensitivity analyses. CONCLUSIONS: This study corroborates the previous microsimulation modeling approach and, using real-life data, confirms RCT-based findings that reduced alcohol consumption is a relevant objective for consideration in alcohol dependence management to improve public health.

Will climate change increase the risk of plant invasions into mountains?

Mountain ecosystems have been less adversely affected by invasions of non-native plants than most other ecosystems, partially because most invasive plants in the lowlands are limited by climate and cannot grow under harsher high-elevation conditions. However, with ongoing climate change, invasive species may rapidly move upwards and threaten mid- then high-elevation mountain ecosystems. We evaluated this threat by predicting current and future potential distributions of 48 invasive plant species distributed in Switzerland (CH) and New South Wales (NSW), two areas where climate interacts differently with the elevation gradient. Using a species distribution modeling approach combining two scales, which builds on high-resolution data (< 250 m) but accounts for the global climatic niche of species, we found that different environmental drivers limit the elevation range of invasive species in the two regions, leading to region-specific species responses to climate change. Whereas the optimal suitability for plant invaders is predicted to markedly shift from the lowland to the montane or subalpine zone in CH, such an upward shift is far less pronounced in NSW where montane and subalpine elevations are currently already suitable. Non-native species able to invade the upper reaches of mountains in a future climate will be cold-tolerant in the Swiss Alps but preferring wet soils in the Australian Alps. Other plant traits were only marginally associated with elevation limits. These results demonstrate that a more systematic consideration of future distributions of invasive species is required in conservation plans of not yet invaded mountainous ecosystems.