Synchronization in a lattice model of pulse-coupled oscillators

We analyze the collective behavior of a lattice model of pulse-coupled oscillators. By means of computer simulations we find the relation between the intrinsic dynamics of each member of the population and their mutual interactions that ensures, in a general context, the existence of a fully synchronized regime. This condition turns out to be the same as that obtained for the globally coupled population. When the condition is not completely satisfied we find different spatial structures. This also gives some hints about self-organized criticality.

NUTRIENT RECOMMENDATION MODEL FOR CARROT CROP – FERTICALC CARROT

ABSTRACT The literature on fertilization for carrot growing usually recommends nutrient application rates for yield expectations lower than the yields currently obtained. Moreover, the recommendation only considers the results of soil chemical analysis and does not include effects such as crop residues or variations in yield levels. The aim of this study was to propose a fertilizer recommendation system for carrot cultivation (FERTICALC Carrot) which includes consideration of the nutrient supply by crop residues, variation in intended yield, soil chemical properties, and the growing season (winter or summer). To obtain the data necessary for modeling nutritional requirements, 210 carrot production stands were sampled in the region of Alto Paranaíba, State of Minas Gerais, Brazil. The dry matter content of the roots, the coefficient of biological utilization of nutrients in the roots, and the nutrient harvest index for summer and winter crops were determined for these samples. To model the nutrient supply by the soil, the literature was surveyed in regard to this theme. A modeling system was developed for recommendation of macronutrients and B. For cationic micronutrients, the system only reports crop nutrient export and extraction. The FERTICALC which was developed proved to be efficient for fertilizer recommendation for carrot cultivation. Advantages in relation to official fertilizer recommendation tables are continuous variation of nutrient application rates in accordance with soil properties and in accordance with data regarding the extraction efficiency of modern, higher yielding cultivars.

Defective tumor necrosis factor-alpha-dependent control of astrocyte glutamate release in a transgenic mouse model of Alzheimer disease.

The cytokine tumor necrosis factor-alpha (TNFalpha) induces Ca2+-dependent glutamate release from astrocytes via the downstream action of prostaglandin (PG) E2. By this process, astrocytes may participate in intercellular communication and neuromodulation. Acute inflammation in vitro, induced by adding reactive microglia to astrocyte cultures, enhances TNFalpha production and amplifies glutamate release, switching the pathway into a neurodamaging cascade (Bezzi, P., Domercq, M., Brambilla, L., Galli, R., Schols, D., De Clercq, E., Vescovi, A., Bagetta, G., Kollias, G., Meldolesi, J., and Volterra, A. (2001) Nat. Neurosci. 4, 702-710). Because glial inflammation is a component of Alzheimer disease (AD) and TNFalpha is overexpressed in AD brains, we investigated possible alterations of the cytokine-dependent pathway in PDAPP mice, a transgenic model of AD. Glutamate release was measured in acute hippocampal and cerebellar slices from mice at early (4-month-old) and late (12-month-old) disease stages in comparison with age-matched controls. Surprisingly, TNFalpha-evoked glutamate release, normal in 4-month-old PDAPP mice, was dramatically reduced in the hippocampus of 12-month-old animals. This defect correlated with the presence of numerous beta-amyloid deposits and hypertrophic astrocytes. In contrast, release was normal in cerebellum, a region devoid of beta-amyloid deposition and astrocytosis. The Ca2+-dependent process by which TNFalpha evokes glutamate release in acute slices is distinct from synaptic release and displays properties identical to those observed in cultured astrocytes, notably PG dependence. However, prostaglandin E2 induced normal glutamate release responses in 12-month-old PDAPP mice, suggesting that the pathology-associated defect involves the TNFalpha-dependent control of secretion rather than the secretory process itself. Reduced expression of DENN/MADD, a mediator of TNFalpha-PG coupling, might account for the defect. Alteration of this neuromodulatory astrocytic pathway is described here for the first time in relation to Alzheimer disease.

A model of cellular dosimetry for macroscopic tumors in radiopharmaceutical therapy.

