Genome scale reconstruction of a Salmonella metabolic model comparison of similarity and differences with a commensal Escherichia coli strain

Salmonella are closely related to commensal Escherichia coli but have gained virulence factors enabling them to behave as enteric pathogens. Less well studied are the similarities and differences that exist between the metabolic properties of these organisms that may contribute toward niche adaptation of Salmonella pathogens. To address this, we have constructed a genome scale Salmonella metabolic model (iMA945). The model comprises 945 open reading frames or genes, 1964 reactions, and 1036 metabolites. There was significant overlap with genes present in E. coli MG1655 model iAF1260. In silico growth predictions were simulated using the model on different carbon, nitrogen, phosphorous, and sulfur sources. These were compared with substrate utilization data gathered from high throughput phenotyping microarrays revealing good agreement. Of the compounds tested, the majority were utilizable by both Salmonella and E. coli. Nevertheless a number of differences were identified both between Salmonella and E. coli and also within the Salmonella strains included. These differences provide valuable insight into differences between a commensal and a closely related pathogen and within different pathogenic strains opening new avenues for future explorations.

Molecular and phenotypic characterization of a mouse model of oculopharyngeal muscular dystrophy reveals severe muscular atrophy restricted to fast glycolytic fibres

Oculopharyngeal muscular dystrophy (OPMD) is an adult-onset disorder characterized by ptosis, dysphagia and proximal limb weakness. Autosomal-dominant OPMD is caused by a short (GCG)8–13 expansions within the first exon of the poly(A)-binding protein nuclear 1 gene (PABPN1), leading to an expanded polyalanine tract in the mutated protein. Expanded PABPN1 forms insoluble aggregates in the nuclei of skeletal muscle fibres. In order to gain insight into the different physiological processes affected in OPMD muscles, we have used a transgenic mouse model of OPMD (A17.1) and performed transcriptomic studies combined with a detailed phenotypic characterization of this model at three time points. The transcriptomic analysis revealed a massive gene deregulation in the A17.1 mice, among which we identified a significant deregulation of pathways associated with muscle atrophy. Using a mathematical model for progression, we have identified that one-third of the progressive genes were also associated with muscle atrophy. Functional and histological analysis of the skeletal muscle of this mouse model confirmed a severe and progressive muscular atrophy associated with a reduction in muscle strength. Moreover, muscle atrophy in the A17.1 mice was restricted to fast glycolytic fibres, containing a large number of intranuclear inclusions (INIs). The soleus muscle and, in particular, oxidative fibres were spared, even though they contained INIs albeit to a lesser degree. These results demonstrate a fibre-type specificity of muscle atrophy in this OPMD model. This study improves our understanding of the biological pathways modified in OPMD to identify potential biomarkers and new therapeutic targets.

Sequential variations of phytoplankton growth and mortality in an NPZ model: a remote-sensing-based assessment

Radiometric data in the visible domain acquired by satellite remote sensing have proven to be powerful for monitoring the states of the ocean, both physical and biological. With the help of these data it is possible to understand certain variations in biological responses of marine phytoplankton on ecological time scales. Here, we implement a sequential data-assimilation technique to estimate from a conventional nutrient–phytoplankton–zooplankton (NPZ) model the time variations of observed and unobserved variables. In addition, we estimate the time evolution of two biological parameters, namely, the specific growth rate and specific mortality of phytoplankton. Our study demonstrates that: (i) the series of time-varying estimates of specific growth rate obtained by sequential data assimilation improves the fitting of the NPZ model to the satellite-derived time series: the model trajectories are closer to the observations than those obtained by implementing static values of the parameter; (ii) the estimates of unobserved variables, i.e., nutrient and zooplankton, obtained from an NPZ model by implementation of a pre-defined parameter evolution can be different from those obtained on applying the sequences of parameters estimated by assimilation; and (iii) the maximum estimated specific growth rate of phytoplankton in the study area is more sensitive to the sea-surface temperature than would be predicted by temperature-dependent functions reported previously. The overall results of the study are potentially useful for enhancing our understanding of the biological response of phytoplankton in a changing environment.

