Systems analysis of host-parasite interactions.

Parasitic diseases caused by protozoan pathogens lead to hundreds of thousands of deaths per year in addition to substantial suffering and socioeconomic decline for millions of people worldwide. The lack of effective vaccines coupled with the widespread emergence of drug-resistant parasites necessitates that the research community take an active role in understanding host-parasite infection biology in order to develop improved therapeutics. Recent advances in next-generation sequencing and the rapid development of publicly accessible genomic databases for many human pathogens have facilitated the application of systems biology to the study of host-parasite interactions. Over the past decade, these technologies have led to the discovery of many important biological processes governing parasitic disease. The integration and interpretation of high-throughput -omic data will undoubtedly generate extraordinary insight into host-parasite interaction networks essential to navigate the intricacies of these complex systems. As systems analysis continues to build the foundation for our understanding of host-parasite biology, this will provide the framework necessary to drive drug discovery research forward and accelerate the development of new antiparasitic therapies

Structural insights into the membrane-extracted dimeric form of the ATPase TraB from the Escherichia coli pKM101 conjugation system

Background: Type IV secretion (T4S) systems are involved in secretion of virulence factors such as toxins or transforming molecules, or bacterial conjugation. T4S systems are composed of 12 proteins named VirB1-B11 and VirD4. Among them, three ATPases are involved in the assembly of the T4S system and/or provide energy for substrate transfer, VirB4, VirB11 and VirD4. The X-ray crystal structures of VirB11 and VirD4 have already been solved but VirB4 has proven to be reluctant to any structural investigation so far. Results: Here, we have used small-angle X-ray scattering to obtain the first structural models for the membrane-extracted, dimeric form of the TraB protein, the VirB4 homolog encoded by the E. coli pKM101 plasmid, and for the monomeric soluble form of the LvhB4 protein, the VirB4 homolog of the T4S system encoded by the Legionella pneumophila lvh operon. We have obtained the low resolution structures of the full-length TraB and of its N- and C-terminal halves. From these SAXS models, we derive the internal organisation of TraB. We also show that the two TraB N- and C-terminal domains are independently involved in the dimerisation of the full-length protein. Conclusions: These models provide the first structural insights into the architecture of VirB4 proteins. In particular, our results highlight the modular arrangement and functional relevance of the dimeric-membrane-bound form of TraB

Predicting the connectivity of primate cortical networks from topological and spatial node properties

The organization of the connectivity between mammalian cortical areas has become a major subject of study, because of its important role in scaffolding the macroscopic aspects of animal behavior and intelligence. In this study we present a computational reconstruction approach to the problem of network organization, by considering the topological and spatial features of each area in the primate cerebral cortex as subsidy for the reconstruction of the global cortical network connectivity. Starting with all areas being disconnected, pairs of areas with similar sets of features are linked together, in an attempt to recover the original network structure. Inferring primate cortical connectivity from the properties of the nodes, remarkably good reconstructions of the global network organization could be obtained, with the topological features allowing slightly superior accuracy to the spatial ones. Analogous reconstruction attempts for the C. elegans neuronal network resulted in substantially poorer recovery, indicating that cortical area interconnections are relatively stronger related to the considered topological and spatial properties than neuronal projections in the nematode. The close relationship between area-based features and global connectivity may hint on developmental rules and constraints for cortical networks. Particularly, differences between the predictions from topological and spatial properties, together with the poorer recovery resulting from spatial properties, indicate that the organization of cortical networks is not entirely determined by spatial constraints

Theoretical basis of the community effect in development

Peer reviewedPublisher PD

The RootScope: A Simple High-Throughput Screening System For Quantitating Gene Expression Dynamics In Plant Roots

