Four-Dimensional Neuronal Signaling by Nitric Oxide: A Computational Analysis

Nitric oxide (NO) is now recognized as a transmitter of neurons that express the neuronal isoform of the enzyme nitric oxide synthase. NO, however, violates some of the key tenets of chemical transmission, which is classically regarded as occurring at points of close apposition between neurons. It is the ability of NO to diffuse isotropically in aqueous and lipid environments that has suggested a radically different form of signaling in which the transmitter acts four-dimensionally in space and time, affecting volumes of the brain containing many neurons and synapses. Although ¿volume signaling¿ clearly challenges simple connectionist models of neural processing, crucial to its understanding are the spatial and temporal dynamics of the spread of NO within the brain. Existing models of NO diffusion, however, have serious shortcomings because they represent solutions for ¿point-sources,¿ which have no physical dimensions. Methods for overcoming these difficulties are presented here, and results are described that show how NO spreads from realistic neural architectures with both simple symmetrical and irregular shapes. By highlighting the important influence of the geometry of NO sources, our results provide insights into the four-dimensional spread of a diffusing messenger. We show for example that reservoirs of NO that accumulate in volumes of the nervous system where NO is not synthesized contribute significantly to the temporal and spatial dynamics of NO spread

Engineering simulations for cancer systems biology

Computer simulation can be used to inform in vivo and in vitro experimentation, enabling rapid, low-cost hypothesis generation and directing experimental design in order to test those hypotheses. In this way, in silico models become a scientific instrument for investigation, and so should be developed to high standards, be carefully calibrated and their findings presented in such that they may be reproduced. Here, we outline a framework that supports developing simulations as scientific instruments, and we select cancer systems biology as an exemplar domain, with a particular focus on cellular signalling models. We consider the challenges of lack of data, incomplete knowledge and modelling in the context of a rapidly changing knowledge base. Our framework comprises a process to clearly separate scientific and engineering concerns in model and simulation development, and an argumentation approach to documenting models for rigorous way of recording assumptions and knowledge gaps. We propose interactive, dynamic visualisation tools to enable the biological community to interact with cellular signalling models directly for experimental design. There is a mismatch in scale between these cellular models and tissue structures that are affected by tumours, and bridging this gap requires substantial computational resource. We present concurrent programming as a technology to link scales without losing important details through model simplification. We discuss the value of combining this technology, interactive visualisation, argumentation and model separation to support development of multi-scale models that represent biologically plausible cells arranged in biologically plausible structures that model cell behaviour, interactions and response to therapeutic interventions

A novel pathway producing dimethylsulphide in bacteria is widespread in soil environments

The volatile compound dimethylsulphide (DMS) is important in climate regulation, the sulphur cycle and signalling to higher organisms. Microbial catabolism of the marine osmolyte dimethylsulphoniopropionate (DMSP) is thought to be the major biological process generating DMS. Here we report the discovery and characterisation of the first gene for DMSP-independent DMS production in any bacterium. This gene, mddA, encodes a methyltransferase that methylates methanethiol (MeSH) and generates DMS. MddA functions in many taxonomically diverse bacteria including sediment-dwelling pseudomonads, nitrogen-fixing bradyrhizobia and cyanobacteria, and mycobacteria, including the pathogen Mycobacterium tuberculosis. The mddA gene is present in metagenomes from varied environments, being particularly abundant in soil environments, where it is predicted to occur in up to 76% of bacteria. This novel pathway may significantly contribute to global DMS emissions, especially in terrestrial environments, and could represent a shift from the notion that DMSP is the only significant precursor of DMS

Essential and checkpoint functions of budding yeast ATM and ATR during meiotic prophase are facilitated by differential phosphorylation of a meiotic adaptor protein, Hop1

A hallmark of the conserved ATM/ATR signalling is its ability to mediate a wide range of functions utilizing only a limited number of adaptors and effector kinases. During meiosis, Tel1 and Mec1, the budding yeast ATM and ATR, respectively, rely on a meiotic adaptor protein Hop1, a 53BP1/Rad9 functional analog, and its associated kinase Mek1, a CHK2/Rad53-paralog, to mediate multiple functions: control of the formation and repair of programmed meiotic DNA double strand breaks, enforcement of inter-homolog bias, regulation of meiotic progression, and implementation of checkpoint responses. Here, we present evidence that the multi-functionality of the Tel1/Mec1-to-Hop1/Mek1 signalling depends on stepwise activation of Mek1 that is mediated by Tel1/Mec1 phosphorylation of two specific residues within Hop1: phosphorylation at the threonine 318 (T318) ensures the transient basal level Mek1 activation required for viable spore formation during unperturbed meiosis. Phosphorylation at the serine 298 (S298) promotes stable Hop1-Mek1 interaction on chromosomes following the initial phospho-T318 mediated Mek1 recruitment. In the absence of Dmc1, the phospho-S298 also promotes Mek1 hyper-activation necessary for implementing meiotic checkpoint arrest. Taking these observations together, we propose that the Hop1 phospho-T318 and phospho-S298 constitute key components of the Tel1/Mec1- based meiotic recombination surveillance (MRS) network and facilitate effective coupling of meiotic recombination and progression during both unperturbed and challenged meiosis

Universally Sloppy Parameter Sensitivities in Systems Biology