Sensitivity analysis of an ocean carbon cycle model in the North Atlantic: an investigation of parameters affecting the air-sea CO2 flux, primary production and export of detritus

The sensitivity of the biological parameters in a nutrient-phytoplankton-zooplankton-detritus (NPZD) model in the calculation of the air-sea CO2 flux, primary production and detrital export is analysed. We explore the effect on these outputs of variation in the values of the twenty parameters that control ocean ecosystem growth in a 1-D formulation of the UK Met Office HadOCC NPZD model used in GCMs. We use and compare the results from one-at-a-time and all-at-a-time perturbations performed at three sites in the EuroSITES European Ocean Observatory Network: the Central Irminger Sea (60° N 40° W), the Porcupine Abyssal Plain (49° N 16° W) and the European Station for Time series in the Ocean Canary Islands (29° N 15° W). Reasonable changes to the values of key parameters are shown to have a large effect on the calculation of the air-sea CO2 flux, primary production, and export of biological detritus to the deep ocean. Changes in the values of key parameters have a greater effect in more productive regions than in less productive areas. The most sensitive parameters are generally found to be those controlling well-established ocean ecosystem parameterisations widely used in many NPZD-type models. The air-sea CO2 flux is most influenced by variation in the parameters that control phytoplankton growth, detrital sinking and carbonate production by phytoplankton (the rain ratio). Primary production is most sensitive to the parameters that define the shape of the photosynthesis-irradiance curve. Export production is most sensitive to the parameters that control the rate of detrital sinking and the remineralisation of detritus.

A metabolic system-wide characterisation of the pig: a model for human physiology

The pig is a single-stomached omnivorous mammal and is an important model of human disease and nutrition. As such, it is necessary to establish a metabolic framework from which pathology-based variation can be compared. Here, a combination of one and two-dimensional 1H and 13C nuclear magnetic resonance spectroscopy (NMR) and high-resolution magic angle spinning (HR-MAS) NMR was used to provide a systems overview of porcine metabolism via characterisation of the urine, serum, liver and kidney metabolomes. The metabolites observed in each of these biological compartments were found to be qualitatively comparable to the metabolic signature of the same biological matrices in humans and rodents. The data were modelled using a combination of principal components analysis and Venn diagram mapping. Urine represented the most metabolically distinct biological compartment studied, with a relatively greater number of NMR detectable metabolites present, many of which are implicated in gut-microbial co-metabolic processes. The major interspecies differences observed were in the phase II conjugation of extra-genomic metabolites; the pig was observed to conjugate p-cresol, a gut microbial metabolite of tyrosine, with glucuronide rather than sulfate as seen in man. These observations are important to note when considering the translatability of experimental data derived from porcine models.

An energy budget agent-based model of earthworm populations and its application to study the effects of pesticides

Earthworms are important organisms in soil communities and so are used as model organisms in environmental risk assessments of chemicals. However current risk assessments of soil invertebrates are based on short-term laboratory studies, of limited ecological relevance, supplemented if necessary by site-specific field trials, which sometimes are challenging to apply across the whole agricultural landscape. Here, we investigate whether population responses to environmental stressors and pesticide exposure can be accurately predicted by combining energy budget and agent-based models (ABMs), based on knowledge of how individuals respond to their local circumstances. A simple energy budget model was implemented within each earthworm Eisenia fetida in the ABM, based on a priori parameter estimates. From broadly accepted physiological principles, simple algorithms specify how energy acquisition and expenditure drive life cycle processes. Each individual allocates energy between maintenance, growth and/or reproduction under varying conditions of food density, soil temperature and soil moisture. When simulating published experiments, good model fits were obtained to experimental data on individual growth, reproduction and starvation. Using the energy budget model as a platform we developed methods to identify which of the physiological parameters in the energy budget model (rates of ingestion, maintenance, growth or reproduction) are primarily affected by pesticide applications, producing four hypotheses about how toxicity acts. We tested these hypotheses by comparing model outputs with published toxicity data on the effects of copper oxychloride and chlorpyrifos on E. fetida. Both growth and reproduction were directly affected in experiments in which sufficient food was provided, whilst maintenance was targeted under food limitation. Although we only incorporate toxic effects at the individual level we show how ABMs can readily extrapolate to larger scales by providing good model fits to field population data. The ability of the presented model to fit the available field and laboratory data for E. fetida demonstrates the promise of the agent-based approach in ecology, by showing how biological knowledge can be used to make ecological inferences. Further work is required to extend the approach to populations of more ecologically relevant species studied at the field scale. Such a model could help extrapolate from laboratory to field conditions and from one set of field conditions to another or from species to species.