PURPOSE: In the radiopharmaceutical therapy approach to the fight against cancer, in particular when it comes to translating laboratory results to the clinical setting, modeling has served as an invaluable tool for guidance and for understanding the processes operating at the cellular level and how these relate to macroscopic observables. Tumor control probability (TCP) is the dosimetric end point quantity of choice which relates to experimental and clinical data: it requires knowledge of individual cellular absorbed doses since it depends on the assessment of the treatment's ability to kill each and every cell. Macroscopic tumors, seen in both clinical and experimental studies, contain too many cells to be modeled individually in Monte Carlo simulation; yet, in particular for low ratios of decays to cells, a cell-based model that does not smooth away statistical considerations associated with low activity is a necessity. The authors present here an adaptation of the simple sphere-based model from which cellular level dosimetry for macroscopic tumors and their end point quantities, such as TCP, may be extrapolated more reliably. METHODS: Ten homogenous spheres representing tumors of different sizes were constructed in GEANT4. The radionuclide 131I was randomly allowed to decay for each model size and for seven different ratios of number of decays to number of cells, N(r): 1000, 500, 200, 100, 50, 20, and 10 decays per cell. The deposited energy was collected in radial bins and divided by the bin mass to obtain the average bin absorbed dose. To simulate a cellular model, the number of cells present in each bin was calculated and an absorbed dose attributed to each cell equal to the bin average absorbed dose with a randomly determined adjustment based on a Gaussian probability distribution with a width equal to the statistical uncertainty consistent with the ratio of decays to cells, i.e., equal to Nr-1/2. From dose volume histograms the surviving fraction of cells, equivalent uniform dose (EUD), and TCP for the different scenarios were calculated. Comparably sized spherical models containing individual spherical cells (15 microm diameter) in hexagonal lattices were constructed, and Monte Carlo simulations were executed for all the same previous scenarios. The dosimetric quantities were calculated and compared to the adjusted simple sphere model results. The model was then applied to the Bortezomib-induced enzyme-targeted radiotherapy (BETR) strategy of targeting Epstein-Barr virus (EBV)-expressing cancers. RESULTS: The TCP values were comparable to within 2% between the adjusted simple sphere and full cellular models. Additionally, models were generated for a nonuniform distribution of activity, and results were compared between the adjusted spherical and cellular models with similar comparability. The TCP values from the experimental macroscopic tumor results were consistent with the experimental observations for BETR-treated 1 g EBV-expressing lymphoma tumors in mice. CONCLUSIONS: The adjusted spherical model presented here provides more accurate TCP values than simple spheres, on par with full cellular Monte Carlo simulations while maintaining the simplicity of the simple sphere model. This model provides a basis for complementing and understanding laboratory and clinical results pertaining to radiopharmaceutical therapy.

Deciphering genetic disease in the genomic era: the model of GnRH deficiency.

Isolated gonadotropin-releasing hormone (GnRH) deficiency is a treatable albeit rare form of reproductive failure that has revealed physiological mechanisms controlling human reproduction, but despite substantial progress in discovering pathogenic single-gene defects, most of the genetic basis of GnRH deficiency remains uncharted. Although unbiased genetic investigations of affected families have identified mutations in previously unsuspected genes as causes of this disease in some cases, their application has been severely limited because of the negative effect of GnRH deficiency on fertility; moreover, relatively few of the many candidate genes nominated because of biological plausibility from in vitro or animal model experiments were subsequently validated in patients. With the advent of exciting technological platforms for sequencing, homozygosity mapping, and detection of structural variation at the whole-genome level, human investigations are again assuming the leading role for gene discovery. Using human GnRH deficiency as a paradigm and presenting original data from the screening of numerous candidate genes, we discuss the emerging model of patient-focused clinical genetic research and its complementarities with basic approaches in the near future.

What drives invasibility? A multi-model inference test and spatial modelling of alien plant species richness patterns in northern Portugal

Understanding and anticipating biological invasions can focus either on traits that favour species invasiveness or on features of the receiving communities, habitats or landscapes that promote their invasibility. Here, we address invasibility at the regional scale, testing whether some habitats and landscapes are more invasible than others by fitting models that relate alien plant species richness to various environmental predictors. We use a multi-model information-theoretic approach to assess invasibility by modelling spatial and ecological patterns of alien invasion in landscape mosaics and testing competing hypotheses of environmental factors that may control invasibility. Because invasibility may be mediated by particular characteristics of invasiveness, we classified alien species according to their C-S-R plant strategies. We illustrate this approach with a set of 86 alien species in Northern Portugal. We first focus on predictors influencing species richness and expressing invasibility and then evaluate whether distinct plant strategies respond to the same or different groups of environmental predictors. We confirmed climate as a primary determinant of alien invasions and as a primary environmental gradient determining landscape invasibility. The effects of secondary gradients were detected only when the area was sub-sampled according to predictions based on the primary gradient. Then, multiple predictor types influenced patterns of alien species richness, with some types (landscape composition, topography and fire regime) prevailing over others. Alien species richness responded most strongly to extreme land management regimes, suggesting that intermediate disturbance induces biotic resistance by favouring native species richness. Land-use intensification facilitated alien invasion, whereas conservation areas hosted few invaders, highlighting the importance of ecosystem stability in preventing invasions. Plants with different strategies exhibited different responses to environmental gradients, particularly when the variations of the primary gradient were narrowed by sub-sampling. Such differential responses of plant strategies suggest using distinct control and eradication approaches for different areas and alien plant groups.

