Probabilistic model checking of complex biological pathways

Probabilistic model checking is a formal verification technique that has been successfully applied to the analysis of systems from a broad range of domains, including security and communication protocols, distributed algorithms and power management. In this paper we illustrate its applicability to a complex biological system: the FGF (Fibroblast Growth Factor) signalling pathway. We give a detailed description of how this case study can be modelled in the probabilistic model checker PRISM, discussing some of the issues that arise in doing so, and show how we can thus examine a rich selection of quantitative properties of this model. We present experimental results for the case study under several different scenarios and provide a detailed analysis, illustrating how this approach can be used to yield a better understanding of the dynamics of the pathway

Improved Continuous Approximation of PEPA Models through Epidemiological Examples

We present two individual based models of disease systems using PEPA (Performance Evaluation Process Algebra). The models explore contrasting mechanisms of disease transmission: direct transmission (e.g. measles) and indirect transmission (e.g. malaria, via mosquitos). We extract ordinary differential equations (ODEs) as a continuous approximation to the PEPA models using the Hillston method and compare these with the traditionally used ODE disease models and with the results of stochastic simulation. Improvements to the Hillston method of ODE extraction for this context are proposed, and the new results compare favourably with stochastic simulation results and to ODEs derived for equivalent models in WSCCS (Weighted Synchronous Calculus of Communicating Systems)

Aldehyde dehydrogenase inhibition blocks mucosal fibrosis in human and mouse ocular scarring

Mucous membrane pemphigoid (MMP) is a systemic mucosal scarring disease, commonly causing blindness, for which there is no antifibrotic therapy. Aldehyde dehydrogenase (ALDH) family 1 is upregulated in both ocular MMP (OMMP) conjunctiva and cultured fibroblasts. Application of the ALDH metabolite, retinoic acid (RA), to normal human conjunctival fibroblasts in vitro induced a diseased phenotype. Conversely, application of ALDH inhibitors, including disulfiram, to OMMP fibroblasts in vitro restored their functionality to that of normal controls. ALDH1 is also upregulated in the mucosa of the mouse model of scarring allergic eye disease, used here as a surrogate for OMMP, in which topical application of disulfiram decreased fibrosis in vivo. These data suggest that progressive scarring in OMMP results from ALDH/RA fibroblast autoregulation, that the ALDH1 subfamily has a central role in immune-mediated ocular mucosal scarring, and that ALDH inhibition with disulfiram is a potential, and readily translatable, antifibrotic therapy

Phase variable DNA repeats in 'Neisseria gonorrhoeae' influence transcription, translation, and protein sequence variation

There are many types of repeated DNA sequences in the genomes of the species of the genus Neisseria, from homopolymeric tracts to tandem repeats of hundreds of bases. Some of these have roles in the phase-variable expression of genes. When a repeat mediates phase variation, reversible switching between tract lengths occurs, which in the species of the genus Neisseria most often causes the gene to switch between on and off states through frame shifting of the open reading frame. Changes in repeat tract lengths may also influence the strength of transcription from a promoter. For phenotypes that can be readily observed, such as expression of the surface-expressed Opa proteins or pili, verification that repeats are mediating phase variation is relatively straightforward. For other genes, particularly those where the function has not been identified, gathering evidence of repeat tract changes can be more difficult. Here we present analysis of the repetitive sequences that could mediate phase variation in the Neisseria gonorrhoeae strain NCCP11945 genome sequence and compare these results with other gonococcal genome sequences. Evidence is presented for an updated phase-variable gene repertoire in this species, including a class of phase variation that causes amino acid changes at the C-terminus of the protein, not previously described in N. gonorrhoeae

Using process algebra to develop predator-prey models of within-host parasite dynamics

As a first approximation of immune-mediated within-host parasite dynamics we can consider the immune response as a predator, with the parasite as its prey. In the ecological literature of predator-prey interactions there are a number of different functional responses used to describe how a predator reproduces in response to consuming prey. Until recently most of the models of the immune system that have taken a predator-prey approach have used simple mass action dynamics to capture the interaction between the immune response and the parasite. More recently Fenton and Perkins (2010) employed three of the most commonly used functional response terms from the ecological literature. In this paper we make use of a technique from computing science, process algebra, to develop mathematical models. The novelty of the process algebra approach is to allow stochastic models of the population (parasite and immune cells) to be developed from rules of individual cell behaviour. By using this approach in which individual cellular behaviour is captured we have derived a ratio-dependent response similar to that seen in previous models of immune-mediated parasite dynamics, confirming that, whilst this type of term is controversial in ecological predator-prey models, it is appropriate for models of the immune system

Influence of FADS Polymorphisms on Tracking of Serum Glycerophospholipid Fatty Acid Concentrations and Percentage Composition in Children

Tracking of fatty acid (FA) contribution to plasma or serum lipids over time was shown in children and adults. However, the potential role of FADS gene variants has not been investigated. Serum GP FA composition of 331 children aged 2 and 6 years, participating in an ongoing birth cohort study, was analyzed. Correlation coefficients were estimated to describe FA tracking over 4 years and to assess the influence of FADS variants on tracking. We found low to moderate tracking (r = 0.12-0.49) of FA compositions and concentration between 2 and 6 years. Concentration changes of total monounsaturated FA and total saturated FA over time correlated closely (r = 0.79) but percentage values were unrelated (r = -0.02). Tracking for n-6 long chain polyunsaturated fatty acid (LC-PUFA) concentrations was lower in subjects homozygous for the major allele of FADS variants and higher in carriers of at least one minor allele, whereas for total n-3 LC-PUFA concentrations and compositions this was vice versa. For individual n-3 PUFA inconsistent results were found. Serum GP FA composition shows low to moderate tracking over 4 years with a higher tracking for LC-PUFA metabolites than for their precursor FA. Serum PUFA levels and their tracking seem to be more influenced by lipid and lipoprotein metabolism than by FA specific pathways