Quantitative computational models play an increasingly important role in modern biology. Such models typically involve many free parameters, and assigning their values is often a substantial obstacle to model development. Directly measuring \emph{in vivo} biochemical parameters is difficult, and collectively fitting them to other data often yields large parameter uncertainties. Nevertheless, in earlier work we showed in a growth-factor-signaling model that collective fitting could yield well-constrained predictions, even when it left individual parameters very poorly constrained. We also showed that the model had a `sloppy' spectrum of parameter sensitivities, with eigenvalues roughly evenly distributed over many decades. Here we use a collection of models from the literature to test whether such sloppy spectra are common in systems biology. Strikingly, we find that every model we examine has a sloppy spectrum of sensitivities. We also test several consequences of this sloppiness for building predictive models. In particular, sloppiness suggests that collective fits to even large amounts of ideal time-series data will often leave many parameters poorly constrained. Tests over our model collection are consistent with this suggestion. This difficulty with collective fits may seem to argue for direct parameter measurements, but sloppiness also implies that such measurements must be formidably precise and complete to usefully constrain many model predictions. We confirm this implication in our signaling model. Our results suggest that sloppy sensitivity spectra are universal in systems biology models. The prevalence of sloppiness highlights the power of collective fits and suggests that modelers should focus on predictions rather than on parameters.Comment: Submitted to PLoS Computational Biology. Supplementary Information available in "Other Formats" bundle. Discussion slightly revised to add historical contex

Computational analysis of the LRRK2 interactome.

LRRK2 was identified in 2004 as the causative protein product of the Parkinson's disease locus designated PARK8. In the decade since then, genetic studies have revealed at least 6 dominant mutations in LRRK2 linked to Parkinson's disease, alongside one associated with cancer. It is now well established that coding changes in LRRK2 are one of the most common causes of Parkinson's. Genome-wide association studies (GWAs) have, more recently, reported single nucleotide polymorphisms (SNPs) around the LRRK2 locus to be associated with risk of developing sporadic Parkinson's disease and inflammatory bowel disorder. The functional research that has followed these genetic breakthroughs has generated an extensive literature regarding LRRK2 pathophysiology; however, there is still no consensus as to the biological function of LRRK2. To provide insight into the aspects of cell biology that are consistently related to LRRK2 activity, we analysed the plethora of candidate LRRK2 interactors available through the BioGRID and IntAct data repositories. We then performed GO terms enrichment for the LRRK2 interactome. We found that, in two different enrichment portals, the LRRK2 interactome was associated with terms referring to transport, cellular organization, vesicles and the cytoskeleton. We also verified that 21 of the LRRK2 interactors are genetically linked to risk for Parkinson's disease or inflammatory bowel disorder. The implications of these findings are discussed, with particular regard to potential novel areas of investigation

Long-term high fat feeding of rats results in increased numbers of circulating microvesicles with pro-inflammatory effects on endothelial cells

Obesity and type 2 diabetes lead to dramatically increased risks of atherosclerosis and CHD. Multiple mechanisms converge to promote atherosclerosis by increasing endothelial oxidative stress and up-regulating expression of pro-inflammatory molecules. Microvesicles (MV) are small ( < 1 μm) circulating particles that transport proteins and genetic material, through which they are able to mediate cell–cell communication and influence gene expression. Since MV are increased in plasma of obese, insulin-resistant and diabetic individuals, who often exhibit chronic vascular inflammation, and long-term feeding of a high-fat diet (HFD) to rats is a well-described model of obesity and insulin resistance, we hypothesised that this may be a useful model to study the impact of MV on endothelial inflammation. The number and cellular origin of MV from HFD-fed obese rats were characterised by flow cytometry. Total MV were significantly increased after feeding HFD compared to feeding chow (P< 0·001), with significantly elevated numbers of MV derived from leucocyte, endothelial and platelet compartments (P< 0·01 for each cell type). MV were isolated from plasma and their ability to induce reactive oxygen species (ROS) formation and vascular cell adhesion molecule (VCAM)-1 expression was measured in primary rat cardiac endothelial cells in vitro. MV from HFD-fed rats induced significant ROS (P< 0·001) and VCAM-1 expression (P= 0·0275), indicative of a pro-inflammatory MV phenotype in this model of obesity. These findings confirm that this is a useful model to further study the mechanisms by which diet can influence MV release and subsequent effects on cardio-metabolic health

An evolutionary perspective on the kinome of malaria parasites

Malaria parasites belong to an ancient lineage that diverged very early from the main branch of eukaryotes. The approximately 90-member plasmodial kinome includes a majority of eukaryotic protein kinases that clearly cluster within the AGC, CMGC, TKL, CaMK and CK1 groups found in yeast, plants and mammals, testifying to the ancient ancestry of these families. However, several hundred millions years of independent evolution, and the specific pressures brought about by first a photosynthetic and then a parasitic lifestyle, led to the emergence of unique features in the plasmodial kinome. These include taxon-restricted kinase families, and unique peculiarities of individual enzymes even when they have homologues in other eukaryotes. Here, we merge essential aspects of all three malaria-related communications that were presented at the Evolution of Protein Phosphorylation meeting, and propose an integrated discussion of the specific features of the parasite's kinome and phosphoproteome