Global parameter identification of stochastic reaction networks from single trajectories

We consider the problem of inferring the unknown parameters of a stochastic biochemical network model from a single measured time-course of the concentration of some of the involved species. Such measurements are available, e.g., from live-cell fluorescence microscopy in image-based systems biology. In addition, fluctuation time-courses from, e.g., fluorescence correlation spectroscopy provide additional information about the system dynamics that can be used to more robustly infer parameters than when considering only mean concentrations. Estimating model parameters from a single experimental trajectory enables single-cell measurements and quantification of cell--cell variability. We propose a novel combination of an adaptive Monte Carlo sampler, called Gaussian Adaptation, and efficient exact stochastic simulation algorithms that allows parameter identification from single stochastic trajectories. We benchmark the proposed method on a linear and a non-linear reaction network at steady state and during transient phases. In addition, we demonstrate that the present method also provides an ellipsoidal volume estimate of the viable part of parameter space and is able to estimate the physical volume of the compartment in which the observed reactions take place.Comment: Article in print as a book chapter in Springer's "Advances in Systems Biology

Scheme for generating entangled states of two field modes in a cavity

This paper considers a two-level atom interacting with two cavity modes with equal frequencies. Applying a unitary transformation, the system reduces to the analytically solvable Jaynes-Cummings model. For some particular field states, coherent and squeezed states, the transformation between the two bare basis's, related by the unitary transformation, becomes particularly simple. It is shown how to generate, the highly non-classical, entangled coherent states of the two modes, both in the zero and large detuning cases. An advantage with the zero detuning case is that the preparation is deterministic and no atomic measurement is needed. For the large detuning situation a measurement is required, leaving the field in either of two orthogonal entangled coherent states.Comment: Accepted in J. Mod. Opt.; 12 pages; Replaced with revised version. Extended discussion of experimental realizations, earlier studies in the field and on the frequency dependence in the adiabatic eliminatio

An IP-10 (CXCL10)-derived peptide inhibits angiogenesis

Angiogenesis plays a critical role in processes such as organ development, wound healing, and tumor growth. It requires well-orchestrated integration of soluble and matrix factors and timely recognition of such signals to regulate this process. Previous work has shown that newly forming vessels express the chemokine receptor CXC receptor 3 (CXCR3) and, activation by its ligand IP-10 (CXCL10), both inhibits development of new vasculature and causes regression of newly formed vessels. To identify and develop new therapeutic agents to limit or reverse pathological angiogenesis, we identified a 21 amino acid fragment of IP-10, spanning the α-helical domain residues 77-98, that mimic the actions of the whole IP-10 molecule on endothelial cells. Treatment of the endothelial cells with the 22 amino acid fragment referred to as IP-10p significantly inhibited VEGF-induced endothelial motility and tube formation in vitro, properties critical for angiogenesis. Using a Matrigel plug assay in vivo, we demonstrate that IP-10p both prevented vessel formation and induced involution of nascent vessels. CXCR3 neutralizing antibody was able to block the inhibitory effects of the IP-10p, demonstrating specificity of the peptide. Inhibition of endothelial function by IP-10p was similar to that described for IP-10, secondary to CXCR3-mediated increase in cAMP production, activation of PKA inhibiting cell migration, and inhibition of VEGF-mediated m-calpain activation. IP-10p provides a novel therapeutic agent that inhibits endothelial cell function thus, allowing for the modulation of angiogenesis. © 2012 Yates-Binder et al

Climate-driven range extension of Amphistegina (protista, foraminiferida) : models of current and predicted future ranges

© The Author(s), 2013. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in PLoS ONE 8 (2013): e54443, doi:10.1371/journal.pone.0054443.Species-range expansions are a predicted and realized consequence of global climate change. Climate warming and the poleward widening of the tropical belt have induced range shifts in a variety of marine and terrestrial species. Range expansions may have broad implications on native biota and ecosystem functioning as shifting species may perturb recipient communities. Larger symbiont-bearing foraminifera constitute ubiquitous and prominent components of shallow water ecosystems, and range shifts of these important protists are likely to trigger changes in ecosystem functioning. We have used historical and newly acquired occurrence records to compute current range shifts of Amphistegina spp., a larger symbiont-bearing foraminifera, along the eastern coastline of Africa and compare them to analogous range shifts currently observed in the Mediterranean Sea. The study provides new evidence that amphisteginid foraminifera are rapidly progressing southwestward, closely approaching Port Edward (South Africa) at 31°S. To project future species distributions, we applied a species distribution model (SDM) based on ecological niche constraints of current distribution ranges. Our model indicates that further warming is likely to cause a continued range extension, and predicts dispersal along nearly the entire southeastern coast of Africa. The average rates of amphisteginid range shift were computed between 8 and 2.7 km year−1, and are projected to lead to a total southward range expansion of 267 km, or 2.4° latitude, in the year 2100. Our results corroborate findings from the fossil record that some larger symbiont-bearing foraminifera cope well with rising water temperatures and are beneficiaries of global climate change.This work was supported by grants from the German Science Foundation (DFG; www.dfg.de) to ML and SL (LA 884/10-1, LA 884/5-1)

The Space Physics Environment Data Analysis System (SPEDAS)

With the advent of the Heliophysics/Geospace System Observatory (H/GSO), a complement of multi-spacecraft missions and ground-based observatories to study the space environment, data retrieval, analysis, and visualization of space physics data can be daunting. The Space Physics Environment Data Analysis System (SPEDAS), a grass-roots software development platform (www.spedas.org), is now officially supported by NASA Heliophysics as part of its data environment infrastructure. It serves more than a dozen space missions and ground observatories and can integrate the full complement of past and upcoming space physics missions with minimal resources, following clear, simple, and well-proven guidelines. Free, modular and configurable to the needs of individual missions, it works in both command-line (ideal for experienced users) and Graphical User Interface (GUI) mode (reducing the learning curve for first-time users). Both options have “crib-sheets,” user-command sequences in ASCII format that can facilitate record-and-repeat actions, especially for complex operations and plotting. Crib-sheets enhance scientific interactions, as users can move rapidly and accurately from exchanges of technical information on data processing to efficient discussions regarding data interpretation and science. SPEDAS can readily query and ingest all International Solar Terrestrial Physics (ISTP)-compatible products from the Space Physics Data Facility (SPDF), enabling access to a vast collection of historic and current mission data. The planned incorporation of Heliophysics Application Programmer’s Interface (HAPI) standards will facilitate data ingestion from distributed datasets that adhere to these standards. Although SPEDAS is currently Interactive Data Language (IDL)-based (and interfaces to Java-based tools such as Autoplot), efforts are under-way to expand it further to work with python (first as an interface tool and potentially even receiving an under-the-hood replacement). We review the SPEDAS development history, goals, and current implementation. We explain its “modes of use” with examples geared for users and outline its technical implementation and requirements with software developers in mind. We also describe SPEDAS personnel and software management, interfaces with other organizations, resources and support structure available to the community, and future development plans