Generation of biological association networks: A novel strategy to detect new targets in cancer therapy

The aim of this work was to design a novel strategy to detect new targets for anticancer treatments. The rationale was to build Biological Association Networks from differentially expressed genes in drug-resistant cells to identify important nodes within the Networks. These nodes may represent putative targets to attack in cancer therapy, as a way to destabilize the gene network developed by the resistant cells to escape from the drug pressure. As a model we used cells resistant to methotrexate (MTX), an inhibitor of DHFR. Selected node-genes were analyzed at the transcriptional level and from a genotypic point of view. In colon cancer cells, DHFR, the AKR1 family, PKC¿, S100A4, DKK1, and CAV1 were overexpressed while E-cadherin was lost. In breast cancer cells, the UGT1A family was overexpressed, whereas EEF1A1 was overexpressed in pancreatic cells. Interference RNAs directed against these targets sensitized cells towards MTX.

Biological monitoring of chemical exposure : toxicokinetic differences due to age and sex

Human biomonitoring is a widely used method in the assessment of occupational exposure to chemical substances and recommended biological limits are published periodically for interpretation and decision-making. However, it is increasingly recognized that a large variability is associated with biological monitoring, making interpretation less efficient than assumed. In order to improve the applicability of biological monitoring, specific factors responsible for this variability should be identified and their contribution quantified. Among these factors, age and sex are easily identifiable, and present knowledge about pharmaceutical chemicals suggests that they play an important role on the toxicokinetics of occupational chemical agents, and therefore on the biological monitoring results.The aim of the present research project was to assess the influence of age and sex on biological indicators corresponding to organic solvents. This has been done experimentally and by toxicokinetic computer simulation. Another purpose was to explore the effect of selected CYP2E1 polymorphisms on the toxicokinetic profile.Age differences were identified by numerical simulations using a general toxicokinetic model from a previous study which was applied to 14 chemicals, representing 21 specific biological entities, with, among others, toluene, phenol, lead and mercury. These models were runn with the modified parameters, indicating in some cases important differences due to age. The expected changes are mostly of the order of 10-20 %, but differences up to 50 % were observed in some cases. These differences appear to depend on the chemical and on the biological entity considered.Sex differences were quantified by controlled human exposures, which were carried out in a 12 m3 exposure chamber for three organic solvents separately: methyl ethyl ketone, 1-methoxy-2-propanol and 1,1,1-trichloroethane. The human volunteer groups were composed 12 of ten young men and fifteen young women, the latter subdivided into those with and without hormonal contraceptive. They were exposed during six hours at rest and at half of the threshold limit value. The kinetics of the parent compounds (organic volatiles) and their metabolite(s) were followed in blood, urine and expired air over time. Analyses of the solvent and their metabolites were performed by using headspace gas chromatography, CYP2E1 genotypes by using PCR-based RFLP methods. Experimental data were used to calibrate the toxicokinetic models developed for the three solvents. The results obtained for the different biomarkers of exposure mainly showed an effect on the urinary levels of several biomarkers among women due to the use of hormonal contraceptive, with an increase of about 50 % in the metabolism rate. The results also showed a difference due to the genotype CYP2E1*6, when exposed to methyl ethyl ketone, with a tendency to increase CYP2E1 activity when volunteers were carriers of the mutant allele. Simulations showed that it is possible to use simple toxicokinetic tools in order to predict internal exposure when exposed to organic solvents. Our study suggests that not only physiological differences but also exogenous sex hormones could influence CYP2E1 enzyme activity. The variability among the urinary biological indicators levels gives evidence of an interindividual susceptibility, an aspect that should have its place in the approaches for setting limits of occupational exposure.