Influence of FADS Polymorphisms on Tracking of Serum Glycerophospholipid Fatty Acid Concentrations and Percentage Composition in Children

BACKGROUND: Tracking of fatty acid (FA) contribution to plasma or serum lipids over time was shown in children and adults. However, the potential role of FADS gene variants has not been investigated. METHODS AND PRINCIPAL FINDINGS: Serum GP FA composition of 331 children aged 2 and 6 years, participating in an ongoing birth cohort study, was analyzed. Correlation coefficients were estimated to describe FA tracking over 4 years and to assess the influence of FADS variants on tracking. We found low to moderate tracking (r = 0.12-0.49) of FA compositions and concentration between 2 and 6 years. Concentration changes of total monounsaturated FA and total saturated FA over time correlated closely (r = 0.79) but percentage values were unrelated (r = -0.02). Tracking for n-6 long chain polyunsaturated fatty acid (LC-PUFA) concentrations was lower in subjects homozygous for the major allele of FADS variants and higher in carriers of at least one minor allele, whereas for total n-3 LC-PUFA concentrations and compositions this was vice versa. For individual n-3 PUFA inconsistent results were found. CONCLUSIONS AND SIGNIFICANCE: Serum GP FA composition shows low to moderate tracking over 4 years with a higher tracking for LC-PUFA metabolites than for their precursor FA. Serum PUFA levels and their tracking seem to be more influenced by lipid and lipoprotein metabolism than by FA specific pathways

Variants of the FADS1 FADS2 Gene Cluster, Blood Levels of Polyunsaturated Fatty Acids and Eczema in Children within the First 2 Years of Life

Association of genetic-variants in the FADS1-FADS2-gene-cluster with fatty-acid-composition in blood of adult-populations is well established. We analyze this genetic-association in two children-cohort-studies. In addition, the association between variants in the FADS-gene-cluster and blood-fatty-acid-composition with eczema was studied. Data of two population-based-birth-cohorts in The Netherlands and Germany (KOALA, LISA) were pooled (n = 879) and analyzed by (logistic) regression regarding the mutual influence of single-nucleotide-polymorphisms (SNPs) in the FADS-gene-cluster (rs174545, rs174546, rs174556, rs174561, rs3834458), on polyunsaturated fatty acids (PUFA) in blood and parent-reported eczema until the age of 2 years. All SNPs were highly significantly associated with all PUFAs except for alpha-linolenic-acid and eicosapentaenoic-acid, also after correction for multiple-testing. All tested SNPs showed associations with eczema in the LISA-study, but not in the KOALA-study. None of the PUFAs was significantly associated with eczema neither in the pooled nor in the analyses stratified by study-cohort. PUFA-composition in young children's blood is under strong control of the FADS-gene-cluster. Inconsistent results were found for a link between these genetic-variants with eczema. PUFA in blood was not associated with eczema. Thus the hypothesis of an inflammatory-link between PUFA and eczema by the metabolic-pathway of LC-PUFAs as precursors for inflammatory prostaglandins and leukotrienes could not be confirmed by these data

Simulations of multiphase turbulence in jet cocoons

M. Krause and P. Alexander, 'Simulations of multiphase turbulence in jet cocoons', Monthly Notices of the Royal Astronomical Society, Vol. 376, pp. 465-478, April 2007, the version of record is available online at doi: 10.1111/j.1365-2966.2007.11480.x. Published by Oxford University Press on behalf of the Royal Astronomical Society. © 2007 The Authors. Journal compilation © 2007 RASThe interaction of optically emitting clouds with warm X-ray gas and hot, tenuous radio plasma in radio jet cocoons is modelled by 2D compressible hydrodynamic simulations. The initial setup is the Kelvin–Helmholtz instability at a contact surface of density contrast 104. The denser medium contains clouds of higher density. Optically thin radiation is realized via a cooling source term. The cool phase effectively extracts energy from the other gas which is both, radiated away and used for acceleration of the cold phase. This increases the system’s cooling rate substantially and leads to a massively amplified cold mass dropout. We show that it is feasible, given small seed clouds of the order of 100 M, that all of the optically emitting gas in a radio jet cocoon may be produced by this mechanism on the propagation time-scale of the jet. The mass is generally distributed as T−1/2 with temperature, with a prominent peak at 14 000 K. This peak is likely to be related to the counteracting effects of shock heating and a strong rise in the cooling function. The volume filling factor of cold gas in this peak is of the order of 10−5–10−3 and generally increases during the simulation time. The simulations tend towards an isotropic scale-free Kolmogorov-type energy spectrum over the simulation time-scale. We find the same Mach-number density relation as Kritsuk & Norman and show that this relation may explain the velocity widths of emission lines associated with high-redshift radio galaxies, if the environmental temperature is lower, or the jet-ambient density ratio is less extreme than in their low-redshift counterparts.Peer reviewe