Applying robust control theory to solve problems in bio-medical sciences: study of an apoptotic model

Biological models of an apoptotic process are studied using models describing a system of differential equations derived from reaction kinetics information. The mathematical model is re-formulated in a state-space robust control theory framework where parametric and dynamic uncertainty can be modelled to account for variations naturally occurring in biological processes. We propose to handle the nonlinearities using neural networks.

Mathematical model for NF-kappaB-driven proliferation of adult neural stem cells

Neural stem cells (NSCs) are early precursors of neuronal and glial cells. NSCs are capable of generating identical progeny through virtually unlimited numbers of cell divisions (cell proliferation), producing daughter cells committed to differentiation. Nuclear factor kappa B (NF-kappaB) is an inducible, ubiquitous transcription factor also expressed in neurones, glia and neural stem cells. Recently, several pieces of evidence have been provided for a central role of NF-kappaB in NSC proliferation control. Here, we propose a novel mathematical model for NF-kappaB-driven proliferation of NSCs. We have been able to reconstruct the molecular pathway of activation and inactivation of NF-kappaB and its influence on cell proliferation by a system of nonlinear ordinary differential equations. Then we use a combination of analytical and numerical techniques to study the model dynamics. The results obtained are illustrated by computer simulations and are, in general, in accordance with biological findings reported by several independent laboratories. The model is able to both explain and predict experimental data. Understanding of proliferation mechanisms in NSCs may provide a novel outlook in both potential use in therapeutic approaches, and basic research as well.

A mathematical model for the dynamics of Banana Xanthomonas Wilt with vertical transmission and inflorescence infection

A mathematical model for Banana Xanthomonas Wilt (BXW) spread by insect is presented. The model incorporates inflorescence infection and vertical transmission from the mother corm to attached suckers, but not tool-based transmission by humans. Expressions for the basic reproduction number R0 are obtained and it is verified that disease persists, at a unique endemic level, when R0 > 1. From sensitivity analysis, inflorescence infection rate and roguing rate were the parameters with most influence on disease persistence and equilibrium level. Vertical transmission parameters had less effect on persistence threshold values. Parameters were approximately estimated from field data. The model indicates that single stem removal is a feasible approach to eradication if spread is mainly via inflorescence infection. This requires continuous surveillance and debudding such that a 50% reduction in inflorescence infection and 2–3 weeks interval of surveillance would eventually lead to full recovery of banana plantations and hence improved production.

Comparison of biological behavior between early-stage adenocarcinoma and squamous cell carcinoma of the uterine cervix