Background: High temperature stress responses are vital for plant survival. The mechanisms that plants use to sense high temperatures are only partially understood and involve multiple sensing and signaling pathways. Here we describe the development of the RootScope, an automated microscopy system for quantitating heat shock responses in plant roots.Results: The promoter of Hsp17.6 was used to build a Hsp17.6(p):GFP transcriptional reporter that is induced by heat shock in Arabidopsis. An automated fluorescence microscopy system which enables multiple roots to be imaged in rapid succession was used to quantitate Hsp17.6p: GFP response dynamics. Hsp17.6(p):GFP signal increased with temperature increases from 28 degrees C to 37 degrees C. At 40 degrees C the kinetics and localization of the response are markedly different from those at 37 degrees C. This suggests that different mechanisms mediate heat shock responses above and below 37 degrees C. Finally, we demonstrate that Hsp17.6(p):GFP expression exhibits wave like dynamics in growing roots.Conclusions: The RootScope system is a simple and powerful platform for investigating the heat shock response in plants

Narrative-based computational modelling of the Gp130/JAK/STAT signalling pathway.

BACKGROUND: Appropriately formulated quantitative computational models can support researchers in understanding the dynamic behaviour of biological pathways and support hypothesis formulation and selection by "in silico" experimentation. An obstacle to widespread adoption of this approach is the requirement to formulate a biological pathway as machine executable computer code. We have recently proposed a novel, biologically intuitive, narrative-style modelling language for biologists to formulate the pathway which is then automatically translated into an executable format and is, thus, usable for analysis via existing simulation techniques. RESULTS: Here we use a high-level narrative language in designing a computational model of the gp130/JAK/STAT signalling pathway and show that the model reproduces the dynamic behaviour of the pathway derived by biological observation. We then "experiment" on the model by simulation and sensitivity analysis to define those parameters which dominate the dynamic behaviour of the pathway. The model predicts that nuclear compartmentalisation and phosphorylation status of STAT are key determinants of the pathway and that alternative mechanisms of signal attenuation exert their influence on different timescales. CONCLUSION: The described narrative model of the gp130/JAK/STAT pathway represents an interesting case study showing how, by using this approach, researchers can model biological systems without explicitly dealing with formal notations and mathematical expressions (typically used for biochemical modelling), nevertheless being able to obtain simulation and analysis results. We present the model and the sensitivity analysis results we have obtained, that allow us to identify the parameters which are most sensitive to perturbations. The results, which are shown to be in agreement with existing mathematical models of the gp130/JAK/STAT pathway, serve us as a form of validation of the model and of the approach itself

A thin layer angiogenesis assay: a modified basement matrix assay for assessment of endothelial cell differentiation

BACKGROUND: Basement matrices such as Matrigel™ and Geltrex™ are used in a variety of cell culture assays of anchorage-dependent differentiation including endothelial cell tube formation assays. The volumes of matrix recommended for these assays (approximately 150 μl/cm(2)) are costly, limit working distances for microscopy, and require cell detachment for subsequent molecular analysis. Here we describe the development and validation of a thin-layer angiogenesis (TLA) assay for assessing the angiogenic potential of endothelial cells that overcomes these limitations. RESULTS: Geltrex™ basement matrix at 5 μl/cm(2) in 24-well (10 μl) or 96-well (2 μl) plates supports endothelial cell differentiation into tube-like structures in a comparable manner to the standard larger volumes of matrix. Since working distances are reduced, high-resolution single cell microscopy, including DIC and confocal imaging, can be used readily. Using MitoTracker dye we now demonstrate, for the first time, live mitochondrial dynamics and visualise the 3-dimensional network of mitochondria present in differentiated endothelial cells. Using a standard commercial total RNA extraction kit (Qiagen) we also show direct RNA extraction and RT-qPCR from differentiated endothelial cells without the need to initially detach cells from their supporting matrix. CONCLUSIONS: We present here a new thin-layer assay (TLA) for measuring the anchorage-dependent differentiation of endothelial cells into tube-like structures which retains all the characteristics of the traditional approach but with the added benefit of a greatly lowered cost and better compatibility with other techniques, including RT-qPCR and high-resolution microscopy. ELECTRONIC SUPPLEMENTARY MATERIAL: The online version of this article (doi:10.1186/s12860-014-0041-5) contains supplementary material, which is available to authorized users

A computational model of the hypothalamic - pituitary - gonadal axis in female fathead minnows (Pimephales promelas) exposed to 17α-ethynylestradiol and 17β-trenbolone