Neutralizing antibodies in Multiple Sclerosis a model-based analysis of Interferon beta signaling pathway in macrophages

Multiple Sclerosis is the most common non-traumatic cause of neurologicaldisability in young people. There is no cure yet, and until recently, few long-termtherapies existed. Interferon beta (IFNβ) was the first treatment, and remains the mostcommonly prescribed. One of the most significant problems of IFNβ therapy is theproduction of drug specific antibodies. Up to 45% of patients develop neutralizingantibodies (NAbs) to IFNβ products. The neutralizing antibody binds to the biologicalagent preventing its interaction with its receptor, inhibiting the biological action of theprotein, which abrogates the clinical efficacy of IFNβ treatment. Interferon-betamediates its response by binding to its high affinity cell surface receptor and initiatingthe JAK/STAT signalling cascade. In this project we have analyzed the IFNβ signalingpathway in macrophages when neutralizing antibodies are present. The response tothis pathway after IFNβ stimulation shows a transient oscillatory rhythm of STAT1phosphorylation, which varies as NAbs concentration increases. To improve ourunderstanding of that behavior, we extended an existing mathematical model based onnonlinear ordinary differential equations of JAK/STAT pathway by including IFN-NAbassociation and IFN-activation receptor. Combining our theoretical model withexperimental data we could study the role of neutralizing antibodies on the molecularresponse and determine its lifetime after cytokine stimulation.

Comprehensive approach for the detection of antifungal compounds using a susceptible strain of Candida albicans and confirmation of in vivo activity with the Galleria mellonella model.

An efficient screening strategy for the identification of potentially interesting low-abundance antifungal natural products in crude extracts that combines both a sensitive bioautography assay and high performance liquid chromatography (HPLC) microfractionation was developed. This method relies on high performance thin layer chromatography (HPTLC) bioautography with a hypersusceptible engineered strain of Candida albicans (DSY2621) for bioactivity detection, followed by the evaluation of wild type strains in standard microdilution antifungal assays. Active extracts were microfractionated by HPLC in 96-well plates, and the fractions were subsequently submitted to the bioassay. This procedure enabled precise localisation of the antifungal compounds directly in the HPLC chromatograms of the crude extracts. HPLC-PDA-mass spectrometry (MS) data obtained in parallel to the HPLC antifungal profiles provided a first chemical screening about the bioactive constituents. Transposition of the HPLC analytical conditions to medium-pressure liquid chromatography (MPLC) allowed the efficient isolation of the active constituents in mg amounts for structure confirmation and more extensive characterisation of their biological activities. The antifungal properties of the isolated natural products were evaluated by their minimum inhibitory concentration (MIC) in a dilution assay against both wild type and engineered strains of C. albicans. The biological activity of the most promising agents was further evaluated in vitro by electron microscopy and in vivo in a Galleria mellonella model of C. albicans infection. The overall procedure represents a rational and comprehensive means of evaluating antifungal activity from various perspectives for the selection of initial hits that can be explored in more in-depth mode-of-action studies. This strategy is illustrated by the identification and bioactivity evaluation of a series of antifungal compounds from the methanolic extract of a Rubiaceae plant, Morinda tomentosa, which was used as a model in these studies.

Specific inhibition of HCN channels slows rhythm differently in atria, ventricle and outflow tract and stabilizes conduction in the anoxic-reoxygenated embryonic heart model.

The hyperpolarization-activated cyclic nucleotide-gated (HCN) channels are expressed in pacemaker cells very early during cardiogenesis. This work aimed at determining to what extent these channels are implicated in the electromechanical disturbances induced by a transient oxygen lack which may occur in utero. Spontaneously beating hearts or isolated ventricles and outflow tracts dissected from 4-day-old chick embryos were exposed to a selective inhibitor of HCN channels (ivabradine 0.1-10microM) to establish a dose-response relationship. The effects of ivabradine on electrocardiogram, excitation-contraction coupling and contractility of hearts submitted to anoxia (30min) and reoxygenation (60min) were also determined. The distribution of the predominant channel isoform, HCN4, was established in atria, ventricle and outflow tract by immunoblotting. Intrinsic beating rate of atria, ventricle and outflow tract was 164+/-22 (n=10), 78+/-24 (n=8) and 40+/-12bpm (n=23, mean+/-SD), respectively. In the whole heart, ivabradine (0.3microM) slowed the firing rate of atria by 16% and stabilized PR interval. These effects persisted throughout anoxia-reoxygenation, whereas the variations of QT duration, excitation-contraction coupling and contractility, as well as the types and duration of arrhythmias were not altered. Ivabradine (10microM) reduced the intrinsic rate of atria and isolated ventricle by 27% and 52%, respectively, whereas it abolished activity of the isolated outflow tract. Protein expression of HCN4 channels was higher in atria and ventricle than in the outflow tract. Thus, HCN channels are specifically distributed and control finely atrial, ventricular and outflow tract pacemakers as well as conduction in the embryonic heart under normoxia and throughout anoxia-reoxygenation.