Objectives: (1) To compare the anatomopathological variables and recurrence rates in patients with early-stage adenocarcinoma (AC) and squamous cell carcinoma (SCC) of the uterine cervix; (2) to identify the independent risk factors for recurrence. Study design: This historical cohort study assessed 238 patients with carcinoma of the uterine cervix (113 and IIA), who underwent radical hysterectomy with pelvic lymph node dissection between 1980 and 1999. Comparison of category variables between the two histological types was carried out using the Pearson`s X-2 test or Fisher exact test. Disease-free survival rates for AC and SCC were calculated using the Kaplan-Meier method and the curves were compared using the log-rank test. The Cox proportional hazards model was used to identify the independent risk factors for recurrence. Results: There were 35 cases of AC (14.7%) and 203 of SCC (85.3%). AC presented lower histological grade than did SCC (grade 1: 68.6% versus 9.4%; p < 0.001), lower rate of lymphovascular space involvement (25.7% versus 53.7%; p = 0.002), lower rate of invasion into the middle or deep thirds of the uterine cervix (40.0% versus 80.8%; p < 0.001) and lower rate of lymph node metastasis (2.9% versus 16.3%; p = 0.036). Although the recurrence rate was lower for AC than for SCC (11.4% versus 15.8%), this difference was not statistically significant (p = 0.509). Multivariate analysis identified three independent risk factors for recurrence: presence of metastases in the pelvic lymph nodes, invasion of the deep third of the uterine cervix and absence of or slight inflammatory reaction in the cervix. When these variables were adjusted for the histological type and radiotherapy status, they remained in the model as independent risk factors. Conclusion: The AC group showed less aggressive histological behavior than did the SCC group, but no difference in the disease-free survival rates was noted. (C) 2006 Elsevier Ireland Ltd. All rights reserved.

Time-lag in biological responses to landscape changes in a highly dynamic Atlantic forest region

Time-lagged responses of biological variables to landscape modifications are widely recognized, but rarely considered in ecological studies. In order to test for the existence of time-lags in the response of trees, small mammals, birds and frogs to changes in fragment area and connectivity, we studied a fragmented and highly dynamic landscape in the Atlantic forest region. We also investigated the biological correlates associated with differential responses among taxonomic groups. Species richness and abundance for four taxonomic groups were measured in 21 secondary forest fragments during the same period (2000-2002), following a standardized protocol. Data analyses were based on power regressions and model selection procedures. The model inputs included present (2000) and past (1962, 1981) fragment areas and connectivity, as well as observed changes in these parameters. Although past landscape structure was particularly relevant for trees, all taxonomic groups (except small mammals) were affected by landscape dynamics, exhibiting a time-lagged response. Furthermore, fragment area was more important for species groups with lower dispersal capacity, while species with higher dispersal ability had stronger responses to connectivity measures. Although these secondary forest fragments still maintain a large fraction of their original biodiversity, the delay in biological response combined with high rates of deforestation and fast forest regeneration imply in a reduction in the average age of the forest. This also indicates that future species losses are likely, especially those that are more strictly-forest dwellers. Conservation actions should be implemented to reduce species extinction, to maintain old-growth forests and to favour the regeneration process. Our results demonstrate that landscape history can strongly affect the present distribution pattern of species in fragmented landscapes, and should be considered in conservation planning. (C) 2009 Elsevier Ltd. All rights reserved.

BnP1, a novel P-I metalloproteinase from Bothrops neuwiedi venom: Biological effects benchmarking relatively to jararhagin, a P-IIISVMP

Snake venom metalloproteinases (SVMPs) have been extensively studied and their effects associated with the local bleeding observed in human accidents by viper snakes. Representatives of P-I and P-III classes of SVMPs similarly hydrolyze extracellular matrix proteins or coagulation factors while only P-III SVMPs induce significant hemorrhage in experimental models. In this work, the effects of P-I and P-III SVMPs on plasma proteins and cultures of muscle and endothelial cells were compared in order to enlighten the mechanisms involved in venom-induced hemorrhage. To reach this comparison, BnP1 was isolated from B. neuwiedi venom and used as a weakly hemorrhagic P-I SVMPs and jararhagin was used as a model of potently hemorrhagic P-III SVMP. BnP1 was isolated by size exclusion and anion-exchange chromatographies, showing apparent molecular mass of approximately 24kDa and sequence similarity with other members of SVMPs, which allowed its classification as a group P-I SVMP. The comparison of local effects induced by SVMPs showed that BnP1 was devoid of significant myotoxic and hemorrhagic activities and jararhagin presented only hemorrhagic activity. BnP1 and jararhagin were able to hydrolyze fibrinogen and fibrin, although the latter displayed higher activity in both systems. Using HUVEC primary cultures, we observed that BnP1 induced cell detachment and a decrease in the number of viable endothelial cells in levels comparable to those observed by treatment with jararhagin. Moreover, both BnP1 and jararhagin induced apoptosis in HUVECs while only a small increase in LDH supernatant levels was observed after treatment with jararhagin, suggesting that the major mechanism involved in endothelial cell death is apoptosis. Jararhagin and BnP1 induced little effects on C2C12 muscle cell cultures, characterized by a partial detachment 24h after treatment and a mild necrotic effect as evidenced by a small increase in the supernatants LDH levels. Taken together, our data show that P-I and P-III SVMPs presented comparable effects except for the hemorrhagic activity, suggesting that hydrolysis of coagulation factors or damage to endothelial cells are not sufficient for induction of local bleeding. (C) 2007 Elsevier Ltd. All rights reserved.

Autolytic Mycobacterium leprae Hsp65 fragments may act as biological markers for autoimmune diseases

Investigating the proteolytic activity of the recombinant Mycobacterium leprae Heat Shock Protein of 65 kDa (rHsp65), chaperonin 2 (cpn2), we observed that it displays high instability. The fragmentation process starts at the C-terminus followed by progressive degradation of the N-terminus, which leads to a stable fragment comprising the middle region of the molecule. Urea was able to prevent autolysis, probably due to its denaturing action, while EDTA increased degradation levels indicating the need for metal ions. Peptides originated from autolysis were purified and analyzed by mass spectrometry, generating a continuous map. Since the bacteria and mammalian Hsp60 are known to be targets of the immune response and have been implicated in autoimmune diseases and chronic inflammation, the in vivo effect of rHsp65 peptides was evaluated in the spontaneous Systemic Lupus Erythematosus (SLE) model developed by the (NZB/NZW)F(1) mouse hybrids, and their individual anti-rHsp65 IgG2a/IgG1 antibody titer ratio was determined. The results showed orientation toward a T(H)1 responsiveness, and the treatment with the rHsp65 peptides diminished the environmental variance of the survival time of treated animals. These results outline the fact that environmental factors may also act through the modified stability expression of Heat Shock Proteins intervening during autoimmune processes. (C) 2011 Elsevier Ltd. All rights reserved.

Melting Regime of the Anionic Phospholipid DMPG: New Lamellar Phase and Porous Bilayer Model

Aqueous dispersions of the anionic phospholipid dimyristoyl phosphatidylglycerol (DMPG) at pH above the apparent pK of DMPG and concentrations in the interval 70-300 mM have been investigated by small (SAXS) and wide-angle X-ray scattering, differential scanning calorimetry, and polarized optical microscopy. The order. disorder transition of the hydrocarbon chains occurs along an interval of about 10 degrees C (between T(m)(on) similar to 20 degrees C and T(m)(off) similar to 30 degrees C). Such melting regime was previously characterized at lower concentrations, up to 70 mM DMPG, when sample transparency was correlated with the presence of pores across the bilayer. At higher concentrations considered here, the melting regime persists but is not transparent. Defined SAXS peaks appear and a new lamellar phase L(p) with pores is proposed to exist above 70 mM DMPG, starting at similar to 23 degrees C (similar to 3 degrees C above T(m)(on)) and losing correlation after T(m)(off). A new model for describing the X-ray scattering of bilayers with pores, presented here, is able to explain the broad band attributed to in-plane correlation between pores. The majority of cell membranes have a net negative charge, and the opening of pores across the membrane tuned by ionic strength, temperature, and lipid composition is likely to have biological relevance.

Density-Profile Processes Describing Biological Signaling Networks: Almost Sure Convergence to Deterministic Trajectories

We introduce jump processes in R(k), called density-profile processes, to model biological signaling networks. Our modeling setup describes the macroscopic evolution of a finite-size spin-flip model with k types of spins with arbitrary number of internal states interacting through a non-reversible stochastic dynamics. We are mostly interested on the multi-dimensional empirical-magnetization vector in the thermodynamic limit, and prove that, within arbitrary finite time-intervals, its path converges almost surely to a deterministic trajectory determined by a first-order (non-linear) differential equation with explicit bounds on the distance between the stochastic and deterministic trajectories. As parameters of the spin-flip dynamics change, the associated dynamical system may go through bifurcations, associated to phase transitions in the statistical mechanical setting. We present a simple example of spin-flip stochastic model, associated to a synthetic biology model known as repressilator, which leads to a dynamical system with Hopf and pitchfork bifurcations. Depending on the parameter values, the magnetization random path can either converge to a unique stable fixed point, converge to one of a pair of stable fixed points, or asymptotically evolve close to a deterministic orbit in Rk. We also discuss a simple signaling pathway related to cancer research, called p53 module.