© 2011 Li et al; licensee BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.Background - Endocrine disrupting chemicals (e.g., estrogens, androgens and their mimics) are known to affect reproduction in fish. 17α-ethynylestradiol is a synthetic estrogen used in birth control pills. 17β-trenbolone is a relatively stable metabolite of trenbolone acetate, a synthetic androgen used as a growth promoter in livestock. Both 17α-ethynylestradiol and 17β-trenbolone have been found in the aquatic environment and affect fish reproduction. In this study, we developed a physiologically-based computational model for female fathead minnows (FHM, Pimephales promelas), a small fish species used in ecotoxicology, to simulate how estrogens (i.e., 17α-ethynylestradiol) or androgens (i.e., 17β-trenbolone) affect reproductive endpoints such as plasma concentrations of steroid hormones (e.g., 17β-estradiol and testosterone) and vitellogenin (a precursor to egg yolk proteins). Results - Using Markov Chain Monte Carlo simulations, the model was calibrated with data from unexposed, 17α-ethynylestradiol-exposed, and 17β-trenbolone-exposed FHMs. Four Markov chains were simulated, and the chains for each calibrated model parameter (26 in total) converged within 20,000 iterations. With the converged parameter values, we evaluated the model's predictive ability by simulating a variety of independent experimental data. The model predictions agreed with the experimental data well. Conclusions - The physiologically-based computational model represents the hypothalamic-pituitary-gonadal axis in adult female FHM robustly. The model is useful to estimate how estrogens (e.g., 17α-ethynylestradiol) or androgens (e.g., 17β-trenbolone) affect plasma concentrations of 17β-estradiol, testosterone and vitellogenin, which are important determinants of fecundity in fish.The Medical Research Foundation of Oregon, U.S. Environmental Protection Agency, and the National Center for Computational Toxicology of the EPA Office of Research and Development

Functionalization of Bacterial Microcompartment Shell Proteins With Covalently Attached Heme.

Heme is a versatile redox cofactor that has considerable potential for synthetic biology and bioelectronic applications. The capacity to functionalize non-heme-binding proteins with covalently bound heme moieties in vivo could expand the variety of bioelectronic materials, particularly if hemes could be attached at defined locations so as to facilitate position-sensitive processes like electron transfer. In this study, we utilized the cytochrome maturation system I to develop a simple approach that enables incorporation of hemes into the backbone of target proteins in vivo. We tested our methodology by targeting the self-assembling bacterial microcompartment shell proteins, and inserting functional hemes at multiple locations in the protein backbone. We found substitution of three amino acids on the target proteins promoted heme attachment with high occupancy. Spectroscopic measurements suggested these modified proteins covalently bind low-spin hemes, with relative low redox midpoint potentials (about -210 mV vs. SHE). Heme-modified shell proteins partially retained their self-assembly properties, including the capacity to hexamerize, and form inter-hexamer attachments. Heme-bound shell proteins demonstrated the capacity to integrate into higher-order shell assemblies, however, the structural features of these macromolecular complexes was sometimes altered. Altogether, we report a versatile strategy for generating electron-conductive cytochromes from structurally-defined proteins, and provide design considerations on how heme incorporation may interface with native assembly properties in engineered proteins

Structurally robust biological networks

Background: The molecular circuitry of living organisms performs remarkably robust regulatory tasks, despite the often intrinsic variability of its components. A large body of research has in fact highlighted that robustness is often a structural property of biological systems. However, there are few systematic methods to mathematically model and describe structural robustness. With a few exceptions, numerical studies are often the preferred approach to this type of investigation. Results: In this paper, we propose a framework to analyze robust stability of equilibria in biological networks. We employ Lyapunov and invariant sets theory, focusing on the structure of ordinary differential equation models. Without resorting to extensive numerical simulations, often necessary to explore the behavior of a model in its parameter space, we provide rigorous proofs of robust stability of known bio-molecular networks. Our results are in line with existing literature. Conclusions: The impact of our results is twofold: on the one hand, we highlight that classical and simple control theory methods are extremely useful to characterize the behavior of biological networks analytically. On the other hand, we are able to demonstrate that some biological networks are robust thanks to their structure and some qualitative properties of the interactions, regardless of the specific values of their parameters