Characterization of cerebral glucose dynamics in vivo with a four-state conformational model of transport at the blood-brain barrier.

Determination of brain glucose transport kinetics in vivo at steady-state typically does not allow distinguishing apparent maximum transport rate (T(max)) from cerebral consumption rate. Using a four-state conformational model of glucose transport, we show that simultaneous dynamic measurement of brain and plasma glucose concentrations provide enough information for independent and reliable determination of the two rates. In addition, although dynamic glucose homeostasis can be described with a reversible Michaelis-Menten model, which is implicit to the large iso-inhibition constant (K(ii)) relative to physiological brain glucose content, we found that the apparent affinity constant (K(t)) was better determined with the four-state conformational model of glucose transport than with any of the other models tested. Furthermore, we confirmed the utility of the present method to determine glucose transport and consumption by analysing the modulation of both glucose transport and consumption by anaesthesia conditions that modify cerebral activity. In particular, deep thiopental anaesthesia caused a significant reduction of both T(max) and cerebral metabolic rate for glucose consumption. In conclusion, dynamic measurement of brain glucose in vivo in function of plasma glucose allows robust determination of both glucose uptake and consumption kinetics.

N-acetylcysteine normalizes neurochemical changes in the glutathione-deficient schizophrenia mouse model during development.

BACKGROUND: Glutathione (GSH) is the major cellular redox-regulator and antioxidant. Redox-imbalance due to genetically impaired GSH synthesis is among the risk factors for schizophrenia. Here we used a mouse model with chronic GSH deficit induced by knockout (KO) of the key GSH-synthesizing enzyme, glutamate-cysteine ligase modulatory subunit (GCLM).¦METHODS: With high-resolution magnetic resonance spectroscopy at 14.1 T, we determined the neurochemical profile of GCLM-KO, heterozygous, and wild-type mice in anterior cortex throughout development in a longitudinal study design.¦RESULTS: Chronic GSH deficit was accompanied by an elevation of glutamine (Gln), glutamate (Glu), Gln/Glu, N-acetylaspartate, myo-Inositol, lactate, and alanine. Changes were predominantly present at prepubertal ages (postnatal days 20 and 30). Treatment with N-acetylcysteine from gestation on normalized most neurochemical alterations to wild-type level.¦CONCLUSIONS: Changes observed in GCLM-KO anterior cortex, notably the increase in Gln, Glu, and Gln/Glu, were similar to those reported in early schizophrenia, emphasizing the link between redox imbalance and the disease and validating the model. The data also highlight the prepubertal period as a sensitive time for redox-related neurochemical changes and demonstrate beneficial effects of early N-acetylcysteine treatment. Moreover, the data demonstrate the translational value of magnetic resonance spectroscopy to study brain disease in preclinical models.

Entropic transport in confined media: a challenge for computational studies in biological and soft-matter systems

Transport in small-scale biological and soft-matter systems typically occurs under confinement conditions in which particles proceed through obstacles and irregularities of the boundaries that may significantly alter their trajectories. A transport model that assimilates the confinement to the presence of entropic barriers provides an efficient approach to quantify its effect on the particle current and the diffusion coefficient. We review the main peculiarities of entropic transport and treat two cases in which confinement effects play a crucial role, with the appearance of emergent properties. The presence of entropic barriers modifies the mean first-passage time distribution and therefore plays a very important role in ion transport through micro- and nano-channels. The functionality of molecular motors, modeled as Brownian ratchets, is strongly affected when the motor proceeds in a confined medium that may constitute another source of rectification. The interplay between ratchet and entropic rectification gives rise to a wide variety of dynamical behaviors, not observed when the Brownian motor proceeds in an unbounded medium. Entropic transport offers new venues of transport control and particle manipulation and new ways to engineer more efficient devices for transport at the nanoscale.